Numerical Investigation of Gob-Side Entry Retaining through Precut Overhanging Hard Roof to Control Rockburst

Gob-side entry retaining through precut overhanging hard roof (GERPOHR) method is one of the commonly used methods for nonpillar mining. However, feasibility studies of controlling rockburst by this method are few. Rockburst occurs in hard thick strata with a higher probability, larger scale, and higher risk. To better understand the GERPOHR method is beneficial for rockburst mitigation. In this paper, the design of GERPOHR was first introduced. And the layout of the working face was optimized. 1en, based on the numerical simulation, the stress and displacement distribution characteristics were compared under the condition of conventional mining and GERPOHRmethod.1e research shows that the intervals of main roof weighting could be decreased through the precut overhanging hard roof method. And the peak value of abutment pressure decreased. Meanwhile, the energy accumulation and the stress fluctuation could be alleviated in roadway surrounding rock.


Introduction
Rockburst is one of the major disasters in underground mining, especially under hard thick strata.e key reason is elastic strain energy that accumulates in overhanging hard roof and coal pillar [1][2][3].Based on incomplete statistics in China, rockburst disaster had happened 107 times in 1998-2014 (more than 30 people died) in Yima coal mine.
e rockburst accident of Shandong Chaoyang Coal Mine led five people to death in November 17, 2012.Since March 8, 1992, Huafeng Coal Mine rockburst has occurred more than 100 times, which has caused many casualties.And rockburst accidents have also occurred many times in Xing'an Coal Mine and Junde Coal Mine.
e common characteristics of these accidents are the coal seam with hard and thick immediate roof.
e safe and reliable control measures are important for rockburst.us, it still needs further study.e main controlling strategies of rockburst are to change the stress state and the pattern of energy release.e common method includes mine layout and pillar design [7,8], deep blasting relieving pressure technology [9], water infusion [10], stress relief boreholes [4], rock support [11], hard rock roof fracturing method [12], etc.
Controlling the rockburst caused by overhanging hard roof, the hard rock roof fracturing and no-pillar mining are efficient and robust methods.Hu et al. [12] proposed that the directional hydraulic fracturing method (deep-hole directional fracturing of thick hard roof) was an effective measure to control rockburst, by means of number simulation, drilling bits, and microseismic system.e long-hole drilling and the blasting methods were used to precrack hard roof so that the panel was smoothly extracted without any further rockburst [13].ere has been a large amount of research on explaining the mechanical mechanisms of rock bursting induced by a coal pillar and analyzing the pressure relief method to control rockburst [14,15].Gob-side entry retaining technology by overlying strata cutting was used as an effective method to decrease the pressure of roadway [16].However, current research was insufficient on the control mechanism of rockburst by gobside entry retaining technology.erefore, this paper introduced the GERPOHR technology and the control mechanism of rockburst.And based on numerical simulation theory, the comparative analysis between the conventional and GERPOHR mining method was investigated.
e contrast indexes include overhanging hard roof movement pattern, stress distribution, and roadways roof stability.
is study can provide guidance for understanding the influences of precut overhanging hard roof technology to control rockburst.

Rockburst Mitigation Mechanism of Gob-Side Entry Retaining through Precut Overhanging Hard Roof Method
Nonpillar mining method of gob-side entry retaining by roof cutting is emerging and developing in China [17][18][19] However, there are few studies on controlling rockburst using this method.In this paper, based on the GERPOHR method, the layout of the working face was improved and the hard roof method with the precut roof was adopted on both sides of the working face.e GERPOHR method is shown in Figure 1.Firstly, the roadway was strengthened to keep the roadway stable in the roof cutting process.e reinforcement measures commonly used have advance support and add constant resistance anchorage cable [20].en, the presplitting blasting was conducted inside the roadway to fracture the overlying thick rock.
e overlying roof will fall under mine pressure.Eventually, a roadway with solid coal on the one side and gangue goaf on the other side was created.e process would be repeated, and finally coal pillar would be nonexistent in the mining area.
e propagation path of stress in hard and thick strata will be cut off, and the roadway will be in the low stress region.Based on this technique, the elastic energy greatly decreased from coal pillar and load transfer of the roof from goaf.And the roof boundary condition with four edges fixed turns into opposite boundary completely clamped support and others free boundary condition.erefore, the first weighting interval of hard roof will decrease, and the bending elastic energy of the roof will decrease.

Numerical Simulation Establishment
Based on the analysis result, the GERPOHR method can decrease roadway surrounding rock stress and the effect of overlying hard roof.Based on the numerical simulation method, this research focuses on comparing with the traditional technology about the roof weighting step, stress distribution characteristics, and roadway roof stability.e simulation was based on the same geological condition.

Engineering Geology Situation. Junde Coal Mine 9103
working face is located in the region of No. 17 mining area, and the inclination is from 27 °to 32 °. e average thickness of coal seam is 12 m.e average depth of coal mining is 500 m. e coal seam is overlain successively by siltstone (7 m), fine sandstone (41 m), grit stone (25 m), fine sandstone (40 m), siltstone (2 m), and coal seam 15(0.9 m) and underlain successively by tuff (4.9 m) and fine sandstone (25.8 m), as shown in Figure 2. Fully mechanized slice mining is adopted in Junde Coal Mine 9103 working face, with the mining height of 3.5 m. e designed face is 1000 m in length and 150 m in width.e roadway was reinforced by "double deck steel mesh + steel anchor + anchor," with 5.6 m × 3.5 m (20 m 2 ) section.Five accidents of rockburst have occurred in the 9103 working face.
Based on the GB/T 25217.2-2010about rockburst tendency classification, the comprehensive judgement of hard roof indicates the rock with strong rockburst tendency.e formula was shown for calculating the bending energy index of single roof.e result of roof rockburst tendency is shown in Table 1: where U WQ is single roof bending energy index (kJ); R t is tensile strength of rock (MPa); h is roof thickness (m); and E is elastic modulus (MPa).
Based on the in situ stress result of No. 17 coal seam in Junde Coal Mine, the maximum principal stress is about 23 MPa on horizontal principal stress.
e direction is northward, and it is nearly parallel to the mining roadway.
e direction of minimum horizontal stress is east, 13 MPa.e vertical principal stress is 16.3 MPa.

Numerical Simulation Model.
is paper adopts 3DEC dynamic module.
e 3DEC numerical simulation software is a discrete element algorithm program.e software can truly demonstrate the nonlinear deformation state of strata movement caused by coal seam excavation.Rock blocks interact with each other through structural planes.When the structural planes exceed the bearing capacity, the shear dislocation or tensile disengagement will produce in rock blocks.Dynamic analysis is required for mining problems involving seismic release of energy (i.e., rockburst) and flow of particles (angular or rounded) under gravity.e auto damping method was applied in this paper, which is proposed by Professor Cundall [21].
e damping increased when the energy ratio is less than 0.5, and the damping decreased when the energy ratio is greater than 0.5.
e energy ratio can be calculated as follows: where R is the energy ratio; E D is the energy is absorbed by damping; and _ E K is system energy gradient.e numerical model is a cuboid of dimensions 410 m × 200 m × 300 m with gravity acting along the negative z axis.
e inclination of working face is 27 °, as shown in Figure 3(a).Figure 3(b) is the numerical calculation model under conventional mining condition.Figure 3(c) is the GERPOHR model.Figure 3(d) shows the location of monitoring point for roadway roof stress and displacement.

Evolution Characteristic of Vertical Displacement in Key
Strata.In order to study the fracture characteristic of overhanging hard roof (key strata) under the different mining methods, the vertical displacement of key stratum was analyzed on the upper working face 150 m. e z-displacement curve of key strata is shown in Figure 4. Before the strata movement, the coordinates of all nodes on the middle line of the hard thick rock stratum were obtained.e numerical software uses the Lagrange solution model.erefore, the final coordinates of the above nodes can be obtained after the stratum movement.
e difference between node coordinates is the distance of strata movement.
e maximum vertical displacement of key strata is 3.3 m.It showed that the fracture interval of key strata is 150 m under the conventional mining condition.e GERPOHR method contributed to decrease the fracture interval, the first weighting step 110 m.
Figures 5(a)-5(b) show the movement of overlying strata and the distribution of plastic zones under the GERPOHR condition.It shows that the rock stratum is asymmetric deformation under the effect of dip angle.In the second working face, the roadway is in the reduced pressure zone, which can effectively alleviate energy accumulation.It plays an important role in controlling rockburst.Figures 5(c)-5(d) show the movement of overlying strata and the distribution of plastic zones under the conventional mining condition.In the second working face, the roadway is in the pressure unrelieved zone.

e Distribution of Surrounding Rock Stress during Two Mining Methods.
e vertical stress will redistribute in the effect of mining disturbance.It may cause stress concentration near roadway and working face.Stress concentration increases the risk of rockburst occurrence.e rock stratum below 0.5 m of coal seam was selected as data extraction plane.
e distribution of vertical abutment pressure was obtained under the condition of conventional mining and GERPOHR method, as shown in Figure 6.Figures 6(a)-6(b) show that the stress concentration zone is located ahead of the coalfaces, surrounding roadways and coal pillar.In the mining process of second working face, the roadway was affected by the stress concentration on the entity coal and the coal pillar.e GERPOHR method can eliminate the stress concentration of coal pillars and decrease advance support pressure of working face, as shown in Figures 6(c)-6(d).
In order to better analyze the influence of two mining methods on roadway stability, the positions of the measuring lines were selected at advanced abutment stress area in the X-direction.Figure 7(a) demonstrates that stress concentration occurs on both sides of the roadway during upper working face mining.e peak value of abutment pressures is 25.1 MPa and 22.8 MPa, respectively, in conventional mining and precut overhanging hard roof mining condition.Because of the GERPOHR method, the gob-side entry retaining roadway is in the plastic zone.
e surrounding rock of roadway would not accumulate more stress.When mining at the second working face, precut overhanging hard roof mining can reduce abutment pressure, as shown in Figures 7(b)-7(c).

Dynamic Evolution Process of Roadway Roof Stress.
Figure 8 illustrates the movement pattern of roof under two different mining conditions.Under conventional mining condition, the hard and thick roof will gather a lot Advances in Civil Engineering of elastic energy before failure.Based on the pattern of main roof collapse [22], the exposed area of main roof strata in the goaf also increases with the continuous advance of working face.Based on the distribution pattern of rectangular plate bending moment, the middle of the long edge of the plate reaches the ultimate strength of rock layer at first to form a fracture when the exposed area increases to a certain value.en, the failure position occurs in the middle of the short edge.And the cracks will form around the plate and destroy gradually.e release of bending elastic energy led to produce the stress wave, and the subsidence of the roof will be rapidly transmitted to the    Advances in Civil Engineering roadway.e stress will fluctuate violently in the roadway roof.It may cause rockburst.Compared with conventional mining method, precut overhanging hard roof mining also produces elastic energy release of hard and thick strata.However, due to the effect of presplitting cut, stress wave transmission will greatly decrease.Meanwhile, the lateral influence caused by the subsidence of hard and thick strata will be avoided.e numerical simulation model can nearly realistically simulate the response of the low and hard thick strata.e detection point was selected for szz analysis (about 1.5 m from roof and the distance from the initial mining position is 160 m).When the hard and thick strata break, the vertical stress of roadway roof fluctuates from 40 MPa to 25 MPa under conventional mining condition.However, the vertical stress of roadway roof fluctuates more smoothly in the GERPOHR method (Figure 9).

Conclusion
Based on the mechanism of rockburst under the condition of hard and thick strata, this paper proposed the idea of preventing rockburst by the GERPOHR method.en, the stress distribution characteristics of conventional mining and precut overhanging hard roof mining were studied by using 3DEC discrete element software.e following conclusions were obtained: (1) e study shows that the GERPOHR method can reduce stress concentration of coal pillar.e roadway was located in the less stress region.us, the possibility of rockburst was reduced.Advances in Civil Engineering interval of main roof was reduced.Under the conventional mining condition, the hard and thick strata break when advancing about 140-150 m.Under the roof cutting condition, the initial pressure interval of hard and thick rock strata is about 110-120 m.After the gob-side entry retaining was completed, the surrounding rock of roadway would be in the plastic zone.us, the accumulation of energy could be avoided in second working face mining.Due to the influence of rock dip angle, the distribution of plastic zone is asymmetric.
(3) e GERPOHR method can reduce the peak value of abutment pressure.Under the traditional mining condition, the peak value of abutment pressure is about 27 MPa.e peak value of abutment pressure is about 21 MPa when the precut overhanging hard roof method is adopted.(4) eoretical and numerical analysis results show that this method can restrain the stress perturbation.Under the conventional mining condition, the stress fluctuation range caused by hard and thick rock fracture is 25-40 MPa.Stress fluctuation is more smooth under the GERPOHR method.

Figure 2 :
Figure 2: Geological and mining condition of selected panel.

Figure 3 :
Figure 3: Numerical model, model dimensions, and monitoring points.(a) Numerical simulation model.(b) Conventional mining model.(c) Precut overhanging hard roof model.(d) Monitoring points for roadway roof.

( 2 )Figure 6 :
Figure 6: Distribution of vertical abutment pressure under the conventional mining and GERPOHR condition.(a) Conventional mining to excavate upper working face.(b) Conventional mining to excavate second working face.(c) Precut overhanging hard roof mining to excavate upper working face.(d) Precut overhanging hard roof mining to excavate second working face.

Figure 7 :Figure 8 :
Figure 7: Distribution of stress around roadway.(a) Excavate upper working face.(b) Excavate second working face by conventional mining.(c) Excavate second working face by precut overhanging hard roof mining.

Table 1 :
Determination of rockburst tendency.In order to absorb the reflection of stress wave, viscous boundary was used in the numerical model.ephysicaland mechanical properties of intact rock used in the model are given in Table2.e models of intact rock and joint rock were the Mohr-Coulomb model and the Coulomb slip model, respectively.