The influence of boundary soil properties on tunneling-induced ground settlement is generally not considered in current analytic solutions, and the hypothesis of equal initial stress in vertical and horizontal makes the application of the above solutions limited. Based on the homogeneous half-plane hypothesis, by defining the boundary condition according to the ground loss pattern in shallow tunnel, and with the use of Mohr-Coulomb plastic yielding criteria and classic Lame and Kiersch elastic equations by separating the nonuniform stress field to uniform and single-direction stress field, a semiempirical solution for ground settlement induced by single shallow circular tunnel is presented and sensitivity to the ground parameters is analyzed. The methods of settlement control are offered by influence factors analysis of semiempirical solution. A case study in Beijing Metro tunnel shows that the semiempirical solution agrees well with the in situ measured results.

1. Introduction

With the rapid development of China’s urban infrastructure construction, the shield tunneling is widely used as a safe and efficient tunnel construction method in urban subway, underground pipeline, and so on. At present, there are still a series of problems needing to be solved in the metro shield tunneling construction; one of the problems that should be urgently solved is the prediction and control of the shield tunneling construction induced ground surface settlement [1, 2]. The research methods of this problem include the observational method, numerical simulations, and analytical methods [3].

In recent years, a large number of analytical solutions for axisymmetric circular tunnels in strain plane state have been presented [311], but those analytic methods are seldom considering the influence of ground soil properties on shield tunnel induced ground settlement, and some methods are based on the assumption of homogeneous semi-infinite plane; the calculated results are suitable for the deep buried tunnel surface settlement calculation but are not fully applicable for most shallow-buried tunnel, which limits the use of those methods. Therefore, it is an urgent problem to find suitable methods to calculate shallow shield tunnel surface settlement.

2. Semianalytical Analysis of Shallow Single Tunnel

As shown in Figure 1, the assumption of shallow tunnel surrounding rock analytic analysis is as follows: the rock is homogeneous and isotropic; tunnel excavation is a plane strain axisymmetry problem; the original stress affecting the surrounding rock is not equal in the vertical and the horizontal direction; the static earth pressure coefficient is 0.5. The other parameters are tunnel buried depth of (m), average overburden soil weight (kN/m3), excavation radius (m), plastic radius (m), and face support reaction stress (kPa).

Equilibrium differential equation for axisymmetric problem is

During the tunnel excavation, the soil deformation is elastic; then the expression of the stress-strain relationship of axisymmetric elasticity is where is the Poisson ratio of surrounding rock of tunnel and is Young’s modulus of surrounding rock of tunnel.

Two-way nonisobaric state could be divided into the external pressure as uniform stress field and the external pressure as one-way stress field. For uniform stress field and one-way stress field, the Lame and Kiersch formula is, respectively, used to calculate the stress field; the elastic zone can be obtained by summing the above calculated uniform stress field and one-way stress field of the surrounding rock:

By (2) and (3), and considering the boundary conditions and symmetry, the elastic zone displacement can be obtained, and the tunnel excavation induced relative displacement can be obtained by subtracting the initial stress generated initial displacement from the calculated elastic zone displacement , as shown inwhere is shear modulus of soil.

The Mohr-Coulomb theory strength criterion is adopted to judge whether the soil is in plastic state. The yield condition iswhere is soil cohesion, kPa, which can be obtained from the shear test and is soil internal friction angle, °, which can be obtained from the shear test.

The stress expression of the plastic zone can be obtained by (1) and (5), shown aswhere is soil cohesion, kPa; is soil internal friction angle, °; is tunnel excavating radius, m; is face supporting pressure, kPa; and is the radius of the calculating stress point, m.

Along the boundary between the elastic and plastic zone, the stress boundary conditions should be met at the same time; that is, there is when the radius ; the plastic radius is shown as follows:

When , the Mohr-Coulomb yield criterion nonassociative plastic flow method is used to solve the displacement [12]. The tunnel excavation is unloading process; (8) can be obtained by nonassociated flow rule.where is dilatancy angle of soil.

Substituting (2) and (3) into (8), (9) can be available.where .

For the soil or rock, the value of is very large, the value to the right of (9) is very small, and then (9) can be simplified to (10) by assuming no deformation in the tunnel axis direction and neglecting the elastic strain in the plastic zone [6].

From the geometric equation , , and (10), the relative displacement of the plastic zone can be solved as follows:

The deformation of (11) is radial uniform, which is different from the actual situation. Park [13] presented that the ground deformation model for shallow-buried tunnel displacement boundary conditions is adopted to modify (11) so that it can better predict the subsidence of shallow-buried tunnels. The displacement boundary expression of Park is shown as (12), where is the uniform radial deformation boundary, which is the value of (11) at the position; the meaning of is shown in Figure 1, , , and is the vertical distance of the forecast point to the ground surface. The deformation pattern of (12) is shown in Figure 2.

3. Parameters Sensitivity Analysis of Semianalytical Solution

For a given depth and radius of the tunnel, taking the soil horizontal pressure coefficient , it is shown that the influence parameters of the ground deformation are , Young’s modulus , Poisson’s ratio , soil cohesion , internal friction angle , and dilatancy . The calculated results of (13) by Matlab software are shown as Figures 37. Figure 3 shows the influence of Young’s modulus on the settlement when  kPa, , , and ; Young’s modulus is inversely proportional to the maximum subsidence of the surface. Figure 4 shows the influence of Poisson’s ratio on the ground settlement, which shows that the maximum settlement increases with increase of Poisson’s ratio .

Figure 5 is the influence of soil internal friction angle on surface settlement when  MPa, ,  kPa, and , Figure 6 is the influence of cohesion on surface settlement when  MPa, , , and , and Figure 7 is the influence of dilatancy on surface settlement when  MPa,  kPa, , and . Figures 57 show that the maximum surface settlement decreases with and increase and increases with , , and increase. In comparison, the influence of , , , and on the maximum surface settlement is larger than .

The above analysis for the influence of surrounding soil properties on the ground subsidence shows that there are two ways to control the shield tunnel induced ground settlement: one way is the active way, which is to reduce the surrounding soil properties weakening degree by optimizing the construction parameters of shield tunneling, such as choosing the rationale supporting face force and excavating speed or shield tail grouting reinforcement; the other way is the passive protection measures, which is to improve soil physical parameters by engineering methods to reduce the surface settlement of shield construction, the feasible methods including bolt reinforcement, partition wall, pile foundation underpinning, or the affected building body reinforcement [14].

4. Case Study

The 11th Project of Beijing Metro Line 10 is underground twin tunnels and is constructed by shield tunnel; the 6 m diameter tunnels are buried 16 m deep with nonsynchronous construction; the left tunnel is excavated more than 200 m ahead of the right tunnel. The ground elevation is 38.15 m. The shield tunnels are mainly through ⑤ the medium-coarse sand, ⑥ silty clay, ⑥1 clay, and ⑥2 silt. The perched groundwater is 2.74–3.58 m deep and the diving groundwater is 7.26–7.52 m deep. The layered and main mechanical parameters of foundation soil are found in Table 1.

As shown in Figure 8, the left shield tunnel machine is nearby underpassing the North Street 8# and South Street 8# residential building; in this underpassing process, the shield machine is just in the transition from circular curve to straight line, and the two tunnels center distance is 11.15 m~15.28 m, which puts great challenge to the control of strata settlement caused by the shield construction. In order to carry out the information construction and to reasonably control the ground surface subsidence for ensuring the safety of the affected building group, the corresponding settlement measurement point is laid out before the construction.

According to the relevant provisions of the national code of “code for design of building foundation (GB50007-2011),” taking into account the importance, the longtime service, and the poor deformation coordination of the upper panel type structure of the tunnels excavation affected buildings, the foundation deformation limit should be more strict than the national code; the previous studies suggest that the whole building tilt value should be controlled within 0.002, but the building tilt value of North Street 8# and South Street 8# residential building should be controlled within 0.001 and the maximum settlement under the wall which is nearby the street should be less than 10 mm because those buildings have tiled about 0.001 [15]. In order to test the impact of shield construction on the settlement of buildings, according to the mechanical properties of various strata and construction parameters, and considering the shield tunneling straightening on the surrounding soil disturbance, the following parameters are used for analytical analysis of the shield tunnel induced ground surface settlement, that is, the tunnel buried depth  m, the static earth pressure coefficient , excavation radius  m, weight of soil  kPa, face support reaction stress  kPa, cohesion  kPa, internal friction angle , Young’s modulus  MPa, Poisson ratio , and dilatancy angle . After the shield tunnel straightened, the shield tunneling causing surrounding soil properties weakening is relatively small, so the above-mentioned parameters were adjusted to  kPa, ,  MPa, , and . The two sets of parameters are taken into (13), and the Matlab7.0 is used for programming; the calculated results are shown in Figure 9. The maximum ground surface settlement is 23.8 mm when shield tunnel is in straightening phase, and the maximum ground surface settlement is 16.4 mm after posture adjustment.

Figure 9 shows that the maximum settlement of the foundation under the external wall nearby street will be more than 3 mm due to the left-line tunnel construction if the corresponding engineering measures are not taken. With the same construction technology and considering the influence of left tunnel excavation, the right tunnel construction induced maximum settlement of the foundation under the external wall nearby street will exceed 10 mm; it is difficult to ensure the safety of nearby buildings. So, the double-row high-pressure jet grouting cement piles were constructed between the external wall nearby street and the right tunnel; before the shield tunneling, the North Street 8# building pile diameter is 500 mm with 1000 mm center distance, and the South Street 8# building pile diameter is 450 mm with 950 mm center interval. Pile length is more than 1 m at the bottom of the shield tunnel, the water-cement ratio of pile is 0.8 ~1.2, the P.O.32.5 Portland cement is used, and the I20~I30 steel I-beam is inserted in piles. In the construction process, the inertia synchronous grouting slurry is replaced as hard synchronous grouting slurry, and the amount of cement is increased, the secondary grouting pressure is 0.3~0.4 MPa, and the grouting amount is 1.5 m3 per liner segment.

Figure 10 is settlement monitoring results of some sections; it shows that the left tunnel excavation induced about 5 mm ground surface subsidence on the right tunnel axis which is 11 m away from left tunnel axis and induced about 2 mm ground surface subsidence on the building which is 18 m away from left tunnel axis; the measured results of shield tunneling causing ground surface settlement are accordant with the analytical calculation results. It is revealed that the foundation reinforcement and grouting slurry adjustment measures are reasonable. By adopting active and passive engineering measures, the twin tunnels successfully underpassed the adjacent buildings, the shield tunneling-induced ground settlement is in standard control, and the nearby buildings safety is ensured.

5. Conclusions

Based on the assumption of homogeneous semi-infinite space, the initial two-way nonisobaric stress field of shallow-buried shield tunnel is divided into the uniform stress field and the one-way stress field; the semiempirical solution for shallow tunnel surface subsidence prediction is presented by using the analytical mechanic theory and definition the surrounding rock displacement boundary conditions of shield tunnel. The Matlab7.0 is chosen for programming to analyze the influence of the surrounding rock and soil mass properties on the ground surface settlement, and the corresponding surface settlement control method is proposed. The results of this research are used in the 11th Project of Beijing Metro Line 10 to predict and control the settlement of adjacent residential buildings affected by shield tunneling; the results show that the semiempirical solution is in good agreement with the in situ measured results; it can be used in the control of shield tunnel induced ground subsidence.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


The support of the Fujian Province Natural Science Fund (Grant no. 2016J01743) and the support of the Sanming University Science Research Development Fund and Fujian Province Young Teacher Education Research Project (Grant nos. JAT160454 and B201603) are gratefully acknowledged.