A fracture mapping and extended finite element scheme for coupled deformation and fluid flow in fractured porous media

This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd.     

A mixed solution for two-dimensional unsteady flow in fractured porous media

A mixed finite element–boundary element solution for the analysis of two-dimensional flow in porous media composed of rock blocks and discrete fractures is described. The rock blocks are modelled implicitly by using boundary elements whereas finite elements are adopted to model the discrete fractures. The computational procedure has been implemented in a hybrid code which has been validated first by comparing the numerical results with the closed-form solution for flow in a porous aquifer intercepted by a vertical fracture only. Then, a more complex problem has been solved where a pervious, homogeneous and isotropic matrix containing a net of fractures is considered. The results obtained are shown to describe satisfactorily the main features of the flow problem under study. © 1997 by John Wiley & Sons, Ltd.     

A finite element method for modeling coupled flow and deformation in porous fractured media

Modeling the flow in highly fractured porous media by finite element method (FEM) has met two difficulties: mesh generation for fractured domains and a rigorous formulation of the flow problem accounting for fracture/matrix, fracture/fracture, and fracture/boundary fluid mass exchanges. Based on the recent theoretical progress for mass balance conditions in multifractured porous bodies, the governing equations for coupled flow and deformation in these bodies are first established in this paper. A weak formulation for this problem is then established allowing to build a FEM. Taking benefit from recent development of mesh‐generating tools for fractured media, this weak formulation has been implemented in a numerical code and applied to some typical problems of hydromechanical coupling in fractured porous media. It is shown that in this way, the FEM that has proved its efficiency to model hydromechanical phenomena in porous media is extended with all its performances (calculation time, couplings, and nonlinearities) to fractured porous media. Copyright © 2015 John Wiley & Sons, Ltd.     

Theoretical and numerical study of the steady‐state flow through finite fractured porous media

This paper investigates the two‐dimensional flow problem through an anisotropic porous medium containing several intersecting curved fractures. First, the governing equations of steady‐state fluid flow in a fractured porous body are summarized. The flow follows Darcy's law in matrix and Poiseuille's law in fractures. An infinite transversal permeability is considered for the fractures. A multi‐region boundary element method is used to derive a general pressure solution as a function of discharge through the fractures and the pressure and the normal flux on the domain boundary. The obtained solution fully accounts for the interaction and the intersection between fractures. A numerical procedure based on collocation method is presented to compute the unknowns on the boundaries and on the fractures. The numerical solution is validated by comparing with finite element solution or the results obtained for an infinite matrix. Pressure fields in the matrix are illustrated for domains containing several interconnected fractures, and mass balance at the intersection points is also checked. Copyright © 2013 John Wiley & Sons, Ltd.     

Hybrid time integration and coupled solution methods for nonlinear finite element analysis of partially saturated deformable porous media at small strain

The goal of the paper is to determine the most efficient, yet accurate and stable, finite element nonlinear solution method for analysis of partially saturated deformable porous media at small strain. This involves a comparison between fully implicit, semi‐implicit, and explicit time integration schemes, with monolithically coupled and staggered‐coupled nonlinear solution methods and the hybrid combination thereof. The pore air pressure p_a is assumed atmospheric, that is, p_a=0 at reference pressure. The solid skeleton is assumed to be pressure‐sensitive nonlinear isotropic elastic. Coupled partially saturated ‘consolidation’ in the presence of surface infiltration and traction is simulated for a simple one‐dimensional uniaxial strain example and a more complicated plane strain slope example with gravity loading. Three mixed plane strain quadrilateral elements are considered: (i) Q4P4; (ii) stabilized Q4P4S; and (iii) Q9P4; “Q” refers to the number of solid skeleton displacement nodes, and “P” refers to the number of pore fluid pressure nodes. The verification of the implementation against an analytical solution for partially saturated pore water flow (no solid skeleton deformation) and comparison between the three time integration schemes (fully implicit, semi‐implicit, and explicit) are presented. It is observed that one of the staggered‐coupled semi‐implicit schemes (SIS(b)), combined with the fully implicit monolithically coupled scheme to resolve sharp transients, is the most efficient computationally. Copyright © 2015 John Wiley & Sons, Ltd.     

采用复合单元法模拟裂隙多孔介质变饱和流动

采用复合单元法建立了模拟裂隙多孔介质变饱和流动的数值模型。该模型具有以下特点:裂隙不需要离散成特定单元,而是根据几何位置插入到孔隙基质单元中形成复合单元;在复合单元中,分别建立裂隙流和孔隙基质流的计算方程,二者通过裂隙-基质界面产生联系并整合成复合单元方程;复合单元方程具有和常规有限单元方程相同的格式,因此,可以使用常规有限单元方程的求解技术。采用欠松弛迭代、集中质量矩阵以及自适应时步调节等技术,开发了裂隙多孔介质变饱和流动计算程序。通过模拟一维干土入渗和复杂裂隙含水层内的流动问题,验证了该模型的合理性和适用性。模拟结果为进一步认识非饱和裂隙含水层地下水流动特性提供了理论依据。     

Transverse and longitudinal fluid flow modelling in fractured porous media with non‐matching meshes

A new discrete fracture model is introduced to simulate the steady‐state fluid flow in discontinuous porous media. The formulation uses a multi‐layered approach to capture the effect of both longitudinal and transverse permeability of the discontinuities in the pressure distribution. The formulation allows the independent discretisation of mesh and discontinuities, which do not need to conform. Given that the formulation is developed at the element level, no additional degrees of freedom or special integration procedures are required for coupling the non‐conforming meshes. The proposed model is shown to be reliable regardless of the permeability of the discontinuity being higher or lower than the surrounding domain. Four numerical examples of increasing complexity are solved to demonstrate the efficiency and accuracy of the new technique when compared with results available in the literature. Results show that the proposed method can simulate the fluid pressure distribution in fractured porous media. Furthermore, a sensitivity analysis demonstrated the stability regarding the condition number for wide range values of the coupling parameter.     

A model for moisture and heat flow in fractured unsaturated porous media

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 John Wiley & Sons, Ltd.     

Embedded discontinuity finite element modeling of fluid flow in fractured porous media

In many geomaterials, particularly rocks and clays, permeability is greatly enhanced by the presence of fractures. Fracture sets create an overall permeability that is anisotropic, enhanced in the directions of the fractures. In modeling the fractures via a finite element method, for example, meshing around these fractures can become quite difficult and result in computationally intensive systems. In this article, we develop a relatively simple method for including the fractures within the elements. Flow through the bulk medium is assumed to be governed by Darcy’s law, and the flow on the fracture by a generalization of the law. This model is embedded in a finite element framework, with the fractures passing through the elements. In this formulation, elements with fractures are given an enhanced permeability in the direction of the fractures. With these enhancements, the material essentially becomes anisotropically more permeable in the direction of fracture sets.     

储层流固耦合的数学模型和非线性有限元方程

根据饱和多孔介质固体骨架的平衡方程和多孔介质中流体的连续性方程,建立了储层流固耦合数学模型。模型中引入了Jaumann应力速率公式描述多孔介质固体骨架的大变形效应,并考虑了地应力、初始孔隙压力、初始流体密度和初始孔隙度对耦合模型的影响。基于与微分方程等价的加权余量公式,在空间域采用有限元离散,对时间域进行隐式差分格式离散,导出了以单元节点位移和单元节点孔隙压力为未知量的储层流固耦合的非线性有限元增量方程。该模型在石油工程中有广泛的应用,为储层流固耦合的数值模拟奠定了理论基础。     

裂隙岩体流固耦合双重介质模型的有限元计算

基于孔隙-裂隙岩体的双重孔隙介质流固耦合计算的微分方程,利用伽辽金有限元法提出的相应有限元公式,并基于岩体分类指标（RQD,RMR）提出了与岩体应力状态相关的渗透系数计算公式。编制了相应的有限元程序并给出了应用算例,将计算结果与相关文献作了比较,得出相关结论。     

模拟裂隙多孔介质中变饱和渗流的广义等效连续体方法

描述了一种计算裂隙多孔介质中变饱和渗流的广义等效连续体方法。这种方法忽略裂隙的毛细作用,设定一个与某孔隙饱和度相对应的综合饱和度极限值,并假定：（1）如果裂隙多孔介质的综合饱和度小于该极限值,水只在孔隙中存在并流动,而裂隙中则没有水的流动;（2）如果综合饱和度等于或大于该极限值,水将进入裂隙,并在裂隙内运动。分析比较了等效连续体模型的不同计算方法,并给出了一个模拟裂隙岩体中变饱和渗流与传热耦合问题的应用算例。结果表明,所述方法具有一般性,可以有效地模拟裂隙多孔介质中变饱和渗流的基本特征。     

A finite element porothermoelastic model for dual-porosity media

An existing dual-porosity finite element model has been extended to include thermo-hydro-mechanical coupling in both media. The model relies on overlapping distinct continua for the fluid and solid domains. In addition, conductive and convective heat transfers are incorporated using a single representative thermodynamics continuum. The model is applied to the problem of an inclined borehole drilled in a fractured formation subjected to a three-dimensional state of stress and, a temperature gradient between the drilling fluid and the formation. A sensitivity analysis has been carried out to study the impact of thermal loading, effect of heat transport by pore fluid flow and, the effect of parameters of the secondary medium used to represent the fractures. Copyright © 2004 John Wiley & Sons, Ltd.     

A finite element method for localized erosion in porous media with applications to backward piping in levees

This paper presents a computational method able to effectively model both the simultaneous processes typically observed in backward erosion piping, ie, the pipe tip propagation and the conduit cross-section enlargement. The numerical method is based on the novel formulation of a problem of localized erosion along a line propagating in a multidimensional porous medium. In this line, a conduit with evolving transverse size is embedded, which conveys a multiphase flow. The two systems, porous medium and pipe, are bridged by exchange terms of multiphase fluid mass and by a shared fluid pressure field. On the contrary, different fields are considered to describe flows, which are assumed as Darcian in the porous medium and turbulent in the conduit. These two flows drive pipe propagation and enlargement, respectively, as modeled by means of proper erosion kinetic laws. The corresponding numerical formulation is based on the combination between one- and multidimensional finite elements, to model the erosion conduit and the porous medium, respectively. Several simulations are proposed to demonstrate the ability of the proposed approach in reproducing available experimental data of real-scale tests on levees. Our results point out the crucial role played by the combined influence of pipe propagation and enlargement, as well as of three-dimensional (3D) effects. We also assess the mesh independence of the proposed numerical solution, particularly as concerns the calculated pipe propagation history.     

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.     

Hydro‐mechanical modeling of cohesive crack propagation in multiphase porous media using the extended finite element method

In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd.     

A pore-scale numerical model for flow through porous media

A pore-scale numerical model based on Smoothed Particle Hydrodynamics (SPH) is described for modelling fluid flow phenomena in porous media. Originally developed for astrophysics applications, SPH is extended to model incompressible flows of low Reynolds number as encountered in groundwater flow systems. In this paper, an overview of SPH is provided and the required modifications for modelling flow through porous media are described, including treatment of viscosity, equation of state, and no-slip boundary conditions. The performance of the model is demonstrated for two-dimensional flow through idealized porous media composed of spatially periodic square and hexagonal arrays of cylinders. The results are in close agreement with solutions obtained using the finite element method and published solutions in the literature. Copyright © 1999 John Wiley & Sons, Ltd.     

Modelling of coupled fluid‐mechanical problems in fractured geological media using enriched finite elements

Geological environments, such as petroleum reservoirs, normally exhibit physical discontinuities, for example, fractures and faults. Because of the reduced thickness of these discontinuities, finite element formulations with strong discontinuity have been applied to the numerical modelling of geological environments. Until now, two relevant characteristics of petroleum reservoirs have not been addressed by these formulations. The first is the pore pressure jump in the direction normal to a discontinuity in a fluid‐mechanical coupling condition, which is present primarily in sealing faults owing to the contrast of permeability with the porous medium. The absence of this jump can affect the prediction of the deformability of a physical discontinuity. Furthermore, reservoir models frequently use coarse meshes. Thus, the method used to evaluate the pore pressure in the discontinuity may exhibit a strong dependence relative to the mesh refinement. Based on these characteristics, in this study, a formulation of an enriched finite element for application to coupled fluid‐mechanical problems with pre‐existing physical discontinuities saturated by a single fluid is presented. The formulation employs discontinuous interpolation functions and enables the reproduction of jumps of displacement and pore pressure associated with a discontinuity inside the element without the need to discretise it. An approximation to estimate the pore pressure in the discontinuity was developed, one which seeks to minimise the influence of refinement. The element's response is verified by comparison with a one‐dimensional analytical solution and simple examples that are simulated using commercial software. Copyright © 2015 John Wiley & Sons, Ltd.     

油藏多相渗流的面向对象有限元程序设计

近年来,面向对象的编程技术越来越受到人们的重视,并用来解决和模拟各种物理问题和工程问题,特别是在有限元方法编程和模拟中得到了广泛地应用。介绍了面向对象的程序设计方法和它与有限元程序相结合的原理和优点。在此基础上提出了面向对象有限元方法在油藏多相渗流中的应用实例,并按照面向对象的程序设计方法,基于Windows XP平台采用Visual C++.NET语言,建立了有关描述多相渗流的有限元模型的类,给出了类的描述和它的实现方法。建立的相关类及其实现方法主要包括处理矩阵计算的类、节点类、单元类、材料类、形函数类等。据此编制了油藏多相渗流的有限元分析的数值计算程序,并给出了一个数值实例。     

A finite element algorithm for parameter identification of material models for fluid saturated porous media

In this contribution an algorithm for parameter identification of geometrically linear Terzaghi–Biot‐type fluid‐saturated porous media is proposed, in which non‐uniform distributions of the state variables such as stresses, strains and fluid pore pressure are taken into account. To this end a least‐squares functional consisting of experimental data and simulated data is minimized, whereby the latter are obtained with the finite element method. This strategy allows parameter identification based on in situ experiments. In order to improve the efficiency of the minimization process, a gradient‐based optimization algorithm is applied, and therefore the corresponding sensitivity analysis for the coupled two‐phase problem is described in a systematic manner. For illustrative purpose, the performance of the algorithm is demonstrated for a slope stability problem, in which a quadratic Drucker–Prager plasticity model for the solid and a linear Darcy law for the fluid are combined. Copyright © 2001 John Wiley & Sons, Ltd.