@Article{ee.2025.066033,
AUTHOR = {Xin Wan, Shuyi Li, Tiankui Guo, Ming Chen, Xing Yang, Guchang Zhang, Zi’ang Wang},
TITLE = {Numerical Simulation of Hydraulic Fracture Propagation in Deep Elasto-Plastic Reservoirs},
JOURNAL = {Energy Engineering},
VOLUME = {122},
YEAR = {2025},
NUMBER = {8},
PAGES = {3013--3039},
URL = {http://www.techscience.com/energy/v122n8/63086},
ISSN = {1546-0118},
ABSTRACT = {Hydraulic fracturing is a key technology for the efficient development of deep oil and gas reservoirs. However, fracture propagation behavior is influenced by rock elastoplasticity and thermal stress, making it difficult for traditional linear elastic models to accurately describe its dynamic response. To address this, this study employs the Continuum-Discontinuum Element Method (CDEM), incorporating an elastoplastic constitutive model, thermo-hydro-mechanical (THM) coupling effects, and cohesive zone characteristics at the fracture tip to establish a numerical model for hydraulic fracture propagation in deep elastoplastic reservoirs. A systematic investigation was conducted into the effects of fluid viscosity, reservoir temperature, injection rate, elastic modulus, and horizontal stress difference on fracture propagation. The findings show that a larger horizontal stress differential results in a more rectangular fracture geometry, a shorter fracture length, and a wider fracture. An increase in elastic modulus has a negligible impact on fracture length but reduces fracture width, resulting in a rounded rectangular morphology. Elevated reservoir temperature induces thermal tensile stress around the fracture, mitigating <i>in-situ</i> stress effects and reducing both breakdown and propagation pressures. Higher injection rates and fluid viscosity increase fracture initiation difficulty, promoting shorter but wider fractures with enhanced height growth beyond interlayer barriers. Additionally, horizontal stress significantly affects near-fracture plastic deformation: when the stress difference increases from 10 to 25 MPa, the maximum cumulative plastic strain in the surrounding rock rises by 66.67%. By integrating elastoplasticity and thermal stress effects, this study overcomes the limitations of conventional hydraulic fracturing simulations, offering novel insights for optimizing extraction strategies in deep unconventional reservoirs.},
DOI = {10.32604/ee.2025.066033}
}