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.

Modeling of multiphase flow in low permeability porous media:Effect of wettability and pore structure properties

Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery.

Flow characteristics and regime transition of aqueous foams in porous media over a wide range of quality,velocity,and surfactant concentration

Aqueous foam is broadly applicable to enhanced oil recovery(EOR).The rheology of foam as a function of foam quality,gas and liquid velocities,and surfactant concentration constitute the foundation of its application.The great variations of the above factors can affect the effectiveness of N2 foam in EOR continuously in complex formations,which is rarely involved in previous relevant studies.This paper presents an experimental study of foam flow in porous media by injecting pre-generated N2 foam into a sand pack under the conditions of considering a wide range of gas and liquid velocities and surfactant concentrations.The results show that in a wide range of gas and liquid velocities,the pressure gradient contours are L-shaped near the coordinate axes,but V-shaped in other regions.And the surfactant concentration is a strong factor influencing the trend of pressure gradient contours.Foam flow resistance is very sensitive to the surfactant concentration in both the high-and low-foam quality regime,especially when the surfactant concentration is less than CMC.The foam quality is an important variable to the flow resistance obtained.There exists a transition point from low-to high-quality regime in a particular flow system,where has the maximum flow resistance,the corresponding foam quality is called transition foam quality,which increases as the surfactant concentration increases.The results can add to our knowledge base of foam rheology in porous media,and can provide a strong basis for the field application of foams.

Dynamic behavior of coalbed methane flow along the annulus of single-phase production

Dynamic behavior of coalbed methane (CBM) flow will provide the theoretical basis to optimize production performance for a given well.A mathematical model is developed to simulate flowing pressures and pressure drops of CBM column from well head to bottom hole.The measured parameters and independent variables of flow rates,flowing pressures and temperatures are involved in CBM producing process along the annulus.The developed relationships are validated against full-scale measured data in single-phase CBM wellbores.The proposed methodology can analyze the dynamic behavior in CBM reservoir and process of CBM flow with an overall accuracy of 2%.The calculating process of flowing pressures involves friction factor with variable Reynolds number and CBM temperature and compressibility factor with gravitational gradients.The results showed that the effect of flowing pressure on CBM column was more obvious than that on CBM and water column accompanied by an increase of dynamic water level.The ratios of flowing pressure on increment of CBM column to the whole column increased with the declined flow rates of water column.Bottom-hole pressure declined with the decreased flowing pressure of CBM column along the annulus.It will lead to the results of the increased pressure drop of CBM column and CBM flow rate in single-phase CBM wellbores.

Flow and heat transfer of a nanofluid through a porous medium due to stretching/shrinking sheet with suction,magnetic field and thermal radiation

This study investigates the suction and magnetic field effects on the two-dimensional nanofluid flow through a stretching/shrinking sheet at the stagnation point in the porous medium with thermal radiation.The governing partial differential equations(PDEs)are converted into ordinary differential equations(ODEs)using the similarity transformation.The resulting ODEs are then solved numerically by using the bvp4c solver in MATLAB software.It was found that dual solutions exist for the shrinking parameter values up to a certain range.The numerical results obtained are compared,and the comparison showed a good agreement with the existing results in the literature.The governing parameters’effect on the velocity,temperature and nanoparticle fraction fields as well as the skin friction coefficient,the local Nusselt number and the Sherwood number are represented graphically and analyzed.The variation of the velocity,temperature and concentration increase with the increase in the suction and magnetic field parameters.It seems that the thermal radiation effect has increased the local Sherwood number while the local Nusselt number is reduced with it.

The capillary pressure curves from upscaling interfacial and unsaturated flows in porous layers with vertical heterogeneity

We provide the capillary pressure curves p_(c)(s)as a function of the effective saturation s based on the theoretical framework of upscaling unsaturated flows in vertically heterogeneous porous layers proposed recently(Z.Zheng,Journal of Fluid Mechanics,950,A17,2022).Based on the assumption of vertical gravitational-capillary equilibrium,the saturation distribution and profile shape of the invading fluid can be obtained by solving a nonlinear integral-differential equation.The capillary pressure curves p_(c)(s)can then be constructed by systematically varying the injection rate.Together with the relative permeability curves k_(rn)(s)that are already obtained.One can now provide quick estimates on the overall behaviours of interfacial and unsaturated flows in vertically-heterogeneous porous layers.

Numerical Approach of a Coupled Pressure-Saturation Model Describing Oil-Water Flow in Porous Media

Two-phase flow in porous media is a very active field of research,due to its important applications in groundwater pollution,CO_(2)sequestration,or oil and gas production from petroleum reservoirs,just to name a few of them.Fractional flow equations,which make use of Darcy's law,for describing the movement of two immiscible fluids in a porous medium,are among the most relevant mathematical models in reservoir simulation.This work aims to solve a fractional flow model formed by an elliptic equation,representing the spatial distribution of the pressure,and a hyperbolic equation describing the space-time evolution of water saturation.The numerical solution of the elliptic part is obtained using a finite-element(FE)scheme,while the hyperbolic equation is solved by means of two dif-ferent numerical approaches,both in the finite-volume(FV)framework.One is based on a monotonic upstream-centered scheme for conservation laws(MUSCL)-Hancock scheme,whereas the other makes use of a weighted essentially non-oscillatory(ENO)reconstruc-tion.In both cases,a first-order centered(FORCE)-αnumerical scheme is applied for inter-cell flux reconstruction,which constitutes a new contribution in the field of fractional flow models describing oil-water movement.A relevant feature of this work is the study of the effect of the parameterαon the numerical solution of the models considered.We also show that,in the FORCE-αmethod,when the parameterαincreases,the errors diminish and the order of accuracy is more properly attained,as verified using a manufactured solution technique.

Soft Template-Induced Porous Polyvinylidene Fluoride Membrane for Vanadium Flow Batteries

Vanadium flow batteries(VFBs)are considered ideal for grid-sc ale,long-duration energy storage applications owing to their decoupled output power and storage capacity,high safety,efficiency,and long cycle life.However,the widespread adoption of VFB s is hindered by the use of expensive Nafion membranes.Herein,we report a soft template-induced method to develop a porous polyvinylidene fluoride(PVDF)membrane for VFB applications.By incorporating water-soluble and flexible polyethylene glycol(PEG 400)as a soft template,we induced the aggregation of hydrophilic sulfonated poly(ether ether ketone),resulting in phase separation from the hydrophobic PVDF polymer during membrane formation.This process led to the creation of a porous PVDF membrane with controllable morphologies determined by the polyethylene glycol content in the cast solution.The optimized porous PVDF membrane enabled a stable VFB performance for 200 cycles at a current density of 80 mA/cm^(2),and the VFB exhibited a Coulombic efficiency of 95.2%and a voltage efficiency of 87.8%.These findings provide valuable insights for the development of highly stable membranes for VFB applications.

P-and SV-wave dispersion and attenuation in saturated microcracked porous rock with aligned penny-shaped fractures

P-and SV-wave dispersion and attenuation have been extensively investigated in saturated poroelastic media with aligned fractures.However,there are few existing models that incorporate the multiple wave attenuation mechanisms from the microscopic scale to the macroscopic scale.Hence,in this work,we developed a unified model to incorporate the wave attenuation mechanisms at different scales,which includes the microscopic squirt flow between the microcracks and pores,the mesoscopic wave-induced fluid flow between fractures and background(FB-WIFF),and the macroscopic Biot's global flow and elastic scattering(ES)from the fractures.Using Tang's modified Biot's theory and the mixed-boundary conditions,we derived the exact frequency-dependent solutions of the scattering problem for a single penny-shaped fracture with oblique incident P-and SV-waves.We then developed theoretical models for a set of aligned fractures and randomly oriented fractures using the Foldy approximation.The results indicated that microcrack squirt flow considerably influences the dispersion and attenuation of P-and SV-wave velocities.The coupling effects of microcrack squirt flow with the FB-WIFF and ES of fractures cause much higher velocity dispersion and attenuation for P waves than for SV waves.Randomly oriented fractures substantially reduce the attenuation caused by the FB-WIFF and ES,particularly for the ES attenuation of SV waves.Through a comparison with existing models in the limiting cases and previous experimental measurements,we validated our model.

Anisotropic dynamic permeability model for porous media

Based on the tortuous capillary network model,the relationship between anisotropic permeability and rock normal strain,namely the anisotropic dynamic permeability model(ADPM),was derived and established.The model was verified using pore-scale flow simulation.The uniaxial strain process was calculated and the main factors affecting permeability changes in different directions in the deformation process were analyzed.In the process of uniaxial strain during the exploitation of layered oil and gas reservoirs,the effect of effective surface porosity on the permeability in all directions is consistent.With the decrease of effective surface porosity,the sensitivity of permeability to strain increases.The sensitivity of the permeability perpendicular to the direction of compression to the strain decreases with the increase of the tortuosity,while the sensitivity of the permeability in the direction of compression to the strain increases with the increase of the tortuosity.For layered reservoirs with the same initial tortuosity in all directions,the tortuosity plays a decisive role in the relative relationship between the variations of permeability in all directions during pressure drop.When the tortuosity is less than 1.6,the decrease rate of horizontal permeability is higher than that of vertical permeability,while the opposite is true when the tortuosity is greater than 1.6.This phenomenon cannot be represented by traditional dynamic permeability model.After the verification by experimental data of pore-scale simulation,the new model has high fitting accuracy and can effectively characterize the effects of deformation in different directions on the permeability in all directions.

A quasi single-phase model for debris flows and its comparison with a two-phase model

A depth-averaged quasi single-phase mixture model is proposed for debris flows over inclined bed slopes based on the shallow water hydrosediment-morphodynamic theory with multi grain sizes. The stresses due to fluctuations are incorporated based on analogy to turbulent flows, as estimated using the depth-averaged k-? turbulence model and a modification component. A fully conservative numerical algorithm, using wellbalanced slope limited centred scheme, is deployed to solve the governing equations. The present quasi single-phase model using four closure relationships for the bed shear stresses is evaluated against USGS experimental debris flow and compared with traditional quasi single-phase models and a recent physically enhanced two-phase model. It is found that the present quasi single-phase model performs much better than the traditional models, and is attractive in terms of computational cost while the two-phase model performs even better appreciably.

Volumetric lattice Boltzmann method for pore-scale mass diffusionadvection process in geopolymer porous structures

Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.

SENSITIVITY COEFFICIENTS OF SINGLE-PHASE FLOW IN LOW-PERMEABILITY HETEROGENEOUS RESERVOIRS

Theoretical equations for computing sensitivity coefficients of wellbore pressures to estimate the reservoir parameters in low-permeability reservoirs conditioning to non-Darcy flow data at low velocity were obtained. It is shown by a lot of numerical calculations that the wellbore pressures are much more sensitive to permeability very near the well than to permeability a few gridblocks away from the well. When an initial pressure gradient existent sensitivity coefficients in the region are closer to the active well than to the observation well. Sensitivity coefficients of observation well at the line between the active well and the observation well are influenced greatly by the initial pressure gradient.

Single-Phase Velocity Determination Based in Video and Sub-Images Processing:An Optical Flow Method Implemented with Support of a Programmed MatLab Structured Script

Important in many different sectors of the industry, the determination of stream velocity has become more and more important due to measurements precision necessity, in order to determine the right production rates, determine the volumetric production of undesired fluid, establish automated controls based on these measurements avoiding over-flooding or over-production, guaranteeing accurate predictive maintenance, etc. Difficulties being faced have been the determination of the velocity of specific fluids embedded in some others, for example, determining the gas bubbles stream velocity flowing throughout liquid fluid phase. Although different and already applicable methods have been researched and already implemented within the industry, a non-intrusive automated way of providing those stream velocities has its importance, and may have a huge impact in projects budget. Knowing the importance of its determination, this developed script uses a methodology of breaking-down real-time videos media into frame images, analyzing by pixel correlations possible superposition matches for further gas bubbles stream velocity estimation. In raw sense, the script bases itself in functions and procedures already available in MatLab, which can be used for image processing and treatments, allowing the methodology to be implemented. Its accuracy after the running test was of around 97% (ninety-seven percent);the raw source code with comments had almost 3000 (three thousand) characters;and the hardware placed for running the code was an Intel Core Duo 2.13 [Ghz] and 2 [Gb] RAM memory capable workstation. Even showing good results, it could be stated that just the end point correlations were actually getting to the final solution. So that, making use of self-learning functions or neural network, one could surely enhance the capability of the application to be run in real-time without getting exhaust by iterative loops.

A fractal approach to low velocity non-Darcy flow in a low permeability porous medium

In this paper, the mechanism for fluid flow at low velocity in a porous medium is analyzed based on plastic flow of oil in a reservoir and the fractal approach. The analytical expressions for flow rate and velocity of non-Newtonian fluid flow in the low permeability porous medium are derived, and the threshold pressure gradient （TPG） is also obtained. It is notable that the TPG （J） and permeability （K） of the porous medium analytically exhibit the scaling behavior J ~ K-D＇r/（l＋Or）, where DT is the fractal dimension for tortuous capillaries. The fractal characteristics of tortuosity for capillaries should be considered in analysis of non-Darcy flow in a low permeability porous medium. The model predictions of TPG show good agreement with those obtained by the available expression and experimental data. The proposed model may be conducible to a better understanding of the mechanism for nonlinear flow in the low permeability porous medium.

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.

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.

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.

Drag reduction of turbulent channel flows over an anisotropic porous wall with reduced spanwise permeability

The direct numerical simulation (DNS) is carried out for the incompressible viscous turbulent flows over an anisotropic porous wall. Effects of the anisotropic porous wall on turbulence modifications as well as on the turbulent drag reduction are investigated. The simulation is carried out at a friction Reynolds number of 180, which is based on the averaged friction velocity at the interface between the porous medium and the clear fluid domain. The depth of the porous layer ranges from 0.9 to 54 viscous units. The permeability in the spanwise direction is set to be lower than the other directions in the present simulation. The maximum drag reduction obtained is about 15.3% which occurs for a depth of 9 viscous units. The increasing of drag is addressed when the depth of the porous layer is more than 25 wall units. The thinner porous layer restricts the spanwise extension of the streamwise vortices which suppresses the bursting events near the wall. However, for the thicker porous layer, the wall-normal fluctuations are enhanced due to the weakening of the wall-blocking effect which can trigger strong turbulent structures near the wall.

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.