# Distribution Features of the Minimum Rock Cover Thickness of the Surrounding Rock Self-Stability of the Metro Tunnel in the Soil-Rock Dualistic Stratum

Authors: Ziguang Zhang; Tao Xu; Guangyong Cao; Xiaofeng Liu
Journal: Advances in Civil Engineering
 (2020)

---

## Abstract

Based on the maximum value of the self-stability of surrounding rock, the key issues affecting the distribution of the minimum mock cover thickness of the surrounding rock self-stability (MRCT-SRS) of the metro tunnel in the soil-rock dualistic stratum are discussed. The MRCT-SRS was studied for different soft stratum thicknessesHs and excavation spans D based on cusp catastrophe theory. Moreover, the 3D spatial distribution was determined, and the fitted curves were constructed according to the MRCT-SRS, Hs, and D. This facilitated the assessment of the stability of the surrounding rock in the excavation of a metro tunnel. Combined with the construction practices used for the Qingdao metro project, the validity of the obtained MRCT-SRS has been verified by the inverse analysis of the monitored data from the metro tunnel. Furthermore, the application of the MRCT-SRS in determining the reasonable buried depth of a metro line has been explored. The research results have provided theoretical support and scientific basis for the preparation, revision, and improvement of the relevant codes, standards, and guidelines for metro tunnel planning and design.


---

## Body

## 1. Introduction

The fundamental goal of an excellent metro tunnel design scheme is to make full use of the self-stability of the surrounding rock. Recently, many subsurface metro tunnels have been carried out in the soil-rock dualistic stratum in cities, such as Qingdao [1], Dalian [2], Guangzhou [3], Chongqing [4], and Guiyang [5]. When a metro tunnel is excavated in the urban area where the buried depth of the bedrock is shallow and the soft soil stratum thickness is thin, the tunnel should be placed in the rock stratum and maintained such that the minimum rock cover thickness meets the requirements of the self-stability of the rock surrounding. In this case, the construction of a metro tunnel can greatly reduce the tunnel support measures, save construction costs, reduce the tunnel construction safety risk and construction difficulty, and reduce the impact on the surrounding environment.Many scholars have studied the relationship between the stability of the surrounding rock and the rock cover thickness. Sun et al. [1] used an example of a metro station in Qingdao to determine the safety factor of a tunnel in surrounding rock at different buried depths by using the strength reduction method. They also considered the safety factor equal to 1 as the basis for judging the surrounding rock self-stability of the metro tunnel and determined the relationship between the minimum rock cover thickness and the soil stratum thickness. Yang and Xiao [6] established the cusp catastrophe model of the safety thickness of a karst tunnel floor and derived the expression for the critical safety thickness. Xu et al. [7] designed a study based on the Qingdao Jiaozhou Bay Subsea Tunnel using two mutually verified discrimination methods: the minimum vertical displacement of the tunnel vault and the maximum safety factor of the tunnel surrounding rock. They then studied the minimum rock cover thickness of the subsea tunnel and considered that the minimum rock cover thickness of a subsea tunnel is the reasonable and effective thickness of the tunnel. Wang et al. [8] pointed out that the reasonable buried depth of the metro station in the soil-rock dualistic stratum should be determined by using the double index system of the overburden thickness-span ratio and the rock cover thickness-span ratio and noted that the rock cover thickness-span ratio is more sensitive to the safety of tunnel construction than the overburden thickness-span ratio. Qiu et al. [9] studied the relationship between the minimum rock cover thickness and the surrounding rock stability of the subsea tunnel by the numerical method. Zhang et al. [10] discussed and expounded upon the issues of the surrounding rock self-stability of the metro tunnel in the soil-rock dualistic stratum from the three aspects of the stress, the displacement, and the safety factor. Wang et al. [11] studied the conversion timing of tunnel excavation method in upper soft and lower hard strata based on displacement direction angle theory under different working conditions.The catastrophe theory was put forward by Thom, a French scholar, in the 1970s [12] and was then expanded into a perfect system by Trotman and Zeeman [13], Zeeman et al. [14], and Poston and Stewart [15]. Catastrophe theory is a mathematical method used to study the phenomenon of discontinuity in nature. The cusp catastrophe model is one of seven primary catastrophe models in catastrophe theory, which has been widely used in the study of tunnel surrounding rock stability, and it has achieved ideal results [16–22]. Xia et al. [16] used cusp catastrophe theory and the discontinuous deformation analysis method to study the stability of the tunnel surrounding rock and obtain the safety factor. Based on the measured deformation data of the surrounding rock, Ren et al. [17] established the catastrophe model of the tunnel surrounding rock cusp by using catastrophe theory and deduced the instability criterion of the surrounding rock. Wang et al. [18] used the cusp catastrophe model to study the stability of gas storage pillars in layered salt caverns. Zhang and Han [19] used catastrophe theory to study the collapse mechanism and possible collapse block shape of a shallow, unlined tunnel. Zhou [20] et al. established a model based on the cusp catastrophe theory and analyzed the stability of the support system in the goaf of room pillar gypsum. Zhang et al. [21] and Huang and Zhang [22] used catastrophe theory to study the instability mechanism of a shallow tunnel.The problem of the minimum mock cover thickness of the surrounding rock self-stability (MRCT-SRS) of the tunnel has been studied by many scholars, who have used numerical calculations to analyze the stress, the deformation, and the safety factor. However, based on cusp catastrophe theory, rarely no relevant report has been found on the MRCT-SRS of the metro tunnel in the soil-rock dualistic stratum. This paper based on the calculation and analysis of the MRCT-SRS under different soft stratum thicknessHs and excavation span D was studied based on cusp catastrophe theory and numerical calculation, the 3D spatial distribution was drawn, and the fitted curves was constructed among MRCT-SRS, Hs, and D, which facilitated the assessment on the stability of the surrounding rock in excavation of metro tunnel. Combined with the construction practice of Qingdao metro project, the validity of the obtained MRCT-SRS has been verified by the inverse analysis on the data monitored for metro tunnel, and the application of MRCT-SRS in determining the reasonable buried depth of metro line has been expounded. The research results can provide theoretical support for the design of reasonable depth of the vertical section of metro tunnel in soil-rock dualistic stratum and further offer a scientific basis for the preparation, revision, and improvement of relevant codes, standards, and guidelines.


## 2. Methodology and Research Process

### 2.1. Cusp Catastrophe Theory

The cusp catastrophe model uses the simplest catastrophe mode with lag, divergence, and other behaviors, which is unique to the high-order catastrophe model. The critical interface is easy to construct and has a strong geometric rationale. Of the seven primary catastrophe models used in catastrophe theory, the cusp catastrophe model is the one with the most extensive application [23–29].The potential function V of the cusp catastrophe model is a 3D space (x, u, v) composed of the control variables u and v and the state variable x as follows:(1)Vx=x4+ux2+vx.An equation for the cusp catastrophe model for the balanced curved surface has been obtained after first deriving the potential function (1) with respect to x as follows:(2)Vx˙=4x3+2ux+v.The cusp catastrophe model uses a continuous balanced surface with a 3D phase space, the interior of which is folded into upper, middle, and lower pages, as shown in Figure1. Different areas have different numbers of balance points. On the middle page, the balance point falls on the maximum value of the potential function, so it is unstable. However, on the upper and lower pages, the potential function takes its minimum value, so the balance point is stable. The boundary between the upper and lower pages, on the one hand, and the middle page, on the other hand, is composed of points with vertical tangents, which is a set of mutation points as follows:(3)Vx¨=12x2+2u.Figure 1
Balance-curved surface and bifurcation set of the cusp catastrophe mode.Combining formulas (2) and (3) results in a bifurcation set, which is the projection of the fold of the balance-curved surface in the u-v plane as follows:(4)Δ=8u3+27v2.When applying catastrophe theory, scholars have used formula (4) as a direct criterion to determine the cusp mutation. When Δ > 0, the system is in a stable state; when Δ = 0, the system is in a critical equilibrium state; when Δ < 0, the system is in an unstable state [6].


### 2.2. Analysis Model

See Figure2, for the analysis model of the MRCT-SRS of the metro tunnel in the soil-rock dualistic stratum. In the figure, H represents the buried tunnel depth; D, the excavation span; h, the excavation height; and Hs and Hr, respectively, the soft soil stratum thickness and the rock cover thickness.Figure 2
Analysis model.A statistical analysis was performed of the geotechnical test results of 3227 geological boreholes in nearly 2000 sets of soil-rock strata along the Qingdao metro tunnel. Accordingly, the values of the physical and mechanical parameters of the soft soil stratum and the rock stratum in the soil-rock dualistic stratum structure are shown in Table1.Table 1
The physical and mechanical parameters.
Types of stratumCohesion (MPa)Internal friction angle (°)Modulus of elasticity (GPa)Poisson’s ratioDensity (kg/m3)Soft soil stratum0.032200.050.3822.5Rock stratum0.60355.00.2524.5


### 2.3. Target

The vault is the point at which the deflection is the largest around the cavern of the shallow metro tunnel. Therefore, the vertical displacement of the vault is the key indicator when evaluating the stability of the surrounding rock of the tunnel [30, 31]. If the rock cover thickness is large over the vault of the tunnel in the soil-rock dualistic stratum and it is gradually reduced, the vertical displacement of the vault at first is slowly reduced. Then, the displacement gradually increases after reaching its minimum as the thickness continues to decrease. When the thickness decreases to a certain value, the displacement sharply rises [7]. Currently, the vertical displacement of the tunnel vault is selected as the target for evaluating the self-stability of the surrounding rock of the tunnel. Furthermore, the corresponding rock cover thickness at the mutation point is used as the MRCT-SRS of the tunnel.


### 2.4. Process

For the MRCT-SRS of the metro tunnel in the soil-rock dualistic stratum based on featured point displacement in cusp catastrophe theory, the calculation is performed as follows:(1)
First, a numerical calculation model is constructed of the metro tunnel in the soil-rock dualistic stratum. In the numerical calculation, the stratum is considered according to the plane strain of ideal elastic-plastic materials, which should conform to the Mohr–Coulomb strength criterion and the large deformation hypothesis. ANSYS numerical calculation software and the inscribed circle DP4 yielded criterion based on the Mohr–Coulomb conditions, which were adopted for the one-time excavation of the whole section without considering the influence of groundwater and other factors. The upper boundary of the calculation model is set at the ground surface, and the lower boundary is set from the vault bottom of the excavation section down to 3.5 times the excavation height. Moreover, the distance between the left and right boundaries is more than 7 times the excavation span. The initial geostress incorporates not only the gravity stress of the rock and soil but also a uniformly distributed load of 20 kPa imposed on the ground surface. A large deformation mode is used during the calculation. The failure of the surrounding rock by large deformation can be reproduced by updating the coordinates. The vertical and horizontal ranges of the unit grid around the tunnel are 0.5 m, and those of the more peripheral unit grid are set at 1.0 m.(2)
By continuous calculations, the rock cover thickness is determined with respect to the smallest vault vertical displacement of the tunnel surrounding rock as the thickness gradually increases. Then, the thickness is reduced gradually at a specific interval (e.g., 1.0 m or 0.5 m) for multiple calculations. The vertical displacementsSi of the tunnel vault in each calculation are recorded to derive a sequential Si-Hr curve between the displacement Si and the rock cover thickness Hr. The number of calculations should increase in the mutation part of the Si-Hr curve.(3)
The calculation should begin starting with the rock cover thickness with respect to the smallest vault vertical displacement of the tunnel surrounding rock as the thickness gradually increases. The end interpolation method is used to fit the data of vertical displacementS and overburden rock thickness Hr, producing a quartic function as follows:(5)S=a0+a1Hr+a2Hr2+a3Hr3+a4Hr4,wherea0,a1,a2,a3,anda4 are undetermined coefficients.(4)
The fitting equation is transformed, and the standard potential function of the cusp mutation is derived asVx=x4+ux2+vx, for which u=a2/a4−3a32/8a42 and v=a1/a4−a2a3/2a42+a33/8a43.(5)
The bifurcation set equation is derived asΔ=8u3+27v2 according to cusp catastrophe theory.In this research,Hr with respect to Δ=8u3+27v2=0 is selected as the MRCT-SRS. Its calculation is shown in Figure 3.Figure 3
MRCT-SRS calculation process for metro tunnel.


### 2.5. Example Analysis

The MRCT-SRS calculation is demonstrated in detail in the following example of a metro tunnel with soft soil stratum thickness Hs of 12 m. The shape and dimensions of the tunnel section are shown in Figure4. The calculation results of the vertical displacement Si under different Hr values are shown in Table 2.Figure 4
Section size.Table 2
Calculation results of the vault vertical displacement.
Hr (m)1211109876S (mm)4.3254.1224.1604.2494.3504.5305.018Hr (m)54321.51.41S (mm)5.4376.1077.28710.22414.66016.28826.149The data ofS and Hr are used, so Hr is set to 11 m to perform the polynomial fitting. When Hr is 1 m, the fitting is carried out on 11 sets of data with Hr vales from 1 m to 11 m, producing the result as follows:(6)S=47.75224+28.81764Hr+6.90942Hr2+0.69678Hr3+0.02495Hr4.Substitute the values of the four parameters ofa1∼a4 from formula (6) into formulas (1) and (4). The result is u=−15.5395 and v=10.69658.Δ=−26929.91153<0, indicating that Hr does not meet the MRCT-SRS requirements.WhenHr is 2 m, Δ=1195349.531>0, indicating that it meets the requirements. Therefore, the MRCT-SRS should be between 1 m and 2 m.To improve further the calculation accuracy of the MRCT-SRS,Hr is set as 1.5 m. The result shows Δ=7211.982>0, indicating that the Hr meets the requirements; when Hr is 1.4, Δ=−3072.25<0, indicating that it does not meet the requirements. Therefore, the MRCT-SRS should be between 1.4 m and 1.5 m. The MRCT-SRS is finally set at 1.45 m when the soft soil stratum thickness is 12 m in this paper.


## 2.1. Cusp Catastrophe Theory

The cusp catastrophe model uses the simplest catastrophe mode with lag, divergence, and other behaviors, which is unique to the high-order catastrophe model. The critical interface is easy to construct and has a strong geometric rationale. Of the seven primary catastrophe models used in catastrophe theory, the cusp catastrophe model is the one with the most extensive application [23–29].The potential function V of the cusp catastrophe model is a 3D space (x, u, v) composed of the control variables u and v and the state variable x as follows:(1)Vx=x4+ux2+vx.An equation for the cusp catastrophe model for the balanced curved surface has been obtained after first deriving the potential function (1) with respect to x as follows:(2)Vx˙=4x3+2ux+v.The cusp catastrophe model uses a continuous balanced surface with a 3D phase space, the interior of which is folded into upper, middle, and lower pages, as shown in Figure1. Different areas have different numbers of balance points. On the middle page, the balance point falls on the maximum value of the potential function, so it is unstable. However, on the upper and lower pages, the potential function takes its minimum value, so the balance point is stable. The boundary between the upper and lower pages, on the one hand, and the middle page, on the other hand, is composed of points with vertical tangents, which is a set of mutation points as follows:(3)Vx¨=12x2+2u.Figure 1
Balance-curved surface and bifurcation set of the cusp catastrophe mode.Combining formulas (2) and (3) results in a bifurcation set, which is the projection of the fold of the balance-curved surface in the u-v plane as follows:(4)Δ=8u3+27v2.When applying catastrophe theory, scholars have used formula (4) as a direct criterion to determine the cusp mutation. When Δ > 0, the system is in a stable state; when Δ = 0, the system is in a critical equilibrium state; when Δ < 0, the system is in an unstable state [6].


## 2.2. Analysis Model

See Figure2, for the analysis model of the MRCT-SRS of the metro tunnel in the soil-rock dualistic stratum. In the figure, H represents the buried tunnel depth; D, the excavation span; h, the excavation height; and Hs and Hr, respectively, the soft soil stratum thickness and the rock cover thickness.Figure 2
Analysis model.A statistical analysis was performed of the geotechnical test results of 3227 geological boreholes in nearly 2000 sets of soil-rock strata along the Qingdao metro tunnel. Accordingly, the values of the physical and mechanical parameters of the soft soil stratum and the rock stratum in the soil-rock dualistic stratum structure are shown in Table1.Table 1
The physical and mechanical parameters.
Types of stratumCohesion (MPa)Internal friction angle (°)Modulus of elasticity (GPa)Poisson’s ratioDensity (kg/m3)Soft soil stratum0.032200.050.3822.5Rock stratum0.60355.00.2524.5


## 2.3. Target

The vault is the point at which the deflection is the largest around the cavern of the shallow metro tunnel. Therefore, the vertical displacement of the vault is the key indicator when evaluating the stability of the surrounding rock of the tunnel [30, 31]. If the rock cover thickness is large over the vault of the tunnel in the soil-rock dualistic stratum and it is gradually reduced, the vertical displacement of the vault at first is slowly reduced. Then, the displacement gradually increases after reaching its minimum as the thickness continues to decrease. When the thickness decreases to a certain value, the displacement sharply rises [7]. Currently, the vertical displacement of the tunnel vault is selected as the target for evaluating the self-stability of the surrounding rock of the tunnel. Furthermore, the corresponding rock cover thickness at the mutation point is used as the MRCT-SRS of the tunnel.


## 2.4. Process

For the MRCT-SRS of the metro tunnel in the soil-rock dualistic stratum based on featured point displacement in cusp catastrophe theory, the calculation is performed as follows:(1)
First, a numerical calculation model is constructed of the metro tunnel in the soil-rock dualistic stratum. In the numerical calculation, the stratum is considered according to the plane strain of ideal elastic-plastic materials, which should conform to the Mohr–Coulomb strength criterion and the large deformation hypothesis. ANSYS numerical calculation software and the inscribed circle DP4 yielded criterion based on the Mohr–Coulomb conditions, which were adopted for the one-time excavation of the whole section without considering the influence of groundwater and other factors. The upper boundary of the calculation model is set at the ground surface, and the lower boundary is set from the vault bottom of the excavation section down to 3.5 times the excavation height. Moreover, the distance between the left and right boundaries is more than 7 times the excavation span. The initial geostress incorporates not only the gravity stress of the rock and soil but also a uniformly distributed load of 20 kPa imposed on the ground surface. A large deformation mode is used during the calculation. The failure of the surrounding rock by large deformation can be reproduced by updating the coordinates. The vertical and horizontal ranges of the unit grid around the tunnel are 0.5 m, and those of the more peripheral unit grid are set at 1.0 m.(2)
By continuous calculations, the rock cover thickness is determined with respect to the smallest vault vertical displacement of the tunnel surrounding rock as the thickness gradually increases. Then, the thickness is reduced gradually at a specific interval (e.g., 1.0 m or 0.5 m) for multiple calculations. The vertical displacementsSi of the tunnel vault in each calculation are recorded to derive a sequential Si-Hr curve between the displacement Si and the rock cover thickness Hr. The number of calculations should increase in the mutation part of the Si-Hr curve.(3)
The calculation should begin starting with the rock cover thickness with respect to the smallest vault vertical displacement of the tunnel surrounding rock as the thickness gradually increases. The end interpolation method is used to fit the data of vertical displacementS and overburden rock thickness Hr, producing a quartic function as follows:(5)S=a0+a1Hr+a2Hr2+a3Hr3+a4Hr4,wherea0,a1,a2,a3,anda4 are undetermined coefficients.(4)
The fitting equation is transformed, and the standard potential function of the cusp mutation is derived asVx=x4+ux2+vx, for which u=a2/a4−3a32/8a42 and v=a1/a4−a2a3/2a42+a33/8a43.(5)
The bifurcation set equation is derived asΔ=8u3+27v2 according to cusp catastrophe theory.In this research,Hr with respect to Δ=8u3+27v2=0 is selected as the MRCT-SRS. Its calculation is shown in Figure 3.Figure 3
MRCT-SRS calculation process for metro tunnel.


## 2.5. Example Analysis

The MRCT-SRS calculation is demonstrated in detail in the following example of a metro tunnel with soft soil stratum thickness Hs of 12 m. The shape and dimensions of the tunnel section are shown in Figure4. The calculation results of the vertical displacement Si under different Hr values are shown in Table 2.Figure 4
Section size.Table 2
Calculation results of the vault vertical displacement.
Hr (m)1211109876S (mm)4.3254.1224.1604.2494.3504.5305.018Hr (m)54321.51.41S (mm)5.4376.1077.28710.22414.66016.28826.149The data ofS and Hr are used, so Hr is set to 11 m to perform the polynomial fitting. When Hr is 1 m, the fitting is carried out on 11 sets of data with Hr vales from 1 m to 11 m, producing the result as follows:(6)S=47.75224+28.81764Hr+6.90942Hr2+0.69678Hr3+0.02495Hr4.Substitute the values of the four parameters ofa1∼a4 from formula (6) into formulas (1) and (4). The result is u=−15.5395 and v=10.69658.Δ=−26929.91153<0, indicating that Hr does not meet the MRCT-SRS requirements.WhenHr is 2 m, Δ=1195349.531>0, indicating that it meets the requirements. Therefore, the MRCT-SRS should be between 1 m and 2 m.To improve further the calculation accuracy of the MRCT-SRS,Hr is set as 1.5 m. The result shows Δ=7211.982>0, indicating that the Hr meets the requirements; when Hr is 1.4, Δ=−3072.25<0, indicating that it does not meet the requirements. Therefore, the MRCT-SRS should be between 1.4 m and 1.5 m. The MRCT-SRS is finally set at 1.45 m when the soft soil stratum thickness is 12 m in this paper.


## 3. Distribution Features of the MRCT-SRS

### 3.1. MRCT-SRSs of the Metro Station Tunnel and the Running Tunnel

The example of the engineering of the Qingdao metro tunnel is used in this research to investigate the MRCT-SRS of the station tunnel and the running tunnel with two main types of sections in a soil-rock dualistic stratum (see Figure5).Figure 5
Metro tunnel section type and size.The results of the MRCT-SRS calculations for the metro station tunnel and the running tunnel under different soft soil stratum thicknessesHs are shown in Table 3 and Figure 6. The MRCT-SRS of the metro tunnel approximately linearly increases or decreases with respect to Hs. The linearity fitting is used to predict the MRCT-SRSs of the metro station and the running tunnel. The fitted equation for the station tunnel MRCT-SRS and Hs is shown in formula (7). The adjusted R-square value is 0.9991, and the MRCT-SRS is approximately linear along with Hs. The fitted equation of the running tunnel MRCT-SRS and Hs is shown in formula (8), and the adjusted R-square value is 0.99258. With formulas (7) and (8), the MRCT-SRSs for the station tunnel and the running tunnel can be derived with respect to different Hs.Table 3
MRCT-SRS calculation results for metro tunnel.
Hs (m)36912151821242730MRCT-SRS (m)Station tunnel0.620.921.171.451.762.022.282.582.923.12Runing tunnel0.200.340.450.520.590.710.780.900.971.13Figure 6
Metro tunnel MRCT-SRSs change with respect to Hs. (a) Station tunnel. (b) Running tunnel.
(a)(b)The MRCT-SRS of the metro tunnel increases approximately linearly with the increase inHs. This result occurs because, when excavating the metro tunnel in the soil-rock dualistic stratum under the requirements of the MRCT-SRS, the overlying soft soil stratum bears a higher load, while the rock cover stratum bears more of the structure. The rock mechanics performance around the tunnel is similar to the force mode of a composite arch or a composite beam. Therefore, when the Hs is greater, the corresponding load and the MRCT-SRS become larger.(7)Station tunnel:MRCT−SRS=0.0935Hs+0.3406,(8)Running tunnel:MRCT−SRS=0.03218Hs+0.128.


### 3.2. Discussion of How MRCT-SRS Works for the Metro Tunnel

The stability characteristics of the tunnel surrounding rock are closely related to the shape and size of the tunnel section. By focusing on the metro station tunnel, the general rule of how the MRCT-SRS varies with the excavation spanD is discussed in the present study while maintaining the original section shape and high-span ratio but using only the scaling tunnel section size. The tunnel section size is shown in Figure 7. The MRCT-SRS calculation results of the metro tunnel are shown in Table 4 when excavation spans D are 10.4 m and 31.2 m with respect to Hs.Figure 7
Metro tunnel excavation section.Table 4
MRCT-SRS for the metro tunnel.
Hs (m)612182430MRCT-SRS (m)D = 10.4 m0.310.781.271.481.83D = 31.2 m1.232.243.063.774.42The 3D spatial distribution of MRCT-SRS ∼ Hs ∼ D among the metro tunnel MRCT-SRS, soft soil stratum thicknessHs, and excavation span D is shown in Figure 8. The figure overly shows the distribution characteristics of the metro tunnel MRCT-SRS with respect to different Hs and D. At the same time, according to the figure, the surrounding rock stability and supplementary measures in the excavation for the metro tunnel in the soil-rock dualistic stratum can be determined as follows:(1)
WhenHr falls below the interface plane, the excavation for the surrounding rock cannot meet the self-stability condition. Strong supplementary engineering measures are needed to prevent its instability to ensure normal tunnel excavation.(2)
WhenHr falls on the interface plane, the excavation of the surrounding rock just meets the self-stability condition. Some supplementary engineering measures are needed to prevent its instability to ensure normal tunnel excavation.(3)
WhenHr falls above the interface plane, the excavation for the surrounding rock can meet the self-stability condition. Only simple supplementary engineering measures are needed to prevent the geometric instability of rocks to ensure normal tunnel excavation.Figure 8
3D MRCT-SRS ∼ Hs ∼ D diagram for metro tunnel.The fitting equation of MRCT-SRS ∼ Hs ∼ Dis shown as formula (9), and the adjusted R-square is 0.99597. By this formula, the MRCT-SRS with respect to different Hs and D can be derived:(9)MRCT−SRS=−0.59406+6.8447exp−D−51.8360221742.17836−Hs−38.3197521037.83038.


## 3.1. MRCT-SRSs of the Metro Station Tunnel and the Running Tunnel

The example of the engineering of the Qingdao metro tunnel is used in this research to investigate the MRCT-SRS of the station tunnel and the running tunnel with two main types of sections in a soil-rock dualistic stratum (see Figure5).Figure 5
Metro tunnel section type and size.The results of the MRCT-SRS calculations for the metro station tunnel and the running tunnel under different soft soil stratum thicknessesHs are shown in Table 3 and Figure 6. The MRCT-SRS of the metro tunnel approximately linearly increases or decreases with respect to Hs. The linearity fitting is used to predict the MRCT-SRSs of the metro station and the running tunnel. The fitted equation for the station tunnel MRCT-SRS and Hs is shown in formula (7). The adjusted R-square value is 0.9991, and the MRCT-SRS is approximately linear along with Hs. The fitted equation of the running tunnel MRCT-SRS and Hs is shown in formula (8), and the adjusted R-square value is 0.99258. With formulas (7) and (8), the MRCT-SRSs for the station tunnel and the running tunnel can be derived with respect to different Hs.Table 3
MRCT-SRS calculation results for metro tunnel.
Hs (m)36912151821242730MRCT-SRS (m)Station tunnel0.620.921.171.451.762.022.282.582.923.12Runing tunnel0.200.340.450.520.590.710.780.900.971.13Figure 6
Metro tunnel MRCT-SRSs change with respect to Hs. (a) Station tunnel. (b) Running tunnel.
(a)(b)The MRCT-SRS of the metro tunnel increases approximately linearly with the increase inHs. This result occurs because, when excavating the metro tunnel in the soil-rock dualistic stratum under the requirements of the MRCT-SRS, the overlying soft soil stratum bears a higher load, while the rock cover stratum bears more of the structure. The rock mechanics performance around the tunnel is similar to the force mode of a composite arch or a composite beam. Therefore, when the Hs is greater, the corresponding load and the MRCT-SRS become larger.(7)Station tunnel:MRCT−SRS=0.0935Hs+0.3406,(8)Running tunnel:MRCT−SRS=0.03218Hs+0.128.


## 3.2. Discussion of How MRCT-SRS Works for the Metro Tunnel

The stability characteristics of the tunnel surrounding rock are closely related to the shape and size of the tunnel section. By focusing on the metro station tunnel, the general rule of how the MRCT-SRS varies with the excavation spanD is discussed in the present study while maintaining the original section shape and high-span ratio but using only the scaling tunnel section size. The tunnel section size is shown in Figure 7. The MRCT-SRS calculation results of the metro tunnel are shown in Table 4 when excavation spans D are 10.4 m and 31.2 m with respect to Hs.Figure 7
Metro tunnel excavation section.Table 4
MRCT-SRS for the metro tunnel.
Hs (m)612182430MRCT-SRS (m)D = 10.4 m0.310.781.271.481.83D = 31.2 m1.232.243.063.774.42The 3D spatial distribution of MRCT-SRS ∼ Hs ∼ D among the metro tunnel MRCT-SRS, soft soil stratum thicknessHs, and excavation span D is shown in Figure 8. The figure overly shows the distribution characteristics of the metro tunnel MRCT-SRS with respect to different Hs and D. At the same time, according to the figure, the surrounding rock stability and supplementary measures in the excavation for the metro tunnel in the soil-rock dualistic stratum can be determined as follows:(1)
WhenHr falls below the interface plane, the excavation for the surrounding rock cannot meet the self-stability condition. Strong supplementary engineering measures are needed to prevent its instability to ensure normal tunnel excavation.(2)
WhenHr falls on the interface plane, the excavation of the surrounding rock just meets the self-stability condition. Some supplementary engineering measures are needed to prevent its instability to ensure normal tunnel excavation.(3)
WhenHr falls above the interface plane, the excavation for the surrounding rock can meet the self-stability condition. Only simple supplementary engineering measures are needed to prevent the geometric instability of rocks to ensure normal tunnel excavation.Figure 8
3D MRCT-SRS ∼ Hs ∼ D diagram for metro tunnel.The fitting equation of MRCT-SRS ∼ Hs ∼ Dis shown as formula (9), and the adjusted R-square is 0.99597. By this formula, the MRCT-SRS with respect to different Hs and D can be derived:(9)MRCT−SRS=−0.59406+6.8447exp−D−51.8360221742.17836−Hs−38.3197521037.83038.


## 4. Case Analysis and Engineering Application

### 4.1. Project Overview

#### 4.1.1. Project Descriptions

The planned metro lines for Qingdao’s rail transit are approximately 807 km long and occur in 16 lines. Qingdao Metro Line 3, the first urban rail transit construction, is approximately 24.8 km long and includes 22 stations. Its construction is followed by Line 2, covering a length of approximately 25.2 km with 22 stations. The two lines have 16 underground excavated stations. Except for Dunhua Road Station in the tower-column structure, the other excavated stations all adopt a single-arch large-span structure with an excavation span of 19.2 m∼23.8 m (see Figure9 and Table 5). The metro running tunnel mainly adopts the horseshoe-shaped subterranean excavation form with a single span of 5.8 m∼6.2 m for a one-way single lane.Figure 9
Qingdao metro lines 2 and 3 plane map.Table 5
Qingdao metro station tunnel in excavation.
S/NStation nameExcavation techniqueD (m) × h (m)H (m)1Qingdao Railway StationOpen-cut + CRD23.45 × 14.582.3 ∼ 8.52Zhongshan ParkArch cover method19.20 × 16.2010.0 ∼ 12.03ZhanshanArch cover + CRD20.52 × 18.1813.0 ∼ 15.14Jiangxi RdArch cover + CRD20.60 × 15.508.9 ∼ 10.55Dunhua RdBench method10.38 × 9.5014.6 ∼ 20.06Qingjiang RdArch cover + CRD20.36 × 17.356.3 ∼ 8.37Wannianquan RdArch cover + CRD20.66 × 18.165.1 ∼ 10.18Junfeng RdArch cover + CRD20.79 × 18.377.4 ∼ 15.69Yan’anlu RdArch cover method23.76 × 19.1313.4 ∼ 16.410Zhiquan RdArch cover method23.76 × 19.1313.2 ∼ 17.011Haichuan RdStep + CD22.66 × 16.8013.3 ∼ 19.012Xujia MaidaoStep + CD22.10 × 16.8013.3 ∼ 18.013Tong’an RdArch cover method19.80 × 16.508.96 ∼ 10.314East Bus StationArch cover method22.96 × 16.9910.8 ∼ 11.7515Zaoshan RdArch cover method19.80 × 16.5010.17 ∼ 16.2216Licun ParkArch cover + CRD24.40 × 16.8010.0 ∼ 13.0NoteCRD means cross diagram method, and CD means center diagram method


#### 4.1.2. Stratigraphic Features

The stratum distribution along Qingdao Metro from top to bottom includes the Quaternary stratum as well as the strongly weathered, the moderately weathered, and the microweathered (or nonweathered) granite strata. To simplify the analysis and highlight the main points, the Quaternary stratum and the strong weathering stratum were considered soft soil stratum, and the strata in the middle weathering strata and below are considered rock stratum, thus constructing a distribution model for the soil-rock dualistic stratum. The structural characteristics of the dualistic stratum along the two metro lines in Qingdao are shown in Figure10.Figure 10
Distribution map of the dualistic stratum structure along the lines.According to the statistical analysis of 3227 geological borehole samples along the two metro lines, the average thickness of the soft soil stratum is 10.2 m; the maximum is 38.0 m; the median is 10.0 m; and the lower and upper quartiles are, respectively, 5.0 m and 14.8 m (as shown in Figure11 and Table 6). Further analysis shows that the number of borehole samples with soft soil stratum thicknesses of 9 m or less along the metro lines accounts for 45.6% of the samples; those with less than 12 m, 61.4%; those with less than 15 m, 77.5%; and those with less than 18 m, 89.7% (as shown in Table 7).Figure 11
Statistical box diagram of the soft soil thickness along lines.Table 6
Statistics of the soft soil thickness.
Min. (m)Low position (m)Lower quartile (m)Median (m)Upper quartile (m)High position (m)Max. (m)Average (m)SD (m)0.00.05.010.014.829.838.010.26.35Table 7
Statistical table of the soft soil thickness along lines.
Soil thickness (m)≤3 m≤6 m≤9 m≤12 m≤15 m≤18 mPercentage15.6%30.7%45.6%61.4%77.5%89.7%


### 4.2. Analysis of Metro Tunnel Monitoring Data

#### 4.2.1. Data Resource

The monitoring data come from six station tunnels of the two metro lines, which have approximately the same section size, excavation depth and technique, support measures, etc. The number of monitoring points, the soft soil stratum thicknessHs, the rock cover thickness Hr, and the MRCT-SRS are shown in Table 8. The statistical data are from a third party. The monitoring points for the ground surface settlement are selected within the main area of the tunnel, among which the abnormal points should be excluded in the statistics.Table 8
Basic overview of metro station tunnel.
Serial numberNameNumber of measuring points (number)Hs (m)Hr (m)MRCT-SRS (m)Distribution rangeMean valueDistribution rangeMean value1Zhanshan Road Station4113.0 ∼ 20.516.80.00.01.80–2.502Jiangxi Road Station689.1 ∼ 1411.50.00.01.12–2.103Zhongshan Park618.5 ∼ 16.511.60.0 ∼ 1.60.81.10–2.254Junfeng Road Station190.7 ∼ 6.83.24.7 ∼ 9.87.30.03–0.805Haichuan Road Station280.8 ∼ 8.34.37.0 ∼ 119.00.03–0.906Zhiquan Road Station196.1 ∼ 16.58.98.2 ∼ 14.311.20.50–0.90


### 4.3. Inversion Analysis

Table9 shows that, with regard to the overburden thicknesses exceeding the MRCT-SRS requirements for the stations of Zhiquan Road, Haichuan Road, and Junfeng Road, the average surface settlement is between 2.29 mm and 6.37 mm, with a maximum value between 6.22 mm and 12.59 mm, which are small in general. Regarding those for the stations of Zhongshan Park, Jiangxi Road, and Zhanshan Road, the average surface settlement is between 37.5 mm and 54.75 mm, with a maximum value between 71.54 mm and 80.96 mm, which are large in general.Table 9
Statistical table of ground surface settlement of the station tunnel.
NameScope (mm)Mean value (mm)Coefficient of variationStandard deviation (mm)The size relationship between Hr and MRCT-SRSZhiquan Road Station1.37 ∼ 6.222.290.361.18Hr > MRCT-SRSHaichuan Road Station0.18 ∼ 9.974.591.872.45Hr > MRCT-SRSJunfeng Road Station2.03 ∼ 12.596.372.093.05Hr > MRCT-SRSZhongshan Park Station17.31 ∼ 71.5437.580.3513.30Hr < MRCT-SRSJiangxi Road Station16.79 ∼ 76.9846.330.2712.69Hr < MRCT-SRSZhanshan Road Station21.19 ∼ 80.9654.750.3116.97Hr < MRCT-SRSSurface settlement can comprehensively reflect the effect of metro construction quality control and the degree of construction influence on the surrounding rock as well as the soil and neighboring environment. For a metro tunnel with shallow buried depth under the approximately same conditions of section size, excavation depth and technique, support measures, etc., the rock cover thickness on the vault meets the MRCT-SRS requirements or not, the surface settlement caused by the tunnel excavation is different by an order of magnitude. This situation fully demonstrates that the self-stability of the surrounding rock has a good correlation with the rock cover thickness (as shown in Figure12). The rationality of the conclusions obtained in this study has been fully verified.Figure 12
The change of ground settlement with different rock cover thicknesses.


### 4.4. Engineering Application

The fundamental goal of an excellent metro line design is to make full use of the self-stability of the surrounding rock and to minimize supplementary engineering measures. Costs should be considered only when engineering security is guaranteed. Furthermore, geological risks and construction difficulty should be reduced. The degree to which the surrounding rock can be self-stabilized depends mainly on the distribution conditions of the strata along the metro lines, the understanding of these conditions, and the ability to execute the required tasks. Therefore, the Geologically Easily Stable Area (GSEA) in this study refers to the area that meets the MRCT-SRS requirements for the excavation of the tunnel along the metro in the soil-rock dualistic stratum, otherwise referred to as the Geologically Difficultly Stable Area (GSDA).The metro station and the running tunnel are considered, respectively, as a single-arch large-span excavation structure with the excavation span of 20.8 m and a horseshoe-shaped structure with the excavation span of 6.2 m. The buried depth is set at equal intervals of 3 m falling within the scope of 6 m∼30 m. Table10 and Figure 13 show the calculation results of the proportions of the GSEA lengths of the station (considering only 16 subterranean excavation stations) and the running tunnels, respectively, as well as the lengths of the station and the running tunnel of Qingdao’s two metro lines under different buried depths.Table 10
The GSEA ratio of metro tunnel under different buried depths.
H (m)6912151821242730Station tunnel (%)49.265.076.484.190.794.398.099.0100.0Running tunel (%)26.340.557.071.687.593.996.698.599.3Figure 13
GSEA proportion variation for metro tunnel under different buried depths.As the buried depth increases, the proportion of the GSEA in the station tunnel rises first rapidly but then slowly. When the buried depth of the station tunnel is at 6 m, the proportion is 49.2%. When the depth increases to 12 m, the proportion increases rapidly to 76.4%. However, when the depth reaches 18 m, the proportion increase to 90.7%. The change of GSEA with buried depth in running tunnel is similar to that in station tunnel. When the buried depth of the running tunnel is 6 m, the proportion is only 26.3%. However, when the depth increases to 15 m, the proportion rapidly increases to 71.6%. As the depth continues to rise to 18 m, the proportion increases to 87.5%. When the depth reaches 21 m, the proportion slowly increases to 93.9%.The metro line is mainly composed of the station and the running tunnel. The buried depth of the line depends on the depth of the station tunnel. The depth of the running tunnel is mainly determined by factors such as the buried depth of the station tunnel and the longitudinal slope. The height difference between the station and the running tunnels is set as 12 m. Therefore, when the buried depth of the Qingdao metro station falls within 12m∼15 m, the GSEA of the station tunnel accounts for 76.4%∼84.1%, and the GSEA of the running tunnel accounts for 96.6%∼98.5%.


## 4.1. Project Overview

### 4.1.1. Project Descriptions

The planned metro lines for Qingdao’s rail transit are approximately 807 km long and occur in 16 lines. Qingdao Metro Line 3, the first urban rail transit construction, is approximately 24.8 km long and includes 22 stations. Its construction is followed by Line 2, covering a length of approximately 25.2 km with 22 stations. The two lines have 16 underground excavated stations. Except for Dunhua Road Station in the tower-column structure, the other excavated stations all adopt a single-arch large-span structure with an excavation span of 19.2 m∼23.8 m (see Figure9 and Table 5). The metro running tunnel mainly adopts the horseshoe-shaped subterranean excavation form with a single span of 5.8 m∼6.2 m for a one-way single lane.Figure 9
Qingdao metro lines 2 and 3 plane map.Table 5
Qingdao metro station tunnel in excavation.
S/NStation nameExcavation techniqueD (m) × h (m)H (m)1Qingdao Railway StationOpen-cut + CRD23.45 × 14.582.3 ∼ 8.52Zhongshan ParkArch cover method19.20 × 16.2010.0 ∼ 12.03ZhanshanArch cover + CRD20.52 × 18.1813.0 ∼ 15.14Jiangxi RdArch cover + CRD20.60 × 15.508.9 ∼ 10.55Dunhua RdBench method10.38 × 9.5014.6 ∼ 20.06Qingjiang RdArch cover + CRD20.36 × 17.356.3 ∼ 8.37Wannianquan RdArch cover + CRD20.66 × 18.165.1 ∼ 10.18Junfeng RdArch cover + CRD20.79 × 18.377.4 ∼ 15.69Yan’anlu RdArch cover method23.76 × 19.1313.4 ∼ 16.410Zhiquan RdArch cover method23.76 × 19.1313.2 ∼ 17.011Haichuan RdStep + CD22.66 × 16.8013.3 ∼ 19.012Xujia MaidaoStep + CD22.10 × 16.8013.3 ∼ 18.013Tong’an RdArch cover method19.80 × 16.508.96 ∼ 10.314East Bus StationArch cover method22.96 × 16.9910.8 ∼ 11.7515Zaoshan RdArch cover method19.80 × 16.5010.17 ∼ 16.2216Licun ParkArch cover + CRD24.40 × 16.8010.0 ∼ 13.0NoteCRD means cross diagram method, and CD means center diagram method


### 4.1.2. Stratigraphic Features

The stratum distribution along Qingdao Metro from top to bottom includes the Quaternary stratum as well as the strongly weathered, the moderately weathered, and the microweathered (or nonweathered) granite strata. To simplify the analysis and highlight the main points, the Quaternary stratum and the strong weathering stratum were considered soft soil stratum, and the strata in the middle weathering strata and below are considered rock stratum, thus constructing a distribution model for the soil-rock dualistic stratum. The structural characteristics of the dualistic stratum along the two metro lines in Qingdao are shown in Figure10.Figure 10
Distribution map of the dualistic stratum structure along the lines.According to the statistical analysis of 3227 geological borehole samples along the two metro lines, the average thickness of the soft soil stratum is 10.2 m; the maximum is 38.0 m; the median is 10.0 m; and the lower and upper quartiles are, respectively, 5.0 m and 14.8 m (as shown in Figure11 and Table 6). Further analysis shows that the number of borehole samples with soft soil stratum thicknesses of 9 m or less along the metro lines accounts for 45.6% of the samples; those with less than 12 m, 61.4%; those with less than 15 m, 77.5%; and those with less than 18 m, 89.7% (as shown in Table 7).Figure 11
Statistical box diagram of the soft soil thickness along lines.Table 6
Statistics of the soft soil thickness.
Min. (m)Low position (m)Lower quartile (m)Median (m)Upper quartile (m)High position (m)Max. (m)Average (m)SD (m)0.00.05.010.014.829.838.010.26.35Table 7
Statistical table of the soft soil thickness along lines.
Soil thickness (m)≤3 m≤6 m≤9 m≤12 m≤15 m≤18 mPercentage15.6%30.7%45.6%61.4%77.5%89.7%


## 4.1.1. Project Descriptions

The planned metro lines for Qingdao’s rail transit are approximately 807 km long and occur in 16 lines. Qingdao Metro Line 3, the first urban rail transit construction, is approximately 24.8 km long and includes 22 stations. Its construction is followed by Line 2, covering a length of approximately 25.2 km with 22 stations. The two lines have 16 underground excavated stations. Except for Dunhua Road Station in the tower-column structure, the other excavated stations all adopt a single-arch large-span structure with an excavation span of 19.2 m∼23.8 m (see Figure9 and Table 5). The metro running tunnel mainly adopts the horseshoe-shaped subterranean excavation form with a single span of 5.8 m∼6.2 m for a one-way single lane.Figure 9
Qingdao metro lines 2 and 3 plane map.Table 5
Qingdao metro station tunnel in excavation.
S/NStation nameExcavation techniqueD (m) × h (m)H (m)1Qingdao Railway StationOpen-cut + CRD23.45 × 14.582.3 ∼ 8.52Zhongshan ParkArch cover method19.20 × 16.2010.0 ∼ 12.03ZhanshanArch cover + CRD20.52 × 18.1813.0 ∼ 15.14Jiangxi RdArch cover + CRD20.60 × 15.508.9 ∼ 10.55Dunhua RdBench method10.38 × 9.5014.6 ∼ 20.06Qingjiang RdArch cover + CRD20.36 × 17.356.3 ∼ 8.37Wannianquan RdArch cover + CRD20.66 × 18.165.1 ∼ 10.18Junfeng RdArch cover + CRD20.79 × 18.377.4 ∼ 15.69Yan’anlu RdArch cover method23.76 × 19.1313.4 ∼ 16.410Zhiquan RdArch cover method23.76 × 19.1313.2 ∼ 17.011Haichuan RdStep + CD22.66 × 16.8013.3 ∼ 19.012Xujia MaidaoStep + CD22.10 × 16.8013.3 ∼ 18.013Tong’an RdArch cover method19.80 × 16.508.96 ∼ 10.314East Bus StationArch cover method22.96 × 16.9910.8 ∼ 11.7515Zaoshan RdArch cover method19.80 × 16.5010.17 ∼ 16.2216Licun ParkArch cover + CRD24.40 × 16.8010.0 ∼ 13.0NoteCRD means cross diagram method, and CD means center diagram method


## 4.1.2. Stratigraphic Features

The stratum distribution along Qingdao Metro from top to bottom includes the Quaternary stratum as well as the strongly weathered, the moderately weathered, and the microweathered (or nonweathered) granite strata. To simplify the analysis and highlight the main points, the Quaternary stratum and the strong weathering stratum were considered soft soil stratum, and the strata in the middle weathering strata and below are considered rock stratum, thus constructing a distribution model for the soil-rock dualistic stratum. The structural characteristics of the dualistic stratum along the two metro lines in Qingdao are shown in Figure10.Figure 10
Distribution map of the dualistic stratum structure along the lines.According to the statistical analysis of 3227 geological borehole samples along the two metro lines, the average thickness of the soft soil stratum is 10.2 m; the maximum is 38.0 m; the median is 10.0 m; and the lower and upper quartiles are, respectively, 5.0 m and 14.8 m (as shown in Figure11 and Table 6). Further analysis shows that the number of borehole samples with soft soil stratum thicknesses of 9 m or less along the metro lines accounts for 45.6% of the samples; those with less than 12 m, 61.4%; those with less than 15 m, 77.5%; and those with less than 18 m, 89.7% (as shown in Table 7).Figure 11
Statistical box diagram of the soft soil thickness along lines.Table 6
Statistics of the soft soil thickness.
Min. (m)Low position (m)Lower quartile (m)Median (m)Upper quartile (m)High position (m)Max. (m)Average (m)SD (m)0.00.05.010.014.829.838.010.26.35Table 7
Statistical table of the soft soil thickness along lines.
Soil thickness (m)≤3 m≤6 m≤9 m≤12 m≤15 m≤18 mPercentage15.6%30.7%45.6%61.4%77.5%89.7%


## 4.2. Analysis of Metro Tunnel Monitoring Data

### 4.2.1. Data Resource

The monitoring data come from six station tunnels of the two metro lines, which have approximately the same section size, excavation depth and technique, support measures, etc. The number of monitoring points, the soft soil stratum thicknessHs, the rock cover thickness Hr, and the MRCT-SRS are shown in Table 8. The statistical data are from a third party. The monitoring points for the ground surface settlement are selected within the main area of the tunnel, among which the abnormal points should be excluded in the statistics.Table 8
Basic overview of metro station tunnel.
Serial numberNameNumber of measuring points (number)Hs (m)Hr (m)MRCT-SRS (m)Distribution rangeMean valueDistribution rangeMean value1Zhanshan Road Station4113.0 ∼ 20.516.80.00.01.80–2.502Jiangxi Road Station689.1 ∼ 1411.50.00.01.12–2.103Zhongshan Park618.5 ∼ 16.511.60.0 ∼ 1.60.81.10–2.254Junfeng Road Station190.7 ∼ 6.83.24.7 ∼ 9.87.30.03–0.805Haichuan Road Station280.8 ∼ 8.34.37.0 ∼ 119.00.03–0.906Zhiquan Road Station196.1 ∼ 16.58.98.2 ∼ 14.311.20.50–0.90


## 4.2.1. Data Resource

The monitoring data come from six station tunnels of the two metro lines, which have approximately the same section size, excavation depth and technique, support measures, etc. The number of monitoring points, the soft soil stratum thicknessHs, the rock cover thickness Hr, and the MRCT-SRS are shown in Table 8. The statistical data are from a third party. The monitoring points for the ground surface settlement are selected within the main area of the tunnel, among which the abnormal points should be excluded in the statistics.Table 8
Basic overview of metro station tunnel.
Serial numberNameNumber of measuring points (number)Hs (m)Hr (m)MRCT-SRS (m)Distribution rangeMean valueDistribution rangeMean value1Zhanshan Road Station4113.0 ∼ 20.516.80.00.01.80–2.502Jiangxi Road Station689.1 ∼ 1411.50.00.01.12–2.103Zhongshan Park618.5 ∼ 16.511.60.0 ∼ 1.60.81.10–2.254Junfeng Road Station190.7 ∼ 6.83.24.7 ∼ 9.87.30.03–0.805Haichuan Road Station280.8 ∼ 8.34.37.0 ∼ 119.00.03–0.906Zhiquan Road Station196.1 ∼ 16.58.98.2 ∼ 14.311.20.50–0.90


## 4.3. Inversion Analysis

Table9 shows that, with regard to the overburden thicknesses exceeding the MRCT-SRS requirements for the stations of Zhiquan Road, Haichuan Road, and Junfeng Road, the average surface settlement is between 2.29 mm and 6.37 mm, with a maximum value between 6.22 mm and 12.59 mm, which are small in general. Regarding those for the stations of Zhongshan Park, Jiangxi Road, and Zhanshan Road, the average surface settlement is between 37.5 mm and 54.75 mm, with a maximum value between 71.54 mm and 80.96 mm, which are large in general.Table 9
Statistical table of ground surface settlement of the station tunnel.
NameScope (mm)Mean value (mm)Coefficient of variationStandard deviation (mm)The size relationship between Hr and MRCT-SRSZhiquan Road Station1.37 ∼ 6.222.290.361.18Hr > MRCT-SRSHaichuan Road Station0.18 ∼ 9.974.591.872.45Hr > MRCT-SRSJunfeng Road Station2.03 ∼ 12.596.372.093.05Hr > MRCT-SRSZhongshan Park Station17.31 ∼ 71.5437.580.3513.30Hr < MRCT-SRSJiangxi Road Station16.79 ∼ 76.9846.330.2712.69Hr < MRCT-SRSZhanshan Road Station21.19 ∼ 80.9654.750.3116.97Hr < MRCT-SRSSurface settlement can comprehensively reflect the effect of metro construction quality control and the degree of construction influence on the surrounding rock as well as the soil and neighboring environment. For a metro tunnel with shallow buried depth under the approximately same conditions of section size, excavation depth and technique, support measures, etc., the rock cover thickness on the vault meets the MRCT-SRS requirements or not, the surface settlement caused by the tunnel excavation is different by an order of magnitude. This situation fully demonstrates that the self-stability of the surrounding rock has a good correlation with the rock cover thickness (as shown in Figure12). The rationality of the conclusions obtained in this study has been fully verified.Figure 12
The change of ground settlement with different rock cover thicknesses.


## 4.4. Engineering Application

The fundamental goal of an excellent metro line design is to make full use of the self-stability of the surrounding rock and to minimize supplementary engineering measures. Costs should be considered only when engineering security is guaranteed. Furthermore, geological risks and construction difficulty should be reduced. The degree to which the surrounding rock can be self-stabilized depends mainly on the distribution conditions of the strata along the metro lines, the understanding of these conditions, and the ability to execute the required tasks. Therefore, the Geologically Easily Stable Area (GSEA) in this study refers to the area that meets the MRCT-SRS requirements for the excavation of the tunnel along the metro in the soil-rock dualistic stratum, otherwise referred to as the Geologically Difficultly Stable Area (GSDA).The metro station and the running tunnel are considered, respectively, as a single-arch large-span excavation structure with the excavation span of 20.8 m and a horseshoe-shaped structure with the excavation span of 6.2 m. The buried depth is set at equal intervals of 3 m falling within the scope of 6 m∼30 m. Table10 and Figure 13 show the calculation results of the proportions of the GSEA lengths of the station (considering only 16 subterranean excavation stations) and the running tunnels, respectively, as well as the lengths of the station and the running tunnel of Qingdao’s two metro lines under different buried depths.Table 10
The GSEA ratio of metro tunnel under different buried depths.
H (m)6912151821242730Station tunnel (%)49.265.076.484.190.794.398.099.0100.0Running tunel (%)26.340.557.071.687.593.996.698.599.3Figure 13
GSEA proportion variation for metro tunnel under different buried depths.As the buried depth increases, the proportion of the GSEA in the station tunnel rises first rapidly but then slowly. When the buried depth of the station tunnel is at 6 m, the proportion is 49.2%. When the depth increases to 12 m, the proportion increases rapidly to 76.4%. However, when the depth reaches 18 m, the proportion increase to 90.7%. The change of GSEA with buried depth in running tunnel is similar to that in station tunnel. When the buried depth of the running tunnel is 6 m, the proportion is only 26.3%. However, when the depth increases to 15 m, the proportion rapidly increases to 71.6%. As the depth continues to rise to 18 m, the proportion increases to 87.5%. When the depth reaches 21 m, the proportion slowly increases to 93.9%.The metro line is mainly composed of the station and the running tunnel. The buried depth of the line depends on the depth of the station tunnel. The depth of the running tunnel is mainly determined by factors such as the buried depth of the station tunnel and the longitudinal slope. The height difference between the station and the running tunnels is set as 12 m. Therefore, when the buried depth of the Qingdao metro station falls within 12m∼15 m, the GSEA of the station tunnel accounts for 76.4%∼84.1%, and the GSEA of the running tunnel accounts for 96.6%∼98.5%.


## 5. Conclusion

(1)
The MRCT-SRS of the metro tunnel approximately linearly increases or decreases with respect to the soft stratum thicknessHs in the soil-rock dualistic stratum. When the Hs is larger, so is the corresponding MRCT-SRS.(2)
The 3D spatial distribution was determined, and the fitted curves were constructed according to the MRCT-SRS,Hs, and D. This facilitated the assessment of the stability of the surrounding rock in the excavation of a metro tunnel.(3)
The variation rule of the Geologically Easily Stable Area with respect to buried depth has been revealed by examining the distribution features of the strata along the Qingdao metro lines. Moreover, the application of the MRCT-SRS in engineering has been expanded to determine the reasonable buried depth of the metro line.



---

*Source: 9120983-2020-07-08.xml*