Propagation Characteristics of Multi-Cluster Hydraulic Fracturing in Shale Reservoirs with Natural Fractures

Hydraulic fracturing of gas and oil reservoirs is the primary stimulation method for enhancing production in the field of petroleum engineering. The hydraulic fracturing technology plays a crucial role in increasing shale gas production from shale reservoirs. Understanding the effects of reservoir and fracturing conditions on fracture propagation is of great significance for optimizing the hydraulic fracturing process and has not been adequately explored in the current literature. In the context of shale reservoirs in Yibin, Sichuan Province, China, the study selects outcrops to prepare samples for uniaxial compression and Brazilian splitting tests. These tests measure the compressive and tensile strengths of shale in parallel bedding and vertical bedding directions, obtaining the shale’s anisotropic elastic modulus and Poisson’s ratio. These parameters are crucial for simulating reservoir hydraulic fracturing. This paper presents a numerical model utilizing a finite element (FE) analysis to simulate the process of multi-cluster hydraulic fracturing in a shale reservoir with natural fractures in three dimensions. A numerical simulation of the intersection of multiple clusters of 3D hydraulic fractures and natural fractures was performed, and the complex 3D fracture morphologies after the interaction between any two fractures were revealed. The influences of natural fractures, reservoir ground stress, fracturing conditions, and fracture interference concerning the spreading of hydraulic fractures were analyzed. The results highlight several key points: (1) Shale samples exhibit distinct layering with significant anisotropy. The elastic compressive modulus and Poisson’s ratio of parallel bedding shale samples are similar to those of vertical bedding shale samples, while the compressive strength of parallel bedding shale samples is significantly greater than that of vertical bedding shale samples. The elastic compressive modulus of shale is 6 to 10 times its tensile modulus. (2) The anisotropy of shale’s tensile properties is pronounced. The ultimate load capacity of vertical bedding shale samples is 2 to 4 times that of parallel bedding shale samples. The tensile strength of vertical bedding shale samples is 2 to 5 times that of parallel bedding shale samples. (3) The hydraulic fractures induced by the injection well closest to the natural fractures expanded the fastest, and the natural fractures opened when they intersected the hydraulic fractures. When the difference in the horizontal ground stress was significant, natural fractures were more inclined to open after the intersection between the hydraulic and natural fractures. (4) The higher the injection rate and viscosity of the fracturing fluid, the faster the fracture propagation. The research findings could improve the fracturing process through a better understanding of the fracture propagation process and provide practical guidance for hydraulic fracturing design in shale gas reservoirs.