Abstract

In order to reveal the failure mechanism of the deep roadway under mining-induced pressure in coal mine, the boundary equations for the plastic zone around the deep roadway were deduced, and then the evolution laws for morphology of the plastic zone and the relationship between the morphological indexes and the stability of surrounding rock were discussed. The results show that, for the deep roadway, the effect of mining on the plastic zone is more sensitive than that on the shallow one. Even if the changes of mining influence are small, they may also cause extremely serious plastic failure of surrounding rock masses, leading to the sudden instability of the roadway. When the plastic wings of the plastic zone are approximately perpendicular to the roof, floor, or sidewall, the large deformation and failure of the deep roadway are very likely to occur. Compared with the index of the uniformity coefficient, the irregular shape coefficient can be used to better characterize the differences in the plastic zone morphology. Finally, a case study was provided to apply the principles for the formation and extension of a butterfly-shaped plastic zone.

1. Introduction

With the increasing shortage of coal resources in shallow depth, the mining depth of coal mine in China is increasing at a speed of 10–25 m per year. Currently, China has 47 coal mines more than 1000 m deep [13]. Deep coal mining at 1000 m in other countries like Poland, Germany, Britain, and Japan is also common [4]. Compared with shallow mining, deep coal mining is confronted with some complicated problems such as large deformation, rockbursts, and water hazards [58]. Among them, the large deformation and failure of deep roadways threaten seriously mining production and safety [9].

In order to minimize the disasters, the first thing to be solved is the clarity of the failure mechanism for deep roadways or deep rock. In recent years, scholars have carried out a lot of research on it by using numerical simulation, model testing, theoretic analysis, etc. For example, Dolezalova et al. [1016] used FEM, FLAC3D, FDEM, and 3EDC to simulate stress distribution and failure evolution of deep engineering, respectively. In view of the nonlinear problems of large deformation occurring in a deep soft rock roadway, He [17] adopted the material point method (MPM) to simulate the large deformation and failure process of deep rocks. Geomechanical model test is also an important way to investigate the failure behavior of deep surrounding rock. Sun et al. [1820] conducted large-scale model tests. For theoretical analysis, Wang et al. [21] developed a dynamic failure constitutive model to computer deep large deformation, and Zareifard and Fahimifar [22] deduced an analytical solution for the stresses and deformations of deep tunnels considering the damaged zone. By field detection, the failure characteristics of the deep rock were also obtained [23, 24]. Zhao et al. [2528] performed rock creep tests under multilevel load and revealed the nonlinear relationships between the instantaneous strain, viscoelastic strain, viscoplastic strain, and high deviatoric stress.

Large deformation and failure of deep roadways are not only related to high stress and weak rock masses but also closely associated with the mining-induced pressure (the side and front abutment pressure); thus, the influence of mining is sometimes a nonignorable factor [29, 30], whereas limited research studies [31, 32] have been carried out on the damage and failure of deep roadways subjected to mining-induced pressure.

Although the deformation and instability of deep roadways have been studied, the failure mechanism is still far from being complete and satisfactory. The previous research [33, 34] shows a close relation between the range and shape of the plastic zone and the stability of the roadway. Moreover, the plastic zone will be extended under the external extra load, which would result in the increase of unrestrained deformations and deterioration of the surrounding rock strength [3540]. Therefore, the distribution and evolution of the plastic zone around the deep roadway under the secondary stress field is key to investigate its failure mechanism. In this work, the morphological evolution of the plastic zone and the relationship between the morphological indexes proposed and the stability of surrounding rock were explored.

2. Failure Characteristics of Deep Roadway under Mining-Induced Pressure in Coal Mine

For deep roadway, the high geostress and the intense mining activities nearby mainly lead to a quite nonuniform, complex stress field, which causes severe failure of the mining-influenced deep roadway. The main failure characteristics of this kind of roadway are as follows:(1)The roadway roof subsides sharply, which can reach 500∼1200 mm, as shown in Figure 1(a). Part of rock masses in the top anchorage region slides down.(2)The deformation of sidewalls is serious. Extruding magnitude of a single sidewall ranges from 400 mm to 800 mm. Besides, the phenomenon of sidewall caving appears in some regions (Figure 1(b)).(3)The obvious floor heave causes cracking of shotcrete lining and serious inclination of the brace and track; thus, the roadway cannot be used properly (Figure 1(c)).(4)Influenced by high geostress, repeated mining activities, and geological structures, the deformation of the roadway surrounding rock appears significantly asymmetry, which hinders the utilization of the roadway (Figure 1(d)).

3. Evolutionary Laws of Morphology for Plastic Zone around Deep Roadways under Mining-Induced Pressure

3.1. Evolution Equation for Boundary of Mining-Induced Plastic Zone
3.1.1. Basic Assumptions and Mechanical Model

The deep gate roadways are often subjected to mining influence from their own working face and adjacent excavation, which contributes to the superposition of the initial stress field, abutment pressure, and other dynamic loads. Therefore, the regional stress field of deep gate roadways is quite complex. In order to solve problems conveniently, the following assumptions are made:(1)The length of the roadway along the horizontal direction is infinite, with a circular cross section, the radius , and the buried depth H ≥ 20.(2)The surrounding rock is considered the isotropic, homogeneous medium and an incompressible material in plasticity.(3)The initial maximum and minimum principal stresses of the roadway are and , respectively, which are parallel to the coordinate axes, without consideration of the supporting force.(4)Two parameters and (i.e., mining coefficients for the maximum and minimum principal stresses) that characterize the effect of mining on the initial principal stress are introduced. The principal stress field of the roadway changes with mining activities except for the stress direction. The mechanical model of the roadway is shown in Figure 2, and the surrounding rock can be divided into the plastic zone (radius equal to ) and elastic zone from inside to outside.

3.1.2. Boundary Equation of Plastic Zone (the Principle Stresses Parallel to Coordinate Axes)

The pressure (in Figure 2) acting on the surrounding rock is decomposed into two cases I and II, as shown in Figure 3.

For Case I, according to the theory of elasticity [17, 21], the elastic stress field of the roadway surrounding rock subjected to uniform pressure is as follows:where and are the radial stress and tangential stress at any point in surrounding rock under polar coordinates, respectively, and r represents the radial coordinate of this point.

For Case II, the stress field of the surrounding rock is as follows:where is the shear stress in polar coordinates and is the tangential coordinate of the point.

After superposition of the stress fields of surrounding rock in two different cases, the elastic stress field of the circular roadway subjected to mining-induced pressure can be obtained as follows:

In terms of the literatures [21, 22], the polar expression of the Mohr–Coulomb yield criterion iswhere denotes the cohesion of the surrounding rock and indicates the internal friction angle.

Assuming that , , , and , substituting (3) into (4), the implicit function referred to and is obtained as follows [41]:

When , it becomes the interface between the elastic and plastic zones of the deep roadway, i.e., the equation for boundary of the plastic zone. Especially for and , the radius of the plastic zone around the circular roadway under uniform stress field can be obtained from formula (5), which is

3.1.3. Boundary Equation of Plastic Zone (the Principle Stresses with a Certain Rotation Angle)

The engineering practice shows that not only the magnitude of principal stress of the roadway but also the direction of that will be changed due to the nearby coal mining. It is assumed that the deflection angle of the maximum principal stress is relative to Figure 2, and the counterclockwise direction is positive. A new coordinate system is established, as shown in Figure 4.

In the new and old coordinates, the polar coordinates of any point A in the surrounding rock are and , respectively. It is seen from Figure 4 that the rotation angle of the new coordinate axes relative to the old ones is .

In the new coordinate system, the boundary equation of the plastic zone around the deep roadway under mining-induced pressure is similar to formula (5), given as follows:where , , , and .

Note that the relationship between the new and old polar coordinates of any point A is and ; therefore, the implicit equation for the boundary of the plastic zone under deflection of the principal stress field can be obtained from formula (7), given as follows:where and .

3.2. Morphological Evolution of Plastic Zone around Roadway under Mining-Induced Pressure

Based on formula (5), the evolution of the plastic zone and the failure of surrounding rock of the deep roadway under mining-induced pressure were discussed. The diagrams of plastic zone distribution were drawn (Figure 5).

In Figure 5, the parameters are set as follows:  = 2 m,  = 3 MPa,  = 25°, and  = 20 MPa. It is shown in Figure 5 that the pressure caused by coal mining has significant effect on the distribution of the plastic zone around the deep roadway. The shape of the plastic zone varies with the increase of the pressure difference between horizontal and vertical direction, from circle to ellipse, rounded rectangle, and butterfly shape. Except for the butterfly shape, the dimension of the other three plastic zones is generally between 2.6 m and 3.2 m, without significant differences. In contrast, the plastic radius of the butterfly-shaped plastic zone increases greatly in the inclined direction between the maximum and minimum principal stress due to the formation of the plastic wing, and the maximum plastic radius occurs nearby the bisector of the two directions.

Due to the influence of mining, the butterfly-shaped plastic zone around the deep roadway is more likely to occur. It is observed from Figure 5(a) that, with the increase of the vertical mining coefficient for the maximum principal stress, the plastic wings form and experience a rapid expansion, which present the slender type. For example, the size of the plastic wing grows exponentially during changing from 2.5 to 2.8, with the maximum plastic radius () of 15 m. This demonstrates that the plastic damage of rock masses within this area is very sensitive to the mining coefficient. Even if the changes of mining influence are small, it may also cause serious localized deformation and failure of the surrounding rock mass. The phenomenon that the shear deformation and failure often occur at the top or bottom corner of the deep roadway under the mining-induced pressure is proved theoretically by the formation and nonlinear extension of the plastic wings, which also shows that the deformation of the mining-influenced roadway is larger than the ordinary one. Thus, it is quite difficult to ensure the stability of the roadway.

If the direction of the principal stress field for the mining-influenced roadway rotates, how the morphology of the plastic zone distributes and what is the relationship between the stability of surrounding rock and the morphological characteristics of the plastic zone should be investigated. The distribution diagram of the plastic-zone boundary around the circular roadway is plotted by formula (8), as shown in Figure 6. The parameters are set as follows:  = 2 m,  = 3 MPa,  = 25°,  = 1,  = 2.8, and  = 20 MPa.

It is found from Figure 6(a) that when the maximum principal stress is perpendicular to the roof stratum of the roadway (i.e., ), the maximum depths of plastic damage mainly exist at the top and bottom corner, while the range of the plastic zone in the roof is quite small. With the increase of the deflection angle β of the maximum principal stress, the plastic damage of the roadway roof is gradually enlarged, and consequently, the stability of the roof is also deteriorated, as shown in Figure 6(b).

In Figure 6(c), the maximum principal stress is deflected at , while the four plastic wings are just perpendicular to the sidewalls, floor, and roof. Under the ultrahigh stress field ( = 42 MPa) caused by strong mining, the depth of the plastic wings can reach up to 14 m; therefore, floor heave, roof sinking, and side wall extrusion often occur in the deep mining-affected roadway. It should be pointed out that if the rock mass structure of the roof is poor or the support is improper, the caving risk of the roof rock mass within the plastic wing will be very high. The plastic wings rotate with the continuous rotation of the maximum principal stress, and the area of the plastic zone in the roof decreases gradually, as shown in Figure 6(d).

3.3. Malignant Expansion of the Plastic Zone around a Roadway

Figure 7 shows the effect of on the plastic radius r, and the necessary parameters are set as follows:  = 2 m,  = 3 MPa,  = 25°, and  = 1. The values of r are all taken from the intersection of the plastic wing in the first quartile and the direction of to represent the maximum plastic radius . It is seen from Figure 7 that whether, for low stress or for high stress, the curves all can be divided into two stages. During the first stage, the plastic radius keeps increasing near-linearly as the mining coefficient for the maximum principal stress raises, and r can be stabilized at a certain value. When Dζ1 exceeds the critical value Dζc, the second stage starts. Clearly, the growth rate of r increases rapidly and r approaches to infinity finally, which shows that the malignant expansion of the plastic zone occurred. Here, Dζc is called the critical mining coefficient. The greater the initial stress, the smaller the value of Dζc, which indicates that the mining-influenced roadway is more likely to produce malignant expansion of the plastic zone under high stress.

Figure 8 shows the relationship between the radius of the plastic zone and the mining coefficients for the maximum and minimum principal stresses under nonuniform high stress, where  = 20 MPa. When Dζ3 takes a larger value (3∼5), the plastic radius r decreases from the infinite value to the minimum value (at the neutral line) with the rise of Dζ1, and then it increases to a relatively small value. The neutral line indicates the congruent relationship between Dζ1 and Dζ3 when r takes the minimum value. It means that the radius of the plastic zone around the roadway reaches the minimum when the pressure distributes uniformly. The existence of the neutral line proves that one of the principles for the design and support of the deep roadway under mining-induced pressure is that the stress environment of surrounding rock should be improved, and its nonuniformity is required to be minimized for reducing the depth of plastic failure.

Based on the above analysis, the following definitions are made:(1)The malignant expansion of the plastic zone (MEPZ): the radius r of the plastic zone increases rapidly with the increase of a certain variable ξ; finally, r cannot be stable at a certain magnitude but grows infinitely. This can be expressed as , where is the function of the plastic radius and denotes the upper limit of the variable .(2)The criticality for the malignant expansion of the plastic zone (): the plastic radius r increases gradually with the growth of the variable ; when increases to a certain value , the growth rate of r increases rapidly, and eventually r is approaching infinity. Here, is called criticality for the malignant expansion of the plastic zone.

4. Index System for Morphological Characteristics of Plastic Zone around a Circular Roadway

In this section, the index system for characterizing the morphological features of the plastic zone around the circular roadway was established using the geometric morphology theory.

4.1. Main Indexes for Morphological Characteristics of the Plastic Zone

Generally, the elastic or plastic solution for the underground roadway is treated as a plane problem, and the plastic failure field obtained is also a 2D figure. Therefore, the indexes established and related parameters are all for the 2D diagrams.

4.1.1. Uniformity Coefficient ()

The uniformity coefficient is used to characterize the differences in plastic failure depths on the plastic zone boundary around the circular roadway along any direction from 0 to 2π. The closer the value of is to 1, the smaller the differences in plastic radius of the plastic zone in all directions are, the closer the shape of the plastic zone is to circle. Generally, the roadway with the circular plastic zone is more stable and has better integrity than those with the heteromorphic plastic zone with the same extension area. can be expressed as follows:where is the perimeter of the plastic zone boundary around the roadway and denotes the area enclosed by the plastic zone boundary.

4.1.2. Extension Factor ()

The extension factor is defined as the ratio of the area of the plastic zone of the surrounding rock to sectional area of the roadway. The extension factor reflects the relative size of the plastic zone extension. The larger the is, the larger the range of plastic failure of the surrounding rock is, i.e., the more unstable the surrounding rock mass is. The index of is not only related to the size of the plastic zone but also associated with the dimension of the tunnel. Hence, it is applicable to comparatively analyze the plastic zones around roadways with different sizes. can take the following form:where is the radius of the roadway section; is the sectional area of the roadway; and represents the area of the plastic zone, and it can be obtained from .

4.1.3. Equivalent Radius ()

The equivalent size and range of the plastic zone can be described by the equivalent radius that is one of the important geometric parameters for intuitively reflecting the macroscopic geometry of the plastic zone. The plastic zone with irregular shape can be transformed into an equivalent circular one in terms of the principle of the same area, shown as the dashed line in Figure 9. The equivalent radius is calculated by the following formula:

4.1.4. Extended Depth Coefficient ()

In order to characterize the relative size of the extended depth for plastic limit, the extended depth coefficient is introduced. It can be written as follows:where is the maximum plastic radius of the roadway. The greater the is, the more serious the plastic failure of surrounding rock is at the direction in which is located. Although the different roadways have the same values, the values are not equal. Thus, the local stability of the surrounding rock is also different.

4.1.5. Irregular Shape Coefficient ()

The shape of the plastic zone is closely related to the stability of the roadway. Usually, the roadways with an irregular plastic zone are less stable than those with the regular one due to weak regions in surrounding rock. Then, the index is used to characterize the degree of irregularity for the shape of the plastic zone. It is defined as the ratio of the maximum plastic radius to the minimum one, given as follows:where is the minimum plastic radius of the roadway surrounding rock. The range of the irregular shape coefficient is ; the closer the value of is to 1, the more regular the shape of the plastic zone in surrounding rock masses is.

4.2. Analysis for Indexes of Morphological Characteristics of the Plastic Zone

In order to clarify the change of indexes with the evolution of the plastic zone, the plastic zones (Figure 10) with different shapes were analyzed. The conditions for calculation are the same as given in Section 3.2. The relationship between the indexes and is shown in Figures 1113. is the ratio of the vertical principal stress to the horizontal one, i.e., .

It can be seen from Figures 10 and 11 that and increase continuously with the increase of . When exceeds 1.9, the growth rate of and becomes larger, which indicates that the differences in the plastic extension depth at various directions are accelerating after the plastic zone evolved into the butterfly shape. Thus, the shape of the plastic zone becomes more irregular, and the stability of the surrounding rock for the roadway is reduced greatly. Additionally, even if the change laws of and are similar, is more sensitive to the morphological changes than , which shows that the differences in morphological changes can be better reflected by the index .

The change process of and with can be all divided into two stages as shown in Figure 12. The first stage is a slow growth period, and then the second stage of rapid growth starts, which shows that the area of the plastic zone and the limit of extended plastic depth will accelerate when the shape of the plastic zone changes into the butterfly shape. For , the growth rate of the area of the plastic zone is much larger than that of the limit of extended plastic depth, and the overall stability of the surrounding rock will be decreased greatly.

In Figure 13, the equivalent radius increases approximately at a constant acceleration during the initial stage (); after that, the accelerating growth stage begins, in which the maximum equivalent radius of the plastic zone can be achieved (8.805 m). It is far beyond the anchorage range of the rock bolt, which indicates that the size of the plastic zone around the deep roadway under the high deviatoric stress becomes quite large and that the size change of the plastic zone is more sensitive when the plastic zone enters into the stage of butterfly-shaped evolution; therefore, the roadway tends to more easily fail.

5. Field Validation

5.1. Project Overview and Failure Status of the Roadway

The #603 tailgate is located at the west second raise mining area of a coal mine in Jiangxi Province, China. B4 coal seam with the average thickness of 2.8 m is the main mining object. The design length of #603 tailgate is 787 m with a buried depth of 806 m. The geological structure of coal-measured strata is relatively simple. The immediate roof of this roadway is siltstone with a thickness of 8∼10 m, and the main roof with a thickness of 2∼4 m is the interbed of siltstone and fine sandstone. While the immediate floor is mudstone with a thickness of 2∼4 m, the main floor of sandstone is 12 m in thickness.

For the deep coal mining, it is very difficult to maintain the roadways, resulting in unbalance of the mining-drifting. The gate roadways have to be tunneled along the unstable gob before finishing the #602 working face. Therefore, the #603 tailgate is not only influenced by the unstable gob but also by the present working face mining. The failure characteristics of #603 tailgate were obtained by field investigation and testing:(1)Serious subsidence of the roof stratum and the extrusion of the sidewall occurred. As shown in Figure 14, during the recovery of #603 working face, the roof had an extremely serious subsidence, leading to large angle lean of the metal pillar, failure of bolts and cables, and localized roof fall. Both sidewalls squeezed out severely, the deformation of a single sidewall was up to 300∼700 mm, and the contraction rate of the roadway section was nearly 40%.(2)Serious floor heave occurred within the 60% length of the roadway, which hindered the use of roadway and the normal production of the coal mine.(3)The distribution characteristics of the failure field in #603 tailgate surrounding rock were detected in boreholes by the recorder for strata detection, as shown in Figure 15. The results show that the damage depth of surrounding rock reaches nearly 8 m and the damage range is abnormally large; thus, the strength and integrity of internal rock masses around the roadway deteriorate. The failure depth of the roof and the sidewalls is significantly larger than that of the two top corners, while the failure depth of the roof is larger than that of the sidewalls.

5.2. Analysis for Failure Mechanism of Roadway Based on the Mining-Induced Plastic Zone

According to the field situations, the parameters for calculation are simplified as follows:  = 2.3 m,  = 3.3 MPa,  = 27°,  = 20.2 MPa, and  = 16 MPa. Additionally, ,  = 1.0, and  = 40° in terms of the past monitoring data for ground pressure. Then, the distribution diagram of the plastic zone around the roadway can be drawn by formula (8), as shown in Figure 16, and the magnitudes of morphological indexes are shown in Table 1.

It can be seen from Figure 16 that the plastic zone of the surrounding rock shows an obvious butterfly shape and the plastic wings extend greatly along the roof, floor, and sidewalls of the roadway. The indexes , , , and in Table 1 indicate that the shape of the plastic zone is quite irregular and that the area and dimension of plastic failure in the surrounding rock are very large, especially for the limit of extended depth for the plastic wings () which is 5.76 times the radius of the roadway so that it is extremely serious for the plastic deformation and failure of rock masses in these plastic wings grown along the roof, floor, and sidewalls.

Based on the above analysis and the field situation, the failure mechanism of #603 tailgate is deemed that the soft surrounding rock produced notable deformation in the early stage due to insufficient strength and stiffness of support under high geostress. With the violent influence of roof caving in the gob and mining for #603 working face, the stress field around the roadway became the ultrahigh, nonuniform stress field; furthermore, the direction of the principal stress deflected, and consequently, a typical butterfly-shaped plastic zone formed (Figure 16). Since the plastic wings with great extended depths were nearly perpendicular to the roof, floor, or sidewalls and the fissures were more developed in surrounding rock, all of them caused nonlinear failure of the surrounding rock of the roadway.

6. Conclusions

(1)The analysis of the boundary equation for the plastic zone shows that when considering the influence of mining, the butterfly-shaped plastic zone in surrounding rock is more likely to occur, and the size of the plastic wing and the range of plastic failure are greater than those without consideration. If the plastic wings are approximately perpendicular to the roof, floor, or sidewall, large deformation and failure are very likely to occur for the deep roadway.(2)The definitions of malignant expansion of the plastic zone and its criticality were made. The effect of mining on the plastic zone around the deep roadway is more sensitive than that on the shallow one. Even if the changes of mining influence are small in short time, it may also cause serious deformation and plastic failure of surrounding rock masses, leading to the sudden instability of the roadway.(3)The index system for morphological characteristics of the plastic zone around the circular tunnel was built, and the relationship between these indexes and the stability of surrounding rock were deeply analyzed. Compared with , is more sensitive to the morphological changes; thus, the differences in morphological changes can be better reflected by the index .(4)The failure characteristics of #603 tailgate were obtained from the field detection and geological conditions. After that, the mechanism of large deformation and failure for #603 tailgate under deep mining-induced pressure was revealed based on the principles for the formation and extension of the butterfly-shaped plastic zone around the mining-influenced roadway.

Data Availability

The data (Table 1 and Figures 1 to 16) used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Natural Science Foundation of China (Nos. 51504091, 51434006, 51774131, and 51774130) and the Natural Science Foundation of Hunan Province, China (No. 2018JJ3166).