# Study on the Mechanism and Control of Rock Burst of Coal Pillar under Complex Conditions

Authors: Xingping Lai; Huicong Xu; Jingdao Fan; Zeyang Wang; Zhenguo Yan; Pengfei Shan; Jie Ren; Shuai Zhang; Yanbin Yang
Journal: Geofluids
 (2020)

---

## Abstract

In order to explore the mechanism of coal pillar rock burst in the overlying coal body area, taking W1123 working face of Kuangou Coal Mine as the engineering background, the full mining stage of W1123 is simulated by FLAC3D. It is found that the high stress concentration area has appeared on both sides of the coal pillar when W1123 does not start mining. With the advance of the working face, the high stress concentration area forms X-shaped overlap. There is an obvious difference in the stress state between the coal pillar under the solid coal and the coal pillar under the gob in W1123. The concrete manifestation is that the vertical stress of the coal pillar below the solid coal is greater than the vertical stress of the coal pillar below the gob. The position of the obvious increase of the stress of the coal pillar in the lower part of the solid coal is ahead of the advancing position of the working face, and the position of the obvious increase of the stress of the lower coal pillar in the gob lags behind the advancing position of the working face. At the same time, in order to accurately reflect the true stress environment of coal pillars, the author conducted a physical similarity simulation experiment in the laboratory to study the local mining process of the W1123 working face, and it is found that under the condition of extremely thick and hard roof, the roof will be formed in the gob, the mechanical model of roof hinged structurer is constructed and analyzed, and the results show that the horizontal thrust of roof structure increases with the increase of rotation angle. With the development of mining activities, the self-stable state of the high stress balance in the coal pillar is easily broken by the impact energy formed by the sudden collapse of the key strata. Therefore, the rock burst of coal pillar in the overlying coal body area is the result of both static load and dynamic load. In view of the actual situation of the Kuangou Coal Mine, the treatment measures of rock burst are put forward from the point of view of the coal body and rock mass.


---

## Body

## 1. Introduction

Rock burst is a common phenomenon of rock dynamic failure in the process of coal mining, which is characterized by a large deformation of rock mass, throwing out rock blocks, producing a loud sound, ejecting gas waves, and so on [1, 2]. Rock burst is studied through various methods such as rock mechanics experiments, numerical calculations, and theoretical analysis. It is generally believed that the cause of rock burst is the accumulation of large amounts of elastic energy in the rock mass, and when the accumulation of elastic energy reaches the load limit of the rock mass, it will instantly release elastic energy [3–9].In the field of rock burst generation mechanism. Lai et al. pointed out that the occurrence of steeply inclined rock burst is caused by the prying effect of rock pillar between coal seams after it enters the depth with mining [10, 11]. Pan et al. put forward that the intrinsic factor of rockburst start-up is the excessive accumulation of static load in the near-field system of mining rock; the possible areas of rockburst start-up are the maximum peak areas of stress in the region of limit equilibrium [12]. Gong et al. improved the prediction accuracy of rock burst by adaptive-grid method and passive seismic tomography [13]. Zhu et al. investigated the influence of overstress state on rockburst in deep underground coal mining and found that compared to triaxial loading, coal samples are more easy to deform and bring to failure under triaxial unloading [14]. For Dai et al., the rock physical and mechanical parameters of 36 near horizontal coal seams are collected and the prediction model of coal pillar width is constructed [15, 16].Scholars at home and abroad have studied rockburst from the relatively independent aspects of rockburst mechanism, unloading, and prevention technology. However, the difference of stress distribution below the overlying coal body is rarely taken into account.This paper considers that due to the sudden thinning of coal seam in the W1145 working face of Kuangou Coal Mine, it is impossible to use existing equipment for mining, so a large number of coal body are left, which has an impact on the mining of W1143 working face which is below gob W1145. Though the rock burst is caused by the dynamic and static binding load mechanisms, the numerical simulation method is used to study the internal stress distribution characteristics of a coal pillar in different environments of overlying coal body area and overlying gob area. The mechanical model of the roof is constructed and analyzed, and the results show that the horizontal thrust of the roof structure increases with the increase of the rotation angle. The engineering treatment of the coal and rock mass is carried out from two aspects of structure and energy. The treatment effect is evaluated microseismic monitoring; microseismic data show that blasting and large-diameter boreholes can effectively reduce the accumulated elastic properties in the rock, which can provide guarantee for the safety of the personnel and property.


## 2. General Engineering Situation

### 2.1. Geology of W1123 Working Face

The Kuangou Mine is under the jurisdiction of Shenxin Company. The administrative division is under the jurisdiction of Queergou Town, Hutubi County, Changji Prefecture, Xinjiang. The coal mine is located in Shawan County-Manas County-Hutubi County earthquake belt.The coal mine main mining coal seam no. B4-1 and coal seam no. B2 are gently inclined coal seams. The average thickness of coal seam no. B4-1 is 3 m, and the two working faces are arranged in W1143 and W1145, in which the total length of the W1143 face is 164 m and the total length of the W1145 face is 175 m. The average thickness of the B2 coal seam is 9.5 m, and the two working faces are W1121 and W1123, respectively.The total length of W1121 is 149 m, and the total length of W1123 is 192 m. The average vertical distance between coal seam no. B4-1 and coal seam no. B2 is about 427.9 m, the thickness of bedrock is 421 m, the thickness of loose soil is 15.9 m, and the buried depth is about 427.9 m. The 15.9 m thick coarse-grained sandstone at 21.9 m above the B4-l coal seam is the main key layer of the overlying rock, and the 14 m thick coarse-grained sandstone at 21.8 m above coal seam no. B2 is the subkey layer of the overlying rock, and the mine coal seam belongs to the group of deep-buried coal seams at close distance. Longwall fully mechanized top coal caving is adopted.Due to the sudden thinning of coal seam no. B4-1, mechanized coal mining cannot be continued, so the remaining coal seam has not been mined. The strike length of W1145 is about 480 m and that of W1123 is about 1000 m. There will be two stages during the mining period of the W1123 working face. At the beginning of W1123, the open cut is under the thinning area of coal seam no.4-1; after mining about 400 m, the working face will be under gob W1145.The structural diagram of the Kuangou Coal Mine is shown in Figure1.Figure 1
Structural diagram of the Kuangou Coal Mine. (a) Plane diagram of the mine. (b) Section diagram of the mine. (c) Schematic diagram of the close-range coal seam.
(a)(b)(c)


### 2.2. Characteristics of Microseismic Distribution

Microseismic monitoring technology, which is based on acoustic emission and seismology, has developed into a new type of high-tech monitoring technology. It is a geophysical technology to monitor the influence, effect, and underground state of production activities by observing and analyzing small seismic events in production activities. When the underground rock breaks and moves due to manmade or natural factors, a weak seismic wave propagates around. By arranging several groups of geophones in the space around the fracture zone and collecting microseismic data in real time, after data processing, using the principle of vibration location, the location of the rupture can be determined and displayed in a three-dimensional space.A rockburst event occurred 60 meters in front of the working face during the mining period of W1123 on March 7, 2018. During the mining of the working face, the coal mine uses the microseismic monitoring system to monitor and locate the microseismic events and energy of the working face. Through the analysis of microseismic records, the energy value is3.1×105 J, and the specific position is on the side of the coal pillar in the headgate (Figure 2). The different colored dots in the picture represent the energy magnitude and location of microseismic events.Figure 2
Microseismic event distribution. (a) Plane diagram. (b) Section diagram.
(a)(b)When the high-energy event occurs, it is accompanied by a loud noise, resulting in a net pocket in the shoulder fossa in the lower section of the leading working face 50-55 m along the headgate, and the secondary small energy event occurs frequently after the occurrence of the high-energy event.We analyze events and energy collected by the microseismic equipment in the working face on March 7 and before extracting them as shown in Figure3. Before the occurrence of rock burst, the total energy and number of microseismic events continue to fluctuate, but three days before the occurrence of shock events, the energy and number of microseismic events begin to increase significantly.Figure 3
Characteristics of microseismic events.


## 2.1. Geology of W1123 Working Face

The Kuangou Mine is under the jurisdiction of Shenxin Company. The administrative division is under the jurisdiction of Queergou Town, Hutubi County, Changji Prefecture, Xinjiang. The coal mine is located in Shawan County-Manas County-Hutubi County earthquake belt.The coal mine main mining coal seam no. B4-1 and coal seam no. B2 are gently inclined coal seams. The average thickness of coal seam no. B4-1 is 3 m, and the two working faces are arranged in W1143 and W1145, in which the total length of the W1143 face is 164 m and the total length of the W1145 face is 175 m. The average thickness of the B2 coal seam is 9.5 m, and the two working faces are W1121 and W1123, respectively.The total length of W1121 is 149 m, and the total length of W1123 is 192 m. The average vertical distance between coal seam no. B4-1 and coal seam no. B2 is about 427.9 m, the thickness of bedrock is 421 m, the thickness of loose soil is 15.9 m, and the buried depth is about 427.9 m. The 15.9 m thick coarse-grained sandstone at 21.9 m above the B4-l coal seam is the main key layer of the overlying rock, and the 14 m thick coarse-grained sandstone at 21.8 m above coal seam no. B2 is the subkey layer of the overlying rock, and the mine coal seam belongs to the group of deep-buried coal seams at close distance. Longwall fully mechanized top coal caving is adopted.Due to the sudden thinning of coal seam no. B4-1, mechanized coal mining cannot be continued, so the remaining coal seam has not been mined. The strike length of W1145 is about 480 m and that of W1123 is about 1000 m. There will be two stages during the mining period of the W1123 working face. At the beginning of W1123, the open cut is under the thinning area of coal seam no.4-1; after mining about 400 m, the working face will be under gob W1145.The structural diagram of the Kuangou Coal Mine is shown in Figure1.Figure 1
Structural diagram of the Kuangou Coal Mine. (a) Plane diagram of the mine. (b) Section diagram of the mine. (c) Schematic diagram of the close-range coal seam.
(a)(b)(c)


## 2.2. Characteristics of Microseismic Distribution

Microseismic monitoring technology, which is based on acoustic emission and seismology, has developed into a new type of high-tech monitoring technology. It is a geophysical technology to monitor the influence, effect, and underground state of production activities by observing and analyzing small seismic events in production activities. When the underground rock breaks and moves due to manmade or natural factors, a weak seismic wave propagates around. By arranging several groups of geophones in the space around the fracture zone and collecting microseismic data in real time, after data processing, using the principle of vibration location, the location of the rupture can be determined and displayed in a three-dimensional space.A rockburst event occurred 60 meters in front of the working face during the mining period of W1123 on March 7, 2018. During the mining of the working face, the coal mine uses the microseismic monitoring system to monitor and locate the microseismic events and energy of the working face. Through the analysis of microseismic records, the energy value is3.1×105 J, and the specific position is on the side of the coal pillar in the headgate (Figure 2). The different colored dots in the picture represent the energy magnitude and location of microseismic events.Figure 2
Microseismic event distribution. (a) Plane diagram. (b) Section diagram.
(a)(b)When the high-energy event occurs, it is accompanied by a loud noise, resulting in a net pocket in the shoulder fossa in the lower section of the leading working face 50-55 m along the headgate, and the secondary small energy event occurs frequently after the occurrence of the high-energy event.We analyze events and energy collected by the microseismic equipment in the working face on March 7 and before extracting them as shown in Figure3. Before the occurrence of rock burst, the total energy and number of microseismic events continue to fluctuate, but three days before the occurrence of shock events, the energy and number of microseismic events begin to increase significantly.Figure 3
Characteristics of microseismic events.


## 3. Evolution Law of Stress Field during the Mining Progress

### 3.1. Numerical Simulating Model

Through on-the-spot investigation and analysis, the numerical simulation model of W1123 is established according to the geological column chart and the thickness and inclination of rock strata in the Kuangou Coal Mine. The model size is420m×525m×326m (length×width×height), and each roadway size is 4.5m×3.5m (width×height) Figure 4 show the three-dimensional model established by FLAC3D numerical simulation software.Figure 4
Numerical model indication.The numerical simulation model is founded by Mohr-Coulomb elements, and the gravity is set at 9.8 m/s2. The horizontal displacement constraint is imposed on both sides of the model, the displacement in the horizontal direction is limited, and the vertical direction is not limited. According to the results of rock stress measurement, as shown in Table 1.After angle conversion, the vertical load of 6.59 MPa is applied on the top of the model, the horizontal load of 11.91 MPa is applied on the front and rear boundaries of the model, the horizontal load of 9.92 MPa is applied on the left and right boundary of the model, and the mechanical parameters of the coal-rock mass are presented in Table 2.Table 1
Results of rock stress measurement.
Maximum principal stressIntermediate principal stressMinimum principal stressValue (MPa)Direction (°)Dip angle (°)Value (°)Direction (°)Dip angle (°)Value (MPa)Direction (°)Dip angle (°)12.8183.1-9.97.543.6-77.06.894.58.3Table 2
Physical and mechanical parameters of the coal-rock mass.
No.Rock typeT (m)γ (kg/m3)G (GPa)K (GPa)C (MPa)φ (°)1Mudstone25973.837.424.3930.412Sandy mudstone13.025464.418.125.4230.413Mudstone10.025973.837.424.3930.414Sandy mudstone14.025464.418.125.4230.415Sandstone13.024678.514.1716.2231.746Sandy mudstone8.025464.418.125.4230.417Packsand7.0130414.0719.5721.3838.868Mudstone6.025973.837.424.3930.419Sandy mudstone8.025464.418.125.4230.4110Coarse sandstone16.0254110.9413.4521.6329.9811Sandy mudstone9.025464.418.125.4230.4112Mudstone8.025973.837.424.3930.4113Coarse sandstone5.0254110.9413.4521.6329.9814Coal seam B4-13.013042.886.253.8137.4915Mudstone8.025973.837.424.3930.4116Coarse sandstone14.0254110.9413.4521.6329.9817Coal seam B31.813031.122.424.530.1018Mudstone3.025973.837.424.3930.4119Packsand16.0261814.0719.5721.3828.8620Coal seam B29.516401.913.914.931.2621Mudstone4.025973.837.424.3930.4122Packsand261814.0719.5721.3828.86Note: the top and bottom thickness of the model have different thickness with the change of dip angle.T is the thickness; γ is the unit weight, G is the shear modulus; K is the bulk modulus; C is the cohesion; and φ is the internal friction.


### 3.2. Stress Distribution Characteristics of the Coal Pillar

Figure5(a) is the distribution of vertical stress after penetration of the headgate, tailgate, and open-off cut in W1123 which shows that the area of high stress concentration has been shown in the coal pillar when the mining of W1123 has not begun, and the stress in the range of W1123 below the gob is obviously less than that below the solid coal. However, there are obvious zoning characteristics in the high stress area inside the coal pillar: the stress concentration area is widely distributed along the side of the gob of the coal pillar and is distributed in the shape of a point on the side of the gob of W1123. Because of the existence of solid coal in the upper layer, the stress is not released, and with the advance of W1123 working face, the scope of the gob is gradually expanded, and the vertical stress is mainly borne by the gob gangue and coal pillar.Figure 5
Vertical stress distribution of coal pillar: (a) advancing 0 m and (b) advancing 80 m.
(a)(b)When the model is used to simulate mining, the high stress area inside the coal pillar is distributed at the edge of the coal pillar at first, and the high stress area gradually moves to the middle of the coal pillar as the mining work continues. When the model mining reaches 80 m, the high stress area on both sides of the coal pillar produces X-shaped superposition in the middle of the coal pillar (Figure5(b)). When the high stress area is superimposed, the three-dimensional spatial distribution of stress field in coal pillar is shown in Figure 6.Figure 6
Relationship between peak stress of coal pillar and mining distance.We can cognize the characteristics of the zoning inside the coal pillar more figurative: (1) the coal pillar is basically divided into two regions. The first area has a higher overall stress below the upper layer solid coal. The lower overall stress in the second area is below the gob of the upper layer. (2) The stress of the coal pillar increases gradually from outside to inside and finally reaches the peak value on the side of W1123 in the middle of the coal pillar.


### 3.3. Spatial Evolution Law of Stress Field in Coal Pillar

Numerical simulation of the peak stress and peak stress increment of the coal pillar in the process of mining is shown in Figure7. By comparing the stress peak value and stress peak increment of the coal pillar under different mining progress, it is found that (1) in the mining process of W1123, the peak stress of coal pillar increases from 28.1094 MPa to 55.1590 MPa; the growth rate is 96.2%. The process of stress peak growth can be divided into three stages: stable zone, surge zone, and slow increase zone. (2) The peak stress increment curve is stable at first in the mining process of W1123 but increases suddenly when the model is excavated to about 80 m, and the high stress areas on both sides of the coal pillar are superimposed. (3) The peak stress increment fluctuates violently in the range of 80 to 150 m, which reflects that the coal pillar is obviously affected by nonlinear loading at this stage, and the stress in the coal pillar changes significantly, which is easy to induce rock burst.Figure 7
Spatial distribution of stress field in coal pillar.


### 3.4. Variation of Stress in Different Positions of Coal Pillar

Seven monitoring points are arranged in the interior of the coal pillar, which are located at 10 m, 70 m, 130 m, 190 m, 250 m, 310 m, and 370 m, respectively, of which 210 m is directly below the interface between the upper layer solid coal and the gob. The stress data of these seven points in the whole mining process are extracted. The stress curves at different positions of the coal pillar are shown in Figure7.By comparing the stress curves of different positions of the coal pillar (Figure8), it is found that there is an obvious difference between the stress of coal pillar under solid coal and that under gob: (1) the position of the obvious increase of the stress of the coal pillar in the lower part of the solid coal is ahead of the advancing position of the working face, that is, when the working face is not mined to the current area, the stress of the coal pillar in this area has begun to increase. (2) The position of the obvious increase of the coal pillar stress in the lower part of the gob lags behind the advancing position of the working face, that is, the coal pillar stress in this area begins to increase obviously after the end of the current mining area of the working face.Figure 8
Relationship between different positions and stress variation in coal pillar.


## 3.1. Numerical Simulating Model

Through on-the-spot investigation and analysis, the numerical simulation model of W1123 is established according to the geological column chart and the thickness and inclination of rock strata in the Kuangou Coal Mine. The model size is420m×525m×326m (length×width×height), and each roadway size is 4.5m×3.5m (width×height) Figure 4 show the three-dimensional model established by FLAC3D numerical simulation software.Figure 4
Numerical model indication.The numerical simulation model is founded by Mohr-Coulomb elements, and the gravity is set at 9.8 m/s2. The horizontal displacement constraint is imposed on both sides of the model, the displacement in the horizontal direction is limited, and the vertical direction is not limited. According to the results of rock stress measurement, as shown in Table 1.After angle conversion, the vertical load of 6.59 MPa is applied on the top of the model, the horizontal load of 11.91 MPa is applied on the front and rear boundaries of the model, the horizontal load of 9.92 MPa is applied on the left and right boundary of the model, and the mechanical parameters of the coal-rock mass are presented in Table 2.Table 1
Results of rock stress measurement.
Maximum principal stressIntermediate principal stressMinimum principal stressValue (MPa)Direction (°)Dip angle (°)Value (°)Direction (°)Dip angle (°)Value (MPa)Direction (°)Dip angle (°)12.8183.1-9.97.543.6-77.06.894.58.3Table 2
Physical and mechanical parameters of the coal-rock mass.
No.Rock typeT (m)γ (kg/m3)G (GPa)K (GPa)C (MPa)φ (°)1Mudstone25973.837.424.3930.412Sandy mudstone13.025464.418.125.4230.413Mudstone10.025973.837.424.3930.414Sandy mudstone14.025464.418.125.4230.415Sandstone13.024678.514.1716.2231.746Sandy mudstone8.025464.418.125.4230.417Packsand7.0130414.0719.5721.3838.868Mudstone6.025973.837.424.3930.419Sandy mudstone8.025464.418.125.4230.4110Coarse sandstone16.0254110.9413.4521.6329.9811Sandy mudstone9.025464.418.125.4230.4112Mudstone8.025973.837.424.3930.4113Coarse sandstone5.0254110.9413.4521.6329.9814Coal seam B4-13.013042.886.253.8137.4915Mudstone8.025973.837.424.3930.4116Coarse sandstone14.0254110.9413.4521.6329.9817Coal seam B31.813031.122.424.530.1018Mudstone3.025973.837.424.3930.4119Packsand16.0261814.0719.5721.3828.8620Coal seam B29.516401.913.914.931.2621Mudstone4.025973.837.424.3930.4122Packsand261814.0719.5721.3828.86Note: the top and bottom thickness of the model have different thickness with the change of dip angle.T is the thickness; γ is the unit weight, G is the shear modulus; K is the bulk modulus; C is the cohesion; and φ is the internal friction.


## 3.2. Stress Distribution Characteristics of the Coal Pillar

Figure5(a) is the distribution of vertical stress after penetration of the headgate, tailgate, and open-off cut in W1123 which shows that the area of high stress concentration has been shown in the coal pillar when the mining of W1123 has not begun, and the stress in the range of W1123 below the gob is obviously less than that below the solid coal. However, there are obvious zoning characteristics in the high stress area inside the coal pillar: the stress concentration area is widely distributed along the side of the gob of the coal pillar and is distributed in the shape of a point on the side of the gob of W1123. Because of the existence of solid coal in the upper layer, the stress is not released, and with the advance of W1123 working face, the scope of the gob is gradually expanded, and the vertical stress is mainly borne by the gob gangue and coal pillar.Figure 5
Vertical stress distribution of coal pillar: (a) advancing 0 m and (b) advancing 80 m.
(a)(b)When the model is used to simulate mining, the high stress area inside the coal pillar is distributed at the edge of the coal pillar at first, and the high stress area gradually moves to the middle of the coal pillar as the mining work continues. When the model mining reaches 80 m, the high stress area on both sides of the coal pillar produces X-shaped superposition in the middle of the coal pillar (Figure5(b)). When the high stress area is superimposed, the three-dimensional spatial distribution of stress field in coal pillar is shown in Figure 6.Figure 6
Relationship between peak stress of coal pillar and mining distance.We can cognize the characteristics of the zoning inside the coal pillar more figurative: (1) the coal pillar is basically divided into two regions. The first area has a higher overall stress below the upper layer solid coal. The lower overall stress in the second area is below the gob of the upper layer. (2) The stress of the coal pillar increases gradually from outside to inside and finally reaches the peak value on the side of W1123 in the middle of the coal pillar.


## 3.3. Spatial Evolution Law of Stress Field in Coal Pillar

Numerical simulation of the peak stress and peak stress increment of the coal pillar in the process of mining is shown in Figure7. By comparing the stress peak value and stress peak increment of the coal pillar under different mining progress, it is found that (1) in the mining process of W1123, the peak stress of coal pillar increases from 28.1094 MPa to 55.1590 MPa; the growth rate is 96.2%. The process of stress peak growth can be divided into three stages: stable zone, surge zone, and slow increase zone. (2) The peak stress increment curve is stable at first in the mining process of W1123 but increases suddenly when the model is excavated to about 80 m, and the high stress areas on both sides of the coal pillar are superimposed. (3) The peak stress increment fluctuates violently in the range of 80 to 150 m, which reflects that the coal pillar is obviously affected by nonlinear loading at this stage, and the stress in the coal pillar changes significantly, which is easy to induce rock burst.Figure 7
Spatial distribution of stress field in coal pillar.


## 3.4. Variation of Stress in Different Positions of Coal Pillar

Seven monitoring points are arranged in the interior of the coal pillar, which are located at 10 m, 70 m, 130 m, 190 m, 250 m, 310 m, and 370 m, respectively, of which 210 m is directly below the interface between the upper layer solid coal and the gob. The stress data of these seven points in the whole mining process are extracted. The stress curves at different positions of the coal pillar are shown in Figure7.By comparing the stress curves of different positions of the coal pillar (Figure8), it is found that there is an obvious difference between the stress of coal pillar under solid coal and that under gob: (1) the position of the obvious increase of the stress of the coal pillar in the lower part of the solid coal is ahead of the advancing position of the working face, that is, when the working face is not mined to the current area, the stress of the coal pillar in this area has begun to increase. (2) The position of the obvious increase of the coal pillar stress in the lower part of the gob lags behind the advancing position of the working face, that is, the coal pillar stress in this area begins to increase obviously after the end of the current mining area of the working face.Figure 8
Relationship between different positions and stress variation in coal pillar.


## 4. Design of Similar Material Model Experiment

The related results show that the similar material simulation experiment can reproduce the structure in the mining process of the working face [17, 18]. In order to further explore the rock burst mechanism under repeated mining in gently inclined extra thick coal seam, the roof hinged structure affected by repeated mining is studied by means of similar material simulation experiment.


### 4.1. Determining the Model Frame and Similarity Coefficient

The model test frame needs to be strong enough and has a certain width to ensure the stability of the model. A large number of researchers have explored the relationship between the ratio of simulated materials and the strength of real rock strata. The experimental model will be completely based on the on-site geological conditions of the Kuangou Coal Mine and the mining situation of working face.In this experiment, the model will be constructed by using the experimental frame of3000mm∗20mm∗1500mm, and the W1123 and its overlying strata will be simulated by using sand as an aggregate and large white powder and gypsum powder as cementitious materials. The simulation of rock formation strength in the field is achieved by adjusting the ratio of different raw materials. Combining the unilateral boundary conditions, then the constant of geometric similarity is as follows:(1)αl=lmlp=1200.According to the properties of the chosen simulation material and the ratio of the material, the constant of bulk density similarity is as follows:(2)αr=rmrp=1.According to similarity principle and dimensional analysis, the constant of geometric similarityαl, the constant of stress similarity αr, the constant of stress similarity ασ, and the constant of time similarity αt, at a certain relationship should be met.(3)ασ=αrαl,αt=αlThen, the constant of stress similarity isασ=0.081, and the constant of time similarity is αt=0.071.


### 4.2. Ratio and Amount of Similar Material

Through the analysis of the experimental data of rock mechanics provided by the coal mine, the ratio of aggregate and cementitious materials is determined, and the model is made.The specific ratio of the model material is as shown in Table3, the ratio of aggregate : lime : large white powder gives the experimental value.Table 3
Experiment ratio.
No.LithologySimulated thickness (cm)Material proportionSand (kg)Gypsum (kg)Flour (kg)1Mudstone40 : 3 : 713.870.521.212Sandy mudstone620 : 1 : 413.60.341.363Mudstone540 : 3 : 713.870.521.214Sandy mudstone720 : 1 : 413.60.341.365Sandstone645 : 1 : 413.30.291.186Sandy mudstone420 : 1 : 413.60.341.367Packsand470 : 9 : 2113.760.591.388Mudstone340 : 3 : 713.870.521.219Sandy mudstone720 : 1 : 413.60.341.3610Coarse sandstone835 : 3 : 1213.340.381.5211Sandy mudstone4.520 : 1 : 413.60.341.3612Mudstone420 : 1 : 413.60.341.3613Coarse sandstone2.535 : 3 : 1213.340.381.5214Coal seam B4-11.520 : 20 : 1 : 1.55.0850.2541.27215Mudstone440 : 3 : 713.870.521.2116Coarse sandstone735 : 3 : 1213.340.381.5217Coal seam B30.920 : 20 : 1 : 1.53.390.1690.84818Mudstone1.540 : 3 : 720.810.781.8219Packsand870 : 9 : 2113.760.591.3820Coal seam B24.7520 : 20 : 1 : 1.54.280.221.0721Mudstone283713.870.521.2122Packsand70 : 9 : 215362354Note: the top and bottom thickness of the model have different thickness with the change of dip angle. The material proportion of rock is aggregate : gypsum : flour. The material proportion of coal seam is aggregate : gypsum : flour : pulverized coal.


### 4.3. Structural Analysis

The results of similar material simulation experiments show that because the roof of the B2 coal seam is composed of a layer of extremely thick and hard fine sandstone; therefore, in the gob, it will not be fully broken under the action of gravity as other roofs, but to form a certain degree of integrity of the block. Because the block near the coal pillar does not have enough space to collapse, a kind of roof hinged structure is formed which is shown in Figures9 and 10.Figure 9
Characteristics of overlying strata in W1121working face mining.Figure 10
Characteristics of overlying strata in W1123 working face mining.The roof hinged structure is shown in Figure11, where α=14°, β is the angle between the rock block and the horizontal direction, and γ is the rotation angle of the rock block. The geometric relationship of the rock block in the process of rotation is shown in Figure 12.Figure 11
Indication of roof hinged structure.Figure 12
Geometric relations when rock rotates.In addition, when the extrusion pressure is linearly distributed along the interface, the position of the resultant force action point isn=p/3 [19]. The thickness and weight of the rock mass are much smaller than the thickness and weight of the overlying strata;therefore, in the process of calculation, the weight of the rotating rock block itself is ignored and the length of the rotating rock block is much smaller than that of the overlying strata; so, the load above the rotating rock block is simplified to a uniformly distributed load, and the force analysis is shown in Figure 13. From the force analysis, the following is obtained.Figure 13
Force analysis and geometric analysis.For mechanical condition,(4)RD=RB,(5)TD+qLtanγ=TB,(6)RD+RB=qL,(7)TDb+RDl=qL22.For geometric condition,(8)b=L′sinγcos14°≈Lsinγcos14°,(9)l≈L′cosγcos14°≈Lcosγcos14°.According to Equations (4)–(9),(10)TD=qL2•cos14°‐cosγsinγ.For the derivation ofγ,(11)TD′=qL2•1−sinγ‐cos14°sin2γ.In the process of roof rotation0<γ<90°, so TD′>0 (0<γ<90°), TD is a monotone increasing function about γ in the process of rotation.When the excavation of W1121 is over, the roof hinged structure appears only on one side of the pillar (Figure14(a)). The weight of some of the overlying strata is borne by the unilateral roof hinged structure formed by the W1121 gob. Therefore, the roof is bound to produce unilateral stress concentration at this time. And the roof hinged structure is a transitional state, and with the passage of time, the gob is gradually compacted under the load of the overlying strata, and the free space formed at this time will make the roof hinged structure rotate along the hinge point. The process of rotation will inevitably squeeze the roof on the coal pillar, which aggravates the degree of unilateral stress concentration. In addition, the coal quality of the B2 coal seam is hard, and the roof of the B2 coal seam also has impact tendency. So, a large number of elastic properties can be stored in the coal pillar, which provides an energy storage basis for the occurrence of rock burst, and the possibility and risk of dynamic disasters are obviously increased.Figure 14
Roof hinged structure in different mining stages. (a) Structure in gob W1121. (b) Structure in gob W1123.
(a)(b)When the excavation of W1123 is over, the roof hinged structure appears on both sides of the coal pillar (Figure14(b)). At this time, the coal pillar is affected by the bilateral roof hinged structure, which will inevitably produce bilateral stress concentration, and the average inclination angle of the B2 coal seam is 14°. The influence of the W1123 side structure on the internal structure of the roof is greater than that of the W1121 side structure.


## 4.1. Determining the Model Frame and Similarity Coefficient

The model test frame needs to be strong enough and has a certain width to ensure the stability of the model. A large number of researchers have explored the relationship between the ratio of simulated materials and the strength of real rock strata. The experimental model will be completely based on the on-site geological conditions of the Kuangou Coal Mine and the mining situation of working face.In this experiment, the model will be constructed by using the experimental frame of3000mm∗20mm∗1500mm, and the W1123 and its overlying strata will be simulated by using sand as an aggregate and large white powder and gypsum powder as cementitious materials. The simulation of rock formation strength in the field is achieved by adjusting the ratio of different raw materials. Combining the unilateral boundary conditions, then the constant of geometric similarity is as follows:(1)αl=lmlp=1200.According to the properties of the chosen simulation material and the ratio of the material, the constant of bulk density similarity is as follows:(2)αr=rmrp=1.According to similarity principle and dimensional analysis, the constant of geometric similarityαl, the constant of stress similarity αr, the constant of stress similarity ασ, and the constant of time similarity αt, at a certain relationship should be met.(3)ασ=αrαl,αt=αlThen, the constant of stress similarity isασ=0.081, and the constant of time similarity is αt=0.071.


## 4.2. Ratio and Amount of Similar Material

Through the analysis of the experimental data of rock mechanics provided by the coal mine, the ratio of aggregate and cementitious materials is determined, and the model is made.The specific ratio of the model material is as shown in Table3, the ratio of aggregate : lime : large white powder gives the experimental value.Table 3
Experiment ratio.
No.LithologySimulated thickness (cm)Material proportionSand (kg)Gypsum (kg)Flour (kg)1Mudstone40 : 3 : 713.870.521.212Sandy mudstone620 : 1 : 413.60.341.363Mudstone540 : 3 : 713.870.521.214Sandy mudstone720 : 1 : 413.60.341.365Sandstone645 : 1 : 413.30.291.186Sandy mudstone420 : 1 : 413.60.341.367Packsand470 : 9 : 2113.760.591.388Mudstone340 : 3 : 713.870.521.219Sandy mudstone720 : 1 : 413.60.341.3610Coarse sandstone835 : 3 : 1213.340.381.5211Sandy mudstone4.520 : 1 : 413.60.341.3612Mudstone420 : 1 : 413.60.341.3613Coarse sandstone2.535 : 3 : 1213.340.381.5214Coal seam B4-11.520 : 20 : 1 : 1.55.0850.2541.27215Mudstone440 : 3 : 713.870.521.2116Coarse sandstone735 : 3 : 1213.340.381.5217Coal seam B30.920 : 20 : 1 : 1.53.390.1690.84818Mudstone1.540 : 3 : 720.810.781.8219Packsand870 : 9 : 2113.760.591.3820Coal seam B24.7520 : 20 : 1 : 1.54.280.221.0721Mudstone283713.870.521.2122Packsand70 : 9 : 215362354Note: the top and bottom thickness of the model have different thickness with the change of dip angle. The material proportion of rock is aggregate : gypsum : flour. The material proportion of coal seam is aggregate : gypsum : flour : pulverized coal.


## 4.3. Structural Analysis

The results of similar material simulation experiments show that because the roof of the B2 coal seam is composed of a layer of extremely thick and hard fine sandstone; therefore, in the gob, it will not be fully broken under the action of gravity as other roofs, but to form a certain degree of integrity of the block. Because the block near the coal pillar does not have enough space to collapse, a kind of roof hinged structure is formed which is shown in Figures9 and 10.Figure 9
Characteristics of overlying strata in W1121working face mining.Figure 10
Characteristics of overlying strata in W1123 working face mining.The roof hinged structure is shown in Figure11, where α=14°, β is the angle between the rock block and the horizontal direction, and γ is the rotation angle of the rock block. The geometric relationship of the rock block in the process of rotation is shown in Figure 12.Figure 11
Indication of roof hinged structure.Figure 12
Geometric relations when rock rotates.In addition, when the extrusion pressure is linearly distributed along the interface, the position of the resultant force action point isn=p/3 [19]. The thickness and weight of the rock mass are much smaller than the thickness and weight of the overlying strata;therefore, in the process of calculation, the weight of the rotating rock block itself is ignored and the length of the rotating rock block is much smaller than that of the overlying strata; so, the load above the rotating rock block is simplified to a uniformly distributed load, and the force analysis is shown in Figure 13. From the force analysis, the following is obtained.Figure 13
Force analysis and geometric analysis.For mechanical condition,(4)RD=RB,(5)TD+qLtanγ=TB,(6)RD+RB=qL,(7)TDb+RDl=qL22.For geometric condition,(8)b=L′sinγcos14°≈Lsinγcos14°,(9)l≈L′cosγcos14°≈Lcosγcos14°.According to Equations (4)–(9),(10)TD=qL2•cos14°‐cosγsinγ.For the derivation ofγ,(11)TD′=qL2•1−sinγ‐cos14°sin2γ.In the process of roof rotation0<γ<90°, so TD′>0 (0<γ<90°), TD is a monotone increasing function about γ in the process of rotation.When the excavation of W1121 is over, the roof hinged structure appears only on one side of the pillar (Figure14(a)). The weight of some of the overlying strata is borne by the unilateral roof hinged structure formed by the W1121 gob. Therefore, the roof is bound to produce unilateral stress concentration at this time. And the roof hinged structure is a transitional state, and with the passage of time, the gob is gradually compacted under the load of the overlying strata, and the free space formed at this time will make the roof hinged structure rotate along the hinge point. The process of rotation will inevitably squeeze the roof on the coal pillar, which aggravates the degree of unilateral stress concentration. In addition, the coal quality of the B2 coal seam is hard, and the roof of the B2 coal seam also has impact tendency. So, a large number of elastic properties can be stored in the coal pillar, which provides an energy storage basis for the occurrence of rock burst, and the possibility and risk of dynamic disasters are obviously increased.Figure 14
Roof hinged structure in different mining stages. (a) Structure in gob W1121. (b) Structure in gob W1123.
(a)(b)When the excavation of W1123 is over, the roof hinged structure appears on both sides of the coal pillar (Figure14(b)). At this time, the coal pillar is affected by the bilateral roof hinged structure, which will inevitably produce bilateral stress concentration, and the average inclination angle of the B2 coal seam is 14°. The influence of the W1123 side structure on the internal structure of the roof is greater than that of the W1121 side structure.


## 5. Analysis of Dynamic Load Source

Based on the above analysis, the internal elastic performance of the coal pillar has accumulated greatly under the influence of the overlying coal body and the roof hinged structure in the gob, and if it is affected by the impact load, the instantaneous storage energy of coal pillar and roof may exceed the limit of its storage and release most of the energy for a short time, resulting in rock burst.


### 5.1. Dynamic Load by Key Stratum

Through the above analysis, the instability of the key strata structure formed during the overlying rock failure caused by coal seam mining is the dynamic cause of rock burst. There is a key strata with the thickness of 14 m above the W1123 working face. When the key strata structure breaks, it will hit the collapsed strata under the action of gravity and produce an elastic wave at the same time. When the elastic wave passes through the coal pillar which has reached the critical state, the coal pillar will instantly release a lot of energy and reflect it in the form of a rock burst.


### 5.2. Dynamic Load by Periodic Pressure

The data of hydraulic support sensor in the W1123 working face, 31 days before and 8 days after the impact event, support pressure distribution as shown in Figure15.Figure 15
Support pressure distribution.It can be seen from the diagram that there is a small periodic weighting in the W1123 working face around February 17. The period ended on February 18 and lasted for three days. The period reappeared around March 5 and ended on March 9 for a total of five days. The scale of the second period is wider, and the strength is larger than that of the previous one. And the impact event happens to be in the time period of large-scale weighting in the working face, so the periodic weighting of the W1123 working face is also one of the causes of the impact event.


### 5.3. Dynamic Load by Earthquake

The Kuangou Coal Mine is located in the Shawan County-Manas County-Hutubi County earthquake belt, where there are frequent earthquakes; this area is an earthquake-prone area. Therefore, the elastic wave caused by the earthquake can apply a dynamic load to working face which may trigger the rockburst to occur, but the earthquake did not occur on the March 7, so the earthquake has nothing to do with the occurrence of March 7 coal pillar rock burst, but enough attention should be paid to it. Because of the high degree of unpredictability of the seismic activity, the prevention and control measures of rock burst must be included in the daily work.


## 5.1. Dynamic Load by Key Stratum

Through the above analysis, the instability of the key strata structure formed during the overlying rock failure caused by coal seam mining is the dynamic cause of rock burst. There is a key strata with the thickness of 14 m above the W1123 working face. When the key strata structure breaks, it will hit the collapsed strata under the action of gravity and produce an elastic wave at the same time. When the elastic wave passes through the coal pillar which has reached the critical state, the coal pillar will instantly release a lot of energy and reflect it in the form of a rock burst.


## 5.2. Dynamic Load by Periodic Pressure

The data of hydraulic support sensor in the W1123 working face, 31 days before and 8 days after the impact event, support pressure distribution as shown in Figure15.Figure 15
Support pressure distribution.It can be seen from the diagram that there is a small periodic weighting in the W1123 working face around February 17. The period ended on February 18 and lasted for three days. The period reappeared around March 5 and ended on March 9 for a total of five days. The scale of the second period is wider, and the strength is larger than that of the previous one. And the impact event happens to be in the time period of large-scale weighting in the working face, so the periodic weighting of the W1123 working face is also one of the causes of the impact event.


## 5.3. Dynamic Load by Earthquake

The Kuangou Coal Mine is located in the Shawan County-Manas County-Hutubi County earthquake belt, where there are frequent earthquakes; this area is an earthquake-prone area. Therefore, the elastic wave caused by the earthquake can apply a dynamic load to working face which may trigger the rockburst to occur, but the earthquake did not occur on the March 7, so the earthquake has nothing to do with the occurrence of March 7 coal pillar rock burst, but enough attention should be paid to it. Because of the high degree of unpredictability of the seismic activity, the prevention and control measures of rock burst must be included in the daily work.


## 6. Pressure Relief Measures of Coal Seam

### 6.1. Rock Mass Engineering

The roof of coal seam is one of the most important factors affecting the occurrence of rock burst.Roof blasting is to break the roof, reduce its strength, release the energy gathered due to pressure, and reduce the impact vibration to the coal seam and support. The rock mass engineering is divided into two parts, namely, the roof advance precracking and the end advance precracking of the end.


#### 6.1.1. Calculation of Blasting Parameters

The average mining height of the W1123 working face is 11 m, the roof caving thickness is set, and the rock swelling coefficient is 1.4. In order to ensure that the caving roof can fully fill the goaf, the following formula is valid [20]:(12)Hx•ξ=HC+Hx.After calculation,Hx=11/1.4−1=27.5m.According to the geological conditions of the Kuangou Coal Mine, the upper 33 m of the B2 coal seam is the B3 coal seam. According to the calculation of roof treatment height, in order to leave a certain surplus coefficient, it is tentatively determined that the roof rock treatment range of the W1123 working face is along the vertical height from the roof of the working face to the upper 27 m.Because there is a certain gap in the lower 27 m collapse area, the pressure caused by the natural collapse of the remaining 6 m rock will be effectively degraded and will not produce a strong impact risk.The roof is weakened by deep-hole blasting, and the diameter of the hole is generally 70 mm~100 mm. When the blasthole diameter is small, the charge is more difficult, and the blasting influence range is proportional to the blasthole blasting diameter. The blasting effect is poor when the blasthole diameter is small, but when the blasthole diameter is too large, it is difficult to seal the mud, which will affect the blasting effect. Therefore, the vertical diameter of the roof blasthole weakened by deep-hole blasting should not be too large or too small. According to the experience of other working faces during mining, the blasthole diameter is determined to be 94 mm.After the explosive explosion, the crushing zone, crack zone, and vibration zone are formed from the explosion source outward.By calculating the range of the crack zone produced by blasting, the reasonable hole spacing can be determined.Because the blasting is carried out without free surface, the range of the crack zone of the pressure relief blasting is calculated according to the explosion stress wave in the case of uncoupled charge.In the blasting of uncoupled charge, the peak value of radial stress acting on the hole wall, that is, the initial impact pressurePr can be expressed as [21, 22]:(13)Pr=18ρDe2dcdb6n,where ρe is the explosive density, taking 1.27×103 kg/m3; D is the explosion velocity, taking 4400 m/s; dc is the explosive diameter, taking 90 mm; db is the blast hole diameter, taking 94 mm; and n is the stress increase multiple caused by explosive gas colliding with the rock, taking 12. After calculation Pr is 29411 MPa.Crack radiusR can be expressed as [23, 24](14)R=ν/1−νPrσt1/arb,where ν is the Poisson ratio of rock, taking 0.2; σt is the tensile strength, taking 7.48 MPa; a is the attenuation index, taking 1.5; and rb is the Blasthole radius; taking 47 mm. After calculation, R is 9081 mm.According to the abovementioned mechanical parameters of the roof and the performance parameters of the no. 3 emulsion explosive, the calculated results show that the diameter of the crack zone is larger than that of 9000 mm.Therefore, under the condition of ensuring enough safety, the roof of the working face is precracked at the same time, so that the roof cannot have a stress concentration, so the blast hole spacing is set as 10000 mm.In cylindrical charge blasting, the explosion resistance of coal and rock mass increases with the increase of blasthole depth, and the blasting ability of explosive is related to the length of plugging hole and the explosion resistance of coal and rock mass.If the plugging length is too short and the explosion resistance is reduced, the throwing funnel will affect the blasting effect during the explosion, especially the stability and support of the trench will be destroyed.If the plugging length is greater than its critical length, the blasting capacity of coal and rock is less than the antiexplosion capacity, and the coal and rock mass of the plugging section cannot form cracks during explosion, which also affects the blasting effect.Therefore, a reasonable plugging length should not only ensure that the coal and rock mass in the plugging section is loose and presplit but also cannot produce a throwing funnel at the same time.According to the general situation, the plugging length of deep-hole blasting should be 25%-30% of the hole depth, but considering that the blasting point is close to the roof of the coal seam, it is easy to cause roof leakage in front of the shelf. According to experience, the blasting point should be controlled above 4 m above the roof of the coal seam. Due to the limited space, the rest of the schemes are also calculated according to this method; it will not be repeated later.


#### 6.1.2. Roof Precracking

Construction begins at a distance of 30 m from the working face, and every 10 m a group of blastholes are arranged perpendicular to the center line of the roadway. The layout of the blasthole is shown in Figure16. The parameters of blasting hole drilling and filling explosive is shown in Table 4.Figure 16
Roof precracking blasthole layout scheme. (a) Plane diagram. (b) Section diagram.
(a)(b)Table 4
Parameters of roof precracking blasthole drilling and filling explosive.
Hole numberBorehole length (m)Horizontal angle (°)Elevation angle (°)Borehole diameter (mm)Filling explosive length (m)Plugging length (m)Explosive filling weight (kg)139904294261313024990259433161653459012943015150447-9038943215160545-9055943015150634-90819422121107319068942110105839903794261313094890149432161601050-9046943020751141-906594261565Note: the horizontal angle is the angle between the drilling direction and the center line of the roadway. Elevation angle is the angle between the drilling direction and the horizontal line of the roadway.


#### 6.1.3. End Precracking

In order to prevent the extremely thick and hard roof from forming a hinged structure on both sides of the coal pillar, it is necessary to blasting the side roof of the coal pillar. The blasting hole at the end is 20 m away from the workingface. There are three boreholes in each group of the headgate, with a depth of 30 m and a step of 5 m, and two boreholes in each group of the tailgate, with a depth of 30 m and a step of 5 m. The layout of the blasthole is shown in Figure17. The parameters of blasting hole drilling and filling explosive is shown in Table 5.Figure 17
End precracking blasthole layout scheme. (a) Plane diagram. (b) Section diagram in the headgate. (c) Section diagram in the tailgate.
(a)(b)(c)Table 5
Parameters of end precracking blasthole drilling and filling explosive.
Hole numberBorehole length (m)Horizontal angle (°)Elevation angle (°)Borehole diameter (mm)Filling explosive length (m)Plugging length (m)Explosive filling weight (kg)130-90679420101002300909420101003309062942010100430-908094201050530908094201050Note: the horizontal angle is the angle between the drilling direction and the center line of the roadway. Elevation angle is the angle between the drilling direction and the horizontal line of the roadway.


### 6.2. Coal Body Engineering

#### 6.2.1. Coal Seam Deep-Hole Blasting

The technological roadway is not only used to deal with the roof but also the coal seam deep-hole blasting is used to presplit the top coal. The hole spacing is 5 m and the hole diameter is 94 mm. The layout of the blasthole is shown in Figure18. The parameters of blasting hole drilling and filling explosive are shown in Table 6.Figure 18
Blasthole layout scheme. (a) Plane diagram. (b) Section diagram.
(a)(b)Table 6
parameters of deep-hole blasting, hole drilling, and filling explosive.
Hole numberBorehole length (m)Horizontal angle (°)Elevation angle (°)Borehole diameter (mm)Filling explosive length (m)Plugging length (m)Explosive filling weight (kg)148-9025943216160231903942011100Note: the horizontal angle is the angle between the drilling direction and the center line of the roadway. Elevation angle is the angle between the drilling direction and the horizontal line of the roadway.


#### 6.2.2. Large Diameter Drilling

The construction of large diameter pressure relief borehole in the coal body can release the elastic property accumulated in the coal body and eliminate stress accumulation. In the impact danger area, the stress value of the roof strata acting on the coal body reaches the critical stress state of the coal body. When the borehole is squeezed and deformed, the stress of the coal body is released and reduced to a new energy balance point in the process of deformation, thus reaching below the limit stress of the coal body impact and avoiding the occurrence of rock burst.The diameter of the drill bit is 153 mm, the drilling depth is 90 m, and the hole spacing is 8 m. Pressure relief boreholes are constructed along the coal seam.


### 6.3. Verification of Pressure Relief Measures

Take the data collected by microseismic monitoring equipment as an example. Analyze the number of events and energy collected by the microseismic equipment in the working face before and after the implementation of the treatment measures and extract them as shown in Figure19.Figure 19
Characteristics of microseismic events after pressure relief measures.On the day before the rock burst, the total energy of microearthquakes has been greatly increased, indicating that the elastic performance of the working face is gathering and finally reaches the cumulative limit at March 7, which leads to the occurrence of impact events. After the above blasting pressure relief measures, the total energy of microseismic events in the working face decreased greatly on March 8, and was lower than that produced in other normal mining periods. Therefore, it can be judged that the degree of elastic agglomeration of the blasting W1123 working face has been greatly reduced, and the pressure relief scheme is effective.


## 6.1. Rock Mass Engineering

The roof of coal seam is one of the most important factors affecting the occurrence of rock burst.Roof blasting is to break the roof, reduce its strength, release the energy gathered due to pressure, and reduce the impact vibration to the coal seam and support. The rock mass engineering is divided into two parts, namely, the roof advance precracking and the end advance precracking of the end.


### 6.1.1. Calculation of Blasting Parameters

The average mining height of the W1123 working face is 11 m, the roof caving thickness is set, and the rock swelling coefficient is 1.4. In order to ensure that the caving roof can fully fill the goaf, the following formula is valid [20]:(12)Hx•ξ=HC+Hx.After calculation,Hx=11/1.4−1=27.5m.According to the geological conditions of the Kuangou Coal Mine, the upper 33 m of the B2 coal seam is the B3 coal seam. According to the calculation of roof treatment height, in order to leave a certain surplus coefficient, it is tentatively determined that the roof rock treatment range of the W1123 working face is along the vertical height from the roof of the working face to the upper 27 m.Because there is a certain gap in the lower 27 m collapse area, the pressure caused by the natural collapse of the remaining 6 m rock will be effectively degraded and will not produce a strong impact risk.The roof is weakened by deep-hole blasting, and the diameter of the hole is generally 70 mm~100 mm. When the blasthole diameter is small, the charge is more difficult, and the blasting influence range is proportional to the blasthole blasting diameter. The blasting effect is poor when the blasthole diameter is small, but when the blasthole diameter is too large, it is difficult to seal the mud, which will affect the blasting effect. Therefore, the vertical diameter of the roof blasthole weakened by deep-hole blasting should not be too large or too small. According to the experience of other working faces during mining, the blasthole diameter is determined to be 94 mm.After the explosive explosion, the crushing zone, crack zone, and vibration zone are formed from the explosion source outward.By calculating the range of the crack zone produced by blasting, the reasonable hole spacing can be determined.Because the blasting is carried out without free surface, the range of the crack zone of the pressure relief blasting is calculated according to the explosion stress wave in the case of uncoupled charge.In the blasting of uncoupled charge, the peak value of radial stress acting on the hole wall, that is, the initial impact pressurePr can be expressed as [21, 22]:(13)Pr=18ρDe2dcdb6n,where ρe is the explosive density, taking 1.27×103 kg/m3; D is the explosion velocity, taking 4400 m/s; dc is the explosive diameter, taking 90 mm; db is the blast hole diameter, taking 94 mm; and n is the stress increase multiple caused by explosive gas colliding with the rock, taking 12. After calculation Pr is 29411 MPa.Crack radiusR can be expressed as [23, 24](14)R=ν/1−νPrσt1/arb,where ν is the Poisson ratio of rock, taking 0.2; σt is the tensile strength, taking 7.48 MPa; a is the attenuation index, taking 1.5; and rb is the Blasthole radius; taking 47 mm. After calculation, R is 9081 mm.According to the abovementioned mechanical parameters of the roof and the performance parameters of the no. 3 emulsion explosive, the calculated results show that the diameter of the crack zone is larger than that of 9000 mm.Therefore, under the condition of ensuring enough safety, the roof of the working face is precracked at the same time, so that the roof cannot have a stress concentration, so the blast hole spacing is set as 10000 mm.In cylindrical charge blasting, the explosion resistance of coal and rock mass increases with the increase of blasthole depth, and the blasting ability of explosive is related to the length of plugging hole and the explosion resistance of coal and rock mass.If the plugging length is too short and the explosion resistance is reduced, the throwing funnel will affect the blasting effect during the explosion, especially the stability and support of the trench will be destroyed.If the plugging length is greater than its critical length, the blasting capacity of coal and rock is less than the antiexplosion capacity, and the coal and rock mass of the plugging section cannot form cracks during explosion, which also affects the blasting effect.Therefore, a reasonable plugging length should not only ensure that the coal and rock mass in the plugging section is loose and presplit but also cannot produce a throwing funnel at the same time.According to the general situation, the plugging length of deep-hole blasting should be 25%-30% of the hole depth, but considering that the blasting point is close to the roof of the coal seam, it is easy to cause roof leakage in front of the shelf. According to experience, the blasting point should be controlled above 4 m above the roof of the coal seam. Due to the limited space, the rest of the schemes are also calculated according to this method; it will not be repeated later.


### 6.1.2. Roof Precracking

Construction begins at a distance of 30 m from the working face, and every 10 m a group of blastholes are arranged perpendicular to the center line of the roadway. The layout of the blasthole is shown in Figure16. The parameters of blasting hole drilling and filling explosive is shown in Table 4.Figure 16
Roof precracking blasthole layout scheme. (a) Plane diagram. (b) Section diagram.
(a)(b)Table 4
Parameters of roof precracking blasthole drilling and filling explosive.
Hole numberBorehole length (m)Horizontal angle (°)Elevation angle (°)Borehole diameter (mm)Filling explosive length (m)Plugging length (m)Explosive filling weight (kg)139904294261313024990259433161653459012943015150447-9038943215160545-9055943015150634-90819422121107319068942110105839903794261313094890149432161601050-9046943020751141-906594261565Note: the horizontal angle is the angle between the drilling direction and the center line of the roadway. Elevation angle is the angle between the drilling direction and the horizontal line of the roadway.


### 6.1.3. End Precracking

In order to prevent the extremely thick and hard roof from forming a hinged structure on both sides of the coal pillar, it is necessary to blasting the side roof of the coal pillar. The blasting hole at the end is 20 m away from the workingface. There are three boreholes in each group of the headgate, with a depth of 30 m and a step of 5 m, and two boreholes in each group of the tailgate, with a depth of 30 m and a step of 5 m. The layout of the blasthole is shown in Figure17. The parameters of blasting hole drilling and filling explosive is shown in Table 5.Figure 17
End precracking blasthole layout scheme. (a) Plane diagram. (b) Section diagram in the headgate. (c) Section diagram in the tailgate.
(a)(b)(c)Table 5
Parameters of end precracking blasthole drilling and filling explosive.
Hole numberBorehole length (m)Horizontal angle (°)Elevation angle (°)Borehole diameter (mm)Filling explosive length (m)Plugging length (m)Explosive filling weight (kg)130-90679420101002300909420101003309062942010100430-908094201050530908094201050Note: the horizontal angle is the angle between the drilling direction and the center line of the roadway. Elevation angle is the angle between the drilling direction and the horizontal line of the roadway.


## 6.1.1. Calculation of Blasting Parameters

The average mining height of the W1123 working face is 11 m, the roof caving thickness is set, and the rock swelling coefficient is 1.4. In order to ensure that the caving roof can fully fill the goaf, the following formula is valid [20]:(12)Hx•ξ=HC+Hx.After calculation,Hx=11/1.4−1=27.5m.According to the geological conditions of the Kuangou Coal Mine, the upper 33 m of the B2 coal seam is the B3 coal seam. According to the calculation of roof treatment height, in order to leave a certain surplus coefficient, it is tentatively determined that the roof rock treatment range of the W1123 working face is along the vertical height from the roof of the working face to the upper 27 m.Because there is a certain gap in the lower 27 m collapse area, the pressure caused by the natural collapse of the remaining 6 m rock will be effectively degraded and will not produce a strong impact risk.The roof is weakened by deep-hole blasting, and the diameter of the hole is generally 70 mm~100 mm. When the blasthole diameter is small, the charge is more difficult, and the blasting influence range is proportional to the blasthole blasting diameter. The blasting effect is poor when the blasthole diameter is small, but when the blasthole diameter is too large, it is difficult to seal the mud, which will affect the blasting effect. Therefore, the vertical diameter of the roof blasthole weakened by deep-hole blasting should not be too large or too small. According to the experience of other working faces during mining, the blasthole diameter is determined to be 94 mm.After the explosive explosion, the crushing zone, crack zone, and vibration zone are formed from the explosion source outward.By calculating the range of the crack zone produced by blasting, the reasonable hole spacing can be determined.Because the blasting is carried out without free surface, the range of the crack zone of the pressure relief blasting is calculated according to the explosion stress wave in the case of uncoupled charge.In the blasting of uncoupled charge, the peak value of radial stress acting on the hole wall, that is, the initial impact pressurePr can be expressed as [21, 22]:(13)Pr=18ρDe2dcdb6n,where ρe is the explosive density, taking 1.27×103 kg/m3; D is the explosion velocity, taking 4400 m/s; dc is the explosive diameter, taking 90 mm; db is the blast hole diameter, taking 94 mm; and n is the stress increase multiple caused by explosive gas colliding with the rock, taking 12. After calculation Pr is 29411 MPa.Crack radiusR can be expressed as [23, 24](14)R=ν/1−νPrσt1/arb,where ν is the Poisson ratio of rock, taking 0.2; σt is the tensile strength, taking 7.48 MPa; a is the attenuation index, taking 1.5; and rb is the Blasthole radius; taking 47 mm. After calculation, R is 9081 mm.According to the abovementioned mechanical parameters of the roof and the performance parameters of the no. 3 emulsion explosive, the calculated results show that the diameter of the crack zone is larger than that of 9000 mm.Therefore, under the condition of ensuring enough safety, the roof of the working face is precracked at the same time, so that the roof cannot have a stress concentration, so the blast hole spacing is set as 10000 mm.In cylindrical charge blasting, the explosion resistance of coal and rock mass increases with the increase of blasthole depth, and the blasting ability of explosive is related to the length of plugging hole and the explosion resistance of coal and rock mass.If the plugging length is too short and the explosion resistance is reduced, the throwing funnel will affect the blasting effect during the explosion, especially the stability and support of the trench will be destroyed.If the plugging length is greater than its critical length, the blasting capacity of coal and rock is less than the antiexplosion capacity, and the coal and rock mass of the plugging section cannot form cracks during explosion, which also affects the blasting effect.Therefore, a reasonable plugging length should not only ensure that the coal and rock mass in the plugging section is loose and presplit but also cannot produce a throwing funnel at the same time.According to the general situation, the plugging length of deep-hole blasting should be 25%-30% of the hole depth, but considering that the blasting point is close to the roof of the coal seam, it is easy to cause roof leakage in front of the shelf. According to experience, the blasting point should be controlled above 4 m above the roof of the coal seam. Due to the limited space, the rest of the schemes are also calculated according to this method; it will not be repeated later.


## 6.1.2. Roof Precracking

Construction begins at a distance of 30 m from the working face, and every 10 m a group of blastholes are arranged perpendicular to the center line of the roadway. The layout of the blasthole is shown in Figure16. The parameters of blasting hole drilling and filling explosive is shown in Table 4.Figure 16
Roof precracking blasthole layout scheme. (a) Plane diagram. (b) Section diagram.
(a)(b)Table 4
Parameters of roof precracking blasthole drilling and filling explosive.
Hole numberBorehole length (m)Horizontal angle (°)Elevation angle (°)Borehole diameter (mm)Filling explosive length (m)Plugging length (m)Explosive filling weight (kg)139904294261313024990259433161653459012943015150447-9038943215160545-9055943015150634-90819422121107319068942110105839903794261313094890149432161601050-9046943020751141-906594261565Note: the horizontal angle is the angle between the drilling direction and the center line of the roadway. Elevation angle is the angle between the drilling direction and the horizontal line of the roadway.


## 6.1.3. End Precracking

In order to prevent the extremely thick and hard roof from forming a hinged structure on both sides of the coal pillar, it is necessary to blasting the side roof of the coal pillar. The blasting hole at the end is 20 m away from the workingface. There are three boreholes in each group of the headgate, with a depth of 30 m and a step of 5 m, and two boreholes in each group of the tailgate, with a depth of 30 m and a step of 5 m. The layout of the blasthole is shown in Figure17. The parameters of blasting hole drilling and filling explosive is shown in Table 5.Figure 17
End precracking blasthole layout scheme. (a) Plane diagram. (b) Section diagram in the headgate. (c) Section diagram in the tailgate.
(a)(b)(c)Table 5
Parameters of end precracking blasthole drilling and filling explosive.
Hole numberBorehole length (m)Horizontal angle (°)Elevation angle (°)Borehole diameter (mm)Filling explosive length (m)Plugging length (m)Explosive filling weight (kg)130-90679420101002300909420101003309062942010100430-908094201050530908094201050Note: the horizontal angle is the angle between the drilling direction and the center line of the roadway. Elevation angle is the angle between the drilling direction and the horizontal line of the roadway.


## 6.2. Coal Body Engineering

### 6.2.1. Coal Seam Deep-Hole Blasting

The technological roadway is not only used to deal with the roof but also the coal seam deep-hole blasting is used to presplit the top coal. The hole spacing is 5 m and the hole diameter is 94 mm. The layout of the blasthole is shown in Figure18. The parameters of blasting hole drilling and filling explosive are shown in Table 6.Figure 18
Blasthole layout scheme. (a) Plane diagram. (b) Section diagram.
(a)(b)Table 6
parameters of deep-hole blasting, hole drilling, and filling explosive.
Hole numberBorehole length (m)Horizontal angle (°)Elevation angle (°)Borehole diameter (mm)Filling explosive length (m)Plugging length (m)Explosive filling weight (kg)148-9025943216160231903942011100Note: the horizontal angle is the angle between the drilling direction and the center line of the roadway. Elevation angle is the angle between the drilling direction and the horizontal line of the roadway.


### 6.2.2. Large Diameter Drilling

The construction of large diameter pressure relief borehole in the coal body can release the elastic property accumulated in the coal body and eliminate stress accumulation. In the impact danger area, the stress value of the roof strata acting on the coal body reaches the critical stress state of the coal body. When the borehole is squeezed and deformed, the stress of the coal body is released and reduced to a new energy balance point in the process of deformation, thus reaching below the limit stress of the coal body impact and avoiding the occurrence of rock burst.The diameter of the drill bit is 153 mm, the drilling depth is 90 m, and the hole spacing is 8 m. Pressure relief boreholes are constructed along the coal seam.


## 6.2.1. Coal Seam Deep-Hole Blasting

The technological roadway is not only used to deal with the roof but also the coal seam deep-hole blasting is used to presplit the top coal. The hole spacing is 5 m and the hole diameter is 94 mm. The layout of the blasthole is shown in Figure18. The parameters of blasting hole drilling and filling explosive are shown in Table 6.Figure 18
Blasthole layout scheme. (a) Plane diagram. (b) Section diagram.
(a)(b)Table 6
parameters of deep-hole blasting, hole drilling, and filling explosive.
Hole numberBorehole length (m)Horizontal angle (°)Elevation angle (°)Borehole diameter (mm)Filling explosive length (m)Plugging length (m)Explosive filling weight (kg)148-9025943216160231903942011100Note: the horizontal angle is the angle between the drilling direction and the center line of the roadway. Elevation angle is the angle between the drilling direction and the horizontal line of the roadway.


## 6.2.2. Large Diameter Drilling

The construction of large diameter pressure relief borehole in the coal body can release the elastic property accumulated in the coal body and eliminate stress accumulation. In the impact danger area, the stress value of the roof strata acting on the coal body reaches the critical stress state of the coal body. When the borehole is squeezed and deformed, the stress of the coal body is released and reduced to a new energy balance point in the process of deformation, thus reaching below the limit stress of the coal body impact and avoiding the occurrence of rock burst.The diameter of the drill bit is 153 mm, the drilling depth is 90 m, and the hole spacing is 8 m. Pressure relief boreholes are constructed along the coal seam.


## 6.3. Verification of Pressure Relief Measures

Take the data collected by microseismic monitoring equipment as an example. Analyze the number of events and energy collected by the microseismic equipment in the working face before and after the implementation of the treatment measures and extract them as shown in Figure19.Figure 19
Characteristics of microseismic events after pressure relief measures.On the day before the rock burst, the total energy of microearthquakes has been greatly increased, indicating that the elastic performance of the working face is gathering and finally reaches the cumulative limit at March 7, which leads to the occurrence of impact events. After the above blasting pressure relief measures, the total energy of microseismic events in the working face decreased greatly on March 8, and was lower than that produced in other normal mining periods. Therefore, it can be judged that the degree of elastic agglomeration of the blasting W1123 working face has been greatly reduced, and the pressure relief scheme is effective.


## 7. Rockburst Control Method and Strategy

According to the analysis above, two aspects affect the internal stress of the coal pillar. On the one hand, the coal pillar carries the static load of the strata above; on the other hand, it is the dynamic load caused by mining in the working face. The dynamic load caused by the sudden collapse of the key strata after reaching the limit span with the increase of the mining space is the main reason for the rock burst. The coal pillar can still maintain its own steady state under the influence of mining, but it is close to the equilibrium limit. At this time, if the energy produced by the structural instability of the key strata will cause the equilibrium limit failure, it is very likely to cause rock burst. Therefore, the prevention and control methods of rock burst in the W1123 working face should be considered from the following two aspects, as shown in Figure20.(i)
The rock mass engineering: the advance precracking of the roof will weaken the impact load caused by the breaking of the key strata, because the roof will play a bearing role on the overlying strata, and the treatment of the roof in time will dissipate actively the elastic energy accumulated above the roof under human intervention. Avoid a sudden release of a large amount of elastic power at one time. Precracking the end is actually a method of cutting the roof, that is, allowing the roof to collapse along the designed position to avoid forming a roof hinged structure to form a high stress area inside the coal pillar(ii)
The coal body engineering: blasting of the coal body can not only improve the resource recovery rate of top coal caving but also cause many cracks in the coal body, and large diameter drilling can cause free space in the coal body. To guide the stress release in the coal body, the above two measures are to actively release the elastic energy in the coal body to avoid the accumulation of elastic energyFigure 20
Prevention and control measures of rock burst.In the effect detection method, the microseismic monitoring system, acoustic emission monitoring system, and electromagnetic emission monitoring system can process the vibration signal, acoustic signal, and electromagnetic radiation signal to reflect the degree of deformation, fracture, and energy release in the coal and rock mass [25–30]. The drilling television system can directly reflect the effect of roof cutting and precracking after blasting [31, 32]. Through these kinds of methods, the effect of treatment is evaluated, and the parameters of treatment measures in the later stage are provided.


## 8. Conclusion

(1)
The numerical simulation shows that the high stress concentration area has appeared on both sides of the coal pillar when the W1123 working face does not start mining. With the advance of the working face, the high stress concentration area forms an X-shaped overlap. The position of the obvious increase of the stress of the coal pillar in the lower part of the solid coal is ahead of the advancing position of the working face, and the position of the obvious increase of the stress of the lower coal pillar in the gob lags behind the advancing position of the working face(2)
The result of a similar material simulation experiment shows that it is easy to form roof hinged structure under the condition of super thick and hard roof, and the existence of structure greatly increases the stress value in the interior of coal pillar(3)
The pillar rock burst in the W1123 working face is the result of both static load and dynamic load, but the fuse causing the rock burst is the energy release caused by the breaking of the key strata. The energy released by the key strata breaks the stress balance limit in the coal pillar and results in the occurrence of rock burst(4)
Considering the various factors of induced rock burst, the source of induced rock burst will be dealt with the strata to avoid the accumulation of elastic properties of roof strata and deal with coal seams to avoid the accumulation of elastic properties of coal seams. The effect of treatment measures is evaluated, and the parameters of treatment measures are adjusted by using microseismic monitoring, acoustic emission monitoring, and other equipment to form a feedback regulation mode



---

*Source: 8847003-2020-10-28.xml*