A Thermal Effect Model for the Impact of Vertical Groundwater Migration on Temperature Distribution of Layered Rock Mass and Its Application

: On the basis of the one-dimensional heat conduction–convection equation, a thermal effect model for vertical groundwater migration in the stratiﬁed rock mass was established, the equations for temperature distribution in layered strata were deduced, and the impacts of the vertical seepage velocity of groundwater and the thermal conductivity of surrounding rocks on the temperature ﬁeld distribution in layered strata were analyzed. The proposed model was employed to identify the thermal convection and conduction regions at two temperature-measuring boreholes in coal mines, and the vertical migration velocity of groundwater was obtained through reverse calculation. The results show that the vertical temperature distribution of the layered rock mass is subject to the migration of the geothermal water; the temperature curve of the layered formation is convex when the geothermal water travels upward, but concave when the water moves downward. The temperature distribution in the stratiﬁed rock mass is also subject to the thermal conductivity of the rock mass; greater thermal conductivity of the rock mass leads to a larger temperature difference among regions of the rock mass, while weaker thermal conductivity results in a smaller temperature difference. A greater velocity of the vertical migration of geothermal water within the surrounding rock leads to a larger curvature of the temperature curve. The model was applied to a study case, which showed that the model could appropriately describe the variation pattern of the ground temperature in the stratiﬁed rock mass, and a comparison between the modeling result and the measured ground temperature distribution revealed a high goodness of ﬁt of the model with the actual situation.


Introduction
Geothermal energy, a renewable energy, has shown great application prospects, and geothermal water is currently the most widely used form of geothermal energy around the world [1][2][3][4]. Coal, a reliable source of energy for many countries including China, accounts for 40% of all electric power generated worldwide. As the exploration and the development of coal intensify, more medium-and low-temperature geothermal water resources have been found in coal measure strata, which provides a basis for further investigation and evaluation of geothermal water resources. In coalfield exploration, temperature measurement is generally carried out along with geological drilling. Migration of geothermal water will affect the temperature of the surrounding rock mass; thus, in areas that witness seepage of geothermal water, the local temperature will differ substantially from that in other areas of the rock mass due to the heat transfer from the water. In this case, the temperature curve can be used to identify where the geothermal water occurs, as well as the temperature and seepage conditions of the water; thus, it can be used to lay a theoretical foundation for further analysis of geothermal water resources. two temperature-measuring boreholes, as well as to reverse-calculate the vertical migration velocity of groundwater. The research result is expected to provide a theoretical basis for future development of geothermal water resources.

One-Dimensional Geothermal Field of Homogeneous Rock Mass
In the present work, it is assumed that the underground rock mass is composed of an ideal water-resisting layer free from water penetration. The rock strata are all isotropic and homogeneous conductors of heat, the total thickness of the rock mass is h, and the ground temperature distribution of surrounding rocks is a thermal conduction temperature field. When the heat transfer is stable, the vertical temperature distribution can be expressed by a one-dimensional mathematical model, as shown below.
Then, the ground temperature at a depth of z is T = T 0 + T 1 −T 0 h z, where T 0 and T 1 represent the known temperatures of the upper and lower boundary of the rock stratum, respectively.
In areas where movements of groundwater are observed, the temperature field is affected by both conduction and convection of heat. When the groundwater moves fast in a geological structure, there is a temperature gap between the water-moving area and the surrounding rock mass, leading to heat exchange in between. In general, a onedimensional mathematical model can be established to describe the distribution of the ground temperature T at any depth z [5,18]. The model is shown below.
where a is heat exchange coefficient between the groundwater and the surrounding rock mass, k is thermal conductivity of the rock (W·(m· • C) −1 ), c w and ρ w are the specific heat capacity (J/(kg· • C)) and density (kg/m 3 ) of the groundwater, respectively, c 1 and ρ 1 are the specific heat capacity (J/(kg· • C)) and density (kg/m 3 ) of the groundwater-rock mixture, respectively, n is porosity of the surrounding rock mass, T c is the temperature of surrounding rock mass ( • C), and t is time (s).
In the actual strata, the ground water moves vertically at a small rate. During the movement, the water is in full contact with the surrounding rock, and the water flow at higher temperatures from deep down slowly warms up the surrounding rock. In this case, the temperature difference between the water flow and the surrounding rock will gradually decrease, and the heat exchange between the two can be ignored at the end. Therefore, we assume a = 0, and Equation (2) can be simplified as follows: When the flow migration time t → ∞ , a thermal balance occurs between the water flow and the surrounding rock, which is considered as a thermal steady state. Thus, Equation (3) can be transformed into Here, the surrounding rock mass is assumed to be an isotropic and homogeneous heat conductor. T 0 and T 1 are known temperatures of the upper and lower boundaries of the rock, respectively. Then, the analytical solution of Equation (4) can be expressed as where β 1 = nc w ρ w v z (1−n)k .

One-Dimensional Geothermal Field of Stratified Rock Mass
According to the one-dimensional mathematical model for the temperature distribution of thermal conduction and conduction-convection of the homogeneous rock mass mentioned above, the thermal effect model of vertical fluid migration of the stratified rock mass is established. There are two types of temperature fields in the model, namely, the thermal conduction (Tc-type) and thermal conduction-convection (Tcc-type) temperature fields. The model consists of n layers, in which the upper layer m is a Tc-type temperature field, the middle layer j is a Tcc-type temperature field, and the lower layer n − m + j is a Tc-type temperature field ( Figure 1). It is assumed that the contact boundaries of adjacent rock strata are S 1 , S 2 , . . . . . . , S n−1 , with temperatures of T 1 , T 2 , . . . . . . , T m (interfacial temperatures within the m layer), T m+1 , T m+2 , . . . . . . , T m+j (interfacial temperatures in the j layer), and T m+j+1 , T m+j+2 , . . . . . . , T n−1 (interfacial temperatures within the n − m + j layer). The temperature of the interface between two adjacent temperature fields (z = h 1 + h 2 + . . . . . . + h n−1 ) is the same, set to T n−1 . The temperature T 0 at z = 0 and the temperature T n at z = h 1 + h 2 + . . . . . . + h n are known. The vertical temperature distribution in the layer i Here, the surrounding rock mass is assumed to be an isotropic and homogeneous heat conductor. T0 and T1 are known temperatures of the upper and lower boundaries of the rock, respectively. Then, the analytical solution of Equation (4) can be expressed as where ( )

One-Dimensional Geothermal Field of Stratified Rock Mass
According to the one-dimensional mathematical model for the temperature distribution of thermal conduction and conduction-convection of the homogeneous rock mass mentioned above, the thermal effect model of vertical fluid migration of the stratified rock mass is established. There are two types of temperature fields in the model, namely, the thermal conduction (Tc-type) and thermal conduction-convection (Tcc-type) temperature fields. The model consists of n layers, in which the upper layer m is a Tc-type temperature field, the middle layer j is a Tcc-type temperature field, and the lower layer n − m + j is a Tc-type temperature field ( Figure 1). It is assumed that the contact boundaries of adjacent rock strata are S1, S2, ……, Sn−1, with temperatures of T1, T2, ……Tm (interfacial temperatures within the m layer), Tm+1, Tm+2, ……Tm+j (interfacial temperatures in the j layer), and Tm+j+1, Tm+j+2, ……Tn−1 (interfacial temperatures within the n − m + j layer). The temperature of the interface between two adjacent temperature fields (z = h1 + h2 + …… + hn−1) is the same, set to Tn−1. The temperature T0 at z = 0 and the temperature Tn at z = h1 + h2 + …… + hn are known. The vertical temperature distribution in the layer i (1 ≤ i ≤ n) is derived below.  The "iterative" and the "chasing" derivation methods are often used for the establishment of models. The "chasing method" is to establish an iterative equation between the temperature T i−1 and T i of the two adjacent strata and the known temperature T 0 of the upper surface; the iteration continues until i = n. According to the obtained iterative equation, the known lower surface temperature T n was used for substitution in the equation to complete the "chasing" process, after which the relation equation for the temperature of one layer in the model (T i ) expressed by the known T 0 and T n was obtained. The results are described below. (1) The temperature of each layer in the Tc-type temperature field of the upper m layer is expressed as The temperature gradient is The boundary between the i − 1 layer and the i layer is the contact surface of two solids. According to the law of heat conduction, the heat flux density of the fluid flowing through the contact surface of two objects is equal, i.e., It can be concluded that the temperature of the contact surface of any two adjacent layers is where (2) The temperature of each layer in the Tcc-type temperature field in the middle layer j is expressed as follows: .
The geothermal gradient is As the temperature T m of the interface between the m layer and m + 1 layer is the temperature of the interface between the Tc-type and Tcc-type temperature fields, the temperatures of these two layers need to be calculated separately. On the contact surface between the m layer and the m + 1 layer, the following equation holds: As per Equation (12), the following equation is obtained: Thus, the temperature of the contact surface is where ψ m = ζ m−1 γ m 1+ζ m−1 , and γ m = k m ·(e β m+1 ·h m+1 −1) k m+1 ·β m+1 ·h m . The m + 1 to m + j layers feature a Tcc-type temperature field, and, on the contact surface of two adjacent layers, the following equation holds: As per Equation (15), the following equation is obtained: Thus, the temperature of the contact surface can be obtained as follows: (3) The temperature of each stratum in the Tc-type temperature field in the lower n-m-j layer is expressed as follows: The geothermal gradient is As the temperature T m+j of the interface between the m + j layer and the m + j + 1 layer is the temperature of the interface between the Tc-type and Tcc-type temperature fields, the temperature for the two layers should be calculated separately. On the contact surface between the m + j layer and the m + j + 1 layer, there is As per Equation (20), the following equation is obtained: The temperature of the contact surface can be obtained as follows: . Similarly, where . The calculation specified in Equations (1)-(3) is a "chasing" process using the upper surface temperature T 0 already given in the model. On that basis, all other unknown variables can be obtained using the chasing method based on the known lower surface temperature T n . The relationship among the temperatures T i , T 0 , and T n of each layer can be described as follows: where, T i = T 0 , T n , when i = 0, n. Upon substituting Equation (24) into the temperature distribution equations in the chasing process specified in Equations (1)-(3), the equation of the temperature distribution in the stratified rock mass can be obtained.
This can be applied to the following formula: where h i is the thickness of layer i (m), k i is the thermal conductivity of the stratum in layer i (W·(m· • C) −1 ), T 0 is the known temperature at the upper boundary of the first layer ( • C), T n is the known temperature at the lower boundary of layer n ( • C), (n i is the porosity of each layer in the convection zone), v zi is the volume velocity of groundwater in the i layer along the z axis (m 3 /s), and v z is positive when the water moves downward, while it is negative if the water moves upward.

Effect of Seepage Velocity on Geothermal Field
In the present work, it is assumed that a geological formation consists of three layers: the upper thermal conduction layer comprising three rock strata, the middle thermal conduction-convection layer comprising three rock strata, and the lower thermal conduction layer also comprising three rock strata. The upper surface temperature of the model (T 0 ) is 35 • C, and the lower surface temperature (T 9 ) is 43.1 • C. The model upper surface elevation (H 0 ) is −500 m, and the lower surface elevation (H 9 ) is −770 m. The specific heat capacity (C w ) and density (ρ) of water are 4200 J/(kg· • C) and 1000 kg/m 3 , respectively. It is assumed that the thermal conductivity of the surrounding rock remains constant, and the rock porosity (n) in the thermal convection area is 0.15. The groundwater seepage velocity (v) was set to ±[0.5, 1, 1.5, 2] × 10 −8 m 3 /s. The relationship between the vertical temperature and the buried depth z in the model was analyzed. Specific parameters of the model are shown in Table 1. According to the above established model, the software Matlab (MathWorks Inc., Natick, MA, USA) was used to program Equation (25) to obtain the temperature-buried depth (z) curve along the vertical direction under different seepage velocities, as shown in Figure 2.   (25) to obtain the temperatureburied depth (z) curve along the vertical direction under different seepage velocities, as shown in Figure 2. As Figure 2 shows, the geothermal curve presents a linear distribution in the Tc-type region. Under the action of groundwater seepage, the curve in the Tcc-type region is convex when the seepage is upward, and it is concave when the seepage is downward. In general, the ground temperature in the model is higher when the groundwater percolates upward than when it percolates downward. In addition, a lower seepage velocity leads to a smaller curvature of the curve in the Tcc-type region. When the seepage velocity is close to 0, the curve in the Tcc-type region presents a linear distribution, which is similar to that in the Tc-type region. This indicates that the presence and migration of groundwater are among the major factors that affect the distribution of the temperature field.
According to the calculation equations for the geothermal gradient, the variation pattern of the geothermal gradient with the buried depth z under different seepage velocities can be obtained. Curves in Figure 3 show the specific variation patterns. As Figure 2 shows, the geothermal curve presents a linear distribution in the Tc-type region. Under the action of groundwater seepage, the curve in the Tcc-type region is convex when the seepage is upward, and it is concave when the seepage is downward. In general, the ground temperature in the model is higher when the groundwater percolates upward than when it percolates downward. In addition, a lower seepage velocity leads to a smaller curvature of the curve in the Tcc-type region. When the seepage velocity is close to 0, the curve in the Tcc-type region presents a linear distribution, which is similar to that in the Tc-type region. This indicates that the presence and migration of groundwater are among the major factors that affect the distribution of the temperature field.
According to the calculation equations for the geothermal gradient, the variation pattern of the geothermal gradient with the buried depth z under different seepage velocities can be obtained. Curves in Figure 3 show the specific variation patterns.  As Figure 3 shows, when there is groundwater seepage in the surrounding rock, the geothermal gradient changes accordingly. When the seepage is upward, the gradient increases obviously, and a greater seepage velocity leads to a larger increase in the gradient. On the contrary, the gradient decreases when the groundwater percolates downward, and a greater seepage velocity leads to a more drastic decrease in the gradient. When the seepage velocity is 0, the curve is a straight line, and, in this case, the Tcc-type temperature field is the only temperature field in the surrounding rock. This indicates that the formation of the Tcc-type temperature field entails not only the presence but also migration of groundwater.

Effect of Thermal Conductivity of Rock on Temperature Field
The rock thermal conductivity, a thermophysical parameter of rocks, mainly reflects the heat transfer characteristics of rocks and is also a necessary parameter for the calculation of the terrestrial heat flow value (Wu et al., 2019). In order to facilitate the analysis of the influence of rock thermal conductivity on the geothermal gradient distribution, the thermal conductivity was set to increase by values of equal differences. The thermal conductivity of each layer in the rock mass along the vertical direction is shown in Table 2. As Table 2 shows, the rock thermal conductivity of the five schemes was set to decrease or increase by an equal difference of 0, ±0.2, or ±0.4 along the vertical direction, and the mean conductivity of all the schemes remained at 2.5 k/W·(m·°C) −1 ; the seepage velocity was set to v = −1 × 10 −8 m 3 /s. The change pattern of the ground temperature with the buried depth z under the five schemes of rock thermal conductivity was obtained, as shown by the curves in Figure 4. As Figure 3 shows, when there is groundwater seepage in the surrounding rock, the geothermal gradient changes accordingly. When the seepage is upward, the gradient increases obviously, and a greater seepage velocity leads to a larger increase in the gradient. On the contrary, the gradient decreases when the groundwater percolates downward, and a greater seepage velocity leads to a more drastic decrease in the gradient. When the seepage velocity is 0, the curve is a straight line, and, in this case, the Tcc-type temperature field is the only temperature field in the surrounding rock. This indicates that the formation of the Tcc-type temperature field entails not only the presence but also migration of groundwater.

Effect of Thermal Conductivity of Rock on Temperature Field
The rock thermal conductivity, a thermophysical parameter of rocks, mainly reflects the heat transfer characteristics of rocks and is also a necessary parameter for the calculation of the terrestrial heat flow value (Wu et al., 2019). In order to facilitate the analysis of the influence of rock thermal conductivity on the geothermal gradient distribution, the thermal conductivity was set to increase by values of equal differences. The thermal conductivity of each layer in the rock mass along the vertical direction is shown in Table 2. As Table 2 shows, the rock thermal conductivity of the five schemes was set to decrease or increase by an equal difference of 0, ±0.2, or ±0.4 along the vertical direction, and the mean conductivity of all the schemes remained at 2.5 k/W·(m· • C) −1 ; the seepage velocity was set to v = −1 × 10 −8 m 3 /s. The change pattern of the ground temperature with the buried depth z under the five schemes of rock thermal conductivity was obtained, as shown by the curves in Figure 4.  As shown in Table 2, the low-thermal-conductivity area in Scheme 1 was located at the lower part of the model, while that in Scheme 5 was located at the upper part of the model. As shown in Figure 4, in Scheme 1 and Scheme 5, the high-geothermal-gradient area corresponded to the low-thermal-conductivity area. This is because the rock mass with low thermal conductivity blocked the transfer of heat that moved up from deep, and the temperature at the upper and lower boundaries of the rock mass varied greatly, resulting in a large geothermal gradient inside the thermal conductivity formation. In addition, comparison of the thermal conductivity variances of all five schemes revealed that a greater variance of thermal conductivity, i.e., a greater difference in the thermal conductivity of the rock mass in the vertical direction, resulted in a greater curvature of the temperature-z curve.

Case Study
Huainan Coalfield is located in the south of the coal accumulation area in north China, with a length of about 100 km, a width of 20-30 km, and an area of 2500 km 2 , as shown in Figure 5. Strata in Huainan coalfield have typical sedimentary characteristics of a north China platform type. Strata mainly developed from a Lower Proterozoic to Quaternary overlay, but the strata of upper Ordovician O1 and Upper and Middle Triassic T1~2 to Middle Jurassic J2 are missing. The geothermal formation of Huainan Coalfield is closely related to the lithological changes [19][20]. This coalfield, with a total thermal energy reserve of 2.32 × 10 16 kJ, features a wide distribution of thermal energy reserves, even thickness, and uniform lithology. Furthermore, with rich geothermal energy, the coalfield is categorized as a medium-low temperature (II) class and II-3 type stratified geothermal reserve according to the national standard of China (GB/T 1161-2010:3-45) [21]. As shown in Table 2, the low-thermal-conductivity area in Scheme 1 was located at the lower part of the model, while that in Scheme 5 was located at the upper part of the model. As shown in Figure 4, in Scheme 1 and Scheme 5, the high-geothermal-gradient area corresponded to the low-thermal-conductivity area. This is because the rock mass with low thermal conductivity blocked the transfer of heat that moved up from deep, and the temperature at the upper and lower boundaries of the rock mass varied greatly, resulting in a large geothermal gradient inside the thermal conductivity formation. In addition, comparison of the thermal conductivity variances of all five schemes revealed that a greater variance of thermal conductivity, i.e., a greater difference in the thermal conductivity of the rock mass in the vertical direction, resulted in a greater curvature of the temperature-z curve.

Case Study
Huainan Coalfield is located in the south of the coal accumulation area in north China, with a length of about 100 km, a width of 20-30 km, and an area of 2500 km 2 , as shown in Figure 5. Strata in Huainan coalfield have typical sedimentary characteristics of a north China platform type. Strata mainly developed from a Lower Proterozoic to Quaternary overlay, but the strata of upper Ordovician O 1 and Upper and Middle Triassic T 1~2 to Middle Jurassic J 2 are missing. The geothermal formation of Huainan Coalfield is closely related to the lithological changes [19,20]. This coalfield, with a total thermal energy reserve of 2.32 × 10 16 kJ, features a wide distribution of thermal energy reserves, even thickness, and uniform lithology. Furthermore, with rich geothermal energy, the coalfield is categorized as a medium-low temperature (II) class and II-3 type stratified geothermal reserve according to the national standard of China (GB/T 1161-2010:3-45) [21].
Guqiao Mine and Dingji Mine ( Figure 5) are two mines located in the Panxie mining area of Huainan Coalfield. Many temperature-measuring boreholes have been deployed at each stage of exploration in these two mines, among which the holes XLZM1 and 16-11 are located to the west of Guqiao Mine and to the north of Dingji Mine, respectively. According to exploration data of the geothermal resources, there are abundant geothermal water resources in the Huainan Coalfield below the buried depth of −600 m. The data measured within the buried depth z intervals of (−800 m, −1500 m) and (−600 m, −1020 m) of the two boreholes were taken as the research data. The XLZM1 hole penetrates the Permian, Carboniferous, Ordovician, and Cambrian strata, among which the main aquifers are Cambrian, Ordovician, and Carboniferous. The 16-11 hole penetrates Permian and Carboniferous strata, the latter of which is the main aquifer. Table 3 shows the specifics of these strata. Guqiao Mine and Dingji Mine ( Figure 5) are two mines located in the Panxie mining area of Huainan Coalfield. Many temperature-measuring boreholes have been deployed at each stage of exploration in these two mines, among which the holes XLZM1 and 16-11 are located to the west of Guqiao Mine and to the north of Dingji Mine, respectively. According to exploration data of the geothermal resources, there are abundant geothermal water resources in the Huainan Coalfield below the buried depth of −600 m. The data measured within the buried depth z intervals of (−800 m, −1500 m) and (−600 m, −1020 m) of the two boreholes were taken as the research data. The XLZM1 hole penetrates the Permian, Carboniferous, Ordovician, and Cambrian strata, among which the main aquifers are Cambrian, Ordovician, and Carboniferous. The 16-11 hole penetrates Permian and Carboniferous strata, the latter of which is the main aquifer. Table 3 shows the specifics of these strata. The parameters of elevation and temperature of the upper and lower boundaries of the two boreholes, porosity, fluid density, and specific heat capacity of the aquifer are shown in Table 4.  The parameters of elevation and temperature of the upper and lower boundaries of the two boreholes, porosity, fluid density, and specific heat capacity of the aquifer are shown in Table 4. Table 4. Other parameters of the strata and the aquifer at the two temperature-measuring holes. All the known parameters presented in Tables 3 and 4 were put into the model for calculation, and the calculation results were compared with the measured results. Given the lithology histogram and measured geo-temperature curve of the formation at these two boreholes, it was confirmed that there are both upper Tc-type and lower Tcc-type temperature fields in the study interval.

Boreholes
By adjusting the elevation of the boundary between the two temperature fields and the velocity and direction of groundwater seepage in the Tcc-type region, we fitted the analytical solution curve and the measured geothermal curve. The elevation value of the boundary between groundwater seepage velocity and temperature field of each aquifer is shown in Table 5. The corresponding fitting of the analytical solution of geothermal curve with the measured value is shown in Figure 6.
All the known parameters presented in Tables 3 and 4 were put into the model for calculation, and the calculation results were compared with the measured results. Given the lithology histogram and measured geo-temperature curve of the formation at these two boreholes, it was confirmed that there are both upper Tc-type and lower Tcc-type temperature fields in the study interval.
By adjusting the elevation of the boundary between the two temperature fields and the velocity and direction of groundwater seepage in the Tcc-type region, we fitted the analytical solution curve and the measured geothermal curve. The elevation value of the boundary between groundwater seepage velocity and temperature field of each aquifer is shown in Table 5. The corresponding fitting of the analytical solution of geothermal curve with the measured value is shown in Figure 6.   As Figure 6 shows, Tc-type and Tcc-type regions are both present in the rock strata at the two boreholes XLZM1 and 16-11. The interfaces of the two regions are shown in the figure. Above the interface, the geothermal temperature curves are inclined straight lines,   Figure 6 shows, Tc-type and Tcc-type regions are both present in the rock strata at the two boreholes XLZM1 and 16-11. The interfaces of the two regions are shown in the figure. Above the interface, the geothermal temperature curves are inclined straight lines, and the temperature field type is the Tc-type; below the interface, the geothermal temperature curves are convex, the temperature field type is the Tcc-type, and vertical geothermal water occurs in the area. The fitting result revealed that the difference between the measured value and the calculated value was less than 0.2 • C, indicating high goodness of fit and that the value of groundwater seepage velocity was set reasonably. In summary, geothermal water exists in the strata of both boreholes, with buried depths ranging from −1020 m to −1500 m and from −890 m to −1020 m, respectively. Geothermal water migrates vertically upward. The vertical seepage velocity of groundwater in the Carboniferous, Ordovician, and Cambrian aquifers at Guqiao Mine was 0.7 × 10 −8 m 3 /s, 1.35 × 10 −8 m 3 /s, and 0.60 × 10 −8 m 3 /s, respectively, and that at the Carboniferous aquifer of Dingji Mine was 0.75 × 10 −8 m 3 /s.
Regarding the geothermal gradient, the geothermal gradient obtained by modeling and the measured gradient were compared [22,23]. As Figure 7 shows, affected by vertical groundwater seepage, the geothermal gradient curves of the two boreholes first increased and then decreased with the increase in buried depth z, indicating the presence of vertical and upward migration of geothermal water in both mines. The variation trend of the geothermal gradient obtained by modeling was nearly the same as that of the measured gradient, and the maximum error was only 0.03 • C/100 m, indicating the high accuracy of the proposed model in reverse-calculating the seepage velocity of groundwater and identifying the occurrence of geothermal water.

Conclusions
On the basis of the one-dimensional heat conduction-convection equation, a thermal effect model of vertical groundwater migration in the stratified rock mass was established, and the equations for temperature distribution in stratified rock mass were deduced. The proposed model was employed to analyze the impact of thermal conductivity of the surrounding rocks and the groundwater seepage velocity on the temperature distribution in geo-formations. Furthermore, the model was applied for identification of the parameters

Conclusions
On the basis of the one-dimensional heat conduction-convection equation, a thermal effect model of vertical groundwater migration in the stratified rock mass was established, and the equations for temperature distribution in stratified rock mass were deduced. The proposed model was employed to analyze the impact of thermal conductivity of the surrounding rocks and the groundwater seepage velocity on the temperature distribution in geo-formations. Furthermore, the model was applied for identification of the parameters of the temperature curve at two temperature-measuring boreholes in coal mines. The following conclusions were drawn: (1) The calculation equations for ground temperature distribution in layered strata were derived, which could be employed to calculate the variation of ground temperature caused by vertical seepage of geothermal water and the thermal conductivity of surrounding rocks at any location under the actual geological formations. The parameters required in the calculation equations are easily accessible, and the calculation results are highly accurate as the equations fully consider the factors in actual scenarios. (2) In the present work, the geothermal heat transfer patterns in the Tc-type and Tcc-type strata in stratified rock masses were summarized. In the Tc-type strata, the ground temperature was mainly affected by the thermal conductivity of the rock; higher thermal conductivity of the rock led to smaller variations of the ground temperature, while lower conductivity led to larger variations of the ground temperature, and the temperature curve resembled a broken line. In the Tcc-type-type strata, the ground temperature was mainly subject to the vertical seepage of geothermal water; the ground temperature curve was convex when the seepage was upward, but concave when the seepage was downward, and the curvature of the curve grew as the seepage velocity increased. In addition, different from the homogeneous model in which the ground temperature curves were straight lines, the layered model fully showed the geothermal variations caused by the difference of thermal conductivity of the rock mass. (3) Application of the model to a study case showed that the difference between ground temperature obtained by modeling and the measured temperature was less than 0.2 • C, the maximum geothermal gradient error was only 0.03 • C/100 m, and the curve fitting degree was high, indicating that the geothermal transfer model of stratified rock masses proposed in this study can appropriately identify the Tc-type and Tcc-type regions in the temperature curve in actual stratified formations. Moreover, the groundwater seepage velocity in the Tcc-type region can be obtained by reverse calculation. The research result is envisioned to provide a geological basis for further development of geothermal water resources.  Informed Consent Statement: Not applicable.

Data Availability Statement:
The data used in this paper can be accessed by contacting the corresponding author directly.