Dynamic capillary effects in heterogeneous porous media

In standard multi-phase flow models on porous media, a capillary pressure saturation relationship developed under static conditions is assumed. Recent experiments have shown that this static relationship cannot explain dynamic effects as seen for example in outflow experiments. In this paper, we use a static capillary pressure model and a dynamic capillary pressure model based on the concept of Hassanizadeh and Gray and examine the behavior with respect to material interfaces. We introduce a new numerical scheme for the one-dimensional case using a Lagrange multiplier approach and develop a suitable interface condition. The behavior at the interface is discussed and verified by various numerical simulations.     

Immiscible two-phase Darcy flow model accounting for vanishing and discontinuous capillary pressures: application to the flow in fractured porous media

Fully implicit time-space discretizations applied to the two-phase Darcy flow problem leads to the systems of nonlinear equations, which are traditionally solved by some variant of Newton’s method. The efficiency of the resulting algorithms heavily depends on the choice of the primary unknowns since Newton’s method is not invariant with respect to a nonlinear change of variable. In this regard, the role of capillary pressure/saturation relation is paramount because the choice of primary unknowns is restricted by its shape. We propose an elegant mathematical framework for two-phase flow in heterogeneous porous media resulting in a family of formulations, which apply to general monotone capillary pressure/saturation relations and handle the saturation jumps at rocktype interfaces. The presented approach is applied to the hybrid dimensional model of two-phase water-gas Darcy flow in fractured porous media for which the fractures are modelled as interfaces of co-dimension one. The problem is discretized using an extension of vertex approximate gradient scheme. As for the phase pressure formulation, the discrete model requires only two unknowns by degree of freedom.     

Sequential implicit vertex approximate gradient discretization of incompressible two-phase Darcy flows with discontinuous capillary pressure

This work deals with sequential implicit schemes for incompressible and immiscible two-phase Darcy flows which are commonly used and well understood in the case of spatially homogeneous capillary pressure functions. To our knowledge, the stability of this type of splitting schemes solving sequentially a pressure equation followed by the saturation equation has not been investigated so far in the case of discontinuous capillary pressure curves at different rock type interfaces. It will be shown here to raise severe stability issues for which stabilization strategies are investigated in this work. To fix ideas, the spatial discretization is based on the Vertex Approximate Gradient (VAG) scheme accounting for unstructured polyhedral meshes combined with an Hybrid Upwinding (HU) of the transport term and an upwind positive approximation of the capillary and gravity fluxes. The sequential implicit schemes are built from the total velocity formulation of the two-phase flow model and only differ in the way the conservative VAG total velocity fluxes are approximated. The stability, accuracy and computational cost of the sequential implicit schemes studied in this work are tested on oil migration test cases in 1D, 2D and 3D basins with a large range of capillary pressure parameters for the drain and barrier rock types. It will be shown that usual splitting strategies fail to capture the right solutions for highly contrasted rock types and that it can be fixed by maintaining locally the pressure saturation coupling at different rock type interfaces in the definition of the conservative total velocity fluxes. The numerical investigation of the sequential schemes is also extended to the widely used finite volume Two-Point Flux Approximation spatial discretization.

     

多孔介质两相系统毛细压力与饱和度关系试验研究

两相系统毛细压力-饱和度(h~S)关系曲线的确定是多孔介质多相流动研究的基础。采用简易试验装置对理想和实际介质中水-气和油-水两相系统中的h~S关系曲线进行了测定。试验结果表明,对于相同两相系统,多孔介质孔隙度愈小,同一毛细压力对应的饱和度相应愈大;对于不同两相系统,理想介质的关系曲线在一定毛细压力以下平缓,较大毛细压力时陡直,实际介质关系曲线走势相对较陡。分析结果表明,水-气和油-水两相系统的实测数据符合Parker等提出的基于van Genuchten(1980)关系式的折算理论;应用折算理论,可以在同一多孔介质某一两相系统h~S关系已知的情况下较好地估计同一孔隙度条件下其它两相系统的h~S关系曲线。     

Numerical modeling of two-phase binary fluid mixing using mixed finite elements

Diffusion coefficients of dense gases in liquids can be measured by considering two-phase binary nonequilibrium fluid mixing in a closed cell with a fixed volume. This process is based on convection and diffusion in each phase. Numerical simulation of the mixing often requires accurate algorithms. In this paper, we design two efficient numerical methods for simulating the mixing of two-phase binary fluids in one-dimensional, highly permeable media. Mathematical model for isothermal compositional two-phase flow in porous media is established based on Darcy’s law, material balance, local thermodynamic equilibrium for the phases, and diffusion across the phases. The time-lag and operator-splitting techniques are used to decompose each convection–diffusion equation into two steps: diffusion step and convection step. The Mixed finite element (MFE) method is used for diffusion equation because it can achieve a high-order and stable approximation of both the scalar variable and the diffusive fluxes across grid–cell interfaces. We employ the characteristic finite element method with moving mesh to track the liquid–gas interface. Based on the above schemes, we propose two methods: single-domain and two-domain methods. The main difference between two methods is that the two-domain method utilizes the assumption of sharp interface between two fluid phases, while the single-domain method allows fractional saturation level. Two-domain method treats the gas domain and the liquid domain separately. Because liquid–gas interface moves with time, the two-domain method needs work with a moving mesh. On the other hand, the single-domain method allows the use of a fixed mesh. We derive the formulas to compute the diffusive flux for MFE in both methods. The single-domain method is extended to multiple dimensions. Numerical results indicate that both methods can accurately describe the evolution of the pressure and liquid level.     

Modelling Two-Phase Incompressible Flow in Porous Media Using Mixed Hybrid and Discontinuous Finite Elements

In this paper, we present a numerical model for simulating two-phase (oil–water and air–water) incompressible and immiscible flow in porous media. The mathematical model which is based on a fractional flow formulation is formed of two nonlinear partial differential equations: a mean pressure equation and a water saturation equation. These two equations can be solved in a sequential manner. Two numerical methods are used to discretize the equations of the two-phase flow model: mixed hybrid finite elements are used to treat the pressure equation, h-based Richards' equation and the diffusion term in the saturation equation, the advection term in the saturation equation is treated with the discontinuous finite elements. We propose a better way to calculate the nonlinear coefficients contained in our equations on each element of the discretized domain. In heterogeneous porous media, the saturation becomes discontinuous at the interface between two porous media. We show in this paper how to use the capillary pressure–saturation relationship in order to handle the saturation jump in the mixed hybrid finite element method. The two-phase flow simulator is verified against analytical solutions for some flow problems treated by other authors.     

Investigation of the effect of specific interfacial area on strength of unsaturated granular materials by X-ray tomography

This paper studies the effect of interfacial areas (air–water interfaces and solid–water interfaces) on material strength of unsaturated granular materials. High-resolution X-ray computed tomography technique is employed to measure the interfacial areas in wet glass bead samples. The scanned 3D images are trinarized into three phases and meshed into representative volume elements (RVEs). An appropriate RVE size is selected to represent adequate local information. Due to the local heterogeneity of the material, the discretized RVEs of the scanned samples actually cover a very large range of degree of saturation and porosity. The data of RVEs present the relationship between the specific interfacial areas and degree of saturation and gives boundaries where the interfacial area of a whole sample should fall in. In parallel, suction-controlled direct shear tests have been carried out on glass beads and the material strength has been corroborated with two effective stress definitions related to the specific air–water interfacial areas and fraction of wetted solid surface, respectively. The comparisons show that the specific air–water interfacial area reaches the peak at about 25% of saturation and contributes significantly to the material strength (up to 60% of the total capillary strength). The wetted solid surface obtained from X-ray CT is also used to estimate Bishop’s coefficient χ based on the second type of effective stress definition, which shows a good agreement with the measured value. This work emphasizes the importance to include interface terms in effective stress formulations of unsaturated soils. It also suggests that the X-ray CT technique and RVE-based multiscale analysis are very valuable in the studies of multiphase geomaterials.

     

A linear constitutive model for unsaturated poroelasticity by micromechanical analysis

This paper extends the Biot theory of poroelasticity from the saturated to unsaturated case. The Biot phenomenological model uses parameters that are easily observable, such as the deformation of porous frame, total stress, pore pressure, and fluid specific discharge. Such model is preferred for engineering applications. At this macroscopic level, the extension of Biot theory from saturated to unsaturated is straightforward. The constitutive constants, however, are combined properties of solid, pore space, and fluids. In the unsaturated case, the constants are functions of the degree of saturation. Their measurements and tabulation over a range of saturation is generally not feasible for practical applications. In this work, a Biot-Willis type analysis is performed for the unsaturated case to provide a theory that the bulk material constants can be evaluated using laboratory measurable micromechanical constants under saturated condition, plus a capillary pressure curve (saturation versus suction pressure) typically available for unsaturated porous medium, without the need of measurement at each state of saturation. In particular, it is demonstrated that the surface energy contained in the meniscus interface manifests as a “capillary modulus,” given by the negative inverse slope of the capillary pressure curve. A rigorous analysis based on the thermodynamic variational energy approach is also conducted to lend theoretical support to the phenomenological approach. The presented model can bring a closure to the practical engineering modeling of the deformation of partially saturated porous medium that lacks the information of material constants over the range of saturation.     

Solving the Buckley–Leverett equation with gravity in a heterogeneous porous medium

Immiscible two‐phase flow in porous media can be described by the fractional flow model. If capillary forces are neglected, then the saturation equation is a non‐linear hyperbolic conservation law, known as the Buckley–Leverett equation. This equation can be numerically solved by the method of Godunov, in which the saturation is computed from the solution of Riemann problems at cell interfaces. At a discontinuity of permeability this solution has to be constructed from two flux functions. In order to determine a unique solution an entropy inequality is needed. In this article an entropy inequality is derived from a regularisation procedure, where the physical capillary pressure term is added to the Buckley‐Leverett equation. This entropy inequality determines unique solutions of Riemann problems for all initial conditions. It leads to a simple recipe for the computation of interface fluxes for the method of Godunov. This revised version was published online in July 2006 with corrections to the Cover Date.     

Accurate seabed modeling using finite difference methods

Finite difference is the most widely used method for seismic wavefield modeling. However, most finite-difference implementations discretize the Earth model over a fixed grid interval. This can lead to irregular model geometries being represented by ‘staircase’ discretization, and potentially causes mispositioning of interfaces within the media. This misrepresentation is a major disadvantage to finite difference methods, especially if there exist strong and sharp contrasts in the physical properties along an interface. The discretization of undulated seabed bathymetry is a common example of such misrepresentation of the physical properties in finite-difference grids, as the seabed is often a particularly sharp interface owing to the rapid and considerable change in material properties between fluid seawater and solid rock. There are two issues typically involved with seabed modeling using finite difference methods: firstly, the travel times of reflections from the seabed are inaccurate as a consequence of its spatial mispositioning; secondly, artificial diffractions are generated by the staircase representation of dipping seabed bathymetry. In this paper, we propose a new method that provides a solution to these two issues by positioning sharp interfaces at fractional grid locations. To achieve this, the velocity model is first sampled in a model grid that allows the center of the seabed to be positioned at grid points, before being interpolated vertically onto a regular modeling grid using the windowed sinc function. This procedure allows undulated seabed bathymetry to be represented with improved accuracy during modeling. Numerical tests demonstrate that this method generates reflections with accurate travel times and effectively suppresses artificial diffractions.     

Finite volume approximation for an immiscible two-phase flow in porous media with discontinuous capillary pressure

We consider an immiscible incompressible two-phase flow in a porous medium composed of two different rocks so that the capillary pressure field is discontinuous at the interface between the rocks. This leads us to apply a concept of multivalued phase pressures and a notion of weak solution for the flow which have been introduced in Cancès and Pierre (SIAM J Math Anal 44(2):966–992, 2012). We discretize the problem by means of a numerical algorithm which reduces to a standard finite volume scheme in each rock and prove the convergence of the approximate solution to a weak solution of the two-phase flow problem. The numerical experiments show in particular that this scheme permits to reproduce the oil-trapping phenomenon.     

考虑毛细滞回效应的非饱和多孔介质渗流与变形耦合本构模型

在Wheeler本构模型框架的基础上,提出了一个水力与力学耦合的本构模型。该模型中的土-水特征曲线采用毛细滞回内变量模型,能够更好地描述水力历史变化下毛细滞回现象对非饱和多孔介质变形的影响,同时也可描述非饱和多孔介质变形对渗流的影响。非饱和土的强度不仅与吸力有关,而且受到饱和度的影响。相同的吸力下,土样经过吸湿和脱湿路径的饱和度不同,因此,非饱和土的强度也不同。此模型以体积含水率的塑性变化和体变的塑性变化为硬化参数,不仅能描述基质吸力对非饱和土的强化作用,而且考虑了饱和度对强度及变形的影响。试验结果与模型预测基本吻合,证明该模型能够模拟非饱和土的主要特性。为了简化,此模型是在各向同性荷载下推得的,有待于推广到一般的应力状态     

Moisture Transport Through Sprayed Concrete Tunnel Linings

Waterproofing of permanent sprayed concrete tunnel linings with sprayed membranes in a continuous sandwich structure has been attempted since 2000 and has seen increased use in some countries. The main function of a sprayed membrane from a waterproofing perspective is to provide crack bridging and hence prevent flow of liquid water into the tunnel through cracks and imperfections in the concrete material. However, moisture can migrate through the concrete and EVA-based membrane materials by capillary and vapor diffusion mechanisms. These moisture transport mechanisms can have an influence on the degree of saturation, and may influence the pore pressures in the concrete material as well as risk of freeze–thaw damage of the concrete and membrane. The paper describes a detailed study of moisture transport material parameters, moisture condition in tunnel linings and climatic conditions tunnels in hard rock in Norway. These data have been included in a hygrothermal simulation model in the software WUFI for moisture transport to substantiate moisture transport and long-term effects on saturation of the concrete and membrane material. The findings suggest that EVA-based membranes exhibit significant water absorption and vapor transport properties although they are impermeable to liquid water flow. State-of-the-art sprayed concrete material applied with the wet mix method exhibits very low hydraulic conductivities, lower than 10^?14 m/s, thus saturated conductive water flow is a very unlikely dominant transport mechanism. Moisture transport through the lining structure by capillary flow and vapor diffusion are calculated to approximately 3 cm³/m² per day for lining thicknesses in the range of 25–35 cm and seasonal Nordic climate variations. The calculated moisture contents in the tunnel linings from the hygrothermal simulations are largely in agreement with the measured moisture contents in the tunnel linings. The findings also indicate that the concrete material exhibits a reduction of saturation on the immediate inside of the membrane. Near the location of the waterproofing membrane on either side, the concrete material exhibits degrees of capillary saturation between 85 and 95 %. Moisture content in the membrane is found to be consistently in the range of 12–17 % by weight, corresponding to a degree of saturation of 30–35 %. Possible effects of such moisture contents are lower risk of freezing degradation, higher tensile bonding strengths at the membrane interfaces, and a reduced risk of pore pressure in the concrete material.     

Numerical methods for flow in a porous medium with internal boundaries

Flow of fluids and transport of solutes in porous media are subjects of wide interest in several fields of applications: reservoir engineering, subsurface hydrology, chemical engineering, etc. In this paper we will study two-phase flow in a model consisting of two different types of sediments. Here, the absolute permeability, the relative permeabilities and the capillary pressure are discontinuous functions in space. This leads to interior boundary value problems at the interface between the sediments. The saturation S_w will be discontinuous or experience large gradients at the interface. A new solution procedure for such problems will be presented. The method combines the modified method of characteristics with a weak formulation where the basis functions are discontinuous at the interior boundary. The modified method of characteristics will provide a good first approximation for the jump in the discontinuous basis functions, which leads to a fast converging iterative solution scheme for the complete problem. The method has been implemented in a two-dimensional simulator, and results from numerical experiments will be presented. This revised version was published online in August 2006 with corrections to the Cover Date.     

岩石裂隙毛管压力-饱和度关系曲线的试验研究

介绍了三峡花岗岩体裂隙毛管压力-饱和度试验。试验采用互不溶混驱替法。试验结果表明,在渗流基本特征方面,裂隙非饱和渗流毛管压力-饱和度关系曲线与空隙介质水分特征曲线具有相似性,如毛管压力-饱和度关系曲线的滞后现象;湿润流体(水)的排泄曲线具有进气压和束缚水饱和度;非湿润流体的吸湿曲线具有残余饱和度。这种相似性表明,孔隙介质非饱和渗流的研究成果可用于裂隙非饱和渗流,孔隙介质水分特征曲线的解析模型,可用于研究裂隙毛管压力-饱和度关系曲线和拟合毛管压力-饱和度排泄曲线的试验数据。     

Variational inequalities for modeling flow in heterogeneous porous media with entry pressure

One of the driving forces in porous media flow is the capillary pressure. In standard models, it is given depending on the saturation. However, recent experiments have shown disagreement between measurements and numerical solutions using such simple models. Hence, we consider in this paper two extensions to standard capillary pressure relationships. Firstly, to correct the nonphysical behavior, we use a recently established saturation-dependent retardation term. Secondly, in the case of heterogeneous porous media, we apply a model with a capillary threshold pressure that controls the penetration process. Mathematically, we rewrite this model as inequality constraint at the interfaces, which allows discontinuities in the saturation and pressure. For the standard model, often finite-volume schemes resulting in a nonlinear system for the saturation are applied. To handle the enhanced model at the interfaces correctly, we apply a mortar discretization method on nonmatching meshes. Introducing the flux as a new variable allows us to solve the inequality constraint efficiently. This method can be applied to both the standard and the enhanced capillary model. As nonlinear solver, we use an active set strategy combined with a Newton method. Several numerical examples demonstrate the efficiency and flexibility of the new algorithm in 2D and 3D and show the influence of the retardation term. This work was supported in part by IRTG NUPUS.     

非离子表面活性剂对两相系统毛细压力和饱和度关系的影响

表面活性剂至水相的引入将对原先包含纯水两相系统的毛细压力和饱和度关系存在影响.对含有非离子表面活性剂Triton X-100的水-气和水-油两相系统中的毛细压力-饱和度关系曲线进行了试验测定.试验结果表明,同不含表面活性剂的纯水系列相比,在同一饱和度情况下,含有表面活性剂的Triton X-100(0.1%)系列对应的毛细压力水头值都有不同程度的减小,说明在Triton X-100存在的情况下,驱替出同样数量的湿润相体积所需的毛细压力值较小.以van Genuchten关系式为基础的拟合结果表明,在已知纯水系列毛细压力饱和度关系的情况下,对于Triton X-100-气系统的毛细压力饱和度关系,考虑界面张力降低作用引入折算系数得到的拟合值更接近于真实值;而对于Triton X-100-油系统拟合值接近真实值的程度则随多孔介质的不同而有所不同.     

考虑毛细管压力的纵波速度衰减特征分析——以球状斑块饱和模型为例

地震波在斑块饱和岩石中传播往往会引起地震波速度的频散与衰减,不同的斑块大小及分布引起速度的频散与衰减不同.毛细管压力作为影响斑块物理特征的主要因素之一,其对速度的频散与衰减的影响知之甚少.为了研究毛细管压力的影响,笔者利用斑块膜刚度来表示毛细管压力的宏观响应,通过改变球状斑块饱和模型的边界条件将毛细管压力考虑到球状斑块饱和模型中,得到了一种新的改进球状斑块饱和模型,并基于该模型对比了考虑毛细管压力前后的地震频段的速度频散与衰减变化.数值模拟结果表明,在地震频段内,与原始球状斑块饱和模型相比,新模型的速度较原始模型速度大,频散降低,衰减减小.除此之外,利用新模型解释了已发表的不同饱和度情况下速度和衰减系数的实验室测试结果.与原始球状斑块饱和模型相比,新模型能够更好的解释不同饱和度下速度与衰减系数的变化趋势,对于斑块饱和岩石的速度饱和度,衰减系数饱和度关系解释具有重要的指示意义.     

Node-centered finite volume discretizations for the numerical simulation of multiphase flow in heterogeneous porous media

Three node-centered finite volume discretizations for multiphase porous media flow are presented and compared. By combination of these methods two additional discretization methods are generated. The ability of these schemes to describe flows at textural interfaces of different geologic formations is investigated. It was found that models with nonzero-entry pressures for the capillary pressure-saturation relationship in conjunction with the Box discretization may give rise to spurious oscillations for flows around low permeable lenses. Furthermore, the applicability and sensitivity of the discretization methods with regard to the used computational grids is discussed. The schemes are used for the numerical study of two-phase flow in porous media with zones of different material properties. This revised version was published online in July 2006 with corrections to the Cover Date.