Thermal-hydraulic-mechanical analysis of enhanced geothermal reservoirs with hybrid fracture patterns using a combined XFEM and EDFM-MINC model

The long-term fluid circulation of Enhanced Geothermal Systems (EGS) involves complex coupled Thermal Hydrological-Mechanical (THM) processes dominated by both hydraulic and activated natural fractures. The simulation of such multiphysical processes is known to be a challenging problem. In this work, we have developed an integrated numerical approach to advance the multiphysical simulation of EGS. In our model, the hydraulic fractures of arbitrary shapes in response to pressure changes and thermal strains are handled by the three-dimensional (3D) eXtended Finite Element Method (XFEM). The induced/natural fractures are incorporated into the model and treated as one continuum of the Multiple INteracting Continua (MINC) for the investigation of their impacts. A TOUGH-based program, TOUGH2-EGS, is utilized to simulate the Thermal Hydrological processes. The 3D Embedded Discrete Fracture Method (EDFM), compatible with the 3D XFEM, is adopted to model the hydraulic fractures. Sensitivity analyses are performed for this model in terms of injection temperature and thermal expansivity. The hybrid EDFM and MINC model is established and analyzed for an EGS with both hydraulic and induced/natural fractures. The convergence performance of the single-fracture model shows that an appropriate stiffness coefficient is essential for this model and different choices of the coefficient values result in distinct behaviors. Our numerical results illustrate that the fracture aperture is opened by the cold fluid injection and the reservoir performance is dominated by the thermal stress/strain, which highlights the necessity of multiphysical simulation. In the hybrid model, the enhancement of the natural fracture permeability notably delays the thermal breakthrough by allowing injected fluid to contact more hot reservoirs. The novelties of this work lie in three aspects. Firstly, the combined 3D XFEM and EDFM is capable of handling arbitrary fracture shapes in a 3D EGS model. Secondly, the hybrid hydraulic and induced/natural fracture model enables us to establish the stimulated reservoir volume of the EGS and investigate the operational and geological parameters. Thirdly, the introduction of fracture stiffness stabilizes the numerical solution.