Prediction of Porous Media Fluid Flow with Spatial Heterogeneity Using Criss-Cross Physics-Informed Convolutional Neural Networks

Recent advances in deep neural networks have shed new light on physics,engineering,and scientific computing.Reconciling the data-centered viewpoint with physical simulation is one of the research hotspots.The physicsinformedneural network(PINN)is currently the most general framework,which is more popular due to theconvenience of constructing NNs and excellent generalization ability.The automatic differentiation(AD)-basedPINN model is suitable for the homogeneous scientific problem;however,it is unclear how AD can enforce fluxcontinuity across boundaries between cells of different properties where spatial heterogeneity is represented bygrid cells with different physical properties.In this work,we propose a criss-cross physics-informed convolutionalneural network(CC-PINN)learning architecture,aiming to learn the solution of parametric PDEs with spatialheterogeneity of physical properties.To achieve the seamless enforcement of flux continuity and integration ofphysicalmeaning into CNN,a predefined 2D convolutional layer is proposed to accurately express transmissibilitybetween adjacent cells.The efficacy of the proposedmethodwas evaluated through predictions of several petroleumreservoir problems with spatial heterogeneity and compared against state-of-the-art(PINN)through numericalanalysis as a benchmark,which demonstrated the superiority of the proposed method over the PINN.

Gas Condensate Two Phase Flow Performance in Porous Media Considering Capillary Number and Non-Darcy Effects

Retrograde condensation frequently occurs during the development of gas condensate reservoirs. The loss of productivity is often observed due to the reduced relative permeability to gas as condensate accumulates near the well bore region. How to describe the condensate blockage effect exactly has been a continuous research topic. However, up to now, the present methods usually over-estimate or underestimate the productivity reduction due to an incorrect understanding of the mechanism of flow in porous medium, which inevitably results in an inaccurate prediction of production performance. It has been found in recent numerous theoretical and experimental studies that capillary number and non-Darcy flow have significant influence on relative permeability in regions near the well bore. The two effects impose opposite impacts on production performance, thus leading to gas condensate flow showing characteristics different from general understanding. It is significant for prediction of performance in gas condensate wells to understand the two effects exactly. The aim of the paper is to describe and analyze the flow dynamics in porous media accurately during the production of gas condensate reservoirs. Based on the description of three-zone flow mechanism, capillary number and non-Darcy effect are incorporated in the analysis of relative permeability, making it possible to describe the effect of condensate blockage. The effect of capillary number and inertial flow on gas and condensate relative permeability is analyzed in detail. Novel Inflow Performance Relation (IPR) models considering high velocity effects are formulated and the contrast analysis of different IPR models is conducted. The result shows that the proposed method can help predict the production performance and productivity more accurately than conventional methods.

MODFLOW-CFPv2 模型在岩溶隧道突涌水及对地下水环境影响中的应用: 以云南鹤庆锰矿沟岩溶水系统为例

滇中地区构造复杂、岩溶发育,隧洞突涌水及泉流量衰减是隧洞施工过程中最棘手的问题之一。锰矿沟岩溶水系统岩溶管道化程度高,岩溶裂隙与岩溶管道2种含水介质差异显著。采用MODFLOW-CFPv2双重介质数值模型对锰矿沟岩溶水系统展开数值模拟研究,精细刻画岩溶管道与引水隧洞,进而掌握隧洞施工对地下水流场影响以及泉流量变化的规律。结果表明:(1)MODFLOW-CFPv2模型能够刻画岩溶地区复杂的地质结构,较好地模拟研究区地下水位的动态特征和岩溶泉流量响应特征。(2)在强排工况下隧洞单位长度最大涌水量为164 m^(3)/(d·m),单位长度稳定涌水量为69 m^(3)/(d·m),锰矿沟岩溶泉流量也出现显著下降的趋势,在模拟期内平均泉流量从天然条件下1578 L/s下降至1098 L/s,总体减少了30.4%;峰值泉流量从2133 L/s下降至1375 L/s,减少了35.5%,强排工况施工会对隧洞工程施工和地下水环境造成显著影响;限排工况下隧洞单位长度最大涌水量为39 m^(3)/(d·m),单位长度稳定涌水量为24 m^(3)/(d·m),隧洞单位长度涌水量显著降低,锰矿沟岩溶泉流量的下降趋势也得到了一定程度的改善,模拟期内平均泉流量降低至1284 L/s,减少了18.6%,峰值泉流量降低至1617 L/s,减少了22.1%。采用的MODFLOW-CFPv2双重介质模型具有较精确刻画岩溶区管道、溶洞、裂隙共存的高度非均质岩溶水系统的能力,能够定量评价香炉山隧洞施工对锰矿沟岩溶水系统地下水流场及泉流量的影响,为香炉山隧洞工程的突涌水灾害防治提供参考依据,也为岩溶地区复杂地质条件下地下水研究提供借鉴经验。

GLOBAL EXISTENCE AND BLOW-UP PHENOMENA OF CLASSICAL SOLUTIONS FOR THE SYSTEM OF COMPRESSIBLE ADIABATIC FLOW THROUGH POROUS MEDIA

By means of maximum principle for nonlinear hyperbolic systems, the results given by HSIAO Ling and D. Serre was improved for Cauchy problem of compressible adiabatic flow through porous media, and a complete result on the global existence and the blow-up phenomena of classical solutions of these systems. These results show that the dissipation is strong enough to preserve the smoothness of ‘small ’ solution.

Evaluation of gas wettability and its effects on fluid distribution and fluid flow in porous media

The special gas wettability phenomenon of reservoir rocks has been recognized by more and more researchers.It has a significant effect on efficient development of unconventional reservoirs.First,based on the preferentially gas-covered ability and surface free energy changes,definition and evaluation methods have been established.Second,a method for altering rock wettability and its mechanisms have been studied,surface oriented phenomena of functional groups with low surface energy are the fundamental reason for gas wettability alteration of rock.Third,the effect of gas wettability on the surface energy,electrical properties and dilatability are investigated.Last,the effects of gas wettability on capillary pressure,oil/gas/water distribution and flow are investigated with capillary tubes and etchedglass network models.The gas wettability theory of reservoir rocks has been initially established,which provides theoretical support for the efficient production of unconventional reservoirs and has great significance.

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.

Numerical simulation of two-phase flow in fractured porous media using streamline simulation and IMPES methods and comparing results with a commercial software

Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consuming. In streamline method, transport equations are solved on one-dimensional streamlines to reduce the computation time with less memory for simulation. First, pressure equation is solved on an Eulerian grid and streamlines are traced. Defining the "time of flight", saturation equations are mapped and solved on streamlines. Finally, the results are mapped back on Eulerian grid and the process is repeated until the simulation end time. The waterflooding process is considered in a fractured reservoir using the dual porosity model. Afterwards, a computational code is developed to solve the same problem by the IMPES method and the results of streamline simulation are compared to those of the IMPES and a commercial software. Finally, the accuracy and efficiency of streamline simulator for simulation of two-phase flow in fractured reservoirs has been proved.

A control volume based finite element method for simulating incompressible two-phase flow in heterogeneous porous media and its application to reservoir engineering

Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.

Interactions between gas flow and reversible chemical reaction in porous media

Taking into consideration the gas compressibility and chemical reaction reversibility, a model was developed to study the interactions between gas flow and chemical reaction in porous media and resolved by the finite volume method on the basis of the gas-solid reaction aA(g)+bB(s)cC(g)+dD(s).The numerical analysis shows that the equilibrium constant is an important factor influencing the process of gas-solid reaction. The stoichiometric coefficients, molar masses of reactant gas, product gas and inert gas are the main factors influencing the density of gas mixture. The equilibrium constant influences the gas flow in porous media obviously when the stoichiometric coefficients satisfy a/c≠1.

A numerical study on heat transfer enhancement and design of a heat exchanger with porous media in continuous hydrothermal flow synthesis system

The aim of this study is to use a new configuration of porous media in a heat exchanger in continuous hydrothermal flow synthesis(CHFS)system to enhance the heat transfer and minimize the required length of the heat exchanger.For this purpose,numerous numerical simulations are performed to investigate performance of the system with porous media.First,the numerical simulation for the heat exchanger in CHFS system is validated by experimental data.Then,porous media is added to the system and six different thicknesses for the porous media are examined to obtain the optimum thickness,based on the minimum required length of the heat exchanger.Finally,by changing the flow rate and inlet temperature of the product as well as the cooling water flow rate,the minimum required length of the heat exchanger with porous media for various inlet conditions is assessed.The investigations indicate that using porous media with the proper thickness in the heat exchanger increases the cooling rate of the product by almost 40% and reduces the required length of the heat exchanger by approximately 35%.The results also illustrate that the most proper thickness of the porous media is approximately equal to 90% of the product tube's thickness.Results of this study lead to design a porous heat exchanger in CHFS system for various inlet conditions.

An overview on nonlinear porous flow in low permeability porous media

This paper gives an overview on nonlinear porous flow in low permeability porous media, reveals the microscopic mechanisms of flows, and clarifies properties of porous flow fluids. It shows that, deviating from Darcy＇s linear law, the porous flow characteristics obey a nonlinear law in a low-permeability porous medium, and the viscosity of the porous flow fluid and the permeability values of water and oil are not constants. Based on these characters, a new porous flow model, which can better describe low permeability reservoir~ is established. This model can describe various patterns of porous flow, as Darcy＇s linear law does. All the parameters involved in the model, having definite physical meanings, can be obtained directly from the experiments.

Measurement of Fluid Flow in Pipe and Porous Media by High-Resolution Magnetic Resonance Imaging

The objective of this study is to understand the process of fluid flow in pipe and porous media with different pore structures. High-resolution Magnetic Resonance Imaging （MRI） technique was used to visualize the pore structure and measure fluid flow. The porous media was formed by packed bed of glass beads. Flow measurement was carried out by a modified spin echo sequence. The results show that the velocity distribution in pipe is annular and the linear relation between MRI velocity and actual velocity is found in pipe flow measurement. The flow distribution in porous media is rather heterogeneous, and it is consistent with heterogeneous pore structure. The flow through pores with the high volume flow rate is determined largely by geometrical effects such as pore size and cross-sectional area.

GLOBAL EXISTENCE AND BLOW-UP PHENOMENA OF CLASSICAL SOLUTIONS FOR THE SYSTEM OF COMPRESSIBLE ADIABATIC FLOW THROUGH POROUS MEDIA

By means of maximum principle for nonlinear hyperbolic systems, the results given by HSIAO Ling and D. Serre was improved for Cauchy problem of compressible adiabatic flow through porous media, and a complete result on the global existence and the blow-up phenomena of classical solutions of these systems. These results show that the dissipation is strong enough to preserve the smoothness of ‘small ’ solution.

TWO PHASE COMPRESSIBLE FLOW IN POROUS MEDIA

We study the mathematical model of two phase compressible flows through porous media. Under the condition that the compressibility of rock, oil, and water is small, we prove that the initial-boundary value problem of the nonlinear system of equations admits a weak solution.

Enduring effect of permeability texture for enhancing accuracy and reducing uncertainty of reservoir fluid flow through porous media

Modeling reservoir permeability is one of the crucial tasks in reservoir simulation studies.Traditionally,it is done by kriging-based methods.More rigorous modeling of the permeability results in more reliable outputs of the reservoir models.Recently,a new category of geostatistical methods has been used for this purpose,namely multiple point statistics(MPS).By this new category of permeability modeling methods,one is able to predict the heterogeneity of the reservoir permeability as a continuous variable.These methods consider the direction of property variation in addition to the distances of known locations of the property.In this study,the reservoir performance of a modified version of the SPE 10 solution project as a pioneer case is used for investigating the efficiency of these methods and paralleling them with the kriging-based one.In this way,the permeability texture concept is introduced by applying some MPS methods.This study is accomplished in the conditions of real reservoir dimensions and velocities for the whole reservoir life.A continuous training image is used as the input of calculation for the permeability modeling.The results show that the detailed permeability of the reservoir as a continuous variable makes the reservoir simulation show the same fluid front movement and flooding behavior of the reservoir similar to the reference case with the same permeability heterogeneity.Some MPS methods enable the reservoir simulation to reproduce the fluid flow complexities such as bypassing and oil trapping during water flooding similar to the reference case.Accordingly,total oil production is predicted with higher accuracy and lower uncertainty.All studied cases are identical except for the permeability texture.Even histograms and variograms of permeabilities for the studied reservoir are quite similar,but the performance of the reservoir shows that kriging-based method results have slightly less accuracy than some MPS methods.Meanwhile,it results in lower uncertainty in outputs for this water flooding case performance.

Digital image processing of saturation for two-phase flow in planar porous media model

In this paper, the accuracy of estimating stained non-wetting phase saturation using digital image processing is examined, and a novel post-processing approach for calculating threshold is presented. In order to remove the effect of the background noise of images and to enhance the high-frequency component of the original image, image smoothing and image sharpening methods are introduced. Depending on the correct threshold, the image binarization processing is particularly useful for estimating stained non-wetting phase saturation. Calculated saturation data are compared with the measured saturation data during the two-phase flow experiment in an artificial steel planar porous media model. The results show that the calculated saturation data agree with the measured ones. With the help of an artificial steel planar porous media model, digital image processing is an accurate and simple method for obtaining the stained non-wetting phase saturation.

Plane Stagnation Double-Diffusive MHD Convective Flow with Convective Boundary Condition in a Porous Media

A numerical analysis has been carried out to study the problem of plane stagnation double-diffusive MHD convective flow with convective boundary condition in a porous media. The governing nonlinear partial differential equations have been reduced to systems of nonlinear ordinary differential equations by the similarity transformations. The transformed equations are solved numerically by using the classical fourth order Runge-Kutta method together with the shooting technique implemented on a computer program. The effects of the physical parameters are examined on the velocity, temperature and concentration profiles. Numerical data for the skin-friction coefficients, Nusselt and Sherwood numbers have been tabulated for various parametric conditions and are also shown graphically and discussed.

Determining factors for carotid mean/max intima-media thickness and brachial flow-mediated dilation in healthy young women

Background: Many factors can contribute to atherosclerotic-type vascular changes in older individuals or men. Thus, confining the investigation to young women with no clinical evidence of the condition could enhance understanding of the early stages of cardiovascular disease. The aim of this study was to determine whether carotid mean/max intima-media thickness (IMT) and brachial flow-mediated dilation (FMD) values, which are well-known event-related indices, are associated with laboratory data and the other vascular indices of atherosclerosis in healthy young women. Methods: Carotid mean/max IMT and brachial FMD were measured in young women with no clinical evidence of atherosclerosis (n = 110;mean age, 39 years) who were instructed not to eat, drink or smoke after 9 PM the evening before testing. All participants also underwent laboratory assessment, including simultaneous measurements of arterial stiffness such as augmentation index (AI), cardioankle vascular index (CAVI) and brachial-ankle pulse wave velocity (baPWV). Results: Mean IMT was signifi-cantly and positively associated with age (p = 0.002), CAVI (p = 0.044), low-density lipoprotein-cholesterol (LDL-C, p = 0.047) and high-sensitive C-reactive protein (hs-CRP, p = 0.002) values but was not related to FMD, AI, baPWV or triglycerides (TG) in the multivariate regression analysis. Similarly, max IMT was positively associated with age (p p = 0.003) and hs-CRP (p = 0.005) values but was not related to FMD, AI, CAVI, baPWV, TG or blood pressure level in the multivariate regression analysis. The association between LDL-C and max IMT was much stronger than that between LDL-C and mean IMT. Brachial FMD was positively associated only with heart rate in the multivariate regression analysis. Conclusions: These results suggest that mean IMT more closely represents the sclerotic aspect of vascular change, whereas max IMT represents the atherotic aspect in healthy young women. Although the relationship between the autonomic nervous system and heart rate is well-known, there may be a complex interaction between the autonomic nervous system and endothelial function.

Plane Transverse MHD Flow through Porous Media

Plane, transverse MHD flow through a porous structure is considered in this work. Solution to the governing equations is obtained using an inverse method in which the streamfunction of the flow is considered linear in one of the space variables. Expressions for flow quantities are obtained for finitely conducting and infinitely conducting fluids.

Unbalanced International News Flow and News Bias in Western Media

With the rapid development of mass media, a new media age comes into being. Accompanied with the rising of advanced capitalism, the print media brings us a new“mass sphere”. In the field of international journalism and communication, it is still exist that international news is unbalanced and biased in first world countries’media. However, the imbalance of international public opinion is just a big problem for third world countries. Therefore, it is very necessary and significant to analyze this subject. Starting from the globalization, this paper argues the imbalance and news bias of reporting upon international issues in western media, so as to provide some ideas and promote the better development of the subject.