*Article* **Study on the Influence and Control of Stress Direction Deflection and Partial-Stress Boosting of Main Roadways Surrounding Rock and under the Influence of Multi-Seam Mining**

**Dongdong Chen, Fangfang Guo, Zijian Li, Xiang Ma, Shengrong Xie \*, Yiyi Wu and Zhiqiang Wang**

School of Energy and Mining Engineering, China University of Mining & Technology-Beijing, Beijing 100083, China

**\*** Correspondence: xsrxcq@cumtb.edu.cn

**Abstract:** A large coal pillar (usually more than 90 m) is generally left in place to ensure the stability of main roadway groups, due to its long service lifespan, which commonly also causes a significant loss of coal resources. The design of the width of the protective coal pillar and the control system of the surrounding rock are directly determined by the characteristics of the stress field and the damage mechanism under the influence of the mining activities. However, there are few studies on the effects of the partial-stress boosting and the direction deflection of the stress field on the failure evolution of the surrounding rock (especially in multi-seam mining). In this paper, theoretical analysis and numerical simulation are used to investigate the direction evolution of the maximum principal stress in front of the working face with malposition distances between the upper and lower working faces during the influence of double coal seams mining. Furthermore, a large-scale numerical model is used to study the deviatoric stress evolution of the surrounding rock and the propagation process of the plastic zone in the main roadway group with different widths of protective coal pillars. Then, an asymmetric cooperative anchoring classification method is proposed to strengthen the roadway support, depending on the critical area of the deviatoric stress in the roadway surrounding rock. The peak zone deflection of the deviatoric stress determines the evolution direction of the plastic area, and the peak value of the deviatoric stress presents a typical asymmetric stress boosting on both sides of the roadway. These findings are validated by the on-site ground pressure monitoring results and the practical failure modes of the surrounding rock.

**Keywords:** main roadway group; mining-induced stress; surrounding rock stability; ground control; deflection of principal stress

#### **1. Introduction**

Energy and mineral resources serve as the economic development drivers for countries around the world [1]. China's proven coal resource reserves rank second in the world, and the total consumption of coal in China's primary energy consumption in 2020 still accounts for more than 50% of the proportion [2,3]. Ensuring the safe production of coal means a stable economic development. Long-wall underground mining arose in Europe and is widely used in China because of its high productivity. About 85 percent of China's coal production comes from underground, and long-wall mining is used for most of it, about 20 percent of the total; and about 5 percent of India's production comes from underground [4–7]. Although long-wall mining has a high productivity, because of the characteristics of roadway layout, large coal pillars should be set in terminal stage of long-wall face to ensure the safety of the roadway group [8]. Due to the long service life of the main roadway group, the mine tends to set up excess coal pillars to maintain the main roadway safety, so the study on the width of the roadway protective coal pillar is important.

The main roadways are the lifeblood of mine transportation, ventilation, pedestrians, etc., and they generally need to serve the whole production process of the mine. Therefore,

**Citation:** Chen, D.; Guo, F.; Li, Z.; Ma, X.; Xie, S.; Wu, Y.; Wang, Z. Study on the Influence and Control of Stress Direction Deflection and Partial-Stress Boosting of Main Roadways Surrounding Rock and under the Influence of Multi-Seam Mining. *Energies* **2022**, *15*, 8257. https://doi.org/10.3390/ en15218257

Academic Editor: Nikolaos Koukouzas

Received: 10 October 2022 Accepted: 2 November 2022 Published: 4 November 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

most of the main roadways are built on the roof or floor of the coal seam in underground mining. Furthermore, the coal seam roadways are easy to excavate, and coal resources can be recovered as soon as possible, reducing the financial pressure on mine production. In recent years, the main roadways are increasingly inclined to be placed in coal seams in some mine areas with relatively good geological conditions, such as in Yulin and Erdos, China.

For the roadway problems of a coal mine, scholars around the world have made many efforts to study the properties of the roadway surrounding rock, the state of the deposit, the influence of mining stress, the mechanism of the surrounding rock deformation, and the role of support, etc., and have accumulated rich theoretical and practical engineering experience. Gao et al. [9–11] conducted a series of numerical simulations to study the roadway's squeezing failure and shear failure, using the UDEC Trigon method and analyzed the effects of rock bolts on the roadway support. Wang et al. [12] tested the properties of surrounding rock considering the influence of the lower coal seam on the stress state by a numerical simulation. The creep characteristics of the surrounding rock of the mining roadway and its effects on the deformation were systematically performed. Pan et al. [13] focused on the deformation and failure process and the mechanism of the surrounding rock induced by the stress adjustment in a deeply arched roadway. A self-developed large-scale simulation test was carried out. Mohammadi et al. [14] developed a geometric computational model to calculate the extension of the excavation-damaged zone above the gate roadways, induced by longwall mining and proposed a new mathematical formulation to determine the face influence coefficient.

Huang et al. [15] clarified the connotation of roadway mining and large deformation under the consideration of the deviatoric stress. They proposed a theoretical framework of a large deformation of surrounding rock rheology and structural instability in deep mining roadways. Y et al. [16] investigated the problems of roadway stability with repeated mining-induced disturbance and focused on the direction and emphasis of the stability control of roadway surrounding rock. Shan et al. [17] summarized the current research status of coal tunnel support technology at home and abroad, in recent years.

Moreover, they proposed that roadway support adopted a more diversified integrated active support technology and showed the trend of gradually moving into intelligent support. Some researchers [18–24] concentrated on the significant deformation of the surrounding rock in deep roadways and targeted control technologies were proposed to solve this problem. These papers are dedicated to discovering the failure mechanism of the roadway or efficient roadway support technology to prevent the instability of the roadway. The Coal pillar width does contribute an essential part to the stability of the roadway. The non-pillar technologies (gob-side entry driving or gob-side entry retaining) have been widely used in China, due to the high coal resource recovery rate [25–32]. Currently, for a coal seam main roadway group in China, the protective coal pillar is usually larger than 90 m for the stability of the roadway, which causes a tremendous waste of coal resources. Some experts have also studied the width of coal pillars to protect the main roadways. For example, some scholars [33,34] studied the stability control of the main deep roadway groups, considering the characteristic rheology coupled with multiple disturbances of mining activities of the adjacent working faces. It determined that it was rational to leave a protective coal pillar for main roadway groups with a width of 100 m.

However, at present, for the close multi-seam mining, when the width of the endmining coal pillar of each coal seam is different, the stress superposition and direction deflection law in the area in front of the working face, caused by the influence of mining and the evolution law of the damage to the main roadway surrounding rock, is rarely studied. Therefore, it is impossible to determine the location of the high deviatoric stress area and the critical control orientation of the main roadway surrounding rock. This also leads to the roadway surrounding rock support design that is not a targeted basis.

The study background of this paper is to examine the specific engineering geological conditions of close two-seam coal mining and the different widths of end-mining coal pillars. On this basis, the following studies are conducted: the law of the principal stress direction deflection and plastic damage evolution in the main roadway surrounding rock, the stress superposition and deflection of the peak deviatoric stress area. Through the study, the location of critical high deviatoric stress zones and support orientation of the roadway surrounding rock are obtained, and the asymmetric reinforcement support technology of the channel steel-anchor cable truss with a bi-directional reinforcement function is proposed. The experimental results show that the technique is effective.

#### **2. Project Overview**

#### *2.1. Geological Conditions*

The No.3 and No.4 coal seams are mainly exploited in the mine, and the No.4 coal seam is located 25 m above the No.3 coal seam. The relationship between the two coal seams is shown in Figure 1. Both are near horizontal coal seams adopting the comprehensive mechanized top coal caving mining method. The thickness of the No.4 coal seam is about 7.5 m, which has been recovered, and a 90 m coal pillar was left to ensure the stability of the main roadway group. The geological histogram of the long-wall working face (LWF) 8216 is shown in Figure 1. The long-wall working face (LWF) 8216 in the No.3 coal seam under the residual coal pillar is the main mining face, with an average coal seam thickness of about 5.5 m and a burial depth of about 400 m. The spacing between the two adjacent roadways is 30 m, and there is slight deformation in the surrounding rock of the roadway group during the upper coal seam recovery. Therefore, during the retreating of the lower coal seam, the mine intended to further reduce the width of the protective coal pillar for the roadway group to maximize the recovery of coal resources. The location of the project is shown in Figure 2.

**Figure 1.** Main roadway group layout of the multiple coal seam and generalized stratigraphic column in the test site.

**Figure 2.** Location of the Project.

#### *2.2. Description of the Mining Conditions*

The No.4 coal seam above the longwall working face 8216 of the No.3 coal seam in the mine has been mined out, and the width of the protective coal pillar of the main roadway group is 90 m. During the mining process of the No.3 coal seam, it usually chooses a giant protective coal pillar (90 m~140 m) to avoid the damaging effect with multi-seam mining stress disturbance to the main roadway. The main roadway is supported by a symmetrical anchor bolting and cable combination using a large coal pillar width. The anchor cables are Φ17.8 × 7000 mm with an inter-row distance of 1500 × 3000 mm and pallet size of 300 × 300 × 16 mm (length × width × height); the anchor bolts are left hand non-longitudinal rebar threaded steel anchor bolts with Φ20 × 2400 mm and inter-row distance of 900 × 1000 mm, and the pallet is a dished pallet with 150 × 150 × 10 mm (length × width × height). The roof, floor, and two ribs of the roadway use a shotcrete support. The roof and two ribs' concrete thickness is 100 mm, and the concrete grade is C20. The support parameters of the main return-air roadway and main haulage roadway are shown in Figure 3.

Although the giant protective coal pillar prevented the main roadway group from destruction and instability, it also caused the loss of coal resources. In addition, the asymmetrical damage to the main roadway can occur under symmetrical anchor rope support conditions after the suspension of mining at specific working faces of the No.3 coal seam. In the ribs of the main return-air roadway and the main haulage roadway, the support structure showed apparent damage, as shown in Figure 4. This is inconsistent with our traditional understanding of the location of the roadway damage. The roadway failure location under the influence of mining stress should first appear in the area near the working face. Therefore, in the study of multi-seam mining conditions, the direction of the evolution of the roadway damage is conducive to the better design of targeted support programs to improve support efficiency.

**Figure 3.** Support parameters of the main roadway group.

**Figure 4.** Asymmetric damage of the main roadway after stopping mining. (**A**) main return-air roadway, and (**B**) main haulage roadway.

#### **3. Theoretical Analysis of the Direction Deflection of the Stress and Partial-Stress Boosting in the Surrounding Rock**

*3.1. Solution for the Direction of Stress Deflection and Partial-Stress Boosting in the Mining Superposition Conditions*

The problems of stress distribution and plastic damage depth of a coal mass have been of great concern. Especially the superimposed loading law of the stress field with multi-seam mining, which is an essential guideline for engineering practice.

The traditional view mainly focuses on the change of lead vertical pressure (vertical direction). However, vertical stress only represents one of the nine stress components. The direction is all vertical, so the simple vertical stress cannot scientifically and comprehensively characterize the change of stress in the over-front area, caused by the mining of the working face, and it cannot effectively propose the damage characteristics of the advanced abutment pressure on the main roadway and the critical direction of control.

Deviatoric stress synthesizes three principal stresses, making the analysis more comprehensive and scientific than the single vertical stress index. According to the elastic-plasticity theory [35], the deviatoric stress plays a dominant role in the plastic deformation of the surrounding rock, so the analysis of the mining superposition main roadway surrounding stress increase law by the deviatoric stress index, can more comprehensively reflect the essence of the surrounding rock deformation damage development.

The *S*<sup>1</sup> is the principal deviatoric stress, which plays a dominant role in the stress tensor, and the commonly referred deviatoric stress refers to the maximum principle of the deviatoric stress, which is also the analytical index of the deviatoric stress in this section, and its formula is as follows.

$$S\_1 = \sigma\_1 - \frac{\sigma\_1 + \sigma\_2 + \sigma\_3}{3} \tag{1}$$

In the three-dimensional stress field, the magnitude of the principal stress and its direction can be solved by the magnitude of nine (six independent) stress components (*σxx*, *σyy*, *σzz*, *τxz*, *τxy*, *τyz*, *τzx*, *τzy*, *τyx*). Take any unitary body for an oblique section ABC, as in Figure 5, and let the outer normal of the plane ABC be *N*, whose direction cosine is

$$\cos(N, x) = l, \cos(N, y) = m, \cos(N, z) = n \tag{2}$$

**Figure 5.** Stress components of the oblique section of the element body.

The full stress *P* of the triangle ABC can be decomposed into *px*, *py*, *pz*, and the following formula is satisfied:

$$\begin{aligned} p\_x &= l\sigma\_x + m\tau\_{yx} + n\tau\_{zx} \\ p\_y &= m\sigma\_y + n\tau\_{zy} + l\tau\_{xy} \\ p\_z &= n\sigma\_z + l\tau\_{xz} + m\tau\_{yz} \end{aligned} \tag{3}$$

Let the positive stress on triangle ABC be *σN*. Then from the projection, we have

$$
\sigma\_N = lp\_x + mp\_y + np\_z \tag{4}
$$

It follows that

$$
\sigma\_N = l^2 \sigma\_x + m^2 \sigma\_y + n^2 \sigma\_z + 2mnr \tau\_z + 2nl \tau\_{zx} + 2lm \tau\_{xy} \tag{5}
$$

$$p^2 = \sigma\_N^2 + \mathbf{r}\_N^2 = p\_x^2 + p\_y^2 + p\_z^2 \tag{6}$$

Suppose there is a principal stress plane at *P*, and the projection of the full stress (principal stress *σ*) on this plane on the coordinate axis is:

$$p\_x = l\sigma, \ p\_y = m\sigma, \ p\_z = n\sigma \tag{7}$$

Satisfies Equation (9)

$$l^2 + m^2 + n^2 = 1\tag{8}$$

Equation (9) can thus be obtained

$$\begin{cases} (\sigma\_x - \sigma)l + \tau\_{yx}m + \tau\_{zx}n = 0\\ \tau\_{xy}l + (\sigma\_y - \sigma)m + \tau\_{zy}n = 0\\ \tau\_{xz}l + \tau\_{yz}m + (\sigma\_z - \sigma)n = 0 \end{cases} \tag{9}$$

Following a simplification, Equation (10) can be derived

$$\begin{cases} \pi\_{yx}\frac{m\_1}{l\_1} + \pi\_{zx}\frac{n\_1}{l\_1} + (\sigma\_x - \sigma\_1) = 0\\ \left(\sigma\_y - \sigma\_1\right)\frac{m\_1}{l\_1} + \pi\_{zy}\frac{n\_1}{l\_1} + \pi\_{xy} = 0 \end{cases} \tag{10}$$

The two equations solved for the two unknowns (*m*1/*l*1) and (*n*1/*l*1) so that *l*1, as well as *m*<sup>1</sup> and *n*1, can be obtained by bringing in the following equation:

$$l\_1 = \frac{1}{\sqrt{1 + \left(\frac{m\_1}{l\_1}\right)^2 + \left(\frac{n\_1}{l\_1}\right)^2}}\tag{11}$$

*l*1, *m*1, and *n*<sup>1</sup> are the cosine of the angle between the principal stress of the principal plane and the three axes. Therefore, based on the nine stress components at any point, the three principal stresses and the direction of the principal stresses can be solved.

#### *3.2. Solution for the Degree of the Stress Direction Deflection and the Partial-Stress Boosting and the Damage Analysis with the Disturbance of Multi-Seam Mining*

As shown in Figure 6, the stress components and the changes of the principal stress during the mining of each coal seam can be extracted through a numerical calculation. Then, the changes in the magnitude and direction of the principal stress before and after being affected by the mining superposition at any location, can be obtained. At the same time, as the magnitude of the principal stress and its direction change, the peak area of the deviatoric stress of the roadway surrounding rock is bound to change, which will inevitably lead to asymmetric damage to the surrounding rock. That is, through the size of the principal stress field of the roadway surrounding rock and its direction, changes to determine the critical damage area of the roadway, to guide the determination of the vital reinforcement direction of the roadway.

**Figure 6.** Expression of the principal stresses for the stress components. (**a**) Arbitrary stress field, (**b**) The direction deflection of the principal stress field and the partial-stress boosting.

According to the principal stress field, to solve the damage range of the roadway, using the stress distribution formula around the hole in the elastoplastic mechanics, the stress expression at any point of the surrounding rock of a circular roadway in polar coordinates [36] can be found as follows.

$$\begin{cases} \sigma\_{r} = \left[3(\lambda - 1)\cos 2\theta \frac{a^{4}}{R\_{\theta}^{4}} + (1 + \lambda)\left(1 - \frac{a^{2}}{R\_{\theta}^{2}}\right) + (\lambda - 1)\cos 2\theta - 4(\lambda - 1)\cos 2\theta \frac{a^{2}}{R\_{\theta}^{2}}\right] \frac{\gamma H}{2} \\\ \sigma\_{\theta} = \left[(1 + \lambda)\left(1 + \frac{a^{2}}{R\_{\theta}^{2}}\right) - \cos 2\theta(\lambda - 1) - 3\cos 2\theta(\lambda - 1)\frac{a^{4}}{R\_{\theta}^{4}}\right] \frac{\gamma H}{2} \\\ \tau\_{r\theta} = \left[\sin 2\theta(\lambda - 1) + 2\sin 2\theta(\lambda - 1)\frac{a^{2}}{R\_{\theta}^{2}} - 3\sin 2\theta(\lambda - 1)\frac{a^{4}}{R\_{\theta}^{4}}\right] \frac{\gamma H}{2} \end{cases} (12)$$

where: *σ<sup>r</sup>* is the radial stress at any point; *σθ* is the circumferential stress at any point; *τr<sup>θ</sup>* is the shear stress at any point; *γ* is the weight density of rock; *H* is the burial depth of the roadway; *λ* is the lateral pressure coefficient; *a* is the radius of the circular roadway; the polar coordinates of any point are set as (*θ*, *Rθ*). In the polar coordinate system, the formula for calculating the principal stress determined by the stress component [36] is:

$$\begin{cases} \sigma\_1 = \frac{\sigma\_r + \sigma\_\theta}{2} + \frac{1}{2} \sqrt{\left(\sigma\_r - \sigma\_\theta\right)^2 + 4\tau\_{r\theta}^2} \\ \sigma\_\theta = \frac{\sigma\_r + \sigma\_\theta}{2} - \frac{1}{2} \sqrt{\left(\sigma\_r - \sigma\_\theta\right)^2 + 4\tau\_{r\theta}^2} \end{cases} \tag{13}$$

The stress state of a point can be decomposed into spherical stress and deviatoric stress, where the deviatoric stress is the dominant factor leading to the plastic deformation of the surrounding rock, so use the deviatoric stress for the theoretical analysis.

When conducting the study of the stress in the roadway enclosure, the problem can be simplified to a plane problem, which results in the principal deviatoric stress calculation formula:

$$S\_1 = \frac{\gamma H}{6}(1+\lambda) + \frac{\gamma H \overline{a}^2 (1-\lambda) \cos 2\theta}{3R\_\theta^2}$$

$$= 1 + \frac{\gamma H}{2} \sqrt{\left[ (1+\lambda)(-\frac{\overline{a}^2}{R\_\theta^2}) + (\lambda - 1)(1 + 3\frac{\overline{a}^4}{R\_\theta^4} - 2\frac{\overline{a}^2}{R\_\theta^2}) \cos 2\theta \right]^2 + \left[ (\lambda - 1) \left( 1 - 3\frac{\overline{a}^4}{R\_\theta^4} + 2\frac{\overline{a}^2}{R\_\theta^2} \right) \sin 2\theta \right]^2} \tag{14}$$

The surrounding rock plasticity criterion satisfies the following relation [37,38].

$$\frac{\sigma\_{\mathfrak{c}}^{2}(\sigma\_{\mathfrak{t}}+\sigma\_{\mathfrak{t}})}{2\left(\sigma\_{\mathfrak{c}}+2\sigma\_{\mathfrak{t}}\right)} + \frac{\sigma\_{\mathfrak{t}}\sigma\_{\mathfrak{c}}^{2}}{\sigma\_{\mathfrak{c}}+2\sigma\_{\mathfrak{t}}} - \frac{\sigma\_{\mathfrak{c}}^{4}}{8\left(\sigma\_{\mathfrak{c}}+2\sigma\_{\mathfrak{t}}\right)^{2}} - \frac{\left(\sigma\_{\mathfrak{t}}-\sigma\_{\mathfrak{t}}\right)^{2}}{2} = 0\tag{15}$$

By combining the deviatoric stress formula and its plasticity criterion, we can analyze the location of the high deviatoric stress zone in the roadway surrounding rock. We can also make it possible to quantify the deflection of the high deviatoric stress zone when the principal stress direction changes. Combined with the numerical simulation tests, the evolution of the peak stress zone of the surrounding rock in the roadway, from no stress deflection and partial-stress boosting to stress deflection and partial-stress boosting conditions, is shown in Figure 7.

**Figure 7.** Evolutionary pattern of surrounding rock damage under different stress direction deflection and partial-stress boosting conditions.

With different degrees and directions of stress deflection and partial-stress boosting, there exist different deviatoric stress critical zones in the roadway surrounding rock, and this zone highlights the roadway surrounding rock's key damage development area and direction: When the degree of stress deflection and partial-stress boosting is small, the roadway surrounding the rock's deviatoric stress critical zone mainly develops symmetrically in two directions of the roadway; when the degree continues to increase, the deviatoric stress critical zone is symmetrically distributed on both sides of the roadway surrounding rock, and the damaged area tends to develop from four directions to the deep part of the surrounding rock. Therefore, under stress deflection and partial-stress boosting conditions, the surrounding rock is not gradually and uniformly damaged to the deep part of the surrounding rock, but it is a non-uniform process. The different stress direction deflection and partial-stress boosting conditions lead to two directions of damage development or four directions of the damage development pattern of the surrounding rock, which has an essential guidance significance for the design of the surrounding rock control scheme.

#### **4. Numerical Simulation**

The influence with multi-seam mining will cause a stress increase effect on the main roadway. Under different principal stress field conditions, the roadway's stress peak area and damage development area will be deflected, forming the phenomenon of stress deflection and partial-stress boosting. There is an essential difference between the development of the surrounding rock deformation and the critical area of its control. The numerical simulation can more intuitively study the law of the deflected loading of the main roadway surrounding rock, caused by mining multiple coal seams. This is essential for analyzing the plasticization law of the main roadway surrounding rock and determining the critical reinforcement location of the main roadway surrounding rock.

The model selects the LWF 8216 as the test working face. The *X*-axis length is 400 m in the length direction, the *Y*-axis length is 100 m and parallel to the main roadway group, and the *Z*-axis length is 90 m in the vertical direction, as shown in Figure 8. The boundary displacement of the model is constrained in the horizontal and bottom. The upper of the model is subjected to a stress of 7.75 MPa, equivalent to the weight of the overburdened rock, and the lateral stress coefficient is 1.2. The model is calculated using the Mohr– Coulomb model. In contrast, the double-yield model is used for the job. The mechanical parameters of the coal and rock mass were very important to the accuracy of the numerical simulation results. Cai et al. [39] suggested that the tensile strength, cohesion and elastic modulus of the coal and rock mass could be estimated in the range of 0.10–0.25 times of the experimental results in the laboratory. Mohammad et al. [40] suggested that the average uniaxial compressive strength of the numerical model should be 0.284 of the laboratory strength and the average stiffness should be 0.469 of the laboratory stiffness. Combined with the generalized Hoke–Brown failure criterion [41–43], the simplified physical parameters of the coal seam are shown in Table 1.

**Figure 8.** Numerical model overview.

**Table 1.** Properties of the rock mass in LWF 8216.


The numerical model is used to study two aspects: (1) The change of the abutment pressure in front of the working face when the upper and lower coal seams are in different positions (outer fault, overlap, internal fault). (2) To investigate the stress distribution and the mechanism of the surrounding rock damage in the surrounding rock of the roadway when the width of the coal pillar is different in the working face of the lower coal seam, when the protective pillar of the upper coal seam is 90 m.

#### *4.1. Stress Field Advance of the Working Face, Caused by Multi-Seam Mining*

#### 4.1.1. Deflection of the Principal Stress

Multi seam mining generally adopts a downward mining practice, so the relative position relationship between the underlying coal seam and the overlying coal seam mainly includes three types. Namely: external staggering (the end-mining line of the lower coal seam is under the gob of the upper coal seam), overlapping (the end-mining line of the upper and lower coal seams overlap), and internal staggering (the end-mining line of the lower coal seam is under the entity coal of the upper coal seam). Suppose the width of the final end-mining coal pillar of the lower coal seam is smaller than that of the upper coal seam. In that case, the working face of the lower coal seam will experience the whole process of the above three types of positional relationships.

The stress component of each point before each working face, is extracted through a numerical calculation. The maximum principal stress direction of the leading working face under three positional relations, is calculated using the stress component. It is specified that the angle of the maximum principal stress line rotating counterclockwise to the vertical direction line is "+" and vice versa is "−". The absolute value of the principal stress direction angle is less than 90◦. When the upper and lower coal seams are at different off-set distances, the deflection characteristics of the principal stress of the element in the center of the lower coal seam with variable distances from the working face are shown in Figure 9.

**Figure 9.** Deflection law of the principal stress in front of the different layout modes of the upper and lower coal seams. (**A**) External-offset of 30 m. (**B**) Overlapping arrangement. (**C**) Internal-offset of 30 m.

The primary conclusions are as follows.


ping type → internal staggered type. It is reasonable to adopt symmetrical support measures for the roadway at the turning point of the principal stress.

(3) Under the condition of multi-seam mining, the maximum principal stress direction within 60 m in front of the lower coal seam is approximately vertical and deflected to the mining void. If the roadway is arranged here, the roadway surrounding rock will be damaged first along the direction of the minimum principal stress, namely the entity coal top angle of the roadway.

#### 4.1.2. Direction Deflection of the Principal Stress

Under the influence of the repeated mining of multiple coal seams, the bearing stress in front of the workings will be subjected to the cyclic loading-unloading effect of mining stress. The geological rock body shows a typical non-homogeneity in the vertical stress direction, and the nature of the different rock seams varies. Therefore, the different relative positions of the upper and lower coal seams will significantly affect the peak and distribution characteristics of the stress field in front of the workings. The vertical stress of the model is selected as the index, and the distribution of the bearing stress field of the working face of the lower coal seam, when the upper and lower coal seams are at different staggered distances, is shown in Figure 10.

**Figure 10.** Abutment stress field distribution in the lower coal seam with different offset distances. (**A**) External-offset of 30 m. (**B**) Overlapping arrangement. (**C**) Internal-offset of 30 m.

As shown in Figure 10:


#### *4.2. Roadway Group Stress Field Characteristics under Various Protective Coal Pillars*

The relative positional relationship of the upper and lower coal seams working face significantly affects the abutment stress distribution condition in front of the working face. Due to the superimposed influence of the mining stress in multiple coal seams, the stress loading law of the main roadway group in the process of the upper and lower coal seams working face going from an outer staggered arrangement → overlapping arrangement → inner staggered arrangement, is significant for the selection of the width of the protection coal pillar. This section explores the stress evolution process of this seam's main roadway group affected by mining when the end-mining coal pillar of the upper coal seam is 90 m, with the gradual mining of the lower coal seam, as shown in Figure 11.

As shown in Figure 11:


surrounding rock is asymmetric, so it is necessary to use targeted support solutions in some places where the surrounding rock is weak.

**Figure 11.** Stress evolution process of the main roadway group during the lower coal seam recovery. (**A**–**E**) with the protective coal pillar widths of ∞ m, 120 m, 90 m, 60 m, and 30 m, respectively, and (**F**) distribution of peak stress value in both ribs.

#### *4.3. Law of Stress Direction Deflection and Partial-Stress Boosting of the Deviatoric Stress Field of the Surrounding Rock with Different End-Mining Coal Pillar Widths*

*Deviatoric stress* is the stress that deviates from hydrostatic stress and causes deformation, reflecting the essence of plastic deformation at a point in the surrounding rock. The high deviatoric stress area in the surrounding rock indicates that the area is in a state of damage development. Mining activities, especially the disturbance of multi-seam mining, will promote the direction deflection of stress and partial-stress boosting effect of the main roadway surrounding rock, further leading to the formation and transfer of the high deviatoric stress zone in the surrounding rock. It is of great significance to investigate the development and evolution of the high deviatoric stress zone in the surrounding rock

to develop the damage mechanism and key reinforcement technology of the roadway surrounding rock. This section explores the evolution of the peak deviatoric stress zone and plastic zone of this seam's main roadway group affected by mining when the end-mining coal pillar width of the upper coal seam is 90 m, with the gradual mining of the lower coal seam, as shown in Figure 12.

**Figure 12.** Development characteristics of the deviatoric stress peak zone and plastic zone. (**A**–**E**) with the protective coal pillar widths of ∞ m, 120 m, 90 m, 60 m, and 30 m, respectively.

As shown in Figures 12 and 13:


**Figure 13.** Stress distribution of the main roadway group in the lower coal seam.

#### **5. Surrounding Rock Control Techniques**

#### *5.1. Principle of the Asymmetric Reinforcement Support Cooperation*

In the process of the lower coal seam mining, with the reduction of the width of the protective coal pillar of the main roadway group, the principal stress magnitude and direction will change in the main roadway surrounding rock. In the process of the mining influence, there is critical area of the direction deflection of stress and partial-stress boosting in the main roadway surrounding rock, and the area of the direction deflection of the stress and partial-stress boosting is the peak area of the deviatoric stress in the main roadway surrounding rock, which is also the precursor and leading area for the development of damage or continuous damage in the surrounding rock. Thus, the principle of the asymmetric directional cooperative anchorage in the critical area of the stress direction deflection and the partial-stress boosting, is proposed in the roadway surrounding rock. Based on the evolution process of the stress direction deflection and the partial-stress boosting area of the roadway surrounding rock, a directional high prestressing anchorage reinforcement technology is adopted to realize the synergistic anchorage with the original support scheme and guarantee the stability of the main roadway group.

Based on the study mentioned above, the multi-coal seam mining disturbance leads to a significant change in the magnitude and directional deflection of the principal stress. In a comprehensive view, there are three main categories of the deflected load critical areas in the upper and lower coal seam main roadway: two ribs area, roof area, floor area, and rib angle area (two types). The critical area of the directional cooperative anchorage is the peak area of the deviatoric stress, the key area of the stress direction deflection and the partial-stress boosting of the main roadway surrounding rock. For this reason, four types of essential roadway surrounding rock reinforcement anchorage schemes are proposed for the different degrees of the direction deflection of the stress and partial-stress boosting, as shown in Figure 14.

**Figure 14.** Principle of the directional cooperative anchoring in the roadway.

In the principle of the high pre-stressed anchor cable reinforcement anchoring effect, the average equivalent mechanical parameters of the surrounding rock in the critical area of the stress direction deflection and the partial-stress boosting can be increased from *C*s0 and *ϕ*s0 to the equivalent *C*s1 and *ϕ*s1, *C*s2 and *ϕ*s2, or *C*s3 and *ϕ*s3, by the high prestress and high shear failure resistance provided by the anchor cable, where *C*s3 > *C*s2 > *C*s01 > *C*s0 and *ϕ*s3 > *ϕ*s02 > *ϕ*s01 > *ϕ*s0. This means that the mechanical properties of the surrounding rocks in the key zone of the direction deflection of the stress and partial-stress boosting are strengthened to improve the overall bearing capacity of the surrounding rocks of the main roadway.

#### *5.2. Support Parameters*

Based on the principle mentioned above of the directional cooperative anchoring of the roadway, a targeted support scheme is put forward for the main return-air roadway of the No.3 coal seam. Based on the original support scheme, two anchor cables are additionally supported on the mining side of the roadway, with an angle of 15◦ with the horizontal, the anchor cable specification of Φ17.8 × 7000 mm, and the row spacing of 1500 mm × 2000 mm. The same row of anchor cables and each row of anchor cables are connected with the W steel strip, respectively, forming a crisscross net-like combined support structure, increasing the surface stress range of the roadway surrounding rock. Add one anchor cable of Φ17.8 × 7000 mm with 2000 mm row spacing on the entity coal side and the W steel strip connects each row of the anchor cables. Detailed support parameters are shown in Figure 15.

**Figure 15.** No. 3 coal seam main return-air roadway reinforcement support program. (**A**) Entry section. (**B**) Support pattern in roof. (**C**) Support patter in coal pillar rib. (**D**) Support pattern in virgin coal rib.

In order to further analyze the rationality of the directional cooperative anchoring support scheme, the support stress field of the No. 3 coal seam main return-air roadway is shown in Figure 16. A deep compressive stress zone and a shallow core compressive stress zone are formed in the surrounding rock after adopting the bolt-cable combined support. The stress field of the anchor rods of the two ribs avoided the fall of the broken surrounding rock of the roadway surface but did not form a practical strengthening effect on the peak area of the deviatoric stress. The principle of support reinforcement is shown in Figure 17. Following the use of the asymmetric reinforcement support with anchor cables, the channel steel anchor cable reinforced mining rib formed two functional compressive stress areas in the center of the deviated stress core area. This can effectively deliver the external mining stress of the surrounding rock to the deep elastic rock, thus inhibiting the damage to the roadway surrounding rock in this area.

**Figure 16.** No.3 coal seam main return-air roadway support stress field.

**Figure 17.** The support reinforcement principle of the No.3 coal seam main return-air roadway. (**A**) No support. (**B**) Original support. (**C**) Strengthen support.

#### *5.3. Ground Pressure Monitoring*

To analyze the mining pressure response of the working face and the return-air roadway of this coal seam during the pushing process of the C3# coal seam, we arranged three measurement stations (each with an interval of 15 m) along the axis of the working face return-air road and the C3# coal seam return-air roadway, respectively. The observation indexes of the stations are mainly: (1) the abutment stress of the mining roadway; (2) the stress of the bolts in the mining roadway and main roadway; (3) the amount of surrounding rock displacement. The changing trend of these indicators is analyzed to reveal further the law of change of the stress field in the surrounding rock of the main roadway when the width of the protective coal pillar is different, to determine the width of the protective coal pillar in the main roadway of LWF 8216. The monitoring principle is as Figure 18.

**Figure 18.** Ground pressure monitoring arrangement.

The location of the drilling stress monitoring is the C3# coal seam working face mining roadway. Each measurement point has two stress gauges, A and B, where A is 14 m deep into the coal body, and B drill holes are 8 m deep into the coal body. The results of the LWF 8216 track roadway's abutment stress monitoring are shown in Figure 19.

**Figure 19.** Monitoring results of the deep hole stress meters.

When the distance between the working face and the stress meter is more than 90 m, the reading of the stress meter is stable, and there is no apparent stress-increasing trend. However, the reading shows a slow and slightly increasing trend when the working face is 90~50 m from the stress meter. When the distance between the working face and the stress meter is less than 50 m, and the face is closer to the monitoring section of the high-sensitivity borehole stress station, the abutment stress value begins to show an apparent acceleration trend in front of the working face.

The anchor cable support resistance can reflect the deformation state of the roadway surrounding rock and the influence of the working face in real time. It can effectively determine the reasonable width of the end-mining coal pillar of the main roadway group. Monitoring stations were set up in the LWF 8216 track roadway, and the C3# coal seam main return-air roadway, respectively, and the stress meters were arranged, as in Figure 18; the monitoring results are shown in Figure 20.

**Figure 20.** Anchor bolts stress meter monitoring results. (**A**) The main haulage roadway. (**B**) The main return-air roadway.

The pre-stress range of the anchor bolts in the LWF 8216 track roadway is 50~53 KN. The values of the stress meters are 163~170 KN after the working face is stopped mining and is stabilized, which is about 90.6~94.4% of its breaking load, and still in a safe working condition; with the decrease of the distance between the working face and the anchor meter, the values of stress meters goes through the process of basic stability (more than 100 m) → slow increase (100~40 m) → accelerated increase (less than 40 m). This indicates that the stress of the advanced mining will be significantly enhanced within 40 m in front of the working face.

In the C3# coal seam main return-air roadway, the anchor pre-stress values range from 49 to 53 KN. Once the working face has stopped mining 56 m from the main roadway, the anchor stress meter of the main roadway showed an apparent asymmetry, and the values of the two ribs were significantly larger than the roof. With the decrease of the distance between the working face and the anchor bolts stress meter, the values of stress meters went through the process of a primary stability (greater than 100 m) → slow rise (100~56 m), and the final value of the stress meters was about 34.4.0~48.3% of its breaking load. The comparison and analysis of the values of the stress meters (G1~G6), G1 > G2 > G5, G6 > G3 > G4, show that the surrounding rock of the two sides of the main roadway, especially the mining side, is more sensitive to the mining effect of the working face. The field monitoring results are consistent with the numerical simulation results of the critical control area of the surrounding rock, indicating that the theory of the stress direction deflection and the partialstress boosting of the surrounding rock is suitable for early warning and strengthening the support of the surrounding rock deformation.

#### **6. Conclusions**

In this paper, the theoretical analysis, numerical simulation, and on-site ground pressure monitoring are used to study the effects of the partial-stress boosting and the direction deflection of the stress field on the failure evolution of the surrounding rock. The conclusions are as follows:


**Author Contributions:** Conceptualization, D.C.; Data curation, Z.L.; Formal analysis, Z.L. and X.M.; Funding acquisition, D.C.; Investigation, X.M., Y.W. and Z.W.; Methodology, D.C., F.G. and S.X.; Project administration, D.C. and S.X.; Software, D.C., F.G., Z.L., X.M. and Y.W.; Supervision, S.X.; Validation, Y.W. and Z.W.; Writing–original draft, D.C., F.G. and Z.L. Writing–Review and Editing, D.C., F.G. and Z.L. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was financially supported by the National Natural Science Foundation of China (Grant No. 52004286), the Fundamental Research Funds for the Central Universities (Grant No. 2022XJNY02), the China Postdoctoral Science Foundation (Grant No. 2020T130701, 2019M650895), all of which were gratefully acknowledged.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**

