ABSTRACT: Creation of a fracture network in a hydraulic fracturing process is essential for subsurface energy extraction and CO2 sequestration. It is facilitated by reactivation of pre-existing intersecting weak layers and cemented cracks in the rock. In this study, a poromechanical model is developed for the hydraulic fracturing process in rocks containing such pre-existing weak layers. Based on the mixture theory, the crack band model is used to simulate the growth of a crack system. The governing equations with the parameters for hydromechanical coupling are derived, to describe the evolution of the opening and branching of cracks caused by water injection. Microplane model M7 is adopted to characterize the deformation and fracturing of the solid skeleton of the rock, and the Poiseuille law is used to characterize fluid flow through the hydraulic fractures. Numerical simulations are performed to reproduce and interpret recently published laboratory-scale hydraulic fracturing experiments conducted at Los Alamos National Laboratory (LANL). In these experiments, the rock was represented by confined plaster slabs containing orthogonal intersecting weak layers of higher porosity. Numerical simulations reveal how poromechanical characteristics such as the Biot coefficient and the fluid injection rate lead to various typical fracture modes observed in the experiments. These modes include formation of one dominant planar crack or various orthogonal fracture networks.
Xu et al. (Sun,) studied this question.