The effectiveness of grouting scheme has been simulated to prevent the coal spontaneous combustion at a goaf in Haizi Colliery, China. The colliery has been operated for long period over 27 years and has a complex ventilation network including airflow leakages which could possibly lead to the spontaneous combustion of coal at goafs. Firstly, the mine ventilation simulator MIVENA was used to analyze the mine ventilation network airflows to control airflows in and out of working faces and goafs. As the second approach, numerical simulations were carried by the simulator FLUENT in order to predict spontaneous combustion of residual coal with leakage flow in the #3205 goaf. It was cleared that the goaf can be divided into three zones based on oxygen concentration in the goaf area. Finally, the numerical simulation results show that the slurry grouting method is able to be an effective and economical method by reducing porosity in the goaf area to prevent spontaneous combustion of residual coal.

1. Introduction

It is well known that the key area relating to underground coal mine fire is goaf areas where there exists a potential to induce spontaneous combustion of residual coal in the goaf areas. According to the statistics of Chinese coal production, 90% of total production coal is produced by underground mines, and there are around 600 state-owned key coalmines of which 25.1% are highly gassy mines and 17% are gas out-burst mines, approximately [1]. In Chinese underground coal mines, serious accidents occurred especially in last ten years. Since they produced a large amount of coal meets domestic demand, it resulted in insufficient effort to keep safety operations at working faces and goaf areas. In particular, most of the accidents were directly or indirectly related to coal spontaneous combustion in goaf areas.

Qi and Chen (2009) [2] studied the distributing condition of abutment pressure and air flowing state in goaf and then presented the theoretical formula to predict airflow velocity in goaf area and the spontaneous combustion zone where airflow velocity is less than 0.24 m/min. Pan et al. (2013) [3] analyzed a goaf and divided it into the three zones related to coal spontaneous combustion using with distribution of O2 and CO concentration. They classified the goaf area to radiating and asphyxiation zones by O2 concentrations larger than 18% and less than 7%, respectively.

In this research, a preventing method for spontaneous combustion of residual coal at a goaf in Haizi Colliery has been studied, because it is a typical lignite coal mine in China and has a potential to get spontaneous combustion due to fusain existing in coal beds [4].

As shown in Figure 1, the Haizi Colliery’s mine area is nearly 6.44 km2 consisting of 1,200 to 5,200 m along the strike and 0 to 1,700 m in tendency. Its coal seams are cleft by three normal faults that formed an isolated block area like a triangle coal-bearing zone. The current level that operates working faces is −320 m from the surface. Average self-ignition point of the produced lignite coal is 275°C. And the elapsed time after contacting air to spontaneous combustion in a pail has been estimated as roughly 1 to 1.5 month [5]. Based on the mine’s operation record during 2002 to 2013, spontaneous combustion phenomenon has been frequently observed at working faces. The characteristics of Haizi coal are listed in Table 1.

In this study, as the first step, authors have analyzed ventilation network airflows in the mine using the simulator MIVENA [6], especially pressure drop at a working face to control leakage airflow into a goaf area. The FLUENT was also used to simulate leakage flow in the goaf area with comparison monitoring data by the tube bundle gas monitoring system. Finally the effectiveness of slurry grout injection into the goaf area has been investigated as an economical method preventing spontaneous combustion of residual coal in the goaf.

2. Mine Ventilation Network Analysis

In the Haizi Colliery, the centralized ventilation system is used with a main intake shaft for fresh ventilation air and an exhausting shaft utilized only for exhausting airflow. The intake shaft is 6 m in diameter (28.3 m2 in sectional area) and also is used for coal elevating and hoisting workers, while the exhaust shaft is 9.6 m2 in sectional area. Two main centrifugal fans, of which one is running and other is standby, were equipped to provide adequate ventilation capabilities. A 220 kW motor operates each fan with 740 rpm in rotating speed at an acceptable efficiency. The total intake airflow = 3564 m3/min operated by total pressure loss  Pa.

In the present study, the mine ventilation network simulator MIVENA [6] has been applied to predict airflows and ventilation characteristics based on the network data (see Figure 2). In particular, the airflow around the #3205 working face, that is fully-mechanized coal cutting faces [7], was focused by changing conditions with working face operation. After that, the simulated data was used to control airflows in the working face considering whole ventilation airflows in the network.

Based on the Chinese safety regulations [8] for coal mine ventilation, airflow velocity must satisfy the range of 0.25 m/s to 4 m/s. Herein, Haizi Colliery evaluated air quantity through the goaf was 480 m3/min that meets the regulation mentioned above, because the air leakage was estimated as 368 m3/min (see Figure 2). The problem which was found was that many bypass or shortcut airways are existing at the bottom of heading with relatively large ventilation resistance in the left flank. The numerical simulations were carried for the ventilation cases with applying several useful methods to solve the problem. Based on the simulation results, the effective solution was to add an additional intake airway newly excavated in this area, because the new airway can contribute to decrease of 30% of airflow resistance at the #3205 working face 5 m in width. It was predicted that the pressure drop decreases from 319 Pa to 270 Pa at the face; therefore it was also expected to be easier to control airflow leakage into the goaf area than before.

By using the procedure mentioned above, air leakage can be reduced to 200 m3/min from 338 m3/min at the #3205 working face. On the other hand, at working faces after excavation (such as #2207 and #1709 working faces), stopping walls were built immediately to prevent air leakage into goaf areas and to improve effectiveness of fresh-air utilization in the mine. After finishing optimizing mine ventilation airflows and pressure drop around working faces, authors focused on analyzing the goaf area and preventing spontaneous combustion of residual coal in the goaf.

3. Three Zones at Goaf

With the working face advancing, the state of goaf area changing step by step. Air leakage at the working faces supplies oxygen to the residual coal in the goaf area and coal seams close to the goaf, and it may lead coal spontaneous combustion in the goaf area classified into three zones that are radiating, oxidized, and asphyxiation zones as shown in Figure 3 [913]. These three zones can be characterized by flow velocity of air leakage (), oxygen concentration in the goaf area, and rate of temperature rise that are evaluated by the method defined in the Coal Mine Safety Regulation of China [8] as in Table 2.

3.1. Oxygen Concentration Detection in Goaf

The tube bundle gas monitoring system, shown in Figure 4, has been used to monitor the coal spontaneous combustion process in the goaf, continuously. It can detect the early stage of coal spontaneous combustion based on monitoring oxidation products, such as CO, and then make an alarm to the mine engineers [1416]. At the #3205 face, gas sample is sucked into the tube bundle with an air pump. Numerous gas contents, such as O2, N2, CO, CH4, CO2, C2H4, C2H6, and so forth, are measured to detect the spontaneous combustion of coal at the monitoring spot. Gas concentration changes and distribution in the goaf area were used to make division into three zones in the goaf area, and its pros and cons were discussed to control the spontaneous combustion of coal [17].

About 50 gas monitoring points are installed in the underground, especially in the goaf areas by the monitoring system shown in Figure 5. Using continuous monitoring of gas compositions during three months, it was analyzed to pick up useful information on detecting spontaneous combustion.

The measured distribution of oxygen concentration in the goaf area is shown in Figure 6 where the horizontal axis is the distance along the #3205 working face from its upstream corner and the vertical axis is distance from working face in the goaf. The contour lines of oxygen concentration were drawn by cubic spline function fitting the limited monitoring data (around 30). There is a relatively higher temperature range in the goaf area where O2 concentration was over 10% compared with others of which oxygen concentration is not up to this critical circumscription. Therefore, the goaf area can be divided into two zones of which boundary is the contour line of 10% concentration of oxygen. The determination of boundary line between the oxidized and asphyxiation zones in the goaf was relatively easy using the oxygen concentration curves, while that of the radiating zone boundary was difficult due to high oxygen concentration range close to 21% of air. The oxidized zone of which O2 concentration is over 10% shows a possibility to start spontaneous combustion in a few days [18]; however, the other zones showed a stable condition in temperature.

It has been confirmed from those results that measurements of oxygen concentration in a goaf are the most effective method to evaluate the potential of spontaneous combustion of residual coal.

3.2. Goaf Air Flow Simulation

The numerical simulator FULENT is an effective tool to investigate spontaneous combustion of coal in the goaf area, because monitoring data are limited to reproduce completely its heating process because spontaneous combustion is comprehensible phenomena including many physical and chemical parameters. As described in previous chapter, the radiating zone boundary is difficult to determine by using the oxygen distribution curve. In this research, the numerical simulation by FLUENT [19, 20] was taken into account to assist in drawing the boundary.

In the model for present 2-dimensional numerical simulations by FLUENT, the goaf was assumed as a zone of porous media composed filled with collapsed rock and unmined coal. At first, values of goaf porosity and differential pressure are required to simulate the airflow leakage along the working face. In the present study, uniform permeability was given for the gob area based on the measured data [21] and the observations at the mine site. The goaf average porosity was set as 0.3 in goaf according to the empirical formula and measured data [22]. On the other hand, oxygen concentration of 10% O2 was used for compacted area based on the gas monitoring data recorded by the tube bundle system. The airflow is leaked through a working face into goaf area [23]. Furthermore, airflow in the goaf was treated as laminar flow expressed by Darcy’s law [24]. Based on the monitoring data, the numerical simulation results by the ventilation simulator MIVENA, inlet airflow, and its velocity into the #3205 face were evaluated as 852 m3/min and 1.3 m/s, respectively. Meanwhile, the differential pressure along the working face was simulated as 319 Pa (see Figure 2).

The numerical simulations on airflow leakage distribution in goaf area were carried out using with the above airflow conditions by the FULENT using the grid block size of 0.5 m × 0.5 m as shown in Figure 7.

In the monitoring result shown in Figure 6, the area where air leakage velocity is over 0.001 m/s corresponding to 7% O2 is located 97 meters away from the #1011 inlet roadway and #2205 outlet airway, while this area was 79 m away from the middle of #3205 working face. The oxidation zone boundary is located at 74 to 92 m far from the working face. The pressure drop in the #3205 working face that installed coal cutting machine, supports, and conveyors is much higher than that of the main airways. The measured pressure loss was 319 Pa along the working face 134 m in length. Airflow in the working face was treated as turbulent flow, because airflow velocity and roughness were large enough. In the FULENT, the standard - model which is a semiempirical model based on model transport equations for the turbulence kinetic energy () and its dissipation rate () was chosen to simulate the airflow in the working face.

Figure 8 shows the simulation result of the airflow velocity distribution in the goaf along the #3205 working face. The flow quantity of air leakage into the goaf is increasing with the pressure drop along the working face. Thus ventilation control along the working face can reduce air leakage providing oxygen to residual coal in the goaf.

As shown in Figure 9, air leakage is observed at the zone located 20 m away from the #3205 working face. And velocity is reduced with increasing distance from the 3205 working face. The working face 134 m in length was modeled between #1011 inlet and #2205 outlet airways.

The position showing the minimum airflow leakage velocity of 0.001 m/s is located at the middle of #3205 working face. Figure 9 also indicates a larger amount of air leakage and higher velocity at inlet and exhaust ends in the #3205 working face. In other words, the air leakage velocity equal to 1.2 m/min (=0.02 m/s) boundary is 26 meters away from #1011 inlet roadway and #2205 outlet airway. Moreover, the gap between the boundary and middle part of #3205 working face is 10 meters. Thus, installing air deflectors at upper-corner and down-corner angles are important to reduce the air leakage and its potential of spontaneous combustion of residual coal in the goaf.

The calculated width of oxidized zone boundary in the #3205 working face is approximately 5 to 21 m. Then, according to the oxygen concentration curves in goaf presented in Figure 6, the air leakage velocity of 0.001 m/s shows approximately 7% O2. The deep oxidation zone boundary is located at 74 to 92 m far from the #3205 working face. Figure 10 shows the farthest boundary of the oxidation zone with 10% O2 is located at 80 meters away from the #3205 working face. However, the nearest one which is next to the radiating zone is about 30 m far from the #3205 working face. By drawing the boundary lines, it is clear that the width of oxidation zone becomes narrower with increasing distance from both ends of the working face. At the same time, the radiating zone expands gradually from the center to both sides of the working face. The maximum temperature zone is located at 21 m away from the working face. These characteristics of three zones are important to design the safety working face and to prevent coal mine fires.

By combining the simulation results of airflows and the tube bundle gas monitoring results, the goaf area can be divided easily into three zones as shown in Figure 10.

4. Preventing Method for Goaf Fire

4.1. Injecting Grouting Slurry into Goaf Area

Some effective methods have been studied to prevent fires in the #3205 working face, such as nitrogen filling, grouting, gelatin injecting, and inhibitor spraying. Their costs and effectiveness were evaluated based on the current conditions of the coal mine, and then the grouting filling method using pulp slurry was chosen as the best method to prevent spontaneous combustion of residual coal in the goaf area. In the grouting filling system, the pulp slurry was transferred by natural hydraulic head between the surface and goafs. The grouting routes are as follows.

Filling station on the surface → Ventilation shaft → Return airway → Concentrated return airway in working face → Return airway in working face → Goaf area.

The slurry pool 200 m3 in volume was constructed on the level +32.7 m. The pipes to transport the slurry were vinyl tubes 2.5 inches in diameter. Total length of the slurry pipeline is  m and the vertical level difference is  m; thus the slurry line coefficient is given by = 6.03.

4.2. Grouting Method and Result

The goaf area was filled with grouting slurry to prevent the coal spontaneous combustion of residual coal by installing the pipe end at the position 40 m away from both ends of the #3205 working face (see Figure 11). And the distance between the filling area and the working face was controlled to be longer than 10 m, because grouting slurry moved from 15 to 20 m toward the working face. When a spontaneous combustion is found in the goaf, the working face and goaf area must be spilled with slurry immediately. According to Chinese regulation [25], the slurry spill must cover the area 5 m away from the working face in the goaf area. To operate grouting slurry to be effective, required amount of solid pulp is given bywhere is amount of solid pulp (m3); is coal seam thickness ( = 4.36 m); is grouting zone strike length ( = 1058 m); is coal recovery ratio, 86%; is grouting zone in clination length ( = 133.7 m); is grouting coefficient (generally 0.1 to 0.5), ( = 0.3).

Also based on the regulation [23], the required amount of water is given bywhere is amount of water, (m3); is coefficient of water reserve (=1.10 to 1.25); is water and soil ratio (= 2 to 5, average ≈ 3).

The amount of solid pulp for the goaf area was roughly calculated from (1) as

By substituting above values into (2), the amount of water, , was also calculated from (2) asThus, total amount of grouting slurry was calculated as

4.3. Grouting Result Simulation Verify

Figure 12 shows the numerical simulation result of velocity contour lines after filling with the grouting slurry by FLUENT. When the goaf is assumed to be filled uniformly with grouting slurry, the oxidation zone exists only in a very small zone just behind the #3205 working face. This result shows that the slurry grouting is a useful and effective method to prevent the coal spontaneous combustion in the goaf. The area of oxidation zone is affected by porosity in the goaf area. Figure 13 shows the numerical simulation result of the relationship between the area of oxidation zone and goaf porosity. The area increases exponentially with increasing goaf porosity. It shows that filling mined-out goaf area after coal cutting by grouting slurry makes reducing area of oxidation zone and less risk of spontaneous combustion of residual coal in goaf areas. It is also expected that reducing air leakage from the working face by control of ventilation pressure and airflow rate around there is an effective method to prevent spontaneous combustion.

Based on the numerical predictions presented in Figures 12 and 13, the Haizi Colliery applied the slurry grouting in the operation, and the accidents related to spontaneous combustion in the goafs have been not reported.

5. Conclusions

In this study, an effective and economical method to prevent spontaneous combustion of residual coal in a goaf area was investigated based on numerical predictions by ventilation simulator MIVENA and flow simulator FLUENT. The characteristics of oxidation and temperature rising in the goaf area in the Haizi Colliery have been studied by combining the tube bundle gas monitoring results and numerical simulation results. Furthermore, dividing the goaf area into three zones that are radiating, oxidized, and asphyxiation zones, were discussed in considerations of spontaneous combustion of residual coal. It has been cleared that filling porosity in the goaf area by injecting grouting slurry has a function to prevent the coal spontaneous combustion of residual coal in the goaf area.

The present research results are summarized as follows.(1)The mine ventilation simulator MIVENA was applied to analyze mine ventilation network in the mine, and the flow resistances at the #3205 working face was controlled to reduce air leakage into goaf area based on the simulation results.(2)The #3205 goaf area has been successfully divided into three zones in considering oxygen concentration and air leakage velocity by combining the tube bundle gas monitoring data and numerical simulation results by flow simulator FLUENT.(3)It has been cleared that slurry grouting is the most economical and effective method to prevent coal spontaneous combustion of residual coal in the goaf area by reducing goaf porosity and also area of oxidation zone behind the working face.

Conflict of Interests

None of the authors have any conflict of interests associated with this study, and financial supports listed in the acknowledgement are also not conflicted, because present research was carried out by international cooperation and financial supports.


This work was supported by JSPS KAKENHI (Grant-in-Aid for Scientific Research-B, Grant no. 25303030). This study was also done as the international cooperated research on coal mining safety, especially spontaneous combustion. Y. Wang thanks Japanese Government (MEXT) and Institute of Engineering and Environment LTU for supporting his study.