Analysis of one-dimensional horizontal two-phase flow in geothermal reservoirs

In the absence of capillarity the single-component two-phase porous medium equations have the structure of a nonlinear parabolic pressure (equivalently, temperature) diffusion equation, with derivative coupling to a nonlinear hyperbolic saturation wave equation. The mixed parabolic-hyperbolic system is capable of substaining saturation shock waves. The Rankine-Hugoniot equations show that the volume flux is continuous across such a shock. In this paper we focus on the horizontal one-dimensional flow of water and steam through a block of porous material within a geothermal reservoir. Starting from a state of steady flow we study the reaction of the system to simple changes in boundary conditions. Exact results are obtainable only numerically, but in some cases analytic approximations can be derived. When pressure diffusion occurs much faster than saturation convection, the numerical results can be described satisfactorily in terms of either saturation expansion fans, or isolated saturation shocks. At early times, pressure and saturation profiles are functionally related. At intermediate times, boundary effects become apparent. At late times, saturation convection dominates and eventually a steady-state is established. When both pressure diffusion and saturation convection occur on the same timescale, initial simple shock profiles evolve into multiple shocks, for which no theory is currently available. Finally, a parameter-free system of equations is obtained which satisfactorily represents a particular case of the exact equations.     

Mathematical model of two-phase fluid nonlinear flow in low-permeability porous media with applications

IntroductionItisasuccessfulexampleinadevelopmentstoryofscienceandtechnologyformechanicsoffluidsinporousmediatocombinewithengineeringtechnology .Fieldsinfluencedbythemechanicsinvolveddevelopmentofoil_gasandgroundwaterresources,controlonseawaterintrusionandsubsidenceandgeologichazards,geotechnicalengineeringandbioengineering ,andairlineindustry[1～ 7].Aproblemonnonlinearflowinlow_permeabilityporousmediaisbutonlyabasiconeindifferentkindsofengineeringfields,butalsooneoffrontlineresearchfieldsofmod…     

Effect of fluids properties on non-equilibrium capillarity effects: Dynamic pore-network modeling

We have developed a Dynamic Pore-network model for Simulating Two-phase flow in porous media (DYPOSIT). The model is applicable to both drainage and imbibition processes. Employing improved numerical and geometrical features in the model facilitate a physically-based pore-scale simulator. This computational tool is employed to perform several numerical experiments (primary and main drainage, main imbibition) to investigate the current capillarity theory. Traditional two-phase flow formulations state that the pressure difference between the two phase is equal to the capillary pressure, which is assumed to be a function of saturation only. Many theoretical and experimental studies have shown that this assumption is invalid and the pressure difference between the two fluids is not only equal to the capillary pressure but is also related to the variation of saturation with time in the domain; this is referred to as the non-equilibrium capillarity effect. To date, non-equilibrium capillarity effect has been investigated mainly under drainage. In this study, we analyze the non-equilibrium capillarity theory under drainage and imbibition as a function of saturation, viscosity ratio, and effective viscosity. Other aspects of the dynamics of two-phase flow such as trapping and saturation profile are also studied.     

Degassing of water-solved gases in thin pipings during fast pressure drops

 A one-dimensional model is presented, which describes the transient two-phase flow in thin pipes during fast pressure drops and degassing by use of Eulerian and Lagrangian systems. The reduction in dimension is obtained by introduction of a geometry model for bubbly and slug flow regimes. The complete model includes the transient two-phase flow, bubble formation and bubble growth. The flow model predicts rising velocities of bubbles and plugs in arbitrary inclined highly accurate pipes. The mass transfer (diffusion) of the dissolved phase is calculated by the bubble growth model. The quality of the model was examined by simulation of experimental series, whereby water was depressurised from the saturation pressure of the dissolved gas mixture (air), by variation of saturation pressure, pressure gradient and pipe geometry. The results of numerical simulation fit the experimental data well. Received on 17 January 2000     

Vertical two-phase flow in porous media

New concepts are introduced to describe single-component two-phase flow under gravity. The phases can flow simultaneously in opposite directions (counterflow), but information travels either up or down, depending on the sign of the wavespeedC. Wavespeed, saturation and other quantities are defined on a two-sheeted surface over the mass-energy flow plane, the sheets overlapping in the counterflow region. A saturation shock is represented as an instantaneous displacement along a line of constant volume fluxJ _Q in the flow plane. Most shocks are of the wetting type, that is, they leave the environment more saturated after their passage. When flow is horizontal all shocks are wetting, but it is a feature of vertical two-phase flow that for sufficiently small mass and energy flows there also exist drying shocks associated with lower final saturations.     

Relative Permeability Analysis of Tube Bundle Models,Including Capillary Pressure

The analytical equations for calculating two-phase flow, including local capillary pressures, are developed for the bundle of parallel capillary tubes model. The flow equations that are derived were used to calculate dynamic immiscible displacements of oil by water under the constraint of a constant overall pressure drop across the tube bundle. Expressions for averaged fluid pressure gradients and total flow rates are developed, and relative permeabilities are calculated directly from the two-phase form of Darcy's law. The effects of pressure drop and viscosity ratio on the relative permeabilities are discussed. Capillary pressure as a function of water saturation was delineated for several cases and compared to a steady-state mercury-injection drainage type of capillary pressure profile. The bundle of serial tubes model (a model containing tubes whose diameters change randomly at periodic intervals along the direction of flow), including local Young-Laplace capillary pressures, was analyzed with respect to obtaining relative permeabilities and macroscopic capillary pressures. Relative permeabilities for the bundle of parallel tubes model were seen to be significantly affected by altering the overall pressure drop and the viscosity ratio; relative permeabilities for the bundle of serial tubes were seen to be relatively insensitive to viscosity ratio and pressure, and were consistently X-like in profile. This work also considers the standard Leverett (1941) type of capillary pressure versus saturation profile, where drainage of a wetting phase is completed in a step-wise steady fashion; it was delineated for both tube bundle models. Although the expected increase in capillary pressure at low wetting-phase saturation was produced, comparison of the primary-drainage capillary pressure curves with the pseudo-capillary pressure profiles, that are computed directly using the averaged pressures during the displacements, revealed inconsistencies between the two definitions of capillary pressure.     

低渗透煤层气藏中气-水两相不稳定渗流动态分析

针对低渗透煤层渗流问题,考虑启动压力梯度及其引起的动边界和动边界内吸附气解吸作用的渗流模型研究目前仅限于单相流,而更符合实际的气-水两相渗流动边界模型未见报道.本文综合考虑了煤层吸附气的解吸作用、气-水两相渗流、非达西渗流、地层应力敏感等影响因素,进行了低渗透煤层的气-水两相渗流模型研究.采用了试井技术中的"分相处理"方法,修正了两相渗流的综合压缩系数和流度,并基于含气饱和度呈线性递减分布的假设,建立了煤层气藏的气-水两相渗流耦合模型.该数学模型不仅可以描述由于低渗透煤层中渗流存在启动压力梯度而产生的可表征煤层有效动用范围随时间变化的移动边界,还可以描述煤层有效动用范围内吸附气的解吸现象以及吸附气解吸作用所引起的煤层含气饱和度的上升;为了提高模型精度,控制方程还保留了二次压力梯度项.采用了稳定的全隐式有限差分方法进行了模型的数值求解,并验证了数值计算方法的正确性,获得了模型关于瞬时井底压力与压力导数响应的双对数特征曲线,由此分析了各渗流参数的敏感性影响.本文研究结果可为低渗透煤层气藏开发的气-水两相流试井技术提供渗流力学的理论基础.     

An Improved IMPES Method for Two-Phase Flow in Porous Media

In this paper we develop and numerically study an improved IMPES method for solving a partial differential coupled system for two-phase flow in a three-dimensional porous medium. This improved method utilizes an adaptive control strategy on the choice of a time step for saturation and takes a much larger time step for pressure than for the saturation. Through a stability analysis and a comparison with a simultaneous solution method, we show that this improved IMPES method is effective and efficient for the numerical simulation of two-phase flow and it is capable of solving two-phase coning problems.     

Geothermal heat pipe stability: Solution selection by upstreaming and boundary conditions

In a geothermal reservoir, the heat pipe mechanism can transfer heat very efficiently, with vapor rising and liquid falling in comparable quantities, driven by gravity. For a given heat and mass flux that is not too large, there are two possible steady solutions with vapor-liquid counterflow, one liquid-dominated, and one vapor-dominated. Numerical solution of the equations for two-phase vertical counterflow displays intriguing stability behaviour. If pressure and saturation are fixed at depth, and heat and mass flux specified at the top, the vapor-dominated solution is almost always obtained. That is, for a variety of boundary values, the solution settles to the vapor-dominated steady-state, and only for very special values is it possible to obtain the liquid-dominated case. Similarly the liquid-dominated solution is almost always obtained if the boundary conditions are reversed, with pressure and saturation fixed at the top and heat and mass flux specified at depth.This behaviour is here explained in two complementary ways. It is shown to be a consequence of upstream differencing of the flow terms in the numerical method. It is also shown to be expected behaviour for wavelike saturation solutions. Hence the observed behaviour is not only a direct consequence of the numerical method used, but is fundamental to geothermal heat pipes.     

Modeling and Numerical Simulations of Immiscible Compressible Two-Phase Flow in Porous Media by the Concept of Global Pressure

A new formulation is presented for the modeling of immiscible compressible two-phase flow in porous media taking into account gravity, capillary effects, and heterogeneity. The formulation is intended for the numerical simulation of multidimensional flows and is fully equivalent to the original equations, contrary to the one introduced in Chavent and Jaffré (Mathematical Models and Finite Elements for Reservoir Simulation, 1986). The main feature of this formulation is the introduction of a global pressure. The resulting equations are written in a fractional flow formulation and lead to a coupled system which consists of a nonlinear parabolic (the global pressure equation) and a nonlinear diffusion–convection one (the saturation equation) which can be efficiently solved numerically. A finite volume method is used to solve the global pressure equation and the saturation equation for the water and gas phase in the context of gas migration through engineered and geological barriers for a deep repository for radioactive waste. Numerical results for the one-dimensional problem are presented. The accuracy of the fully equivalent fractional flow model is demonstrated through comparison with the simplified model already developed in Chavent and Jaffré (Mathematical Models and Finite Elements for Reservoir Simulation, 1986).     

ON A NEW KIND OF METHODS TO SOLVE THE PLANE PROBLEMS OF TWO-PHASE FLOW THROUGH POROUS MEDIA

This paper presents a new kind of method for solving the planeproblems of two-phase flow in porous media.The ellipticalpartial differential equation for pressure distribution is sol-ved by the finite element method,and then the semi-analyticalsolution for pressure gradient is used to determine the new sa-turation field according to the existing exact formula describ-ing the saturation propagation along the streamlines.The maindistinguishing feature and advantage of this kind of method arethe ability to overcome the numerical dispersion which is inhe-rent in the ordinary numerical simulation methods.and thereby,to give a precise and clear-cut position of the saturation dis-continuity in the water-oil displacement front.moreover.thesaturation equation,which should commonly be solved simultan-eously or alternatively with the pressure equation,is complete-ly avoided,so that the computing time is greatly reduced.     

非牛顿流体两相渗流问题的摄动解

本文用奇异摄动法结合正则摄动法求解了考虑毛管力因素时多孔介质中弱非牛顿流体的两相驱替问题,得到了分流函数和湿相饱和度的渐近解析解。所得结果同数值解和经典的牛顿流体两相渗流结果进行了比较,并着重讨论了非牛顿因素的影响。     

Constant rate production of geothermal fluid from a two-phase vertical column I: Theory

Previous theoretical results on geothermal two-phase flows in porous media are extended and applied to the case of withdrawal of fluid at a constant rate from a vertical column. Dimensional considerations show that pressure and saturation behaviour is controlled by a single parameter which is the ratio of the withdrawal speed to buoyancy speed. For large flows (large) fluid withdrawal is a mining process, and a vapour dominated zone spreads out from the production level. Production enthalpies tend towards steam values. However if is small then gravity dominates, and buoyancy forces can lead to the formation of a steam bubble which escapes from the production boundary and rises towards the surface. Production enthalpy may then remain at the liquid value over long periods. In addition certain saturation ranges at the sink may be forbidden as a consequence of the constant rate boundary condition. Then saturation shocks will form at the production boundary and travel out from the sink. Internally generated shocks may also occur. Pressure and saturation response to a steady withdrawal of fluid is more complicated than in a two-phase gravity-free situation. Since gravity is an essential component of even horizontal two-phase flow this suggests that two-phase studies which ignore the role of gravity may be too simplistic.     

Two-phase flow pattern identification in a tube bundle based on void fraction and pressure measurements,with emphasis on churn flow

Two-phase flow are frequently encountered in the industry. In particular, in steam generators of nuclear plants, water is heated so that at the top of the generator an important fraction of water flows as vapor. In this upper part, a rising co-current two phase flow transverse to the tube bundle takes place. Fluids exert significant forces on the tubes in this area which highly depend on the two-phase flow pattern. Thus, as a prerequisite, it is essential to gather information on the flow conditions associated with the different two-phase flow patterns, which can be bubbly, intermittent, or annular. Then we must analyze the potentially dangerous flow patterns. This paper presents an experimental campaign aimed at characterizing those flow patterns for a rising co-current transverse flow in a tube bundle representative of the geometry in a steam generator. A new methodology based on the understanding of key contributions to vertical two-phase flow pattern maps in tube bundles is proposed that leads to a more complete flow pattern map. Finally, the paper focuses on the churn flow, which is the flow pattern for which significant pressure fluctuations occur. For this pattern, important damages could be expected on the tubes of a steam generator. Different kinds of pressure fluctuations are observed at different frequencies depending on the flow rates and the location in the test section.     

Experimental and numerical analysis of steam jet pump

Steam jet pump is the best choice for pumping radioactive and hazardous liquids because it has no moving parts and so no maintenance. However, the physics involved is highly complicated because of the mass, momentum and energy transfer between the phases involved. In this study the characteristics of SJP are studied both experimentally and numerically to pump water using saturated steam. In the experimental study the static pressure, temperature along the length of the steam jet pump and the steam and water flow rates are recorded. The three dimensional numerical study is carried out using the Eulerian two-phase flow model of Fluent 6.3 software and the direct-contact condensation model developed previously. The experimental and CFD results, of axial static pressure and temperature, match closely with each other. The mass ratio and suction lift are calculated from experimental data and it is observed that the mass ratio varies from 10 to 62 and the maximum value of suction lift is 2.12 m under the conditions of the experiment.     

Interfaces of Phase Transition and Disappearance and Method of Negative Saturation for Compositional Flow with Diffusion and Capillarity in Porous Media

The specific case of interfaces separating a single-phase fluid and a two-phase continuum appears in the theory of compositional flow through porous media. They are usually called the interfaces of phase transition (PT-interfaces) or the interfaces of phase disappearing (PD-interfaces). The principle of equivalence is proved which shows that a single-phase multi-component fluid may be replaced by an equivalent fictitious two-phase fluid having specific properties. The equivalent properties are such that the extended saturation of a fictitious phase is negative. This principle enables us to develop the uniform system of two-phase equations in the overall domain in terms of the extended saturation (the NegSat model), and to apply the direct numerical simulation. In the case with diffusion, the uniform NegSat model contains a new term proportional to the gradient of saturation in the relation for flow velocity. The canonical NegSat model represents a transport equation with discontinuous nonlinearities. The qualitative analysis of this model shows that the PT-interfaces represent the shocks of the extended saturation, or, in some cases, can transform into weak shocks. The diffusion and capillarity do not destroy necessarily the shocks, but change their velocity. The analytical technique is developed which allows capturing PT-shocks. The method is illustrated by several examples of miscible gas injection in oil reservoir. In two-dimensional case, the effects of multiple shock collisions in heterogeneous media are automatically modeled. In the case of immiscible fluids and a classic interface, the suggested method converges to the VOF-method.     

Implicit finite element methods for two-phase flow in oil reservoirs

The equations governing immiscible, incompressible, two-phase, porous media flow are discretized by generalized streamline diffusion Petrov–Galerkin methods in space and by implicit differences in time. Systems of non-linear algebraic equations are solved by Newton–Raphson iteration employing ILU-preconditioned conjugate-gradient-like methods to the non-symmetric matrix system in each iteration. The resulting solution methods are robust, enable complex grids with irregular nodal orderings and allow capillary effects. Several numerical formulations are tested and compared for one-, two- and three-dimensional flow cases, with emphasis on problems involving saturation shocks, heterogeneous media and curved boundaries. For reservoirs consisting of multiple rock types with differing capillary pressure properties, it is shown that traditional Bubnov-Galerkin methods give poor results and the new Petrov–Galerkin formulations are required. Investigations regarding the behaviour of several preconditioned conjugate-gradient-like methods in these type of problems are also reported.     

Numerical simulation of 1-D unsteady two-phase flows with shocks

In the present paper, random-choice method (RCM) and second-order GRP difference method, which are high resolution methods used for pure gas flows with shocks, are extended and employed to study the problem of one-dimensional unsteady two-phase flows. The two-phase shock wave and the flow field behind it in a dusty gas shock tube are calculated and the time-dependent change of the fiow parameters for the gas antiparticle phase are obtained. The numerical results indicate that both the two methods can give the relaxation structure of the two-phase shocks with a sharp discontinuous front and that the GRP method has the advantages of less time-consuming and higher accuracy over the RCM method.     

A Mathematical Model for Hysteretic Two-Phase Flow in Porous Media

We develop a mathematical model for hysteretic two-phase flow (of oil and water) in waterwet porous media. To account for relative permeability hysteresis, an irreversible trapping-coalescence process is described. According to this process, oil ganglia are created (during imbibition) and released (during drainage) at different rates, leading to history-dependent saturations of trapped and connected oil. As a result, the relative permeability to oil, modelled as a unique function of the connected oil saturation, is subject to saturation history. A saturation history is reflected by history parameters, that is by both the saturation state (of connected and trapped oil) at the most recent flow reversal and the most recent water saturation at which the flow was a primary drainage. Disregarding capillary diffusion, the flow is described by a hyperbolic equation with the connected oil saturation as unknown. This equation contains functional relationships which depend on the flow mode (drainage or imbibition) and the history parameters. The solution consists of continuous waves (expansion waves and constant states), shock waves (possibly connecting different modes) and stationary discontinuities (connecting different saturation histories). The entropy condition for travelling waves is generalized to include admissible shock waves which coincide with flow reversals. It turns out that saturation history generally has a strong influence on both the type and the speed of the waves from which the solution is constructed.     

Analytical Solution to the Riemann Problem of Three-Phase Flow in Porous Media

In this paper we study one-dimensional three-phase flow through porous media of immiscible, incompressible fluids. The model uses the common multiphase flow extension of Darcys equation, and does not include gravity and capillarity effects. Under these conditions, the mathematical problem reduces to a 2 × 2 system of conservation laws whose essential features are: (1) the system is strictly hyperbolic; (2) both characteristic fields are nongenuinely nonlinear, with single, connected inflection loci. These properties, which are natural extensions of the two-phase flow model, ensure that the solution is physically sensible. We present the complete analytical solution to the Riemann problem (constant initial and injected states) in detail, and describe the characteristic waves that may arise, concluding that only nine combinations of rarefactions, shocks and rarefaction-shocks are possible. We demonstrate that assuming the saturation paths of the solution are straight lines may result in inaccurate predictions for some realistic systems. Efficient algorithms for computing the exact solution are also given, making the analytical developments presented here readily applicable to interpretation of lab displacement experiments, and implementation of streamline simulators.