Discrete fracture model for coupled flow and geomechanics

We present a fully implicit formulation of coupled flow and geomechanics for fractured three-dimensional subsurface formations. The Reservoir Characterization Model (RCM) consists of a computational grid, in which the fractures are represented explicitly. The Discrete Fracture Model (DFM) has been widely used to model the flow and transport in natural geological porous formations. Here, we extend the DFM approach to model deformation. The flow equations are discretized using a finite-volume method, and the poroelasticity equations are discretized using a Galerkin finite-element approximation. The two discretizations—flow and mechanics—share the same three-dimensional unstructured grid. The mechanical behavior of the fractures is modeled as a contact problem between two computational planes. The set of fully coupled nonlinear equations is solved implicitly. The implementation is validated for two problems with analytical solutions. The methodology is then applied to a shale-gas production scenario where a synthetic reservoir with 100 natural fractures is produced using a hydraulically fractured horizontal well.