Abstract

In coal mines, a reasonable design of the widths of coal pillars is critically important, particularly for rockburst mines. This is because the frequent occurrence of rock burst in coal mines often arises from the inappropriate widths of the coal pillars. To address this problem, we first review two recent incidents of rockburst occurring and thus present a theoretical calculation of lateral bearing pressure distribution in goaf of inclined thick coal seam to illuminate the occurrence mechanism of rockburst induced by inappropriate width of sectional coal pillars. Based on the mechanism model, we then propose a design criterion for sectional coal pillar widths in inclined thick coal seams that can effectively reduce the risk of the induced rock burst. Our theoretical calculations, field stress monitoring, numerical simulation and field investigation all demonstrate the validity of the proposed design criterion in preventing the induced rock burst and large deformation of surrounding rock. The results from this paper may be used as a theoretical guidance of sectional coal pillar designs in inclined thick coal seams.

1. Introduction

Inappropriate width of coal pillars is an important factor that induces rockbursts in the working face along the goaf [13]. At present, in the central and eastern part of China, narrow coal pillars are used to protect roadways along the gob in mines with a depth of more than 600 m. However, In Xinjiang, Inner Mongolia, and Shanxi provinces, mines with depths between 300 m and 500 m still protect roadways with 15 m to 30 m wide sectional coal pillar, which are typically used in shallow mines. In the process of mining the working face, the mine pressure of roadway along the goaf is strong, and frequent rockburst and dynamic roof-fall accidents occur [4, 5], which causes many casualties and serious economic losses. On January 17, 2017, a rockburst accident occurred along the goaf of fully mechanized working face 4203 in Danshuigou Coal Mine of China National Coal Group, Shanxi Province; 9 people were killed and the 115 m long roadway was severely damaged (Figure 1). The average mining depth of working face 4203 was about 400 m, and the inappropriate design of a coal pillar of 20 m wide led to concentration of stress, which was the main cause of the rockburst accident.

If the design of the sectional coal pillar width along goaf working face in the rockburst mine well is inappropriate, rockbursts and large deformations in the surrounding rock may occur at the same time [69]. Many scholars have studied the proper width of coal pillars. Liu et al. [10] utilized field monitoring, theoretical calculation, and numerical simulation to investigate the distribution characteristics of lateral bearing pressure of fully mechanized working face in goaf of deep mine thick coal seams and proposed that four factors needed to be taken into account when determining sectional coal width. Hauquin T et al. [11] analyzed how an explicit numerical modelling method may be used to calculate and locate the damping of this kinetic energy during pillar failure. Cao et al. [12] studied the microseismic multidimensional information for the identification of rock bursts and spatial–temporal pre-warning in a specific coalface which suffered high rock burst risk in a mining area near a large residual coal pillar. Li et al. [13] studied the characteristic strata movement and mechanism underlying fault–pillar induced rock bursts (FPIRBs) and proposed that the most important factors affecting FPIRB are the static stress in the pillar and the dynamic stress induced by fault slides. Das AJ et al. [14] developed a generalised analytical solutions to estimate the strength of the coal pillars which can be applied for both the inclined and flat coal pillars. Jaiswal A et al. [15] describes an approach to estimation of strain-softening constitutive behaviour of coal-mass through calibration of a numerical model with field cases. Qin et al. [16] analyzed the relationship between the width of coal pillar and its permeability under the premise of stability and proposed an equation for the calculation of proper coal pillar width for preventing spontaneous combustion of coal pillar. Zhao et al. [17] carried out research and showed that the stability of coal pillar was mainly affected by the top main roof and multiple mining activities. With the increase of coal pillar width, the stress concentration zone shifted from the solid coal laterally to the middle of coal pillar. Li et al. [18] studied the stress distribution rules of sectional coal pillars under various buried depth and dip angles. When the dip angle of coal seam increased, the extent of stress concentration was more apparent. These scholars studied the distribution rules of bearing pressure in working face and sectional coal pillar width, but most were based on numerical simulation and field monitoring.

This paper is based on the coal pillar in the sectional fully mechanized working face (4-5)04 in an inclined thick coal seam in Liuhuanggou coal mine of Xinjiang Autonomous Region, China. Based on the rock mass test [19, 20], through theoretical calculation [2123], field monitoring [24], and numerical simulation [25], we determined the appropriate width of coal pillars, and this work may provide reference information for determining the sectional coal width of mine wells like the engineering conditions in the Xinjiang area.

2. Conditions of setting wide sectional coal pillars in Liuhuanggou coal mine

2.1. General conditions of the working face

Working face (4-5)04 of Liuhuanggou coal mine of Xinjiang is a gob-side working face with an average strike length of 2,636 m and incline width of 180 m. The main mining seam, seam 4-5, has a thickness of 6.06 m to 7.52 m and average thickness of 6.15 m. The dip angle is 21–26° with an average value of 25°. Figure 2(a) shows the plane position map of working face (4-5)04, which sets a 40 m wide sectional coal pillar to isolate goal (4-5)02. Figure 2(b) shows the section plan. The rock structure of the working face is provided in Table 1. As per the appraisal conclusion of rockburst tendency in seam 4-5 of Liuhuanggou coal mine, seam 4-5 was found to have a weak rockburst tendency, as shown in Table 2.

2.2. Mining conditions of wide coal pillars in the working face

After the working face (4-5)04 entered the stage of gob-side mining, the mining pressure in the tailgate became very strong due to the high stress concentration of the 40 m wide coal pillar, and frequent coal bursts resulted in a larger mining roof subsidence (Figure 3(a)). The large-diameter boreholes on the side of the coal pillar tended to close (Figure 3(b)). For this reason, establishing a properly sized sectional coal pillar, dispersing the coal pillar stress, and the deformation of surrounding rock in gob-side entry could significantly improve mining safety of the working face in seam 4-5 in Liuhuanggou coal mines.

2.3. Study concept

Our observations of the patterns of pressure in the Liuhuanggou coal mine show that the recovery of coal resources, the stability of roadway surrounding rock, and the prevention of rockburst and secondary disaster (such as flood, fire, and marsh gas) in adjacent goaf should be taken into account when determining the proper sectional coal pillar width. Of these, the stability and rockburst risk of surrounding rock in gob-side entry are related to the stress condition, and the recovery of coal resource and the prevention of secondary disaster in adjacent goaf are related to the width and integrity of coal pillar. This paper first studies the characteristics of lateral bearing pressure distribution in goaf of a fully mechanized working face in inclined thick coal seams, to determine the range of the low stress zone; then, the integrity of coal body in low stress zone is divided to ensure that the gob-side entry is not affected by high lateral stress and a secondary disaster in the goaf.

3. Theoretical calculation of lateral bearing pressure distribution in goaf of inclined thick coal seam

3.1. Theoretical calculation model

According to the monitoring results from the ground subsidence in Liuhuanggou coal mine, the ground subsidence coefficient after mining of working face (4-5)04 is only 0.1. This meant that it was in the incomplete mining stage, and the overlying rock in goaf is still in an overhanging state. Based on the monitoring results from ground subsidence, this research established a lateral bearing pressure estimation model in goaf of fully mechanized working face in inclined thick coal seam, as shown in Figure 4.

When the rock is in incomplete mining stage after mining of working face, the rock rupture height is approximately half of the width of goaf [26]. Because the overlying strata of both sides of the goaf are in the state of overhanging, half of the weight of the overhanging rock will act as a load on one side of the coal body in the goaf. The angles β and θ between the horizontal direction and the lines that connect the goaf and the stratum end are the rock stratum angular displacements.

The later bearing pressure σ on one side of coal body in goaf consists of self-weight stress σq and stress increment Δσ, i.e.,where Δσ equals to the sum of the pressure transferred to one side of the coal body by the hanging parts of the key layers above the goaf, i.e., Δσ = Σσi; σi is the pressure transferred to one side of coal body from the hanging part of the ith layer; i = 1∼n, n is the number of key layers above goaf.

The weight of each key layer of suspended coal transferring to one side of the goaf is half its weight, and the stress increment transferred to one side of the goaf is approximately in isosceles triangle distribution as shown in Figure 4. Therefore, the stress increment transferred to one side of coal body in goaf from the ith key layer is as follows:

Here, σmaxi is the maximum bearing pressure on one side of coal body in goaf produced by the ith key layer; σmaxi = 2Qi/Si; Mi is the thickness of the ith key layer; and Si is the influential distance of transferring stress from the ith key layer.

The weight Qi of each key layer is the sum of its self-weight and the weight of its controlled weak rock layer. The weight transferred from the key layer to one side of coal body is half of its weight, i.e.,

Here, Qi represents the weight transferred from key layer to one side of coal body in goaf; i = 1 – n, n is the number of key layers above coal layer; γ represents the bulky density of rock; Li is the overhang length for the thickness center of the ith key layer locating in goaf, and Li = 2I + Hicotθ − Hitan (α + β − π/)2. Hi is the distance from the thickness center of the ith key layer to the coal floor, and Hi = I+Mi/2+ΣMj(j = 1…i − 1); 2I is the width of goaf; Mi and mi are the thickness of the ith key layer and its controlled rock layer, respectively; γ is the bulky density of rock.

Because the △BCG part of the overlying rock in the inclined coal seam cannot form the rock beam structure to transfer the load to the front of the coal wall, the weight of the rock beam structure should be borne by the uppermost key layer of the working face. Then the weight transferred onto the coal body from the nth key layer is as follows:

The area of ΔBCG part from the overly rock layer is as follows:

According to the trigonometric function relation in Figure 4, the range of influence Si of the key layer from the ith layer on the lateral coal seam in goaf is obtained by the following:

The stress increment from n critical layers is superposed and the stress increment Δσ is obtained.

The stress σq produced by self-weight is as follows:where H is mining depth.

3.2. Calculation results

Because the overlying rock of seam 4-5 was mostly hard sandstone, according to the statistical results from borehole 27-2 near the working face, from seam 4-5 to the surface of the rock roof, the thickness of sandstone formation accounted for 96% of overall rock thickness. Considering the larger hardness and thickness of the overlying rock, to simplify the calculation process, the rock layer above the fracture range in goaf was calculated as one rock formation, and its thickness M1 was 309 m. According to the practical conditions of working face (4-5)04, the dip angle of coal seam was selected as α = 25°, the rock stratum angular displacement under the working face was β = 85°, the rock stratum angular displacement above the working face was θ = 80°, the mining depth of the tailgate in working face was 440 m, the mining depth of headgate was 500 m, the inclined length of working face was 2I=180 m, and the bulky density of rock γ is 25 kN/m3. Substitute these parameters into Equation (1) through Equation (7) and calculate the lateral bearing pressure distribution curves for the lower side of goaf in working face (4-5)04 as shown in Figure 5.

The red dotted line in Figure 5 represents strong rockburst risk. The determination basis was that the vertical stress in coal should be 2.0 times greater than the uniaxial compressive strength of coal [27]. As shown, the range of influence of lateral bearing pressure in goaf was approximately 52 m (OD section); the peak value of lateral bearing pressure in goaf was about 26 m (OB section) away from the goaf; the location 12 m to 45 m (AC section) away from the goaf was the high rockburst risk zone, and mining should not be performed in this range; OA section was a low rockburst risk zone (relative low stress zone), and its width is about 12 m (arranging mining roadway within this range needs to consider the dynamic stress issues during mining).

4. Field measurement of lateral bearing pressure distribution in goaf of inclined thick coal seam

To understand the characteristics of lateral bearing pressure of coal body in goaf during mining, the real monitoring alert system of rockburst was adopted to monitor the characteristics of the variations in bearing pressure of the heading side coal body in headgate of working face (4-5)04. One group of stress measurement stations was arranged in the headgate heading side from 46 m of advanced working face (4-5)04. Four stress gauges were installed at 4 m, 7 m, 12 m, and 16 m from the roadway (Figure 1). The boreholes were parallel to dip angle of coal seam.

These stress measurement stations started their monitoring service on April 12, 2017 (from 46 m of advanced working face). The monitoring ended on May 16, 2017 (entering 85 m of goaf) after 35 days of continuous monitoring. Figure 6 shows the relative stress variation curves of coal body. The distance below date showed the distance between working face and stress measurement station. Positive values indicate that the stress measurement station was located ahead of the working face and negative values indicate the stress measurement station was located behind the working face and had entered the goaf.

It was found from Figure 6 that the stress gauge at the 4 m deep measurement point began to continuously drop from 7.3 MPa to the minimum 0.42 MPa after the working face entered 13.5 m in goaf (April 25). The main reason was that the roof rotary compression caused the yield failure of the shallow coal seam, and the slow rise of stress later was due to the recovery of the stress in shallow coal body caused by the continuous rotary compression of the roof in the goaf. After the stress values at the 7 m and 12 m deep measurement points stabilized at the initial installation period, both stress curves exhibited a steady increasing tendency. This meant that two measurement points were located within the rotary influential range of lateral rock failure in goaf; the stress curve at the 16 m deep measurement point kept slowly decreasing tendency after installation. This meant that the lateral rock fracture position was in 16 m. Therefore, the lateral fracture position of the main roof of working face (4-5)04 was located between 12 m and 16 m, i.e. average lateral low stress range in goaf was 14 m.

5. Numerical simulation of lateral bearing pressure distribution in goaf of inclined thick coal seam

A numerical calculation model was established based on the geological and mining conditions of the working face (4-5)04 in Liuhuanggou mine, and the model size was length × width × height = 571 × 500 × 415 m with 161,387 nodes and element cells of 152,550 as shown in Figure 7. The four sides and bottom were displacement boundary. The 90 m high rock weight between model top to ground was substituted by applying 2.25 MPa distribution loading. The model calculation utilized Mohr-Coulomb criterion. The physical and mechanical parameters of the rock were obtained from the geological report and borehole data of the mine.

Mining of working faces (4-5)02 and (4-5)04 was simulated to resemble the actual mining situation of Liuhuanggou mine. The vertical stress distribution cloud map and plastic zone after mining completion of two working faces were presented in Figure 8 (partial enlarged detail).

As shown in Figure 8, after mining completion of working faces (4-5)02 and (4-5)04, 40 m wide sectional coal pillar exhibited a higher stress concentration, and there was apparent elastic core, which was the main source of mining pressure in the track grooves in the working face (4-5)04 and the headgates in (4-5)06 and large deformation of surrounding rock. The low stress zone width of the down side in the goaf of working face (4-5)04 was approximately 10 m, and the results of numerical simulation were consistent with those of theoretical calculation.

6. Determination of proper width of coal pillar in fully mechanized working face of inclined thick coal seam

6.1. Determination of lateral low stress zone range

The range of influence and peak position of bearing pressure obtained by theoretical calculation, stress monitoring, and numerical simulation are summarized in Table 3. By summarizing stress monitoring, theoretical calculation, and numerical simulation results, the width of low stress zone of the lateral bearing pressure in the goaf of the upper working face was determined to be 12 m.

6.2. Determination of sectional coal pillar width

The determination of sectional coal pillar width must consider the stability of gob-side entry. For this reason, the gob-side entry should be arranged within 12 m of goaf sides. From the viewpoints of controlling flooding, fire, gas, and other secondary disasters, the width of sectional coal pillar must be greater than 3 m. With consideration of the effectiveness of anchoring bearing (anchoring rod has a length of 2.2 m in field), the width of coal pillar should not be less than 3 m. From the viewpoint of preventing rockbursts during mining of working face, the width of low stress protection zone for gob-side entry solid side should not be less than 4 m, and given a 4 m wide gob-side entry, the width of sectional coal pillar should not be less than 5 m. Through comprehensive analysis, a proper coal pillar width for roadway protection should be between 3 m and 5 m.

6.3. Construction validation

Working face (4-5)06 of Liuhuanggou mine was mined forward along gob-side entry. Considering the fact that there was no experience in forwarding mining along gob-side entry for this mine well and adjacent mine wells, and severe disasters of flooding, fire and gas in goaf of mine wells, the sectional coal pillar width at the initial period of forwarding mining in gob-side entry was enlarged to 7 m. During mining, there was intense mining pressure. There were head-on entry explosions and frequent breakage of the anchor cable. The displacement on the two sides of entry, roof and floor all exceeded 500 mm (Figure 9(a)). The main cause was that a larger sectional coal pillar width resulted in the gob-side entry entering the lateral high stress range of goaf. The mining construction party later accepted the advice to reduce the sectional coal pillar width to 4 m, and both mining pressure and surrounding rock deformation greatly decreased (Figure 9(b)).

Therefore, according to the field pressure behavior of roadway along the gob and the control requirement of surrounding rock deformation, it is appropriated to have the width of coal pillar around 4 m.

7. Conclusion

Based on pre-setting the sectional coal pillar in fully mechanized working face of an inclined coal seam of Liuhuanggou coal mine, this research studied the proper sectional coal pillar width in fully mechanized working face in inclined thick coal seam by means of theoretical analysis, stress monitoring, and numerical simulation and reached the following main conclusions:(1)Through establishing the theoretical calculation model of lateral bearing pressure in goaf of inclined thick coal seam, this research concluded that the influence range of lateral bearing pressure in the goaf of working face (4-5)04 was about 52 m; the peak value position of lateral bearing pressure in goaf was about 26 m away from the goaf; the location 12 m to 45 m away from the goaf was high rockburst risk zone, and mining construction cannot be arranged within this range.(2)By summarizing the stress monitoring, theoretical calculation, and numerical simulation results, we here concluded that the width of low stress zone for the lateral bearing pressure in upper working face was 12 m. With consideration of prevention of secondary disasters in goaf, roadway bearing, and rockburst prevention, the proper width of sectional coal pillar for working face (4-5)06 was finally found to be between 3 m to 5 m.(3)Field practice showed that the determined sectional coal pillar was appropriate. When the sectional coal pillar was too large, rockburst and large deformation of surrounding rock could occur simultaneously; so field construction should strictly control the sectional coal pillar width.

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 there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge financial support from the National Natural Science Foundation of China (Grant No. 51874133, 51674014, 51904017, 51634001) and Engineering Research Center of Geothermal Resources Development Technology and Equipment, Ministry of Education, Jilin University (No. 21013). The authors also thank Liuhuanggou Coal Mine for support and assistance during field verification and application.