Supercritical carbon dioxide (ScCO2) fracturing is a promising waterless stimulation technology for unconventional reservoirs, but the mechanisms controlling fracture-network evolution in naturally fractured formations remain unclear. In this study, a stress–flow–damage coupled numerical model based on the continuum–discontinuum element method (CDEM) was established to investigate hydraulic fracture propagation and natural fracture activation during ScCO2 fracturing. Water-based fracturing and ScCO2 fracturing were first compared under an identical symmetric natural fracture configuration. The effects of horizontal principal stress difference and natural fracture approach angle on fracture morphology, injection-pressure response, and fracture-width evolution were then analyzed. Finally, ScCO2-induced fracture-network development in a complex natural fracture network was examined under different injection rates and stress differences. The results show that water-based fracturing is dominated by rapid pressure buildup and relatively straight main-fracture propagation, with a peak injection pressure of 22.9 MPa at 0.075 s and a fracture width of approximately 1.0 mm at 8 s. In contrast, ScCO2 fracturing shows delayed pressure buildup, stronger pressure fluctuations, more pronounced natural fracture activation, and tortuous step-like fracture paths, while producing a smaller fracture width of approximately 0.44 mm at 120 s. A lower stress difference and lower approach angle favor fracture deflection and natural fracture activation, whereas a higher stress difference strengthens stress-guided main-fracture propagation. In complex fracture networks, a lower injection rate promotes more sufficient ScCO2–natural fracture interaction, while a higher injection rate enhances main-fracture extension. These findings provide a numerical basis for optimizing ScCO2 fracturing in naturally fractured unconventional reservoirs.



