Abstract
Thermal insulation gunite (TIG) in roadways is an effective method for regional thermal hazard control in mines. The development of mine TIG materials is the foundation of thermal insulation technique. However, some conventional and advanced insulation materials are inapplicable to deep mines which are rather humid with high in situ stress and high geo-temperature. In this study, a kind of fly ash-inorganic mineral TIG material was developed and applied to the modelling of a high geo-temperature roadway. Moreover, the thermal insulation effect of the TIG layer was analyzed, and the temperature field characteristics of the TIG surrounding rock were discussed. Results reveal that (1) the TIG layer has a significant impact on the heat release of the wall and stability of the surrounding rock temperature field; (2) the initial temperature disturbance times, temperature disturbance ranges, and temperature drop rates differ with whether a TIG layer exists or not; (3) after the TIG roadway starts to be ventilated, the thermal flux densities tend to be consistent, which indicates the end of temperature disturbance; besides, the dimensionless temperature shares an exponential relation with the dimensionless radius; and (4) the characteristics of temperature drops vary with the radial positions of the surrounding rock. The research results provide a certain reference for thermal hazard control, temperature prediction, and ventilation network adjustment.
1. Introduction
Thermal hazards, one of the common problems in modern mines, are an inevitable geological disaster in the process of mine development to the deep [1–3]. The problem of thermal hazards appeared earlier in developed mining countries such as the United States [4], Australia [5], Poland [6], and South Africa [7]. To this day, there have been about 140 high geo-temperature mines in China [8, 9]. High geo-temperature is the main reason for the formation of thermal hazards which are inevitably brought about by deep mining [10–12]. By 2016, there have been 47 coal mines whose depths exceed 1,000 m in China [13, 14], such as the Suncun Coal Mine (1,501 m, geo-temperature 48°C), Sanhejian Coal Mine (1,010 m, geo-temperature 46.8°C), and Zhangji Coal Mine (1,260 m, geo-temperature 51.5°C) [8], and over 150 coal mines whose depths exceed 1,000 m abroad, such as the Pniowek Coal Mine (1,000 m, geo-temperature 45°C) in Poland [6]. In addition, some coal mines which are 600–1,000 m deep have high geo-temperatures as a result of local heat accumulation, such as the Pingdingshan No. 8 Coal Mine (660 m, geo-temperature 43°C) [1]. Domestic and foreign researches on cooling of high geo-temperature mines are mainly focused on mechanical cooling technique [7, 15], including central air conditioning cooling, ice cooling, air cooling [16], electric-thermal-glycol cooling, and deep well HEMS cooling [17]. However, the high cost, large cooling capacity losses, and large energy consumption of mechanical cooling bring considerable burden to enterprises. The thermal insulation and cooling of a roadway is to shoot the mine thermal insulation materials on the roadway surface with high surrounding rock temperature, so as to reduce heat transfer from high-temperature rock to air flow. This method, in combination with the adjustment of ventilation parameters, can achieve a good cooling effect with a relatively low cost [9]. Therefore, the research on the thermal insulation and cooling of roadways is of great significance for the control of regional thermal hazards and the efficient exploitation of resources in mines.
The development of mine thermal insulation gunite material (TIG material) is the basis for the application of thermal insulation and cooling to roadways. Mine TIG materials are required not only to feature light weight, thermal insulation, and flame retardance but also to possess high mechanical properties. These requirements make certain materials that are commonly used in construction and chemical fields unsuitable for use in mine roadways, such as form materials (EPS, XPS, phenolic foam, and polyurethane), fibrous materials (stone wool and glass wool), and advanced materials (vacuum insulation panels, gas filed panels, and aerogel) [18, 19]. Foreign scholars rarely conducted researches on TIG materials, and Chinese scholars mainly focused on inorganic mineral materials such as expanded perlite, vitrified microsphere, and vermiculite. Fly ash, a solid waste produced by the rapid oxidation of pulverized coal during thermal power generation [20], is enormously produced in China (686 million t, 2017). It features multiple holes, low density and low thermal conductivity, and is pretty economical [21, 22]. Therefore, fly ash-inorganic mineral insulation materials are advantageous in mining. The previous literature pays more attention to the development of materials and takes thermal conductivity as the evaluation standard, but the developed materials are seldom applied in the field. Li et al. [23] developed a kind of TIG material by taking vitrified microsphere, cement, sand, and aluminum powder as the raw materials. Then, through an industrial application of the material, they found that the heat release from the shotcrete roadway was reduced by 64%. By studying the composition and ratio of raw materials, Zhang [24], Yang [25], and Zhu [26] obtained thermal insulation materials with a better performance. Liu [27] and Wang [28] simulated the temperature distribution of a TIG roadway with the aid of Fluent software. Zhou [29] discussed the influence of thermophysical parameters of the TIG layer on its thermal insulation effect. It can be concluded that scholars mostly evaluated the thermal insulation effect of TIG materials in high geo-temperature roadways by means of numerical simulation. Nevertheless, since the unstable heat transfer coefficients of the roadways are difficult to determine, the simulation results differ notably from the actual situation.
Therefore, despite some previous achievements made in the development of TIG materials, there is still a lack of researches on the thermal insulation effect of mine TIG materials and the influence of TIG layers on the temperature field of surrounding rock. Physical simulation experiments provide us with solutions. Wang [30], Zhang [31], Zhu [32], and Wang [33] analyzed the temperature and humidity fields of air flow and surrounding rock with the aid of similar simulation devices, but their conclusions were all targeted at non-TIG roadways. On the basis of previous studies, the influence of the TIG layer on the geo-temperature field was investigated by using a similarity simulation device (Figure 1) in this paper. The investigation can serve to evaluate the thermal insulation effect of TIG materials. Besides, it is meaningful for the prediction of airflow and rock temperature and the adjustment and optimization of the ventilation network. Ultimately, it conduces to the control of thermal hazards in mines.

This paper is focused on the application of a kind of fly ash-inorganic mineral thermal insulation material to the high geo-temperature roadway. The specific objectives of this study are (a) to develop a kind of mine TIG material which achieves both roadway support and thermal insulation; (b) to compare and verify the thermal insulation effect of the TIG layer; and (c) to analyze the temperature field characteristics of the TIG roadway.
2. Materials and Methods
2.1. Experimental Apparatus
As displayed in Figure 2, the physical similarity simulation experimental apparatus for a high geo-temperature TIG roadway mainly consists of a surrounding rock heating system, a modelling roadway, a constant temperature and humidity ventilation system, and a data acquisition system. The apparatus has a geometric similarity constant Cl = 20, a time similarity constant Cτ = 400, and a wind speed similarity constant = 0.05. Detailed design of the experimental apparatus can be found in the literature [9].

The surrounding rock heating system, which contains a closed loop constituted by a rectifier, electric heating electrodes, an alloy heating net, copper plates, and other components, works by outputting low voltage and high current through the thyristor rectifier. With the aid of the alloy heating net, the system can heat the cylinder through stepless adjustment of output current in the range of 0–1000 A. The constant temperature and humidity ventilation system can provide air flow conditions in the temperature range of 5°C–45°C, the humidity range of 20%–90%, and a maximum air volume of 5000 m3/h. The modelling roadway is composed of six cylinders, each of which has a length of 1 m, an inner diameter of 0.8 m, and an outer diameter of 1 m. In the 0.1 m space between the inner wall and the outer wall, mica paper, alloy heating net, mica paper, and thermal insulation asbestos are arranged, respectively. Temperature sensors are arranged at the positions of 0.9 R, 1 R, 2 R, 3 R, and 4 R.
2.2. Materials of the Surrounding Rock and TIG Layer
Surrounding rock materials (the range of 1 R–4 R) and TIG layer materials (the range of 0.9 R–1 R) were poured or built into the cylinders, as shown in Figure 2.
The surrounding rock material was made of high-purity quartz sand and river sand as the fine aggregates and iron powder (to improve thermal conductivity) as the additive. The ratio and thermophysical properties are given in Table 1 [9, 31]. For the material, the thermal conductivity λ2 was 1.2 W/(m·K) and the thermal diffusivity α2 was 0.8 × 10–6 m2/s. Three layers of surrounding rock material were arranged within the range of 1 R–4 R in the cylinders in Figure 1, of which the material within 1 R–3 R was built and cured by the mould while that within 3 R–4 R was directly poured in the cylinders.
The TIG material was composed of 425 ordinary Portland cement (28 d uniaxial compressive strength ≥ 36 MPa), river sand (particle size < 1 mm), and fly ash (particle size <0.5 mm). Besides, it was added with a certain amount of vitrified microsphere (particle size = 0.5–1 mm; thermal conductivity < 0.043 W/(m·K)) for enhancing its thermal insulation performance. The TIG layer was arranged within the range of 0.9 R–1 R in the cylinders, as in Figure 2.
2.3. Experimental Scheme
2.3.1. Development of the TIG Material
The TIG material should meet the requirements of thermal insulation, light weight, and excellent load-bearing capacity. To explore the material performance under different ratios, a three factor and three-level orthogonal experiment was designed. As listed in Table 2, the three factors are fly ash, vitrified microsphere, and river sand, and the three levels are three ratios corresponding to each raw material. Since cement can improve the mechanical properties of the material significantly, it is not researched on any more here.
The nine groups of experimental samples with varying ratios were numbered 1–9, representing level combinations A1B1C1, A1B2C2, A1B3C3, A2B1C3, A2B2C1, A2B3C2, A3B1C2, A3B2C3, and A3B3C1, respectively. In each group, three 70.7 mm × 70.7 mm × 70.7 mm cubic samples and one 400 mm × 400 mm × 30 mm plate sample were prepared. During sample preparation, the materials at the corresponding ratio were added with a certain amount of water and stirred evenly. Next, the evenly stirred mixture was injected into three-groove moulds, vibrated, and compacted. After 24 h, the samples were demoulded and placed in the curing room to be cured standardly (temperature 20 ± 2°C; humidity over 95%) for 28 d.
Three basic physical and mechanical parameters, namely, apparent density, uniaxial compressive strength (UCS), and thermal conductivity, of the nine groups of samples were measured. The apparent density was obtained by measuring the mass and volume; the compressive strength was measured by using a compression testing machine (MTS/SANS CMT5105); and the thermal conductivity was measured by the method of a plate heat flux meter (DRPL-400 thermal conductivity tester produced in Xiangtan, China).
2.3.2. Application of the TIG Material to the High Geo-Temperature Roadway
To investigate the application effect of the mine TIG material to the high geo-temperature roadway, two groups of experiments were performed with the aid of the physical similarity simulation experimental apparatus. The first group aims to investigate the temperature field of surrounding rock of a non-TIG roadway, namely, the roadway without a TIG layer (see Scheme 1 in Table 3 for the working conditions), while the second group aims to explore that of a TIG roadway, namely, the roadway with a TIG layer (see Schemes 1–4 in Table 3 for the working conditions).
2.4. Method for Data Analysis
The experimental data were subjected to dimensionless processing. After the processing, the dimensionless temperature, radius, and time expressions become T = (t − tf)/(t0 − tf), R = r/r1, and F0 = α2 ∗ r/(r12), respectively. In the expressions, T, R, and F0 represent the dimensionless temperature, radius, and time, respectively; t0 is the temperature of virgin rock, °C; tf is the temperature of inlet air, °C; t is the temperature of the measuring point, °C; r is the distance from a certain point to the center of roadway, m; r1 is the inner radius of roadway, 0.1 m; α2 is the thermal diffusivity of surrounding rock, m2/s; and τ is the ventilation time, s. Because the relation between dimensionless quantities is universally applicable, the experimental data need not be calculated into the real model according to the similar constant relations.
3. Results and Discussion
3.1. The Optimum Ratio of the TIG Material
Figure 3 illustrates the apparent density, thermal conductivity, and UCS data of nine groups of compounded samples. As displayed in Figure 3, the physical and mechanical properties of samples differ notably but vary in basically consistent patterns overall. For the samples, the apparent density, thermal conductivity, and UCS lie in the ranges of 1200–1630 kg/m3, 0.119–0.181 W/(m·K), and 4.8–11.7 MPa, respectively.

Table 4 gives the range analysis results of apparent density, thermal conductivity, and UCS of the nine groups of compounded samples (the data are derived from Figure 3). It can be known from Table 4 that the three factors are vitrified microsphere (B), fly ash (A), and river sand (C) in descending order of their degrees of influence on the apparent density. The optimal level combination for apparent density is A2B3C1, i.e., fly ash (50%), vitrified microsphere (30%), and river sand (60%). The three factors are vitrified microsphere (B), river sand (C), and fly ash (A) in descending order of their degrees of influence on the thermal conductivity. The optimal level combination for thermal conductivity is also A2B3C1. The three factors are river sand (C), fly ash (A), and vitrified microsphere (B) in descending order of their degrees of influence on the UCS. The optimal level combination of 28 d UCS is A1B1C2, i.e., fly ash (40%), vitrified microsphere (10%), and river sand (80%).
From the analysis of data in Figure 3 and Table 4, the apparent density, thermal conductivity, and UCS of the samples change synchronously, except some data points that vary due to experimental errors. That is, the sample with a small apparent density generally possesses a small thermal conductivity and a low UCS. The range analysis demonstrates that vitrified microsphere remarkably reduces the weight of the material and improves the thermal insulation performance, but it notably causes damage to the load-bearing capacity of the material. Hence, fly ash and cement become important factors balancing the relationship among weight reduction, thermal insulation, and load-bearing capacity.
To determine the optimal ratio of the material, the method of efficacy coefficient was adopted. It is assumed that n indexes (three indexes in this experiment, namely, apparent density, thermal conductivity, and UCS) are examined and that the efficiency coefficient of each index is di (0 ≤ di ≤ 1). If the ith examined index is the optimal index, its efficiency coefficient is 1, i.e., di = 1. Efficiency coefficients of the other indexes are ratios of these indexes to the optimal index. Then, the total efficiency coefficient is . The greater the total efficacy coefficient is, the better the ratio is. The results are listed in Table 5 The data are derived from Figure 3).
As can be observed from Table 5, d = 0.86 (fly ash 40%, vitrified microsphere 10%, river sand 60%) in Group 1 is the largest efficiency coefficient. In other words, the ratio of Group 1 is the optimal ratio. From the perspective of material properties, the thermal insulation material with this ratio can balance roadway support and thermal insulation. With a density of 1480 kg/m3, it meets the requirement for lightweight concrete 1950 kg/m3; its thermal conductivity 0.15 W/(m·K) meets the requirement for thermal insulation material 0.23 W/(m·K); and it boasts a UCS of 11.7 MPa which exceeds the strength (10 MPa) of C10 concrete.
3.2. The Heat Reduction of the TIG Layer
The experimental data of the TIG roadway and non-TIG roadway were compared under the conditions of the virgin rock temperature °C, inlet temperature °C, and average wind speed . Figures 4 and 5 exhibit the temperature distribution characteristics of radial points (1 R–4 R) in the surrounding rock and their dynamic variations with time in the two experimental cases.


When dimensionless temperature T in the surrounding rock falls to 0.99, it is considered that thermal disturbance occurs in the surrounding rock. In the non-TIG roadway, thermal disturbance occurs at 1 R and 2 R of the surrounding rock when the dimensionless time ; it occurs at 3 R when and at 4 R when . Radial points (1 R–4 R) in the surrounding rock undergo thermal disturbance earlier for the non-TIG roadway than the TIG one.
For the non-TIG roadway wall, the initial and end dimensionless temperatures T are 0.73 and 0.51 , respectively, while for the TIG roadway wall, the initial and end dimensionless temperatures T are 0.79 and 0.6 , respectively. The temperature drop range and rate of the surrounding rock (2 R–4 R) are both greater for the non-TIG roadway than for the TIG roadway. Obviously, the TIG layer prevents heat in the deep surrounding rock from dissipating to the airflow in the roadway and plays a significant role in cooling the roadway.
To interpret the thermal insulation effect of the TIG layer, it is necessary to calculate the amount of heat reduction. In engineering, the amount of heat release (or heat absorption) from surrounding rock can be calculated in accordance with the following formula:where is the amount of heat released from surrounding rock, W; is the coefficient of unstable heat transfer between surrounding rock and airflow, ; is the perimeter of the roadway interface, ; and L is the length of the roadway, .
The value of the coefficient of unstable heat transfer is difficult to obtain by this experimental system. The analysis indicates that it is through the roadway wall that heat dissipates from the surrounding rock in the high geo-temperature roadway to the airflow. The heat loss equals the integral of heat flow rate in the unit area to the time , i.e.,
In this experiment, the homogeneous surrounding rock material has the same thermal conductivity . The slope of curve in Figure 4 numerically equals the temperature gradient of surrounding rock (sharing a linear conversion relationship with ). The shadow area in Figure 5 physically means the integral term of the temperature gradient between 1 R and 2 R of surrounding rock (approximately representing the wall) over a period of time. The ratio represents the heat reduction per unit wall in a period of time after the high geo-temperature roadway is equipped with a TIG layer. As shown in Table 6, under the experimental conditions, the heat reduction of the roadway wall is 64.2% in the initial stage of ventilation, whereas it is greatly reduced in the later stage, namely, after . Thus, the TIG layer should be constructed at the early stage of ventilation.
3.3. The Temperature Field Characteristics of the TIG Roadway
In the similarity experiment of the second group of the TIG roadway, the experimental results reveal that, after the excavation and ventilation of the TIG roadway, the temperature distributions of radial points in the surrounding rock in Schemes 1–4 develop in similar trends and so do their dynamic variations with time in Schemes 1–4, as illustrated in Figures 6–9. In Figure 6(a), Zone I is the TIG layer zone (0.9 R–1 R) and Zone II is the surrounding rock zone (1 R–4 R). Here, only the experimental data under the working conditions of Scheme 1 are analyzed.

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)
3.3.1. Relationship between Dimensionless Temperature and Dimensionless Radius
As can be observed from Figure 6(a), an obvious exponential function relationship exists between and in the surrounding rock zone (Zone II), and the expression is . That is, the temperature of surrounding rock rises gradually with the increase in depth and approximates the temperature of the virgin rock. Partial fitting parameters are displayed in Table 7. The analysis shows that and are both the functions of and (Biot number) under the condition of constant airflow state [31].
It can also be seen from Figure 6(a) that ① the temperature gradient between the TIG layer and the surrounding rock experiences an obvious abrupt change due to the great difference in thermal conductivity of materials; ② assuming that thermophysical properties inside the surrounding rock are constant, the thermal flux density of the surrounding rock near the roadway wall is larger in the early stage of ventilation; the various thermal flux densities inside the surrounding rock tend to be consistent as the ventilation time increases from 0 to 1. This indicates that the temperature field of the surrounding rock regains a balance, which marks the end of temperature disturbance.
3.3.2. Thermal Disturbance of the TIG Roadway
From Figure 6, the temperature of the surrounding rock of the TIG roadway is slowly disturbed in the beginning, and the range of disturbance increases gradually with dimensionless time .
According to the data, when , the dimensionless temperatures T of 0.9 R, 1 R, and 2 R are disturbed immediately. Among them, the temperature of the TIG layer wall (0.9 R) drops to 30.6% of the initial dimensionless temperature; that of roadway wall (1R) declines to 78.8% of the initial dimensionless temperature; and that of internal surrounding rock (2 R) begins to be slightly disturbed. The dimensionless temperatures of 3 R and 4 R in the surrounding rock do not begin to be disturbed until and , respectively.
3.3.3. Characteristics of Temperature Drop at Different Radial Positions of the TIG Roadway
The temperatures at 0.9 R and 1 R on both sides of the TIG layer drop suddenly at the beginning of ventilation, but the degrees of drops differ, leading to the formation of a relatively large temperature difference. The temperature outside the TIG layer (0.9 R) drops in an equal slope; the temperature inside the TIG layer (1 R) declines at a gradually reducing rate; and the temperature inside the surrounding rock (2 R) falls at an increasing rate first and at a decreasing rate later.
(1) TIG Layer Zone (Zone I). The temperature on the TIG layer wall falls sharply at the beginning of ventilation. As presented in Figure 10, after , the temperature here (0.9 R) decreases in an equal slope with the increase of under the four schemes.

(2) Surrounding Rock Zone (Zone II). According to the temperature data at 1 R and 0.9 R on both sides of the TIG layer, the temperature at 1 R drops sharply at the beginning of ventilation when , but the range of drop is much smaller than that of the TIG layer wall, leading to the formation of a relatively large temperature difference on the two sides. The temperature at 1 R declines at a gradually decreasing rate with the increase of and tends to level off when .
The temperatures at 2 R in the surrounding rock decreases with the increase in , but their temperature variation rates differ from that of the wall. The temperature drop rates keep growing before , and they tend to level off after , as shown in Figure 6(b). Since 3 R and 4 R are closer to the outer boundary, the abovementioned law is not evidently exhibited at the two positions. However, it can be inferred that for the position father away from the wall, the turning point of the temperature drop rate appears later temporally. This indicates that the inside of the surrounding rock exhibits a lagging response time for the adjustment of boundary temperature disturbance, and the longer the distance is, the longer the response time is. Compared with the non-TIG roadway, the TIG roadway presents more obvious lagging response time.
By comparing Figure 6 with Figures 7–9, it is found that the variations of ventilation parameters (wind speed and wind temperature) and virgin rock temperature do not alter the variation trends of the curve and curve. Besides, variations of these parameters exert no significant effect on the internal temperature field of the TIG roadway in this experimental scope. As presented in Figure 10, the variations of ventilation parameters (wind speed and wind temperature) and virgin rock temperature directly affect the wall temperature, which demonstrates that the TIG layer has a significant effect on the stability of the temperature field of internal surrounding rock.
Figure 11 is a comparison between the temperature drop and the heat reduction at 2 R of the surrounding rock under Scheme 1. From Figure 11, the heat reduction of the TIG roadway is significantly lowered after , which corresponds to the gentle temperature drop rates at 2 R, 3 R, and 4 R, as in Figure 6(b). This suggests that the TIG structure should be constructed in the early stage of ventilation in order to achieve an excellent effect.

4. Conclusions
In this paper, a kind of mine TIG material that boasts light weight, thermal insulation, and large load-bearing capacity was developed by taking cement, sand, fly ash, and vitrified microsphere as raw materials and density, thermal conductivity, and UCS as evaluation indexes. Furthermore, the thermal insulation effect of the TIG layer was analyzed with the aid of a physical similarity simulation experimental apparatus independently developed by our research team. Besides, the temperature field characteristics of the surrounding rock were researched on. The following main conclusions were drawn:(1)The thermal insulation effect of the TIG material in the high geo-temperature roadway was analyzed by means of physical similarity simulation. The temperature disturbance time, temperature drop rate, and temperature drop range of radial points in the surrounding rock of the TIG roadway and the non-TIG roadway were compared. The results show that, in the early stage of TIG layer construction, the heat reduction of the roadway wall reaches over 60%, which directly demonstrates the excellent thermal insulation performance of the TIG material. The changes in ventilation parameters (wind speed and wind temperature) and virgin rock temperature do not notably influence the internal temperature field of the surrounding rock of the TIG roadway, which suggests that the TIG layer has an important influence on the stability of the surrounding rock temperature field. This indirectly demonstrates the good thermal insulation performance of the TIG material.(2)Characteristics of temperature field in the surrounding rock of the TIG roadway were researched on. The thermal flux density of surrounding rock near the roadway wall is larger in the early stage of ventilation; then, the thermal flux densities tend to be consistent as the ventilation time increases from 0 to 1. This indicates that the temperature field of surrounding rock regains a balance, which marks the end of temperature disturbance. An obvious exponential function relationship exists between and during the entire time period.(3)Different radial points of the surrounding rock in the TIG roadway exhibit varying characteristics of temperature drop. The temperature outside the TIG layer drops sharply as ventilation begins, and then, the temperature drop develops in an equal slope. The temperature of the inner wall of the TIG layer declines in a smaller range than that of the outer wall at the beginning of ventilation, after which the temperature drop rate decreases gradually. The temperature drop inside the surrounding rock increases first and, then, decreases. The inside of the surrounding rock exhibits a lagging response time for the adjustment of boundary temperature disturbance; and the longer the distance is, the longer the response time is. Compared with the non-TIG roadway, the TIG roadway presents more obvious lagging response time. The stage of decreasing temperature drop rate of the surrounding rock corresponds to the stage of gentle heat reduction of the wall.
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 in the language proofreading. This research was funded by the National Natural Science Foundation of China, grant no. 51674242.