Phase-Field Modeling of Thermal Fracturing Mechanisms in Reservoir Rock Under High-Temperature Conditions

Thermal stimulation represents an effective method for enhancing reservoir permeability, thereby improving geothermal energy recovery in Enhanced Geothermal Systems (EGS). The phase-field method (PFM) has been widely adopted for its proven capability in modeling the fracture behavior of brittle solids. Consequently, a coupled thermo-mechanical phase-field model (TM-PFM) was developed in COMSOL 6.2 Multiphysics to probe thermal fracturing mechanisms in reservoir rocks. The TM-PFM was validated against the analytical solutions for the temperature and stress fields under steady-state heat conduction in a thin-walled cylinder, three-point bending tests, and thermal shock tests. Subsequently, two distinct thermal fracturing modes in reservoir rock under high-temperature conditions were investigated: (i) fracture initiation driven by sharp temperature gradients during instantaneous thermal shocks, and (ii) crack propagation resulting from heterogeneous thermal expansion of constituent minerals. The proposed TM-PFM has been validated through systematic comparison between the simulation results and the corresponding experimental data, thereby demonstrating its capability to accurately simulate thermal fracturing. These findings provide mechanistic insights for optimizing geothermal energy extraction in EGS.