Performance of Heat-Insulating Materials Doped with Basalt Fibres for Use in Mines

To solve high-temperature-induced hazards in mines, heat-insulating materials were prepared by utilising basalt fibres and high-strength ceramsite combined with cementing materials. Through orthogonal tests and data analyses, the optimal combination of the heat-insulating materials doped with basalt fibres was determined as A1B1C1, that is, doping with 45% basalt fibres, a length of the basalt fibres of 6 mm, and doping with 20% ceramsite. The performance indices corresponding to the optimal comprehensive combination of the heat-insulating materials doped with basalt fibres included a density of 1200 kg/m3, thermal conductivity of 0.151 W/(mK), compressive strength of 9.7 MPa, flexural strength of 3.6 MPa, and a water-seepage depth of 25.4 mm. Numerical simulations verified that the materials presented favourable thermal insulation performance.


Introduction
Coal is the primary energy source in China, and the coal mining industry promotes the socio-economic development of the nation. With increases in mining depth, the temperature of the rock increases, triggering high-temperature hazards. These hazards affect workers' physical and mental health, and cause a reduction in their labour output; more seriously, scientific research has shown that high-temperature-induced hazards can also trigger rock collapse and gas explosion accidents, threatening mine safety [1].
At present, high-temperature hazards in mines are ameliorated by applying a mechanical cooling refrigeration system combined with cooling devices set in local areas underground; however, this method relies on the reasonable distribution and application of large-scale refrigeration equipment, which will increase the cost of mine construction. In the trend of pursuing environmental protection and energy saving benefits, it is thus deemed necessary to explore and implement new cooling measures for mines to improve the quality of such cooling systems. In this context, heat-insulating materials provide a new idea for controlling heat hazards in mines.
Heat-insulating materials refer to raw materials characterised by favourable thermal insulation performance, low material density, and suitability for application as spraying and support materials underground, aiming at high-temperature hazards in deep mines [2]. In many deep mines in South Africa, engineers have applied organic heat-insulating materials (e.g., polyethylene foam, polyurethane, and polyurethane foam) to control high-temperature hazards [3]. These materials can be used as investigated and compared with the behaviour of 100% flax and 100% basalt composites. After the optimisation of the epoxy resin curing cycle, tension-tension fatigue tests were performed on samples with two orientations, i.e., 0/90 • and ±45 • . The results showed that hybridization was able to produce a positive effect on the fatigue resistance of basalt laminates. Chlup et al. [14] characterised both representative materials with the aim to determine similarities and differences in the fracture processes. The microstructural, elastic, and fracture properties were also examined. The fracture resistance was obtained in two typical directions, i.e., along and across the fibres. Based on this, this paper studies the influence of basalt fibre length on thermal insulation materials, and determines its influence on thermal insulation and mechanical properties.
Basalt fibres feature many excellent physical properties such as low cost, high temperature resistance, stretching resistance, corrosion resistance, and safety during production. They can be degraded into soil parent materials after being abandoned and are, effectively, green industrial materials [15][16][17]. By comparing and analysing three inorganic heat-insulating materials experimentally, Jiang et al. [18] found an optimal test material blended with basalt fibres with a thermal conductivity coefficient of 0.1323 W/(mK) and a compressive strength of 10.98 MPa. The material was able to be used as a heat-insulating material in mines. Heat-insulating materials for mines were prepared by selecting basalt fibres as the main raw materials, and their mechanical properties and thermal insulation performance were experimentally explored to reveal their comprehensive superiority in preparing heat-insulating materials for mines.

Test Materials
In light of relevant stipulations in current Chinese standards for cements, 425# Portland cement produced in Qingdao, Shandong Province, China was used to satisfy requirements for sprayed concretes in terms of strength, durability and workability. The fly ash (ceramsite) was prepared by adding a certain amount of fly ash and doping a certain amount of binder and water in the mixing process of raw materials to process the mixture into particles with approximately equal sizes (the diameter of the ceramsite selected in this experiment was 3−5 mm). The ceramsite was prepared through high-temperature calcination, showing a series of advantages, e.g., it was light-weight, of high strength, prepared from abundant raw materials, environmentally-friendly, and cheap. Fly ash ceramsite produced in a construction material factory in Jinan, Shandong Province, was employed in the tests. Basalt fibres were first developed in the Soviet Union in the 20th century. At present, basalt fibres are regarded as the ideal filler for polymer composites and concretes due to their mechanical and thermodynamic performances [19]. The basalt fibres (thermal conductivity was between 0.031-0.038 W/(mK) and the breaking strength was 3500-4300 Mpa) produced by a material factory in Jinan were applied in these tests, and natural medium-grained sands were used as the test aggregates.

Design of the Basic Mixture Proportions
According to the single factor test performed initially, the basic proportions of the heat-insulating materials for mines doped with basalt fibres were set to 1:3.58:0.5 (cement: Aggregate: Water). Some aggregates were replaced by ceramsite.

Orthogonal Test Design
According to the design method of orthogonal tests, three-factor, three-level tests (a total of nine groups of tests) were designed. As shown in Table 1, based on existing research, three factors (i.e., dosage and length of basalt fibres, as well as the dosage of ceramsite) were designed; three levels represented three corresponding dosages of raw materials. The dosages of basalt fibres (percentage of basalt fibres in the total mass of cement) were 45%, 50%, and 55%, respectively; the lengths of the basalt fibres were (separately) 6, 12, and 18 mm, respectively; the dosages of ceramsite were 40%, 45%, and 50%. A total of three factors were considered, and an L 9 (3 4 ) orthogonal table was used in the present experimental design.

Preparation and Maintenance of Specimens
The proportions of raw materials in various groups were designed according to the orthogonal tests, and further raw materials were separately weighed. The high-strength ceramsite, sands, and cements were poured into a stirrer, in which basalt fibres were added after sufficient stirring. To avoid agglomeration of the basalt fibres, manual dispersion was necessary. Afterwards, the stirred concrete specimens were reinforced with basalt fibres and placed into the mould required for different tests; subsequently, the moulds were placed on a horizontal shaking table and vibrated. After preparation, the specimens with different dimensions should be placed apart according to strength class for about 24 h; next, the specimens were demoulded and labelled.
All specimens were placed into a Ca(OH) 2 saturated solution to be cured. In the curing box, the specimens were not allowed to come into contact with each other, i.e., they were separated by at least 10 mm and always submerged with the solution level over 50 mm above their upper surfaces. The specimens were cured for 28 days indoors at a temperature of (20 ± 2) • C to avoid exposure to sunlight. After 28 days, the specimens were removed from the curing tank and air-dried in shade under indoor conditions at constant temperature ( Figure 1).

Preparation and Maintenance of Specimens
The proportions of raw materials in various groups were designed according to the orthogonal tests, and further raw materials were separately weighed. The high-strength ceramsite, sands, and cements were poured into a stirrer, in which basalt fibres were added after sufficient stirring. To avoid agglomeration of the basalt fibres, manual dispersion was necessary. Afterwards, the stirred concrete specimens were reinforced with basalt fibres and placed into the mould required for different tests; subsequently, the moulds were placed on a horizontal shaking table and vibrated. After preparation, the specimens with different dimensions should be placed apart according to strength class for about 24 h; next, the specimens were demoulded and labelled.
All specimens were placed into a Ca(OH)2 saturated solution to be cured. In the curing box, the specimens were not allowed to come into contact with each other, i.e., they were separated by at least 10 mm and always submerged with the solution level over 50 mm above their upper surfaces. The specimens were cured for 28 days indoors at a temperature of (20 ± 2) °C to avoid exposure to sunlight. After 28 days, the specimens were removed from the curing tank and air-dried in shade under indoor conditions at constant temperature ( Figure 1).

Compressive Strength Testing
The compressive strengths of nine groups of specimens (nonstandard specimens measuring 100 mm × 100 mm × 100 mm) were tested. To improve the accuracy and reliability of the test data, three replicates were tested. An AGX-250 electronic universal testing machine (Shimadzu Corporation, Kyoto, Japan) was used to measure the compressive strength.

Flexural Strength Testing
According to the relevant stipulations of the standard on the test method of the mechanical properties of sprayed concretes and the test for flexural strength, the flexural strength was measured by using beam-type standard specimens measuring 150 mm × 150 mm × 550 mm. Three such specimens were designed as a group, and a total of nine groups were conducted. A YDW-10 microcomputer-controlled electric flexural testing machine (Xingaoweiye Scientific, Tianjin, China) was used.
The specimens should be used timeously after curing. The specimens were wiped and cleaned, and then installed according to "Standard for Test Methods of Mechanical Properties of Ordinary Concrete" (GB/T 50081-2002), with the deviation of the installation dimensions being no greater than 1 mm. The flank of a formed specimen was supposed to act as the load-bearing surface thereof. The pedestal and the support surface were steadily and uniformly brought into contact with the cylinder before loading at 0.02 to 0.05 MPa/s. When a specimen was almost damaged, it was necessary to stop adjusting the accelerator of the testing machine until the specimen was damaged; afterwards, the failure load was recorded. The reading on the testing machine under the failure load and the fracture position of the lower edge of the specimen were recorded.

Permeability Testing
A HP-4.0 automatic pressure regulating concrete permeability meter (Jingrui Instrument, Cangzhou, China), was used. The specimen measured 175 mm × 185 mm × 150 mm, and three replicates were tested.
The cured specimens were taken out one day before testing. The surface of the specimens was dried in the sun and the testing mould was heated to 40 • C. Moreover, after heating the sealing wax until it was completely melted, the specimens were rolled through a circle in the melted wax (without getting wax on the two end faces of the specimens, which could be protected by pasting newspapers thereon). The specimens coated with wax were pressed into the testing mould preheated by the drying oven. After being cooled slightly, the pressure could be removed and the specimens together with the testing mould were nested on the permeability meter. The nut of the watering gun was screwed off and various valves were opened; subsequently, a funnel was connected to the water-injection nozzle to inject water into the water tank (the chassis of the testing mould was also filled with water to displace air in the pipeline system). After they were filled with water, six valves feeding water to the testing mould were closed. The packed testing mould with its specimens was fixed onto the instrument. At first, the small valve was opened and the No. 0 control valve was closed. Until the linear water flow was generated from the small water nozzle, the six control valves were opened and the small water valve was closed. The initial water pressure was set as 0.1 MPa, and then water pressure was increased by 0.1 MPa at intervals of 8 h. The water seepage of the upper end face of the specimens was observed until the test was ceased when the end faces of three of six specimens showed seepage, and the water pressure at that time was recorded. The water-seepage height can be measured by fracturing the water-seeping specimens at the middle using the fracturing equipment. The impermeability is calculated as follows: where S and H refer to the degree of impermeability and water seepage pressure (MPa), respectively.

Thermal Conductivity Testing
A DRPL-1 thermal conductivity coefficient instrument (Xiangtan Instrument, Hunan, China), was applied in the test. The measuring range of instrument thermal conductivity was 0.015-5 W/(mK), and heat-insulating materials measuring 100 mm × 100 mm × 100 mm were tested using the steady-state plate method. The instrument was divided into the upper (hot) and lower (cold) surfaces. After setting the temperature, the heat was transferred to the cold surface through the specimens; afterwards, according to the thickness (100 mm) and the heat transfer area (10,000 mm 2 ) of the specimens, the thermal conductivity of the heat-insulating materials was calculated.

Analysis of the Apparent Density
Apparent density refers to the mass of material particles of unit apparent volume (including material solid volume and volume of closed pores) in a natural state, which is calculated as follows: where ρ , m and v denote the apparent density (kg/m 3 ), the mass (kg), and volume (m 3 ) of materials, respectively. It can be seen from the statistical results in Table 2 (A-dosage of basalt fibres, B-length of basalt fibres, C-dosage of ceramsite) that the apparent density of the specimens in the test groups ranged from 1170 to 1480 kg/m 3 . Compared with the sprayed concrete (with an apparent density of 2100 to 2600 kg/m 3 ) commonly used in underground coal mines, the apparent density of the heat-insulating materials was lower, showing the advantage of being light.

Data Analysis: Compressive Strength
An L 9 (3 4 ) orthogonal table was employed and three tests were performed at each factor level. K i represents the sum of the three test results at the ith level; on this basis, the mean k i was attained by dividing K i by the number of levels. After preliminary processing of the data, the result is displayed in Table 3.  To explore the trend in the compressive strength, the range results during the orthogonal tests were plotted and analysed, as shown in Figure 2.
Polymers 2020, 12, x FOR PEER REVIEW 7 of 18 Table 3. Statistics pertaining to the measured compressive strengths. To explore the trend in the compressive strength, the range results during the orthogonal tests were plotted and analysed, as shown in Figure 2. According to Table 3 and Figure 2, it can be seen that in terms of the importance, the compressive strength was significantly affected by the dosage of basalt fibres, i.e., the larger the dosage, the lower the compressive strength. The dosage of ceramsite was negatively related to the compressive strength, and the length of fibres had no obvious influence on the compressive strength.

Test Number Level Combination Experimental Variables Compressive Strength (MPa) A (%) B (mm) C (%)
The three factors influencing the heat-insulating materials were (in descending order): (A) dosage of basalt fibres, (C) dosage of ceramsite and (B) length of basalt fibres. As for the strength of the materials, the higher the compressive strength, the better the mechanical properties thereof. Therefore, we should select those factors with the maximum K value in various columns.
Owing to K1 > K2 > K3 in the column for factor A; K2 > K3 > K1 in the column for factor B, and K3 > K2 > K1 in the column for factor C, the optimal scheme combination was A1 B2 C3 in terms of the compressive strength, that is, the optimal proportion for the compressive strength of the heat- According to Table 3 and Figure 2, it can be seen that in terms of the importance, the compressive strength was significantly affected by the dosage of basalt fibres, i.e., the larger the dosage, the lower the compressive strength. The dosage of ceramsite was negatively related to the compressive strength, and the length of fibres had no obvious influence on the compressive strength.
The three factors influencing the heat-insulating materials were (in descending order): (A) dosage of basalt fibres, (C) dosage of ceramsite and (B) length of basalt fibres. As for the strength of the materials, the higher the compressive strength, the better the mechanical properties thereof. Therefore, we should select those factors with the maximum K value in various columns.
Owing to K 1 > K 2 > K 3 in the column for factor A; K 2 > K 3 > K 1 in the column for factor B, and K 3 > K 2 > K 1 in the column for factor C, the optimal scheme combination was A 1 B 2 C 3 in terms of the compressive strength, that is, the optimal proportion for the compressive strength of the heat-insulating

Data Analysis: Flexural Strength
The measured flexural strengths are listed in Table 4. To investigate the trend in the flexural strength, the range results during the orthogonal tests were drawn and analysed ( Figure 3). insulating materials was listed as follows: dosage of basalt fibres of 45%, length of basalt fibres of 12 mm, and dosage of ceramsite of 60%.

Data Analysis: Flexural Strength
The measured flexural strengths are listed in Table 4. To investigate the trend in the flexural strength, the range results during the orthogonal tests were drawn and analysed ( Figure 3). As shown in Table 4 and  As shown in Table 4 and properties of the developed materials. Thus, it was necessary to select the factors with the maximum K value in each column.
Due to K 3 > K 2 > K 1 in the column for factor A; K 2 > K 3 > K 1 in the column for factor B, and K 1 > K 2 > K 3 in the column for factor C, A 3 B 1 C 2 were regarded as the optimal scheme combination in terms of the flexural strength, that is, the optimal proportion for the flexural strength of the heat-insulating materials was described as follows: doping with 55% basalt fibres, using basalt fibres with a length of 6 mm, and doping with 40% ceramsite.

Data Analysis: Impermeability
The impermeability test was conducted on the basis of the water-seepage height. The test data are presented in Table 5. To investigate the trend in the impermeability, the range results found during orthogonal tests were plotted and analysed, as shown in Figure 4.
Polymers 2020, 12, x FOR PEER REVIEW 9 of 18 developed materials. Thus, it was necessary to select the factors with the maximum K value in each column. Due to K3 > K2 > K1 in the column for factor A; K2 > K3 > K1 in the column for factor B, and K1 > K2 > K3 in the column for factor C, A3 B1 C2 were regarded as the optimal scheme combination in terms of the flexural strength, that is, the optimal proportion for the flexural strength of the heat-insulating materials was described as follows: doping with 55% basalt fibres, using basalt fibres with a length of 6 mm, and doping with 40% ceramsite.

Data Analysis: Impermeability
The impermeability test was conducted on the basis of the water-seepage height. The test data are presented in Table 5. To investigate the trend in the impermeability, the range results found during orthogonal tests were plotted and analysed, as shown in Figure 4.  According to Table 5 and Figure 4, an increase in the dosage of ceramsite, as the aggregate in the material, would increase the void ratio in the finished concrete. Although this meant that cementing materials were combined with the aggregates to form a skeletal void structure to some extent, contributing to certain enhancement of the strength of the materials, it significantly affected the permeability of the materials. Due to the growth of the dosage of basalt fibres, a crossed network structure was formed in the heat-insulating materials, which blocked the pore water and prevented capillary water flow. Therefore, increasing the dosage of the basalt fibres in the mixture within a certain range was conducive to improving the permeability of the material. It can be seen that in terms of the importance, the three factors influencing the impermeability of the heat-insulating materials were (in descending order) as: (C) dosage of ceramsite, (A) dosage of basalt fibres, and (B) length of basalt fibres. It was better to have a lower permeability of materials. Thus, it was necessary to select the factors with the minimum K value in each column.
Due to K 3 < K 2 < K 1 in the column for factor A; K 3 < K 1 < K 2 in the column for factor B, and K 1 < K 2 < K 3 in the column for factor C, the optimal scheme combination was A 3 B 3 C 1 in terms of the impermeability, that is, the optimal proportion for the impermeability of the heat-insulating materials was as follows: doping with 55% basalt fibres, using basalt fibres with a length of 18 mm, and doping with 20% ceramsite.

Data Analysis: Thermal Conductivity
The measured thermal conductivity is listed in Table 6. Table 6. Statistics pertaining to the thermal conductivity. To investigate the trend in the thermal conductivity, the range results during the orthogonal tests were plotted and analysed, as shown in Figure 5.

Experimental Variables Thermal Conductivity (W/(mK)) A (%) B (mm) C (%)
As shown in Table 6, the ranges among three columns for factors A, B, and C were different, implying that the change of the three factor levels exhibited different influences on the thermal conductivity. According to the importance, the three factors are (in descending order): (A) (dosage of basalt fibres), (B) (length of basalt fibres), and (C) (dosage of ceramsite).
According to Figure 5, it can be found that the thermal conductivity varied significantly with the change in the dosage of basalt fibres. There was a negative correlation between the two, and the dosage presented a significant influence on the index. From the perspective of the length of basalt fibres, the thermal conductivity was low when the length of the selected materials was 6 mm. As the length of the basalt fibres varied from 6 mm to 18 mm, the thermal conductivity coefficient gradually rose. In the stirring process, it was difficult to mix the basalt fibres with the other materials after their length exceeded a certain value, and therefore, the mixture failed to be sufficiently stirred. Through analysis according to the change of the dosage of ceramsite in the third column, it was found that the thermal conductivity coefficient was marginally changed after ceramsite was added into the developed materials. Thus, ceramsite exerted a slight influence on the thermal conductivity of the heat-insulating materials. As shown in Table 6, the ranges among three columns for factors A, B, and C were different, implying that the change of the three factor levels exhibited different influences on the thermal conductivity. According to the importance, the three factors are (in descending order): (A) (dosage of basalt fibres), (B) (length of basalt fibres), and (C) (dosage of ceramsite).
According to Figure 5, it can be found that the thermal conductivity varied significantly with the change in the dosage of basalt fibres. There was a negative correlation between the two, and the dosage presented a significant influence on the index. From the perspective of the length of basalt fibres, the thermal conductivity was low when the length of the selected materials was 6 mm. As the length of the basalt fibres varied from 6 mm to 18 mm, the thermal conductivity coefficient gradually rose. In the stirring process, it was difficult to mix the basalt fibres with the other materials after their length exceeded a certain value, and therefore, the mixture failed to be sufficiently stirred. Through analysis according to the change of the dosage of ceramsite in the third column, it was found that the thermal conductivity coefficient was marginally changed after ceramsite was added into the developed materials. Thus, ceramsite exerted a slight influence on the thermal conductivity of the heat-insulating materials.
According to the test results, due to K1 > K2 > K3 in the column for factor A; K3 > K2 > K1 in the column for factor B, and K1 > K2 > K3 in the column for factor C, the optimal scheme combination for the thermal conductivity was A3B1C3, that is, the optimal proportion of the heat-insulating materials was as follows: the optimal dosage of the basalt fibres was 55%, the optimal length of the basalt fibres was 6 mm, and the dosage of ceramsite was 60%.

Comprehensive Analysis: The Efficacy Coefficient Method
By measuring the compressive strength, flexural strength, degree of impermeability, and thermal conductivity of these materials, the optimal scheme combinations in terms of the compressive strength, flexural strength, degree of impermeability, and thermal conductivity were A1B2C3, A3B1C2, A3B3C1, and A3B1C3, respectively. To determine the optimal mixture design, these were also analysed by applying the efficacy coefficient method.
The efficacy coefficient method is a method for solving the multiobjective decision problem [20]. At first, the number n of the final indices for data analysis in the design test was determined, with the total of five indices (n = 5) in the test: apparent density, compressive strength, flexural strength, degree of impermeability, and thermal conductivity coefficient of materials. Among the nine groups of test data for a single index, the efficacy coefficient of an optimal combination for a performance According to the test results, due to K 1 > K 2 > K 3 in the column for factor A; K 3 > K 2 > K 1 in the column for factor B, and K 1 > K 2 > K 3 in the column for factor C, the optimal scheme combination for the thermal conductivity was A 3 B 1 C 3 , that is, the optimal proportion of the heat-insulating materials was as follows: the optimal dosage of the basalt fibres was 55%, the optimal length of the basalt fibres was 6 mm, and the dosage of ceramsite was 60%.

Comprehensive Analysis: The Efficacy Coefficient Method
By measuring the compressive strength, flexural strength, degree of impermeability, and thermal conductivity of these materials, the optimal scheme combinations in terms of the compressive strength, flexural strength, degree of impermeability, and thermal conductivity were A 1 B 2 C 3 , A 3 B 1 C 2 , A 3 B 3 C 1 , and A 3 B 1 C 3 , respectively. To determine the optimal mixture design, these were also analysed by applying the efficacy coefficient method.
The efficacy coefficient method is a method for solving the multiobjective decision problem [20]. At first, the number n of the final indices for data analysis in the design test was determined, with the total of five indices (n = 5) in the test: apparent density, compressive strength, flexural strength, degree of impermeability, and thermal conductivity coefficient of materials. Among the nine groups of test data for a single index, the efficacy coefficient of an optimal combination for a performance index was set to 1, that is, d i = 1. The efficacy coefficient of any other test combination was equal to the ratio of test data of the combination to those of the optimal combination. The total efficacy of each test group was calculated using Equation (3) and the proportion of the test group with the maximum total efficacy coefficient was considered to be the optimal mixture design.
As shown in Table 7, the total efficacy coefficient (0.83) of the first group of data was found to be the largest after processing the nine groups of data. Therefore, the optimal combination was A 1 B 1 C 1 , that is, doping with 45% basalt fibres, using basalt fibres with a length of 6 mm, and doping with 20% ceramsite. The corresponding performance indices of this combination of the heat-insulating materials were listed as follows: an apparent density of 1200 kg/m 3 , a thermal conductivity of 0.151 W/(mK), a compressive strength of 9.7 MPa, a flexural strength of 3.6 MPa, and a water-seepage height of 25.4 mm, respectively.

Simulation of the Thermal Insulation
The related simulations of the heat-insulating materials doped with basalt fibres were conducted by applying COMSOL Multiphysics software [21]. By comparing the simulation effects of the heat-insulating materials and common concretes with a same thickness, the thermal insulation effect of the heat-insulating materials in high-temperature mine roadways was verified.
A model was established based on the practical parameters of an auxiliary uphill haulage-way of a certain mine at the depth of 900 m underground. The specific parameters are listed in Table 8. A physical simulation model was established based on the aforementioned scheme design and parameters; thereafter, grid partition was applied to the established physical model for the roadway. Not every domain is shown due to a large aspect ratio, and therefore, the vertical coordinate (length direction of the roadway) is lowered by a factor of 10 to show the geometric model ( Figure 6). In the process of grid partition, it is necessary to refine the grids of the roadway, as they are the most important part of the mesh. Corresponding grids for the roadway are further refined; afterwards, the grid quality is examined to ensure that all elements satisfy the computational demand imposed during subsequent simulation. parameters; thereafter, grid partition was applied to the established physical model for the roadway. Not every domain is shown due to a large aspect ratio, and therefore, the vertical coordinate (length direction of the roadway) is lowered by a factor of 10 to show the geometric model ( Figure 6). In the process of grid partition, it is necessary to refine the grids of the roadway, as they are the most important part of the mesh. Corresponding grids for the roadway are further refined; afterwards, the grid quality is examined to ensure that all elements satisfy the computational demand imposed during subsequent simulation. A numerical simulation was performed on the roadway by separately taking the common sprayed concretes and sprayed heat-insulating materials A1B1C1 as supporting materials in the ventilation area. The nephograms of the temperature field of the two materials are displayed in Figure 7. A numerical simulation was performed on the roadway by separately taking the common sprayed concretes and sprayed heat-insulating materials A 1 B 1 C 1 as supporting materials in the ventilation area. The nephograms of the temperature field of the two materials are displayed in Figure 7.
For common concretes and heat-insulating materials A 1 B 1 C 1 (when sprayed to the same thickness), the external boundary conditions were consistent and the temperature of the airflows at the air inlet was the same. In Figure 7a, the temperature of surrounding rocks close to the profile of the roadway changed: The thermal conductivity of the common sprayed concretes used as a support material was relatively high; the heat dissipation in surrounding rocks would exert a larger influence on the airflow temperature when cold airflows constantly went through the roadway, and more cold energy was used to balance the heat dispersed from the surrounding rocks. Therefore, more heat was dispersed from surrounding rocks to the roadway through the common sprayed concrete (support materials) per unit time, which led to an increased airflow temperature. Due to the gradual increase of the airflow temperature, the regulating effect of airflows on surrounding rocks along the direction of the length of the roadway was gradually weakened and their zone of influence also shrank. In Figure 7b, the characteristic of the nephogram formed due to the heat exchange between airflows and surrounding rocks when air flowed through the roadway was less significant than that in Figure 7a. The main reason for this was that the thermal conductivity (0.151 W/(mK)) of the heat-insulating materials for mines doped with basalt fibres was relatively low. This impaired the dissipation of heat from surrounding rocks to the roadway, thus improving the efficiency of utilisation of the cold energy carried by the airflow, so that the heat-insulating materials exerted a significant thermal-insulating effect.
It can be seen from Figure 8a that the area of the circles of the same colour gradually shrank through analysis according to the temperature distribution of three transverse sections at the air inlet, middle part, and air outlet. This indicated that the temperature difference between airflows and surrounding rocks of the roadway gradually decreased as the surrounding rocks constantly transferred heat to the interior of the roadway; therefore, the effect of thermal insulation of the airflows on surrounding rocks was gradually diminished. In Figure 8b, although the change was relatively insignificant, it was evident that there was a similar trend throughout. For common concretes and heat-insulating materials A1B1C1 (when sprayed to the same thickness), the external boundary conditions were consistent and the temperature of the airflows at the air inlet was the same. In Figure 7a, the temperature of surrounding rocks close to the profile of the roadway changed: The thermal conductivity of the common sprayed concretes used as a support material was relatively high; the heat dissipation in surrounding rocks would exert a larger influence on the airflow temperature when cold airflows constantly went through the roadway, and more cold energy was used to balance the heat dispersed from the surrounding rocks. Therefore, more heat was dispersed from surrounding rocks to the roadway through the common sprayed concrete (support materials) per unit time, which led to an increased airflow temperature. Due to the gradual increase of the airflow temperature, the regulating effect of airflows on surrounding rocks along the direction of the length of the roadway was gradually weakened and their zone of influence also shrank. In Figure 7b, the characteristic of the nephogram formed due to the heat exchange between airflows and surrounding rocks when air flowed through the roadway was less significant than that in Figure 7a. The main reason for this was that the thermal conductivity (0.151 W/(mK)) of the heat-insulating materials for mines doped with basalt fibres was relatively low. This impaired the dissipation of heat from surrounding rocks to the roadway, thus improving the efficiency of utilisation of the cold energy carried by the airflow, so that the heat-insulating materials exerted a significant thermal-insulating effect.
It can be seen from Figure 8a that the area of the circles of the same colour gradually shrank through analysis according to the temperature distribution of three transverse sections at the air inlet, middle part, and air outlet. This indicated that the temperature difference between airflows and surrounding rocks of the roadway gradually decreased as the surrounding rocks constantly transferred heat to the interior of the roadway; therefore, the effect of thermal insulation of the airflows on surrounding rocks was gradually diminished. In Figure 8b, although the change was relatively insignificant, it was evident that there was a similar trend throughout. Through comparison, it can be seen that as support materials, these new heat-insulating materials present favourable thermal insulation at the same sprayed thickness. The temperature change in the materials can be analysed. Figure 9 shows the nephograms of the change of the internal temperature fields of the common sprayed concrete and heat-insulating material A1B1C1 containing basalt fibres. Through comparison, it can be seen that as support materials, these new heat-insulating materials present favourable thermal insulation at the same sprayed thickness.
The temperature change in the materials can be analysed. Figure 9 shows the nephograms of the change of the internal temperature fields of the common sprayed concrete and heat-insulating material A 1 B 1 C 1 containing basalt fibres.
(b) Figure 8. Comparison of effects of thermal insulation of support materials (longitudinal section): (a) Nephogram of the temperature field of the roadway supported using common sprayed concrete; (b) Nephogram of the temperature field of the roadway supported by using the new heat-insulating materials with basalt fibres.
The temperature change in the materials can be analysed. Figure 9 shows the nephograms of the change of the internal temperature fields of the common sprayed concrete and heat-insulating material A1B1C1 containing basalt fibres.
(a) (b) Figure 9. Nephograms of the internal temperature field of support materials: (a) Nephogram of the internal temperature field of the common sprayed concretes; (b) Nephogram of the internal temperature field of the heat-insulating materials doped with basalt fibres.
As shown in Figure 9, the temperature distribution at different positions of the materials was irregular due to being influenced by the airflow temperature at the air inlet; however, the internal temperature of materials along the airflow direction from the air inlet (1000 m) to the air outlet (0 m) gradually increased. In addition, owing to the temperature of surrounding rocks being relatively high and that of the airflows being relatively low, heat exchange occurred between the surrounding rocks and airflows in the roadway through the materials. This caused that the temperature of the heat-insulating materials to fall towards the interior of the roadway on the transverse section of the materials. Compared with the heat-insulating materials A 1 B 1 C 1 with basalt fibres, the common sprayed concrete had a higher thermal conductivity and its heat exchange with the airflows was greater. Thus, the temperature increased rapidly along the air flow. As shown in Figure 9, the thermal insulation performance of the heat-insulating materials A 1 B 1 C 1 was superior to that of the common sprayed concrete. It can be also seen from Figure 10 that the internal temperature of the common sprayed concrete varied between 290 K to 293.2 K, and that of the heat-insulating materials varied between 293.2 K and 295.2 K. These ranges also reflected (to some extent) that the barrier performance of the heat-insulating materials for heat exchange was better than that of the common sprayed concrete.
A straight line perpendicular to the centreline of the roadway was separately drawn at the inlet (1000 m), middle (500 m), and outlet (0 m) of the model. The temperature data were subjected to statistical analysis along the straight lines from the end of the lines close to the wall surface of the roadway. The statistical results of the temperatures of the common sprayed concrete and heat-insulating materials doped with basalt fibres are displayed in Figure 11.
insulation performance of the heat-insulating materials A1B1C1 was superior to that of the common sprayed concrete. It can be also seen from Figure 10 that the internal temperature of the common sprayed concrete varied between 290 K to 293.2 K, and that of the heat-insulating materials varied between 293.2 K and 295.2 K. These ranges also reflected (to some extent) that the barrier performance of the heat-insulating materials for heat exchange was better than that of the common sprayed concrete. of the heat-insulating materials for heat exchange was better than that of the common sprayed concrete.
(a) (b) Figure 10. Changes of the internal temperatures of support materials along their centreline: (a) Change of the internal temperature of the common sprayed concrete; (b) Change of the internal temperature of the heat-insulating materials doped with basalt fibres.
A straight line perpendicular to the centreline of the roadway was separately drawn at the inlet (1000 m), middle (500 m), and outlet (0 m) of the model. The temperature data were subjected to statistical analysis along the straight lines from the end of the lines close to the wall surface of the roadway. The statistical results of the temperatures of the common sprayed concrete and heatinsulating materials doped with basalt fibres are displayed in Figure 11. It can be seen from Figure 11 that the temperatures at the air inlet and outlet of the roadway were separately the lowest and highest. The temperature was about 314 K at the position with a horizontal coordinate of 50 m, that is, the external surface of the model. By taking the temperature at the air outlet as an example, as shown in Figure 11a, the heat from surrounding rocks was gradually transferred to the interior of the roadway along the direction from 50→0 (horizontal coordinate, m), and the whole curve smoothly declined, to about 291.7 K at the wall surface (2 m) of the roadway. The reduction in temperature slowed under the influence of the sprayed concrete. In Figure 11b, the heat from surrounding rocks was also transferred to the interior of the roadway along the direction from 50→0; however, when the heat was transferred to the wall surface (2 m) of the roadway, the temperature of surrounding rocks decreased linearly to about 290.2 K due to the presence of the heat-insulating materials doped with basalt fibres, showing a temperature difference of 1.2 K compared to the common sprayed concrete at the air outlet. This suggested that the heat-insulating materials exhibited satisfactory thermal insulation.

Conclusions
(1) By performing orthogonal tests and data analyses, the optimal proportion of materials for a single index was separately analysed. The apparent densities of the nine groups of specimens were between 1170 and 1480 kg/m 3 , indicating that the heat-insulating materials doped with basalt fibres were light; the optimal scheme combinations for the compressive strength, degree of impermeability, and flexural strength were A 1 B 2 C 3 , A 3 B 3 C 1 , and A 3 B 1 C 2 , respectively. This provides a basis for further multi-index comprehensive analyses of heat-insulating materials doped with basalt fibres. (2) The thermal conductivity of the materials gradually decreased with increasing proportion of basalt fibres therein, the dosage of which greatly influenced the indices assessed here. Basalt fibres with three different specifications were selected for testing. The thermal conductivity was low when the length of the selected materials was 6 mm. The test result showed that shorter basalt fibres were more beneficial in improving the thermal insulation performance of the materials. By analysing the change of the dosage of ceramsite, it can be seen that the thermal conductivity changed slightly with the increase in the dosage of ceramsite. (3) By comparing the simulations of the heat-insulating materials doped with basalt fibres and common sprayed concrete in terms of their relative cooling effects, it can be found that the difference of the internal temperatures was 3 K and the temperature difference at the air outlet was 1.5 K. This verifies the effect of thermal insulation of materials doped with basalt fibres.
Author Contributions: Y.J. and S.X. conducted the main work and wrote the paper; H.L. and L.Z. completed basic experiments; C.H. completed the data processing and analysis; Z.Z. and J.G. designed the numerical simulation experiment and operation; they all provided insightful suggestions and revised the paper. All authors have read and agreed to the published version of the manuscript.
Funding: This work has been funded by the National Natural Science Foundation of China, grant numbers 51774197.

Conflicts of Interest:
The authors declare no conflict of interest.