A Mechanics-Based Recursive Propagation Framework for Modeling Complex Hydraulic Fracture Networks in Naturally Fractured Shale Reservoirs

Hydraulic fracturing in naturally fractured shale reservoirs commonly generates complex mesh-like fracture networks governed by hydraulic fracture–natural fracture interactions, which strongly affect stimulated volume, fracture connectivity, and early-time production. Existing simulation and monitoring-based methods often cannot simultaneously capture interaction mechanisms, rapidly generate field-scale fracture networks, and validate production responses. This study proposes a mechanics-constrained recursive propagation framework. A field-constrained stochastic natural-fracture model is first constructed, an explicit hydraulic fracture–natural fracture interaction criterion is incorporated to identify penetration, opening, and shear slipping, and a fully vectorized bidirectional recursive algorithm is developed to efficiently generate complex fracture networks. The method is applied to a 40-stage fractured horizontal well in the Changqing Oilfield, where the target interval has a porosity of 6.1%, a permeability of 0.1 mD, and a horizontal stress contrast of 7.0 MPa. The simulated network reproduces crossing, arrest, unilateral diversion, and bilateral diversion, and agrees well with microseismic observations. EDFM-based fully implicit flow simulation further shows early-time production deviations of 2–10%. These results demonstrate that the proposed framework can efficiently generate physically plausible field-scale fracture networks for fracturing design, post-fracturing evaluation, and short-term production forecasting.