Study on the Evolution Law of Four-Dimensional In Situ Stress During Hydraulic Fracturing of Deep Shale Gas Reservoir

The increasing burial depth of deep shale formations in the southern Sichuan leads to more complex in situ stresses and geological structures, which in turn raises the challenges of hydraulic fracturing. Although enlarging the treatment scale and injection rate can enhance reservoir stimulation, the intensive development of faults and fractures in deep shale formations aggravates stress instability, inducing casing deformation, fracture communication, and other engineering risks that constrain efficient shale gas production. In this study, a cross-scale geomechanical model linking the regional to near-wellbore domains was constructed. A dynamic evolution equation was established based on flow–stress coupling, and a numerical conversion from the geological model to the finite element model was implemented through self-programming, thereby developing a simulation method for dynamic geomechanical evolution during hydraulic fracturing. Results indicate that dynamic variations in pore pressure dominate stress redistribution, while near-wellbore heterogeneity and mechanical property distribution significantly affect prediction accuracy. The injection of fracturing fluid generates a high-pressure gradient that drives pore pressure diffusion along fracture networks and faults, exhibiting a strong near-wellbore but weak far-field non-steady spatial attenuation. As the pore pressure increases, the peak value reaches 1.4 times the original pressure. The triaxial stress shows a negative correlation and decreases. The horizontal minimum principal stress shows the most significant drop, with a reduction of 15.79% to 20.68%, while the vertical stress changes the least, with a reduction of 5.7% to 7.14%. Compared with the initial stress state, the horizontal stress difference increases during fracturing. Rapid pore-pressure surges and fault distributions further trigger stress reorientation, with the magnitude of rotation positively correlated with the intensity of pore-pressure variation. The high porosity and permeability characteristics of the initial fracture network lead to a rapid attenuation of the stress around the wellbore. In the middle and later stages, as the pressure balance is achieved through fracture filling, the pore pressure rises and the stress decline tends to stabilize. The findings provide significant insights into the dynamic stress evolution of deep shale reservoirs during fracturing and offer theoretical support for optimizing fracturing design and mitigating engineering risks.