Abstract

Liquid nitrogen (LN2), which can greatly improve the efficiency of hot dry rock (HDR) mining, is commonly used as a cooling material in the enhanced geothermal system (EGS). Physical property, triaxial compression, and permeability tests were undertaken on treated granite samples, for a better scientific understanding of the effect of the LN2 cooling method on the mechanical and permeability properties of the rocks after heat treatment. The experimental results indicated that the physical properties of the treated granite change significantly, such as the density and wave velocity are substantially reduced. Meanwhile, with the increase of treatment temperature, the macroscopic cracks on its surface are gradually generated and the volume is expanded clearly. In addition, the surface wettability of granite gradually increases with increasing temperature. Compared with the air/water cooling methods, under LN2 cooling condition, the mechanical properties decrease markedly. When the temperature exceeds 600°C, the granite strength decreases significantly to only 56.16% of the reference value. The deformation properties also change significantly, with a final strain of about 3% at failure for a sample at 800°C, showing an obvious ductile deformation characteristic. Further, an appreciable correlation also exists between the initial permeability of granite and temperature. Once the temperature exceeds 200°C, the increase in temperature contributes to the increase in initial permeability. In addition to the effect of temperature, the increase in load also leads to a change in the permeability coefficient. When the temperature reaches 600°C, the permeability of granite first decreases and then increases with the increases in axial stress. The results of this paper are valuable in understanding the effect of thermal shock by LN2 on the fracturing efficiency and permeability characteristics of dry hot rocks.

1. Introduction

Development of hot dry rock (HDR) resources plays an important role in meeting the demand for energy [15]. Generally, the geothermal energy is mainly stored in a hot granite reservoir, which is distributed at the depth of 2–6 km with a temperature above 150°C [6]. Due to the low permeability of reservoir matrices, the permeability enhancement technologies are often suggested to improve the thermodynamic efficiency and productivity of the enhanced geothermal system [7, 8]. Liquid nitrogen (LN2) cryogenic fracturing has been proven to be an efficient stimulation method, which has a great potential in establishing pathways with high permeability [9, 10]. During the process of cryogenic fracturing, the use of LN2 can produce a high thermal gradient and a rapid temperature change in the rock inner and outer body, leading to a greater thermal stress and more thermal cracks in HDRs compared with traditional water/air cooling [11]. Water was often used as a coolant in previous works, but water has the disadvantage that can even cause important reductions in the mechanical properties of nonheat-treated rocks [12, 13]. Consequently, the physicomechanical properties of the HDR reservoirs are greatly influenced by LN2 cooling treatments. The effects of different thermal shocks by water cooling on physical properties, mechanical behaviors, and permeability parameters of rocks are critical knowledge for the successful reformation of deep geothermal reservoirs [14].

In recent years, plentiful in-depth studies on the physical and mechanical properties of granites after high-temperature treatment have been conducted, and many valuable conclusions have been obtained by various investigators. Previous experimental studies have indicated that thermal-induced microcracks are formed in granite after high-temperature treatment, which was caused by thermal expansion and uneven temperature distribution of minerals [13, 14]. Once the critical temperature is reached, macroscopic cracks developed inside the rock by expansion of mineral particles further contributing to the high-temperature erosion of granite [15, 16]. Zhang et al. [17] analyzed in detail the effect of temperature variations on the distribution of several variables in the physical field and the cracking by an improved hydrothermal-salt-mechanical (IHTSM) model. Shao et al. [18, 19] found that rapid cooling would instigate more thermal cracking within the granite, leading to more pronounced deterioration of mechanical properties. Subsequently, they considered the potential causes of the thermal sensitivity discrepancy among the three rocks based on the result of scanning electron microscope (SEM) and X-ray diffraction (XRD) [8]. However, Wu et al. [20] indicated that the uniaxial compressive strength (UCS) of granites cooled in water did not drop a lot until the rock temperature reached to 400°C. At higher cooling rates, the strength and elasticity of the granite might be further reduced under strong thermal shock. Chen et al. [21] indicated that the peak stress of heated granite decreases with the heating temperature increases (from 200°C to 1,000°C). Brotóns et al. [22] found that the thermal stress produced by the temperature gradient can degrade physical properties of the stone as the temperature increases (from 105°C to 600°C). Liu and Xu [23] investigated physical and mechanical properties of granite and sandstone after high-temperature treatment (from 100°C to 1,000°C) and found that 400°C is the sensitivity point for strength.

Recently, researchers have shown an increased interest in the relationship between high-temperature treatment and permeability, because the permeability of granite may increase in thermal-induced microcrack formation and propagation. Chen et al. [24] carried out compression tests and found the change in permeability was negligible before microfracture coalescence, and when the crack was in the growth region, it increased dramatically. Subsequently, to understand the effect of temperature on the permeability properties of granite, Chen et al. [21, 25] found that the thermal sensitivity temperature of permeability is 500°C. Then, Tian et al. [26] found a significant effect on initial and residual permeabilities of granite with increases of 3-4 orders of magnitude under temperatures from 25°C to 750°C by air cooling.

However, in contrast to the study of the mechanical properties of rock, there is much less information about the effects of LN2 cooling on permeability characteristics. Despite the importance of permeability for the extraction of HDR, there remains a paucity of evidence on the permeability evolution pattern of granites after high-temperature heating and LN2 cooling treatment. Thus, the aim of this essay is to explore the effect of different heating temperatures on the mechanical and permeability properties of granite under the action of LN2 cooling. The granite samples were subjected to heating-cooling in four levels, i.e., a slow heating to 200°C, 400°C, 600°C, and 800°C followed by a rapid cooling with LN2. A comparative analysis was conducted to study the variations in the physicomechanical and permeability behavior of granite subjected to different heating-cooling thermal shocks.

2. Materials and Methods

2.1. Sample Preparation

Granite, as a typical rock type in HDR [17], was selected as experimental materials in our study. The granite specimens were collected from an outcrop located at Hunan Province, China. According to the mineralogical analyses using X-ray diffraction (XRD), the granite is primarily composed of feldspar (47.6%), quartz (15.7%), and biotite (23.8%). The rock samples were shaped into the standard cylinders of in the laboratory according to the International Society for Rock Mechanics (ISRM) standards [27]. Figure 1(a) shows the rock samples and test setup used in our experiment. The thermal loading schedule is shown in Figure 1(b). In the first step, 12 granite samples were divided into four groups, which were all first heated to different target temperatures (200°C, 400°C, 600°C, and 800°C) using a high-temperature furnace with a heating rate of 2°C/min and then were kept at the preset temperature for 2 h. Secondly, the heated specimens were taken out of the high-temperature furnace and placed in a tank with 25 L of liquid nitrogen. Finally, granite specimens with treatment completed were installed on the designated device, as shown on the right side of Figure 1(a).

2.2. Testing Methods
2.2.1. Physical Parameter Test Procedure

Prior to mechanical testing, the average and standard deviations of physical properties of the sample, such as density, size (i.e., height and diameter), and P-wave velocity, were first measured before heating and after cooling. The contact angle of droplets is an important indicator of rock wettability which reflects the hydrophilic properties of the rock, as shown in Figure 2(a). Wanniarachchi [28] provided the wettability test method which is using the high-definition image technique to obtain the contact angle of water droplets on the granite surface and the method schedule as follows. Firstly, the syringe was cleaned using distilled water and filled with water. Secondly, a clean surface specimen was placed stably on a horizontal table. Subsequently, a water droplet was slowly injected by the syringe on the surface of the granite, as shown in Figure 2(a). Finally, a high-definition camera was used to record interfacial surface tension around the droplet. The contact angle of the rock sample surface is the angle between the contour line of water droplet with the horizontal of rock surface, as shown in Figure 2(b).

2.2.2. Mechanical Test Procedure

Triaxial compression tests were performed using the MTS815.04 rock mechanics test system on the thermal shock-treated granite specimens, as shown in Figure 3. The system can perform triaxial deformation experiments at confining stress up to 140 MPa, axial load up to 4600 kN, and hydraulic pressure up to 140 MPa and can measure the permeability of rock samples down to . This study used deionized water as the pore fluid and silicone oil as the confining pressure fluid. The confining and inlet pressures were controlled by an oil pressure intensifier and a gas pressure pump, respectively. The triaxial compression strength was obtained with a confining pressure of 5 MPa, a confining pressure rate of 0.1 MPa/s, and an axial displacement loading rate of 0.001 mm/s.

2.2.3. Permeability Test Procedure

Considering that granite has ultralow permeability, the transient-pulse permeability test method was used to guarantee the measurement accuracy [29], as shown in Figure 4. To investigate the permeability evolution of thermal shock-treated samples with increasing deviatoric stress, the following test procedures were designed: (i)The treated specimens were saturated in water for 48 h and then jacketed with 0.5 mm thick Teflon heat shrink tubing, sealed at both ends with several stainless-steel wires(ii)The assembled jacketed sample was then placed into a triaxial pressure cell with the hydrostatic pressure set at 5 MPa(iii)The circumferential pressure was held at 5 MPa and the axial stress was increased by a displacement-controlled manner with a deformation rate of 0.001 mm/s(iv)From this stage, the sample permeability was measured until the peak stress was reached. For each permeability test, 3 MPa of hydraulic pressure was applied at both the top and bottom ends(v)The hydraulic pressure at the top surface was decreased to 1 MPa in order to create a pressure gradient between the two end surfaces. Water flowed from the bottom end of the sample to the top end and through the interior of the sample

It is important to note that this water pressure should be lower than the confining pressure to avoid lateral leakage.

The - curves recorded were used to calculate the sample permeability values according to Equation (1). It was noted that the limiting load and axial load were kept constant while the permeability measurements were performed. The permeation test was terminated when the pressure gradient reached or approached 0.4 MPa, the permeation time reached 3 h, or the permeation pressure difference stabilized. In the permeability test measured by the transient method, the water pressure difference () between the two ends of the specimen decreased gradually, and the rate of decrease depended on the rock type, rock configuration, specimen length, specimen cross-sectional area, water density and viscosity, and stress state [29].

where is the permeability of the rock sample, is the compressibility of the water, is the volume of the upper and lower reservoirs, is the length of the specimen, is the viscosity of the water, is the duration of the test, is the cross-sectional area of the specimen, and is the difference in water pressure between the upper and lower ends of the rock specimen. Well-known parameter values are , , , , , , and .

3. Results and Discussion

3.1. Physical Properties
3.1.1. Surface Morphology, Density, and P-Wave Velocity

The surface morphology graphs of granite samples subjected to different thermal shock treatments (200, 400, 600, and 800°C) are shown in Figure 5. According to previous investigations, the escape of adsorbed water, crystalline water, and tectonic water existed in the rocks during the high-temperature treatment [30]. The gravity and hydroscopic water escaped completely at about 150°C and 200°C, respectively [31]. As a result, as shown in Figure 5, there is no visible difference between the surface colors of granite samples at 25°C (with no thermal shock treatment) and those of granite samples at 200°C. With further increase in temperature, when at 600°C, the quartz experienced an α-β phase transition, which caused a significant fading of the surface color [32, 33].

The basic physical and mechanical properties of granite samples subjected to LN2 cooling under different temperatures are summarized in Table 1. The relationships between P-wave velocity and density and temperature for granite samples heated to different temperatures followed by LN2-coling treatment are depicted in Figure 6. It can be seen from Figure 6 that the variation in P-wave velocity exhibits a noteworthy downward trend with a relatively small initial reduction at 200°C and then an approximately linear rapid decrease with further rise of temperature. However, the density apparently does not change before 600° and is reduced by approximately 20% at 800°C. These results indicate that the P-wave velocity, compared with the density, is more sensitive and more suitable for characterizing the damage caused by thermal shock to granite.

3.1.2. Wettability

The results of wettability quantified by contact angle for different thermal shock-treated granite samples are shown in Figure 7. It can be found that the contact angles of the water droplet are increased with the heating temperature, and the angles are 63° at 25°C, 147° at 200°C, 165° at 400°C, 170° at 600°C, and 178° at 800°C (Figures 7(a)7(f)). The results suggest that the thermal shock-treated granite surface is more hydrophilic than the untreated granite surface, and water is more readily absorbed on the treated granite surface. It is noteworthy that the surface of the granite after the 800°C treatment has obvious macroscopic cracks and a large number of fine cracks. This allowed the water droplets to be immediately absorbed into the interior as soon as they made contact with the rock surface. The water stains produced by water penetration can be seen in the top view, as shown in Figure 7(f).

3.2. Mechanical Behavior
3.2.1. Stress-Strain Curves

The triaxial compressive stress-strain curve of the granite sample with no thermal treatment at room temperature (25°C) under the confinement pressure of 5 MPa is illustrated in Figure 8. It can be seen that the curve can be divided into four stages, i.e., (i) compaction stage (O-A), (ii) elastic stage (A-B), (iii) yield stage (B-C), and (iv) failure stage (after the C point) [34]. The characteristic stresses corresponding to each stage in the triaxial compression curve are point A (crack closure stress), point B (damage threshold stress), and point C (peak stress). Each stress levels are selected for permeation tests to observe the permeation characteristics of granite under different stages, and the experimental results will be presented in the next section.

The stress-strain curves under hydronic-mechanical coupling condition of granite samples after exposure to different thermal shock treatments are plotted in Figure 9. It can be seen that the stress-strain responses of five specimens have obvious differences, suggesting that thermal shock has a significant effect on the mechanical behavior. With the increase in heating temperature, the initial nonlinear deformation segment induced by the closure of preexisting microcracks becomes more and more pronounced, especially if the temperature is over 600°C. Then, in the linear growth segment, as the temperature increases, the slope of the curve becomes lower and lower. Meanwhile, the decline of the curve in the postpeak segment gradually changes from a sudden drop to a slow drop, especially at 800°C, which shows a clear transformation from brittleness to ductility. These changes all can be attributed to more thermal cracks occurring in the sample led by the rapid release of thermal stress under the action of extreme temperature difference.

3.2.2. Triaxial Compressive Strength

The changes of the triaxial compressive strength of granite samples after heating and LN2 cooling treatment under a confining pressure of 5 MPa are shown in Figure 10. It can be seen from Figure 10 that the triaxial compressive strength of granite does not change significantly during temperatures of 25~400°C. This phenomenon indicates that when the heating temperature is relatively low, only physical morphology variation and slight microcrack initiation occurred in granite with the escape of absorbed, bound, and crystalline water [14]. As the temperature increases, further and exceeding 573 degrees, more thermal-induced microcracking transferred into macrocracks as a result of the phase change of quartz [30]. Due to the formation of macrocracks, thermal stresses were released, and the continuous increase in temperature promoted the formation of more microcracks [33]. In particular, when the temperature reached 800°C, a large amount of crystal and structural water escaped, and Mg2+ oxidation in the black mica led to a large amount of thermal change, which eventually induced macroscopic crack formation. Therefore, the peak strength remained relatively constant at the temperature of less than 600°C. When the temperature reached 600°C, some metal bonds (such as Al-O, K-O, Na-O, and Ca-O) broke, and some minerals (such as calcium carbonate and calcium montmorillonite) melted or decomposed [35]. More cracks were induced by heating, and the peak strength of granite decreased significantly at 800°C.

To further comparatively analyze the effect of the cooling method on the fracturing efficiency of dry hot rock, the relationships between triaxial compressive strength and temperature under air and water cooling conditions reported in previous studies are also plotted together in Figure 10. It can be found that the variation trend of granite strength with temperature can be divided into three stages, which are (I) 25°C-400°C, (II) 400°C-600°C, and (III) above 600°C. Under either air or water cooling condition, the granite strength shows a trend of fluctuation with the change of temperature [35, 36], and there is a critical temperature of transition from strengthening to weakening. Under the LN2 cooling condition, however, the strength of granite showed a sustained decreasing trend with increasing temperature, which conformed to the quadratic function. When the temperature reached 800°C, the granite strength decreases significantly and reaches 80.29 MPa. These differences in results indicate that at a certain heating temperature range, the relatively low cooling rate is more likely to lead to the closure or healing of preexisting microcracks than to the expansion and initiation of cracks.

3.2.3. Failure Pattern

Final failure patterns of granite samples subjected to different thermal shock treatments under triaxial compression are presented in Figure 11. When the temperature is at 25°C, the angle of the main macroscopic fracture with respect to the horizontal direction is about 60°, as shown in Figures 11(a) and 11(b). When the temperature rises to 200°C, 400°C, and 600°C, the angle of the macroscopic fracture is dominated by large-angle tensile fractures of 70°-80°, as shown in Figures 11(c)–11(e). The difference of macrocracks is owing to the temperature field heterogeneous distribution and the thermal stress in minerals (tensile) under rapid quenching. At failure, there are many large-angle macroscopic tensile fractures, which is due to the decrease of tensile resistance with the increase of temperature [37]. As a result, the granite transforms from brittle to ductile deformation.

3.3. Permeability Parameters

During the process of hydraulic-mechanical coupling loading, the pores and microcracks in the treated granite can produce some complex dynamic volume changes, resulting in the permeability evolution characteristics of granite. The time taken for the permeability pressure difference () to reach stability and to decrease below 0.5 MPa at the samewhich is related to the stress state of the granite can be seen in Figure 12(a). In addition, it is noteworthy that with the time increase, the initialdecreases gradually and cannot reach the preset 2 MPa as presented in Figure 12(a), because of the development of internal cracks in the granite with loading increase. In particular, at the yield stage, the initial is below to 0.5 MPa, as shown in Figure 12(a).

After the logarithm fetch on the , the typical evolution law of the logarithm of permeability pressure with loading time under percolation is shown in Figure 12(b). The typical example of experimental data are linearly fitted to obtain slope , fitted correlation coefficient , , , , , and , which are finally calculated based on Equation (1): .

To examine the permeability evolution during triaxial compression, the curves of deviatoric stress and sample permeability as a function of axial strain for each heating temperature are plotted in Figure 13. The granite samples under different thermal shock treatments exhibit a similar mechanical behavior: (1) A certain degree of stress recession occurs during each percolation, which is caused primarily by grain slipping during percolation. (2) In the process of percolation, the fluid under the action of passed through the crack channel which may cause further damage to rock and also reduce the stress.

As shown in Figures 13(a)13(c), the permeability evolutions of granite samples at 200°C and 400°C are similar to the reference sample (25°C). Before reaching 0.8% of the axial strain, the permeability is relatively stable and floats around . It is due to the fact that the crack channels are not fully formed during this stage. When the axial strain continues to increase and the axial stress approaches the peak value of stress, the permeability of the sample increases significantly by nearly two orders of magnitude to , manifesting that liquid seepage channels have been fully developed inside the sample. In addition, the permeability of samples exhibits an interesting phenomenon at 600°C, as shown in Figure 13(d). The permeability decreased first and then increased with the increase of the axial strain, and the turning point of the axial strain is 0.8%. It is due to the closure of a large number of thermally induced cracks during the initial stage, leading to a decrease in permeability reaching . Thereafter, the incremental generation of new cracks induced by compression results in the increase of permeability. Moreover, as exhibited in Figure 13(e), the permeability evolution pattern of the granite after being treated at 800°C shows a distinct difference with other temperature-treated samples. The permeability changes in all loading stages are relatively slight and always fluctuate around , which is significantly larger than the permeability at low heating temperature. This can be contributed to a large number of macroscopic thermal damage cracks existing inside the granite, and the thermal-induced cracks are hard to close even during the action of compression.

To further analyze the effect of treatment temperature on permeability characteristics under different stages, typical permeability data (average values) were selected to obtain the correlation between permeability and temperature, as shown in Figure 14. It can be seen that the permeability of thermal-treated granite during the elastic stage is positively correlated with temperature. In addition, the permeability coefficient in the elastic stage is lowest for each treated granite because of the initial crack closure caused by compression. Thereafter, with the increase of stress, the initiation and coalescence of microfissures and then transformation into macroscopic fractures provide a rapid flow passage for fluids.

Moreover, there is a significant quadratic relationship between temperature and permeability, as illustrated in Figure 14. In the compression-density stage, temperature is greatly affected the permeability of the sample. Permeability of the rock sample increased with the change of temperature in this stage. A similar characteristic has been exhibited in the elastic stage, which is weakly influenced by temperature at 200°C. During the yield stage, there is a phenomenon that the permeability variations show decrease at first and then increase with the increase of temperature. The permeability at the compression stage reached the minimum when the granite is exposed to temperatures between 200°C and 400°C. At the compression-density stage, the permeability is lowest when the rock was under thermal treatment. When the treatment temperature is below 400°C, flow channels can remain open and improve the permeability of the rock which is attributed to the loss of absorbed and bound water [38]. Furthermore, the loss of mineral water may damage the mineral skeleton, and this mineral damage will cause the rock to be denser and the brittleness of the rock increases, which leads to the additional improvements of permeability [39, 40]. When the temperature exceeds 400°C, the stress state takes effect on the compacting of the initial fracture, resulting in a decrease in permeability. When the treatment temperature is exactly 400°C, the change of permeability with the stress state is not obvious indicating that the 400°C is a turning point.

4. Conclusion

Rapid thermal shock can greatly improve the efficiency of HDR fracturing, especially with the use of LN2. In this study, the influence of different thermal shocks (slow heating to different temperatures followed by rapid cooling with LN2) on the physical, mechanical, and permeability properties of granite samples was systematically studied. Nondestructive tests (surface morphology, density, P-wave velocity, and wettability) and destructive tests (triaxial compression and seepage) were adopted to investigate these properties. Furthermore, based on the experimental results, a comparative analysis was conducted to study the variations in the physicomechanical and permeability behavior of granites after different thermal shock treatments. The main conclusions are summarized as follows. (1)The LN2, an efficient coolant, can change the physical properties of HDR by rapid cooling. Notably, the P-wave velocity of granite is reduced to nearly 80%, when the treatment temperature is at 800°C. Additionally, with the increase of temperature, a large number of thermal damage cracks are produced inside granite specimens, and the wettability of granite is enhanced. The fluid directly penetrates the granite interior and cannot form a complete water droplet under high-temperature treatment, which leads to the granite wettability being stronger, especially at 800°C treatments(2)The effect of LN2 cooling on the mechanical property of granite is significant. When the temperature exceeds 600°C, the strength of granite exhibits a distinct decrease, only 56.16% of the reference sample. Moreover, the deformation properties also change apparently; especially at 800°C, the final strain at the time of failure reaches about 3%, showing an obvious ductile deformation characteristic(3)A correlation also exists between the initial permeability of granite and temperature. The permeability evolution pattern of granite after liquid nitrogen treatment shows a quadratic relationship with temperature. When the treatment temperature is lower than 400°C, the increase of stress enhances the impermeability of granite and makes the permeability coefficient decrease. When the temperature reaches 600°C, the permeability of granite first decreases and then increases with the increase of axial stress. With the temperature of treatment further increased, a large number of connected fractures are sprouted inside the rock, leading to an abnormally rapid seepage process and a larger calculated permeability coefficient

Data Availability

The data used to support the findings of this study are included in the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to acknowledge the financial support for this study supported by the Second Tibetan Plateau Scientific Expedition and Research Program (grant No. 2019QZKK0904), the National Natural Science Foundation of China (grant Nos. 42007254 and 41831290), and the Natural Science Foundation of Zhejiang Province (grant No. LQ20E080006).