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

SUMMARY 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.

[1]  Nam H. Tran,et al.  Coupled Fluid Flow Through Discrete Fracture Network: A Novel approach , 2007 .

[2]  D. F. McTigue,et al.  Thermoelastic response of fluid‐saturated porous rock , 1986 .

[3]  J. E. Warren,et al.  The Behavior of Naturally Fractured Reservoirs , 1963 .

[4]  Bernhard A. Schrefler,et al.  Mapped infinite elements in soil consolidation , 1987 .

[5]  T. Belytschko,et al.  MODELING HOLES AND INCLUSIONS BY LEVEL SETS IN THE EXTENDED FINITE-ELEMENT METHOD , 2001 .

[6]  T. Rabczuk,et al.  Three-dimensional crack initiation, propagation, branching and junction in non-linear materials by an extended meshfree method without asymptotic enrichment , 2008 .

[7]  T. Belytschko,et al.  Extended finite element method for three-dimensional crack modelling , 2000 .

[8]  M. Biot.  General Theory of Three‐Dimensional Consolidation , 1941 .

[9]  J. Mandel.  Consolidation Des Sols (Étude Mathématique) , 1953 .

[10]  Bernhard A. Schrefler,et al.  Non‐isothermal consolidation of unbounded porous media using mapped infinite elements , 1987 .

[11]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[12]  Ted Belytschko,et al.  Modelling crack growth by level sets in the extended finite element method , 2001 .

[13]  Jean-Herve Prevost,et al.  MODELING QUASI-STATIC CRACK GROWTH WITH THE EXTENDED FINITE ELEMENT METHOD PART II: NUMERICAL APPLICATIONS , 2003 .

[14]  Milind Deo,et al.  Finite element, discrete‐fracture model for multiphase flow in porous media , 2000 .

[15]  Derek Elsworth,et al.  FLOW-DEFORMATION RESPONSE OF DUAL-POROSITY MEDIA , 1992 .

[16]  T. Belytschko,et al.  Arbitrary branched and intersecting cracks with the eXtended Finite Element Method , 2000 .

[17]  B. Schrefler,et al.  The Finite Element Method in the Static and Dynamic Deformation and Consolidation of Porous Media , 1998 .

[18]  T. Belytschko,et al.  Non‐planar 3D crack growth by the extended finite element and level sets—Part I: Mechanical model , 2002 .

[19]  Ted Belytschko,et al.  An extended finite element method with higher-order elements for curved cracks , 2003 .

[20]  小山 毅,et al.  拡張有限要素法(XFEM)・一般化有限要素法(GFEM)を用いた材料モデリングのレビュー Ted Belytschko,Robert Gracie and Giulio Ventura:A Review of Extended/Generalized Finite Element Methods for Material Modeling [Modeling and Simulations in Materials Science and Engineering, Vol.17, 043001, June 2009](構造,文献抄録) , 2010 .

[21]  Ted Belytschko,et al.  Elastic crack growth in finite elements with minimal remeshing , 1999 .

[22]  T. Belytschko,et al.  A review of extended/generalized finite element methods for material modeling , 2009 .

[23]  F. Brezzi.  On the existence, uniqueness and approximation of saddle-point problems arising from lagrangian multipliers , 1974 .

[24]  R. Al-Khoury,et al.  A computational model for fracturing porous media , 2007 .

[25]  C. Cryer,et al.  A COMPARISON OF THE THREE-DIMENSIONAL CONSOLIDATION THEORIES OF BIOT AND TERZAGHI , 1963 .

[26]  J. S. Y. Wang,et al.  Validity of cubic law for fluid flow in a deformable rock fracture. Technical information report No. 23 , 1979 .

[27]  T. Belytschko,et al.  A First Course in Finite Elements: Belytschko/A First Course in Finite Elements , 2007 .

[28]  Julien Réthoré,et al.  A Numerical Approach for Arbitrary Cracks in a Fluid-Saturated Medium , 2006 .

[29]  L. Durlofsky,et al.  An Efficient Discrete-Fracture Model Applicable for General-Purpose Reservoir Simulators , 2004 .

[30]  Derek Elsworth,et al.  Coupled Processes in Subsurface Deformation, Flow, and Transport , 2000 .

[31]  T. Belytschko,et al.  Non‐planar 3D crack growth by the extended finite element and level sets—Part II: Level set update , 2002 .

[32]  Ted Belytschko,et al.  A finite element method for crack growth without remeshing , 1999 .

[33]  T. Belytschko,et al.  A method for multiple crack growth in brittle materials without remeshing , 2004 .

[34]  G. I. Barenblatt,et al.  Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata] , 1960 .

[35]  Roland W. Lewis,et al.  A FINITE ELEMENT DOUBLE POROSITY MODEL FOR HETEROGENEOUS DEFORMABLE POROUS MEDIA , 1996 .

[36]  Roland W. Lewis,et al.  A novel finite element double porosity model for multiphase flow through deformable fractured porous media , 1997 .

[37]  Olivier Gosselin,et al.  Coupled Deformation And Fluid Flow In Fractured Porous Media Using Dual Permeability And Explicitly Defined Fracture Geometry , 2010 .

[38]  René de Borst,et al.  A discrete model for the dynamic propagation of shear bands in a fluid‐saturated medium , 2007 .

[39]  N. Moës,et al.  Improved implementation and robustness study of the X‐FEM for stress analysis around cracks , 2005 .

[40]  Michel Salaün,et al.  High‐order extended finite element method for cracked domains , 2005 .

[41]  Malgorzata Peszynska,et al.  Coupled fluid flow and geomechanical deformation modeling , 2003 .

[42]  Mark A Fleming,et al.  ENRICHED ELEMENT-FREE GALERKIN METHODS FOR CRACK TIP FIELDS , 1997 .

[43]  T. Belytschko,et al.  A first course in finite elements , 2007 .

[44]  Sung-Uk Choi,et al.  A comparative study of dual-porosity model and discrete fracture network model , 1999 .

[45]  T. Belytschko,et al.  Vector level sets for description of propagating cracks in finite elements , 2003 .

[46]  I. Babuska.  The finite element method with Lagrangian multipliers , 1973 .

[47]  Benoit Prabel,et al.  Level set X‐FEM non‐matching meshes: application to dynamic crack propagation in elastic–plastic media , 2007 .

[48]  I. Babuska.  Error-bounds for finite element method , 1971 .

[49]  J. Prévost,et al.  Modeling quasi-static crack growth with the extended finite element method Part I: Computer implementation , 2003 .

[50]  M. Duflot.  A study of the representation of cracks with level sets , 2007 .