Multi-Factor Coupling Simulation and Mechanism Investigation of Hydraulic Fracture Propagation in Hard Roof
Abstract
The propagation behavior of hydraulic fractures in deep hard roofs is governed by the multi-factor coupling of injection rate, in situ stress, and fluid viscosity. Taking the Cuimu coal mine as the engineering background, this study systematically investigates the effects of injection rate, lateral pressure coefficient, and fluid viscosity on fracture propagation through numerical simulation. Theoretical derivations based on the KGD (Kristianovich-Geertsma-de Klerk) model are further integrated and validated by field tests. The results indicate a critical injection rate of 6 × 10−8 m3/s, above which the marginal increase in acoustic emission events declines significantly. Increasing the lateral pressure coefficient from 1.0 to 3.5 shifts the fracture pattern from relatively simple to increasingly complex and interwoven, accompanied by a logarithmic increase in fractal dimension from 1.38 to 1.80. The total acoustic emission count rises to a peak of 125,910 as viscosity increases from 0.001 Pa·s to 0.5 Pa·s, then drops to 44,061 at 1.0 Pa·s, showing a unimodal trend. Theoretical analysis shows that during the propagation stage, the fracture length follows L∝Q1/2, and the maximum fracture opening follows wmax∝Q1/3. The lateral pressure coefficient controls the complexity of the fracture network through the directional distribution of stress intensity factors. Field tests at the Cuimu coal mine adopted a combination of stepwise injection rate and low-viscosity fluid, together with borehole densification and interval-skipping fracturing sequences. The effective fracturing radius reached 25~30 m, roof convergence was reduced by 31%, and the proportion of high-energy microseismic events decreased from 12% to 4%. This study establishes a complete theoretical framework from initiation theory to propagation dynamics and then to multi-crack competition, providing both a theoretical basis and engineering example for optimizing fracturing parameters in hard roofs under high stress anisotropy.