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 ef...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.展开更多
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 physicsi...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.展开更多
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.T...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.展开更多
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 o...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.展开更多
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 o...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.展开更多
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...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.展开更多
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 ...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.展开更多
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 consum...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.展开更多
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 permeabili...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.展开更多
In this paper, we have discussed the linear stabil- ity analysis of the electrified surface separating two coaxial Oldroyd-B fluid layers confined between two impermeable rigid cylinders in the presence of both interf...In this paper, we have discussed the linear stabil- ity analysis of the electrified surface separating two coaxial Oldroyd-B fluid layers confined between two impermeable rigid cylinders in the presence of both interfacial insoluble surfactant and surface charge through porous media. The case of long waves interfacial stability has been studied. The dispersion relation is solved numerically and hence the ef- fects of various parameters are illustrated graphically. Our results reveal that the influence of the physicochemical pa- rameterβ is to shrink the instability region of the surface and reduce the growth rate of the unstable normal modes. Such important effects of the surfactant on the shape of in- terfacial structures are more sensitive to the variation of the βcorresponding to non-Newtonian fluids-model compared with the Newtonian fluids model. In the case of long wave limit, it is demonstrated that increasing r, has a dual role in- fluence (de-stabilizing effects) depending on the viscosity of the core fluid. It has a destabilizing effect at the large values of the core fluid viscosity coefficient, while this role is ex- changed to a regularly stabilizing influence at small values of such coefficient.展开更多
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 preferent...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.展开更多
We apply a proper orthogonal decomposition(POD)to data stemming from numerical simulations of a fingering instability in a multiphase flow passing through obstacles in a porous medium,to study water injection processe...We apply a proper orthogonal decomposition(POD)to data stemming from numerical simulations of a fingering instability in a multiphase flow passing through obstacles in a porous medium,to study water injection processes in the production of hydrocarbon reservoirs.We show that the time evolution of a properly defined flow correlation length can be used to identify the onset of the fingering instability.Computation of characteristic lengths for each of the modes resulting from the POD provides further information on the dynamics of the system.Finally,using numerical simulations with different viscosity ratios,we show that the convergence of the POD depends non-trivially on whether the fingering instability develops or not.This result has implications on proposed methods to decrease the dimensionality of the problem by deriving reduced dynamical systems after truncating the system’s governing equations to a few POD modes.展开更多
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 po...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.展开更多
This paper presents the analytical solutions in Laplace domain for two-dimensional nonsteady flow of slightly compressible liquid in porous media with double porosity by using the methods of integral transforms and va...This paper presents the analytical solutions in Laplace domain for two-dimensional nonsteady flow of slightly compressible liquid in porous media with double porosity by using the methods of integral transforms and variables separation. The effects of the ratio of storativities to , interporosity flow parameter on the pressure behaviors for a vertically fractured well with infinite conductivity are investigated by using the method of numerical inversion. The new log-log diagnosis graph of the pressures is given and analysed.展开更多
A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function w...A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function with three parameters, a mass conservation law and a concept of turbulent ellipses. A solution to the model was obtained by using a finite difference method and an extrapolation method. Formulas of calculating development index not only before but also after water breaks through an oil well in the condition of two-phase fluid nonlinear flow in the media were derived. An example was discussed. Water saturation distribution was presented. The moving law of drainage front was found. Laws of change of pressure difference with time were recognized. Results show that there is much difference of water saturation distribution between nonlinear flow and linear flow; that drainage front by water moves faster, water breaks through sooner and the index gets worse because of the nonlinear flow; and that dimensionless pressure difference gets larger at the same dimensionless time and difficulty of oil development becomes bigger by the nonlinear flow. Thus, it is necessary that influence of nonlinear flow on development indexes of the oil fields be taken into account. The results provide water-flooding development of the oilfields with scientific basis.展开更多
Based on the mathematical model of one dimension transient flow of the polymer foam in porous media, the numerical calculation method of the flow mentioned above by using the finite difference method is given. Through...Based on the mathematical model of one dimension transient flow of the polymer foam in porous media, the numerical calculation method of the flow mentioned above by using the finite difference method is given. Through the experiments of one dimension transient flow of HPAM (Hydrolytic Polyacrylamide) foam in the artificial sandstone core, the HPAM foam generation and coalescence coefficient of the mathematical model mentioned above are determined. The profiles of the liquid phase saturation, the pressure drop and the number density of one dimension transient flow of HPAM foam with the dimensionless time in artificial sandstone core are numerically calculated and analyzed by using the numerical calculation method.展开更多
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...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.展开更多
Microbial fuel cell (MFC) is a novel environmental friendly energy device which has received great attention due to its technology for producing electricity directly fi-om organic or inorganic matter by using bacter...Microbial fuel cell (MFC) is a novel environmental friendly energy device which has received great attention due to its technology for producing electricity directly fi-om organic or inorganic matter by using bacteria as catalyst. To date, many experiments have been carried out to achieve the maximum power output with advective flow through porous anode to the cathode in the MFC. However, the precise mechanical mechanism of flow through anode and the quantified relationship between electrode spacing and MFC performance are not yet clearly understood. It hasbeen found experimentally that the power output can be increased apparently at certain electrode spacing configuration. Based on these available experimental data, this paper investigates the effect of spacing between electrodes, the Darcy number of porous anode and the Reynolds number on the power production performance of MFC by using lattice Boltzmann method. The numerical simulation results present that the distance between electrodes significantly influences the flow velocity and residence time of the organic matter attached to the anode in the MFC. Moreover, it is found that the Darey number of porous anode and the Reynolds number can regulate the output efficiency of MFC. These results perform better understanding of the complex phenomena of MFC and will be helpful to optimize MFC design.展开更多
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...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.展开更多
基金supported by National Natural Science Foundation of China(Grant No.42172159)Science Foundation of China University of Petroleum,Beijing(Grant No.2462023XKBH002).
文摘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.
基金the National Natural Science Foundation of China(No.52274048)Beijing Natural Science Foundation(No.3222037)+1 种基金the CNPC 14th Five-Year Perspective Fundamental Research Project(No.2021DJ2104)the Science Foundation of China University of Petroleum,Beijing(No.2462021YXZZ010).
文摘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.
基金financially supported by National Natural Science Foundation of China(No.U20B6003).
文摘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.
文摘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.
文摘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.
文摘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.
基金Supported by the National Natural Science Foundation of China(52274048)Beijing Natural Science Foundation Project of China(3222037)Shaanxi Provincial Technical Innovation Project of China(2023-YD-CGZH-02).
文摘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.
文摘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.
基金supported by the National Natural Science Foundation of China(11102237)Program for Changjiang Scholars and Innovative Research Team in University(IRT1294)+1 种基金Specialized Research Fund for the Doctoral Program of Higher Education(20110133120012)China Scholarship Council(CSC)
文摘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.
文摘In this paper, we have discussed the linear stabil- ity analysis of the electrified surface separating two coaxial Oldroyd-B fluid layers confined between two impermeable rigid cylinders in the presence of both interfacial insoluble surfactant and surface charge through porous media. The case of long waves interfacial stability has been studied. The dispersion relation is solved numerically and hence the ef- fects of various parameters are illustrated graphically. Our results reveal that the influence of the physicochemical pa- rameterβ is to shrink the instability region of the surface and reduce the growth rate of the unstable normal modes. Such important effects of the surfactant on the shape of in- terfacial structures are more sensitive to the variation of the βcorresponding to non-Newtonian fluids-model compared with the Newtonian fluids model. In the case of long wave limit, it is demonstrated that increasing r, has a dual role in- fluence (de-stabilizing effects) depending on the viscosity of the core fluid. It has a destabilizing effect at the large values of the core fluid viscosity coefficient, while this role is ex- changed to a regularly stabilizing influence at small values of such coefficient.
基金supported by the Basic Research on Drilling & Completion of Critical Wells for Oil & Gas (Grant No. 51221003)National Science Fund for Petrochemical Industry (Project No. U1262201)+2 种基金"863" National Project (Project No. 2013AA064803)National Science Fund for Distinguished Young Scholars (Project No. 50925414)National Natural Science Foundation (Project No. 51074173)
文摘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.
基金support from YPF-Tecnología(YTEC)support from PICT Grant No.2015-3530.
文摘We apply a proper orthogonal decomposition(POD)to data stemming from numerical simulations of a fingering instability in a multiphase flow passing through obstacles in a porous medium,to study water injection processes in the production of hydrocarbon reservoirs.We show that the time evolution of a properly defined flow correlation length can be used to identify the onset of the fingering instability.Computation of characteristic lengths for each of the modes resulting from the POD provides further information on the dynamics of the system.Finally,using numerical simulations with different viscosity ratios,we show that the convergence of the POD depends non-trivially on whether the fingering instability develops or not.This result has implications on proposed methods to decrease the dimensionality of the problem by deriving reduced dynamical systems after truncating the system’s governing equations to a few POD modes.
基金financially supported by the Major State Basic Research Development Program of China(973 Program,Grant No.2011CB707304)the National Natural Science Foundation of China(Grant Nos.51006016,51006017,51106018 and 51106019)
文摘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.
文摘This paper presents the analytical solutions in Laplace domain for two-dimensional nonsteady flow of slightly compressible liquid in porous media with double porosity by using the methods of integral transforms and variables separation. The effects of the ratio of storativities to , interporosity flow parameter on the pressure behaviors for a vertically fractured well with infinite conductivity are investigated by using the method of numerical inversion. The new log-log diagnosis graph of the pressures is given and analysed.
文摘A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function with three parameters, a mass conservation law and a concept of turbulent ellipses. A solution to the model was obtained by using a finite difference method and an extrapolation method. Formulas of calculating development index not only before but also after water breaks through an oil well in the condition of two-phase fluid nonlinear flow in the media were derived. An example was discussed. Water saturation distribution was presented. The moving law of drainage front was found. Laws of change of pressure difference with time were recognized. Results show that there is much difference of water saturation distribution between nonlinear flow and linear flow; that drainage front by water moves faster, water breaks through sooner and the index gets worse because of the nonlinear flow; and that dimensionless pressure difference gets larger at the same dimensionless time and difficulty of oil development becomes bigger by the nonlinear flow. Thus, it is necessary that influence of nonlinear flow on development indexes of the oil fields be taken into account. The results provide water-flooding development of the oilfields with scientific basis.
文摘Based on the mathematical model of one dimension transient flow of the polymer foam in porous media, the numerical calculation method of the flow mentioned above by using the finite difference method is given. Through the experiments of one dimension transient flow of HPAM (Hydrolytic Polyacrylamide) foam in the artificial sandstone core, the HPAM foam generation and coalescence coefficient of the mathematical model mentioned above are determined. The profiles of the liquid phase saturation, the pressure drop and the number density of one dimension transient flow of HPAM foam with the dimensionless time in artificial sandstone core are numerically calculated and analyzed by using the numerical calculation method.
文摘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.
基金supported by the National Natural Science Foundation of China(Grant No51079043)the Special Fund for Public Welfare Industry of Ministry of Water Resources of China(Grants No200901064 and 201001020)the Research Innovation Program for College Graduates of Jiangsu Province(Grant No CXZZ11_0450)
基金supported by the National Natural Science Foundation of China (10932010,10972208 and11072220)
文摘Microbial fuel cell (MFC) is a novel environmental friendly energy device which has received great attention due to its technology for producing electricity directly fi-om organic or inorganic matter by using bacteria as catalyst. To date, many experiments have been carried out to achieve the maximum power output with advective flow through porous anode to the cathode in the MFC. However, the precise mechanical mechanism of flow through anode and the quantified relationship between electrode spacing and MFC performance are not yet clearly understood. It hasbeen found experimentally that the power output can be increased apparently at certain electrode spacing configuration. Based on these available experimental data, this paper investigates the effect of spacing between electrodes, the Darcy number of porous anode and the Reynolds number on the power production performance of MFC by using lattice Boltzmann method. The numerical simulation results present that the distance between electrodes significantly influences the flow velocity and residence time of the organic matter attached to the anode in the MFC. Moreover, it is found that the Darey number of porous anode and the Reynolds number can regulate the output efficiency of MFC. These results perform better understanding of the complex phenomena of MFC and will be helpful to optimize MFC design.
文摘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.