Simulation of Complex Hydraulic Fracture Propagation in Shale with Interlayers

Shale gas, as an unconventional resource, requires hydraulic fracturing to create complex fracture networks due to its low porosity and permeability. However, the presence of interlayers significantly affects fracture propagation, leading to highly complex fracture morphologies. This study focuses on the interbedded shale of the WJP Formation in southern China. A three-dimensional block discrete element method (BDEM) was employed to establish a hydraulic fracture propagation model, systematically investigating the effects of geological parameters (stress difference, interlayer thickness), engineering parameters pumping rate, fluid volume, viscosity), and perforation parameters (cluster number, cluster spacing, perforation location) on fracture network morphology. The results indicate that: (1) Among geological parameters, interlayer thickness is the key factor inhibiting vertical fracture propagation. Due to the influence of interlayers, an increase in stress difference promotes fracture length but suppresses fracture height and stimulated reservoir volume (SRV); (2) For engineering parameters, there exists a “threshold effect” for pumping rate and fluid volume, with 16 m3/min and 2000 m3 identified as the critical thresholds for interlayer breakthrough. Low viscosity (1 mPa·s) is conducive to forming complex fracture networks, while high viscosity extends fracture length but reduces SRV; (3) Regarding perforation parameters, the optimal stimulation effect is achieved with 6–7 clusters, a cluster spacing of 10 m, and perforation locations in the center of the main shale layer (19.85–21.6 m); (4) By introducing grey relational analysis, the degree of correlation between various influencing factors and the response to interlayer breakthrough is systematically evaluated based on the breakthrough conditions under different factors. Thin interlayers or low stress differences can reduce the critical pumping rate, whereas thick interlayers (≥3 m) become the primary constraint, making breakthrough difficult even at high pumping rates. Reliable interlayer breakthrough requires the simultaneous satisfaction of Δσ ≤ 16 MPa, h < 1 m, and Q ≥ 16 m3/min. The reliability of the model was verified by comparing numerical simulation results with field microseismic data. This study reveals the extension laws of complex fracture networks in interbedded shale, providing a theoretical basis for fracturing design and development optimization.