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 ne...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.展开更多
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.展开更多
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 e...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.展开更多
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 ...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.展开更多
This paper applies the convolutional differentiator method, based on generalized Forsyte orthogonal polynomial （CFPD）, to simulate the seismic wave propagation in two-phase media. From the numerical results we can s...This paper applies the convolutional differentiator method, based on generalized Forsyte orthogonal polynomial （CFPD）, to simulate the seismic wave propagation in two-phase media. From the numerical results we can see that three types of waves, fast P-waves, S-waves and slow P-waves, can be observed in the seismic wave field. The experiments on anisotropic models demonstrate that the wavefront is elliptic instead of circular and S-wave splitting occurs in anisotropic two-phase media. The research has confirmed that the rules of elastic wave propagation in fluid-saturated porous media are controlled by Biot＇s theory. Experiment on a layered fault model shows the wavefield generated by the interface and the fault very well, indicating the effectiveness of CFPD method on the wavefield modeling for real layered media in the Earth. This research has potential applications to the investigation of Earth＇s deep structure and oil/gas exploration.展开更多
Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical sol...Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.展开更多
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 remov...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.展开更多
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.展开更多
The upwind scheme is very important in the numerical approximation of some problems such as the convection dominated problem, the two-phase flow problem, and so on. For the fractional flow formulation of the two-phase...The upwind scheme is very important in the numerical approximation of some problems such as the convection dominated problem, the two-phase flow problem, and so on. For the fractional flow formulation of the two-phase flow problem, the Penalty Discontinuous Galerkin (PDG) methods combined with the upwind scheme are usually used to solve the phase pressure equation. In this case, unless the upwind scheme is taken into consideration in the velocity reconstruction, the local mass balance cannot hold exactly. In this paper, we present a scheme of velocity reconstruction in some H(div) spaces with considering the upwind scheme totally. Furthermore, the different ways to calculate the nonlinear coefficients may have distinct and significant effects, which have been investigated by some authors. We propose a new algorithm to obtain a more effective and stable approximation of the coefficients under the consideration of the upwind scheme.展开更多
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.展开更多
二氧化碳(CO_(2))捕集与封存技术有利于减少CO_(2)的排放量,近年来针对CO_(2)地质封存形成了从纳米尺度到油气藏尺度的大量研究成果,大多数研究只针对单一维度多孔介质中流动行为开展研究,且物理实验方法受许多不确定性因素影响,十分耗...二氧化碳(CO_(2))捕集与封存技术有利于减少CO_(2)的排放量,近年来针对CO_(2)地质封存形成了从纳米尺度到油气藏尺度的大量研究成果,大多数研究只针对单一维度多孔介质中流动行为开展研究,且物理实验方法受许多不确定性因素影响,十分耗费时间和成本。为了从微观角度深入理解CO_(2)地质封存过程中的渗流行为,提高CO_(2)地质埋存量,基于追踪两相界面动态变化的VOF(Volume of Fluid)方法,分别建立了2D和3D模型,开展了超临界CO_(2)-水两相流动数值模拟研究,对比了不同润湿性、毛细管数、黏度比条件下的CO_(2)团簇分布特征、CO_(2)饱和度变化规律,揭示了孔隙尺度CO_(2)埋存的内在机理。研究结果表明:①随着岩石对CO_(2)润湿性增加,CO_(2)波及范围扩大,同时CO_(2)团簇的卡断频率减少,CO_(2)埋存量增加;②随着毛细管数的增加,驱替模式由毛细指进转变为稳定驱替,CO_(2)埋存量增加;③随着注入超临界CO_(2)黏度逐渐接近水的黏度,两相流体之间的流动阻力降低,促进了“润滑效应”,CO_(2)相的渗流能力提高,CO_(2)埋存量增加;④润湿性、毛细管数、黏度比在不同维度多孔介质模型中对CO_(2)饱和度的影响程度不同。结论认为,基于VOF方法的CO_(2)-水两相渗流模拟研究在孔隙尺度上揭示了CO_(2)地质封存过程中的渗流机理,对CCUS技术的发展有指导意义,也为更大尺度的CO_(2)地质封存研究提供了理论指导和技术支撑。展开更多
基金Project“973",a national fundamental research development program
文摘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.
文摘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.
基金Iranian Offshore Oil Company (IOOC) for financial support of this work
文摘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.
基金supported by the China State Major Key Project for Basic Researches
文摘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.
基金supported by the National Natural Science Foundation of China(Grant No.40874045)Special Funds for Sciences and Technology Research of Public Welfare Trades(Grant Nos. 200811021 and 201011042)
文摘This paper applies the convolutional differentiator method, based on generalized Forsyte orthogonal polynomial （CFPD）, to simulate the seismic wave propagation in two-phase media. From the numerical results we can see that three types of waves, fast P-waves, S-waves and slow P-waves, can be observed in the seismic wave field. The experiments on anisotropic models demonstrate that the wavefront is elliptic instead of circular and S-wave splitting occurs in anisotropic two-phase media. The research has confirmed that the rules of elastic wave propagation in fluid-saturated porous media are controlled by Biot＇s theory. Experiment on a layered fault model shows the wavefield generated by the interface and the fault very well, indicating the effectiveness of CFPD method on the wavefield modeling for real layered media in the Earth. This research has potential applications to the investigation of Earth＇s deep structure and oil/gas exploration.
文摘Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.
基金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)
文摘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.
文摘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.
文摘The upwind scheme is very important in the numerical approximation of some problems such as the convection dominated problem, the two-phase flow problem, and so on. For the fractional flow formulation of the two-phase flow problem, the Penalty Discontinuous Galerkin (PDG) methods combined with the upwind scheme are usually used to solve the phase pressure equation. In this case, unless the upwind scheme is taken into consideration in the velocity reconstruction, the local mass balance cannot hold exactly. In this paper, we present a scheme of velocity reconstruction in some H(div) spaces with considering the upwind scheme totally. Furthermore, the different ways to calculate the nonlinear coefficients may have distinct and significant effects, which have been investigated by some authors. We propose a new algorithm to obtain a more effective and stable approximation of the coefficients under the consideration of the upwind scheme.
基金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.
文摘二氧化碳(CO_(2))捕集与封存技术有利于减少CO_(2)的排放量,近年来针对CO_(2)地质封存形成了从纳米尺度到油气藏尺度的大量研究成果,大多数研究只针对单一维度多孔介质中流动行为开展研究,且物理实验方法受许多不确定性因素影响,十分耗费时间和成本。为了从微观角度深入理解CO_(2)地质封存过程中的渗流行为,提高CO_(2)地质埋存量,基于追踪两相界面动态变化的VOF(Volume of Fluid)方法,分别建立了2D和3D模型,开展了超临界CO_(2)-水两相流动数值模拟研究,对比了不同润湿性、毛细管数、黏度比条件下的CO_(2)团簇分布特征、CO_(2)饱和度变化规律,揭示了孔隙尺度CO_(2)埋存的内在机理。研究结果表明:①随着岩石对CO_(2)润湿性增加,CO_(2)波及范围扩大,同时CO_(2)团簇的卡断频率减少,CO_(2)埋存量增加;②随着毛细管数的增加,驱替模式由毛细指进转变为稳定驱替,CO_(2)埋存量增加;③随着注入超临界CO_(2)黏度逐渐接近水的黏度,两相流体之间的流动阻力降低,促进了“润滑效应”,CO_(2)相的渗流能力提高,CO_(2)埋存量增加;④润湿性、毛细管数、黏度比在不同维度多孔介质模型中对CO_(2)饱和度的影响程度不同。结论认为,基于VOF方法的CO_(2)-水两相渗流模拟研究在孔隙尺度上揭示了CO_(2)地质封存过程中的渗流机理,对CCUS技术的发展有指导意义,也为更大尺度的CO_(2)地质封存研究提供了理论指导和技术支撑。
