Experimental and Numerical Investigation of the Effect of Bedding Layer Orientation on Fracture Toughness of Shale Rocks

Yu Suo,Zhixi Chen,Sheikh S. Rahman,Huifang Song
DOI: https://doi.org/10.1007/s00603-020-02131-1
IF: 6.518
2020-04-23
Rock Mechanics and Rock Engineering
Abstract:Fracture toughness is a critical parameter responsible for fracture initiation and propagation. Fracture toughness in shale gas reservoir, however, is highly variable because of its anisotropy and spatial variation of clay content. In this paper, a series of laboratory and numerical experiments are carried out to estimate shale fracture toughness at different bedding plane orientations with respect to loading direction (angle of bedding layer with respect to Cracked Chevron notched Brazilian disc (CCNBD) test samples) and in different brine solutions. The CCNBD test was conducted on 12 cylindrical samples for fracture toughness. Shale samples were prepared at four different angles (0°, 30°, 45° and 90°) relative to the bedding plane. The prepared specimens were saturated in potassium chloride (KCL) solutions of different concentrations. The laboratory results of toughness have shown to be highly variable with respect to both bedding plane and brine concentration and that the sample at the angle of 90° in 4% KCL concentration exhibited the highest fracture toughness. Numerical simulations based on extended finite element method (XFEM) were also carried out to simulate fracture evolution and propagation in CCNBD samples with different bedding planes. The results have shown that the bedding layers caused the fracture path to deflect. The deviation from straight crack path is caused by mixed mode fracture initiation and propagation instead of tension mode. The mixed mode fracture propagation behavior was verified by analyzing fracture propagation path using both laboratory experiments and numerical simulation.
geosciences, multidisciplinary,engineering, geological
What problem does this paper attempt to address?
The problems that this paper attempts to solve are: how does the fracture toughness of shale rocks change under different bedding plane directions and different brine solution conditions, and what impact does this change have on crack initiation and propagation during the hydraulic fracturing process? ### Problem Background As an important natural resource, the extraction efficiency of shale gas is significantly affected by geological conditions and rock mechanical properties. Shale rocks have obvious anisotropy, especially large differences in mechanical properties in the bedding plane direction. In addition, the clay mineral content in shale rocks and its interaction with fluids also affect the mechanical properties of the rocks. Therefore, understanding the fracture toughness of shale rocks under different conditions is crucial for optimizing the hydraulic fracturing process. ### Research Objectives 1. **Experimental Research**: - Measure the fracture toughness of shale rocks at different bedding plane angles (0°, 30°, 45° and 90°). - Analyze the influence of potassium chloride (KCl) solutions with different concentrations on the fracture toughness of shale rocks. 2. **Numerical Simulation**: - Use the extended finite element method (XFEM) to simulate the crack propagation paths and behaviors at different bedding plane angles. - Study the influence of bedding planes on crack propagation paths, especially the mixed - mode crack propagation behavior. ### Key Findings - **Experimental Results**: - The bedding plane angle has a significant impact on the fracture toughness of shale rocks. When the bedding plane is perpendicular to the loading direction (90°), the fracture toughness is the highest; when it is parallel to the loading direction (0°), the fracture toughness is the lowest. - KCl solutions with different concentrations also have an impact on the fracture toughness. Under the condition of 4% KCl solution, the fracture toughness of shale rocks is the highest. - **Numerical Simulation Results**: - The numerical simulation verifies the experimental results, indicating that the bedding plane will cause the crack path to deflect, and the crack propagation is mainly in the mixed - mode rather than the pure tensile mode. - The simulation results show that when the Young's modulus of the bedding plane is different from that of the matrix material, the crack propagation path and stress distribution will also change. ### Formula Summary The calculation formula for fracture toughness is: \[ K_{IC} = \frac{F_{max}}{B \sqrt{R Y^*}} \] where: - \( F_{max} \) is the maximum axial load when the specimen fails; - \( B \) is the thickness of the specimen; - \( R \) is the radius of the specimen; - \( Y^* \) is a dimensionless stress intensity factor, defined as: \[ Y^* = u e^{v / \alpha_0} \] where \( u \) and \( v \) are determined by \( \alpha_0 \) and \( \alpha_B \) respectively. Through these studies, the authors hope to provide theoretical basis and technical support for the hydraulic fracturing design of shale gas fields, thereby improving the extraction efficiency of shale gas.