1. Introduction
The tunnel construction was prone to rockburst [
1,
2,
3], large deformation [
4,
5,
6], and other diseases under the condition of deep-buried depth and high tectonic stress. Large deformation of tunnels in the case of soft rock and high in situ stress was more common than rockburst. The large deformation of the tunnel would aggravate the instability of the surrounding rock and bring about the distortion of the steel frame, the cracking of the initial support, the invasion limit, the loss and cracking of the secondary lining, the inverted arch uplift, and other related diseases [
7,
8], which would seriously affect the construction safety and progress and bring great risks to the construction of the tunnel. To control large deformation disasters, strengthening pre-supporting, improving support rigidity, adopting double-layer initial support, grouting reinforcement, and closing and forming a ring as soon as possible were the commonly used methods [
9,
10,
11]. These methods had achieved good results in many cases. Based on the deformation characteristics of the surrounding rock, Zhang et al. [
12] studied the reasonable tunnel shape of the high in situ stress and steeply inclined surrounding rock and believed that the optimization of the cavity type could fundamentally reduce the stress and deformation of the tunnel structure. Xu et al. [
13] applied theories and tests to propose a steel grid concrete core tube support system, which could improve the support quality and prevent large deformation in tunnels. Qiu et al. [
14] optimized the grille steel frame through element tests and refined the model calculation. This result did not significantly increase the number of steel bars but increased the bearing capacity of the initial support, which could effectively resist the deformation of the tunnel. For tunnels with huge deformations, Wang et al. [
15] applied high-stiffness steel tube concrete as the initial support to resist the deformation of the surrounding rock, which achieved good results. In addition to the large-rigidity support method for resisting the deformation, Qiu et al. [
16] invented the resistance limiter for the initial support of the tunnel by using different stiffnesses between the materials. In addition to the large-rigidity support method for resisting the deformation, Qiu et al. [
16] invented the support resistant limiting damper (SRLD) for the initial support by using the different stiffnesses between the materials. SRLD controlled the deformation by first letting down and then resisting, which had been used in the Menghua Railway. Various types of anchors had also been used in tunnels, such as NRP anchors [
17], high constant resistance anchor cables [
18], prestressed anchors (anchor cables) [
17], etc. These absorbed the stress release caused by tunnel excavation in the surrounding rock and had a good effect in high in situ stress tunnels. In addition to the need for reasonable support parameters and construction methods, it was of great significance to grasp the deformation characteristics of the tunnel and select a reasonable amount of reserved deformation [
19,
20,
21,
22].
The current “Specifications for Design of Higbway Tunnels
Section 1 Civil Engineering” (JTG3370.1-2018) [
23] and “Code for Design of Railway Tunnels” (TB 10003-2016) [
24] both had instructions for the reserved deformation for the classification of surrounding rock. They were the main basis for the design of the reserved deformation of Chinese tunnels. Meanwhile, based on the deformation field survey and on-site monitoring and measurement data of the Lianchengshan tunnel, Chen et al. [
19] statistically analyzed the deformation of the tunnel under different chlorite schist rock mass conditions, studied the defamation law, and reasonably reserved the amount of deformation. Deformation control benchmarks for deep-buried and large-span chlorite schist tunnel construction were established. For different rock mass states, four different reserved deformations of 70~95 cm, 50~70 cm, 30~50 cm, and 15~30 cm were given. Liu et al. [
25] monitored a total of 471 sections in 14 tunnels of the Yuchu Expressway. He used the interval statistical method to calculate the reserved deformation of tunnels at different levels of surrounding rock and at different buried depths. It is recommended that when the buried depth of Grade IV surrounding rock is less than 50 m, the reserved deformation is set to 17 cm. It is recommended that when the buried depth is 50~300 m, the reserved deformation is set to 50 cm. It is recommended that when the buried depth of Grade V surrounding rock is less than 50 m, the reserved deformation is set to 18 cm. It is recommended that when the buried depth is 50~300 m, the reserved deformation is set to 32 cm. Zhao et al. [
26] proposed a method of optimizing the reserved deformation based on the whole-section analysis and on the 3D laser scanning technology. The step-like method of setting aside the amount of deformation was adopted (the dome was 7 cm, and the rest area was 5 cm). Li et al. [
27] used the theory of elastoplastic mechanics to derive the analytical solutions of the surrounding stress and displacement at different stages and the radius of the crack and plastic zone. In response to the problems in a tunnel in Xuzhou, 20 cm reserved deformation layers were proposed. Luo et al. [
28] analyzed the stress and deformation characteristics of the tunnel, using numerical calculation software, and obtained the reserved deformation of the loess tunnel with different moisture contents. Cui et al. [
29] determined the asymmetric reserved deformation value of the loess tunnel through the method of a model test and believed that the reserved deformation of the vault should be two to four times that of the sidewall. Wang et al. [
30] gave the respective reserved deformations under the three conditions of hard plastic, soft plastic, and flow plastic for large-section double-track loess tunnels. They were hard plastic (vault 55–70 mm, sidewall 15–20 mm), soft plastic (vault 166–180 mm, sidewall 40–50 mm), and flow plastic (vault 290–300 mm, sidewall 125–140 mm).
Many scholars have carried out in-depth research on the reserved deformation in tunnel construction and put forward the corresponding reasonable reserved deformation according to different soft rock geological conditions. These research results are mainly related to the strength of the surrounding rock, buried depth (i.e., geostress level), and other key factors. For soil tunnels, especially loess tunnels, the distribution of reserved deformation is obviously uneven, and the important variable of water content should be additionally considered. However, at present, most of the research on reserved deformation focuses on tunnels with buried depths less than 500 m. With the development of tunnel engineering technology to a deeper level, especially for tunnels with a depth of 1000 m, the in situstress level is much higher than that of conventional tunnels, which makes the prediction and control of tunnel deformation more complicated and arduous. Under such extreme conditions, how to determine the reserved deformation scientifically and reasonably to ensure the stability and safety of the tunnel structure has become a technical problem to be solved urgently. In addition, for tunnels with a large cross-section, high ground stress, soft rock, and other special working conditions, the current design specifications often lack specific guidance on reserved deformation, which further increases the construction difficulty and risk. Therefore, future research should focus on the reserved deformation law of tunnels under these complex conditions and explore more accurate and reliable prediction and control methods of reserved deformation through theoretical analysis, numerical simulation, field tests, and other means so as to provide strong support for the safe and efficient construction of tunnel projects. In this paper, the deformation characteristics of a diversion tunnel of Jinping II Hydropower Station are studied. Based on field statistics, theoretical calculation, and numerical simulation, and considering the size, in situ stress, and grade of the surrounding rock, the reasonable reserved deformation of the deep soft rock tunnel is studied.
2. Project Overview and Geological Conditions
The Jinping II Hydropower Station is located on the Jinping River Bay, the mainstream of the Yalong River in Sichuan Province. It was the largest hydropower station in the cascade development of the Yalong River and one of the backbone power stations for power transmission from west to east. The terrain along the tunnel is undulating. The height of the mountain at the western end of the tunnel entrance is about 500 m, and then the buried depth gradually increases. The project location and cave entrance photos are shown in
Figure 1. The diversion system of Jinping II Hydropower Station is meticulously built, including Jinping auxiliary tunnel lines A and B, four diversion tunnels, and drainage tunnels, for a total of seven major tunnels placed in parallel, as illustrated in
Figure 2a. The Jinping auxiliary tunnel was built using efficient drilling and blasting techniques, and its standard section was planned as a distinctive gate-shaped section with an excavation area of 45 m
2. The horseshoe-shaped section facing the emergency stop area has been adjusted for safety, and the excavation area has been increased to 55 m
2. The drilling and blasting process and advanced tunnel boring machine (TBM) technology are used in a unique way to construct a diversion tunnel. For the section excavated by drilling and blasting, advanced bench construction technology was used, resulting in a horseshoe-shaped section with an excavation diameter of 13 m and an upper excavation area of 110 m
2, significantly improving the construction efficiency and quality. To address the project’s water inrush issue and ensure the smooth construction of the diversion tunnel, the project team expertly built a critical drainage tunnel between auxiliary tunnel B and diversion tunnel 4. This drainage tunnel fully utilizes the full-face TBM excavation technology. The excavation diameter is precisely controlled at 7.2 m, and the excavation area is efficiently maintained at 41 m
2, ensuring effective water flow management and material transportation during the construction period while demonstrating the superior engineering design and construction level.
The Jinping II Hydropower Station was affected by the collision of the Eurasian plate and the Indian Ocean plate and mainly composed of middle and upper Triassic marble, limestone, crystalline limestone and sandstone, and slate. The diversion tunnel also includes lower Triassic green sandstone and chlorite schist. From east to west, the tunnels passed through various strata including marble of the Yantang (T
2y) group, marble of the Baishan group (T
2b), upper Triassic (T
3) marble, marble of the Zagunao group, lower Triassic (T
1) chlorite schist, and metamorphic medium-to-fine sandstone, as shown in
Figure 2b. According to the tectonic history and orientation, the faults in the project area could be divided into four structural groups in the NNE direction, NNW direction, NE-NEE direction, and NW-NWW direction (all have a steep inclined angle), respectively. The surrounding rock is dominated by class II and III (hydraulic classification), and class IV and V account for about 23% to 30%, mainly concentrated in the T1 strata. The buried depth of tunnels was generally 1500–2000 m, with the maximum buried depth being 2525 m. The maximum in situ stress in the engineering area is 80 MPa, the uniaxial compressive strength of chlorite schist is 19.47 to 38.80 MPa, and the strength–stress ratio of the surrounding rock is 0.24 to 0.49, which has strong and extremely strong large deformation conditions. The prediction and prevention of large deformation is the key to ensuring construction safety.
5. Reserved Deformation
Because the initial support is closely attached to the surrounding rock, the local cracking and damage of the initial support will not cause the collapse of the entire tunnel. With the increase in the support displacement, the surrounding rock exhibits plastic deformation behavior, which is accompanied by the increase in and development of damage. When the damage develops to a certain extent, the initial support will not play the role of support. This is the excessive destruction of the initial support, and the stable state where the initial support and the surrounding rock deform together is the limit state of the initial support system.
5.1. Real Data Statistics
According to the statistics of the cross-sectional scan data of the deformed tunnel section (Tunnel 1 + 536~759 and Tunnel 2 + 613~643), when different reserved deformations were given, the corresponding guarantee rate relationship was as shown in
Figure 13 and
Table 3.
Figure 14a shows the statistics of the section scanning data of Tunnel 1. When the design reserved deformation was 40 cm, 50 cm, and 60 cm, the guarantee rates were 59.05%, 67.62%, and 78.01%, respectively. The guarantee rates were the probability that there would be no deformation greater than the reserved deformation amount. It could be seen that the intrusion limit in the range of 1 + 660~720 m was significantly larger (generally above 60 cm), and the intrusion limit in other parts of the cave was about 40 cm. Tunnel 1’s design may have a limited deformation range of 50–60 cm due to discrete on-site measurement data and a higher guarantee rate.
Figure 14b shows the statistics of the section scanning data of Tunnel 2. When the design reserved deformation was 40 cm, the guarantee rate reached 100%. Tunnel 1 and Tunnel 2 were parallel and separated by 60 m, whose stratum lithology and in situ stress were comparable. Based on statistical data, Tunnel 2’s reserved deformation may be compared to that of Tunnel 1, and a thorough examination may reveal that Tunnel 2’s design reserved deformation range was between 30 and 40 cm.
5.2. Analytical Method
Within a deeply underground tunnel, a portion of the rock mass may transition into a plastic state due to either the local shear stress exceeding the shear strength of the rock mass or the secondary stress state of the surrounding rock exceeding the compressive strength of the surrounding rock. The Jinping diversion tunnel’s substantial buried depth led to the assumption that the lateral pressure coefficient was 1. The ultimate displacement around the tunnel for three different
pa values (
pa is in situ stress)—0, 1, and 2 MPa—was examined using the elastoplastic method.
where
is the displacement around the tunnel(m),
μ is Poisson’s ratio,
E is the elastic Modulus (MPa),
σZ is the initial in situ stress (MPa),
r0 is the tunnel radius (m),
R0 is the plastic zone radius (m), and σ
R0 is the radial stress on the boundary of the plastic zone.
The radius of the plastic zone (
R0) was calculated using the following formula:
where
c is the cohesion (MPa),
φ is the friction (°), and
pa is the supporting reaction force (MPa).
The radial stress
on the boundary of the plastic zone is calculated using the following formula:
The parameters in
Table 2 were used to calculate the displacement around the tunnel when the supporting resistance
pa is 0 MPa, 1 MPa, and 2 MPa under different in situ stress conditions under Type III and IV surrounding rocks.
Figure 15 showed the relationship between the displacement of the unsupported tunnel and the initial stress.
Figure 16 and
Figure 17 are the curves when the supporting reaction force was 1 MPa and 2 MPa.
The limit displacement around the tunnel under the condition of Type IV surrounding rock was obviously greater than that of Type III surrounding rock. The larger the hole diameter, the greater the displacement. With the increase in in situ stress, the trend of displacement increases is more obvious.
According to the disclosed T1 stratum (Type IV surrounding rock), it was known that the buried depth was 1550~1850 m, and the in situ stress was about 45 MPa. When the tunnel diameter was 13.4 m, the displacement around the tunnel was 70 cm without support. When the support resistance pa = 1 MPa, it was 45 cm. When the support resistance pa = 2 MPa, it was 33 cm. Considering the hysteresis of the supporting effect, the reserved deformation should have been greater than the displacement when there was supporting resistance and less than the displacement around the tunnel when there was no supporting resistance. Considering the support resistance of 1 MPa, the reserved deformation range can be 45~70 cm. Considering the support resistance of 2 MPa, the reserved deformation range can be 33~70 cm.
5.3. Numeral Calculations
The maximum deformation, which depended on the strength of the surrounding rock, the amount of in situ stress, the stiffness of the support, and the duration of the support, was intimately linked to the reserved deformation. This part used the limit displacement of the tunnel without support as a guide to determine the reserved deformation, taking into consideration the uncertainty of the supporting function.
The mechanical calculation parameters of the surrounding rock of the tunnel are also shown in
Table 1. The deformation of the surrounding rock under 40 working conditions and under four conditions of 10 MPa, 20 MPa, 30 MPa, 40 MPa, and 50 MPa and cavern diameters of 12 m, 13 m, 14 m, 15 m in grade III and IV surrounding rock conditions were simulated. The calculation model and related boundary conditions and parameters were the same as in
Section 4.2. The calculation results are shown in
Figure 18.
It is evident that the surrounding rock grade has the greatest impact on the surrounding rock deformation. With the improvement in the surrounding rock grade, the integrity of the rock decreases and the development of cracks increases, which leads to a significant decrease in the strength and stiffness of the rock, which makes the surrounding rock more prone to deformation under the same stress conditions. The displacement of the surrounding rock increases by 5~10 times with each increase in the surrounding rock level. For example, the maximum deformation of Grade III surrounding rock is 12–17 cm, while the maximum deformation of Grade IV surrounding rock is 95–115 cm. This sharp increase in deformation not only increases the burden of the supporting structure but also may pose a threat to the overall safety of the project. In the relatively low stress range of 10–50 MPa, there is an approximate linear relationship between surrounding rock displacement and in situ stress. This means that with the gradual increase in in situ stress, the deformation of surrounding rock will also increase in proportion. However, when the local stress exceeds a certain critical value, this linear relationship may change, because the rock may enter a plastic or failure state at this time, resulting in a more complicated deformation mechanism. Although the effect of the cave diameter on the displacement of the adjacent rock is minor, it cannot be overlooked. As a special geological structure in underground space, the existence of a karst cave will change the stress distribution and deformation mode of the surrounding rock mass. When the diameter of the cave increases, the stress concentration around the cave intensifies, which may lead to local rock mass destruction and a displacement increase. In addition, factors such as the nature of the filler in the cave and the relative position between the cave and the tunnel may also have some influence on the deformation of the surrounding rock. When comparing the theoretical calculation results with the actual situation, the simulation results are often too large, and the arch waist displacement is generally greater than the vault displacement. This deviation comes from many aspects: first, theoretical models are usually based on a series of simplifications and assumptions, and it is difficult to fully reflect the complexity and variability of actual projects; second, the uncertainty of field geological conditions, the disturbance in the construction process, and the error of monitoring data may also lead to the difference between theory and practice; finally, the arch waist position is more prone to large deformation because of its complex stress and relatively weak supporting conditions.
Based on the above calculation results, the reserved deformation of the tunnel can be estimated in a more detailed macro way. This estimation not only depends on the specific surrounding rock grade, cave diameter, and in situ stress conditions but can also be combined with mathematical tools such as an interpolation method to accurately predict the maximum deformation under different conditions. According to the existing calculation results, including the surrounding rock deformation data under different surrounding rock grades, cave diameters, and in situ stresses, a multi-dimensional deformation prediction model is constructed. Check the corresponding calculation results and make predictions. It can effectively deal with continuously changing geological conditions and engineering parameters and improve the accuracy and reliability of reserved deformation estimation. When the in situ stress of the Class IV surrounding rock section is 45 MPa and the tunnel diameter is 13 m, these specific conditions can be input into the prediction model. The maximum deformation of the vault is 82 cm and the maximum horizontal displacement of the arch waist is 84 cm. This shows that under complex geological conditions, the actual deformation may significantly exceed the theoretical calculation results based on simplified assumptions. To confirm the validity of the prediction results, we compare and analyze the simulation results against the theoretical calculation values. As can be seen in
Table 2, the simulation results (such as vault 82 cm and arch waist 84 cm) have a higher assurance rate than the theoretical calculation value of 70 cm, which further proves the importance and application value of numerical simulation in complex geological engineering.
6. Conclusions
In this paper, the deformation characteristics of the diversion tunnel of Jinping II Hydropower Station are studied. Based on field statistics, theoretical calculation, and numerical simulation, considering the size, in situ stress, and grade of the surrounding rock, the reasonable reserved deformation of the deep soft rock tunnel is studied.
(1) The initial support deformation of the diversion tunnel is serious, although it does not directly cause large-scale cracking of the sprayed concrete layer, but its manifestations—the emergence of a cave wall bulge and local “prick” shape—fully show that the stability of the surrounding rock has been significantly affected. This deformation mode not only increases the stress burden of the supporting structure but also poses a threat to the subsequent construction safety and long-term stability of the tunnel. The intrusion of the pure green schist karst cave section into the design lining space is as high as 21.7% to 176.7%, which leads to the deformation of tunnel surrounding rock accounting for 20 to 60 cm of the design lining gap thickness. This intrusion and deformation not only reduce the effective thickness of the lining but also affect the interaction between the lining and surrounding rock, reducing the bearing capacity of the whole structure.
(2) The tunnel diameter, surrounding rock grade, in situ stress, and other factors will affect the tunnel deformation. On the premise of the reasonable selection of calculation parameters, the numerical calculation method is closer to the actual situation and can be used as the calculation method of reserved deformation. Various factors influence the reserved deformation, including the surrounding rock grade (lithology), in situ stress, and cave diameter.
(3) Considering the actual situation of the chlorite schist section at the west end of the Jinping diversion tunnel and the influence of support, in order to ensure the construction safety and structural stability, it is suggested to set the reserved deformation space of the tunnel as a flexible interval of 30 cm to 60 cm. In view of the direct influence of the distribution of the chlorite schist section on the deformation allowance, this paper puts forward some suggestions on optimizing the tunnel diameter design: adjusting the original planned diameter range from 13.4 m to 13.8 m to 14.6 m. Specifically, when the tunnel passes through many sections of chlorite schist, the reserved deformation should be appropriately increased to the upper limit of the section or higher to fully absorb the deformation of the surrounding rock. On the other hand, if the chlorite schist section accounts for a small proportion, the reserved deformation can be controlled near the lower limit of the section to realize a more economical and reasonable tunnel size design.