An XFEM implementation in Abaqus to model intersections between fractures in porous rocks

Abstract This paper presents the key aspects of the implementation of an XFEM element for hydraulic fracturing in porous medium in Abaqus. The implemented element has the capability of intersecting fractures inside an element and includes coupled hydro-mechanical behaviour and fluid exchange between fracture and the porous medium. The algorithm has three main steps: fracture geometry pre-processing, local matrix assembling, and post-processing. The first computes all the geometrical parameters of the XFEM. The second computes the jacobian and the right-hand-side matrices for each mesh element. The third checks for fracture propagation and the direction and length for propagating segments. Considering that Abaqus is a commercial software, adaptations were necessary, such as an algorithm for consideration of in-situ stress states, positioning of the degrees of freedom and numerical integration. A discussion on the limitations of the tool and of the use of the User Subroutine feature closes the implementation description. Then, we present the consolidation of a multi-fractured medium, comparing the results with a model with interface elements. In a second example, we simulate an injection in a fracture that propagates and intersects a naturally fractured model. The simulations show that the capabilities of Abaqus were considerably extended with this implementation.

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

[2]  P. Hansbo,et al.  A finite element method for the simulation of strong and weak discontinuities in solid mechanics , 2004 .

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

[4]  M. Wheeler,et al.  An augmented-Lagrangian method for the phase-field approach for pressurized fractures , 2014 .

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

[6]  T. Belytschko,et al.  Extended finite element method for cohesive crack growth , 2002 .

[7]  Jmrj Jacques Huyghe,et al.  The enhanced local pressure model for the accurate analysis of fluid pressure driven fracture in porous materials , 2015 .

[8]  D. Roehl,et al.  Hydro-mechanical modeling of hydraulic fracture propagation and its interactions with frictional natural fractures , 2019, Computers and Geotechnics.

[9]  Panos Papanastasiou,et al.  The Effective Fracture Toughness in Hydraulic Fracturing , 1999 .

[10]  D. Roehl,et al.  An XFEM element to model intersections between hydraulic and natural fractures in porous rocks , 2018, International Journal of Rock Mechanics and Mining Sciences.

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

[12]  T. Belytschko,et al.  New crack‐tip elements for XFEM and applications to cohesive cracks , 2003 .

[13]  T. Belytschko,et al.  A method for growing multiple cracks without remeshing and its application to fatigue crack growth , 2004 .

[14]  Xiaowei Weng,et al.  Hydraulic Fracture Crossing Natural Fracture at Nonorthogonal Angles: A Criterion and Its Validation , 2012 .

[15]  Amir R. Khoei,et al.  An enriched FEM technique for modeling hydraulically driven cohesive fracture propagation in impermeable media with frictional natural faults: Numerical and experimental investigations , 2015 .

[16]  D. Roehl,et al.  Integrated discrete fracture and dual porosity - Dual permeability models for fluid flow in deformable fractured media , 2019, Journal of Petroleum Science and Engineering.

[17]  Yan Jin,et al.  A criterion for identifying hydraulic fractures crossing natural fractures in 3D space , 2014 .

[18]  Jon E. Olson,et al.  Numerical Modeling of Multistranded-Hydraulic-Fracture Propagation: Accounting for the Interaction Between Induced and Natural Fractures , 2009 .

[19]  John McLennan,et al.  A 3D peridynamic simulation of hydraulic fracture process in a heterogeneous medium , 2016 .

[20]  T. Ishida,et al.  The distinct element analysis for hydraulic fracturing in hard rock considering fluid viscosity and particle size distribution , 2011 .

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

[22]  Amir R. Khoei,et al.  Modeling the interaction between fluid-driven fracture and natural fault using an enriched-FEM technique , 2015, International Journal of Fracture.

[23]  Ahmad Ghassemi,et al.  Simulation of hydraulic fracture propagation near a natural fracture using virtual multidimensional internal bonds , 2011 .

[24]  Ted Belytschko,et al.  A method for dynamic crack and shear band propagation with phantom nodes , 2006 .

[25]  Ercan Gürses,et al.  A robust algorithm for configurational‐force‐driven brittle crack propagation with R‐adaptive mesh alignment , 2007 .

[26]  Eugenio Giner,et al.  An Abaqus implementation of the extended finite element method , 2009 .

[27]  Gen Li,et al.  Numerical Simulation of 3D Hydraulic Fracturing Based on an Improved Flow-Stress-Damage Model and a Parallel FEM Technique , 2012, Rock Mechanics and Rock Engineering.