A comparative study on simulating flow-induced fracture deformation in subsurface media by means of extended FEM and FVM

Accurate and efficient simulation on the fluid flow and deformation in porous media is of increasing importance in a diverse range of engineering fields. At present, there are only several methods can be used to simulate the deformation of fractured porous media. It is very important to know their application scopes, advantages, and disadvantages for solving the practical problems correctly. Therefore, in this paper, we compared two numerical simulation methods for flow-induced fracture deformation in porous media. One is the Extended Finite Element Method (XFEM), which is based on the classical finite element method and can simulate strong or weak discontinuous problems. The other is developed within the finite-volume framework, termed Extended Finite Volume Method (XFVM). We designed three test cases, including single fracture, cross fractures and eight discrete fractures, to investigate the accuracy and efficiency of XFEM and XFVM. The reference solutions were provided by the commercial software, COMSOL , where the standard finite element method is implemented. The research findings showed that the accuracy of the XFEM was slightly higher than that of the XFVM, but the latter was more efficient. These results are likely to be useful in decision making regarding choice of solving methods for the multi-field coupling problem in fractured porous media.

[1]  B. Valley,et al.  THERMAID - A matlab package for thermo-hydraulic modeling and fracture stability analysis in fractured reservoirs , 2018, 1806.10942.

[2]  Abimael F. D. Loula,et al.  On stability and convergence of finite element approximations of biot's consolidation problem , 1994 .

[3]  Peng Wang,et al.  Implementation of the IDEAL algorithm for complex steady-state incompressible fluid flow problems in OpenFOAM , 2017 .

[4]  Fusheng Yang,et al.  A novel POD reduced-order model based on EDFM for steady-state and transient heat transfer in fractured geothermal reservoir , 2020 .

[5]  Philip Cardiff,et al.  On finite volume method implementation of poro‐elasto‐plasticity soil model , 2015 .

[6]  Craig T. Simmons,et al.  Finite volume coupling strategies for the solution of a Biot consolidation model , 2014 .

[7]  R. Deb.  Numerical modeling of fluid injection induced shear failure, tensile opening and flow-mechanics coupling , 2018 .

[8]  Jan Martin Nordbotten.  Finite volume hydromechanical simulation in porous media , 2014, Water resources research.

[9]  Eirik Keilegavlen,et al.  A finite-volume discretization for deformation of fractured media , 2018, Computational Geosciences.

[10]  Xavier Raynaud,et al.  Comparison between cell-centered and nodal-based discretization schemes for linear elasticity , 2018, Computational Geosciences.

[11]  Yu-Shu Wu,et al.  Parallel simulation of fully-coupled thermal-hydro-mechanical processes in CO2 leakage through fluid-driven fracture zones , 2015 .

[12]  Patrick Jenny,et al.  Modeling of Failure along Predefined Planes in Fractured Reservoirs , 2014 .

[13]  K. Terzaghi.  Theoretical Soil Mechanics , 1943 .

[14]  P. Jenny,et al.  Modeling of shear failure in fractured reservoirs with a porous matrix , 2017, Computational Geosciences.

[15]  P. Jenny,et al.  Finite volume–based modeling of flow‐induced shear failure along fracture manifolds , 2017 .

[16]  Gunnar Jansen,et al.  On the Role of Thermal Stresses during Hydraulic Stimulation of Geothermal Reservoirs , 2017 .

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

[18]  Arnaud G. Malan,et al.  An enhanced finite volume method to model 2D linear elastic structures , 2014 .

[19]  D. W. Peaceman.  Interpretation of well-block pressures in numerical reservoir simulation with nonsquare grid blocks and anisotropic permeability , 1983 .

[20]  Ruben Juanes,et al.  A locally conservative finite element framework for the simulation of coupled flow and reservoir geomechanics , 2007 .

[21]  Nadine Lanteri,et al.  Gas Seepage along the Edge of the Aquitaine Shelf (France): Origin and Local Fluxes , 2017 .

[22]  He Liu,et al.  An XFEM-based method with reduction technique for modeling hydraulic fracture propagation in formations containing frictional natural fractures , 2017 .

[23]  Kamy Sepehrnoori,et al.  Coupled Geomechanics and Flow Simulation for an Embedded Discrete Fracture Model , 2013, ANSS 2013.

[24]  Yang Li,et al.  An efficient hydro-mechanical model for coupled multi-porosity and discrete fracture porous media , 2018 .

[25]  Dongxiao Zhang,et al.  A fully coupled thermo-hydro-mechanical, three-dimensional model for hydraulic stimulation treatments , 2016 .

[26]  D. W. Peaceman.  Interpretation of Well-Block Pressures in Numerical Reservoir Simulation(includes associated paper 6988 ) , 1978 .

[27]  Stéphane Bordas,et al.  An extended finite element library , 2007 .

[28]  Olav Møyner,et al.  The multiscale restriction smoothed basis method for fractured porous media (F-MsRSB) , 2016, J. Comput. Phys..

[29]  H. Tchelepi,et al.  Stability and convergence of sequential methods for coupled flow and geomechanics: Fixed-stress and fixed-strain splits , 2011 .

[30]  Mark D. Zoback,et al.  Fully Coupled Nonlinear Fluid Flow and Poroelasticity in Arbitrarily Fractured Porous Media: A Hybrid‐Dimensional Computational Model , 2017 .

[31]  A. Cheng,et al.  Mandel's problem revisited , 1996 .

[32]  M. F. Lough,et al.  Hierarchical modeling of flow in naturally fractured formations with multiple length scales , 2001 .

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

[34]  Patrick Jenny,et al.  Adaptive Hierarchical Fracture Model for Enhanced Geothermal Systems , 2016, Multiscale Model. Simul..

[35]  Jihoon Kim,et al.  Constrained energy minimization based upscaling for coupled flow and mechanics , 2018, J. Comput. Phys..