Influence of Key Strata on the Evolution Law of Mining-induced Stress in the Working Face of Deep and Large-scale Mining

1. Introduction

The movement of overlying strata after coal mining will exert some impact on the stress evolution of underground working face. The key stratum plays a major role in bearing the overlying strata [1]. Accordingly, it is of particular importance to master the influence of the key stratum distribution in the overlying strata on the evolution of mining-induced stress in the working face.

Mining-induced stress distribution of the working face in a kilometer-deep coal mine has been studied [2]. Different key stratum thicknesses and heights have been numerically simulated by the Universal Distinct Element Code (UDEC) software in order to study the effect of key stratum on the mining abutment pressure of a coal seam [3]. The distribution of energy accumulation and fracture positions before and after the fracture of overlying key strata are derived, and the energy release of fractures in each stratum is calculated [4]. The mechanism and evolution control of pillar bursts in multithick key strata are studied using field investigation, theoretical analysis, and numerical simulation [5]. The spatial fracture characteristic of overlying strata was analyzed by Winkler elastic foundation beam theory [6]. Based on the Key strata theory an overburden caving model is proposed to predict the multilayered hard strata behaviour [7]. A new damage model based on the modified thermomechanical continuum constitutive model in coal mass and the contact layers between the rock and coal mass is proposed [8]. Fracture failure analysis of hard and thick key layer and its dynamic response characteristics were analyzed [9]. Strata movement and stress evolution when mining two overlapping panels affected by hard stratum has been studied [10]. The analysis of the distribution characteristics of the primary key stratum in the coal mine reveals the bow-shaped structural characteristics of the overlying thick primary key stratum [11]. The mechanism, prevention measures, and control methods for earthquake disasters typically occurring in mines with thick and hard rock strata were investigated [12]. A simplified mechanical model for the analysis of dynamic destabilization of the overlying strata during underground mining was constructed [13]. A study investigated the spatiotemporal effect of mining-induced stress-fracture-seepage field coupling, characterized by ladder key stratum and composite aquifer [14]. Wilson’s equations for the vertical stress distribution in the vicinity of a single longwall panel after it has been mined have been used in conjunction with finite element modelling to evaluate vertical stresses in the underlying strata [15]. An estimation method for cover pressure re-establishment distance and pressure distribution in the goaf of longwall coal mines [16]. Gravitational rules are used as the basis for the prognosis of vertical stress distribution in exploited rock masses [17]. To analyze the law of movement and caving of the roof rock stratum, the roof subsidence displacement, rock stratum stress, and the rock stratum movement law were analyzed by using the methods of the particle discrete element and similar material simulation test [18]. Based on the analysis of subsidence and overlying structural characteristics, the influencing factors and the control effect laws of strata movement were investigated [19].

A working face in Binchang mining area of China is buried nearly a kilometer deep, with a mining height of over 10m. After coal mining, the surface subsidence is generally slight, and the underground rock burst accidents occur frequently, indicating that the movement law of the overlying strata is of certain particularity. Through field monitoring and numerical simulation, this study explored the strata movement law of the mine after mining and analyzed the influence of the key stratum of overlying strata on the evolution of mining-induced stress in the working face, laying the foundation for revealing the mechanism of rock burst.

2. Project overview

The plane layout of the working face in the mine under study is shown in Figure 1. First, the four working faces LW101~ LW104 in panel 1 were mined; then LW201~ LW 205 in panel 2 were mine successively; then LW302 in panel 3 were mined; and LW 301 in panel 3 was mined finally. This study mainly analyzed the evolution law of mining-induced stress after mining LW204 and LW205, and then LW302 and LW301.

3. Field monitoring

In order to master the internal movement law of overlying strata, an observation borehole ZY1 was arranged on the ground surface of LW205 with a depth of 940m (the layout position is shown in Figure 1). A dual-channel DOFS with a depth of 930m was arranged inside the borehole (the scheme is shown in Figure 2a), and the on-site installation process is shown in Figure 2b. Meanwhile, a set of GNSS system was also arranged on the surface to monitor the changes in surface subsidence in real time (Figure 2c).

After the drilling of observation borehole ZY1 was completed, the corresponding columnar distribution results of strata were obtained by logging in the hole. By means of the key stratum discrimination software, the ZY1 column was judged and the key stratum distribution was obtained, as shown in Figure 3. It can be seen from Figure 3 there are 14 key strata in the whole overburden above the coal seam.

With ZY1 being 292m ahead of the working face, the installation of DOFS in the borehole and GNSS monitoring system for surface subsidence was completed on November 27, 2020 and strain measuring was conducted for the first time. From then on, strain measuring was conducted regularly until October 26, 2021 when the mining in the working face was over. At that time, the working face had advanced beyond ZY1 by 380.4m. Even after the mining in the working face was over, strain measuring was still conducted for many times; and the last measuring was conducted on December 6, 2021. The strain value measured at the first time was set as the benchmark; then the strain data obtained from each subsequent measuring was conducted overall differential processing with the benchmark. The strain difference curves obtained are shown in Figure 4. As can be seen, the negative value in the vertical coordinate represents that the rock formation is under compressive strain, while the positive value indicates that the rock formation is under tensile strain.

As can be seen from Figure 4, when ZY1 was ahead of the working face (the corresponding mining time was before May 10,2021), the overlying strata were affected by the advance stress. On May 6, 2021, when ZY1 was 12.6m ahead of the working face, the DOFS at a depth of 900.5m in the borehole was broken because the compression was beyond its strain limit.

On May 10, 2021, when the working face advanced past ZY1 by 0.5m, the strain form of the DOFS in the borehole changed from compressive strain to tensile strain and the DOFS broke at the depth of 887.4m, which corresponded to KS1 in the borehole.

On May 17, 2021, when the working face advanced past ZY1 by 18.4m, the DOFS broke at the depth of 853.7m, which corresponded to KS2 in the borehole.

On May 27, 2021, when the working face advanced past ZY1 by 47.6m, the DOFS broke at the depth of 838.7m, which corresponded to KS3 in the borehole.

On June 1, 2021, when the working face advanced past ZY1 by 61.3m, the DOFS broke at the depth of 818.3m, which corresponded to KS4 in the borehole.

On June 5, 2021, when the working face advanced past ZY1 by 71.3m, the DOFS broke at the depth of 772.1m, which corresponded to KS5 in the borehole.

On June 9, 2021, when the working face advanced past ZY1 by 81.5m, the DOFS broke at the depth of 609.8m, which corresponded to the range between KS6 and KS7 in the borehole. Since then, as the working face continued to advance, the breakpoint height of theDOFS tended to be stable.

On October 26, 2021, when the mining on the working face was over, the working face had advanced past ZY1 by 380.4m. On December 6, 2021, the strain was measured for the last time, which revealed that the stain value increased slightly but the position of breakpoint remained unchanged. This indicates that the high overlying strata above KS7 remained a basically stable state.

In addition to the above regular monitoring of DOFS in the borehole, real-time monitoring of the surface subsidence near the mouth of ZY1 was conducted and the corresponding curve was obtained, as shown in Figure 5. As can be seen, when the working face advanced past ZY1 by 380m, the corresponding surface subsidence data is approximately 588mm. Since the mining height of LW205 is 10m, the corresponding subsidence ratio is about 0.058, indicating that the overlying strata is not broken yet and in a state of bending deformation.

4. Numerical simulation

4.1. Simulation schemes

Based on the columnar distribution of rock strata in the three panels and actual geological data, the numerical model was established by using 3DEC discrete element software. By using discrete element method, the rock mass is treated as blocks which can produce displacement and torsion motion. In this way, 3DEC discrete element software can effectively simulate the phenomena of fracture, separation, deformation and failure of surrounding rock. According to the key stratum theory, the key stratum controls the deformation and fracture of the whole overlying rock. Therefore, the key stratum in Figure 3 is retained in the model, while other rock strata are replaced by homogenized soft rock. In order to facilitate the simulation, the strike dimensions of the working face in tpanel 2 and panel 3 in Figure 1 are uniformly processed, and the plane layout of the working face in the model is shown in Figure 6.

A 3D numerical model (4100m×2600m×500m) was established (as shown in Figure 7a) to simulate the mining process of three panels and analyze the evolution law of mining-induced stress in the coal seam. In the mining process, one working face is excavated each time, and the excavation sequence is as follows: panel 1, panel 2, LW302, and LW301. After the excavation is over, the model will be calculated to reach a state of equilibrium again. In order to conduct a comparative analysis of the influence law of key stratum on the evolution of mining-induced stress, the block and joint strength of key stratum is weakened on the basis of the model in Figure 7a, and a weakened model without key stratum is obtained, as shown in Figure 7b.

Since the purpose of this study is to obtain the stress evolution law of the surrounding rock and the failure law of roof strata in the mining process, the thickness of floor in the model is slightly smaller. At the same time, in order to eliminate the boundary effect, a distance of more than 600m is reserved on both sides of the working face. Given that the burial depth of the working face is nearly one kilometer, in order to facilitate the model calculation, the top boundary of the model is constrained by uniform load on the premise that the simulation results are not to be affected. To be specific, a vertical stress of 13.2 MPa is applied in accordance with the burial depth of the model. Meanwhile, the horizontal and vertical constraints are applied to the surrounding boundary and the bottom boundary respectively, and their displacement is fixed to 0. Based on the geological histogram of the area under study, the numerical model is built layer by layer from the bottom boundary according to the actual thickness and lithology of the rock (coal) layer. This model design is more complex in structure, but it can simulate the movement of the roof of the coal seam more realistically and accurately.

In order to make the simulation results more aligned with the actual situation, the inversion modeling is carried out based on the field measured surface subsidence data. The simulated surface subsidence curve of LW205 after mining in Model 7(a) is shown in Figure 8. When the working face is pushed 380m past the observation point, the surface subsidence is 590mm, which is basically consistent with the measured results in Figure 5, indicating that the actual model established is in line with the actual mining conditions, thus ensuring the reliability of subsequent simulation results.

4.2. The evolution law of mining-induced stress in panel 1 and panel 2

With the progress of working face, the equilibrium state of rock mass around the stope is broken, and the gravity of overlying rock strata in the goaf is transferred to the surrounding coal body. Within a certain range around the stope, the vertical stress of the coal seam increases significantly, forming stress concentration. In the actual mining process, when the working faces of LW201~LW204 were mined, many rock burst accidents occurred in the main roadway in the isolated coal pillar area between panel 1 and panel 2.

In order to analyze the stress distribution of the surrounding rock of the stope in the mining process, the vertical stress distribution cloud map (Figure 9) of the coal seam was derived for analysis after the mining of LW204.

(As is shown in Figure 9a, the bearing stress around the goaf increases to a large extent after the working face in the first panel and LW201~204 working faces are mined. Comparative speaking, the stress increases mildly around the goaf in the first panel, whereas the stress increases significantly in the second panel. Figure 9b shows that the stress increase in the coal pillar area between the first panel and the second panel is not noticeable after the key stratum is weakened.

In order to further explain the distribution law of mining-induced stress in the first and second panels, measuring line A-A was set in the first panel and the second panel to monitor the vertical stress data. The vertical stress curve of LW204 after mining working face was obtained, as shown in Figure 10.

As is shown in Figure 10a, when panel 1 and LW201~204 in panel 2 are mined, the vertical stress in the isolated coal pillar area between panel 1 and panel 2 appears trapezoidal distribution with one side high and the other side low; and the stress at the edge of the goaf in panel 2 is obviously higher than that at the edge of panel 1. In the isolated coal pillar area between panel 1 and panel 2, the stress increases significantly. The average stress is about 43MPa, and the stress concentration coefficient reaches 1.79. When the key stratum is weakened, the stress increase in the isolated coal pillar area between panel 1 and panel 2 is not obvious (as shown in Figure 10b), the average stress in this area is about 32Mpa, and the stress concentration coefficient reaches 1.33.

The above simulation results show that the key stratum existing in the overlying strata has a significant influence on the mining-induced stress evolution of the coal seam. Moreover, the results also reveal the mechanism of rock burst accidents in the main roadway in the isolated coal pillar area between panel 1 and panel 2 after the mining of LW201-204 working faces.

4.3. The evolution law of mining-induced stress in panel 2 and panel 3

With the increase of goaf area, the stress of coal seam around the goaf further increases. Especially when panel 2 is mined, the goaf is virtually one kilometer wide and the whole panel is basically in a full-mining state. As the mining continues in panel 3, the stress in the isolated coal pillar area between panel 2 and panel 3 witnesses a significant change. Therefore, the vertical stress distribution cloud map and stress curve of coal seam are derived for further analysis after LW205 and LW302 are mined. Figure 11 and Figure 12 show the cloud map and stress curve corresponding to the two working faces respectively.

Theoretically, the stress influence boundary should be divided according to the standard that the bearing stress should be more than 1.05 time of the original stress. As can be seen from Figure 11, when LW205 in panel 2 is mined, the influence range of the lateral bearing stress is 660m. However, the influence range of lateral bearing stress in the weakened model is reduced to 400m.

As can be seen from Figure 12, due to the mining of panel 3, the vertical stress in the isolated coal pillar area between panel 2 and panel 3 rises and a high stress zone is formed noticeably. As LW302 is mined, stress concentration occurs in the isolated coal pillar area between panel 2 and panel 3. The stress in the area is generally over 40 MPa, with the peak value reaching 60.24 MPa and the stress concentration coefficient being as high as 2.51. This indicates that when LW302 is mined, its adjacent working face LW301 is in a state of high stress. In other words, the risk of rock burst is rather high in the mining process. However, in the weakened model, no obvious stress concentration occurs in LW301 and the stress in the central part of isolated coal pillar area between panle 2 and panel 3 does not increase significantly.

With the mining of LW301, the stress in the isolated coal pillar area between panle 2 and panel 3 would display a further rising trend. After LW301 reaches a state of equilibrium, the vertical stress distribution cloud map and stress curve of coal seam are derived, as shown in Figure 13.

As is shown in Figure 13a, due to the mining of LW301, the stress in the isolated coal pillar area between panle 2 and panle 3 displays a relatively remarkable rising trend. The maximum stress in the central part of isolated coal pillar area between panle 2 and panle 3 reaches up to 72 MPa, with the stress concentration coefficient being as high as 2.92. This peak value is far higher than the vertical stress in other areas around the goaf in panel 3. As can be seen from Figure 13b, when the key stratum is weakened, no obvious high stress zone is observed near the coal pillar side in LW301.

Figure 13c shows that the stress in the coal pillar area increases further after LW301 is mined. The average vertical stress upon the isolated coal pillars between panel 2 and panel3 increases by 17MPa. In particular, the stress of LW301 return airway near the coal pillar side increases by about 30MPa compared toLW302 after mining. As can be known from Figure 13d, when the key stratum is weakened, stress concentration is not obvious in LW301 return airway. This indicates that due to the influence of the key stratum of overlying strata, the rock burst risk of the return airway along the coal pillar side is greater in the mining process of LW301. Accordingly, the anti-shock and pressure relief measures should be taken in advance to ensure the safety during the mining period. As is revealed through field investigation, during the mining process of LW301, several roof-fall events occurred in the return airway, which verifies the correctness of the simulation results.

5. Conclusions

(1) The ground observation borehole was arranged in the test mine; and DOFS and the surface subsidence GNSS monitoring system were installed inside and at the mouth of the borehole respectively. According to the monitoring data on strain obtained by DOFS, the height of the broken stratum in the overlying rock after mining is much lower than that of the main key stratum, indicating that through fracture does not occur in the high key strata. As the monitoring results of the GNSS monitoring system show, the surface subsidence volume is 0.058m, and the subsidence coefficient is less than 0.1. This verifies again that the high key strata don’t break and enters a state of bending subsidence instead.

(2) By comparing the simulation results, it is found that when there are multiple key strata in the overlying strata, the stress concentration on the isolated coal pillar between panel 1 and panel 2 is greatly affected after the mining of LW201-LW204. This reveals the working mechanism of rock burst accidents in the main roadway in the isolated coal pillar area between panel 1 and panel 2.

(3) By comparing the simulation results, it is predicted that when there are multiple key strata in the overlying strata, the stress concentration risk of LW301 return airway along the coal pillar between panel 2 and panel 3 is relatively higher compared to LW302 after mining. After the actual mining, several roof-fall events occurred in the return airway, which verifies the accuracy of the predicted results.

(4) The paper reveals the mechanism of rock burst accidents caused by the significant influence of multiple key strata in the overlying strata on stress concentration in coal seam. Moreover, it predicts the potential hazard areas after mining. The research results are of great guiding significance for the prevention and control of rock burst accidents.