Study on the Large Deformation Characteristics and Disaster Mechanism of a Thin-Layer Soft-Rock Tunnel

State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China School of Mechanics and Civil Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China Deformation Survey Center, Yunnan Earthquake Agency, Kunming 650041, China


Introduction
e Muzhailing Tunnel in China is a critical control project along the Weiwu Expressway from Lanzhou to Haikou (G75), a national expressway. Because of its special geological conditions including deep burial, high ground stress, and a thin-layered slate structure, the tunnel deformation is extremely serious and has restricted the opening of the Weiwu Expressway. During the construction process, tunnel deformation and collapse occur frequently because of disturbance from excavation and the complex geological structure, seriously affecting the safety of the tunnel project [1][2][3][4]. erefore, it is urgent to investigate the deformation mechanism of the surrounding rock of the thin-layered slate tunnel to provide an important theoretical basis for support and control measures after future deformation of the surrounding rock.
Many researchers have investigated the deformation of deep-buried tunnels and have obtained useful findings [5][6][7][8]. In 1999, Brox and Hagedorn [9] studied the deformation characteristics of a tunnel crossing a fault along the Anatolian Highway in Turkey.
e study revealed that structural stresses, weak surrounding rock, and stress concentrations after excavation were the main causes of the deformation of the tunnel. In 2000, Hoek and Marinos [10] identified through a study of rock and rock mass strength that rock mass strength and ground stress were the key factors affecting the stability of the surrounding rock of the tunnel. In 2005, Liu et al. [11] studied the damage characteristics and causes of the large deformation of the Muzhailing Tunnel and suggested that the deformation and failure of the tunnel had serious, long-term and asymmetric characteristics. e deformation was found to be mainly caused by the plastic flow and expansion of the excavated unloaded surrounding rock under the combined action of lithology, groundwater, and ground stress field.
Apart from experimental investigations, many studies have numerically analyzed the deformation of deeply buried tunnels. In 2006, Meguid and Rowe [12] studied the deformation and stress distribution of the surrounding rock of a tunnel under different stress conditions by means of a three-dimensional numerical simulation method. Guo et al. [13] analyzed the anisotropy of the strength and deformation of layered surrounding rock through an experiment of basic mechanics and studied the deformation and failure mechanism of the highly stressed and steeply dipping layer of the Yangjiaping Tunnel by numerical analysis. Yu [14] in 2017 studied the deformation mechanism and control measures of the Shiziyuan Tunnel by means of theoretical analysis and on-site measurement and found that the main factors of the deformation were ground stresses, mechanical characteristics of surrounding rocks, type of tunnel-lining structure, and construction disturbance (blasting). In 2018, Cao et al. [15] analyzed the deformation and failure characteristics of the initial support of a typical weak tunnel, Dongshan Tunnel, and obtained the factors influencing the stability of surrounding rocks in the construction process. Sun et al. [16] in 2019 studied the application of a constant resistance large deformation anchor cable support system by performing a field test and numerical simulation. e results showed that the anchor cable effectively controlled the large deformation of the surrounding rock of the tunnel. e above studies have revealed many beneficial findings on deep-buried soft-rock tunnels. However, there are few studies on the deformation and failure mechanisms of the surrounding rock in the complex engineering geological environment of the Muzhailing deep-buried thin-layered slate tunnel. e factors affecting large deformations of the surrounding rock and the degree of influence of each factor on the deformation of the surrounding rock of the tunnel are not adequately understood. erefore, it is necessary to postulate relevant theories and perform numerical analysis.
is study considers the No. 2 inclined shaft of the Muzhailing Tunnel, investigating the deformation characteristics of the thin-layered slate tunnel by statistically analyzing the deformation of the tunnel under the original supporting conditions. Combined with laboratory experiments, the main factors affecting tunnel deformation are then analyzed. At the same time, the effects of different joint inclinations of the surrounding rock, buried depths, water absorption, and softening on the large deformation of the surrounding rock are analyzed by the 3DEC-Tragon numerical simulation method.

Project Overview
e Muzhailing Tunnel is located in the low and middle mountainous area of the West Qinling Mountains (Figure 1). It is located at the junction of several plates, crosses the Muzhailing Watershed between the Zhanghe River and Taohe River and spans Zhangxian County and Minxian County of Gansu Province. e total length of the tunnel is 15.15 km, the maximum buried depth of the tunnel body is 629.1 m, and the minimum buried depth is 8 m. e tunnel is composed of three inclined shafts, of which the No. 2 inclined shaft is located in the Muzhailing Tunnel (middle section, red framed position in Figure 1). e starting and ending pile numbers are (z) K216 + 380 − (z) K220 + 300, and the length of the Muzhailing Tunnel is 3920 m. e opening of the No. 2 inclined shaft is located on the north side slope of Luzagou. e pile number in the main tunnel at its intersection with the inclined shaft is K218 + 400. e intersection angle with the main tunnel plane is 90°. For the inclined shaft, the starting pile number is K0 + 000 and the end pile number is K1 + 813 m, with the length of the shaft being 1813 m.

Formation Lithology Characteristics.
According to the field investigation and exploration, the bank slope is covered with quaternary slope gravel and the underlying bedrock is the lower Permian slate. e main components of the surrounding rock are carbonaceous slate with local intercalations of sandstone, limestone, and shale. Quartzite is sporadically distributed in the rock formations, and the slate in this area is mainly grey-black. Joints are relatively developed, the rock mass is relatively broken, and the strength markedly decreases under the action of water-rock coupling. e surrounding geological structure largely affects the surrounding rock, and the whole area is essentially grade V rock.

Hydrogeological Characteristics.
e distribution of groundwater was relatively uneven and was strongly related to the geological structure, lithology, and topography of the region. As the area was strongly affected by geological actions in the past, the joints and fissures in the area were well developed and the groundwater mainly was present in the form of fissure water. According to the geological survey, the valley area was relatively rich in water content and the water collection capacity was strong, so the water pouring was more serious when the tunnel passes through the area. Besides fissure water, karst water and Quaternary loose rock void water were encountered in this area.

Ground Stress Characteristics.
According to previous studies [17][18][19], hydraulic fracturing was used to measure the ground stress at two measurement points: AK214 + 738 m and AK219 + 422.8 m, at test depths of 492.00 and 499.00 m, respectively. According to the test results, the stress condition was such that the horizontal maximum principal stress S H was greater than the horizontal minimum principal stress S h , which was greater than the straight principal stress S v . e ground stress was mainly tectonic stress. e ratio of stress size to R c in the Muzhailing Tunnel area is shown in Table 1.

Advances in Civil Engineering
According to the Standard for Engineering Classification of Rock Mass (GB/50218-94), R c /S H < 4 is extremely high stress, 4 < R c /S H < 7 is high stress, and 7 < R c /S H is low stress. Table 1 shows that the value of R C /S H in the Muzhailing Tunnel is < 4, which is categorized as extremely high ground stress. Moreover, the stratum rock mass is relatively broken and the possibility of large deformation in the surrounding rock of the tunnel is relatively high.

Properties of Rock Mass.
e surrounding rock inclination angles and surrounding rock integrity of the No. 2 inclined shaft of the Muzhailing Tunnel in the large deformation area vary widely.
e dominant joint angle is between 30°and 85°, and the thickness of the slate layer is 5-30 cm. To understand the mineral composition and mechanical properties of the slate, an X-ray diffraction experiment, a point-load experiment, and a uniaxial compression experiment were performed.

Mineral Composition Analysis.
Six positions (pile number K1430/K1460/K1490/K1520/K1550/K1570) in the No. 2 inclined shafts of the Muzhailing Tunnel were selected as sampling points, and the selected carbonaceous slate samples were quantitatively analyzed by using a TTR III multifunctional X-ray diffractometer. e mineral composition and relative content of rock and clay were obtained as shown in Tables 2 and 3. A small difference was observed in the mineral composition among the six groups, particularly in the quartz and clay minerals, which mainly included illite, chlorite, illite mixed layer, and kaolinite. e montmorillonite content was low, and the samples exhibited low water absorption and expansion.

Effect of Water on the Strength of Slate.
Because of the fissured joints and severe weathering of the surrounding rock of the No. 2 inclined shaft of the Muzhailing Tunnel, the lithology is mainly carbonaceous slate with scattered sandstones. e rock samples were significantly affected by the disturbance during sample acquisition, and the processing of standard samples was found to be challenging. erefore, a field point load intensity meter was selected for a rapid experimental analysis. During the selection of specimens, the integrity of the selected rock should be ensured (avoiding rocks with cracks). A total of 82 specimens of carbonaceous slate from K1740 to K1800 in the No. 2 inclined shaft were selected in the field, of which 41 specimens were tested for uniaxial compressive strength in the natural state, and the other 41 specimens were tested for uniaxial compressive strength after 24 hours of water absorption. e experimental results for the 82 specimens are shown in Figure 2. e uniaxial compressive strength was estimated by the point-load method, and the average uniaxial compressive strength of the carbonaceous slate in the natural state was found to be 29.6 MPa. e average uniaxial compressive strength of the carbonaceous slate after water absorption for 24 hours was 12.7 MPa. e uniaxial compressive strength after water absorption was reduced by 57% when compared with the natural state.

Uniaxial Compressive Strength
Test. Slate with visible joints was processed into standard V50 mm × 100 mm specimens. Specimens with 30°, 45°, and 60°bedding inclinations were tested for their uniaxial compressive strength by using a 2000 kN capacity mono-triaxial test system. is test revealed the influence of different joint surface   Table 4. According to the above results, the uniaxial compressive strength and elastic modulus of the rock specimens were relatively low when the joint surface of the three groups with different joint inclination angles was 45°(that is, the angle between the joint surface and the horizontal plane was 45°). For inclination angles between 30°and 60°, the surrounding rock essentially failed in shear.

Failure Characteristics of the Tunnel.
To further study the large deformation characteristics of the surrounding rock of the Muzhailing Tunnel, we used a three-dimensional laser scanner and a total station together ( Figure 3) to monitor eight sections of the No. 2 inclined shaft of the Muzhailing Tunnel, including XK0 + 060-XK0 + 250, XK0 +     e results of the monitoring are shown in Table 2. Figure 4 shows the initial support damage as photographed on-site. According to Table 5 and Figure 4, the initial deformation and cracking of the tunnel support were severe, the steel frame was distorted, the deformation rate was high, and the amount of deformation was large.

Deformation Characteristics of Surrounding Rock.
On the basis of the deformation of the surrounding rock of the No. 2 inclined shaft, the deformation of the surrounding rock of the tunnel was classified according to the standard of tunnel deformation. e length of the tunnel was K1713 m, which is divided into four zones: basic stable zone (maximum deformation of less than 100 mm), slight deformation zone (maximum deformation of 100-250 mm), larger deformation zone (maximum deformation of 250-1000 mm), and severe deformation zone (deformation of more than 1000 mm).
rough the above analysis, 41% of the total length of the tunnel was found to be in the large deformation zone and the buried depth of the severe deformation zone was more than 500 m. e length of the basic stable, slight deformation, larger deformation, and severe deformation zones were 407, 250, 717, and 338 m, respectively. e deformation of the surrounding rock in the essentially stable area was analyzed in this study. After the initial support of the tunnel was completed, the deformation of the surrounding rock increased sharply and the arch sank by 200-600 mm. Generally, the tunnel reached basic stability within 14 to 18 days, while the deformation continued by 1-2 mm per day. From the displacement distribution along the entire section, the deformation of the tunnel in a given section was asymmetric: the deformation on the left side of the tunnel was significantly larger than that on the right side. Figure 5 shows the typical deformation curve of the surrounding rock.
From the field investigation, the laboratory experiments, and the analysis of the field-monitoring data, we found that the surrounding rock joints of the No. 2 inclined shaft of the Muzhailing Tunnel were developed, the slate layer spacing was 5-30 cm, and the structure largely influenced the strength of the rock mass. e mineral content of the slate clay was approximately 50%, indicating that the rock strength had a certain influence on the damage of the surrounding rock and that water largely influenced the rock strength.
e stress distribution of the surrounding rock varied with an increase in buried depth. e area around the

Advances in Civil Engineering
Muzhailing Tunnel experienced extremely high ground stress, and the deformation of the surrounding rock was significant. erefore, three main factors influencing tunnel deformation were determined from the experimental study: joint inclination angle variation of the surrounding rock, attenuation of the strength of the surrounding rock owing to water absorption, and the buried depth of the tunnel. To further reveal the deformation and failure mechanism of the surrounding rock of the Muzhailing Tunnel, a numerical simulation was performed by considering these three aspects: joint inclination angle, burial depth, and softening caused by water absorption.

Numerical Analysis of the Deformation
Mechanism of the Thin-Layer Slate Tunnel 4.1. Modelling Approach. e software 3 Dimension Distinct Element Code (3DEC) is a three-dimensional discrete element program that represents the discontinuous discrete characteristics of a material. It can explain the interaction between rock blocks and structural surfaces in rock masses,   [21] used the 3DEC joint constitutive model to study the deformation characteristics of layered rock mass in a roadway under the action of the supporting body. Hou [22] used the 3DEC numerical analysis software to analyze the deformation and failure characteristics of a roadway tunnel supported by an anchor mesh and obtained the stress distribution characteristics, displacement characteristics, and plastic zone distribution characteristics of the surrounding rock.
To determine the characteristics of prominent dominant joints and coexistence of multiple joint forms in the No. 2 inclined shaft of the Muzhailing Tunnel, a simulation method using equivalent ubiquitous joints of dominant joint groups was proposed in this study, as shown in Figure 6. Considering the computational cost, we did not individually simulate all joints in the numerical model and we generalized the distribution of surrounding rock joints as shown in Figure 6(b). e joints were grouped into dominant joint groups (Figure 6(c)) and other joints. Locally, other joints (except for superior joints) had anisotropic characteristics, while macroscopically they were regarded as isotropic.
erefore, the 3DEC-Trigon method was used to simulate the rock mass as continuous media equivalence. e simulation considered the representative elementary volume (REV) characteristics of the rock mass and the transformation of rock-to-rock mass parameters in the process of equivalence [23,24]. Compared with other methods, this method can better adapt to the actual characteristics of the surrounding rock of the Muzhailing Tunnel.

Determination of Numerical Simulation Parameters.
According to the principle of joint parameter selection [25], joints were modelled as uniformly distributed. e mechanical properties of the rock were obtained from equations (1) to (4), and the edge length of the irregular blocks was selected as Δz min � 0.8.
where K is the bulk modulus, G is the shear modulus, K n is the normal stiffness of the contact model parameters, and K s is the shear stiffness.
According to the geological strength index (GSI) rock classification standard [26][27][28][29][30], for the mechanical properties as obtained from equations (1) to (4), the numerical inversion analysis was used during the simulation process. Accordingly, the strengths of carbonaceous slate in the natural state and under the saturated condition were obtained as shown in Tables 6 and 7. On the basis of multiple numerical analyses, the REV with a stable structural plane combination was found to be 8 × 8 × 16 m. e results of the numerical test on the uniaxial compressive strength of the rock mass are shown in Figure 7.

Advances in Civil Engineering
After several inversion analyses combined with field measured deformation data of the surrounding rock, the dominant joint parameters selected were jkn � 2e9, jks � 1e9, jcoh � 1e5, jfr � 25, and jten � 5e4. Figure 8. e model includes the tunnel section, uniform joints, and dominant joints. e ground stress was obtained from the field investigation as discussed in Section 2.3. e size of the model was 60 m × 30 m × 60 m (length × width × height). e reverse arch at the bottom of the tunnel was filled with concrete. Fixed constraints were imposed at the bottom of the model, and other boundary conditions were subjected to stress constraints based on the ground stress.

Effect of Joint Inclination on the Deformation of the Tunnel.
Using the numerical model presented in Section 4.3, we estimated the ground stress when the buried depth was 500 m. e vertical stress was 12.5 MPa, the maximum horizontal stress was 21 MPa, and the minimum horizontal stress was 12.5 MPa. is simulation did not consider the tunnel support. By changing the dominant joint angle, we analyzed the deformation characteristics, stress distribution, and plastic zone expansion of the surrounding rock under different joint inclination angles. Figure 9, bending deformation and block slip deformation were caused by tunnel excavation, with bending deformation being the primary deformation encountered by the tunnel. e maximum bending deformation occurred at the left shoulder socket of the tunnel. e bending deformation was mainly due to an increase in the deflection of the plate-like surrounding rock under the action of the ground stresses, and the slip deformation was mainly the result of the shear slip of the structural surface. e deformation of the tunnel was influenced by the presence of joints and was consequently asymmetric. Among the joint inclination angles of 30°, 45°, and 60°, the most evident asymmetry was observed for a joint inclination angle of 45°. As shown in Figure 9, the displacement distribution is in the form of asymmetric butterfly displacement. e maximum value of bending deformation corresponding to the 30°joint inclination was 0.7 m, and the displacement and deformation essentially stopped when the calculation reached 8000 steps. A separation layer developed at a depth of approximately 3.5 m in the surrounding rock, and sliding failure occurred in the surrounding rock of the tunnel. e maximum value of bending deformation corresponding to the 45°joint inclination was 1.2 m, and the displacement deformation essentially stopped after 10000 steps. A separation layer developed at a depth of approximately 7 m in the surrounding rock, the mass spalling caused by tunnel slip failure gradually increased, and the deformation of the surrounding rock of the tunnel presented strong

Saturated
Step  asymmetry.
e maximum bending deformation corresponding to the 60°joint inclination was 0.8 m, and the displacement deformation essentially stopped when the calculation reached 9000 steps. A separation layer developed at a depth of approximately 4 m in the surrounding rock.

Trend in the Maximum Principal Stress Expansion.
e original stress field is rearranged after the excavation of a tunnel. e distribution of the maximum principal stress is an important basis for judging the stability of the geotechnical engineering, and it is an important index for understanding the range of influence of the excavation. As shown in Figure 10, the maximum principal stress transferred to a greater depth into the surrounding rock with the deformation of the surrounding rock, and the range of influence increased. After the tunnel deformation became stable, the maximum principal stress was concentrated in the vault and the bottom plate of the tunnel, and the maximum principal stress of the vault reached 32 MPa (showing a stress

30° inclination 45° inclination 60° inclination
Step 2000 Step 4000 Step 6000 Step 8000 Step 10000  e maximum major principal stress of the vault occurred at 11 m from the tunnel wall. At increasing depth of the surrounding rock, the major principal stress decreased. At a distance of 16 m from the tunnel wall, the major principal stress transitioned to the original rock stress state, indicating that the range of influence of tunnel excavation was up to 16 m.

Impact of Different Buried Depth Conditions on the Surrounding Rock of the Tunnel.
e above analysis shows that the deformation of the surrounding rock is most evident when the joint inclination angle is 45°under the same buried depth. erefore, this section considers the dominant joint with a 45°inclination angle as the structural condition to study the variation in the stress, strain, and plastic zone of the surrounding rock caused by tunnel excavation under different buried depths.
Combined with the site-tunnel buried depth conditions, three tunnel excavation depths H were selected: 250, 330, and 500 m. e displacement distribution of the surrounding rock after tunnel excavation under the three buried depth conditions is shown in Figure 11. e simulation results show that (1) the deformation of the surrounding rock of the tunnel under the influence of the dominant joint is asymmetric under different buried depth conditions, showing that the left side deformation is larger than the right side, the maximum deformation position is located at the left shoulder socket of the tunnel, and the maximum deformation amount and deformation range of the tunnel gradually expand with an increase in buried depth and (2) for H � 250 m, the tunnel does not have obvious damage, the tunnel section shape is essentially good, and the maximum displacement is 350 mm. When H � 330 m, the left shoulder socket of the tunnel appears to have a detached layer, the surrounding rock of the tunnel has bending damage, and the maximum displacement is 520 mm. When H � 500 m, thin-plate bending and damage occurred at the left side and the left shoulder fossa of the tunnel, and the maximum deformation is 1100 mm. e distribution characteristics of the plastic zone of the surrounding rock excavated by tunnels with different depths are shown in Figure 12. e simulation results show the following. (1) e distribution of the plastic zone is crossshaped, with vertical dominant joints and parallel dominant joints in the direction of the cross axis. As the burial depth increases, the plastic zone expands gradually and the plastic zone of the surrounding rock mass has shear failure. Primarily, the shear displacement direction is closely related to the direction of the joint surface. e bottom shear zone gradually transitions to the "W" shape under the action of in situ stress. Both ends of the "W" shape intersect the tunnel bottom angle and expand upward along the side, gradually closing up with the upper "M"-type plastic zone. (4) When H � 500 m, the plastic zone inherits the trend when H � 350 m, the two plastic zones at the top and the shoulder form a large-scale penetration, and the bottom "W" plastic zone completely closes with the upper "M" model. Eventually, a larger "¤" plastic zone forms. e maximum principal stress distribution of the surrounding rock under different buried depths is shown in Figure 13. e simulation results demonstrate the following.

Effect of Softening of the Surrounding Rock on the
Deformation of the Tunnel. During the construction of the No. 2 inclined shaft of the Muzhailing Tunnel, water gushing from the fissure was often encountered, and the deformation of the tunnel was found to be more significant in areas where water gushing was large. According to Section 2.4.2, the uniaxial strength of carbonaceous slate after water absorption was 0.4 times the uniaxial strength in the natural state, and the rock strength significantly decreased. erefore, it is necessary to analyze the influence of water absorption and softening on the deformation of the tunnel. Combined with the difficulty of on-site engineering operations, the influence of water absorption and softening of surrounding rock on the surrounding rock of the tunnel was analyzed under working conditions of a tunnel burial depth of 500 m and a dominant joint angle of 45°.
By reducing the strength of surrounding rock to 40% of the strength before water absorption, we compared the deformation field, plastic zone, and stress field distribution characteristics of the surrounding rock before and after water absorption ( Figure 14). (1) e characteristics of the deformation distribution of the surrounding rock before and after water absorption change greatly. e deformation of the tunnel before water absorption is mainly distributed in the left shoulder socket with a maximum deformation of 1200 mm. e left shoulder socket tends to be unstable, while the deformation of other parts is small and there is no instability. After water absorption and softening, the deformation of the surrounding rock at the left shoulder socket is the largest, with a maximum of 1900 mm, while the tunnel floor and the right side also have large deformation, and the deformation range of the surrounding rock extends to a greater depth. e hole has a symmetrical "C" shape, which indicates that the decrease of the surrounding rock strength after the water absorption softening helps control the deformation of the tunnel, while the dominant joints have an effect on the deformation direction. (2) e plastic zone of the surrounding rock is also evidently expanded after water absorption softening, which is affected by the decrease in the strength of the surrounding rock. Furthermore, the upper plastic zone completely joins and the dominant joints still affect the distribution characteristics of the plastic zone. e radius of the plastic zone after water absorption softening can reach 14.5 m, which is approximately 4.2 m more than before water absorption, indicating that the distribution range is mainly controlled by the strength of the surrounding rock. (3) e stress distribution of the surrounding rock also changes greatly before and after water absorption softening. e maximum principal stresses before and after water absorption are concentrated in the deep surrounding rock at the top and bottom of the tunnel, the maximum          Advances in Civil Engineering 13 principal stress concentration area after water absorption moves to 0.5 times the diameter of the hole to a greater depth, and the maximum principal stress value decreases by approximately 4 MPa after water absorption softening.

Conclusions
In view of the large deformation of the surrounding rock during the construction of the No. 2 inclined shaft of the Muzhailing Tunnel, the deformation characteristics and influence of the surrounding rock of the tunnel under the original supporting form were analyzed. e following conclusions are drawn: (1) e initial support of the No. 2 inclined shaft of the Muzhailing Tunnel was seriously damaged, and the deformation of the surrounding rock of the tunnel was asymmetric. e deformation on the left side was significantly larger than that on the right side, and the area with a deformation greater than 0.25 m accounted for 41% of the total length of the tunnel. In addition, when the buried depth was more than 500 m, serious deformation of the surrounding rock exceeding 1 m was observed. (2) According to the field investigation, the laboratory experiment, and the field-monitoring data analysis, the surrounding rock of the No. 2 inclined shaft of the Muzhailing Tunnel was relatively broken, the lithology was mainly carbonaceous slate, the joints and fractures were well developed, and the area was subjected to extremely high stress. e uniaxial compressive strength after 24 hours of water absorption was markedly lower than that of carbonaceous slate in the natural condition. e uniaxial compressive strength and elastic modulus of the rock with a 45°joint inclination were relatively low when compared with inclinations of 30°and 60°. erefore, the main factors affecting the deformation of the thin-layer rock of the Muzhailing Tunnel included the joint inclination angle, buried depth of the tunnel, and softening of the surrounding rock due to water absorption.
(3) rough the numerical simulation results, it can be concluded that under the same buried depth, the deformation of surrounding rock was maximum for a 45°joint inclination. At the 45°joint inclination, the deformation and failure of the tunnel was intensified with an increase in burial depth. When the burial depth was 500 m, thin-plate bending and damage occurred at the left side and the left shoulder fossa of the tunnel; after the surrounding rock was softened by water absorption, the floor of the tunnel, the left shoulder socket, and the right side of the tunnel deformed greatly and the deformation range of the surrounding rocks extended to a greater depth.

Data Availability
All the data included in this study are available upon request by contact with the corresponding author.  Figure 14: Deformation characteristics of the surrounding rock of the tunnel under water absorption and softening: (a) change in the surrounding rock displacement before and after water absorption, (b) changes in the plastic zone of the tunnel before and after water absorption, and (c) changes in the maximum principal stress field of the tunnel before and after water absorption.