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.

The three-point bending test is a key method for determining parameters related to the mechanical fracture properties of rocks. In this study, shale outcrops from Changning County, Sichuan Province, China, were selected. Three-point bending experiments were performed on shale semi-circular disks with a central straight crack, tested both perpendicular and parallel to the bedding direction. The corresponding load–displacement curves and crack opening displacements were obtained. The opening displacements of the specimens were measured through digital image technology, and the tensile strength and stiffness of the specimens were further calculated. A finite element model of the three-point bending test was developed. By integrating the finite element model with the experimentally obtained load–displacement curves, the anisotropic elastic moduli of the shale were inversely determined. Fracture toughness was calculated using two approaches: a formula from the International Society for Rock Mechanics and numerical methods using the finite element model, which was appropriately configured with the previously obtained elastic modulus values. The stress intensity factors for each specimen were calculated and compared. The energy release rate of shale was computed based on the fracture toughness. Results showed that both the fracture toughness and energy release rate of shale were greater in the perpendicular bedding direction than in the parallel direction. As an example, one specimen’s elastic modulus, opening displacement, and energy release rate obtained from experiments were input into the numerical simulation of the three-point bending test. The simulated load–displacement curve matched the experimental results well. This study provides a comprehensive approach to evaluating the anisotropic mechanical fracture properties of shale formations, which is essential for improving the accuracy of hydraulic fracture prediction models and enhancing the efficiency of shale gas extraction.

Welcome to e-content platform of John Benjamins Publishing Company. Here you can find all of our electronic books and journals, for purchase and download or subscriber access.