Study on the Instability Characteristics and Bolt Support in Deep Mining Roadways Based on the Surrounding Rock Stability Index: Example of Pansan Coal Mine

In view of the influence of mining stress on the stability of the surrounding rock of inclined roof mining roadways in deep mines, the surrounding rock stability index is defined and solved based on the rock strength criterion and the stress distribution. /e mining roadway of the 17102(3) working face of the Pansan Coal Mine is used as the engineering background and example. /e surrounding rock’ stabilities under the conditions of no support and bolt support are analyzed according to the surrounding rock’s stability index and the deformation data. /e results show that the areas of low wall and high wall instability are 1.68m and 2.12m, respectively, and the low wall is more stable than the high wall; the areas of the roof and floor instability are 0.33m and 0.35m, respectively, and the roof and floor are more stable than the two sides. During mining, the area of instability greatly increases at first, then decreases to 0, and reaches a maximum value at the peak of the abutment pressure. /e stability of the surrounding rock decreases first and then increases. Compared with the end anchoring bolt support, the full-length anchoring bolt support reduces the area of instability to a greater extent, and the full-length anchoring bolt support effect is better. /e surrounding rock in the end anchoring zone and the full-length anchoring zone began to deform significantly at 200m and 150m from the working face, respectively. /is indicates that the control effect of the full-length anchoring bolt support is better and verifies the rationality of the surrounding rock stability index to describe the instability characteristics. /is research method can provide a theoretical reference for analysis of the stability characteristics and support design of different cross-section roadways.


Introduction
With the increasing depletion of shallow coal resources, deep mining in various mining areas has become the norm in China. In high ground stress environments, a deep mining roadway often shows the characteristics of soft rock. It is difficult to control the deformation of the mining roadway when it is affected by mining stress. Effective control of the surrounding rock deformation of the mining roadway is the premise to ensure smooth mining of the working face. In deep coal resources, the inclined coal seam is widely distributed and occupies a large proportion of the area. When mining an inclined coal seam, to facilitate the management of the roof of the mining roadway, the mining roadway is often driven along the roof of the coal seam, forming an inclined roof mining roadway (IRMR). With the characteristics of supporting fast driving speeds and simple construction technology, IRMR is widely used in inclined coal seams. e stability characteristics of the surrounding rock of IRMR have their own features, which has attracted the attention of the field engineers and rock mechanics workers.
As the core content of roadway support theory, the stability of a roadway's surrounding rock is always one of the key research directions in mining engineering. Many scholars at home and abroad, using field engineering practices, experimental results, and theoretical reasoning analyses, have carried out in-depth research on the stability of the surrounding rock and bolt support mechanism from different angles and have established many classic theories, such as the combined arch theory [1], the plasticity zone theory [2,3], and the broken rock zone theory [4]. In IRMR, because of the existence of corners, when affected by mining stress, there is a large concentration of stress in these corners, resulting in a large difference between the surrounding rock stability characteristics of IRMR and that of regular crosssection roadways such as circular, elliptical, and straight wall arches.
is study on the stability characteristics of surrounding rock and the bolt supporting mechanism of IRMR affected by mining stress is of great significance to the surrounding rock control of IRMR in deep mining. Based on this, this paper proposes a surrounding rock stability index (SRSI), which can quantitatively describe the stability of the surrounding rock based on the stress distribution of IRMR and the rock strength criterion. SRSI provides the theoretical basis for the study of the evolution law of the surrounding rock stability characteristics and the bolt support mechanism of different anchorage forms of IRMR under the influence of mining stress.
In terms of research on the surrounding rock stability of roadways with irregular cross sections, such as inclined roof cross sections, the numerical simulation method is mainly used to carry out relevant research. Some of the research results are as follows: G. Li and others used FLAC3D numerical simulation software to study the distribution of the surrounding rock plastic zone, the principal stress difference, and the deformation characteristics of six tunnel section shapes, namely, rectangular, straight walls with arched tops, horseshoe shaped, three centered arches, circles, and ellipse, and proposed the concept of equivalent excavation and void reinforced areas [5]. Q. Meng and others used FLAC3D numerical simulation software to study the surrounding rock deformation characteristics and the surrounding rock distribution of the plastic zone of six kinds of cross-section shape roadways, namely, rectangular, trapezoid, straight wall arches, horseshoes, ovals, and round shapes. e results show that the round and oval shapes are the optimal cross-section shapes for deep roadway under high ground stress [6]. L. Hu and others used ANSYS software to simulate and analyze the influence of roadway section structure on the surrounding rock deformation. e results show that the section form has an appreciable impact on both the horizontal displacement and the vertical displacement of the surrounding rock under the same conditions, and the rectangular section has the greatest impact [7]. Although numerical simulation can be used to study and solve a wide range of engineering technical problems, numerical simulation can only obtain quantitative conclusions, and some problems cannot be judged qualitatively. At the same time, the research method of numerical simulation is still based on a theoretical basis, and the imperfection of theoretical research limits the application and development of numerical simulations.
In view of the research on the bolt support mechanism with different anchorage forms, scholars at home and abroad have conducted thorough studies and achieved many results. For the end anchored bolt support, there are mainly classical theories such as the combined arch theory [1], the broken rock zone theory [3], the suspension theory [8,9], and the dynamic response theory [10], which play an irreplaceable role in the design of bolt support. For the full-length anchoring bolt support, there have been numerous studies, and some important achievements are listed as follows: J. Chen and others used an analytical model to study the load transfer behavior of fully grouted rock bolts. e results show that the axial load in the rock bolt decayed from the loaded end to the free end independent of the pull-out load. However, the trend of the load distribution curve was influenced by the pull-out load [11]. M. H. Aghchai and others presented an analytical method to use this curve for determining the bond (bolt-grout and grout-rock interface) shear strength parameters [12]. Cui and others aimed to propose new numerical procedures for elastoplastic analysis of the interactions between fully grouted bolt and rock with slip and nonslip cases at a potential-failure interface [13]. Zou and Zhang investigated the dynamic evolution characteristic of the bond strength at the interface of a bolt and a rock mass under an axial tensile load and the mechanical behavior of fully grouted bolts in situ considering the nonuniform stress of the surrounding rock along the bolts [14]. C. Wu and others studied the tensile behavior of unreinforced and reinforced rock with fully grouted bolts and an empirical approach was proposed to predict the strength of rock reinforced with fully grouted bolts [15]. Liu and others established an analytical model for the interaction between the bolt and surrounding rock based on the bearing mechanism of fully grouted rock bolts, deduced the corresponding controlled differential equation for load transfer, and obtained the stress distributions of the anchorage body [16]. Chen and others employed the elastic theory of a semi-infinite length beam and approximate differential equation of the deflection curve (elastic and elastoplastic) to analyze the lateral behavior of a bolt under surrounding rock pressures [17]. Although there is much research on the mechanisms of bolt support with different anchorage forms in the existing literature, most of them are based on the condition that the mechanical state of the surrounding rock remains unchanged. However, owing to the influence of mining stress, the mechanical state of the surrounding rock is obviously changed, so it is necessary to consider the influence of the change of the surrounding rock's mechanical state when studying bolting mechanisms.
To summarize, this paper proposes a method for the surrounding rock stability analysis with reference to previous research results. Based on the strength criterion and the stress distribution of the surrounding rock, SRSI is defined to measure the surrounding rock stability characteristics. Using the Mohr-Coulomb strength criterion and the complex function method proposed by Muskhelishvili [18][19][20][21][22][23], the analytical solution of SRSI is obtained. Based on the study of the 17102(3) working face of the Pansan Coal Mine (PCM) in the Huainan, China, the area of the unstable area (AUA) of the surrounding rock is calculated according to SRSI, and then the evolution law of surrounding rock stability characteristics is analyzed when IRMR is affected by mining stress in deep mining. By changing the boundary conditions, the distribution of the SRSI with supporting conditions and the change law of AUA with the mining of the working face are obtained. Combined with the observation data of the IRMR deformation, the supporting effect and mechanism of the end anchoring bolt support and the full-length anchoring bolt support are analyzed. Finally, this paper summarizes the research results and gives relevant conclusions.

Definition and Solution of SRSI.
e strength criterion of rock can be expressed by an equation describing the function relationship between the principal stresses: (1) In equation (1), σ 1 represents the maximum principal stress of the critical failure state; σ 2 represents the intermediate principal stress; σ 3 represents the minimum principal stress; function f represents the strength criterion. When the real maximum principal stress σ 1 is less than the critical failure maximum principal stress σ 1 , the surrounding rock is in a stable state. e smaller σ 1 is, the more stable the surrounding rock is. On the contrary, when the surrounding rock is in a state of instability, the larger σ 1 is, the worse the stability of the surrounding rock is. According to this, SRSI can be defined as In equation (2), K represents SRSI. SRSI reflects the relationship between the actual maximum principal stress σ 1 and the critical failure maximum principal stress σ 1 . When K is equal to 1, the surrounding rock is in a critical state of failure; when K is greater than 1, the surrounding rock is in a stable state; when K is less than 1, the surrounding rock is in a state of failure, that is, an unstable state, so this area is called an unstable area (UA). According to this, AUA can be defined to evaluate the stability of the surrounding rock and the bolt support effect. e smaller AUA is, the more stable the surrounding rock is. e definition of AUA is as follows: In equation (3), Ω represents the specific area in the surrounding rock, which can be the entire surrounding rock, the roof and floor, the two sides, or other designated areas; H represents the Heaviside step function, and the Heaviside function H is as follows: In this paper, the basic form of the Mohr-Coulomb rock strength criterion is used to calculate K. When the Mohr-Coulomb strength criterion is expressed by principal stress, it is as follows: In equation (5), ϕ represents the internal friction angle and c represents cohesion. By substituting equation (5) into equation (2), the expression of K under the Mohr-Coulomb criterion can be obtained as follows: e complex function method proposed by Nusyfmjzcjmj is used to solve the principal stress distribution of the surrounding rock of IRMR [18][19][20][21][22][23]. e specific ideas are as follows: first, the outer area of IRMR is mapped to the unit circle on the complex plane with a conformal mapping function; second, the complex function of the representative stress distribution is solved by using relevant knowledge of the complex analysis according to the boundary conditions; finally, the principal stress distribution is calculated with the complex function.
e deep buried IRMR can be regarded as an infinite plane orifice. e vertical pressure σ V , horizontal pressure σ H , and shear stress τ are used to represent the total pressure on the surrounding rock. e stress diagram of IRMR is shown in Figure 1.
According to the Riemann mapping theorem [18,24], the surrounding rock area of IRMR can be mapped into a unit circle via the conformal mapping function ω(ξ). In general, it is very difficult to find the exact conformal mapping function to use. An effective method is to expand it into a Laurent series and take the finite term as the approximation of the conformal mapping function ω(ξ). In this case, the conformal function ω(ξ) can be expressed as In equation (7), z � x + iy represents the coordinate of the point on the z plane where the roadway section is located; ξ � ρe iθ represents the coordinate of the point on the ξ plane where the unit circle is located; C n � A n + iB n represents the complex coefficient, whose value is related to the shape of the roadway section; N represents the order of the Laurent series. en the search for an exact conformal mapping function is transformed into a search for a group of the complex coefficients C n for approximation. According to the algorithm proposed in [24][25][26][27][28][29], this paper uses scientific calculation software to program and solve the complex coefficient C n . e corresponding relationship between the coordinates of the points on the z plane and the ξ plane is as follows: x � Re ω ρe iθ , y � Im ω ρe iθ .
In equation (8), the function Re represents the real part of the complex number; the function Im represents the imaginary part of the complex number.
Advances in Civil Engineering e stress of the surrounding rock can be determined with two complex functions, φ(z) and ψ(z), in the complex plane. To solve the distribution of the principal stress in the surrounding rock of IRMR, the form of the complex function is as follows: In equation (9), the complex functions φ 0 (ξ) and ψ 0 (ξ) represent the analytic function satisfying the Cauchy-Riemann condition in the definition domain, which can be expanded to Taylor series, as shown in equation (10). And the α value is as shown in equation (11). One has φ 0 (ξ) � n�+∞ n�0 a n ξ n , In equation (10), a n and b n represent the complex coefficients. e complex functions φ 0 (ξ) and ψ 0 (ξ) can be obtained using the boundary conditions of the roadway surface. According to Harnack's theorem [18], the boundary conditions of the roadway surface can be equivalent to two functional equations as In equation (12), Γ represents the boundary of the roadway section; ζ represents the point on the boundary Γ; ; f 0 and f 0 represent the reduced boundary condition as follows: By substituting equation (10) into equation (12), the integral is solved by using the Cauchy integral formula. en, the complex coefficients a n and b n can be obtained by comparing the coefficients on both sides of equation (12), and the complex functions φ 0 (ξ) and ψ 0 (ξ) can be obtained. en, the stress distribution of surrounding rock can be obtained with the following equation: In equation (14), σ ρ represents the normal stress in the direction of the coordinate line ρ; σ θ represents the normal stress in the direction of the coordinate line θ; τ ρθ represents the shear stress in the direction of the coordinate lines ρ and θ.
Because the stress calculated by equation (14) is in an orthogonal curvilinear coordinate system determined by the By substituting the main principle into equation (6), the analytical solution of SRSI can be obtained.

Effect of Mining Stress.
Mining stress mainly affects the surrounding rock stress conditions and the strength parameters of the surrounding rock. According to the research results on the abutment pressure in [30][31][32][33][34], the exponential model is used to describe the abutment pressure, which is used to estimate the vertical pressure σ V (l) of the surrounding rock: In equation (16), l represents the distance from the working face; σ 0 V represents the vertical pressure of the original rock; k s represents the vertical stress concentration coefficient at the peak of the abutment pressure; l s represents the peak position of the abutment pressure; coefficients a and b are used to describe the change severity of the abutment pressure curve. e lateral pressure coefficient λ and shear stress coefficient η are defined as Assuming that the lateral pressure coefficient λ and shear stress coefficient η are constant in the mining process, the horizontal pressure σ H (l) and shear stress τ(l) of the surrounding rock of the roadway are It is assumed that the cohesion of the surrounding rock is constant before the peak of the abutment pressure, decreases linearly after the peak of the abutment pressure, and is 0 at the working face. e internal friction angle of the surrounding rock changes little in the mining process, which is assumed to be constant. e cohesion c(l) and the internal friction angle ϕ(l) are as follows: In equation (19), c 0 represents the cohesion of the protolith, and ϕ 0 represents the internal friction angle of the protolith.
By substituting equations (16) and (18) into equation (9), the distribution of the principal stress of the surrounding rock of IRMR under the influence of mining can be obtained. By substituting the principal stress distribution and equation (19) into equation (6), the analytical solution of K under the influence of mining can be obtained. To facilitate the analysis of the stability of the surrounding rock of the roadway, the areas of the roof and floor and the two sides are divided. In the middle of the roof, the floor, and the two sides, a radial survey line is arranged in the vertical rock wall. A radial survey line is arranged in the diagonal direction at the corners of the roadway. A circular survey line is arranged on the roadway surface and 0.3 m into the surrounding rock. And the corners of roadway are numbered. e roadway section shape, the survey line position, the number of corners, the two-side area, and the roof and floor area are shown in Figure 2. In Figure 2, the orthogonal curvilinear coordinate system determined by the conformal mapping function ω(ξ) is also drawn. It is illustrated as a spider network coordinate system composed of circular coordinate line ρ and radial coordinate line θ in the figure.

Distribution Law of Mining Stress in Roadway.
e roof and floor of the 17102(3) working face are mudstone and sandy mudstone, which are similar to the mechanical properties of the coal seam and thus are combined. According to the in situ stress test results, σ 0 V is 16.8 MPa; σ 0 H is 13.3 MPa; τ 0 is 0.5 MPa; λ is 0.79; η is 0.03. According to [30][31][32][33][34], as well as field mine pressure observation results, mining height, mining speed, hydraulic support parameters, and other relevant data, let k s � 1.9, a � 0.1, and b � 0.03. When the distance between the roadway and the working face is more than 250 m, the roadway is not affected by mining stress. Substitute the above parameters into equation Advances in Civil Engineering (16) to obtain the abutment pressure curve, as shown in Figure 3. It can be seen from Figure 3 that σ V increases first and then decreases with mining, reaching a maximum value k s σ 0 V at l � l s . After the laboratory test and using different weights to determine the reduction coefficient, we obtain the rock mechanics parameters, which are more consistent with the actual site. When not affected by mining stress, let the cohesion c 0 be 6.2 MPa and let the internal friction angle ϕ 0 be 28°. According to equation (19), the change in the surrounding rock cohesion of the roadway with the mining can be obtained as shown in Figure 4. It can be seen from Figure 4 that c is a constant before the peak of the abutment pressure, and after the peak, c decreases linearly, reaching a minimum value of c � 0 at the working face.

Distribution Law of SRSI in IRMR without Support.
Without considering roadway support, SRSI in the vertical section of the 17102(3) working face is obtained at different positions (l) in front of the advancing direction, and SRSIs of radial survey line and the circumferential line are obtained.
ey are as shown in Figures 5-7, respectively. It is generally believed that when l > 250 m, it is not affected by mining. It can be seen from Figure 5 that the distribution law of SRSI is similar when the roadway is affected by mining stress and when it is not affected by mining stress. In the roof and floor and the two sides of the roadway, there are unstable areas of the surrounding rock. It can be seen from Figure 6(a) that the minimum value of SRSI at the roof and floor radial survey lines is greater than 0, as shown by arrows b and d in Figure 6(a), indicating that the UAs of the roof and floor are not connected, and the minimum value of SRSI at the two-side survey lines is less than 0, as shown by arrows a and c in Figure 6(a), indicating that the UAs of the two sides are connected. From Figure 5, it can be clearly observed that there are two blue areas in the roof and floor, which are located at the two corners, respectively. With the mining of the working face, there is a trend of connectivity, indicating that the UAs of the roof and floor are located at the two corners, respectively, and with the mining of the working face, there is also a trend of connectivity. ere is a blue area in each of the two sides, which is the arched distribution, and the arch corners are located at the roadway corners, indicating that the UAs of the two sides are an arched distribution and the arch corners are located at the roadway corners. e minimum value of SRSI is 0, and the minimum value is located at the roadway corners, indicating that the stability of the roadway corners is worst in the surrounding rock of the roadway.
It can be seen from Figure 6(b) that, on the diagonals of the roadway, SRSI tends to be constant when it approaches the corners from the depth of the surrounding rock, as shown by arrows e, f, g, and h in Figure 6(b). It can be seen from Figure 7(a) that, on the roadway surface, SRSI tends to be 0 when approaching from the middle of the rock wall to the corners, as shown by arrows i, j, k, and l in Figure 7(a). It can be seen from Figure 7(b) that SRSI near the corners increases suddenly, as shown by arrows m, n, o, and p in Figure 7(b). Based on this, it can be found that the distribution of SRSI at the corners of IRMR has variability. So first, UA appears at the corners of IRMR and then develops to the middle of rock wall.
According to equation (3), when the distance from the working face is different, the AUAs of the roof and floor and the two sides are calculated, as shown in Figure 8. It can be seen from Figure 8 that the AUAs of the roof and floor and the two sides first increase and then decrease with mining of the working face. e maximum values of the AUAs of the roof and floor and the two sides appear at the peak of the abutment pressure, as shown by arrow q in Figure 8, indicating that the stability of the surrounding rock is the worst at the peak of the abutment pressure. Before the peak of the abutment pressure, with the mining, the AUAs of the roof and floor and the two sides first increase slowly, as shown in part (a) in Figure 8, and then increase rapidly, as shown in part (b) in Figure 8, which shows that the stability of the surrounding rock decreases slowly and then rapidly decreases before the peak of the abutment pressure. e AUA changes of the roof and floor are greater than that of the two sides, and the AUA change of the high wall area is greater than that of the low wall, indicating that the roof and floor are more affected by mining stress than the two sides, and the high wall is more affected by mining stress than the low wall. After the peak of the abutment pressure, the AUAs of the roof and floor and the two sides decrease rapidly and are reduced to 0, as shown in part (c) in Figure 8, indicating that the stability of the surrounding rock of the roadway increased.
It can be seen from Figure 8 that the AUA curve of the low wall is below the AUA curve of the high wall. When the roadway is not affected by mining stress, the AUAs of the low wall and the high wall are 1.68 m 2 and 2.12 m 2 , respectively. At the peak of the abutment pressure, the AUAs of the low wall and the high wall are increased to 2.81 m 2 and 4.24 m 2 . is shows that the AUA of the low wall is smaller than that of the high wall; thus, the low wall is more stable. ere is little difference between the AUA curves of the instability areas of the roof and floor. When the roadway is not affected by mining stress, the AUAs of the roof and floor are 0.33 m 2 and 0.35 m 2 , respectively. At the peak of the abutment pressure, the AUAs of the roof and floor increase to 4.11 m 2 and 3.97 m 2 .
e results show that the AUA difference between the roof and floor is small, and the stabilities of the roof and floor are equal. Before the peak of the abutment pressure and being less affected by the mining stress, the AUA curve of the roof and floor is located below the AUA curve of the two sides, and the AUA of the roof and floor is less than that of the two sides, indicating that the roof and floor are more stable than the two sides. e roadway section parameters and support scheme are shown in Figure 9.

Distribution Law of SRSI in IRMR with
In this paper, only the influence of bolt resistance on the stability of the surrounding rock is considered. According to the calculation method of the working resistance of bolts with different anchorage lengths in [11,13], the variations of working resistance of end anchorage and full-length anchorage bolts during the mining are obtained. By changing  Advances in Civil Engineering the boundary conditions, the working resistance of the bolt is superposed onto the stress distribution of the surrounding rock of the roadway, and the distribution of SRSI with support conditions is obtained. e contour map of SRSI with different anchoring methods is shown in Figures 10 and  11. e distribution of SRSI at 0.3 m from the roadway surface is shown in Figure 12. Equation (3) shows calculation of the AUAs of the roof and floor and the two sides under the conditions of the end anchoring bolt support and the fulllength anchoring bolt support when the distances from the working face are different, as shown in Figure 13.
Compared with Figures 5, 10, and 11, it can be seen that the SRSI of the shallow area is greater than that of the nosupport area when there is support and that of the full-length anchoring bolt support is greater than that of the end anchoring bolt support, indicating that both the end anchoring bolt support and the full-length anchoring bolt support can improve the stability of the surrounding rock. Moreover, the control effect of the full-length anchoring bolt support on the surrounding rock is stronger than that of the end anchoring bolt support. Compared with Figures 12(a) and  12(b), it can be clearly observed that, at the corner, the SRSI is higher with the full-length anchoring bolt support than that with the end anchoring bolt support, as shown by ar- rows a, b, c, d, e, and f and a′, b ′ , c ′ , d ′ , e ′ , and f ′ in Figure 12.
e SRSI in the middle of the rock wall is approximately the same as that when adopting the full-length anchoring or end anchoring bolt support, as shown by arrows g, h, and l and g ′ , h ′ , and l ′ in Figure 12. is shows that the full-length anchoring bolt support is better than the end anchoring bolt support in controlling the surrounding rock deformation at the corners of the roadway. Meanwhile, the full-length anchoring bolt support or the end anchoring bolt support has the same control effect on the middle of the rock wall. e main reason for this is that the working resistance of the bolt with full-length anchoring bolt support at the corners is higher than that of the bolt with end anchoring bolt support, while the working resistance of the bolt at the middle of the rock wall is similar to the roadway with end anchoring or full-length anchoring bolt support.
It can be seen from Figure 13 that, under the conditions of end anchoring bolt support, full-length anchoring bolt support, or no support, the AUAs of the roof and floor and the two sides tend to be the same as that with the mining. Before the peak of the abutment pressure and being less affected by the mining stress, the AUA curve under the condition of the full-length anchoring bolt support is located below the AUA curve under the condition of the end anchoring bolt support, as shown in part (a) in Figure 12. is shows that the AUA with the full-length anchoring bolt support is significantly smaller than that of the end anchoring bolt support, and the supporting effect of the fulllength anchoring bolt support is better than that of the end anchoring bolt support. At the peak of the abutment pressure and after the peak of the abutment pressure, the AUA curves with the full-length anchoring bolt support or the end anchoring bolt are basically coincident, as shown in part (b) in Figure 12. is shows that, at the peak of the    Figure 14 that the deformation curve in the full-length anchoring bolt support zone is below the deformation curve in the end anchoring bolt support zone. e deformation in the fulllength anchoring bolt support zone is less than that in the end anchoring bolt support zone, and the full-length anchoring bolt support effect is better than that of the end anchoring bolt support. According to Figures 13 and 14,   when the stations in the end anchoring bolt support zone are about 200 m away from the working face, the deformation begins to increase significantly, as shown by arrows a, b, c, and d in Figures 13 and 14. However, the displacement in the full-length anchoring bolt support zone increases significantly when stations are about 150 m away from the working face, as shown by arrows e, f, g, and h in Figures 13 and 14.
At the peak of the abutment pressure, the displacement rate of the full-length anchoring bolt support zone is obviously lower than that of the full-length anchoring bolt support zone, as shown by sections i and j in Figure 15. e main reason for this is that, before the peak of the abutment pressure, compared with that of the end anchoring bolt support, the full-length anchoring bolt support reduces the AUA to a greater extent, as shown in part (a) in Figures 13; thus, the surrounding rock is more stable when adopting    full-length anchoring bolt support. Near and after the peak of the abutment pressure, the displacement rates of the surrounding rock in the full-length anchoring bolt support zone and the end anchoring bolt support zone are slightly different as shown in section k and l in Figure 15. e main reason for this is that, near and after the peak of the abutment pressure, there is little difference in the AUAs of the surrounding rock under the conditions of the full-length anchoring bolt support and the end anchoring bolt support. e full-length anchoring bolt support of this part is equivalent to the end anchoring bolt support. In conclusion, it can be concluded that the full-length anchoring bolt support effect is significantly better than the end anchoring bolt support, and the theoretical analysis results are verified.

Conclusions
In order to study the stability characteristics of the surrounding rock and the control effect of bolt support with different anchoring methods under the influence of mining stress, the SRSI is defined and solved based on the rock strength criterion and the stress distribution.
e Mohr-Coulomb strength criterion and a complex function method are used to solve the analytical solution of the SRSI. Combined with the actual conditions, through the analysis of the changes of the SRSIs and the AUAs under the conditions of no support, the end anchoring bolt support, and the full-length anchoring bolt support, the evolution law of the stability characteristics with the mining is obtained. Combined with the observation data of the deformation, the supporting effects of the full-length and the end anchoring bolt support are analyzed. e conclusion is as follows: (1) Based on the rock strength criterion and stress distribution of the surrounding rock, the SRSI is defined, which is a more reasonable and effective method with which to analyze the stability characteristics of the surrounding rock. Based on the Mohr-Coulomb strength criterion and a complex function method, the analytical solution of the SRSI of IRMR is obtained. e SRSI can comprehensively reflect the influence of the shape of the roadway section, the stress conditions, and the strength parameters on the stability characteristics. (2) Under the condition of no support and without the influence of mining stress, the AUAs of low wall and high wall are 1.68 m 2 and 2.12 m 2 , respectively. us, the AUA of the high wall is larger than that of the low wall, and the low wall is more stable. e AUAs of the roof and floor are 0.33 m 2 and 0.35 m 2 . Compared with that of the two sides, the roof and floor are more stable. Before the peak of the abutment pressure, the SRSI decreases, the AUA increases, and the stability decreases as a whole. After the peak of the abutment pressure, the AUA decreases to 0 and the stability of the surrounding rock increases.
(3) e theoretical analysis shows that, compared with the end anchoring bolt support, the full-length anchoring bolt support method reduces the AUA to a greater extent. e surrounding rock in the end anchoring zone and the full-length anchoring zone began to deform significantly at 200 m and 150 m from the working face, respectively. e control effect of the full-length anchoring bolt support is better than that of the end anchoring bolt support and verifies the rationality of SRSI to describe the instability characteristics.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.  14 Advances in Civil Engineering