Comparison of two approaches for modelling fracture opening due to cold water injection in geothermal reservoir
Authors: Luo Wen , Josselin Ouf , Vardon Philip J. , Dieudonné Anne-Catherine , Liaudat Joaquín , Khaledi Kavan , Jalali Reza , Amann Florian 
Published: October 2023

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 [1]. 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 (Fig.1 (a)). 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. Then, an equivalent continuum permeability  [m2] is obtained as follows:

where  [m] is the elastic opening, br and bmax [m] are the residual and maximum apertures, respectively, a is the stress dependency coefficient;  [Pa] is effective normal stress, and [m]the fracture spacing. The combination of the continuum elements to model an implicit fracture embedded in a fracture zone and the evolving permeability has been used successfully used to model elastic opening in volcanic rock [2] and granodiorite rock [3, 4].

In the second approach, zero-thickness interface elements are used to explicitly model the opening of an existing fracture (Fig.1 (b)), following the work by Liaudat et al. [5]. 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 this approach, the cubic law is used to estimate the longitudinal transmissivity  [m3] as a function of the normal fracture aperture:

where   [m]  is the normal separation and  [m3] is a constant value which makes it possible to assign an initial longitudinal transmissivity even if the fracture is mechanically closed.

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.

In order to compare both approaches, a simple synthetic model is studied. The model with its boundary and initial conditions are shown in Fig.2. Water with a temperature of 20 °C and pressure varying from 4 MPa – 10 MPa is injected at the centering point of the fracture. Figure 3 shows the injection pressure profile. Fracture permeability evolution  is compared in Fig.4, which shows the two different methods capture different characteristics of the natural fracture behaviour. In the continuum method, the fracture is modelled as a continuum element, thus any pressure or stress change will result in corresponding changes in the fracture permeability, which is reflected in Fig.4 showing that the permeability changes with the injection pressure changing, In contrast, in the interface element, Fig.4 shows only when the injection pressure reaches 10 MPa, the permeability increases sharply from the initial value to the peak value. This is because penalty method is used to allow the small overlapping of the two faces of the discontinuity, thus keeping negligible permeability variation, before it opens. The result in Fig.4 also illustrates the effect of thermal stress on the fracture opening, considering the total vertical stress is 10 MPa while injection pressure is below or equal to 10 MPa.

In conclusion, though the peak permeability predicted from two methods are consistent, both around 1.5 x 10 -12 D, the behaviour before fracture opens in these two approaches are quite different. For the interface method, permeability variations keep negligible before fracture opens, while for the continuum method, permeability changes with corresponding pressure or stress changes without dependence on the opening of the fracture.