A Phase-Field-Based Approach for Modeling Flow and Geomechanics in Fractured Reservoirs

Modeling the geomechanical deformations of fracture networks has become an integral part of designing enhanced geothermal systems and recovery mechanisms for unconventional reservoirs. Stress changes in the reservoir can cause variations in the apertures of fractures resulting in large changes in their transmissivities. At the same time, sustained high-injection pressures can induce shear slipping along existing fractures and faults and trigger seismic activity. In this work, we extend the phase-field method to solve for flow and geomechanical deformations in naturally fractured reservoirs. The framework can predict the opening/closure of fractures as well as their shear slipping because of induced stresses and poromechanical effects. The flow through fractures is modeled on spatially nonconforming grids using the enhanced velocity mixed finite element method. The geomechanics equations are discretized using the continuous Galerkin (CG) finite element method. The flow and mechanics equations are decoupled using the fixed stress iterative scheme. The implementation is validated against the analytical solutions of Mandel's problem and Sneddon's benchmark test. Two synthetic examples are presented to illustrate the impact of poroelastic deformations and the accompanying dynamic behavior of fractures on the safety and productivity of subsurface projects.