Abstract

Deep-circling thermal water upwelling and trickling to high geo-temperature roadway obviously alter the microclimate in mines, which brings difficulty to the prediction of airflow temperature and humidity. This is the basis of air-conditioning cooling load calculation. The heat and mass transfer between trickling water and airflow is rather complicated. Moreover, humid air exhibits the accumulation effect of heat and humidity in the long-distance flow process. In this paper, an apparatus was designed and developed to explore the influence of thermal water trickling on the airflow thermal parameters of a section of roadway (1L–39L, in which 1L–9L is the trickling section). The results show the following (1) With the rise of trickling water temperature, the total enthalpy difference of dry air in the roadway increases within a small range and that of humid air goes up nonlinearly. Besides, the increase of trickling water flow rate has an insignificant effect on the sensible heat of the airflow, while it plays a notable role in increasing the latent heat of the airflow. (2) High trickling water temperature results in a higher growth rate of humidity ratio at 19L than those at 29L and 39L in the early stage of thermal water trickling. Meanwhile, sensible heat exchange, which becomes strong after thermal water trickles for over 30 min, complicates the enthalpy difference variation rates of wet air at the three measuring points. (3) The three measuring points in the 19L–39L section all display a process of enthalpy growth with time. In the case of point 39L, the enthalpy difference of humid air surges sharply when the trickling water temperature is or the flow rate is 200 ml/min. The research results boast some reference value for thermal water management and microclimate change forecasting after the airflow passes through a trickling roadway.

1. Introduction

Mine thermal hazards are one of the geological disasters on a global scale [14]. During deep mining, roadways are often excavated in surrounding rock with high geo-temperature field, which will inevitably bring about thermal hazards because heat dissipation of the rock is one of the main heat sources of roadways [5, 6]. Mine microclimate parameters, including temperature, humidity, and enthalpy, are the essential indicators for evaluating the severity of thermal hazards. There are 47 coal mines whose depths exceed 1,000 m in China [7, 8], such as Jiahe Coal Mine (1,200 m), Suncun Coal Mine (1,500 m), and Tangkou Coal Mine (1,025 m). Airflow temperatures of these mines are mostly higher than 30°C (Table 1) so that the operating environment for workers is rather harsh. Slight thermal hazards can usually be solved by means of thermal insulation and ventilation adjustment, while serious thermal hazards have to be dealt with by the mine air-conditioning system that cools the airflow directly [9]. The required cooling amount and the cooling capacity of conditioners are calculated based on the actual measurement and dynamic prediction of airflow temperature, which is a complex system engineering. Thermal water is common in some mines. On the one hand, it is a feasible resource and a potential approach to the utilization of geothermal energy under specific hydrogeothermal conditions. On the other hand, it can cause remarkable underground microclimate changes [10]. Thus, it needs to be managed rationally [11, 12]. As shown in Table 1, the thermal water (62°C) in Pingdingshan No. 8 Coal Mine notably heats the airflow and raises the temperature in working face to 35°C. In addition, in Hujiawan Coal Mine, when the airflow passes through the trickling section, 82% of its heat increment is attributed to the thermal water (40°C). Thus, the research on the heat and humidity transfer between airflow and high geo-temperature surrounding rock with thermal water trickling is of great significance for the prediction of mine microclimate and the design of air-conditioning.

As the basis of airflow temperature prediction, heat and mass transfer between airflow and surrounding rock is a complex dynamic process. А. Н. Щербань defined this process as unstable heat transfer and used the dimensionless physical parameter to represent the process, but he did not give the theoretical solution of [13]. Numerous scholars have solved the coefficient of unstable heat transfer by means of the analytical solution method (Cen [14], Zhou [15]), the Laplace transform method (Sun [16], Yakovenko [17]), the separated variable method (Wang [18]), the numerical method (Qin [19], Roy [20]), and so forth. Nevertheless, these solution methods are of complicated processes because of the nonlinearity of thermal differential equations, and the results are lack of practicability. Wang et al. [21] obtained an empirical formula for by performing regression analysis on measured data. Cheng et al. [22] proposed a calculation model for prediction based on an improved momentum BP algorithm. In order to investigate the effect of roadway wall moisture on surrounding rock heat dissipation, Chow and Chung [23] adopted an experimental method to analyze the process of water evaporation from the roadway wall into the air, and they focused on the evaporation rate and conversion temperature. Hou and Shen [24] put forward the calculation formula for heat transfer between a partial wet roadway wall and the airflow. Liu et al. [25] and Gao et al. [26] calibrated by introducing the coefficient “roadway surface humidity.” It can be seen that these studies mostly adopted theoretical analysis or numerical calculation methods, yet they failed to analyze the variation law of in depth. Consequently, the actual effect of changes in boundary conditions on heat and mass transfer is still difficult to be fully grasped. Thermal water trickling forms a liquid film in the roadway, and direct contact between air and water will cause heat and mass transfer. Scholars have conducted extensive research on spray room [27], cooling tower [28], steam cooler [29], and so forth in light of the film theory, but their research results are not quite applicable to roadways with trickling thermal water because airflow that has passed through the trickling section exhibits the accumulation effect of heat and humidity in the process of long-distance flow.

Despite the abundant progress made in the research on heat and moisture transfer between surrounding rock and airflow in dry or wet roadways, there are a few related studies on roadways with trickling thermal water. To study roadways with trickling thermal water, theoretical and numerical calculations can hardly be carried out because of the complex boundary conditions and the coefficient of heat and mass transfer that are difficult to determine. In contrast, similarity simulation experiments are more advantageous in the study on hydrodynamics and heat transfer. Wang et al. [30], Zhang et al. [31], Zhu et al. [32], and Wang [33] have designed experimental apparatuses for heat transfer between surrounding rock and airflow, but these apparatuses are just applicable to dry roadways. On the basis of previous studies, an improvement was made to the similar simulation apparatus independently developed by our team, so that the apparatus can be used for exploring the microclimate oscillation with trickling water. The improvement is of great significance for the prediction of underground climate conditions, the adjustment of the ventilation network, and the design of the air-conditioning system. Ultimately, it serves the treatment of mine thermal hazards.

In this paper, we focused on investigating the effects of thermal water trickling on thermal environment change in the high geo-temperature roadway. Firstly, the geological genesis of thermal hazards in Sanhejian Coal Mine was analyzed based on geological survey and field measurement. Then, with Sanhejian Coal Mine taken as the research background, an apparatus was designed and improved to explore the effect of thermal water trickling on a section of roadway. Finally, the temperatures, relative humidity (RH), and enthalpy differences of the measuring points were presented and analyzed.

2. Geological Genesis of Thermal Hazards: A Case Study

2.1. Geological Structure and Strata of Sanhejian Coal Mine

Sanhejian Coal Mine is located in the southwest part of the secondary folds of Tengyu Anticline in the middle of Fengpei Coalfield, as shown in Figure 1. The regional faults developed in Fengpei Coalfield include the nearly EW-trending Yutai Fault, the NS-trending Jiaxiang Fault, and Sunshidian Fault (shattered zone width 38.5–51.53 m). These faults, combined with other subordinated faults, cause the stepped fault blocks [34]. In the south boundary of Sanhejian Coal Mine, the east section is F1 Fault (shattered zone width 31–32.29 m) and the west section is Liuwuzhuang Fault. The northeast boundary is the Sunshidian Fault, and the northwest is a monoclinic structure tilting to the deep part. A secondary fault F2 exists inside the minefield.

The stratigraphic conditions in Sanhejian Coal Mine are given in Table 2. The main mining seams of Sanhejian Coal Mine are 7# and 9# coal seams of Permian Period Shanxi Formation and 17# and 21# coal seams of Carboniferous Period Taiyuan Formation. The main aquifer is the limestone of the Ordovician Period Majiagou Formation.

2.2. Geological Genesis of Thermal Hazards in Sanhejian Coal Mine

Sanhejian Coal Mine is one of the mines suffering from thermal hazards in China [11]. In March 2018, the roadway airflow parameters in the ventilation route of A-I working face were measured (Figure 2). The temperature of the airflow increased from to (Figure 3(a)). The RH at each measuring point was maintained above 85%, and the RH in A-I working face approximated saturation. Thermal hazards in Sanhejian Coal Mine are mainly caused by high geo-temperature (HGT) field and thermal water upwelling (TWU).

The HGT field is mainly attributed to the regional geothermal background and the changes in heat transfer conditions in the shallow crust. (1) The average terrestrial heat flow in Sanhejian Coal Mine is (Table 2) which is higher than the global average . (2) Affected by F2 Fault and Longgu Anticline near A-I working face, the crystalline basement and coal basement of the area have been uplifted, resulting in the heat flow accumulation at the top of the anticline or at the uplifted part. Therefore, a large severe-anomaly zone of geo-temperature has been formed in this region (Zone II in Figure 1(b) [11]). (3) The Quaternary, Cretaceous-Jurassic, and Carboniferous-Permian strata jointly form a good high thermal resistance cover (thickness 600–1200 m) with poor aquosity and thermal conductivity (Table 2).

The additional heat source caused by the TWU along with the HGT field makes the mine climate more complex. It was a process of warming, humidifying, and enthalpy increasing along the ventilation route of A-I working face. As shown in Figures 2 and 3, the problem of thermal water trickling or seeping existed in the ventilation route. Consequently, the relative humidity and the humidity ratio (HR) are of large values (e.g., measuring point 6: 87% RH, 0.015 kg/kg HR); the ratio of water vapor enthalpy to humid air enthalpy reached 65%–71% (Figure 3(b)); the latent heat exchange of airflow accounted for a large part of the heat source.

The problem of TWU also exists in other areas of Sanhejian Coal Mine. For example, a water inrush with a temperature over and water pressure of 7.6 MPa occurred in A-II working face in Figure 1 [36]. Where does thermal water in the roadway come from? Sunjiadian Fault is a large regional fault, but it has a small gap in Sanhejian Mine Field () and is cut through by neotectonic faults (Figure 4). As a result, the Ordovician limestone aquifer is still capable of conducting water laterally. There is an Ordovician suboutcrop area in the axis of Tengyu Anticline. The water head difference between the water supply end and the water trickling end is as high as 980 m, which provides dynamic for the radial flow and vertical upwelling of the Ordovician water. The neotectonic faults cut through the Ordovician strata and coal strata. Meanwhile, fracture propagation is promoted by mining activities. Then, the deep-circle water gets heated by the rock, upwells along the fault, and finally trickles to the roadway.

3. Methodology

As described in Section 2, the HGT field and the TWU have a significant impact on the mine climate. Airflow absorbs the heat released from surrounding rock along the roadway. However, the temperature and humidity of airflow will alter notably if thermal water trickling occurs at multiple positions in the route, which brings trouble to climate forecasting and ventilation management. In Section 3, an apparatus for heat and mass transfer in the modelling roadway was designed and developed, and the precast fracture representing the fracture zone connected with the roadway is supplied with water whose flow rate and temperature are controllable (Figure 4).

3.1. Experimental Apparatus

The experimental apparatus consists of a modelling HGT roadway, a pump, a water flow thermostat, a Surrounding Rock Heating System (SRHS), a Constant Temperature and Humidity Ventilator (CTHV), and a data acquisition system (Figures 4 and 5). The geometric similarity constant of the device , the time similarity constant , and the wind speed similarity constant . Detailed design of the experimental apparatus can be found in the literature [31]. The modelling HGT roadway was composed of six cylinders (1 m/cylinder) that were poured (or filled) with materials similar to the surrounding rock. Temperature sensors were set at 1R, 2R, 3R, and 4R (R = 0.1 m) along the radial direction of the cylinder, and temperature and humidity sensors were set at 1L, 9L, 19L, 29L, and 39L (L = 0.1 m) of the airflow axial centres. The pump and water flow thermostat provided the required water flow temperature and rate. The SRHS ensured the constant temperature of the virgin rock. The CTHV was responsible for keeping the inlet airflow steady.

The 1L–9L section in Figures 5 and 6 was the trickling water section. The casing (at 3R in the cylinder) and the precast fracture surface were preembedded in the surrounding rock (Figure 5). The designed casing was a coaxial double-layered steel tube whose internal casing was slotted upward and external casing was slotted downward. Having overcome the surface tension, the surfactant-containing liquid flowed from the internal casing to the external one and finally entered the roadway along the fracture surface to form a water film.

The fracture surface (width 2 mm, length 800 mm) was equivalent to a 4 cm fracture zone of the mine, and the total length of the roadway and the length of the trickling section were 80 m and 16 m, respectively.

3.2. Experimental Contents and Procedures
3.2.1. Experimental Contents

The virgin rock temperature and inlet airflow parameters in each group were kept unchanged. Since the rock temperature in Sanhejian Coal Mine at the depth of is , in this experiment, the virgin rock temperature was set as . In addition, the airflow speed was ; the temperature of the inlet airflow was ; and the humidity of the inlet airflow was .

In Sanhejian Coal Mine, since the flow rate of trickling water within the range of 16 m shall not exceed , the flow rate of fracture water shall not exceed 250 ml/min in this experiment according to the similarity ratio. The highest temperatures of trickling water are 60– in coal mines (e.g., , Pingdingshan No. 8 Mine) and 70– in metal mines (e.g., , Jiudian Gold Mine). Based on this fact, five experimental schemes were designed, as listed in Table 3.

3.2.2. Experimental Procedures

The tests were conducted in accordance with the following steps: (1) to heat surrounding rock to by the SRHS; (2) to regulate parameters of inlet airflow by the CTHV (, , ) and maintain ventilation until the surrounding rock-airflow system reach heat balance—airflow parameters at this moment are disclosed in Figure 6—(3) to provide the required trickling conditions by the pump and to keep the ventilation and trickling state for 180 min; and (4) to record the temperature and RH of airflow in the 1L–39L section. By this above procedure, tests were conducted under varying trickling conditions in Schemes 1–5. In the tests, the effect of thermal water on the surrounding rock temperature was ignored.

3.3. Calculation of Enthalpy

During 180 min water flow supplying to the 1L–39L section, heat and mass transfer promoted by the enthalpy difference occurs between airflow and trickling thermal water [37]. Humid air absorbs/releases heat from/to the surrounding medium when the enthalpy difference is positive/negative. The enthalpy of humid air, expressed based on 1 kg of dry air, is the sum of the enthalpies of 1 kg of dry air and d kg of water vapor. The enthalpy of humid air iswhere is the enthalpy of 1 kg of dry air, ; is the humidity ratio in air, ; is the latent heat of vaporization of water vapor, ; 1.85 is the average specific heat at a constant pressure of water vapor at normal temperature, ; and is dry bulb temperature, .

The humidity ratio in air iswhere is the relative humidity of air, ; is the atmospheric pressure, ; and is the saturated vapor pressure of air at the temperature , .

4. Results and Discussion

4.1. Temperature and Relative Humidity Variations

The airflow temperature rarely varies when the trickling water is of low temperature and low flow rate, while it rises remarkably when the trickling water is of high temperature and high flow rate. Under the five trickling conditions, the temperature at 1L along the axial direction of the cylinder decreases by about . When the temperature is lower () or the flow rate is lower (), the temperatures at 9L, 19L, 29L, and 39L rise slightly within the range of . In contrast, the maximum temperature rise can be (Figure 7(a) under the trickling condition of 100–80 (100–80 refers to , hereinafter described in the same form).

As shown in Figure 7(b), after thermal water supplying, the relative humidity in the 9L–39L section is all notably higher than the initial value (Figure 6(b)), and the value at 19L is the greatest, followed by those at 29L and 39L in turn all the time. Before thermal water supplying, the relative humidity in the 1L–39L section declines and is lower than 20% of the initial value. This is because the humidity ratios in this section are basically equal (0.003–0.004 kg/kg), but the saturated steam pressure along the axial direction rises. After water flow supplying, the relative humidity in the 9L–39L section grows (the growth rates gradually decrease). Among the five measuring points, the relative humidity of 19L is the largest, reaching up to 51.7% under the condition of 200–60 (Figure 7(b)). The reason for this distribution characteristic is that when the airflow passes through the 1L–9L water trickling section, the evaporation of water film results in the maximum value of the water vapor pressure at 19L. However, the relative humidity keeps falling in the 9L–39L section for two reasons: first, the water vapor pressures at 29L and 39L are lower than those at 19L. Second, the saturated steam pressure is higher than that at 19L.

4.2. Enthalpy Difference Variation
4.2.1. Total Enthalpy Difference Variation of Trickling Roadway

The total enthalpy difference of airflow, which represents the variation of airflow energy, is the impetus of heat and mass transfer. It can be calculated through the following equation:where is the enthalpy difference, ; is the enthalpy difference at 39L during water flow supplying; is the enthalpy value at 39L at the end of water flow supplying; is the enthalpy value at 39L at the beginning of water flow supplying; is the enthalpy difference at 1L during water flow supplying; is the enthalpy value at 1L at the end of water flow supplying; and is the enthalpy value at 1L the beginning of water flow supplying.

The curves of total enthalpy difference variations under the five trickling conditions are presented in Figure 8. The rise of trickling water temperature causes an increase in both the sensible heat and the latent heat of the airflow (Figure 8(a)). At the same trickling water flow rate (), the total enthalpy difference of dry air increases within a small range, while that of humid air goes upward nonlinearly with the rise of trickling water temperature. The increase of trickling water flow rate has an insignificant effect on the sensible heat of the airflow, while it plays a notable role in increasing the latent heat of the airflow (Figure 8(b)). At the same temperature (), the enthalpy difference of dry air fluctuates, while that of humid air goes upward linearly with the rise of trickling water flow rate.

4.2.2. Enthalpy Difference Variation Rate at Each Point under the Condition of High Trickling Temperature and High Flow Rate

Figures 9(a) and 9(c) display the curves of enthalpy difference variations of humid air and dry air at different times with position under the condition of 100–80 (here, the enthalpy difference refers to the difference between the initial enthalpy value and the final one at different times at the same position). Figure 9(b) shows the variations of humidity ratios at various measuring points. It can be concluded that the humidity ratio at 19L changes more and more slowly with time, and the increases of humidity ratios at 29L and 39L experience a hysteretic process. At the end of water flow supplying, the humidity ratios at the three points become basically the same. The hysteretic process causes the higher relative humidity growth rate at 19L than those at 29L and 39L in the first 30 min (Figure 7(b)). Figure 9(d) discloses the change curve of dry bulb temperature during water flow supplying at 19L, 29L, and 39L. It can be seen that the temperature difference between 39L and 19L gradually increases.

High trickling water temperature (e.g., ) significantly affects the enthalpy difference variation rates of wet air at 19L, 29L, and 39L with time. The variation rates at the three points differ remarkably before the surrounding rock-airflow system reaches heat balance (Figure 9(a)). The enthalpy difference of wet air comes from both sensible heat and latent heat. The change in sensible heat, positively correlated with the temperature difference, is represented by the enthalpy difference of dry air. In contrast, the change in latent heat, positively correlated with the humidity ratio, is represented by the enthalpy difference of water vapor. (1) In the first 30 min of thermal water trickling, the enthalpy difference of humid air at 19L is the largest. In this stage, due to the high growth rate of humidity ratio at 19L, the enthalpy difference of water vapor surges rapidly. As shown in Figure 9(b), the humidity ratio at 19L is 0.003 kg/kg greater than that at 39L at 30 min. Nevertheless, the enthalpy differences of dry air at the three points differ slightly. (2) After thermal water trickles for 30–90 min, the enthalpy differences of humid air at 29L and 39L increase, and the differences at the three points become basically equal at 90 min. In this stage, 29L and 39L have already passed the hysteretic period. Their humidity ratios grow at higher rates than that at 19L, and their enthalpy differences of dry air surge. As shown in Figure 9(d), the temperature at 39L is higher than that at 19L at 90 min. (3) After thermal water trickles for over 90 min, the enthalpy difference of humid air at 39L becomes the largest. In this stage, the amounts of latent heat exchange at the three points all reach their high values. Meanwhile, their humidity ratios rise synchronously at 180 min. However, because of the continuously increasing amounts of sensible heat exchange, the enthalpy difference of humid air at 39L is greater than that at 19L after the end of water flow supplying.

Figure 10 presents curves of enthalpy difference variations of humid air and dry air at different times during water flow supplying under the condition of 200–60. High trickling water flow rate (e.g., ) will not remarkably alter the enthalpy difference variation rates of humid air at 19L, 29L, and 39L with time (Figure 10(a)). A high trickling water flow rate leads to small ranges of temperature rise at the three points, and the amounts of sensible heat at the three points do not differ much (Figure 10(b)). In other words, the enthalpy differences of humid air at 19L, 29L, and 39L are mainly influenced by latent heat. (1) In the first 30 min of thermal water trickling, the enthalpy difference of humid air at 19L is the largest. In this stage, due to the high growth rate of humidity ratio at 19L, the enthalpy difference of water vapor surges rapidly. A comparison between Figure 11 and Figure 9(b) reveals that the increases in the humidity ratio at 29L and 39L present no obvious hysteretic processes. (2) After thermal water trickles for over 30 min, the enthalpy difference of humid air at 19L remains the largest. In addition, the enthalpy difference growth rates of humid air at 19L, 29L, and 39L do not differ notably at a certain moment because the humidity ratios at the three points rise synchronously after 30 min (Figure 11).

Evidently, high trickling water temperature results in a higher growth rate of humidity ratio at 19L than those at 29L and 39L in the early stage of thermal water trickling. Meanwhile, sensible heat exchange complicates the enthalpy difference variation rates of humid air at the three points, which is clearly exhibited after thermal water trickles for over 30 min. These facts bring difficulty to the temperature and humidity prediction of the trickling roadway.

4.2.3. Enthalpy Difference Variations at 39L with Different Trickling Temperatures and Flow Rates

After the airflow passes through the trickling section, enthalpy difference variations at 39L with trickling water temperature and flow rate are shown in Figures 12 and 13, respectively. The enthalpy difference at 39L grows with time during water flow supplying. Besides, the enthalpy difference of humid air at 39L rises as well with the increases of trickling water temperature and flow rate.

When the trickling water flow rate is 100 ml/min, the dry air enthalpy difference at 39L is 0.4– under the conditions of and . When the trickling water temperature is , the enthalpy difference of dry air at 39L surges sharply to , and that of humid air at 39L jumps significantly as well (Figure 12(a)).

When the trickling water temperature is , the dry air enthalpy difference at 39L is under the conditions of and . When the trickling water flow rate is , the enthalpy difference of dry air at 39L increases by , and that of humid air at 39L soars (Figure 13(a)).

5. Conclusions

Thermal hazards in Sanhejian Coal Mine are mainly caused by high geo-temperature field and thermal water upwelling. Field measurement results suggest that thermal water trickling to roadway has a significant impact on the mine microclimate. In this paper, a similarity experiment was performed for investigating the influence of thermal water trickling on airflow parameters of a section of roadway. Several beneficial conclusions were drawn. The research results allow the prediction of temperature and RH to be deliberated:(1)The total enthalpy difference variation of the trickling roadway was analyzed. With the rise of trickling water temperature, the total enthalpy difference of dry air in the roadway increases within a small range, and that of humid air goes up nonlinearly. Besides, the rise of trickling water flow rate has an insignificant effect on the sensible heat of the airflow, while it plays a notable role in increasing the latent heat of the airflow.(2)The enthalpy difference variation rates at different measuring points were compared under the condition of high trickling water temperature and flow rate. High trickling water flow rate (e.g., ) will not cause notably varying enthalpy difference variation rates of humid air with time at 19L, 29L, and 39L, but the variation rates at the three points differ remarkably before the surrounding rock-airflow system reach heat balance under the condition of high trickling water temperature (e.g., ). Evidently, high trickling water temperature results in a higher growth rate of humidity ratio at 19L than those at 29L and 39L in the early stage of thermal water trickling. Meanwhile, sensible heat exchange, which becomes strong after thermal water trickles for over 30 min, complicates the enthalpy difference variation rates of humid air at the three points.(3)The curves of temperature and RH variations with time at 19L, 29L, and 39L under the five schemes were recorded. The three points all display a process of enthalpy growth with time. For example, for 39L, the enthalpy difference of humid air surges sharply when the trickling water temperature is or the flow rate is .

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 no conflicts of interest.

Acknowledgments

The authors are grateful to Dr. Zhaoyang Ma from Monash University for his support and contribution to language proofreading. This research was funded by National Natural Science Foundation of China, under Grant no. 51674242.