Hydraulic Fracturation Simulated by a 3D Coupled HM-DEM Model

"Numerical methods used to study hydraulic fracturing in low permeability rock masses need to take into account the interaction between the progressive failure mechanism of the intact rock matrix and the pre-existing natural fractures. A fully coupled hydro-mechanical 3D discrete elements model dealing with such interaction is presented here. The ultimate objective is to grasp the role of the rock-matrix failure on the propagation of the main fluid driven fracture under different stress fields and pre-existing fracture orientations. The proposed method provides the spatiotemporal distribution of the induced crack events which can be tracked during the simulations, making possible the identification of the local mechanisms leading to macro-failure in the rock matrix or reactivation of the discontinuity planes. Fluid flow inside a fracture, pressurized crack and hydro-fracking experimental tests were simulated as part of the validation process of the model and, are presented in this paper.INTRODUCTIONWith the expansion of shale gas and geothermal reservoirs exploitation, the technique of Hydraulic Fracturing (HF) has been widely used over the past years (McClure, 2012; Britt, 2012). Despite recent advances in micro-seismicity analysis or numerical modeling (Nagel et al., 2013), the way progressive failure of rocks affects the fluid flow depending on the in situ stress conditions or the rock heterogeneities is not fully controlled yet (Zoback, 2010).In this work we focused on setting up a three-dimensional discrete elements model to study fluid driven fractures propagation. The intact rock is modeled by a discrete element method (DEM) able to reproduce fracture initiation and propagation (Scholtès and Donzé, 2013, Garcia et al., 2013). The pre-existing fractures are explicitly modeled and can be individually defined or be part of a discrete fracture network (DFN) (Harthong et al., 2012). Their mechanical response is fully coupled with the fluid within a HM-DEM framework. The fluid flow and it’s interaction with the intact matrix and the pre-existing discontinuities is modeled through a Pore-scale Finite Volume (PFV) scheme specially enhanced for fractured rock modeling (Chareyre et al., 2012, Catalano et al., 2014) ."