Abstract

The practice shows that deep strip mining induces rockburst disaster easily. Accurately evaluating rockburst risk of the strip coal pillar is of great significance for ensuring the safety of deep strip mining. In this paper, the catastrophe mechanics model was used to analyze the abrupt instability condition of strip coal pillar. And the three indicators that are the medium stiffness ratio () of the elastic and plastic zone in the coal pillar, the plastic zone width ratio (), and the elastic deformation index () of core zone were put forward with considering the geometry size of coal pillar. Based on the 3202 panel of Gucheng Coal Mine, the evolution characteristics of rockburst risk of coal pillar under different mining widths and coal pillar widths were studied by numerical simulation. The evaluation result shows that the strip coal pillar of the 3202 panel is in danger of strong rockburst, which is more in line with the actual situation than the results of the traditional rockburst tendency identification test and comprehensive index method. These three indicators can be regarded as important indicators to evaluate the rockburst risk in the strip mining engineering field. Based on that, the design principle of strip mining in Gucheng Coal Mine was put forward, which is considered an important reference for similar cases.

1. Introduction

In recent years, with the increasing depletion of shallow coal resources in China, in order to improve mining rate of coal resources and prolong mine life, many mining areas increase the exploitation of coal under buildings, railways, and water bodies, especially in the eastern part of China [15]. Strip mining is one of the common mining methods of coal under buildings, railways, and water bodies, which can effectively alleviate the surface subsidence of the coal mining area. Strip mining is widely used in the old eastern mining areas of China, which is one of the research hotspots in the Green Mining of Coal Resources in China. However, under mining disturbance, the high concentrated stress of the strip coal pillar is easy to induce dynamic disasters such as rockburst, which seriously affects the safety production of the working face [610].

As a complex nonlinear dynamic disaster, the occurrence of rockburst in coal and rock mass needs to meet certain physical conditions, such as the burst risk of coal and rock mass [11, 12]. It is of great significance for the prevention and control of rockburst to accurately evaluate the risk of rockburst [13]. In order to quantitatively measure the burst risk of coal and rock, domestic and overseas scholars have proposed a variety of impact propensity indexes based on energy storage, dynamic failure time, deformation, and stiffness of coal and rock and put forward corresponding discrimination indexes [1416]. Pan et al. [17] put forward three new indexes: impact energy velocity index considering time effect, critical softening region coefficient, and critical stress coefficient. Jiang et al. [1822] studied the influence of mining supporting pressure and structural stress on the rockburst and adopted a multiparameter early warning method to evaluate the rockburst risk of the working face. Lu et al. [23, 24] revealed the mining-induced fault slip mechanism and proposed the multiparameter precursors of a rockburst induced by the fault slip. Zhao et al. [25] established the prediction model of rockburst induced by roof. Yang et al. [26] studied the mechanism of rockburst from the vertical stress of coal pillar. Tan et al. [27] established a cusp catastrophe model for the instability of narrow coal pillars between mining tunnels and strip coal pillars and analyzed the instability mechanism and engineering stability conditions of the coal pillar system of Wongawilli mining. The above studies were based on the characteristics of stress, structure, and roof to establish the rockburst risk assessment model, mainly considering the influence of geological and mining factors on the occurrence of rockburst. However, there are relatively few studies on the dynamic instability characteristics and the main control factors of the strip coal pillar. The existing evaluation methods do not consider the special stress characteristics and instability conditions of the strip coal pillar and lack of pertinence in the evaluation of the rockburst risk of strip mining.

In this study, the medium stiffness ratio of the elastic and plastic zone in coal pillar, the plastic zone width ratio, and the elastic deformation index of kernel zone were put forward based on catastrophe theory, which could be used to guide the prevention and control of rockburst of deep strip working face.

2. Cusp Catastrophe Model of Strip Coal Pillar

Rockburst in the coal mine is characterized by abruptness, concealment, no forewarning, and heavy causalities and economic loss [28]. Rockburst is essentially a discontinuous process of sudden instability of the system. Therefore, catastrophe theory can be used to study the process of stope impact instability [29]. For strip coal pillar, with the increase of the width of goaf, the side abutment pressure increases significantly. Under the concentrated load, the coal pillar shows strain softening and remains in a stable position at the initial stage. With the compression deformation of the coal pillar, the deformation energy in the coal pillar gradually accumulates and gradually reaches the unstable critical state. Under the external disturbance, the critical system suddenly jumps to the corresponding stable equilibrium point. The jump phenomenon of the coal pillar system can be considered as the occurrence of rockburst [30].

In order to explore the process of coal pillar rockburst, Tan et al. [27] established the mechanical model of the coal pillar system as shown in Figure 1. In Figure 1, m is the mining thickness of the coal seam, is the reserved width of the coal pillar, and b is the mining width of the working face, all of which are in meters.

According to the stress-strain curve of the coal pillar, the constitutive relationship between elastic core and plastic zone of coal pillar is different. In order to simplify the calculation, the elastic core of the coal pillar is regarded as linear elastic, the stress-strain curve is a straight line, and the slope of the straight line is the elastic modulus of the coal seam. The plastic zone constitutive relation is approximately replaced by the strain-softening constitutive relation of coal. The stress-strain relationship of coal pillar can be expressed aswhere E is the elastic modulus of the coal seam, is the strain of coal pillar, and is the strain corresponding to the peak stress.

The width of the plastic zone of the coal pillar is defined as Y, and then the width of the elastic core of the coal pillar is . Considering the influence of coal pillar size, the medium stiffness of the elastic core and the plastic zone of the coal pillar can be expressed as follows:where is the medium stiffness in the elastic core of the strip coal pillar, is the medium stiffness in the plastic zone of the strip coal pillar, is the deformation value of the coal pillar, and is the deformation value of the coal pillar under peak load. Since the height of the coal pillar is , there is .

According to the research results of Tan et al. [27], the expression of total potential energy function of the system, denoted by , can be expressed as follows:where is the strain energy of the plastic zone of the strip coal pillar, is the elastic potential energy in the elastic nuclear region, and is the gravitational potential energy of the overlying strata.

The first derivative of equation (3) is expanded by Taylor series at , the cubic term is extracted, and the following is obtained:

Then, the standard equation of the cusp catastrophe equilibrium surface is obtained and can be expressed as equation (5). Point is the mutation point of the following equation [27]:

In equation (5), x is the state variable, p and q are control variables, and thuswhere is the ratio of the medium stiffness in the elastic core to the medium stiffness in the plastic zone of the coal pillar, which is stiffness ratio for short.

According to the mutation theory, only when , the system will mutate [31]. So the necessary condition for the system mutation is

The critical compression of coal pillar before rockburst is denoted by and can be expressed as follows:

According to equations (9) and (10), the internal cause of strip pillar rockburst is the stiffness ratio of the elastic-plastic zone which is less than 1.46 caused by strain softening. And the external cause of the coal pillar rockburst is that the deformation energy accumulation caused by the compression deformation of the coal pillar reaches the instability, and the compression of the elastic core of coal pillar reaches or approaches to the critical value. The strain-softening state reflects the failure degree of the coal pillar under load, while the elastic deformation index reflects the magnitude of energy accumulation. Therefore, it is feasible to evaluate the rockburst risk of the coal pillar by using the indexes of strain-softening degree and compression deformation.

3. Rockburst Risk Evaluation Index of Strip Coal Pillar

3.1. Medium Stiffness Ratio of Elastic and Plastic Zones in Strip Coal Pillar

At present, the commonly used stiffness ratio index for rockburst risk evaluation of coal and rock is the preyield to postyield stiffness ratio on the stress-strain curve of the coal and rock specimen, denoted by , which reflects the strain-softening degree of the specimen after yield [32]. According to the literature [32], when , there is rockburst tendency, and when , there is no rockburst tendency, and the detailed division is as follows: , no rockburst tendency; , weak and medium rockburst tendency; , strong rockburst tendency. However, this stiffness ratio index ignores the influence of strip coal pillar size effect, and the evaluation result of rockburst risk is not accurate.

Based on the old stiffness ratio index of rockburst risk evaluation and considering the influence of coal pillar geometry size, this paper proposes a new index, the medium stiffness ratio of the elastic and plastic zone in strip coal pillar, which is denoted by . The physical significance of this index is the ratio of the medium stiffness in the elastic core and the medium stiffness in the plastic zone of the loaded strip coal pillar, which indicates the degree of strain-softening during the loading process of the strip coal pillar. The smaller the of coal pillar, the greater the rockburst risk. The can be expressed as follows:

According to calculation results of equations (9) and (10), referring to index, the classification standard of plasticity zone width ratio index is as follows: , strong rockburst tendency; , medium rockburst tendency; , no rockburst tendency.

3.2. Plastic Zone Width Ratio of Strip Coal Pillar

According to equation (11), it can be predicted that the field application of the stiffness ratio index in the elastic-plastic zone of strip coal pillar is complicated.

By changing the form of equation (11), the necessary condition for rockburst of the strip coal pillar can be obtained as follows:where is the plastic zone width of strip coal pillar and is the width of strip coal pillar.

In the field application, the width of the plastic zone of the strip coal pillar can be measured by the drilling cuttings method. The larger the width ratio of the plastic zone is, the greater the rockburst risk is. The width ratio of the plastic zone is denoted by aY, which can be expressed as follows:

The classification standard of the plasticity zone width ratio index is as follows: , strong rockburst tendency; , medium rockburst tendency; , no rockburst tendency.

3.3. Elastic Deformation Index of the Strip Coal Pillar Core Zone

The elastic deformation index is also one of the commonly used evaluation indexes of rockburst risk of coal and rock. The old elastic deformation index is obtained by the ratio of the elastic deformation to the total deformation under the cyclic loading and unloading test, and the test loading is not less than 80% of the peak strength of coal and rock [33]. This index reflects the influence of the elastic deformation degree of the specimen on the burst failure. The detailed classification standard of the elastic deformation index which is usually denoted by is as follows: , strong rockburst tendency; , medium rockburst tendency; , no rockburst tendency. However, this index is aimed at small size specimens of coal and rock, and the impact of strip coal pillar geometric size on the rockburst risk is also not considered, resulting in its little significance in the application of strip coal pillar rockburst risk assessment.

Based on the old elastic deformation index and considering the influence of coal pillar geometry size on rockburst risk, this paper proposes a new elastic deformation index of strip coal pillar core zone, which is denoted by . is the ratio of elastic core compression to critical core area compression of loaded strip coal pillar. The elastic deformation index is directly related to the accumulation of deformation energy, which indirectly reflects the relationship between the deformation energy accumulation of strip coal pillar and the critical state of rockburst. The larger elastic deformation index of strip coal pillar core zone, the greater the rockburst risk of strip coal pillar. When the compression amount of strip coal pillar core zone exceeds the critical compression amount of impact instability, it can be considered that the rockburst has occurred. The new elastic deformation index can be expressed aswhere is the compression of the elastic core of strip coal pillar and is the rockburst critical compression of strip coal pillar core zone.

According to the old elastic deformation index of coal and rock, the classification standard of the new elastic deformation index of strip coal pillar core zone is as follows: , no rockburst tendency; , medium rockburst tendency; , strong rockburst tendency; , rockburst has occurred.

4. Application of Rockburst Risk Evaluation

4.1. Mining and Geological Condition

The Gucheng Coal Mine is located in Jining city, Shandong Province, China. The average overburden depth of Gucheng Coal Mine is 1070 m, which is a typical deep mine. The current mining seam is #3 coal seam. The average thickness and dip angle of the coal seam are 8.83 m and 7.5°, respectively. The Protodyakonov scale of hardness, Poisson’s ratio, and elastic modulus of the coal seam are 3.0, 0.22, and 6.72 GPa, respectively. The rock strata above the coal seam are, in order, sandy mudstone with 1.88 m thick and medium sandstone with 18.7 m thick. This study involves mining panel 3202, which is approximately 116 m in a strike. A fully mechanized longwall top coal caving method is used for #3 coal extraction in this working face, with a mechanized mining height of 3.0 m. The top-coal caving height is 5.83 m. 3202 panel is located at the east of the mining area 3203 panel with the width 120 m, where the coal has been extracted. The width of the strip coal pillar between the 3203 panel and the 3202 panel is 120 m, as shown in Figure 2.

4.2. Evaluation of Rockburst Risk

According to Section 3 of this paper, rockburst risk evaluation indexes of strip coal pillar are listed in Table 1.

According to A. H. Wilson progressive failure theory, the plastic zone width of strip coal pillar can be expressed aswhere is the plastic zone width of strip coal pillar, is the height of the strip coal pillar, and is the mining depth.

For the 3202 panel, the height of the strip coal pillar is 8.83 m, the mining depth is 1070 m, while the width of the strip coal pillar is 120 m. Therefore, the plastic zone width of the coal pillar calculated by equation (15) is 46.3 m, and the is 38.6%. According to the evaluation index of in Table 1, the coal pillar of the 3202 panel is of medium rockburst tendency.

Further, and calculated by equation (9) are 2.11 × 1010 and 1.018 × 1011 N/m, respectively, and the medium stiffness ratio is 0.207. According to the evaluation index of in Table 1, the coal pillar of the 3202 panel is of strong rockburst tendency.

According to the test results, the strain corresponding to the peak load of the coal seam under uniaxial compression is 0.013. So is 0.115 m. Further, the critical compression of the strip coal pillar calculated by equation (8) is 0.086 m. However, according to equation (12), is still unknown, which needs to be obtained by numerical simulation.

4.3. Numerical Simulation for the Elastic Deformation Index
4.3.1. Numerical Model Generation

The FALC3D was used to construct a numerical model. According to the geological conditions of the 3202 working face, the FALC3D was used to construct the numerical calculation model. Considering the influence of boundary effect, coal pillars of 30 m in width are reserved at the boundary along the direction of length and width, and the dimensions of the model are 222 m × 158 m × 69.5 m. The horizontal displacements of the four vertical planes of the model were restricted in the normal direction, and the vertical displacement at the base of the model was set to zero. The top surface of the model is allowed to move freely. Vertical stress of 26.875 MPa is applied to the top boundary of the model to simulate the self-weight of the overburden, and the horizontal side pressure coefficient is 0.5.

The model uses Mohr–Coulomb constitutive relations, and the mechanical properties applied in the model are listed in Table 2.

4.3.2. Evaluation of the Plastic Zone Width Ratio

When the width of the strip coal pillar is 120 m and the mining width of the 3202 panel varies from 76 m to 136 m, the plastic zone width ratio of the coal pillar is shown in Figure 3.

As shown in Figure 3, when the mining width of the working face is 76 m and 96 m, the plastic zone width of the strip coal pillar is 40 m and is 33.3%. When the mining width of the working face is 116 m and 136 m, the plastic zone width of the coal pillar is 44 m and 56 m, respectively, and is 36.7% and 46.7%, respectively.

The relationship between and mining width is shown in equation (16). When the width of the strip coal pillar is fixed, increases exponentially with the increase of mining width:

According to equation (16), when the mining width is less than 112 m, there is no rockburst tendency for the strip coal pillar; when the mining width is between 112 m and 128 m, there is medium rockburst tendency for the strip coal pillar; and when the mining width is greater than 128 m, there is strong rockburst tendency for the strip coal pillar.

When the mining width of the 3202 panel is 116 m and the width of the strip coal pillar varies from 60 m to 120 m, the plastic zone width ratio of the coal pillar is shown in Figure 4. When the width of the strip coal pillar is 60 m, the whole coal pillar has a plastic failure and is 100%. When the mining width is fixed, the plastic zone of strip coal pillar decreases with the increase of the width of the coal pillar. When the width of the strip coal pillar is 80 m, 100 m, and 120 m, the plastic zone width of the coal pillar is 60 m, 44 m, and 44 m, respectively, and is 75.0%, 44.0%, and 36.7%, respectively.

As shown in Figure 4, the relationship between and the width of the strip coal pillar is shown in equation (17). When the mining width of the 3202 panel is fixed, decreases exponentially with the increase of the coal pillar width:

According to equation (17), when the width of strip coal pillar is less than 108 m, there is strong rockburst tendency for the strip coal pillar; when the width of strip coal pillar is between 108 m and 120 m, there is medium rockburst tendency for the coal pillar; and when the width of strip coal pillar is greater than 120 m, there is no rockburst tendency for the coal pillar.

4.3.3. Evaluation of the Elastic Deformation Index

In the numerical simulation, the compression of the strip coal pillar is obtained by measuring the subsidence of the coal pillar roof. When the width of the strip coal pillar is 120 m, the compression of the coal pillar with different mining widths is shown in Figure 5.

It can be seen from Figure 5 that with the increase of mining width, the compression of the coal pillar keeps increasing gradually.

When the mining width is 76 m, 96 m, and 116 m, respectively, is 0.061 m, 0.066 m, and 0.072 m, respectively, and is 0.709, 0.767, and 0.837, respectively. When the mining width is 136 m, reaches 0.129 m, and increases to 1.5. As shown in Figure 6, the relationship between the and the mining width is shown in

According to equation (18), when the mining width is less than 76 m, there is no rockburst tendency for the coal pillar; when the mining width is between 76 m and 116 m, there is medium rockburst tendency for the coal pillar; when the mining width is between 116 m and 126 m, there is strong rockburst tendency for the coal pillar; and when the mining width is greater than 126 m, the rockburst has occurred, and increases rapidly which can be shown in Figure 6.

When the mining width of the 3202 panel is 116 m and the width of the strip coal pillar varies from 60 m to 120 m, the compression of the coal pillar is shown in Figure 7.

It can be seen from Figure 7 that with the decrease of the width of strip coal pillar, the compression of the coal pillar keeps increasing gradually. When the width of the strip coal pillar is 120 m and 100 m, respectively, is 0.072 m and 0.083 m, respectively, and is 0.837 and 0.965, respectively. When the width of the coal pillar is 80 m and 60 m, respectively, u0 is 0.128 m and 0.160 m, respectively, and is 1.49 and 1.86, respectively. As shown in Figure 8, the relationship between and the width of the coal pillar is shown in

According to equation (19), when the width of strip coal pillar is greater than 114 m, there is medium rockburst tendency for the coal pillar; when the width of strip coal pillar is between 98 m and 114 m, there is strong rockburst tendency for the coal pillar; and when the width of strip coal pillar is less than 98 m, the rockburst has occurred, and increases rapidly.

4.4. Comprehensive Evaluation Results

According to the above simulation results, the rockburst risk of strip coal pillar under different mining widths and the width of coal pillar conditions is shown in Figure 9. The results of rockburst risk evaluation in this paper are compared with those of the traditional rockburst tendency identification test in the laboratory, as shown in Table 3.

It can be seen from Figure 9 and Table 3 that under the same mining width and width of the coal pillar, the rockburst risk obtained by is higher than the evaluation result of . For the sake of safety, the evaluation result of is adopted. So, the strip coal pillar of the 3202 panel has a strong rockburst tendency.

As shown in Table 3, the result of the traditional rockburst tendency identification test of #3 coal in Gucheng Coal Mine is a weak rockburst tendency. Moreover, the comprehensive index method [34] is used to evaluate the rockburst risk of the 3202 panel, and the evaluation result shows that the 3202 panel is of medium rockburst tendency. However, the field monitoring results [35] showed that the deformation of the roadway on the side of the strip coal pillar was very serious, the maximum deformation was more than 1100 mm, large energy microseismic events gather relatively, and the rockburst risk of the coal pillar was high. It can be seen that the rockburst risk evaluation result in this paper is more in line with the actual situation. Therefore, the rockburst risk evaluation indexes of strip coal pillar proposed in this paper are more accurate and reasonable.

According to Figure 9, from the perspective of reducing the rockburst risk, for the strip mining of Gucheng Coal Mine, when the mining width of the working face is 116 m, the width of strip coal pillar cannot be less than 98 m. The width of strip coal pillar of the 3202 panel is 120 m, which meets the above requirements. Therefore, the rockburst risk of the coal pillar of the 3202 panel is relatively low.

5. Conclusions

(1)The traditional test method of coal and rockburst risk ignores the influence of coal pillar geometry size, and the evaluation results have little significance for the field guidance of deep strip mining. Considering the influence of mining width and coal pillar width in strip mining, three new indexes were put forward to evaluate the rockburst risk of deep strip mining, namely, the medium stiffness ratio () of the elastic and plastic zone in coal pillar, the plastic zone width ratio (), and the elastic deformation index () of the core zone.(2)The stiffness ratio index () of the elastic-plastic zone indicates the degree of strain-softening during the loading process of strip coal pillar. The plastic zone width ratio () of strip coal pillar can be measured by the drilling cutting method, which has strong field operability. And the elastic deformation index () indirectly reflects the relationship between the deformation energy accumulation of strip coal pillar and the critical state of rockburst. For these three indexes, the larger the index value is, the greater the rockburst risk is.(3)According to the new evaluation indexes of this paper, the strip coal pillar of the 3202 panel of Gucheng Coal Mine is in danger of strong rockburst. Compared with the results of the traditional rockburst tendency identification test and comprehensive index method, the evaluation result in this paper is more in line with the actual situation.(4)From the point of view of prevention and control of rockburst, the design principle of strip mining in Gucheng Coal Mine is put forward; that is, when the mining width of the working face is 116 m, the width of the strip coal pillar cannot be less than 98 m.

Data Availability

All data supporting this study are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This work was supported by the Shandong Provincial Natural Science Foundation, China (ZR2019BEE001, ZR2018MEE009) and Key Research and Development Projects of Shandong Province, China (2019GSF111038).