Abstract
It is of great significance for the sustainable development of global energy to develop hot dry rock (HDR) geothermal resources by using enhanced geothermal system (EGS) technology. Different working fluids in EGS have different heat recovery efficiencies. Therefore, this paper takes water and CO2 as the heat-carrying media and establishes a thermal hydraulic mechanical coupling model to simulate the heat recovery process in high-temperature rock mass. By considering the different confining pressures, rock temperature, and injection pressure, the advantages of H2O-EGS and CO2-EGS are obtained. The results show that with the increase of confining pressure, the heat recovery efficiency of water is significantly higher than that of CO2, but at higher reservoir temperature, CO2 has more advantages as a heat-carrying medium. The net heat extraction rate will increase with the increase of injection pressure, which indicates that the mass flow rate plays a leading role in the heat recovery process and increases the injection pressure of fluid which is more conducive to the thermal recovery of EGS. This study will provide a technical guidance for thermal energy exploitation of hot dry rock under different geological conditions.
1. Introduction
The huge consumption of coal, oil, and other fossil fuels not only leads to energy crisis but also causes significant environmental and climate problems [1, 2]. Therefore, the global efforts are to seek green and sustainable alternative energy. As a kind of low-carbon renewable energy, hot dry rock (HDR) geothermal resources have attracted worldwide attention in recent years [3–5]. In the process of developing hot dry rock, it is necessary to form a fracture network connecting injection well and production well by some means, such as hydraulic fracturing, so as to establish a “heat recovery and electricity generation” closed cycle-enhanced geothermal system, as shown in Figure 1 [6, 7].

Conventional EGS uses water as the heat-carrying medium, and there are a lot of researches on this technology at home and abroad [8–11]; it is found that the heat extraction of H2O-EGS has great advantages in the presence of in situ fluid in the reservoir [12]. However, some scholars pointed out that the use of H2O-EGS will cause a lot of water resource loss [13] and also cause a series of problems that are not conducive to the operation and maintenance of geothermal system [14–17]. To solve these problems, in 2000, Brown first proposed that it can use supercritical CO2 (SCCO2) instead of water as a heat-carrying medium [18]. Since then, many studies have been carried out on the feasibility of using CO2 as heat extraction fluid [12, 14, 19–23]. In addition, many studies on the factors affecting the thermal extraction of CO2-EGS (such as CO2 injection rate, permeability, purity of CO2, and injection temperature of CO2) have been carried out [24–27]. CO2 has low density and viscosity, which can obtain large mobility ratio and buoyancy and reduce the pumping power in the heat recovery cycle [18, 28]. However, low specific heat capacity of CO2 leads to lower effective energy of heat extraction than that of water, and the mass flow rate of CO2 is easily affected by reservoir temperature and pressure. At present, domestic and foreign experts and scholars have analyzed the influence of reservoir temperature, gas injection pressure, temperature, water injection temperature, and hot dry rock fracture morphology on CO2-EGS and H2O-EGS through theoretical analysis and numerical simulation, but there is little research on the difference of thermal recovery efficiency of water injection and gas injection in high-temperature rock mass under the same conditions. Therefore, this paper will carry out relevant research on this issue.
This study focuses on the influence of the change of reservoir seepage field, reservoir temperature field, and reservoir stress field on the seepage and heat transfer of a heat-carrying medium, establishes a thermal-fluid-solid coupling model, and, respectively, analyzes the influence of injection pressure, sample temperature, and the change of confining pressure on the heat recovery efficiency when H2O and CO2 are used as working fluids. This study will provide technical guidance for improving hot dry rock productivity under different geological conditions.
2. Governing Equations
The coupling relationship among temperature field, seepage field, and stress field in geothermal exploitation is shown in Figure 2. The coupling of the three fields is realized by the seepage movement of the heat-carrying medium, the deformation of the fractured rock mass, and the heat transfer. Therefore, the first step is to determine the governing equations of each physical field.

2.1. Rock Deformation
The constitutive equation based on Biot’s theory [29] is as follows: where is the solid strain tensor component, is the stress tensor component, is the pore fluid pressure, is the normal stress component, is the Biot coefficient of rock, is the shear modulus, and is the Poisson ratio.
According to Equation (1), the elastic constitutive equation of stress component is obtained as follows: where is the volumetric strain, is the elastic modulus, is the coefficient of thermal expansion of rock mass, and is the temperature of rock mass.
Relationship between strain and displacement under small deformation and deformation geometry equation is as follows:
According to Terzaghi’s [30] effective stress principle, the total stress of the reservoir is as follows: where is the total stress, is the effective stress, and is the equivalent tensor.
Without considering inertial force effect, the momentum conservation equation is as follows:
Through combining Equations (2), (3), (4), and (5), the equilibrium differential equation expressed by displacement is obtained as follows: where is the displacement component, is the Poisson ratio, is the volume force component, and is the lame elastic constant.
2.2. Fluid Flow
This study is based on the geothermal reservoir dominated by tight crystalline rocks such as granite. Here, only the fluid flow in fractures is considered. Fluid flow in the fracture is a saturated single-phase flow process, and the mass conservation equation of fluid flow is expressed by the following formula: where is the fluid density, is the source (sink) term of the fluid, is the dynamic viscosity of the fluid, is the fluid pressure, is the porosity, and is the permeability.
Considering the roughness of fracture, the permeability of fluid in fracture is as follows [31]: where is the roughness coefficient of fracture,
Under the influence of pore pressure and thermal stress, the fracture width is as follows [32]: where is the normal stress and is the normal stiffness of the crack.
The change of porosity under the influence of stress field, seepage field, and temperature field is as follows:
2.3. Heat Transfer
The heat transfer process in fractured rock mass can be expressed by the following formula: where is the temperature of the rock sample, is the density of the fluid, is the specific heat of the fluid, is the flow velocity, is the source term, is the heat conduction flux, is the effective thermal conductivity, is the porosity, is the density of the rock sample, is the thermal conductivity of rock sample, and is the thermal conductivity of fluid.
3. Numerical Model and Simulation Scheme
3.1. Model Description
Previous laboratory tests have studied the thermal recovery efficiency of alternating injection of H2O and CO2 into standard granite samples with single fracture; combined with the experimental results of H2O-SCCO2 alternate injection, the model of rock seepage heat extraction under THM coupling is established [33]. Therefore, the size of the model is consistent with the standard rock sample (), and there is a thorough crack in the axial direction of the model; the geometric model is shown in Figure 3. It contains 3324 free triangle units, the minimum unit size is 1 mm, and the maximum unit size is 6 mm. The finite element software COMSOL Multiphysics is used to solve the problem. The solving module includes structural mechanic module, groundwater flow module, and heat transfer module.

3.2. Initial and Boundary Conditions
The boundary conditions of stress field include the following: the upper boundary applies axial force, the right boundary applies lateral force, and the left and lower boundaries apply roller support.
Boundary conditions of temperature field include the following: the upper boundary applies axial force, the right boundary applies lateral force, and the left and lower boundaries apply roller support.
Boundary conditions of seepage field include the following: the upper boundary and the lower boundary are constant pressure boundaries, and the left boundary and right boundary are both zero flux boundaries.
After reading the relevant literature [34, 35], the thermophysical parameters of CO2 and other fluids are determined, and the relevant parameters in the model are shown in Table 1.
3.3. Simulation Scheme and Data Process
The simulation scheme is shown in Table 2.
The measuring points (25 mm and 50 mm) selected in the model are shown in Figure 4.

The heat energy of heat-carrying fluid is the product of instantaneous sensible enthalpy and mass flow rate, and the net heat productivity of the production well is as follows: where is the heat energy extracted by the heat-carrying fluid, is the mass flow rate of the heat-carrying fluid flowing out of the rock mass (kg/s), is the mass flow rate of the heat-carrying fluid flowing into the rock mass (kg/s), is the specific enthalpy of the fluid at the injection end (J/kg), and is the specific enthalpy of the fluid at the outflow end (J/kg).
4. Results and Discussions
4.1. Effect of Confining Pressure
In the simulation of the influencing factors of confining pressure, the initial temperature of the model is 473 K, and the boundary heat source is 0.001 W/mm2. The fluid injection pressure is 2 MPa, the axial pressure is 35 MPa, and the confining pressure is 26 MPa, 28 MPa, 30 MPa, and 35 MPa, respectively. After, respectively, injecting water and CO2 under different confining pressures, the temperature change at the measuring point is shown in Figure 5.

(a)

(b)
The results show that the temperature at the measuring point increases with the increase of confining pressure after water injection (or CO2), and there are three stages of temperature change: initial stable stage (short time), rapid decrease stage, and final stable stage. When the confining pressure is low, the first stage is not obvious. With the increase of confining pressure, the duration of the initial stable stage is prolonged, and the start of the second stage is delayed. In the third stage, the higher the confining pressure, the higher the temperature is when the pressure is stable. Comparing the temperature after water injection and the temperature after CO2 injection, it is found that the temperature at the measuring point after water injection is slightly higher than that after CO2 injection, especially under the confining pressure of 35 MPa, the temperature difference between the two conditions in the third stage at the measuring point is about 10 K.
Figures 6 and 7 show the stress distribution of rock samples after 600 s water injection and 600 s CO2 injection under 26-35 MPa confining pressure, respectively. It can be seen that when the confining pressure increases, the internal stress of the model increases. According to Equation (9), the model fracture opening decreases with the increase of stress, which reduces the permeability and affects the heat recovery efficiency.


Figure 8 shows the change of net heat with time after water and CO2 injection under different confining pressures. It can be seen that the net heat decreases with the increase of confining pressure in the process of water injection and CO2 injection. Increasing confining pressure reduces fluid velocity, prolongs heat exchange time between fluid and rock mass, but also reduces permeability and fluid mass flux in corresponding time; under the combined influence of the two factors, the net heat decreased, indicating that the latter factor had a greater influence. Comparing with the previous experimental research, the results are generally consistent with the experimental results [36].

(a)

(b)
Comparing the net heat results after water injection and CO2 injection, it is found that the net heat extracted after water injection is slightly higher than that of CO2 injection under the same confining pressure. This difference is not obvious under low confining pressure, but with the increase of confining pressure, the heat recovery efficiency of water is significantly higher than that of CO2.
Figures 9 and 10, respectively, show the temperature distribution of rock samples after water injection for 600 s and CO2 injection for 600 s under confining pressure of 26-35 MPa. It can be seen that the heat loss of rock sample under low confining pressure is obviously lower than that under high confining pressure, and the heat loss of rock sample after water injection is slightly lower than that of CO2 injection under the same confining pressure. It can be seen that the thermophysical properties of fluid have an obvious effect on the heat recovery rate. Although the low mass flow rate of CO2 leads to lower thermal recovery under the same conditions, the heat loss of rock samples by CO2 is higher than that by water.


4.2. Effect of Initial Temperature
In the simulation of the influence factors of initial temperature of rock sample, boundary heat source is 0.001 W/mm2, fluid injection pressure is 2 MPa, axial compression is 35 MPa, and the confining pressure is 26 MPa. We simulated the thermal recovery process of water injection and CO2 injection at the initial temperature of 373 K, 403 K, 433 K, and 493 K, respectively. Figure 11 shows the rule of net heat variation with water and CO2 injection time at different temperatures.

(a)

(b)
It can be seen that in the process of water injection heat recovery, the net heat increases with the increase of the initial temperature of the rock sample. In the process of CO2 injection heat recovery, the net heat decreases with the increase of initial temperature of rock sample. Comparing with the previous experimental research, the results are generally consistent with the experimental results [37]. Due to the small size of rock sample, the temperature of fracture surface decreases rapidly after fluid injection, resulting in the rapid decrease of initial net heat. In the actual EGS reservoir, the heat recovery area will be compensated by the nearby high-temperature rock mass. When boundary heat source is added to the model, the net heat curve fluctuates slightly, which is the result of the interaction between heat recovery and heat compensation. Figures 12 and 13 show the stress distribution of rock samples after 600 s water injection and 600 s CO2 injection at the initial temperature of 373 K-493 K, respectively. It can be seen that when the initial temperature increases, the internal stress of the model increases. Generally, the density viscosity ratio of CO2 is larger than that of water, and the sensitivity of the two fluids to temperature and pressure is quite different. This ratio of water is mainly a function of temperature and less affected by pressure. However, for CO2, this ratio is significantly affected by temperature and pressure. When the temperature of rock sample increases from 373 K to 493 K, the density viscosity ratio of water increases by 0.5, and the specific enthalpy increases by 120 kJ/kg, and the net heat of water increases with the increase of rock temperature. However, the above changes in the rock sample temperature decrease the density viscosity ratio of CO2 by 1 and increase the specific enthalpy by only 30 kJ/kg. The former plays a leading role in the influence of the net heat, while the net heat of CO2 decreases with the increase of rock temperature.


Figure 14 shows the temperature distribution characteristics of rock samples after water injection and CO2 injection for 600 s at 373-493 K initial temperature. It can be seen that the heat loss of rock sample after CO2 injection is slightly higher than that of water injection under the same initial temperature of rock sample. This is consistent with the results of the influence of confining pressure. Therefore, the thermal recovery efficiency of CO2 can be higher than that of water by reasonably selecting the injection conditions of CO2.

4.3. Effect of Injection Pressure
In the simulation of factors affecting CO2 injection pressure, the initial temperature of rock sample is 473 K, the boundary heat source is 0.001 W/mm2, axial pressure is 35 MPa, and the confining pressure is 26 MPa. The model simulates the thermal recovery process of CO2 under high injection pressure of 8 MPa, 9 MPa, 10 MPa, and 11 MPa, respectively. Figure 15 shows the rule of net heat variation with CO2 injection time under different injection pressures.

It is found that the net heat extracted by CO2 increases with the increase of injection pressure in the same amount of time. Comparing with the previous experimental research, the results are generally consistent with the experimental results [37]. The increase of gas injection pressure leads to the increase of gas volume flow in unit time and the heat exchange of gas in rock fracture. The increase of gas injection pressure leads to the increase of rock permeability, resulting in the increase of CO2 mass flow rate at the outflow end. At the same time, it also reduces the heat exchange time of CO2 in fracture. Under the combined influence of the two factors, the net heat increases with the increase of injection pressure, which indicates that the mass flow rate plays a leading role.
Figure 16 shows the stress distribution of rock sample after CO2 injection for 600 s under different injection pressures. It can be seen that the increase of CO2 injection pressure has little effect on the internal stress of the model, and the crack opening will not change much. Compared with the effect of CO2 mass flow rate increasing, the effect of stress can be neglected.

The temperature distribution of rock sample after CO2 injection for 600 s under different injection pressures is shown in Figure 17. It can be found that the heat loss of rock sample gradually increases with the increase of injection pressure in a certain period of time. Therefore, increasing the injection pressure of the fluid is beneficial to the thermal recovery of EGS, but suppresses the thermal compensation effect, and significantly reduces the production life of the reservoir.

5. Conclusions
In this paper, a numerical model is established to simulate the influence of confining pressure and initial temperature of rock sample on the thermal recovery rate during water injection and CO2 injection of fractured rock sample, and the influence of CO2 injection pressure on the thermal recovery rate. This paper compares and analyzes the advantages of heat recovery by water injection and CO2 injection and obtains the following main conclusions: (1)After water injection (or CO2), the temperature at the measuring point increases with the increase of confining pressure, and there are three stages of temperature change: initial stable stage, rapid decrease stage, and final stable stage. After thermal recovery by water injection, the temperature of the measuring point is slightly higher than that of CO2 injection, especially under the confining pressure of 35 MPa. The temperature difference of the measuring point between the two conditions in the third stage is about 10 K. The net heat decreases with the increase of confining pressure, and under the same confining pressure, the net heat extracted after water injection is slightly higher than that of CO2 injection. This difference is not obvious under low confining pressure, but with the increase of confining pressure, the thermal recovery efficiency of water is significantly higher than that of CO2(2)The net heat of water increases with the increase of rock sample temperature, while the net heat of CO2 decreases with the increase of rock sample temperature. At the temperature of 393 K rock sample, the maximum net heat ratio of CO2 to water is about 3. With the increase of sample temperature, the maximum net heat ratio of CO2 to water decreases to 0.4 at 493 K. Therefore, with the increase of rock temperature, the heat extraction rate of CO2 is gradually lower than that of water(3)In the same amount of time, the net heat extracted by CO2 increases with the increase of injection pressure. With the increase of injection pressure, the velocity of CO2 in fracture increases, which increases the mass flow rate of CO2, but also reduces the heat exchange time of CO2 in fracture. The net heat increases with the increase of injection pressure, which indicates that the mass flow rate plays a leading role; increasing the injection pressure of the fluid is beneficial to the thermal recovery of EGS, but suppresses the thermal compensation effect and significantly reduces the production life of the reservoir
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the Fundamental Research Funds for the Central Universities, China University of Mining and Technology (grant number 2015XKZD06).