Comparison of two approaches for modelling fracture opening due to cold water injection in geothermal reservoir

Abstract

Heat and electricity production from deep hot reservoirs through hydrothermal or petrothermal system requires to inject cold fluid in a naturally or artificially fractured medium. Cold water injection in a hot reservoirs causes thermo-hydro-mechanical (THM) coupled processes that may have several influences on operations and long term production. Cold water has a higher viscosity, which means injection can become difficult due to higher flow impedance, while thermal diffusion in the rock matrix causes shrinkage and thus an increasing fracture aperture over time that can lead to flow channeling and a reduction of operation time. A deep understanding of the contribution of coupled THM processes to injection into fractured media is thus important to predict the long-term performance of a geothermal power plant. In this study, coupled THM processes in a single fracture are numerically investigated with two different approaches to model the discontinuity. The first numerical approach proposes modelling an implicit fracture in a fracture zone using solid elements. The impact of the fracture is introduced by means of a stress dependent fracture aperture, which in turn determines the fracture transmissivity. This function allows recovery of the opening, and thus the model replicates opening and closure of the fracture. In the second approach, zero-thickness interface elements are used to explicitly model the opening of an existing fracture, following the work by Liaudat et al. A constitutive law that is capable of describing the fracture initiation, propagation, closing and opening is implemented in this approach, considering heat transfer in the discontinuity. In both approaches, the simulation of cold water injection requires to consider the heat and hydraulic flows in the fracture alongside mechanical behaviour, since changes in pore water pressure and temperature influence the fracture aperture, thus modifying the fracture transmissivity.