An extended finite volume method and fixed‐stress approach for modeling fluid injection–induced tensile opening in fractured reservoirs

Fluid injection–induced tensile opening is modeled using an extended finite volume method (XFVM). An embedded fracture strategy is used for the flow problem, that is, the fractures are discretized using finite volume segments without resolving the grid around them. Further, the discontinuities across fractures are modeled using special basis functions. The fracture openings due to enhanced fluid pressure and the associated shear slip due to traction free boundary condition on the fracture segments are both modeled using these special discontinuity basis functions. Mass transfer between fractures and matrix is modeled using the pressure difference. The enhancement of fracture storativity due to tensile opening leads to stronger coupling between flow and mechanics. An iterative scheme relying on the fixed‐stress approach for fractures, which conserves the stress dependent terms over each iteration of the flow problem, has been introduced. Tensile opening has been simulated for single fractures embedded in two‐ and three‐dimensional matrices. The convergence criterion for sequentially implicit fixed‐stress scheme for fractures embedded in elastic media is established and has been validated numerically. Further, for 2D simulations, the effect of the matrix permeability for fracture propagation due to tensile opening has been studied.

[1]  J. Geertsma,et al.  A Rapid Method of Predicting Width and Extent of Hydraulically Induced Fractures , 1969 .

[2]  S. L. Crouch,et al.  Boundary element methods in solid mechanics , 1983 .

[3]  M. F. Lough,et al.  Efficient Finite-Difference Model for Flow in a Reservoir With Multiple Length-Scale Fractures , 2000 .

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

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

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

[7]  M. Belayneh,et al.  Control-volume finite-element two-phase flow experiments with fractured rock represented by unstructured 3D hybrid meshes , 2005 .

[8]  S. G. Roberts,et al.  Numerical simulation of multi-phase fluid flow in structurally complex reservoirs , 2007 .

[9]  M. Belayneh,et al.  Finite Element - Node-Centered Finite-Volume Two-Phase-Flow Experiments With Fractured Rock Represented by Unstructured Hybrid-Element Meshes , 2007 .

[10]  R. Borja Assumed enhanced strain and the extended finite element methods: A unification of concepts , 2008 .

[11]  Sebastian Geiger,et al.  A parallel FE-FV scheme to solve fluid flow in complex geologic media , 2008, Comput. Geosci..

[12]  T. Belytschko,et al.  The extended/generalized finite element method: An overview of the method and its applications , 2010 .

[13]  R. Horne,et al.  Investigation of injection-induced seismicity using a coupled fluid flow and rate/state friction model , 2011 .

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

[15]  Patrick Jenny,et al.  A hierarchical fracture model for the iterative multiscale finite volume method , 2011, J. Comput. Phys..

[16]  A. Ghassemi,et al.  A 3D Thermal-Poroelastic Model For Naturally Fractured Geothermal Reservoir Stimulation , 2012 .

[17]  D. Karvounis Simulations of Enhanced Geothermal Systems with an Adaptive Hierarchical Fracture Representation , 2013 .

[18]  Roland N. Horne,et al.  An investigation of stimulation mechanisms in Enhanced Geothermal Systems , 2014 .

[19]  J. Nordbotten Cell‐centered finite volume discretizations for deformable porous media , 2014 .

[20]  R. Horne,et al.  An embedded fracture modeling framework for simulation of hydraulic fracturing and shear stimulation , 2016, Computational Geosciences.

[21]  Joshua A. White,et al.  Accuracy and convergence properties of the fixed‐stress iterative solution of two‐way coupled poromechanics , 2015 .

[22]  Ahmad Ghassemi,et al.  Fast multipole displacement discontinuity method (FM‐DDM) for geomechanics reservoir simulations , 2015 .

[23]  M. Wheeler,et al.  Convergence of iterative coupling of geomechanics with flow in a fractured poroelastic medium , 2016, Computational Geosciences.

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

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

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

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