Abstract

In order to understand the instability characteristics of surrounding rock in the process of deep roadway excavation, a three-dimensional numerical model was established by FLAC3D to systematically analyze the influence of roadway surrounding rock stability under different in situ stress distribution forms, and the environmental coefficient of mining-induced stress η was defined, the larger the environmental coefficient of mining-induced stress is, the larger the surrounding rock stress environment is, and the range where the η coefficient is greater than 0.2 is called with the destruction-danger zone. When the initial vertical stress is maximum principal stress and minimum principal stress, by comparing the roadway along the middle ground stress direction and minimum or maximum in-situ stress direction, the variation characteristics of displacement, failure zone and failure hazard zone of roadway surrounding rock are obtained, which provides theoretical basis for the treatment of disaster accidents such as roadway surrounding rock instability and rock burst caused by deep high in-situ stress.

1. Introduction

All underground engineering activities must be carried out in the earth’s crust. Coal mining is to take the crustal rock mass and coal seam as the main research object, so it is extremely important to study the natural properties of underground rock mass and the in situ stress distribution [1, 2]. In 1912, the concept of geostress was first proposed by Heim, a Swiss geologist [3]. It is pointed out that the rock is in a state of approximate hydrostatic pressure due to the rheological action for a long time under the continuous action of large stress [46]. The stress is equal to the gravity of the overlying rock mass, i.e., . Under hydrostatic pressure, there is no shear stress, and any direction is the principal stress. In 1926, Soviet scholars proposed the following theory [7]: assuming that the rock mass is uniform and continuous elastic medium and according to the theory of elasticity, it is proposed that the vertical stress of the rock mass is equal to the pressure of the rock strata above it, i.e., . And the horizontal stress . For the gravitational field with the ground level, the lateral pressure coefficient , and when the hydrostatic pressure. However, the measured value may be greater than 1.0 or less than 1.0 depending on the area and position of the measured result of [810]. Literature [11, 12] showed that the anchor bolt not only provides supporting force on the roadway surface but essentially forms a load-bearing structure with certain strength through the coupling action of anchor bolt and surrounding rock. In literature [1315], the stress evolution, displacement field, local deformation, overall distribution, and destroy characteristics of surrounding rock anchoring structures were studied with different bolt spacing through model tests. In literature [16, 17], the influence of bolt preloading force and bolt spacing on the supporting strength of surrounding rock anchoring structure through numerical simulation software was analyzed, the supporting scheme of excavation roadway was designed, and the effectiveness of the supporting scheme by adopting the roadway displacement measurement scheme was verified. Literature [1820] showed that under the action of high stress, the two-side coal seam is soft, the stress is uneven, the range of loosening circle is large, and the surface deformation shows the characteristics of “two-side displacement amount greater than roof subsidence amount greater than floor heave amount.”

With the deepening of coal mining depth in recent years, the original rock stress of underground rock mass is also increasing, and the anisotropy is significant [21, 22]. As a result, the excavation roadway deformation speed is accelerated; roof caving, side deflection, and floor heave are increasingly serious; and the bursting liability is more intense. The number of roadway repair is also increasing correspondingly, and the repair rate of some soft rock roadway has reached more than 90% [2325]. Redistribution of in situ stress caused by roadway excavation is closely related to the original rock stress distribution state. In order to accurately grasp the influence of in situ stress distribution on the stability of roadway surrounding rock [2628], based on the discussion of the distribution form of in situ stress, this paper systematically analyzes the displacement curve of surrounding rock, the range of plastic zone, and the volume of destruction-danger zone to judge the influence on the stability of surrounding rock [29, 30].

2. The Engineering Background and Numerical Model

2.1. The Engineering Background

This paper takes 7703 bottom pumping roadway in Zouzhuang Mine, Huaibei, Anhui Province, as the engineering background. The construction horizon of the bottom pumping roadway is in sand and mudstone interbed and 16–30 m away from the coal seam. The lithology is mainly purple piebald mudstone and aluminum mudstone, followed by siltstone and fine sandstone. The occurrence of strata is relatively stable. The direct roof is mudstone, gray to dark gray, massive, with slippery surface, a large number of plant fossil fragments, flat fracture, brittle, fragile, thickness of 1.91∼4.13 m, and average of 3.15 m. The old top is fine siltstone, light gray, fine grain, massive, with thin siltstone bands, a few cracks filled by calcite veins, medium sorting, calcareous cementation, parallel bedding, dense, thickness of 4.10∼90.2 m, and average of 6.12 m. The direct bottom is mudstone, dark gray, argilaceous structure, local silty content is high, relatively dense, flat fracture, a large number of plant rhizome fossils, thickness of 1.02∼3.95 m, and average of 2.39 m. The old bottom is fine silty sandstone, gray, silty structure, the local with high fine sand content, banded, flat and dense fracture, a large number of plant rhizome fossils in the lower part, longitudinal fractures, parallel bedding, thickness of 2.90∼5.61 m, and average of 4.47 m. The comprehensive histogram of roof and floor within the area of working face 7703 is shown in Figure 1.

2.2. The Establishment of Numerical Model

Roadway surrounding rock stability of the numerical model is set up by using FLAC3D, as shown in Figure 2; the model size is 80 m ∗ 60 m ∗ 80 m, and the actual depth of the coal seam is 720-800 m. The numerical model includes 9 layers from top to bottom. The sandstone is 11 meters, mudstone is 9 meters, siltstone is 7 meters, conglomerate is 4 meters, mudstone is 5 meters, siltstone is 4 meters, mudstone is 23 meters, sandstone is 5 meters, mudstone is 12 meters, respectively, and the cumulative thickness is 80 meters. The dip angle of each rock stratum is 15°. Physical and mechanical parameters of rock strata are shown in Table 1.

Excavation roadway is located in the middle of the model 23 meters mudstone environment, occurrence depth of 760 m; roadway cross section adopts the design of arch section, the section size of 5 meters wide and 4 meters high, along the Y direction for excavation, excavation length of 60 m, model size greater than the roadway section size, thereby eliminating the boundary effect on the result of numerical calculation. The Mohr–Coulomb model is used to simulate the mechanical response characteristics of the surrounding rock mass in the process of excavation. A fixed displacement boundary condition is adopted at the bottom and around the model, and a vertical stress is applied at the top of the model to simulate other overburden strata that have not been built into the numerical model. The step-by-step characteristics of surrounding rock stress, displacement, plastic zone, and failure risk zone caused by roadway excavation are analyzed to evaluate the stability of surrounding rock under the influence of mining.

3. Influence of In Situ Stress Distribution on the Stability of Roadway Surrounding Rock

3.1. Influence of In Situ Stress Distribution on Roadway Surrounding Rock Displacement

In order to quantitatively analyze the influence of different forms of in situ stress distribution on the displacement of surrounding rock, a vertical survey line was arranged on the roof and floor of the roadway, respectively, and a horizontal survey line was arranged on the left and right sides of the roadway, respectively. According to the elastic mechanics knowledge, the excavation of the chamber in the undisturbed rock will inevitably cause the redistribution of the original rock stress, resulting in the production of plastic zone (and plastic zone). According to the analytical formula of the Gilch solution, the range of the plastic zone is generally 3–5 times of the length of the chamber rock mass. And the width of the section size in the numerical model is 5 m, so the range of the plastic zone caused by excavation is about 15–25 m. Therefore, the maximum length of 25 m is adopted here as the analysis area. The numerical model of the roadway surrounding rock cross survey line was established, as shown in Figure 3.

The vertical displacement of the survey line of the top and bottom was extracted and curves were drawn, as shown in Figure 4. When the initial vertical stress is the maximum in situ stress, the deformation degree of surrounding rock is higher when the roadway is tunneling along the middle in situ stress direction than when the roadway is tunneling along the minimum in situ stress direction. Under the condition of hydrostatic pressure, the surface deformation of surrounding rock of roadway is higher than the first two types of in situ stress distribution and deep roadway surrounding rock deformation is lower than the first two types of ground stress. For the above three kinds of in situ stress condition, vertical displacement shows slowly lowering trend from the roadway surface to the deep roadway direction, and surface reduced to 0 about 15 m away from the roadway. When the initial vertical stress is the minimum in situ stress, the degree of vertical displacement of surrounding rock on the roadway surface is higher than that in the direction of maximum in situ stress when the roadway is excavated along the middle in situ stress direction, and the degree of vertical displacement of surrounding rock on the roadway surface is smaller than that in the case of tunneling along the maximum in situ stress direction when the initial vertical stress is the minimum in situ stress. Compared with the case where the initial vertical stress is the maximum in situ stress, when the initial vertical stress is the minimum in situ stress, the vertical stress decreases faster from the roadway surface to the depth of the surrounding rock.

The horizontal displacement of the left and right sides of the survey line was extracted and curves were drawn, as shown in Figure 5. When the initial vertical stress is the maximum in situ stress, the horizontal deformation decreases rapidly from the roadway surface to the deep surrounding rock. Surrounding rock of roadway is excavated along minimum in situ stress direction of horizontal deformation degree which is higher than the tunneling of the intermediate principal stress, and hydrostatic pressure is bigger than the horizontal deformation of roadway under the conditions of the first two types of in situ stress distribution. It is important to note that the roadway is excavated along the middle ground stress direction, horizontal smaller range of compression zone. When the initial vertical stress is the minimum in situ stress, the decrease rate of horizontal deformation from the roadway surface to the deep surrounding rock is significantly less than that of the initial vertical stress which is the maximum in situ stress. The horizontal deformation of the tunnel is obviously higher than that of the tunnel with the maximum ground stress. According to the geological conditions, the rock strata have an inclination angle of 15°, which leads to the different gravity of the overlying rock strata at the positions of the two sides of the roadway after excavation. Therefore, the stress redistribution process is not symmetrical. However, according to the influence curve of ground stress on the horizontal displacement curve of surrounding rock in Figure 3, the difference between the left and right sides of the deformation displacement is small, so the support of the roadway can still adopt the symmetric support method for surrounding rock reinforcement.

3.2. Influence of In Situ Stress Distribution on Plastic Zone of Roadway Surrounding Rock

The influence of different in situ stress distribution forms on the plastic zone of roadway surrounding rock is shown in Figure 6. Under the condition that the initial vertical stress is the maximum in-situ stress, when the roadway is extracted along the middle in-situ stress direction, the failure zone of the surrounding rock is mainly concentrated on the two sides of the roadway and expands to four shoulder angles at the same time. The failure of surrounding rock only occurs on the roadway surface of the roof and floor, and the failure range is obviously smaller than that of the two sides of the roadway. The failure zone of the two sides extends to the interface between mudstone and upper siltstone and stops expanding. When the roadway was tunneling along the direction of the minimum in situ stress, the plastic zone of the surrounding rock mainly developed on the two sides of the roadway and no longer expanded to the four shoulder angles of the roadway. The plastic zone of the surrounding rock of the roof and floor was similar to that of the roadway tunneling along the direction of the intermediate principal stress but only developed in the surface rock. Compared with the roadway tunneling along the intermediate principal stress direction, the failure range of surrounding rock of the roadway is obviously reduced, and the spatial distribution of the roadway is contracted to the two sides of the roadway. Under the condition of hydrostatic pressure, the plastic zone of surrounding rock around the roadway basically presents uniform distribution. Compared with the previous two types of in situ stress distribution, the range of plastic zone of surrounding rock on both sides of the roadway decreases, while the range of plastic zone of surrounding rock on the roof and floor increases.

When the initial vertical stress is the minimum in situ stress and the roadway is excavated along the direction of maximum in situ stress, the failure range of surrounding rock of the roadway roof and floor is slightly larger than that of the roadway sides, but the failure area of surrounding rock still presents uniform distribution characteristics. When the roadway is extracted along the middle stress direction, the failure phenomenon of the two sides of the roadway only appears in the shallow surrounding rock. The failure zone is mainly concentrated in the roof and floor surrounding rock of the roadway, and the roof failure zone stops developing at the boundary of mudstone and upper siltstone. Compared with the roadway extracting along the direction of maximum in-situ stress, the failure area of surrounding rock on both sides of the roadway is significantly contracted, and the failure area of surrounding rock on top and bottom is expanded. Compared with hydrostatic pressure conditions, the latter two types of in situ stress lead to lower damage degree of surrounding rock on both sides of roadway and higher damage degree of surrounding rock on roof and floor.

The in situ stress distribution not only has a significant influence on the spatial failure modes of surrounding rock but also leads to significant differences in the range of plastic zones. In order to quantitatively analyze the failure range of surrounding rock, the total volume of failure units was calculated by accumulating the failure units, and a histogram was drawn, as shown in Figure 7. It can be seen from the figure that when the initial vertical stress is the maximum in situ stress, the plastic zone of surrounding rock when the roadway was tunneled along the middle in situ stress is significantly larger than that when the roadway was tunneled along the minimum in situ stress. When the initial vertical stress is the minimum in situ stress, the range of surrounding rock plastic zone when the roadway is driven along the middle in situ stress direction is slightly higher than that along the maximum in situ stress direction.

3.3. Influence of In Situ Stress Distribution on Roadway Surrounding Destruction-Danger Zone

The rock fracture which is within the scope of the roadway destruction area grows highly, and destruction area is the key of support and reinforcement. The neighboring area is a part of the destruction rock though it is within range of the elastic state, but its mining stress environment is very bad, it is easily affected by external disturbances and enters a state of destruction, and it will become a damage zone of surrounding rock of the region. When the bolt or anchor cable support is carried out, it is necessary to pay attention to that the anchor end cannot be completely in the destruction-danger zone. Otherwise, once the destruction-danger zone enters the failure state under the external disturbance, the overall failure of the supporting structure will be caused and the reliability of roadway support will be reduced. In order to characterize the advantages and disadvantages of the stress environment in the surrounding rock of the roadway, the mining-induced stress environmental coefficient is defined aswhere is the first invariant of the principal stress and J2 is the second invariant of deviatoric stress, and those can be expressed as follows:where , , and are the maximum, intermediate, and minimum principal stresses, respectively, and S1, S2, and S3 are maximum, intermediate, and minimum deviatoric stresses, respectively, , , and (i = 1, 2, 3) are hydrostatic stress, deviatoric stress, and principal stress, respectively, and and , .

Since the first invariant of principal stress I1 represents the compression degree of the surrounding rock in three directions, the larger its value is, the better the stress environment the surrounding rock is. The second invariant of deviatoric stress, J2, represents the shear degree of surrounding rock. And the larger the value is, the worse the stress environment of surrounding rock is. And the mining stress environment factor is defined the ratio of the second invariant deviatoric stress to first principal stress invariants. Due to rock compressive shear, the greater the value of mining stress environment factor, the easier the rock mass is to be destroyed. And the rock shear strength is about 10%–40% of rock compressive strength, so the scope of mining stress environment factor greater than 0.2 is defined as a destruction-danger zone. The spatial distribution morphology of surrounding rock in the destruction-danger zone under different in situ stress conditions was extracted, as shown in Figure 8. When the initial vertical stress is the maximum in-ground stress, the destruction-danger zone around the roadway takes on a butterfly shape, and its range is obviously larger than that of the plastic zone. When the roadway is excavated along the direction of the middle in-ground stress, the destruction-danger zone has the largest range, and when the roadway is excavated along the direction of the middle in-ground stress, the destruction-danger zone presents a tendency of contraction. Under the condition of hydrostatic pressure, the destruction-danger zone around the roadway is basically uniformly distributed, and its range is slightly larger than that of the plastic zone. When the initial vertical stress is the minimum in situ stress and the roadway is excavated along the direction of the maximum in situ stress, the failure danger area is basically uniform in the two sides of the roadway and the surrounding rock of the floor, and the scope of the failure danger area in the floor increases obviously. When the roadway was excavated along the direction of the middle ground stress, the destruction-danger zone of the surrounding rock on both sides of the roadway contracted obviously, and the destruction-danger zone of the surrounding rock on the roof and floor expanded obviously, and the spatial distribution of the destruction-danger zone was like a wine glass. Compared with the hydrostatic pressure condition, the latter two types of in situ stress lead to the increase in the uneven distribution of the failure risk zones around the roadway and the expansion of the development range.

The volume of all units within the scope of the destruction-danger zone was accumulated, and the column diagram of the destruction-danger zone range of the surrounding rock of the roadway under different in-ground stress distribution was obtained, as shown in Figure 9. Under the condition of hydrostatic pressure, the stress environment of surrounding rock is relatively good, and the range of failure danger of surrounding rock is obviously smaller than the other four types of in situ stress distribution. When the initial vertical stress is the maximum in situ stress, the range of destruction-danger zone is higher when the roadway is tunneling along the middle in situ stress direction than when the roadway is tunneling along the minimum principal stress direction. When the initial vertical stress is the minimum in situ stress, the range of surrounding rock destruction-danger zone is higher when the roadway is tunneling along the middle in situ stress direction than when the roadway is tunneling along the maximum in situ stress direction.

4. Conclusion

(1)When the initial vertical stress is the maximum in situ stress, the roof subsidence and floor bulge of the roadway are significantly greater than those when the initial vertical stress is the minimum in situ stress, but the displacement of the two sides of the roadway is significantly higher than that when the initial vertical stress is the minimum in situ stress.(2)When the initial vertical stress is the maximum in situ stress, the range of plastic zone of surrounding rock is significantly higher than that of the minimum principal stress when the roadway is excavated along the direction of intermediate principal stress. When the initial vertical stress is the minimum in situ stress and the roadway advances along the direction of the maximum in situ stress, the scope of the plastic zone around the roadway is slightly smaller than that along the direction of the middle in situ stress.(3)In order to characterize the advantages and disadvantages of the stress environment in the surrounding rock of the roadway, the environmental coefficient of mining-induced stress is defined. The larger the environmental coefficient of mining-induced stress is, the larger the surrounding rock stress environment is. The range of this coefficient greater than 0.2 is defined as the destruction-danger zone.

Data Availability

The data 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 financially supported by National Natural Science Foundation of China (Grant no. 51934008) and Natural Science Foundation of Heilongjiang (Grant no. LH2021E104).