Abstract

Coal and gas burst is one of the significant and catastrophic hazards in underground longwall operations. To date, the protective coal seam mining has been recognized as the most effective mining method for minimizing or even avoiding the effect of the coal and gas burst. In this paper, numerical modelling and field test were carried out for the longwall operation in Qidong Coal Mine in order to investigate the induced stress and coal seam gas drainage performance in the protected coal seam after the complete extraction of the protective coal seam. It was found that four stress zones can be classified in the protected coal seam being the original stress zone, stress concentration zone, stress relief zone, and recompaction zone. In addition, the monitoring data of gas concentration and volume change in the field agree well with the numerical modelling results.

1. Introduction

China has one of the largest mining industries in the world and due to the increasing complex geological conditions and mining depth, the frequency of the coal and gas burst occurrence tends to increase. The number of incidents resulted by the coal and gas burst in China accounts for 40% globally. As the mining depth increases, the geological stress increases and so does the concentration of the coal seam gas. The existing mines that used to have low gas concentration are now challenged by the increasing gas concentration as the mining operation moves deeper. In addition, the new mines face the challenge by high gas concentration from the beginning of the operation. The incidents caused by the coal and gas burst have been occurring from time to time [13]. The protective coal seam mining has been recognized as one of the most efficient mining methods in order to minimize or even prevent the coal and gas burst hazards and hence it has been investigated extensively globally. Wang et al. [4] confirmed that the protective coal seam mining is the most economical efficient method to prevent the coal and gas burst through analytical and field testing research. Wang et al. [5] pointed out that the coal seam gas predrainage is the most efficient way to eliminate the risk associated with coal and gas burst. Chen et al. [68] carried out research to investigate the permeability of the coal seam during mining operation and hence revealed that such a permeability would increase significantly due to the effect of stress relief and fractures in the seam resulted by mining activities. Nima et al. [9] identified the creep behaviour of the coal seam playing a significant role in the permeability of coal seam and gas drainage. Jin et al. [10] investigated the effect of the protective coal seam mining on the gas management through numerical modelling, statistical analysis, and residual gas analysis. Yuan et al. [11] investigated the pillar stability and change of permeability in coal seam during mining operation and hence concluded that the distance between adjacent coal pillars is a significant factor in stress relief in the protective coal seam and permeability change of coal seam. Li [12] pointed out compared to the protective coal seam mining, the extraction of the thin coal seam would minimize the coal and gas burst hazards and it is costly and timely economical. Wang et al. [13] used the soft coal seam adjacent to the thick seam as the protective coal seam. He et al. [14] investigated the permeability change in the upper coal seam during the extraction of the lower protective coal seam and developed an analytical model to simulate such a change during mining operation. Yuan et al. [15] conducted a series of laboratory tests simulating the filed condition to investigate the role of the distance between adjacent coal seams in protective coal seam mining. Shi and Liu [16] concluded that the distance between the protective and protected coal seams could affect the stress relief and deformation of the protected coal seam significantly. The deformation would tend to decrease as the distance increases and hence the crack initiation and propagation would be limited. Fang et al. [17] investigated the effect of the dip angle on the stress relief and deformation of the coal seam during the extraction of the protective coal seam and they concluded that as the dip angle increases, the deformation and damage in the overlying coal seam tend to decrease. Qin et al. [18] believed the computational fluid dynamics (CFD) is the optimized method for coal seam gas drainage design. Cao et al. [19] used a numerical model to simulate the stress-strain relationship in the lower protective coal seam during the protective coal seam mining operation and hence mapped out the stress relief zone which was further validated by theoretical approaches. Xue et al. [20] investigated the stress relief and permeability increase in the protective coal seam during mining operation and hence developed physical models that capture subsidence zone and fracture distribution zone in the protective coal seam. Wang et al. [21] conducted parametric studies to analyse the effect of a range of parameters associated with longwall panel design on the permeability change in the coal seam. Zhang et al. [22] conducted numerical simulation to study the influence of the overlying coal seam’s stress and strain during the mining of the lower coal seam and then concluded that the stress distribution would be categorized into three zones being the original stress zone where no deformation was detected, stress concentration zone where compression occurs, and stress relief zone where tensile failure occurs. It should be noted that the large deformation occurs in the stress relief zone.

Overall, it can be seen that extensive research has been carried out to investigate the effect of the protective coal seam mining on prevention of the coal and gas burst in underground longwall mining. However, the research on the stress distribution in the coal seams overlying and underlying the protective coal seam is limited, either theoretical or numerical work. In this research, through the field testing and numerical modelling, the stress distribution in the protective and protected coal seams was analysed. It was found that the numerical model results agree well with field monitoring data.

2. Stress Distribution in Overlying Protective Coal Seam

2.1. Input Parameters

Qidong Coal Mine is located in Suzhou City, Anhui Province, China (see Figure 1). There are seven seams that are economically feasible to be mined including seams 61, 62, 63, 71, 81, 82, and 9. At present, seams 61, 71, 82, and 9 are being extracted with average seam thicknesses 32.13, 31.95, and 10.29 m. The mining direction is downwards where seam 61 is the first seam to be mined. All the input parameters were obtained from longwall panel 6135 in the protective coal seam 61. The average seam thickness is 1.93 m. The average dip angle is 13°. The depth of cover ranges from 520 to 598 m. The panel length is 1770 m and face width is 200 m. The protected coal seam 71 has an average sea thickness of 1.42 m. The gas concentration is as high as 1.69 m3/t and outburst pressure can be up to 1.9 MPa. It can be confirmed that seam 61 is less likely to have gas burst occurred than seam 71, and therefore, seam 61 is extracted first to protect seam 71 such that the elastic energy can be reduced and permeability can be increased in seam 71. Such a strategy can aid the further gas drainage in seam 71 and hence reduce its gas concentration. The distance between seams 61 and 71 ranges from 40 to 46 m with an average 44 m.

Based on the aforementioned parameters for longwall panel 6135, the numerical model can be established to simulate the field condition as seen in Figure 2. Based on such a model, the stress distribution, stress relief area, and deformation in the overlying and underlying coal seams can be simulated.

In the model, X axis is along the dip direction and Y axis is along the longwall panel length. The size of the model is X × Y × Z = 400 m × 400 m × 289.7 m. The longwall panel width and length in the model are 100 m and 200 m, respectively. Constrain the velocities in the x, y, and z directions of the bottom boundary nodes, and constrain the horizontal velocity on the boundaries on both sides of x. The panel retreat starts from where Y = 100 m, and the longwall face moves 5 m at each step during modelling, and therefore, in total 40 steps, it would be simulated. Fixed boundary was assigned for bottom and sides, whereas free boundary was for the top of the model. The vertical stress from overburden is mainly resulted by its weight and hence 10 MPa was assigned at the top boundary considering the overburden thickness is around 400 m.

In total, 326,400 meshes and 341,172 nodes were generated. It is noteworthy that the mesh around the roof and floor of seam 61 was finer for pursuing the simulation. All the input parameters in the numerical model were determined by measurement and monitoring on the site (see Table 1).

2.2. Stress Distribution in the Roof and Floor of Protective Coal Seam

Goaf was formed after the retreat of the protective coal seam started. The roof would collapse and floor would heave significantly. Such failures at the roof and floor would have effect on the stress redistribution, permeability, outburst pressure, and deformation at the roof and floor of the seam. Therefore, investigation in these characteristics of the seam is very important.

The distribution of the vertical stress at the roof and floor of the seam after panel retreat is shown in Figures 3 and 4.

It can be seen from Figures 3 and 4 that the stress concentration tends to increase as the longwall face moves forward, and eventually, it tends to stabilize after longwall retreating for 120 m. It should be noted that the vertical stress along the dip direction is higher than that along the retreat direction. The vertical stress distribution along the retreat direction is symmetric, whereas due to the dip angle, the stress distribution along the dip direction is asymmetric. In addition, the stress concentration at the two sides along the retreat direction and lower side along the dip direction tends to increase and the concentration area tends to expand during panel retreating.

As the face moves forward, the middle part at the floor tends to release the stress. It is subjected to a high tensile stress leading to heave behaviour on the floor. The extent of the stress relief at the floor tends to decrease as the depth of cover of the coal seam increases. The stress relief area tends to increase as the longwall face moves forward and eventually reaches a threshold value.

2.3. Displacement in the Protective Coal Seam

Figure 5 shows the vertical displacement at the roof and floor of protected coal seam 71 as the longwall face in the protective coal seam advances up to 200 m. It can be seen from Figure 5(a) that the deformation distribution at the roof and floor at each retreating level is symmetric. The maximum vertical displacement at the roof in the goaf is 38.8 mm, whereas the maximum displacement at the floor is 27.6 mm. It can be seen from Figure 5(b) that the maximum tensile displacement in the protected coal seam 71 is 11.16 mm with the corresponding strain 7.86%. During the protective coal seam mining operation, the roof above pillars tends to move downwards and at the vicinity of pillar, the displacement tends to move horizontally towards goaf area. On the contrary, the floor underneath pillars tends to downwards, and at the vicinity of the pillar on the floor, the rock on the floor tends to move horizontally towards goaf area due to the compression. In the goaf area, the underlying coal and rock tend to move upwards due to the compression and tension on the floor. This leads to the formation of the new cracks on the floor and increase in the openness of the existing cracks. As a result, the permeability of the coal seam increases and such a change would aid in the stress relief and gas drainage in the protected coal seam.

2.4. Division of Four Zones of the Trend of the Protected Coal Seam

Based on the analysis of the stress distribution in longwall panel 6135 of the protective coal seam during mining operation, the stress distribution in the underlying protected coal seam 71 can be classified into four zones including original stress zone, stress concentration zone, stress relief zone, and recompaction zone (see Figure 6).

The area that is located over 38 m ahead of the longwall face is in the original stress zone. In such a zone, mining activity has negligible effect on the stress and strain in the protected coal seam 71. The area that is located within 8 to 38 m ahead of the longwall face is in the stress concentration zone. In such a zone, the protected coal seam is subjected to high stress leading to the close of the fractures as well as reduction in permeability and gas flow. The area that is located from 8 m ahead of the longwall face to 70 behind the longwall face is in the stress relief zone. In such a zone, the protected coal seam is subjected to a lower stress leading to a large displacement as well as increase in the permeability and gas flow in the seam. As a result, the risk of coal and gas burst in this zone is low such that it would be more efficient to carry out the gas drainage work in this zone. The area that is located 70 m behind the longwall face is the recompaction zone. In such a zone, the fracture rocks that collapsed from the roof tend to compact and consolidate, leading to the stress recovering towards the original level. However, the stress after the recompaction is still less than the original level and it would take extensive time for the stress to come back to the original level.

In summary, during the retreat of the longwall panel 6135, the vertical displacement at the floor would follow three stages in sequence including significant increase and significant decrease followed by returning to the original level. This can be observed by compressive and tensile deformation at the floor of the roadway.

3. Field Monitoring

Four stress zones have been classified in the aforementioned numerical modelling work. In order to validate such a hypothesis, field monitoring was conducted in Qidong Coal Mine. Figure 7 shows the pattern of the borehole in protected coal seam 71 for gas monitoring.

The stress evolution in the protected coal seam 71 follows four stages including original stress stage, stress concentration stage, stress relief stage, and recompaction stage. As a result, the gas concentration during the gas drainage would also vary correspondingly. The examples of gas concentration and flow quantity data in a number of selected boreholes are shown in Figures 812.

A few remarks can be concluded based on Figures 8 and 9.(1)For the boreholes located over 20 m ahead of the longwall face of the longwall panel 6135 in the protective coal seam, the gas concentration monitored fluctuated significantly in the range between 5% and 85%, whereas the flow quantity remains low in the range of 0.01 to 0.063 min. It could be observed that the gas drainage performance in this area is not optimized and the protected coal seam 71 was subjected to original stress condition.(2)For the boreholes located between 8 to 20 m ahead of the longwall face of the longwall panel 6135 in the protective coal seam, the gas concentration monitored reduced significantly to 15% and the flow quantity remains low in the range of 0.03 to 0.053 min. It could be observed that the protected coal seam 71 was subjected to stress concentration condition.(3)For the boreholes located between 8 ahead of the longwall face and 55 m behind the longwall face in the protective coal seam, the gas concentration monitored increased significantly up to 100% and the flow quantity also increased significantly to the range of 0.2 to 0.33 min. It could be observed that the protected coal seam 71 was subjected to stress relief condition.(4)For the boreholes located over 55 m behind the longwall face in the protective coal seam, the gas concentration monitored was lower than 40% and the flow quantity reduced to the range of 0.02 to 0.063 min. It could be observed that the protected coal seam 71 was subjected to recompaction condition.

Three more sets of gas monitoring data shown in Figures 1012 were used to further validate the gas monitoring data seen in Figures 8 and 9. As a result, the area of each of the four stress distribution zones was further calibrated based on the field monitoring data and numerical modelling results.(1)For the boreholes located 20 to 25 m ahead of the longwall face in the protective coal seam 61, the protected coal seam 71 was subjected to original stress condition. Such an area is only 13 m different to that determined by numerical modelling. In this area, the gas concentration fluctuated significantly in the range of 1% to 85%, whereas the flow quantity remained low in the range of 0.01 to 0.05 m3/min. It could confirm that the gas drainage performance would not be optimized in this area.(2)For the boreholes located 5 to 20 m ahead of the longwall face in the protective coal seam 61, the protected coal seam 71 was subjected to stress concentration condition. Such an area is only 2 to 10 m different to that determined by numerical modelling. In this area, the gas concentration reduced significantly to 15% and the flow quantity remained low in the range of 0.03 to 0.06 m3/min. It could confirm that the gas drainage performance would not be optimized in this area.(3)For the boreholes located between 5 m ahead of the longwall face and 65 m behind the longwall face in the protective coal seam 61, the protected coal seam 71 was subjected to stress relief condition. Such an area is only 5 m different to that determined by numerical modelling. In this area, the gas concentration surged significantly up to 100% and the flow quantity also increased significantly up to the range of 0.1 to 0.3 m3/min. It could confirm that the gas drainage performance in this area is optimal.(4)For the boreholes located over 65 m behind the longwall face in the protective coal seam 61, the protected coal seam 71 was subjected to recompaction condition. Such an area is only 5 m different to that determined by numerical modelling. In this area, the gas concentration was lower than 30% and the flow quantity also reduced to the range of 0.01 to 0.05 m3/min.

4. Conclusions and Outlook

Using numerical simulation to get the upper protective layer under mining effect, it is found that under the action of upper protective layer mining, the strike upper web stress contour is symmetrical saddle-shaped distribution. The stress of the bottom plate is distributed in an “eight” shape with a small top and a big bottom, and the vertical stresses of the top and bottom plates are distributed asymmetrically.

After the protection layer is mined, the coal and rock under the goaf moved upward as a whole, and the maximum swelling deformation reached 4.29‰, which expanded the original fissures and formed new fissures, creating conditions for the pressure relief gas drainage of the protected layer.

Based on the stress distribution results of the numerical simulation, the range of the strike “four zones” of the protected coal seam 71 was divided, and field verification was carried out. The range of the strike “four zones” was revised according to the gas drainage effect of the protected layer.

As an economical and effective technical means to prevent coal and gas outburst, protective layer mining has been studied more and more in recent years, and research methods have become more and more abundant, which has promoted the rapid development of protective layer mining theory. However, there is currently little research on the combined mining of the upper protective layer and the lower protective layer. With the continuous maturity of protective layer mining technology and the development of theory, under the superimposed effect of the combined mining of the lower protective layer and the upper protective layer, the temporal and spatial changes of coal seams, surrounding rock stress, and interlayer deformation will gradually become the focus of future research.

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

Acknowledgments

This work was financially supported by the National Key Research and Development Program of China (nos. 2017YFC0804206 and 2017YFC0804209) and the Anhui Provincial Natural Science Foundation (no. 2008085ME147).