Orthogonal Numerical Analysis of Deformation and Failure Characteristics of Deep Roadway in Coal Mines: A Case Study

: With the development of deep, underground coal mines in China, the failure mechanism of the rocks surrounding roadways is becoming increasingly complicated and the deformational control is also signiﬁcantly difﬁcult. In this study, based on the temporal and spatial deformational distribution of the deep roadway area in the 2233 working face of Fuxin Hengda Coal Mine, factors affecting the deformation and failure mechanism of deep-buried roadways, such as cohesion ( c ), tensile strength ( σ t ), internal friction angle ( ϕ ), vertical ground stress ( p ), and the horizontal-to-vertical stress ratio ( λ ), were analyzed using orthogonal numerical experiments. The stress and electromagnetic radiation monitoring data were used to locate areas of highly concentrated deformation in the roadway and surrounding rocks. The results show that the order of the degree of inﬂuence of the surrounding rock and geometric parameters on the deformation of the deep-buried roadway is ϕ > p > σ t > λ > c . The vertical stress of the roof and the horizontal stress of the two sides are negatively correlated with the tensile strength and horizontal-to-vertical stress ratio, respectively, and mainly shear failure occurs in the area. The higher the level of the surrounding rock, the more serious is the deterioration and deformation. The electromagnetic radiation reﬂects the distribution range of the high-stress concentration area and strength deterioration area. The test results accurately describe the deformation–deterioration-failure laws of rock surrounding deep-buried roadways inﬂuenced by different factors. The results are of great signiﬁcance for analyzing the deformation and failure characteristics of rocks surrounding roadways, preventing rockburst, and supporting the parameter optimization of roadways.


Introduction
Coal will continue to be one of the main sources of energy in China for a long time and it has been a major strategic deployment of the country over the past few centuries. It is an inevitable trend; considering the growth in Chinese industrial development, the demand for resources from the deep pits is bound to rise. According to statistics on the rate of increase in mining depth across the country, the old mines in the eastern provinces are developing at a mining speed of 50-125 m every five years. By the end of 2017, there were more than 50 mines with a buried depth of more than 1 km. The geological environment in the shallow buried roadway is simple in stope structures, low in ground stress, and less affected by geological faults [1][2][3][4][5]. Deep mining faces many problems, such as significantly increased in situ stress, complex geological environments, and the coupling of multiphase rock masses. Owing to these factors, the physical and mechanical properties of deep rock masses and the instability mechanism of the surrounding rocks of the roadway are becoming increasingly complicated. As a result, it is difficult to accurately estimate the deformation of rocks surrounding the deep-buried roadway and the difficulty of controlling the structural stability of the rocks continues to increase [6][7][8][9][10][11]. They may even cause dynamic disasters, such as roof collapse, rockbursts, coal and gas outbursts, and large-scale floor subsidence, that could seriously jeopardize the safe and efficient operation of the coal mines.
Scholars have done a lot of outstanding work in the past on the deformation characterization and the temporal and spatial evolution of the fractural process, key parameters for the deformation control of the surrounding rocks and related supporting technologies of the roadway. Zhang et al. (2018) established the mechanics and structural model of the roadway and proposed that the roof collapse deformation is caused by buckling when the roof stratum reaches the failure limit and, as a consequence, sliding shear failure is likely to occur on the sidewalls; a decrease in bearing capacity is the root cause of floor undulation [12]. Khalymendyk et al. (2021) analyzed the relationship between the mutual restriction of the "support-rock mass" system and the particularity of the deformation mechanism of the roadway, and a linear relationship between the subsidence of the roof and the uplift of the bottom drum was established [9]. Pan et al. (2020) established a mechanical model of the "coal wall-hydraulic support-gangue" support system from which the cyclical failure span of the stope and the modified expression of the impact kinetic energy were derived [13]. An innovative control technology for fracturing through vertical walls was also proposed. Zhao et al. (2021) found that high in situ stress in deep mines, local stress concentrations in the fault, shear-slip failure, and insufficient support strength, etc., are the main reasons for the deformation and failure of the surrounding rock of the roadway crossing the fault [14]. Ding et al. (2021) found that a weak layered roadway weakened the stress transfer from the deep part of the surrounding rock to the coal bank, promoted the stress transfer from the roof to the coal bank, and intensified the shear failure and horizontal slip of the surrounding rock [15]. Sun et al. (2020) proposed a new coal-roadway pre-reinforcement method with jet grouting to strengthen the surrounding soft-coal wall and form a stable support structure [16]. Jia et al. (2021) found that the maximum depth of the plastic zone of the roadway usually faces the roof, and plastic compression-shear failure occurs when the plastic failure-affected zone contains a weak layer [17]. Until now, the mechanical properties of deep rock masses and the evolution characteristics of the original rock stress fields have not been determined with accuracy. Similarly, the distribution of stress fields in the rock mass around the deep-buried roadway and the evolution characteristics of rock movement under the influence of mining could not be determined with a fair degree of accuracy, either through mathematical calculations or a study of the mechanics involved.
Therefore, an in-depth study of the deformation characteristics of the surrounding rocks and influencing factors in deep-buried roadways is useful for controlling deformation, degradation, and destruction. It can help reveal the mechanism of instability and optimize the support design parameters of deep-buried roadways. Considering the great significance of ensuring and maintaining safe production in deep-buried mines, this study has observed the 2233 working face of the Fuxin Hengda Coal Mine as the engineering research background and adopts the method of combining simulation analysis with orthogonal experiment. Five mechanical factors on the influence of the deformation and deterioration mechanism were analyzed, including the cohesion (c), tensile strength (σ t ), internal friction angle (ϕ), vertical stress (p), and horizontal-to-vertical stress ratio (λ) of the surrounding rock of the deep-buried roadway. The stress, support resistance, and electromagnetic radiation monitoring data collected in the roadway before and after the occurrence of rockbursts in the 2233 working face of the Fuxin Hengda Coal Mine were analyzed. Furthermore, the time and space characteristics of the surrounding rock-destruction of the roadway were corroborated. The research results can provide a theoretical reference for the optimization of the deformation control, reveal the failure mechanism of the surrounding rock of deep-buried roadways, and lead to the optimization of the stability control parameters. Additionally, it could shed new light on the prevention and control of rockburst disasters in roadways.

Failure Mechanism of Surrounding Rock of Roadway
The results of roadway deterioration and instability in the Fuxin Hengda Coal Mine are shown in Figure 1. The elastic potential energy originally stored in the surrounding rock that is naturally formed by geological processes is in a state of dynamic equilibrium [18][19][20]. However, external engineering elements, such as pressure relief owing to blasting, largearea collapse of the roof, etc., disrupt the original stress balance and the redistribution of the stress in the rock mass. This leads to the accumulation and development of the two defects of coal block yield deformation, and the accumulated elastic energy is released sharply, causing the surrounding rock to rapidly transition from the deterioration stage to instability failure (Figure 1a,b). In the evolution process of rock mass deterioration, the block primitives with a certain strength evolve to the main structure with a high degree of damage and the initial strength of the actual material is significantly reduced [21][22][23]. Among them, the Mohr-Coulomb destruction criterion is a constitutive model with strong applicability and the widest application range in geotechnical structure engineering. The Mohr-Coulomb criterion is used to describe the strength characteristics of the plastic zone of the surrounding rock, and its form is: where f s is shear failure criterion, f t is tensile failure criterion, σ 1 is the maximum axial principal stress of the coal rock unit body, σ 3 is the lateral pressure, σ t is the tensile strength, is the internal friction angle, and c is the cohesion.
the surrounding rock of deep-buried roadways, and lead to the optimization of the stabil ity control parameters. Additionally, it could shed new light on the prevention and contro of rockburst disasters in roadways.

Failure Mechanism of Surrounding Rock of Roadway
The results of roadway deterioration and instability in the Fuxin Hengda Coal Min are shown in Figure 1. The elastic potential energy originally stored in the surrounding rock that is naturally formed by geological processes is in a state of dynamic equilibrium [18][19][20]. However, external engineering elements, such as pressure relief owing to blast ing, large-area collapse of the roof, etc., disrupt the original stress balance and the redis tribution of the stress in the rock mass. This leads to the accumulation and developmen of the two defects of coal block yield deformation, and the accumulated elastic energy i released sharply, causing the surrounding rock to rapidly transition from the deteriora tion stage to instability failure (Figure 1a,b). In the evolution process of rock mass deteri oration, the block primitives with a certain strength evolve to the main structure with high degree of damage and the initial strength of the actual material is significantly re duced [21][22][23]. Among them, the Mohr-Coulomb destruction criterion is a constitutiv model with strong applicability and the widest application range in geotechnical structur engineering. The Mohr-Coulomb criterion is used to describe the strength characteristic of the plastic zone of the surrounding rock, and its form is: where is shear failure criterion, is tensile failure criterion, 1 is the maximum axia principal stress of the coal rock unit body, 3 is the lateral pressure, t is the tensil , Φ is the internal friction angle, and is the cohesion. The surrounding rock of the underground roadway is in a three-dimensional stres state [24][25][26][27]. When the local stress is large, the lateral pressure increases beyond the coa rock tensile strength threshold; the stress distribution of the surrounding rock satisfie Equation (2) but does not meet the failure conditions of Equation (1) at a point within th influence range of the spatial constraint effect. Tensile instability deformation occurs in side the rock mass. Although the surface of the surrounding rock will not be destroyed immediately under this condition, dynamic disasters, such as rockbursts, could occur. Th macroscopic appearance is a constant isokinetic creep, and as time progresses, the yield Where the annular region of radius r o ∼ R S is plastic flow zone, R S ∼ R P is plastic softening zone, R P ∼ R c is plastic hardening zone.
The surrounding rock of the underground roadway is in a three-dimensional stress state [24][25][26][27]. When the local stress is large, the lateral pressure increases beyond the coal rock tensile strength threshold; the stress distribution of the surrounding rock satisfies Equation (2) but does not meet the failure conditions of Equation (1) at a point within the influence range of the spatial constraint effect. Tensile instability deformation occurs inside the rock mass. Although the surface of the surrounding rock will not be destroyed immediately under this condition, dynamic disasters, such as rockbursts, could occur. The macroscopic appearance is a constant isokinetic creep, and as time progresses, the yield of the surrounding rock surface causes changes in the properties of the coal rock mass because of the creep. The strength of the coal rock is controlled by parameters such as the internal friction angle and cohesion. Under this condition, the internal point stress of the surrounding rock within a certain depth reaches the Equation (1) requirement, leading to the accelerated failure of the surrounding rock and the expansion of the deformation failure area to the deep part of the surrounding rock. Consequently, this causes the surrounding rock to produce a loose zone and completely collapse, resulting in roadway failure. This phenomenon is manifested in the destruction of the roadway in Fuxin, either as the surface deformation of the surrounding rock-roof collapse or the bottom heave severe compression failure of the roadway; eventually, the status shown in Figure 1 is attained. From the above analysis it can be seen that the main factors affecting the deformation of the surrounding rock include the surrounding rock cohesion (c), tensile strength (σ t ), internal friction angle (ϕ), vertical ground stress (p), and horizontal-to-vertical stress ratio (λ).

Mine Geology and Numerical Model
Six layers of coal are available for mining at the Fuxin Hengda Coal Mine that is located in Fuxin City, Liaoning Province, China. The 2233 working face in the 126-mining area is located in the Sunjiawan coal seam, whose average thickness is 4.8 m, inclination angle is 4 • , and burial depth is approximately 900 m. The designed mining area of the working face was 700 × 210 m. The structure of the coal seam is complex, but the nature of the surrounding rock of the working face roadway is good. The main components of the working face roof consist of thick, hard, and coarse sandstone, whereas the floor consists of siltstone and fine sandstone. The positions of the slope of the working face and the geological histogram are shown in Figures 2 and 3.
of the surrounding rock surface causes changes in the properties of the coal rock mas because of the creep. The strength of the coal rock is controlled by parameters such as th internal friction angle and cohesion. Under this condition, the internal point stress of the surrounding rock within a certain depth reaches the Equation (1) requirement, leading to the accelerated failure of the surrounding rock and the expansion of the deformation fail ure area to the deep part of the surrounding rock. Consequently, this causes the surround ing rock to produce a loose zone and completely collapse, resulting in roadway failure This phenomenon is manifested in the destruction of the roadway in Fuxin, either as th surface deformation of the surrounding rock-roof collapse or the bottom heave severe compression failure of the roadway; eventually, the status shown in Figure 1 is attained From the above analysis it can be seen that the main factors affecting the deformation o the surrounding rock include the surrounding rock cohesion ( ), tensile strength (σ t ), in ternal friction angle ( ), vertical ground stress ( ), and horizontal-to-vertical stress ratio ( ).

Mine Geology and Numerical Model
Six layers of coal are available for mining at the Fuxin Hengda Coal Mine that is lo cated in Fuxin City, Liaoning Province, China. The 2233 working face in the 126-mining area is located in the Sunjiawan coal seam, whose average thickness is 4.8 m, inclination angle is 4°, and burial depth is approximately 900 m. The designed mining area of th working face was 700 × 210 m. The structure of the coal seam is complex, but the nature of the surrounding rock of the working face roadway is good. The main components o the working face roof consist of thick, hard, and coarse sandstone, whereas the floor con sists of siltstone and fine sandstone. The positions of the slope of the working face and the geological histogram are shown in Figures 2 and 3. To better observe the deformation of the roadway combined with the calculation ac curacy requirements of the simulation, a simplified numerical model sized 100 × 50 × 100 m with a composite structure of coal and rock in nine layers corresponding to the geolog ical conditions of the Fuxin Hengda Coal Mine is considered. The roadway, whose mode is shown in Figure 4, is located in the middle of the coal seam; the lower half of the road way section is rectangular with a length of 4 m and width of 2 m, while the upper part i designed as a semicircular arch with a radius of 2 m. A uniform load of 30 MPa was ap plied to the top of the model to compensate for the overlying coal and rock layer's own weight stress, and the horizontal-to-vertical stress ratio was set to 1, i.e., 30 MPa stres around the model. The surrounding boundary uses a fixed constraint method to limit the displacement of the boundary of the model and the Mohr-Coulomb constitutive mode To better observe the deformation of the roadway combined with the calculation accuracy requirements of the simulation, a simplified numerical model sized 100 × 50 × 100 m with a composite structure of coal and rock in nine layers corresponding to the geological conditions of the Fuxin Hengda Coal Mine is considered. The roadway, whose model is shown in Figure 4, is located in the middle of the coal seam; the lower half of the roadway section is rectangular with a length of 4 m and width of 2 m, while the upper part is designed as a semicircular arch with a radius of 2 m. A uniform load of 30 MPa was applied to the top of the model to compensate for the overlying coal and rock layer's own weight stress, and the horizontal-to-vertical stress ratio was set to 1, i.e., 30 MPa stress around the model. The surrounding boundary uses a fixed constraint method to limit the displacement of the boundary of the model and the Mohr-Coulomb constitutive model has been selected. The heights of the coal and rock layers in this model and their related rock and soil mechanical parameters are listed in Table 1. Based on the geological survey report of the Fuxin Hengda Coal Mine, under the premise of ensuring accuracy and computing efficiency, we established a roadway model to describe the real geological form of the mine.
Minerals 2022, 12,185 has been selected. The heights of the coal and rock layers in this model and their r rock and soil mechanical parameters are listed in Table 1. Based on the geological s report of the Fuxin Hengda Coal Mine, under the premise of ensuring accuracy and puting efficiency, we established a roadway model to describe the real geological fo the mine.

Determine Evaluation Indicators
The analytic hierarchy process refers to the establishment of an evaluation system by linking the elements with the overall goal through quantitative standards, thus, providing a reasonable decision-making framework for the problem at hand. On the basis of this

Determine Evaluation Indicators
The analytic hierarchy process refers to the establishment of an evaluation system by linking the elements with the overall goal through quantitative standards, thus, providing a reasonable decision-making framework for the problem at hand. On the basis of this framework, using the method which solves the eigenvectors of the judgment matrix and passing qualitative indicators, the fuzzy quantification method ranks the relative pros and cons of the levels [28][29][30][31]. The analysis shows that the roadway deformation under the influence of mining includes roof collapse, subsidence, the bottom heave of the floor, and the displacement of the sides. The deformation of the roadway can be accurately described through an analysis of these deformations.
This study uses roadway deformation D to characterize the influence of the physical and mechanical parameters of coal on the deformation and failure of the roadway. The more serious the degree of deformation and failure of the surrounding rock, the larger the roadway deformation.
where D i : D 1 is the maximum deformation of the roadway roof, D 2 is the maximum deformation of the roadway sides, D 3 is the maximum deformation of the roadway floor, and ω i is the corresponding deformation weight coefficient of D i .
The deformation weight coefficient ω i can be determined using the analytic hierarchy process. Based on long-term process practice and production experience, using the 1-9 scale matrix method [32][33][34], a third-order judgment matrix of the degree of deformation and weight in the roadway can be defined as follows: Through calculation, the matrix P A−B has a maximum eigenvalue λ max = 3.009, then the eigenvector (0.848, 0.428, 0.3113) corresponding to the maximum eigenvalue is unitized and the following weight vector can be obtained: where ω A−B is the corresponding weights for the deformation of the roof, two sides, and floor, respectively. Then, Equations (6) and (7) are used to test the consistency ratio of the weight vector of the single layer and determine whether the test is passed according to the ratio value: where CR is the eigenvalue of the single-layer weight vector ratio used to analyze and determine the construction matrix. When the judgment result CR is less than 0.1, it can be considered that the single-layer weight vector ranking result of the constructed matrix has satisfactory statistical consistency; otherwise, the construction matrix must be re-analyzed and corrected again. The consistency test eigenvalues of the third-order judgment matrix of the roadway deformation are used to judge whether the degree of inconsistency of the matrix is acceptable; when CI, the characteristic value, is equal to 0, it is considered that the constructed third-order roadway deformation matrix is completely consistent and can accurately reflect the roadway deformation characteristics; when CI is close to zero, it is considered that the matrix shows a high degree of consistency and can be used to study roadway deformation analysis; when CI deviates from 0, the inconsistency is more serious and it cannot accurately reflect the deformation characteristics of the roadway. The values of RI, the random consistency parameters, are listed in Table 2, where n is the number of factors. For the construction of the third-order judgment matrix of roadway deformation, n = 3, RI is the average random consistency index and RI = 0.52 corresponding to the number of factors n = 3. It can be seen from the calculation results that the value of CI is close to zero and that of CR is much less than 0.1, and based on these results the constructed roadway deformation matrix is considered to have satisfactory consistency. When evaluating the influence of different factors on the deformation of deep-buried roadways, the corresponding weights for the deformation of the roof, two sides, and floor were 53.4%, 27%, and 19.6%, respectively.

Orthogonal Experimental Design
The orthogonal experimental design is an experimental design theory and technique proposed by studying the multifactor and multilevel [35][36][37]. Consequent to analyses of the related research of scholars, this study selected five physical, mechanical, and geological factors concerning rocks, viz., tensile strength σ t , cohesion c, internal friction angle ϕ, vertical ground stress p, and horizontal-to-vertical stress ratio λ to study the influence of the deformation and degradation mechanism of deep-buried roadways. The internal friction angle ϕ actual effect in the Mohr-Coulomb equation is determined by the tangent results. The internal friction angle is used instead of the friction coefficient, because it is the internal friction angle input in the FLAC3D 6.0 software, (Itasca Regions, MN, USA), instead of the friction coefficient which facilitates understanding and application [38][39][40][41]. Each factor took four parameter levels, selected L 16 4 5 on the orthogonal table, and conducted related experiments. The grouping value of each level and factor is determined based on the analysis of a large amount of mine data. The specific experimental plan is presented in Table 3. The selection of the parameters mainly considers the equal difference increment-  Table 3, 16 groups of FLAC3D orthogonal experiment schemes were designed. The initial model was used as the prototype to carry out the simulation calculation of the roadway under different parameter combinations and the data results of the roadway surrounding rock deformation under different conditions were obtained [42][43][44][45]. The roadway deformation index results are presented in Table 4. The experimental results were calculated and processed using normalization. After normalization, the influence of the multiplier on the deterioration and deformation analysis of the roadway can be reduced, and an influence effect of different parameters on the deterioration and deformation of the roadway can be obtained more objectively and accurately.

Range Analysis
Range analysis refers to the method of extracting the mean and the range corresponding to each level of various factors from the results of orthogonal experiments and analyzing the changes in various factors that affect the target fluctuation interval [46][47][48]. The range characterizes the range of influence of the corresponding factors on the evaluation index of the orthogonal experiment (the amount of roadway deformation). The range analysis and calculation of the experimental statistical results are shown in Table 4.
Through range analysis, it can be determined that the roadway deformation reaches the maximum value of 540.80 mm when the internal friction angle ϕ is 20 • and the minimum value of 86.55 mm at 50 • internal friction angle; the range is the largest at 454.25 mm when the vertical ground stress p is 15 MPa, the minimum roadway deformation is 50.74 mm, and when it is 45 MPa, the maximum value is 489.11 mm, with a large range of 438.37 mm; when the tensile strength σ t is 0.2 MPa, the minimum roadway deformation is 206.96 mm, and when it is 0.2 MPa, the maximum is 603.64 mm, with a large range of 406.68 mm. Through the sorting of the range, it can be observed that different physical and mechanical properties of the rock have varying effects on the deformation and failure of the deep-buried roadways, and the degree of influence in descending order is: ϕ > p > σ t > λ > c. This shows that within a certain range, ϕ has the greatest impact on the deformation of the deep-buried roadways, followed by p and σ t , whereas λ and c have a relatively minor influence.
Based on the average value obtained from the range analysis in Table 5, the corresponding intuitive analysis chart that describes the deformation trend of the roadway corresponding to changes in different factors can be obtained. It can be seen from Figure 5 that the tensile strength σ t , cohesion c, and internal friction angle have a negative influence on the roadway deformation. When the internal friction angle ϕ is less than 30 • , the slope of 29.7 mm/ • is higher and the internal friction angle has a significant influence on the deformation of the roadway. In contrast, when the internal friction angle is greater than 30 • , the slope reduces and its influence on the deformation of the roadway weakens: the reduction speed gradually slows and the degree of roadway deformation is lower. When the tensile strength σ t increases from 0.3 to 0.4 MPa, and the deformation of the roadway drops from 502.19 to 225.55 mm, it indicates a decrease of 55.08%, which is substantial. However, with a further increase in the tensile strength, the downward trend becomes lower. The cohesion c causes a decrease in the deformation, particularly in the range between 2 and 3.0 MPa, and the change in the deformation of the roadway is significant. The vertical ground stress and horizontal-to-vertical stress ratios are positively correlated with the deformation of the roadway. The vertical ground stress p causes an increase in the road deformation and this relationship is near-linear in fashion: the curve has a slope of 14.62 mm/MPa. When the horizontal-to-vertical stress ratio λ increases from 1 to 2 it has a significant impact on the deformation of the roadway: the curve is nearly horizontal when it is between 2 and 2.5.

Results of Variance Analysis
The method of variance analysis can obtain a quantitative estimate of the significance of the decision-making effect by studying the contributions of different types of decisions to the target [49][50][51]. This is a common method for analyzing the results of orthogonal experiments. The variance analysis table of the roadway deformation is obtained by calculating the method of decomposing the sum of squared deviations, as shown in Table 6. The five parameters have different significances in the influence of the deformation of the deep-buried roadway. As shown in Table 6, the F value of the modified model is 436.237 and the significance is 0.000, which indicates that the modified model is very significant and can accurately represent the characteristics of the roadway deformation. According to the comparison and analysis of the F value, the order of the deformation contribution degree of the various influencing factors in the roadway deformation can be calculated as: > > > > . Among them, the F value of the internal friction angle and vertical ground stress are higher, i.e., 665.806 and 608.174, respectively, whereas the F value of cohesion is lower-only 140.103. Cohesion c is more effective when rock is under the low-stress condition and at a shallow depth, but cohesion c has little influence on the deformation of the roadway under the high-stress state, indicated through the range and variance analysis of the results of the orthogonal simulation experiment. In

Results of Variance Analysis
The method of variance analysis can obtain a quantitative estimate of the significance of the decision-making effect by studying the contributions of different types of decisions to the target [49][50][51]. This is a common method for analyzing the results of orthogonal experiments. The variance analysis table of the roadway deformation is obtained by calculating the method of decomposing the sum of squared deviations, as shown in Table 6. The five parameters have different significances in the influence of the deformation of the deep-buried roadway. As shown in Table 6, the F value of the modified model is 436.237 and the significance is 0.000, which indicates that the modified model is very significant and can accurately represent the characteristics of the roadway deformation. According to the comparison and analysis of the F value, the order of the deformation contribution degree of the various influencing factors in the roadway deformation can be calculated as: ϕ > p > σ t > λ > c. Among them, the F value of the internal friction angle ϕ and vertical ground stress σ t are higher, i.e., 665.806 and 608.174, respectively, whereas the F value of cohesion is lower-only 140.103. Cohesion c is more effective when rock is under the low-stress condition and at a shallow depth, but cohesion c has little influence on the deformation of the roadway under the high-stress state, indicated through the range and variance analysis of the results of the orthogonal simulation experiment. In summary, the internal friction angle ϕ and the strength of the local vertical ground stress p can significantly affect the degree of deformation of a deep-buried roadway. The internal friction angle ϕ represents the sliding friction and occlusal friction caused by the interaction between fine particles inside the coal rock, and under high-stress conditions, the internal friction angle dominates the shear strength of the surrounding rock. The internal friction angle simulated in this paper increases linearly, while the corresponding tangent value friction coefficient changes nonlinearly in the Mohr-Coulomb equation. The larger the internal friction angle ϕ, the greater the shear strength of coal and rock, the stronger the ability of the surrounding rock to resist shear sliding, the smaller the volume and distribution of the plastic zone of shear failure, and the greater the obstacle to the deformation of the surrounding rock of the roadway. Therefore, there is a negative correlation between the road deformation and the internal friction angle ϕ. The vertical ground stress is determined by the buried depth of the roadway. The greater the vertical ground stress, the greater the vertical load on the roadway, and the greater the promotion of roadway deformation. Therefore, roadway deformation and vertical ground stress are positively correlated. The purpose of this paper is mainly to perform an orthogonal analysis, considering the influence on the deterioration and deformation of the surrounding rock of the roadway in the process of depth change. Currently, there are many coal and metal mines with burial depths ranging from 500 to 1500 m, and more and more coal and metal mines will reach this depth in the near-future. In addition, the development trend of China's mining industry is to go deeper. The purpose is to better and more accurately predict the deformation characteristics of the surrounding rock in the deep roadway, and to explore the instability and failure characteristics of the roadway under the state of deep burial and high in situ stress.

Characteristics of Roadway Stress
Stress change is an inherent dynamic cause of the deformation and failure of deepburied roadways. Analyzing the stress trend under different factors is helpful for studying the characteristics of roadway deformation [52][53][54]. In this study, the maximum vertical stress on the roof and floor, and the maximum horizontal stress on the two sides of the roadway, were analyzed. The results are shown in Figure 6. summary, the internal friction angle and the strength of the local vertical ground stress p can significantly affect the degree of deformation of a deep-buried roadway. The internal friction angle φ represents the sliding friction and occlusal friction caused by the interaction between fine particles inside the coal rock, and under high-stress conditions, the internal friction angle dominates the shear strength of the surrounding rock. The internal friction angle simulated in this paper increases linearly, while the corresponding tangent value friction coefficient changes nonlinearly in the Mohr-Coulomb equation. The larger the internal friction angle φ, the greater the shear strength of coal and rock, the stronger the ability of the surrounding rock to resist shear sliding, the smaller the volume and distribution of the plastic zone of shear failure, and the greater the obstacle to the deformation of the surrounding rock of the roadway. Therefore, there is a negative correlation between the road deformation and the internal friction angle φ. The vertical ground stress is determined by the buried depth of the roadway. The greater the vertical ground stress, the greater the vertical load on the roadway, and the greater the promotion of roadway deformation. Therefore, roadway deformation and vertical ground stress are positively correlated. The purpose of this paper is mainly to perform an orthogonal analysis, considering the influence on the deterioration and deformation of the surrounding rock of the roadway in the process of depth change. Currently, there are many coal and metal mines with burial depths ranging from 500 to 1500 m, and more and more coal and metal mines will reach this depth in the near-future. In addition, the development trend of China's mining industry is to go deeper. The purpose is to better and more accurately predict the deformation characteristics of the surrounding rock in the deep roadway, and to explore the instability and failure characteristics of the roadway under the state of deep burial and high in situ stress.

Characteristics of Roadway Stress
Stress change is an inherent dynamic cause of the deformation and failure of deepburied roadways. Analyzing the stress trend under different factors is helpful for studying the characteristics of roadway deformation [52][53][54]. In this study, the maximum vertical stress on the roof and floor, and the maximum horizontal stress on the two sides of the roadway, were analyzed. The results are shown in Figure 6.  For the vertical stress on the floor, the internal friction angle ϕ responds more obviously to the floor stress: the floor stress increases with an increase in the internal friction angle. The changing trend of the other factors presents a horizontal straight line with no obvious change. The tensile strength σ t has the greatest impact on the roof for the vertical stress of the roadway roof. The tensile strength σ t increases from 0.2 to 0.3 MPa and the stress decreases from 18.89 to 8.02 MPa-a decrease of 57.5%. The cohesion c and internal friction angle ϕ have similar influence trends on roof stress. The roof stress presents an "N-shaped" change trend that first increases, then decreases, and finally increases. The horizontal-tovertical stress ratio λ has a significant impact on the vertical stress of the roof of the roadway and has the greatest impact on the horizontal stress of the two sides: the two curves are relatively consistent. When the horizontal-to-vertical stress ratio λ changes from 1.5 to 2, the stress on the two sides decreases from 15.58 to 5.82 MPa-a decrease of 62.6%. When the cohesion c and the vertical ground stress p increase, the stresses of the two sides show an inverted-V-shaped change trend: it first increases and then decreases. In comparison, when the cohesion c changes, the degree and range of the stress on the two sides of the roadway are relatively large and the internal friction angle ϕ has little effect on the stress value of the two sides of the roadway. As the internal friction angle ϕ increases from 20 to 50 • , the stress on the two sides of the roadway is maintained at approximately 11 MPa.

Analysis of Roadway Plastic Zone during Mining
The plastic zone of the roadway refers to the yield zone, where the load exceeds the ultimate bearing capacity of coal and rock, causing a local irreversible deformation [55][56][57]. On completion of the calculation with the FLAC3D software, the unique stress state of the element could be determined. When the stress state reaches the strength envelope, it is considered that the element enters the plastic state and a plastic zone is generated in the model. Otherwise, no plastic zone is generated. With the plastic potential surface used to represent the shear plastic flow (g s ) and tensile plastic flow (g t ), it is determined whether the unit is in the shear plastic zone or the tensile failure zone.
To study the deformation and failure of the roadway during the mining of the straight wall semicircular arch roadway, the plastic zone cloud diagrams of the five stages, viz., mining to 10, 20, 30, 40, and 50 m, were selected for analysis, as shown in Figure 7. The plastic zone is mainly centered on the roadway and extends in an "X" shape. The deformation of the roadway is mainly due to the influence of shear damage. Among them, the shoulders and floor of the two sides of the roadway are sensitive areas that are extremely vulnerable to damage, and the plastic zone first appears and concentrates here. The two shoulders of the straight wall semicircular arch roadway first exhibited shear failure. With the progress of excavation, the shear failure of the roof of the roadway gradually expanded, while the shear failure and tensile failure on the two sides evolved simultaneously. Finally, because the two sides of the roof are mainly subjected to shear stress, the roof sinks and the two sides move inward; the surrounding rock withstands shear failure, while the floor is mainly subjected to a horizontal pulling force, the bottom heave swells, and the floor undergoes tensile failure.

Comprehensive Analysis of Surrounding Rock Classification
The composition parameters of the strength of the surrounding rock were complex. The main parameters include cohesion, internal friction angle, and tensile strength. These different surrounding rock mechanical parameters do not exist in isolation, and they are interrelated and mutually restricted. The influence of various parameters on the strength of the surrounding rock was determined in a previous study [58][59][60][61][62]. Taking into consideration the findings and results of the above research in combination with the characteristics of surrounding rock grades, this article discusses the deformation of roadways under common surrounding rock grades II-IV. The relevant parameters are listed in Table 7 and the simulation results are shown in Figures 7 and 8.

Comprehensive Analysis of Surrounding Rock Classification
The composition parameters of the strength of the surrounding rock were complex. The main parameters include cohesion, internal friction angle, and tensile strength. These different surrounding rock mechanical parameters do not exist in isolation, and they are interrelated and mutually restricted. The influence of various parameters on the strength of the surrounding rock was determined in a previous study [58][59][60][61][62]. Taking into consideration the findings and results of the above research in combination with the characteristics of surrounding rock grades, this article discusses the deformation of roadways under common surrounding rock grades II-IV. The relevant parameters are listed in Table 7   It can be observed from Figure 8 that with a continuous increase of the level of the surrounding rock, the growth rate of roadway deformation gradually increases and the deformation deterioration becomes more serious. Additionally, the horizontal stress on the two sides is significantly higher than the vertical stress on the roof and both gradually increase with the grade of the surrounding rock, whereas the vertical stress on the floor does not change significantly and remains nearly constant at approximately 0.1 MPa. It can be seen from Figure 9 that the plastic zone is still dominated by shear failure and the roadway floor mainly exhibits tensile failure. Evidently, the distribution area of the plastic zone continues to expand with the grade of the surrounding rock, and the shear plastic zone continues to extend in an "X" shape. The plastic zone of the grade II surrounding rock has the smallest area and is also the most stable, whereas the plastic zone of the grade IV surrounding rock has the widest distribution area, a loose structure, and poor resistance to deformation. Furthermore, as the level of the surrounding rock increases, its stability worsens and greater is the impact on the degraded deformation of the roadway. It can be observed from Figure 8 that with a continuous increase of the level of the surrounding rock, the growth rate of roadway deformation gradually increases and the deformation deterioration becomes more serious. Additionally, the horizontal stress on the two sides is significantly higher than the vertical stress on the roof and both gradually increase with the grade of the surrounding rock, whereas the vertical stress on the floor does not change significantly and remains nearly constant at approximately 0.1 MPa. It can be seen from Figure 9 that the plastic zone is still dominated by shear failure and the roadway floor mainly exhibits tensile failure. Evidently, the distribution area of the plastic zone continues to expand with the grade of the surrounding rock, and the shear plastic zone continues to extend in an "X" shape. The plastic zone of the grade II surrounding rock has the smallest area and is also the most stable, whereas the plastic zone of the grade IV surrounding rock has the widest distribution area, a loose structure, and poor resistance to deformation. Furthermore, as the level of the surrounding rock increases, its stability worsens and greater is the impact on the degraded deformation of the roadway.

On-Site Verification
During the stopping period, the evolutionary process of mine pressure is shown as a cycle of "accumulation-transformation-release." In this process, a part of the energy is released in the form of electromagnetic waves [63][64][65]. Therefore, analyzing the distribution characteristics of the electromagnetic radiation energy in the mining process plays an important role in revealing the roadway deformation characteristics and evolution process. In this study, a KJ550 stress-monitoring system, resistance monitoring system of hydraulic support, and KBD5 portable electromagnetic radiation instrument were used to analyze the fracture characteristics of coal and rock in front of the working face, and the  It can be observed from Figure 8 that with a continuous increase of the level of the surrounding rock, the growth rate of roadway deformation gradually increases and the deformation deterioration becomes more serious. Additionally, the horizontal stress on the two sides is significantly higher than the vertical stress on the roof and both gradually increase with the grade of the surrounding rock, whereas the vertical stress on the floor does not change significantly and remains nearly constant at approximately 0.1 MPa. It can be seen from Figure 9 that the plastic zone is still dominated by shear failure and the roadway floor mainly exhibits tensile failure. Evidently, the distribution area of the plastic zone continues to expand with the grade of the surrounding rock, and the shear plastic zone continues to extend in an "X" shape. The plastic zone of the grade II surrounding rock has the smallest area and is also the most stable, whereas the plastic zone of the grade IV surrounding rock has the widest distribution area, a loose structure, and poor resistance to deformation. Furthermore, as the level of the surrounding rock increases, its stability worsens and greater is the impact on the degraded deformation of the roadway.

On-Site Verification
During the stopping period, the evolutionary process of mine pressure is shown as a cycle of "accumulation-transformation-release." In this process, a part of the energy is released in the form of electromagnetic waves [63][64][65]. Therefore, analyzing the distribution characteristics of the electromagnetic radiation energy in the mining process plays an important role in revealing the roadway deformation characteristics and evolution process. In this study, a KJ550 stress-monitoring system, resistance monitoring system of hydraulic support, and KBD5 portable electromagnetic radiation instrument were used to analyze the fracture characteristics of coal and rock in front of the working face, and the

On-Site Verification
During the stopping period, the evolutionary process of mine pressure is shown as a cycle of "accumulation-transformation-release." In this process, a part of the energy is released in the form of electromagnetic waves [63][64][65]. Therefore, analyzing the distribution characteristics of the electromagnetic radiation energy in the mining process plays an important role in revealing the roadway deformation characteristics and evolution process. In this study, a KJ550 stress-monitoring system, resistance monitoring system of hydraulic support, and KBD5 portable electromagnetic radiation instrument were used to analyze the fracture characteristics of coal and rock in front of the working face, and the data of the transportation along the channel before the occurrence of the "6 August rockburst accident" at the 2233 working face of the Fuxin Hengda Coal Mine was selected for analysis. After the rockburst accident occurred on 6 August, three impact damage sections appeared and the section with a distance of 68 m was the most severe impact damage section. The electromagnetic radiation data were received at 20 measuring points in front of the working face on 9 August, with a test interval of 10 m (Figure 9).
The stress and electromagnetic radiation intensity can reflect the local stress level. High-stress areas will significantly reduce the physical and mechanical properties of coal and rock mass, which promotes roadway damage and sudden instability [66][67][68]. It can be seen from Figure 10 that before the rockburst accident the deep coal seam stress in front of the working face increased sharply from 9.2 to 19 MPa, leading to the occurrence of the rockburst: the stress dropped linearly to approximately 3 MPa after the impact. Simultaneously, the resistance of the corresponding hydraulic support also increased sharply before dropping significantly after the impact, because it did not cause any major damage and can still play a supporting role later. Concurrently, the electromagnetic radiation intensity showed an inverted-N-shaped fluctuation trend, i.e., down-up-down, in the monitoring interval. It reached its lowest value on 5 August when the impact risk was the greatest, according to literature records: this is a typical strong shock signal law. On 6 August, the roadway experienced rockbursts and the site survey revealed that the main damage area of the roadway was 68 m outside the front of the working face because of the impact. This corresponds to the relatively high value of the electromagnetic radiation in the early stage, which is a good precursor phenomenon, and the other two damaged areas correspond to relatively low points of electromagnetic radiation. After the impact, the stress dropped sharply and the electromagnetic radiation intensity presented a relatively stable wave-like distribution in the direction of the working surface. There was no obvious high-value concentration area, indicating that after the impact the coal in the broken zone is completely unstable and the plastic zone and elastic zone are again affected by the supporting stress. Accurate identification of the high-stress zone of the roadway is of great significance for strengthening the support to prevent roadway deformation, as well as roadway pressure shock. and rock mass, which promotes roadway damage and sudden instability [66][67][68]. It can be seen from Figure 10 that before the rockburst accident the deep coal seam stress in front of the working face increased sharply from 9.2 to 19 MPa, leading to the occurrence of the rockburst: the stress dropped linearly to approximately 3 MPa after the impact. Simultaneously, the resistance of the corresponding hydraulic support also increased sharply before dropping significantly after the impact, because it did not cause any major damage and can still play a supporting role later. Concurrently, the electromagnetic radiation intensity showed an inverted-N-shaped fluctuation trend, i.e., down-up-down, in the monitoring interval. It reached its lowest value on 5 August when the impact risk was the greatest, according to literature records: this is a typical strong shock signal law. On 6 August, the roadway experienced rockbursts and the site survey revealed that the main damage area of the roadway was 68 m outside the front of the working face because of the impact. This corresponds to the relatively high value of the electromagnetic radiation in the early stage, which is a good precursor phenomenon, and the other two damaged areas correspond to relatively low points of electromagnetic radiation. After the impact, the stress dropped sharply and the electromagnetic radiation intensity presented a relatively stable wave-like distribution in the direction of the working surface. There was no obvious high-value concentration area, indicating that after the impact the coal in the broken zone is completely unstable and the plastic zone and elastic zone are again affected by the supporting stress. Accurate identification of the high-stress zone of the roadway is of great significance for strengthening the support to prevent roadway deformation, as well as roadway pressure shock.

Limitations and Perspectives
The limitations of this paper are mainly reflected by these points: (1) For the plastic zones, qualitative analysis is only carried out on the distribution area of shear failure and tensile failure. In the future, the volume of the plastic zones of tensile stress and shear stress failure can be further obtained for quantitative analysis; (2) In this paper, the internal friction angle is used because the FLAC3D software simulation needs to input the internal friction angle, but the Mohr-Coulomb modular configuration is to calculate the tangent of the internal friction angle, that is, the friction coefficient. The internal friction angle simulated in this paper increases linearly, while the corresponding friction coefficient changes nonlinearly; (3) The orthogonal test has the characteristics of high efficiency, convenience, and few tests, but the method obtained by orthogonality is only a series of combinations of representative test level parameters, and the results will not exceed the range of the selected level.
The suggestion for future work is to continue to improve the model and reduce the limitations, so that the model can more accurately describe the deformation and deterioration characteristics of the surrounding rock in the deep roadway [69].
(1) In a future model, we aim to simulate a more realistic deep mine environment, improve the geological conditions models (stratigraphic dip, fault occurrence, etc.), simulate the optimization of roadway support (bolt length, row spacing, and pre-tightening force, etc.), and determine the best way of support control to keep the roadway in a stable state; (2) The evaluation index of the stability of the surrounding rock of the roadway is not only limited to the displacement and deformation of the roadway, but also should integrate many factors, such as the extensive damage-depth of the surrounding rock caused by the deterioration of the roadway; the volume of the plastic zones; and the accumulation, dissipation, and release of elastic energy, etc.; (3) It is better to check the accuracy of the results of the orthogonal test through a comprehensive test. In this way, the modelling results can make up for the limitations of the small amount of data in the orthogonal test, obtain a larger quantity and wider range of information, and more accurately describe the influence of the surrounding rock under different parameters.

Conclusions
In this study, the FLAC3D software and orthogonal experiment method were used to analyze the influence of multiple factors, such as cohesion c, tensile strength σ t , internal friction angle ϕ, vertical ground stress p, and horizontal-to-vertical stress ratio λ, on the deformation of a deeply buried roadway. Simultaneously, the characteristics of the temporal and spatial evolution of roadway deformation during the stopping period are discussed through the characteristics of the displacement, stress, and plastic zone. Research results show that: (1) The physical and mechanical parameters of a rock mass play a major role in the roadway deformation. The order of influence degree under different parameters is: ϕ > p > σ t > λ > c. The changes in the internal friction angle ϕ and vertical ground stress p have the most significant effects. The internal friction angle is negatively correlated with roadway deformation, but the vertical ground stress is positively correlated; (2) The tensile strengt σ t and horizontal-to-vertical stress ratio λ are most sensitive to the stress on the roof and the two sides, respectively, and both have a negative correlation. When cohesion c and vertical ground stress p increase, the stress on the two sides first increases and then decreases. The vertical stress on the floor changes very little; (3) The tensile failure zones are distributed on the two sides and floor on the roadway, extending inward in the direction of excavation. Besides this, the scope of shear failure is most widely distributed and concentrated on the shoulders of the two sides and roof, and shear failure gradually extends from the roof in the left and right directions; (4) With the increase of the surrounding rock grade, the deformation and the stress of the two sides and roof rise, whereas the stress acting on the floor is not obvious; the stability continues to weaken, the deterioration and deformation are more serious, the distribution of the plastic zone is wider, and more shear failure extends in the "X" shape; (5) On-site electromagnetic radiation can effectively characterize the stress distribution in the roadway; when the electromagnetic radiation detected the overall high-intensity area, it reflects the distribution of the high-stress concentration area and the strength reduction area, which strongly supports the characteristics of coal and rock failure.