The Impacts of Increased Ductility in Organic-Rich Shale on Fracture Network Growth by Hydraulic Fracturing

The difficulty of hydraulic fracturing in organic-rich shale caused by the increased ductility has not been well interpreted quantitatively, although it is well perceived that the increased shale ductility can impede the propagation of hydraulic fractures and enhance the healing of created fractures upon injection shutdown. This study aims to quantitatively study the impacts of increased ductility on the stimulated reservoir volume (SRV) using an advanced XFEM-based simulator. To achieve this goal, a modified cohesive zone model has been integrated into an in-house fully coupled poroelastic XFEM framework. The study continues by extending the functionality of the numerical framework to simulating multiple interacting fractures. The utilization of the object-oriented programming paradigm in the development of the framework makes it an easy extension to include the multi-fracture network by creating more instances of crack segments. A main hydraulic fracture with an arbitrary number of intersected branches can thus be modeled. A series of parametric studies will be conducted to investigate the impacts of increased ductility on the induced SRV by varying four involved material parameters individually. The modified cohesive zone model, which is essentially a traction-separation law (TSL), is characterized by four parameters: the initial tensile strength Tini, ultimate tensile strength Tkrg, the critical separation Dc, and the final crack separation Dmax. It can flexibly model different crack opening scenarios and simulate more realistically the increased shale ductility. The fully coupled poroelastic XFEM framework has been comprehensively verified against the latest semi-analytical solutions on the four well-known propagation regimes. The numerical results show that the shape of TSL does affect the main hydraulic fracture growth as well as the evolvement of the fracture network, given the same cohesive crack energy and tensile strength. It infers that ductility is not only controlled by cohesive crack energy and tensile strength, which further indicates the necessity of the newly proposed cohesive zone model. The magnitude of the initial tensile strength, controlling when the cohesive crack starts propagating, is found to have the greatest impacts on the fracture length, and SRV, among all four TSL parameters. The novelty of this study is two-fold. First, the newly modified cohesive zone model can more realistically represent the increased shale ductility. Second, the advanced XFEM framework that enables the simulation of a fracture network can study the impacts of increased ductility on the whole SRV but not a single crack.