Research Article  Open Access
Jifeng Liu, Huizhi Zhang, "Semianalytical Solution and Parameters Sensitivity Analysis of Shallow Shield TunnelingInduced Ground Settlement", Advances in Materials Science and Engineering, vol. 2017, Article ID 3748658, 8 pages, 2017. https://doi.org/10.1155/2017/3748658
Semianalytical Solution and Parameters Sensitivity Analysis of Shallow Shield TunnelingInduced Ground Settlement
Abstract
The influence of boundary soil properties on tunnelinginduced 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 halfplane hypothesis, by defining the boundary condition according to the ground loss pattern in shallow tunnel, and with the use of MohrCoulomb plastic yielding criteria and classic Lame and Kiersch elastic equations by separating the nonuniform stress field to uniform and singledirection 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 [3–11], 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 semiinfinite plane; the calculated results are suitable for the deep buried tunnel surface settlement calculation but are not fully applicable for most shallowburied 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/m^{3}), 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 stressstrain 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.
Twoway nonisobaric state could be divided into the external pressure as uniform stress field and the external pressure as oneway stress field. For uniform stress field and oneway 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 oneway 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 MohrCoulomb 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 MohrCoulomb 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 shallowburied tunnel displacement boundary conditions is adopted to modify (11) so that it can better predict the subsidence of shallowburied 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 3–7. 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 5–7 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 mediumcoarse 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 (GB500072011),” 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 abovementioned 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 leftline 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 doublerow highpressure 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 watercement ratio of pile is 0.8 ~1.2, the P.O.32.5 Portland cement is used, and the I20~I30 steel Ibeam 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 m^{3} 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 tunnelinginduced ground settlement is in standard control, and the nearby buildings safety is ensured.
5. Conclusions
Based on the assumption of homogeneous semiinfinite space, the initial twoway nonisobaric stress field of shallowburied shield tunnel is divided into the uniform stress field and the oneway 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.
Acknowledgments
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.
References
 L. G. Tao, B. Liu, C. G. Ding et al., “Numerical simulation of ground settlement due to constructing metrostation based on shield tunneling,” Journal of China University of Mining & Technology, vol. 32, no. 3, pp. 236–240, 2003. View at: Google Scholar
 M. Y. Fattah, K. T. Shlash, and N. M. Salim, “Prediction of settlement trough induced by tunneling in cohesive ground,” Acta Geotechnica, vol. 8, no. 2, pp. 167–179, 2013. View at: Google Scholar
 Y. F. Xu, J. S. Chen, and D. M. Fu, “Effect of shield tunneling on mechanical properties of soils,” Chinese Journal of Rock Mechanics and Engineering, vol. 22, no. 7, pp. 1174–1179, 2003. View at: Google Scholar
 A. Verruijt, “A complex variable solution for a deforming circular tunnel in an elastic halfplane,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 21, no. 2, pp. 77–89, 1997. View at: Google Scholar
 W.I. Chou and A. Bobet, “Predictions of ground deformations in shallow tunnels in clay,” Tunnelling and Underground Space Technology, vol. 17, no. 1, pp. 3–19, 2002. View at: Publisher Site  Google Scholar
 A. Bobet, “Analytical solutions for shallow tunnels in saturated ground,” Journal of Engineering Mechanics, vol. 127, no. 12, pp. 1258–1266, 2001. View at: Publisher Site  Google Scholar
 C. González and C. Sagaseta, “Patterns of soil deformations around tunnels: application to the extension of Madrid metro,” Computers & Geosciences, vol. 28, no. 67, pp. 445–468, 2001. View at: Publisher Site  Google Scholar
 P.F. Li, Q. Fang, and D.L. Zhang, “Analytical solutions of stresses and displacements for deep circular tunnels with liners in saturated ground,” Journal of Zhejiang University—Science A (An International Applied Physics & Engineering Journal), vol. 15, no. 6, pp. 395–404, 2014. View at: Google Scholar
 H. N. Wang, Y. Li, Q. Ni et al., “Analytical solutions for the construction of deeply buried circular tunnels with two liners in rheological rock,” Rock Mechanics & Rock Engineering, vol. 46, no. 6, pp. 1481–1498, 2013. View at: Google Scholar
 K. Zeng and J. J. Xu, “Unified semianalytical solution for elasticplastic stress of deep circular hydraulic tunnel with support yielding,” Journal of Central South University, vol. 20, no. 6, pp. 1742–1749, 2013. View at: Google Scholar
 H. N. Wang, S. Utili, and M. J. Jiang, “An analytical approach for the sequential excavation of axisymmetric lined tunnels in viscoelastic rock,” International Journal of Rock Mechanics and Mining Sciences, vol. 68, pp. 85–106, 2014. View at: Publisher Site  Google Scholar
 H. S. Yu and R. K. Rowe, “Plasticity solutions for soil behaviour around contracting cavites and tunnels,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 23, no. 12, pp. 1245–1279, 1999. View at: Google Scholar
 K.H. Park, “Analytical solution for tunnellinginduced ground movement in clays,” Tunnelling and Underground Space Technology, vol. 20, no. 3, pp. 249–261, 2005. View at: Publisher Site  Google Scholar
 Z. S. Wang and M. S. Wang, “Effects of shielddriven tunneling on the satety of adjacent buildings and its countermeasures,” China Safety Science Journal, vol. 4, no. 2, pp. 45–49, 2002. View at: Google Scholar
 C. Y. Heng, Y. J. Teng, and X. Q. Chen, “Research of allowable deformation value of the building foundation soils above metro shield tunnel,” Chinese Journal of Underground Space and Engineering, vol. 12, no. 8, pp. 1336–1340, 2006. View at: Google Scholar
Copyright
Copyright © 2017 Jifeng Liu and Huizhi Zhang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.