# Key Strata Inducing Dynamic Disasters Based on Energy Condition: Criterion and Application

Authors: Chao Xu; Gang Yang; Haoshi Sun; Liangliang Qin; Kai Wang; Bo Ren; Zhijie Wen
Journal: Geofluids
 (2021)

---

## Abstract

The thick and hard rock strata (THRS) exist widely in coal measure strata, which control the movements of overlying rock strata in stopes. When THRS break, great energy is released, which could aggravate the risks of coal and gas outburst, rock burst, and other dynamic disasters. Therefore, the foundation and key of preventing dynamic disasters are to distinguish the THRS that could induce coal-rock dynamic disasters and to analyze the laws of rock stratum breaking and energy releasing. The paper proposed the theoretical calculation methods of the energy accumulation and attenuation of rock breaking which is greatly affected by the hanging length of rock strata and the spreading distance. One or more roof strata that play a leading role in inducing dynamic disasters of the underlying coal mass are defined as the key disaster-inducing strata (KDIS). The disaster-inducing coefficient (DIC) is defined and used as the criterion of KDIS. The greater the source energy, the shorter the spreading distance, and the smaller the attenuation coefficient are, the easier the roof strata are to become KDIS. The disaster-inducing ability of the main THRS was analyzed, and the igneous sill was judged as KDIS, taking the Yangliu Coal Mine as project background. The breaking laws of the igneous sill were obtained by the methods of UDEC numerical simulation and microseismic monitoring, which verified the criterion of KDIS.


---

## Body

## 1. Introduction

With the huge demand of economic development for energy and the gradual exhaustion of shallow coal resources, the depth of coal mining is increasing by 10-25 m per year in China [1, 2]. The deep mining conditions are more complex with complicated geological structure, higher in situ stress, higher gas pressure, and lower permeability and strength of coal [3]. Therefore, the typical dynamic disasters such as rock burst and coal and gas outburst in deep mines are more and more frequent and serious, compared with those in shallow mines [4–8]. Thick and hard rock strata (THRS) exist widely in coal measure strata, which control the movements of overlying rock strata in stopes because of the great thickness, hardness, and integrity. In case of THRS failure and instability, huge elastic energy is released and transmitted to the stopes, which aggravates the risk of coal and rock dynamic disasters (Figure 1) [9, 10].Figure 1
Schematic diagram of dynamic disasters induced by key strata.In recent years, a series of coal and rock dynamic disasters occurred in the environment of THRS and caused huge casualties and property losses. For example, a rock burst accident that killed six miners was triggered by the fracturing of a 550 m-thick conglomerate rock stratum in the Qianqiu coal mine on March 27, 2014 [11]. One person died in a coal and gas outburst accident under a 120 m-thick igneous rock stratum in the Haizi coal mine on April 25, 2009 [12]. A water-gas ejection accident occurred with about 7800 m3 of water and 166400 m3 of gas which were ejected from the bed-splitting space under the 43.6 m-thick igneous rock stratum in the Yangliu coal mine on July 17, 2011 [13]. These examples highlight the need to investigate the effects of the THRS failure on the coal-rock dynamic disasters.The strata are composed of coal and rock layers with uneven thickness and different strength. An increasing number of studies have focused on the recognizing and failure behaviors of THRS [14]. The key strata [15], which are usually thick-hard strata, play an important role in overlying movement and may influence the mining-induced strata behaviors in the working face using top-coal caving method [16]. Based on the key strata theory, an overburden caving model is proposed to predict the multilayered hard strata behavior [17]. The concept of key strata in immediate roof (KSIR) was put forward, and the dynamic effect and control mechanism of KSIR were studied [18]. Three kinds of structural model, which are affected by the relative position of key strata in the overburden, are found and defined [19]. A roof structural model for large mining-height stopes was established based on their overlying rock structure [20]. In situ measurements via vertical boreholes were performed to determine the key strata, and an innovative solution to fracture high-level hard strata by ground hydraulic fracturing was proposed [21, 22]. A physical modeling system was established to explore the fracture mechanism of the hard roof, and the characteristics of acoustic emission signals during the process of hard roof failure were also studied [23]. The key strata breakage could bring great effects on the stress redistribution and energy storage of mining coal seams. The cause and mechanism of localized stress concentration and rock failure behavior in rock interlayer after protective seam exploitation were explored by numerical simulation, mechanical model, and field investigation [24]. The microseismic effects of hard and thick igneous strata separation and fracturing and the corresponding evaluation index of fracturing intensity were analyzed [9]. For controlling the behavior of the high-located main key stratum (HMKS), a new technique is proposed to add grout between the separated beds during mining [25]. However, there are less literatures on how to distinguish the key strata inducing coal-rock dynamic disasters.In this paper, the laws of energy accumulation, spreading, and attenuation for the THRS instability was studied theoretically, the sensitivity of those influencing factors was analyzed, and then the concept and criterion of key disaster-inducing strata (KDIS) were put forward. Afterwards, the KDIS in Yangliu Coal Mine, Huaibei Coalfield, China, was distinguished, in which results were verified by the numerical simulations and microseismic monitoring.


## 2. Energy Spreading Laws of THRS Instability

The elastic energy accumulated in THRS will be released rapidly and then spread and attenuated in the surrounding rock mass when THRS breaks. The residual energy could have a certain degree of impact action on the underlying coal under the state of ultimate stress, easily causing the irreversible damage and even the dynamic disasters.


### 2.1. Theoretical Analysis for Energy Spreading Laws

According to the stress state of roof strata, the total of the accumulated energyE0 is composed of the bending elastic energy, the volumetric strain energy, and the kinetic energy generated by rock strata movement [26], which is expressed by
(1)E0=q2a05576EI+1−2μ1+2k2q26E+12mdudt2,where q is the overburden load, Pa; a0 is hanging length of roof strata, m; E is the elasticity modulus, Pa; μ is Poisson’s ratio; k is the ratio of horizontal stress to vertical stress; m is the mass of broken roof strata, kg; u is the displacement of roof strata movement, m; and I is the inertia moment of roof strata, m4. For rectangular cross-section, I=bh3/12, in which b and h are the width and height of the rectangle, respectively.It is known that the total of the accumulated energy before the failure of THRS is positively correlated with the hanging length, overburden load, and kinematic velocity of THRS. The larger the hanging length, overburden load, and kinematic velocity, the greater the accumulated energy.The energy will be released rapidly and then spread and attenuated in the surrounding rock mass when THRS breaks. However, only about 1%~10% of the energy is spread in the form of seismic waves while most of the energy is converted to heat during collisions of THRS with other rocks [27]. The energy of seismic waves decreases gradually in the form of power with the increase of spreading distance due to the heterogeneity and damping effect of rock strata. From the perspective of safety protection, it could be assumed that 10% of the energy is converted into seismic waves, and then the residual energy after attenuation E0′ is
(2)E0′=0.1E0l−λ,where l is the spreading distance, m, and λ is the energy attenuation coefficient that is related to the properties of the medium.The residual energy decreases with the increasing of the energy spreading distance and attenuation coefficient. Assuming that the seismic waves are spread inm strata with the thickness of h1, h2, h3, …, hm, in turn, the residual energy Em′ would be
(3)Em′=0.1E0∏1mhi−λi.It indicates that the THRS with the hanging lengtha0 accumulates the energy E0; when THRS breaks, the energy is spread in m strata with the thickness of h1, h2, h3, …, hm in turn and attenuated to Em′.


### 2.2. Sensitivity of Influencing Factors for Energy Spreading

The energy accumulation and spread laws are related to many influencing factors, such as the overburden load, the hanging length, thickness, elasticity modulus and mass of roof strata, energy spreading distance, and attenuation coefficient based on the above formula. The sensitivity of influencing factors was obtained when they were assigned the different values (1 time, 2 times, and 4 times), assuming that the seismic waves were spread in the isotropous rock strata (Figure2). When the load is increased by 2 times and 4 times, the initial seismic wave energy increases by 1.22 times and 2.08 times, respectively; however, when the hanging length is increased by 2 times and 4 times, the initial seismic wave energy increases by 2.55 times and 23.04 times, respectively. When the thickness and elasticity modulus of rock strata increase, the initial seismic wave energy decreases gradually, which is contrary to the common understanding. It is because that only the single influence of rock thickness or elasticity modulus on initial seismic wave energy is considered. In fact, the rock strata are more difficult to be broken and form larger hanging length when the thickness and elasticity modulus of rock strata increase and so the greater seismic wave energy will be released when the rock strata are broken [28]. The attenuation velocity of seismic wave energy increases with the increase of energy attenuation coefficient. If the energy attenuation coefficients are set as 1.1, 1.3, and 1.5, the seismic wave energy will be attenuated to 1/14, 1/22.6, and 1/36.5 of the initial values, respectively, when the seismic waves are spread 10 m. The seismic wave energy is attenuated in the form of powers with the increasing of spreading distance. The seismic wave energy will be attenuated to 1/33.8, 1/86.8, and 1/209.4 of the initial values, respectively, when the spreading distances increase by 2 times and 4 times from 15 m. The comprehensive analysis shows that the spread of seismic wave energy is greatly affected by the hanging length of rock strata and the spreading distance.Figure 2
Sensitivity of influencing factors for energy spreading and attenuation. (a) Overburden load. (b) Hanging length of roof strata. (c) Thickness of roof strata. (d) Elasticity modulus of roof strata. (e) Mass of roof strata. (f) Energy attenuation coefficient.
(a)(b)(c)(d)(e)(f)


## 2.1. Theoretical Analysis for Energy Spreading Laws

According to the stress state of roof strata, the total of the accumulated energyE0 is composed of the bending elastic energy, the volumetric strain energy, and the kinetic energy generated by rock strata movement [26], which is expressed by
(1)E0=q2a05576EI+1−2μ1+2k2q26E+12mdudt2,where q is the overburden load, Pa; a0 is hanging length of roof strata, m; E is the elasticity modulus, Pa; μ is Poisson’s ratio; k is the ratio of horizontal stress to vertical stress; m is the mass of broken roof strata, kg; u is the displacement of roof strata movement, m; and I is the inertia moment of roof strata, m4. For rectangular cross-section, I=bh3/12, in which b and h are the width and height of the rectangle, respectively.It is known that the total of the accumulated energy before the failure of THRS is positively correlated with the hanging length, overburden load, and kinematic velocity of THRS. The larger the hanging length, overburden load, and kinematic velocity, the greater the accumulated energy.The energy will be released rapidly and then spread and attenuated in the surrounding rock mass when THRS breaks. However, only about 1%~10% of the energy is spread in the form of seismic waves while most of the energy is converted to heat during collisions of THRS with other rocks [27]. The energy of seismic waves decreases gradually in the form of power with the increase of spreading distance due to the heterogeneity and damping effect of rock strata. From the perspective of safety protection, it could be assumed that 10% of the energy is converted into seismic waves, and then the residual energy after attenuation E0′ is
(2)E0′=0.1E0l−λ,where l is the spreading distance, m, and λ is the energy attenuation coefficient that is related to the properties of the medium.The residual energy decreases with the increasing of the energy spreading distance and attenuation coefficient. Assuming that the seismic waves are spread inm strata with the thickness of h1, h2, h3, …, hm, in turn, the residual energy Em′ would be
(3)Em′=0.1E0∏1mhi−λi.It indicates that the THRS with the hanging lengtha0 accumulates the energy E0; when THRS breaks, the energy is spread in m strata with the thickness of h1, h2, h3, …, hm in turn and attenuated to Em′.


## 2.2. Sensitivity of Influencing Factors for Energy Spreading

The energy accumulation and spread laws are related to many influencing factors, such as the overburden load, the hanging length, thickness, elasticity modulus and mass of roof strata, energy spreading distance, and attenuation coefficient based on the above formula. The sensitivity of influencing factors was obtained when they were assigned the different values (1 time, 2 times, and 4 times), assuming that the seismic waves were spread in the isotropous rock strata (Figure2). When the load is increased by 2 times and 4 times, the initial seismic wave energy increases by 1.22 times and 2.08 times, respectively; however, when the hanging length is increased by 2 times and 4 times, the initial seismic wave energy increases by 2.55 times and 23.04 times, respectively. When the thickness and elasticity modulus of rock strata increase, the initial seismic wave energy decreases gradually, which is contrary to the common understanding. It is because that only the single influence of rock thickness or elasticity modulus on initial seismic wave energy is considered. In fact, the rock strata are more difficult to be broken and form larger hanging length when the thickness and elasticity modulus of rock strata increase and so the greater seismic wave energy will be released when the rock strata are broken [28]. The attenuation velocity of seismic wave energy increases with the increase of energy attenuation coefficient. If the energy attenuation coefficients are set as 1.1, 1.3, and 1.5, the seismic wave energy will be attenuated to 1/14, 1/22.6, and 1/36.5 of the initial values, respectively, when the seismic waves are spread 10 m. The seismic wave energy is attenuated in the form of powers with the increasing of spreading distance. The seismic wave energy will be attenuated to 1/33.8, 1/86.8, and 1/209.4 of the initial values, respectively, when the spreading distances increase by 2 times and 4 times from 15 m. The comprehensive analysis shows that the spread of seismic wave energy is greatly affected by the hanging length of rock strata and the spreading distance.Figure 2
Sensitivity of influencing factors for energy spreading and attenuation. (a) Overburden load. (b) Hanging length of roof strata. (c) Thickness of roof strata. (d) Elasticity modulus of roof strata. (e) Mass of roof strata. (f) Energy attenuation coefficient.
(a)(b)(c)(d)(e)(f)


## 3. Concept and Criterion of KDIS

### 3.1. Concept of KDIS

Although most of the released energy has been dissipated, the residual energy may still induce dynamic disasters when the seismic waves were spread to the mining coal seam. Therefore, one or more roof strata that play a leading role in inducing dynamic disasters of the underlying coal mass are defined as the key disaster-inducing strata (KDIS). The residual energy not only depends on the size and mechanical properties of rock strata but also relates to the energy spreading distance and attenuation coefficient. The greater the source energy, the shorter the spreading distance, and the smaller the attenuation coefficient are, the stronger the disaster-inducing ability of the roof strata is, and the easier the roof strata are to become KDIS. However, the key stratum of ground control mainly depends on the size and mechanical properties of the rock stratum itself [15]. Thus, the key stratum of ground control is not exactly equal to KDIS.Figure3 provides a better understanding of KDIS. When the rock stratum #1 that is close to the mining coal seam is broken, the residual energy may be larger than the minimum energy inducing dynamic disasters although the released energy is less. When the thicker rock stratum #2 is broken, more seismic wave energy is released but cannot induce dynamic disasters due to the relatively large distance from the mining coal seam. But for the thick rock stratum #3 which is the key stratum of ground control, great seismic wave energy is released when broken, and the residual energy is so large that could induce dynamic disasters although it is far from the mining coal seam. Therefore, the rock strata #1 and #3 are the KDIS, and the key stratum of ground control may be the KDIS under certain conditions.Figure 3
Energy attenuation laws for roof rock breakage against distance.


### 3.2. Criterion of KDIS

During the process of dynamic disaster occurring, the crushing work is needed for crack propagation, debris stripping, and new surface forming of coal mass; the kinetic energy is needed for throwing the broken coal into the excavation space; the friction dissipation energy is also needed for overcoming the inner friction during the broken coal peeling and throwing. Thus, the total energy consumption is composed of the crushing work, the kinetic energy, and the friction dissipation energy. If the accumulated energy of mining coal mass is larger than the total energy consumption, the dynamic disasters will occur.They are assumed that there aren overlying rock strata with the distances l1, l2, l3, …, ln from the mining coal seam, respectively, and that the seismic wave energy E1, E2, E3, …, En will be released and attenuated with the average attenuation coefficient λ0 when they are broken. The residual energy E1′, E2′, E3′, …, En′ when spreading to the mining coal seam is, respectively
(4)E1′=E1l1−λ0E2′=E2l2−λ0E3′=E3l3−λ0⋯En′=Enln−λ0.The ratio of the residual energyEi′ and the minimum energy inducing dynamic disasters Emin is defined as the disaster-inducing coefficient (DIC) ξ, that is ξ=Ei′/Emin. Based on DIC, THRS are divided into three types: KDIS, the weak disaster-inducing stratum (WDIS), and the non-disaster-inducing stratum (NDIS). If ξ≥1, the broken rock strata could induce dynamic disasters, which are distinguished to be KDIS; if 0.1≤ξ<1, the broken rock strata may induce dynamic disasters considering in situ stress, mechanical properties of coal, and other factors, which are distinguished to be WDIS; if ξ<0.1, the broken rock strata could not induce dynamic disasters, which are distinguished to be NDIS.


## 3.1. Concept of KDIS

Although most of the released energy has been dissipated, the residual energy may still induce dynamic disasters when the seismic waves were spread to the mining coal seam. Therefore, one or more roof strata that play a leading role in inducing dynamic disasters of the underlying coal mass are defined as the key disaster-inducing strata (KDIS). The residual energy not only depends on the size and mechanical properties of rock strata but also relates to the energy spreading distance and attenuation coefficient. The greater the source energy, the shorter the spreading distance, and the smaller the attenuation coefficient are, the stronger the disaster-inducing ability of the roof strata is, and the easier the roof strata are to become KDIS. However, the key stratum of ground control mainly depends on the size and mechanical properties of the rock stratum itself [15]. Thus, the key stratum of ground control is not exactly equal to KDIS.Figure3 provides a better understanding of KDIS. When the rock stratum #1 that is close to the mining coal seam is broken, the residual energy may be larger than the minimum energy inducing dynamic disasters although the released energy is less. When the thicker rock stratum #2 is broken, more seismic wave energy is released but cannot induce dynamic disasters due to the relatively large distance from the mining coal seam. But for the thick rock stratum #3 which is the key stratum of ground control, great seismic wave energy is released when broken, and the residual energy is so large that could induce dynamic disasters although it is far from the mining coal seam. Therefore, the rock strata #1 and #3 are the KDIS, and the key stratum of ground control may be the KDIS under certain conditions.Figure 3
Energy attenuation laws for roof rock breakage against distance.


## 3.2. Criterion of KDIS

During the process of dynamic disaster occurring, the crushing work is needed for crack propagation, debris stripping, and new surface forming of coal mass; the kinetic energy is needed for throwing the broken coal into the excavation space; the friction dissipation energy is also needed for overcoming the inner friction during the broken coal peeling and throwing. Thus, the total energy consumption is composed of the crushing work, the kinetic energy, and the friction dissipation energy. If the accumulated energy of mining coal mass is larger than the total energy consumption, the dynamic disasters will occur.They are assumed that there aren overlying rock strata with the distances l1, l2, l3, …, ln from the mining coal seam, respectively, and that the seismic wave energy E1, E2, E3, …, En will be released and attenuated with the average attenuation coefficient λ0 when they are broken. The residual energy E1′, E2′, E3′, …, En′ when spreading to the mining coal seam is, respectively
(4)E1′=E1l1−λ0E2′=E2l2−λ0E3′=E3l3−λ0⋯En′=Enln−λ0.The ratio of the residual energyEi′ and the minimum energy inducing dynamic disasters Emin is defined as the disaster-inducing coefficient (DIC) ξ, that is ξ=Ei′/Emin. Based on DIC, THRS are divided into three types: KDIS, the weak disaster-inducing stratum (WDIS), and the non-disaster-inducing stratum (NDIS). If ξ≥1, the broken rock strata could induce dynamic disasters, which are distinguished to be KDIS; if 0.1≤ξ<1, the broken rock strata may induce dynamic disasters considering in situ stress, mechanical properties of coal, and other factors, which are distinguished to be WDIS; if ξ<0.1, the broken rock strata could not induce dynamic disasters, which are distinguished to be NDIS.


## 4. Engineering Application

### 4.1. Project Background

Yangliu Coal Mine, the coal and gas outburst mine, is located in the north Anhui province, China. The coal seam #10 with the average thickness of 3.19 m is the main mining seam adopting the comprehensive mechanized mining method. The magmatic activities in Yangyan Coal Mine were quite intense, and almost all coal seams were affected by magmatic intrusion. The igneous sill, 43.6 m thick averagely and 102 m above the coal seam #10, is the key stratum of ground control (Figure4). The panel #10414 and the adjacent panel #10416 of the coal seam #10 were mined in turn. According to the statistics, multiple dynamic disasters were induced by the movement and failure of the igneous sill during the two panels mining. The large amount of water and gas were ejected from the surface well when the panel #10414 was advanced 527 m [13]; the abnormal pressure on working face occurred when the panel #10416 was advanced 320 m [28].Figure 4
Project background of Yangliu Coal Mine.
(a)(b)


### 4.2. Criterion Results of KDIS

It can be preliminarily analyzed from the column chart of coal-bearing strata that the four thick rock strata, 12 m, 21 m, 72 m, and 102 m above the coal seam #10, respectively, may have the ability to induce dynamic disasters. The released energy of the four rock strata breaking and the residual energy spreading to the coal seam #10 were calculated based on the above formula (Table1 and Figure 5). Researches show that the minimum energy inducing dynamic disasters Emin is 104 orders of magnitude [29]. Assuming that the minimum energy inducing dynamic disasters in the coal seam #10 of Yangliu Coal Mine is 10 kJ for security, the DIC of the four rock strata were 0.155, 0.162, 0.028, and 1.598, respectively. So the rock strata #1 and #2 were judged to be WDIS, the stratum #3 to be NDIS, and the stratum #4 to be KDIS. The breaking of the stratum #4 could induce the dynamic disasters in coal seam #10.Table 1
Criterion results of KDIS in Yangliu Coal Mine.
No.#1#2#3#4LithologySandstoneSiltstoneSandstoneIgneous rockThickness (m)4.811.59.343.6Distance from coal seam #10 (m)122172102Extreme hanging length (m)16.720.018.7184.6Released energy (MJ)0.310.630.4841.10Residual energy spreading to coal seam #10 (kJ)1.551.620.2815.98DIC0.1550.1620.0281.598Criterion resultWDISWDISNDISKDISFigure 5
Energy attenuation laws of rock strata in Yangliu Coal Mine.


### 4.3. Result Validation

The methods of UDEC numerical simulation and microseismic monitoring were adopted for verifying the THRS failure and energy releasing laws.


#### 4.3.1. UDEC Numerical Simulation

The numerical model is 190 m high from 50 m below the coal seam #10 to the top of the igneous sill. The model size is500∗190m while the panel size is 300∗3m and the coal pillars with 100 m width were set at each end (Figure 6). The Mohr-Coulomb criterion was selected, and the mechanical parameters of the coal and rock masses and joint surface were set as shown in the reference [26]. The bottom, left, and right boundaries of the model were set to obey the zero-displacement constraint, and the top boundary was loaded at 10 MPa. The excavation of the panel took place in six steps, each of which was 50 m.Figure 6
UDEC numerical model.The movement velocity and acceleration of the igneous sill and the kinetic energy in the system were recorded by the ordersyvel, yacc, and energy, respectively (Figure 7). The movement velocity and acceleration of the igneous sill at the central position (250, 155) increased gradually with the workface advancing. When the workface was advanced 250 m, the movement velocity increased to the maximum of -2.6 m/s, the acceleration increased to about 3000 m/s2, and the kinetic energy extended to the maximum of 59.3 MJ, which were consistent with the previous calculated results. When the workface advanced 300 m, the maximum of the movement velocity was 1.2 m/s and the acceleration increased to -9600 m/s2 while the kinetic energy decreased to the 56.2 MJ due to the collision and rebound between the igneous sill and the underlying rock strata. It was indicated that the igneous sill had been broken when the workface was advanced 250~300 m, and the released kinetic energy was much larger than that of the common rock strata. The igneous sill would force greater impact on the underlying coal seam and conformed to the characteristics of KDIS, which verified the above research results.Figure 7
Movement laws of igneous sill.


#### 4.3.2. Microseismic Monitoring

The microseismic monitoring technology is a method to monitor the stability of engineering rock mass using the vibration waves emitted by the rock mass deformation and failure [30–32]. The SOS microseismic monitoring system was used to monitor the microseismic signals in real time in Yangliu Coal Mine. The daily frequency and maximum energy of seismic waves during the workface #10416 advancing are shown in Figure 8. When the workface #10416 was advanced 300~380 m, the igneous sill was broken with square form [28], and thus, the daily frequency and maximum energy of seismic waves increased obviously, the maximum of which were 121.1 kJ and 23 times, relatively, according with the characteristics of KDIS. The daily frequency and maximum energy returned to normal levels gradually with the workface #10416 advancing further. A large number of microseismic events with greater magnitude mainly occurred on the side of the intake roadway (Figure 9), which verified the laws of the igneous sill breaking and energy releasing.Figure 8
Daily frequency and maximum energy of the panel #10416.Figure 9
Microseismic event distribution.


## 4.1. Project Background

Yangliu Coal Mine, the coal and gas outburst mine, is located in the north Anhui province, China. The coal seam #10 with the average thickness of 3.19 m is the main mining seam adopting the comprehensive mechanized mining method. The magmatic activities in Yangyan Coal Mine were quite intense, and almost all coal seams were affected by magmatic intrusion. The igneous sill, 43.6 m thick averagely and 102 m above the coal seam #10, is the key stratum of ground control (Figure4). The panel #10414 and the adjacent panel #10416 of the coal seam #10 were mined in turn. According to the statistics, multiple dynamic disasters were induced by the movement and failure of the igneous sill during the two panels mining. The large amount of water and gas were ejected from the surface well when the panel #10414 was advanced 527 m [13]; the abnormal pressure on working face occurred when the panel #10416 was advanced 320 m [28].Figure 4
Project background of Yangliu Coal Mine.
(a)(b)


## 4.2. Criterion Results of KDIS

It can be preliminarily analyzed from the column chart of coal-bearing strata that the four thick rock strata, 12 m, 21 m, 72 m, and 102 m above the coal seam #10, respectively, may have the ability to induce dynamic disasters. The released energy of the four rock strata breaking and the residual energy spreading to the coal seam #10 were calculated based on the above formula (Table1 and Figure 5). Researches show that the minimum energy inducing dynamic disasters Emin is 104 orders of magnitude [29]. Assuming that the minimum energy inducing dynamic disasters in the coal seam #10 of Yangliu Coal Mine is 10 kJ for security, the DIC of the four rock strata were 0.155, 0.162, 0.028, and 1.598, respectively. So the rock strata #1 and #2 were judged to be WDIS, the stratum #3 to be NDIS, and the stratum #4 to be KDIS. The breaking of the stratum #4 could induce the dynamic disasters in coal seam #10.Table 1
Criterion results of KDIS in Yangliu Coal Mine.
No.#1#2#3#4LithologySandstoneSiltstoneSandstoneIgneous rockThickness (m)4.811.59.343.6Distance from coal seam #10 (m)122172102Extreme hanging length (m)16.720.018.7184.6Released energy (MJ)0.310.630.4841.10Residual energy spreading to coal seam #10 (kJ)1.551.620.2815.98DIC0.1550.1620.0281.598Criterion resultWDISWDISNDISKDISFigure 5
Energy attenuation laws of rock strata in Yangliu Coal Mine.


## 4.3. Result Validation

The methods of UDEC numerical simulation and microseismic monitoring were adopted for verifying the THRS failure and energy releasing laws.


### 4.3.1. UDEC Numerical Simulation

The numerical model is 190 m high from 50 m below the coal seam #10 to the top of the igneous sill. The model size is500∗190m while the panel size is 300∗3m and the coal pillars with 100 m width were set at each end (Figure 6). The Mohr-Coulomb criterion was selected, and the mechanical parameters of the coal and rock masses and joint surface were set as shown in the reference [26]. The bottom, left, and right boundaries of the model were set to obey the zero-displacement constraint, and the top boundary was loaded at 10 MPa. The excavation of the panel took place in six steps, each of which was 50 m.Figure 6
UDEC numerical model.The movement velocity and acceleration of the igneous sill and the kinetic energy in the system were recorded by the ordersyvel, yacc, and energy, respectively (Figure 7). The movement velocity and acceleration of the igneous sill at the central position (250, 155) increased gradually with the workface advancing. When the workface was advanced 250 m, the movement velocity increased to the maximum of -2.6 m/s, the acceleration increased to about 3000 m/s2, and the kinetic energy extended to the maximum of 59.3 MJ, which were consistent with the previous calculated results. When the workface advanced 300 m, the maximum of the movement velocity was 1.2 m/s and the acceleration increased to -9600 m/s2 while the kinetic energy decreased to the 56.2 MJ due to the collision and rebound between the igneous sill and the underlying rock strata. It was indicated that the igneous sill had been broken when the workface was advanced 250~300 m, and the released kinetic energy was much larger than that of the common rock strata. The igneous sill would force greater impact on the underlying coal seam and conformed to the characteristics of KDIS, which verified the above research results.Figure 7
Movement laws of igneous sill.


### 4.3.2. Microseismic Monitoring

The microseismic monitoring technology is a method to monitor the stability of engineering rock mass using the vibration waves emitted by the rock mass deformation and failure [30–32]. The SOS microseismic monitoring system was used to monitor the microseismic signals in real time in Yangliu Coal Mine. The daily frequency and maximum energy of seismic waves during the workface #10416 advancing are shown in Figure 8. When the workface #10416 was advanced 300~380 m, the igneous sill was broken with square form [28], and thus, the daily frequency and maximum energy of seismic waves increased obviously, the maximum of which were 121.1 kJ and 23 times, relatively, according with the characteristics of KDIS. The daily frequency and maximum energy returned to normal levels gradually with the workface #10416 advancing further. A large number of microseismic events with greater magnitude mainly occurred on the side of the intake roadway (Figure 9), which verified the laws of the igneous sill breaking and energy releasing.Figure 8
Daily frequency and maximum energy of the panel #10416.Figure 9
Microseismic event distribution.


## 4.3.1. UDEC Numerical Simulation

The numerical model is 190 m high from 50 m below the coal seam #10 to the top of the igneous sill. The model size is500∗190m while the panel size is 300∗3m and the coal pillars with 100 m width were set at each end (Figure 6). The Mohr-Coulomb criterion was selected, and the mechanical parameters of the coal and rock masses and joint surface were set as shown in the reference [26]. The bottom, left, and right boundaries of the model were set to obey the zero-displacement constraint, and the top boundary was loaded at 10 MPa. The excavation of the panel took place in six steps, each of which was 50 m.Figure 6
UDEC numerical model.The movement velocity and acceleration of the igneous sill and the kinetic energy in the system were recorded by the ordersyvel, yacc, and energy, respectively (Figure 7). The movement velocity and acceleration of the igneous sill at the central position (250, 155) increased gradually with the workface advancing. When the workface was advanced 250 m, the movement velocity increased to the maximum of -2.6 m/s, the acceleration increased to about 3000 m/s2, and the kinetic energy extended to the maximum of 59.3 MJ, which were consistent with the previous calculated results. When the workface advanced 300 m, the maximum of the movement velocity was 1.2 m/s and the acceleration increased to -9600 m/s2 while the kinetic energy decreased to the 56.2 MJ due to the collision and rebound between the igneous sill and the underlying rock strata. It was indicated that the igneous sill had been broken when the workface was advanced 250~300 m, and the released kinetic energy was much larger than that of the common rock strata. The igneous sill would force greater impact on the underlying coal seam and conformed to the characteristics of KDIS, which verified the above research results.Figure 7
Movement laws of igneous sill.


## 4.3.2. Microseismic Monitoring

The microseismic monitoring technology is a method to monitor the stability of engineering rock mass using the vibration waves emitted by the rock mass deformation and failure [30–32]. The SOS microseismic monitoring system was used to monitor the microseismic signals in real time in Yangliu Coal Mine. The daily frequency and maximum energy of seismic waves during the workface #10416 advancing are shown in Figure 8. When the workface #10416 was advanced 300~380 m, the igneous sill was broken with square form [28], and thus, the daily frequency and maximum energy of seismic waves increased obviously, the maximum of which were 121.1 kJ and 23 times, relatively, according with the characteristics of KDIS. The daily frequency and maximum energy returned to normal levels gradually with the workface #10416 advancing further. A large number of microseismic events with greater magnitude mainly occurred on the side of the intake roadway (Figure 9), which verified the laws of the igneous sill breaking and energy releasing.Figure 8
Daily frequency and maximum energy of the panel #10416.Figure 9
Microseismic event distribution.


## 5. Conclusions

In this paper, the laws of energy accumulation, spreading, and attenuation for the THRS instability were studied theoretically, the sensitivity of those influencing factors was analyzed, and then the concept and criterion of KDIS were put forward. Afterwards, the KDIS in Yangliu Coal Mine, Huaibei Coalfield, China, was distinguished, in which results were verified by the numerical simulations and microseismic monitoring. The main conclusions are as follows:(1)
The total of the accumulated energy is composed of the bending elastic energy, the volumetric strain energy, and the kinetic energy generated by rock strata movement, which is positively correlated with the hanging length, overburden load, and kinematic velocity. The energy of seismic waves decreases gradually in the form of power with the increase of spreading distance, the residual of which is negatively correlated with the energy spreading distance and attenuation coefficient. The spread of seismic wave energy is greatly affected by the hanging length of rock strata and the spreading distance, relatively(2)
One or more roof strata that play a leading role in inducing dynamic disasters of the underlying coal mass are defined as the KDIS. The greater the source energy, the shorter the spreading distance, and the smaller the attenuation coefficient are, the easier the roof strata are to become KDIS. The key stratum of ground control is not exactly equal to KDIS but may be the KDIS under certain conditions. The DIC is defined and used as the criterion of KDIS, based on which the THRS are divided into three types of KDIS, WDIS, and NDIS(3)
The disaster-inducing ability of the main THRS was analyzed and the igneous sill was judged as KDIS when the overlying coal seam #10 was mined, taking the Yangliu Coal Mine as project background. By the methods of UDEC numerical simulation and microseismic monitoring, it was revealed that the igneous sill was broken and greater energy than that of common rock strata was released resulting to the occurrence of abnormal pressure when the workface of the coal seam #10 was advanced 250~300 m, which verified the criterion of KDIS.



---

*Source: 6672020-2021-02-20.xml*