m (LUZhongxiang moved page Draft Zhongxiang 374164140 to Review 537806263752)
(References)
 
(55 intermediate revisions by 3 users not shown)
Line 4: Line 4:
 
==Analysis of deformation of levee caused by large diameter shield tunneling in river crossing tunnel==
 
==Analysis of deformation of levee caused by large diameter shield tunneling in river crossing tunnel==
  
LU Zhongxiang<sup>1</sup>,SHOU Lingchao<sup>2</sup>,CHEN Yintao<sup>1</sup>,SHAO Xin<sup>3</sup>,FAN Wenju<sup>1</sup>,Amber Islam<sup>1</sup>
+
LU Zhongxiang<sup>1</sup>,SHOU Lingchao<sup>2</sup>,CHEN Yintao<sup>1</sup>,SHAO Xin<sup>3</sup>,FAN Wenju<sup>1</sup>,Amber Islam<sup>1</sup>,ZHANG Hengxu<sup>4</sup>
  
 
1 Universiti Teknologi Malaysia
 
1 Universiti Teknologi Malaysia
Line 12: Line 12:
 
3 Zhejiang Academy of Transportation Science
 
3 Zhejiang Academy of Transportation Science
  
<span id='_GoBack'></span>
+
4 China Railway 14th Bureau Group Corporation Limited
 
-->
 
-->
 +
==Abstract==
  
'''Abstract''':This study utilizes ABAQUS finite element software to analyze embankment deformation during shield tunneling. Results reveal that shield tunneling effects are manageable.The left tunnel line successfully tunnels beneath the embankment, causing a maximum settlement of 8.7 mm, meeting regulations. The ground surface exhibits a "V"-shaped lateral settlement trough, approximately five times the tunnel centerline width. Conversely, the right tunnel line induces a more extensive impact, with a 12.2 mm maximum settlement.Different lateral settlement patterns emerge in soil at varying depths. The left tunnel creates a "V" shape, with slightly increasing settlement above the tunnel axis as depth increases, accompanied by narrower troughs. The right tunnel line results in a pronounced "W" shape, especially at the grouting layer depth.Horizontal displacement forms an "S" shape, with maximum displacement at positions ±i from the tunnel axis. The right tunnel line induces greater maximum horizontal displacement than the left.In summary, this numerical analysis provides insights into embankment deformation during shield tunneling, aiding in assessing the impact and guiding safe tunneling beneath embankments.
+
This study utilizes ABAQUS finite element software to analyze embankment deformation during shield tunneling. Results reveal that shield tunneling effects are manageable. The left tunnel line successfully tunnels beneath the embankment, causing a maximum settlement of 8.7 mm, meeting regulations. The ground surface exhibits a “V” shaped lateral settlement trough, approximately five times the tunnel centerline width. Conversely, the right tunnel line induces a more extensive impact, with a 12.2 mm maximum settlement. Different lateral settlement patterns emerge in soil at varying depths. The left tunnel creates a “V” shape, with slightly increasing settlement above the tunnel axis as depth increases, accompanied by narrower troughs. The right tunnel line results in a pronounced “W” shape, especially at the grouting layer depth. Horizontal displacement forms an “S” shape, with maximum displacement at positions ±i from the tunnel axis. The right tunnel line induces greater maximum horizontal displacement than the left. In summary, this numerical analysis provides insights into embankment deformation during shield tunneling, aiding in assessing the impact and guiding safe tunneling beneath embankments.
  
'''Keywords: '''shield tunnel;finite element simulation;Levee settlement;data fitting
+
'''Keywords:''' Shield tunnelfinite element simulation,  levee settlementdata fitting
  
<br/>
+
==1.  Introduction==
  
With the rapid development of urban transportation in China, the number and scale of urban cross-river tunnel constructions are also expanding. These tunnels, primarily built using large-diameter mud shield methods, have seen their outer diameters grow from approximately 11 meters <sup>[1]</sup> to about 15 meters <sup>[2]</sup>. Notably, these river-crossing tunnels must traverse through river embankments, inevitably causing soil disturbance, particularly during the construction of super-large diameter mud shield tunnels.Given the intricacy of embankment structures and the vital role that river embankments play in safeguarding the lives and property of communities on both sides of the river, it becomes crucial to adopt rational construction parameters and measures. This approach aims to minimize the impact of shield tunneling on embankment structures.
+
With the rapid development of urban transportation, the number and scale of urban cross-river tunnel constructions are also expanding. These tunnels, primarily built using large-diameter mud shield methods, have seen their outer diameters grow from approximately 11 metersto about 15 meters [1,2]. Notably, these river-crossing tunnels must traverse through river embankments, inevitably causing soil disturbance, particularly during the construction of super-large diameter mud shield tunnels. Given the intricacy of embankment structures and the vital role that river embankments play in safeguarding the lives and property of communities on both sides of the river, it becomes crucial to adopt rational construction parameters and measures. This approach aims to minimize the impact of shield tunneling on embankment structures.
  
In recent years, the engineering field has accumulated some experience on the disturbance effect of shield tunneling through levees. Zhang Yazhou et al.<sup> [3] </sup>used finite-difference software to establish a three-dimensional numerical model of Nanjing Weisan Road super-diameter shield tunnel crossing the Yangtze River embankment, and calculated the settlement law and overall stability of the embankment. With the Chengjiang West Road tunnel in Jiangyin as the background, Han Lei et al. <sup>[4]</sup> used the displacement control finite element method based on the formation loss rate to calculate the embankment settlement under different formation loss rates and put forward control suggestions. Based on monitoring data and Peck empirical formula, Li Zongliang et al. <sup>[5-6]</sup> analyzed the disturbance and settlement law of the Qiantang River levee during shield crossing construction of Qingchun Road Tunnel in Hangzhou City, and proposed settlement control measures when mud shield crossing the levee. Wu Weiyi et al. [7] applied FLAC3D finite difference program and combined with field measurement to calculate and analyze the deformation law of the levee caused by the mud-water balanced shield with a diameter of 11.68m penetrating the Qiantang River levee, and found that the transverse settling trough at the top of the levee was similar to that of the subway shield with a diameter of 6.34m, showing Gaussian normal distribution. Peck formula can still be used to estimate the settlement, but the variation of settlement amount and trough width with depth is not as obvious as that of 6.34m diameter subway shield. Huang Hai [8] studied the embankment surface settlement caused by the large-diameter shield tunneling under Qiantang River, and found that the transverse settlement trough of the levee basically accords with normal distribution, the settlement is the largest at the center of the tunnel axis, and the surface settlement away from the tunnel center gradually decreases, and the settlement is symmetrically distributed along the tunnel axis, which accords with the normal distribution curve proposed by Peck. Zhu Muyuan [9] et al., relying on the Jilou Road Tunnel project, analyzed the settlement rules and characteristics of the large embankment section during the process of shield tunnel through the levee through the field measured data, and adopted a fine dynamic simulation of the shield tunnel driving process, and found that the initial settlement rate of the large embankment section during the shield tail ejection stage was larger and more susceptible to disturbance than that of the ordinary flat section. And the top of the levee is more sensitive to the disturbance of deep soil. Xie Xiongyao et al. [10] used finite element software to establish a three-dimensional numerical model for the super-diameter shield tunneling levee across the Yangtze River, and carried out a fine simulation of the construction process on the basis of considering construction elements such as soil excavation, mud-water support, segment lining, grouting behind the wall and formation loss. This paper summarizes the technical measures of settlement control in slurry shield tunneling from 7 aspects, such as slurry pressure control and real-time measurement of synchronous grouting quantity. Lin, QB et al. <sup>[11] </sup>takes Liuyang River levee engineering section of Changsha Metro Line 6 as the engineering background, and studies the levee deformation under shield tunnel excavation condition. Wu Shiming et al.<sup> [12]</sup> analyzed the risk sources of mud shield tunneling through levee engineering, discussed the causes, hazards and avoidance measures of the risks, and verified the rationality of the risk control measures by combining the engineering example of mud shield tunneling through levee in Qingchun Road Tunnel in Hangzhou. Most of the existing research results are the analysis of the settlement of levees, but there are few studies on the deep soil settlement and horizontal displacement characteristics of levees. Therefore, the analysis of three-dimensional deformation characteristics of levees after shield tunneling is helpful to understand the law and control the construction risk.
+
In recent years, the engineering field has accumulated significant experience regarding the disruptive effects of shield tunneling through levees. Zhang et al. employed finite-difference software to construct a three-dimensional numerical model of the Nanjing Weisan Road super-diameter shield tunnel, which crosses the Yangtze River embankment. They calculated the settlement pattern and overall stability of the embankment [3]. Building upon the context of the Chengjiang West Road tunnel in Jiangyin, Han Lei et al. employed the displacement control finite element method, which is based on the formation loss rate, to compute the embankment settlement under various formation loss rates. They also proposed recommendations for control [4]. Utilizing monitoring data and the Peck empirical formula, Li and Huang examined the disruption and settlement patterns of the Qiantang River levee during the shield crossing construction of the Qingchun Road Tunnel in Hangzhou City [5]. They also suggested settlement control measures for the mud shield crossing the levee [5,6]. Wu et al. utilized the FLAC3D finite difference program in conjunction with field measurements to compute and analyze the deformation patterns of the levee resulting from the penetration of the Qiantang River levee by the mud-water balanced shield, which has a diameter of 11.68m [7]. Their study revealed that the transverse settling trough at the top of the levee resembled that of a subway shield with a diameter of 6.34m, displaying a Gaussian normal distribution. The Peck formula remains applicable for estimating settlement; however, the change in settlement volume and trough width with depth is not as pronounced as observed with the 6.34m diameter subway shield. Huang Hai investigated the surface settlement of the embankment caused by large-diameter shield tunneling under the Qiantang River. The study revealed that the transverse settlement trough of the levee generally follows a normal distribution pattern, with the greatest settlement occurring at the center of the tunnel axis. Settlement diminishes gradually as distance from the tunnel center increases, and settlement distribution along the tunnel axis is symmetrical, in line with the normal distribution curve proposed by Huang [8]. Relying on the Jilou Road Tunnel project, Zhu et al. analyzed the settlement patterns and characteristics of the large embankment section during the shield tunneling process through the levee using field-measured data. They employed a detailed dynamic simulation of the shield tunnel driving process and observed that the initial settlement rate of the large embankment section during the shield tail ejection stage was greater and more susceptible to disturbance compared to that of an ordinary flat section [9]. Moreover, the top of the levee is particularly sensitive to the disturbance caused by deep soil. Xie et al. utilized finite element software to construct a three-dimensional numerical model for the super-diameter shield tunneling across the Yangtze River levee [10]. They conducted a detailed simulation of the construction process, considering construction elements such as soil excavation, mud-water support, segment lining, grouting behind the wall, and formation loss. Lin et al. based their study on the Liuyang River levee engineering section of Changsha Metro Line 6, They investigated the deformation of the levee under shield tunnel excavation conditions [11]. Wu et al. analyzed the risk factors associated with mud shield tunneling through levee engineering [12]. They examined the causes, hazards, and avoidance measures of these risks and validated the effectiveness of the risk control measures by referencing the engineering example of mud shield tunneling through the levee in the Qingchun Road Tunnel in Hangzhou. Tamir et al. simplified the modeling of the tunnel boring machine shield while incorporating the hardened soil constitutive model [13]. This approach enabled them to capture crucial aspects of the shield tunneling process, such as joint segment lining, grouting pressure, time-varying grout hardening, and nonlinear soil behavior. Hyobum et al. conducted laboratory pressure blade shear tests to examine the rheological properties of sandy soils under foam conditions [14]. They investigated how varying water content, foam injection ratio (FIR), and confining pressure affected the rheological properties of soil under foam conditions. Milad et al. utilized the finite element limit analysis method to assess the limit support pressure in undrained clay [15]. Their evaluation considered horizontal pseudo-static seismic forces directed outward from the face. Zdenek et al. presented a new approach for incorporating the formation process into the analysis of pressure transfer [16,17]. This method was specifically designed for shallow slurry penetration depths smaller than the cutting depth of the tools. Hoonil  et al. reported on the significant ground subsidence that ensued after the structural deformation of tunnel segments during the excavation of a cross-passage between two main tunnels [18]. These main tunnels were excavated by an earth pressure balance shielded tunnel boring machine (EPB shield TBM) and were constructed in the deltaic region near the Nak-Dong River in Korea.
  
In this paper, the deformation (settlement and displacement) of the levee caused by the large diameter shield machine under the construction of the Jiangshan Road tunnel in Hangzhou is numerically calculated, and the three-dimensional deformation characteristics of the levee are studied.
+
Most of the existing research focuses on the analysis of levee settlement, with few studies addressing the settlement and horizontal displacement characteristics of deep soil in levees. Therefore, analyzing the three-dimensional deformation characteristics of levees after shield tunneling is beneficial for understanding the patterns and managing construction risks. In this paper, the deformation (settlement and displacement) of the levee caused by the large-diameter shield machine during the construction of the Jiangshan Road tunnel in Hangzhou is numerically calculated, and the three-dimensional deformation characteristics of the levee are studied.
  
===1. Project overview and geological conditions===
+
==2. Project overview and geological conditions==
  
====1.1 Project Overview====
+
===2.1 Project overview===
  
"Genshandong Road Cross-River Tunnel" is a significant infrastructure project located between Xiisha Bridge and Jiangdong Bridge, spanning a total length of 4462.26m. This tunnel plays a vital role in the region's transportation network, and its technical aspects are noteworthy for a comprehensive understanding.This shield tunnel features a single-layer lining with segment assembly. It boasts an impressive outer diameter of 14.5m, an inner diameter of 13.3m, a segment thickness of 0.6m, and a ring width of 2m. The tunnel comprises 10 "9+1" sections, and these sections are assembled with staggered joints, tightly secured by high-strength bolts. The construction of the shield tunnel begins at the eastern starting well of the Big River and proceeds westward. Upon crossing the mighty Qiantang River, it reaches its destination at the Xiasha receiving well. It's worth noting that when the shield tunneling machine pierces through the Qiantang River levee, it faces a significant challenge due to the levee's high deformation requirements caused by tidal influence. For a visual representation, please refer to Fig 1, which illustrates the current status of the levee.
+
The Genshandong Road Cross-River Tunnel is a significant infrastructure project situated between the Xisha Bridge and Jiangdong Bridge, spanning a total length of 4462.26m. This tunnel plays a vital role in the region's transportation network, and its technical aspects are noteworthy for a comprehensive understanding. This shield tunnel features a single-layer lining with segment assembly. It boasts an impressive outer diameter of 14.5m, an inner diameter of 13.3m, a segment thickness of 0.6m, and a ring width of 2m. The tunnel comprises 10 “9+1” sections, and these sections are assembled with staggered joints, tightly secured by high-strength bolts. The construction of the shield tunnel, a significant infrastructure project in the region, begins at the eastern starting well of the Big River and proceeds westward, reaching its destination at the Xisha receiving well. However, the tunnel faces a major challenge when it crosses the mighty Qiantang River due to the high deformation requirements of the levee caused by tidal influence. For a visual representation of the current status of the levee, as shown in [[#img-1|Figure 1]].
  
{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;"  
+
<div id='img-1'></div>
 +
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"  
 +
|-style="background:white;"
 +
|align="center" |
 +
{|style="margin: 0em auto 0.1em auto;width:auto;"
 +
|+
 
|-
 
|-
| style="border: 1pt solid black;text-align: center;vertical-align: top;"|''' [[Image:Draft_Zhongxiang_374164140-image1.png|222px]] '''
+
|style="text-align: center;padding:10px;"| [[Image:Review_537806263752-image1.png|222px]]
| style="border: 1pt solid black;text-align: center;vertical-align: top;"|''' [[Image:Draft_Zhongxiang_374164140-image2.png|210px]] '''
+
|style="text-align: center;padding:10px;"| [[Image:Review_537806263752-image2.png|210px]]
| style="border: 1pt solid black;text-align: center;vertical-align: top;"|''' [[Image:Draft_Zhongxiang_374164140-image3.png|192px]] '''
+
|style="text-align: center;padding:10px;"| [[Image:Review_537806263752-image3.png|192px]]
 +
|}
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 1'''. Present situation of levees on the east bank of Qiantang River
 
|}
 
|}
  
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
===2.2  Hydrogeological conditions and monitoring arrangement===
Fig 1. Present situation of levees on the east bank of Qiantang River</div>
+
  
====1.2 Hydrogeological conditions and monitoring arrangement====
+
The shield tunnel passes through various soil layers, including sandy silt, silt, silty clay, mucky clay, clay, Fully weathered rock, and Moderately weathered rock. Despite the levee not being strengthened before shield construction, strategic monitoring points were installed on the levee to ensure safety. These measures are crucial for maintaining project stability and safety. For detailed information about the tunnel's soil layers and monitoring setup, as shown in [[#img-2|Figure 2]].
  
The shield tunnel passes through various soil layers, including sandy silt, silty clay, and sandy silty clay. It's worth noting that the levee wasn't strengthened before shield construction. However, monitoring points were strategically placed on the levee to ensure safety. Let's delve into the specifics of the tunnel's soil layers and monitoring setup in Figure 2. Safety is a top priority in this project. Despite the absence of levee strengthening before shield construction, we've taken precautions by installing settlement monitoring points on the levee. These measures are crucial to ensure the stability and safety of the project. Let's now examine the tunnel's soil layers and the monitoring section in Fig 2.
+
<div id='img-2'></div>
 +
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"
 +
|-style="background:white;"
 +
|style="text-align: center;padding:10px;"| [[Image:Review_537806263752-image4.png|550px]]
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 2'''. Monitoring site layout profile
 +
|}
  
[[Image:Draft_Zhongxiang_374164140-image4-c.png|600px]]
+
==3. Measured analysis of surface settlement of levees==
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
===3.1  Lateral surface settlement of levee top===
Fig.2monitoring site layout profile</div>
+
  
===2. Measured analysis of surface settlement of levees===
+
In [[#img-3|Figure 3]], we observe the settlement patterns above and around the tunnel axis. The most significant settlement occurs directly above the tunnel axis. For the left-line shield tunnel, the transverse surface settlement generally forms a “V” shape distribution, with a maximum settlement of 8.7mm directly above the left line.
  
====2.1 Lateral surface settlement of levee top====
+
Conversely, after the right-line shield tunnel crosses the levee, the transverse surface settlement generally exhibits a “W” shape distribution. The maximum settlement reaches 12.2mm, which is still directly above the left line of the shield tunnel and greater than the surface settlement above the right line. Additionally, the settlement in the middle of the left and right lines is slightly less than the settlement value above the tunnel.
  
In Fig 3, we observe the settlement patterns above and around the tunnel axis. The most significant settlement occurs directly above the tunnel axis. For the left-line shield tunnel, the transverse surface settlement generally forms a 'V' shape distribution, with a maximum settlement of 8.7mm, directly above the left line. Conversely, after the right-line shield tunnel crosses the levee, the transverse surface settlement generally exhibits a 'W' shape distribution. Here, the maximum settlement reaches 12.2mm, which is still directly above the left line of the shield tunnel and greater than the surface settlement above the right line. Additionally, the surface settlement in the middle of the left and right lines is slightly less than the settlement value above the tunnel.
+
<div id='img-3'></div>
 +
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"
 +
|-style="background:white;"
 +
|style="text-align: center;padding:10px;"| [[Image:Review_537806263752-image5.png|414px]]
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 3'''. Transverse surface settlement curve of levee top after the left and right lines cross the levee
 +
|}
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
[[Image:Draft_Zhongxiang_374164140-image5.png|468px]] </div>
 
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
For the prediction of ground displacement caused by single-line tunnel excavation, the horizontal distribution estimation formula of ground settlement proposed by Peck is generally adopted in current engineering practice [19]:
Fig.3 Transverse surface settlement curve of levee top after the left and right lines cross the levee</div>
+
  
For the prediction of ground displacement caused by single-line tunnel excavation, the horizontal distribution estimation formula of ground settlement proposed by Peck is generally adopted in current engineering practice <sup>[13]</sup> :
+
<div style="text-align: right; direction: ltr; margin-left: 1em;"></div>
  
<div style="text-align: right; direction: ltr; margin-left: 1em;">
+
{| class="formulaSCP" style="width: 100%; text-align: center;"
  [[Image:Draft_Zhongxiang_374164140-image6.png|162px]] (2-1)</div>
+
|-
 +
|
 +
{| style="text-align: center; margin:auto;"   
 +
|-
 +
| <math>S\left( x\right) ={S}_{\max}\exp(-\frac{{x}^{2}}{2{i}^{2}})</math>
 +
|}
 +
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
 +
|}
  
<div style="text-align: right; direction: ltr; margin-left: 1em;">
+
{| class="formulaSCP" style="width: 100%; text-align: center;"
[[Image:Draft_Zhongxiang_374164140-image7.png|150px]] (2-2)</div>
+
|-
 +
|
 +
{| style="text-align: center; margin:auto;"  
 +
|-
 +
| <math>{S}_{\max}=\frac{{V}_{s}}{i\sqrt{2\pi }}=\frac{\pi {R}^{2}\eta }{i\sqrt{2\pi }}</math>
 +
|}
 +
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
 +
|}
  
[[Image:Draft_Zhongxiang_374164140-image8.png|30px]] : Ground subsidence caused by ground loss; [[Image:Draft_Zhongxiang_374164140-image9.png|12px]] : the distance from the tunnel center line; [[Image:Draft_Zhongxiang_374164140-image10.png|30px]] : the maximum settlement caused by the formation loss at the center line of the tunnel; [[Image:Draft_Zhongxiang_374164140-image11.png|6px]] : the width coefficient of surface settling trough; [[Image:Draft_Zhongxiang_374164140-image12.png|18px]] : Stratum loss per unit length of tunnel; [[Image:Draft_Zhongxiang_374164140-image13.png|12px]] : the formation loss rate, is the ratio of the surface subsidence trough area to the tunnel excavation area; [[Image:Draft_Zhongxiang_374164140-image14.png|12px]] : the radius of tunnel excavation.
+
where <math display="inline">S(x)</math> is the ground subsidence caused by ground loss; <math display="inline">x</math> the distance from the tunnel center line; <math display="inline">{S}_{\max}</math> the maximum settlement caused by the formation loss at the center line of the tunnel; <math display="inline">i</math> the width coefficient of surface settling trough; <math display="inline">{V}_{s}</math> the stratum loss per unit length of tunnel; <math display="inline">\eta</math> the formation loss rate, is the ratio of the surface subsidence trough area to the tunnel excavation area; and <math display="inline">R</math> is the radius of tunnel excavation.
  
When the Peck formula is used to predict the transverse surface settlement, the most important thing is to determine the formation loss rate V<sub>1</sub> and the width coefficient of sedimentation tank [[Image:Draft_Zhongxiang_374164140-image11.png|6px]] . O'Reilly &New (1982) <sup>[13]</sup> proposed the following values in accordance with the experience of London:
+
When the Peck formula is used to predict the transverse surface settlement, the most important thing is to determine the formation loss rate <math display="inline">{V}_{1}</math> and the width coefficient of sedimentation tank <math display="inline">i</math>. O'Reilly and New [20] proposed the following values in accordance with the experience of London:
  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 82: Line 112:
 
{| style="text-align: center; margin:auto;"  
 
{| style="text-align: center; margin:auto;"  
 
|-
 
|-
| [[Image:Draft_Zhongxiang_374164140-image15.png|48px]]
+
| <math>i=K{z}_{0}</math>
 
|}
 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2-3)
+
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
 
|}
 
|}
  
 +
where <math display="inline">K</math> is the sedimentation tank width parameter; and <math display="inline">{z}_{0}</math>the  buried depth for the tunnel axis.
  
[[Image:Draft_Zhongxiang_374164140-image16.png|18px]] :sedimentation tank width parameter; [[Image:Draft_Zhongxiang_374164140-image17.png|18px]] :buried depth for the tunnel axis.
+
In the paper, the Peck formula is utilized to fit the measured transverse surface settlement curve of the embankment top after the left line crosses the embankment. This analysis aims to compare the surface settlement of the embankment top and evaluate the applicability of the Peck formula in the project, as shown in [[#img-4|Figure 4]]. The findings from this analysis are crucial for understanding the behavior of the embankment top and its response to the tunneling process.
  
In the paper, Peck formula was used to fit the measured transverse surface settlement curve of the embankment top after the left line crossed the embankment, further compared the surface settlement of the embankment top, and evaluated the applicability of Peck formula in the project, as shown in Fig .4.
+
<div id='img-4'></div>
 +
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"
 +
|-style="background:white;"
 +
|style="text-align: center;padding:10px;"| [[File:Review 537806263752-image7.png|455x455px]]
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 4'''. Comparison between the measured value of land surface settlement and the Peck fitting value
 +
|}
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
[[Image:Draft_Zhongxiang_374164140-image18.png|600px]] </div>
 
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
Upon analyzing [[#img-4|Figure 4]], we observe that the variation in the surface settlement curve closely aligns with the Peck formula. However, it's noteworthy that the measured values in the field are slightly larger than those fitted by the Peck formula. The discrepancy in measured values may arise due to the greater structural weight of the levee compared to typical strata. This difference is reflected in the higher soil and water pressure, as well as settlement, beneath the levee excavation surface at the same burial depth compared to shield tunneling in other general strata. The large diameter of the shield used influences these factors. Despite the challenges posed by large-diameter shield tunnels, the Peck formula remains applicable, allowing engineers to predict the distribution of transverse settling troughs at the top of the levee. These predictions, which consider factors such as formation loss (Vs) and the width coefficient of the settling trough (i), are invaluable for assessing and managing levee settlement and deformation during construction.
Fig. 4 Comparison between the measured value of land surface settlement and the Peck fitting value</div>
+
  
Upon analyzing Fig.4, we can observe that the variation in the surface settlement curve closely aligns with the Peck formula. However, it's noteworthy that the measured values in the field are slightly larger than those fitted by the Peck formula. This discrepancy may arise due to the greater structural weight of the levee compared to typical strata. Compared with shield tunneling in other general strata, the soil and water pressure, as well as settlement, beneath the levee excavation surface are greater at the same burial depth. This is influenced by the large diameter of the shield used. Given that the Peck formula remains applicable even in situations involving large-diameter shield tunnels, engineers can use it to quickly predict the distribution of transverse settling troughs at the top of the levee. This prediction considers the formation loss (Vs) caused by shield construction and the width coefficient of the settling trough (i). Such predictions are invaluable for assessing whether levee settlement values and uneven settlement slopes exceed deformation control thresholds, allowing timely adjustments to construction parameters.
+
===3.2  Surface settlement of shield through various stages===
  
====2.2 Surface settlement of shield through various stages====
+
These stages provide a framework for understanding the settlement process during shield tunnel construction and are crucial for monitoring and managing potential impacts on surrounding structures.
  
The surface settlement caused by shield tunnel construction can be divided into five stages based on the relative position and formation mechanism of the shield tunnel and monitoring section. The following is a detailed description of each stage:
+
# The shield approaches: the shield excavation face has not reached the monitoring section, and the distance is about 2 times the diameter of the shield.
 +
# The shield arrives: the excavation face of the shield has not reached the monitoring section, and the distance is 0~2 times the diameter of the shield.
 +
# The shield goes through: when the shield excavation face passes through the monitoring section, the distance is 1 times the diameter of the shield.
 +
# The shield leaves: when the shield excavation face leaves the monitoring section about 1 times the diameter of the shield.
 +
# Final stable times: about 10 times the diameter of the shield when the shield excavation face leaves the monitoring section.
  
Stage 1: When the shield tunnel approaches, and the shield excavation face has not yet reached the monitoring section, the distance is approximately 2 times the tunnel diameter (D :the tunnel diameter).
 
  
Stage 2: Prior to the arrival of the shield, when the shield excavation face has not yet reached the monitoring section, and the distance is within 0~2 times the tunnel diameter (D).
+
As shown in [[#img-5|Figure 5]], this is the lateral surface settlement curve of the levee at each stage of left-line shield construction. By comparing these curves, it can be observed that as the shield advances, surface settlement gradually increases. The width of the transverse settlement trough in the same section remains relatively constant, and the settlement curve becomes steeper. When the shield reaches the monitoring fault front (approximately 1D~2D away from the section, D: the tunnel diameter), a noticeable transverse settlement trough with a width of about 15~20 meters is generated. The maximum land settlement occurs at the center of the tunnel, with settlement decreasing as one moves farther away from the tunnel center.
  
Stage 3: When the shield passes through, and the shield excavation face crosses the monitoring section, with the tunnel diameter being D.
+
<div id='img-5'></div>
 +
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"
 +
|-style="background:white;"
 +
|style="text-align: center;padding:10px;"| [[File:Review 537806263752-image8.png|426x426px]]
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 5'''. Surface settlement curve of the left line shield crossing the levee at each stage
 +
|}
  
Stage 4: When the shield exits, and the shield excavation face moves away from around the monitoring section, with the tunnel diameter as D.
 
  
Stage 5: The final stable stage, where surface settlement for this section is measured for the last time, at approximately 10 times the tunnel diameter (D) when the shield excavation face leaves the monitoring section.
+
[[#tab-1|Table 1]] provides the percentage distribution of settlement in different stages of shield tunnel excavation. In the early settlement stage of the monitoring point, the settlement value accounts for 7% of the total settlement, which, while relatively small, should not be overlooked. Moving on to the arrival stage of the shield, we see that it contributes 14% of the total settlement. However, the phases of shield passage and post-shield passage are particularly critical, with these two stages contributing to a larger proportion of the total settlement, specifically 43% and 23%, respectively. In fact, settlements in these stages collectively account for approximately 66% of the total settlement value, making this phase the most concerning due to its rapid and substantial settlement changes. Finally, the late sedimentation stage produces around 13% of the total settlement value. To summarize, the percentage of settlement in various stages of shield tunnel excavation varies significantly. The early and late stages contribute relatively small percentages, while the arrival stage accounts for 14% of the total settlement. However, it's during the phases of shield passage and post-shield passage that settlement percentages are most substantial, combining to make up approximately 66% of the total settlement value. This phase is of particular concern due to its rapid and pronounced settlement changes.
  
"As shown in Fig 5, this is the lateral surface settlement curve of the levee at each stage of left-line shield construction. By comparing these curves, it can be observed that as the shield advances, surface settlement gradually increases. The width of the transverse settlement trough in the same section remains relatively constant, and the settlement curve becomes steeper. When the shield reaches the monitoring fault front (approximately 1D~2D away from the section), a noticeable transverse settlement trough with a width of about 15~20 meters is generated. The maximum land settlement occurs at the center of the tunnel, with settlement decreasing as one moves farther away from the tunnel center."
+
<div class="center" style="font-size: 75%;">'''Table 1'''.  Percentage of the total settlement at each stage of the monitoring point (%)</div>
  
Table 1 provides valuable insights into the percentage distribution of settlement in different stages of shield tunnel excavation. Understanding these percentages is crucial for assessing settlement patterns during the construction process. In the early settlement stage of the monitoring point, the settlement value accounts for 7% of the total settlement, which, while relatively small, should not be overlooked. Moving on to the arrival stage of the shield, we see that it contributes 14% of the total settlement. However, the phases of shield passage and post-shield passage are particularly critical, with these two stages contributing to a larger proportion of the total settlement, specifically 43% and 23%, respectively. In fact, settlements in these stages collectively account for approximately 66% of the total settlement value, making this phase the most concerning due to its rapid and substantial settlement changes. Finally, the late sedimentation stage produces around 13% of the total settlement value. To summarize, the percentage of settlement in various stages of shield tunnel excavation varies significantly. The early and late stages contribute relatively small percentages, while the arrival stage accounts for 14% of the total settlement. However, it's during the phases of shield passage and post-shield passage that settlement percentages are most substantial, combining to make up approximately 66% of the total settlement value. This phase is of particular concern due to its rapid and pronounced settlement changes.
+
<div id='tab-1'></div>
 
+
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
|-style="text-align:center"
[[Image:Draft_Zhongxiang_374164140-image19.png|480px]] </div>
+
!Pre-arrival settlement of shield !! Phase settlement of shield arrival !! Phase settlement of the shield as it passes through !! Phase settlement after shield passage !! Long-term subsequent stage settlement
 
+
Fig. 5 Surface settlement curve of the left line shield crossing the levee at each stage
+
 
+
Table 1  Percentage of the total settlement at each stage of the monitoring point (%)
+
 
+
{| style="width: 100%;border-collapse: collapse;"  
+
|-
+
style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Pre-arrival settlement of shield
+
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Phase settlement of shield arrival
+
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Phase settlement of the shield as it passes through
+
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Phase settlement after shield passage降
+
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Long-term subsequent stage settlement沉降
+
 
|-
 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|7
+
|  style="text-align: center;"|7
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|14
+
|  style="text-align: center;"|14
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|43
+
|  style="text-align: center;"|43
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|23
+
|  style="text-align: center;"|23
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|13
+
|  style="text-align: center;"|13
 
|}
 
|}
  
 +
==4.  Numerical analysis of levee deformation==
  
===3.Numerical analysis of levee deformation===
+
To gain insights into the soil deformation characteristics of the levee during sequential passage by the left and right-line shields, we conducted an investigation using the advanced finite element software ABAQUS. This analysis aims to calculate and compare the dynamic construction processes of the shields to understand how they impact the soil and the overall construction of the tunnel.
  
To investigate the soil deformation characteristics of the levee during sequential passage by the left and right-line shields, we have utilized the advanced finite element software ABAQUS to calculate and compare the dynamic construction processes.
+
===4.1  Numerical calculation model===
  
====3.1 Numerical calculation model====
+
To gain a comprehensive understanding of the real-world engineering scenario and design parameters, we have developed a numerical calculation model using ABAQUS finite element software. This model is designed to simulate the effects of the shield tunnel's construction, which comprises two parallel lines with a centerline separation of 35m. One of the key features of the model is its representation of the grouting filler within the shield tunnel gap as a uniform and equally thick grouting layer.
  
To better understand the real-world engineering scenario and design parameters, we have created a numerical calculation model using ABAQUS finite element software. This model aims to simulate the effects of the shield tunnel's construction. The tunnel consists of two parallel lines, with a centerline separation of 35m. We've represented the grouting filler within the shield tunnel gap as a uniform and equally thick grouting layer.
+
In our simulation, we have incorporated the thrust generated by the shield machine and the pressure exerted by cement in the form of nodal concentrated forces. Our model encompasses various elements, including the surrounding rock, pipe segments, grouting material, and the levees, all represented using solid elements.
  
In our simulation, we've applied the thrust generated by the shield machine and the pressure exerted by cement in the form of nodal concentrated forces. Our model encompasses various elements, including the surrounding rock, pipe segments, grouting material, and the levees, all represented using solid elements.
+
The numerical model’s dimensions are as follows: a longitudinal length of 56m, a transverse width of 180m, and a vertical height ranging from 76.1m to 84.6m. This extensive model is divided into 140,208 units and comprises 148,071 nodes, ensuring a comprehensive representation of the system. For a visual representation, as shown in [[#img-6|Figure 6]].
  
The numerical model's dimensions are as follows: a longitudinal length of 56m, a transverse width of 180m, and a vertical height ranging from 76.1m to 84.6m. This extensive model is divided into 140,208 units and comprises 148,071 nodes, ensuring a comprehensive representation of the system. For a visual representation, please refer to Fig 6.
+
<div id='img-6'></div>
 +
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"
 +
|-style="background:white;"
 +
|style="text-align: center;padding:10px;"| [[File:Review 537806263752-image9.png|426x426px]]
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 6'''. Numerical calculation model
 +
|}
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
[[Image:Draft_Zhongxiang_374164140-image20.png|336px]] </div>
 
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
In shield construction, a ring gap will be formed between the tunnel excavation wall and the outer ring of the segment when the segment is detached from the shield tail. If the action of grouting and pressure is ignored, the soil mass will move to this gap and produce a large displacement.In the actual grouting process, the slurry fills the gap of the shield tail and penetrates into the soil body, and then gradually hardens. If the hardening of the grouting material is not considered, the settlement value will also deviate from the actual value by using the isogeneration layer simulation. Because the interaction of shield tail grouting will form a mixture of soil and cement slurry, the grouting layer can be simulated by elastic model, and the hardening of grouting material can be simulated by the variable stiffness of grouting material [21]. And an elastic model is utilized for soil. The specific material parameters for these models are detailed in [[#tab-2|Table 2]].
Fig. 6 Numerical calculation model</div>
+
  
In our analysis, we employed the Mohr-Coulomb elastoplastic model for simulating soil behavior, while an elastic model was utilized for both the grouting layer and segment materials. The specific material parameters for these models are detailed in Table 2.
+
<div class="center" style="font-size: 75%;">'''Table 2'''. Material parameters of the model</div>
  
The boundary conditions and deformation constraints play a critical role in accurately representing the real-world scenario. At the top and surface of the levee, we have implemented free boundaries to allow for unrestricted movement. Around the levee, deformation constraints have been applied to capture the complex interactions. Additionally, at the bottom surface, we've limited deformation in the X, Y, and Z directions to reflect the ground's response. Furthermore, lateral constraints have been introduced on the sides of the model.
+
<div id='tab-1'></div>
 
+
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
|-style="text-align:center"
Table 2 Material parameters of the model</div>
+
! style="text-align: left;"| Material name !! Weight (kN/m³) !! Cohesion (kPa) !! Internal friction Angle (°) !! Poisson's ratio  !! Elastic modulus  (MPa)
 
+
{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;"  
+
 
|-
 
|-
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|Material name
+
|  style="text-align: left;"|Sandy silt
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|weight(kN/m³)
+
|  style="text-align: center;"|19.8
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|cohesion(c)kPa
+
|  style="text-align: center;"|12
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|internal friction Angle φ(°)
+
|  style="text-align: center;"|32
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|Poisson's ratio υ
+
|  style="text-align: center;"|0.2
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|elastic modulus
+
|  style="text-align: center;"|19
 
+
E(MPa)
+
|-
+
|  style="border-top: 2pt solid black;text-align: center;"|Sandy silt
+
|  style="border-top: 2pt solid black;text-align: center;"|19.8
+
|  style="border-top: 2pt solid black;text-align: center;"|12
+
|  style="border-top: 2pt solid black;text-align: center;"|32
+
|  style="border-top: 2pt solid black;text-align: center;"|0.2
+
|  style="border-top: 2pt solid black;text-align: center;"|19
+
 
|-
 
|-
|  style="text-align: center;"|silt
+
|  style="text-align: left;"|Silt
 
|  style="text-align: center;"|19.9
 
|  style="text-align: center;"|19.9
 
|  style="text-align: center;"|13.5
 
|  style="text-align: center;"|13.5
Line 191: Line 214:
 
|  style="text-align: center;"|20
 
|  style="text-align: center;"|20
 
|-
 
|-
|  style="text-align: center;"|Silty clay
+
|  style="text-align: left;"|Silty clay
 
|  style="text-align: center;"|17.6
 
|  style="text-align: center;"|17.6
 
|  style="text-align: center;"|19
 
|  style="text-align: center;"|19
Line 198: Line 221:
 
|  style="text-align: center;"|8
 
|  style="text-align: center;"|8
 
|-
 
|-
|  style="text-align: center;"|mucky clay
+
|  style="text-align: left;"|Mucky clay
 
|  style="text-align: center;"|18.4
 
|  style="text-align: center;"|18.4
 
|  style="text-align: center;"|27
 
|  style="text-align: center;"|27
Line 205: Line 228:
 
|  style="text-align: center;"|4
 
|  style="text-align: center;"|4
 
|-
 
|-
|  style="text-align: center;"|clay
+
|  style="text-align: left;"|Clay
 
|  style="text-align: center;"|19.6
 
|  style="text-align: center;"|19.6
 
|  style="text-align: center;"|45
 
|  style="text-align: center;"|45
Line 212: Line 235:
 
|  style="text-align: center;"|9
 
|  style="text-align: center;"|9
 
|-
 
|-
|  style="text-align: center;"|Silty clay
+
|  style="text-align: left;"|Fully weathered rock
|  style="text-align: center;"|19.8
+
|  style="text-align: center;"|21
+
|  style="text-align: center;"|33
+
|  style="text-align: center;"|0.26
+
|  style="text-align: center;"|7.5
+
|-
+
|  style="text-align: center;"|Fully weathered rock
+
 
|  style="text-align: center;"|22.3
 
|  style="text-align: center;"|22.3
 
|  style="text-align: center;"|50
 
|  style="text-align: center;"|50
Line 226: Line 242:
 
|  style="text-align: center;"|320
 
|  style="text-align: center;"|320
 
|-
 
|-
|  style="text-align: center;"|Moderately weathered rock
+
|  style="text-align: left;"|Moderately weathered rock
 
|  style="text-align: center;"|22.6
 
|  style="text-align: center;"|22.6
 
|  style="text-align: center;"|200
 
|  style="text-align: center;"|200
Line 233: Line 249:
 
|  style="text-align: center;"|550
 
|  style="text-align: center;"|550
 
|-
 
|-
|  style="text-align: center;"|Grouting layer
+
|  style="text-align: left;"|Grouting layer
 
|  style="text-align: center;"|21
 
|  style="text-align: center;"|21
|  style="text-align: center;"|/
+
|  style="text-align: center;"|-
|  style="text-align: center;"|/
+
|  style="text-align: center;"|-
 
|  style="text-align: center;"|0.2
 
|  style="text-align: center;"|0.2
 
|  style="text-align: center;"|15
 
|  style="text-align: center;"|15
 
|-
 
|-
|  style="border-bottom: 2pt solid black;text-align: center;"|segment
+
|  style="text-align: left;"|Segment
|  style="border-bottom: 2pt solid black;text-align: center;"|25
+
|  style="text-align: center;"|25
|  style="border-bottom: 2pt solid black;text-align: center;"|/
+
|  style="text-align: center;"|-
|  style="border-bottom: 2pt solid black;text-align: center;"|/
+
|  style="text-align: center;"|-
|  style="border-bottom: 2pt solid black;text-align: center;"|0.17
+
|  style="text-align: center;"|0.17
|  style="border-bottom: 2pt solid black;text-align: center;"|34500
+
|  style="text-align: center;"|34500
 
|}
 
|}
  
  
====3.2 Simulation of shield excavation process====
+
The boundary conditions and deformation constraints are crucial for accurately representing the real-world scenario in our simulation. To achieve this, we have implemented free boundaries at the top and surface of the levee to allow unrestricted movement. This setup enables us to capture the natural behavior of the levee’s surface under various loading conditions. Around the levee, we have applied deformation constraints to account for the complex interactions that occur in this region. These constraints help us model the interaction between the levee and its surroundings more accurately. Additionally, at the bottom surface, we have limited deformation in the <math display="inline">X, Y</math>, and <math display="inline">Z</math> directions to reflect the response of the ground to the construction process. This ensures that our simulation captures the realistic behavior of the ground beneath the levee. Furthermore, we have introduced lateral constraints on the sides of the model to account for the lateral support provided by the surrounding soil. These constraints play a crucial role in simulating the lateral stability of the levee under different loading conditions.
 +
 
 +
===4.2 Simulation of shield excavation process===
  
 
To simulate the step excavation process of the shield accurately, we utilize the element life and death function to manage the elements representing the shield machine, the grouting layer, and the segment. This process involves several key steps:
 
To simulate the step excavation process of the shield accurately, we utilize the element life and death function to manage the elements representing the shield machine, the grouting layer, and the segment. This process involves several key steps:
  
Initial Mud Support (NTH Analysis Step):In the NTH analysis step, we initiate the simulation by applying boundary conditions of mud pressure and pore pressure to the corresponding position of the palm surface, mimicking the mud support provided by the shield machine. This support pressure is continuously applied for one analysis step.
+
Geostatic (Initial Step): In the initial step, We first remove the elements of the shield machine, the grouting layer and the segment, replace them with corresponding displacement boundary conditions, and carry out static analysis of the obtained model under the state of self-weight, so as to obtain a stress field result, and then assign the obtained stress field as the initial stress field to the original model for in-situ stress balance.In this way, the influence of the initial soil stress on the subsequent calculation results can be eliminated [22,23,24].
  
Soil Excavation (n+1 Analysis Step):In the n+1 analysis step, we replicate the excavation of soil as the shield machine advances. Here, we employ the life-death unit method to deactivate the soil elements in front of the original palm face. The shield unit in front of the machine is activated, while the shield unit in the rear is deactivated. We apply the same mud pressure and pore pressure conditions to the next palm face to simulate the shield's progress. Additionally, we activate the grouting layer and segment units at this position and apply radial grouting pressure and pore pressure boundary conditions to simulate the grouting process at the shield tail.
+
Initial Mud Support (N Analysis Step): In the N analysis step, we initiate the simulation by applying boundary conditions of mud pressure and pore pressure to the corresponding position of the palm surface, mimicking the mud support provided by the shield machine. This support pressure is continuously applied for one analysis step.
  
Grout Hardening and Dissipation (n+6 Analysis Step):In the n+6 analysis step, we conclude the process by closing the boundary conditions of grouting pressure and pore pressure at this position. Furthermore, we enhance the modulus of the isoplasmic layer material to model grout hardening and grouting pressure dissipation.
+
Soil Excavation (N+1 Analysis Step): In the N+1 analysis step, we replicate the excavation of soil as the shield machine advances. Here, we employ the life-death unit method to deactivate the soil elements in front of the original palm face. The shield unit in front of the machine is activated, while the shield unit in the rear is deactivated. We apply the same mud pressure and pore pressure conditions to the next palm face to simulate the shield's progress. Additionally, we activate the grouting layer and segment units at this position and apply radial grouting pressure and pore pressure boundary conditions to simulate the grouting process at the shield tail.
 +
 
 +
Grout Hardening and Dissipation (N+6 Analysis Step): In the N+6 analysis step, we conclude the process by closing the boundary conditions of grouting pressure and pore pressure at this position. Furthermore, we enhance the modulus of the isoplasmic layer material to model grout hardening and grouting pressure dissipation.
  
 
These steps are cyclically executed to achieve dynamic tunnel shield excavation.
 
These steps are cyclically executed to achieve dynamic tunnel shield excavation.
  
====3.3 Calculation results and analysis====
+
===4.3 Calculation results and analysis===
  
3.3.1 Model verification
+
====4.3.1   Model verification====
  
Fig.7 presents a comparison between the measured settlement curve and the calculated settlement curve of the embankment top following the shield tunneling through the levee. This visual representation allows us to assess the accuracy of our calculations and their compliance with specified requirements.As depicted in Fig.7, the maximum calculated settlement value at the top of the levee, after the left hole shield tunnel's completion, is 8.7mm. Importantly, this value meets the specified requirements outlined in the project specifications. Furthermore, it is worth noting that the calculated settlement curve closely aligns with the measured settlement curve, indicating a high level of accuracy in both our numerical calculation model and the associated processes outlined in this paper.
+
In [[#img-7|Figure 7]], we compare the measured settlement curve with the calculated settlement curve of the embankment top following the shield tunneling through the levee. This visual representation allows us to assess the accuracy of our calculations and their compliance with specified requirements. The comparison reveals that the maximum calculated settlement value at the top of the levee, after the completion of the left hole shield tunnel, is 8.7mm. This value meets the specified requirements outlined in the project specifications. Furthermore, the calculated settlement curve closely aligns with the measured settlement curve, indicating a high level of accuracy in both our numerical calculation model and the associated processes outlined in this paper.
  
[[Image:Draft_Zhongxiang_374164140-image21.png|528px]]
+
<div id='img-7'></div>
 +
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"
 +
|-style="background:white;"
 +
|style="text-align: center;padding:10px;"| [[File:Review 537806263752-image10.png|467x467px]]
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 7'''. Measured settlement curve and calculated settlement curve of embankment top after <br> shield tunneling through the levee
 +
|}
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
Fig7. measured settlement curve and calculated settlement curve of embankment top after shield tunneling through the levee</div>
 
  
Fig.8 displays the calculated results of roof settlement in the levee during the penetration of the right-line shield. This visual representation allows us to assess the accuracy of our calculations and their compliance with project specifications.As depicted in the figure, the maximum measured settlement of the levee after the completion of the double-hole shield is 12.2mm. It's noteworthy that this value aligns with the specifications' requirements. While the calculated value is slightly smaller than the measured one, the agreement in the curve further substantiates the precision of our numerical calculation model.
+
[[#img-8|Figure 8]] presents the calculated results of roof settlement in the levee during the penetration of the right-line shield. In the figure, the maximum measured settlement of the levee after the completion of the double-hole shield is 12.2mm. Importantly, this value aligns with the requirements specified in the project's specifications. While the calculated value is slightly smaller than the measured one, the agreement in the curve further substantiates the precision of our numerical calculation model.
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
<div id='img-8'></div>
[[Image:Draft_Zhongxiang_374164140-image22.png|456px]] </div>
+
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"
 +
|-style="background:white;"
 +
|style="text-align: center;padding:10px;"| [[File:Review 537806263752-image11.png|402x402px]]
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 8'''. Calculation results of the roof settlement of the levee when the right line shield<br> penetrates through the levee
 +
|}
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
Fig. 8 Calculation results of the roof settlement of the levee when the right line shield penetrates through the levee</div>
 
  
====3.3.2 Deep soil settlement of levee====
+
Combined with [[#img-3|Figure 3]] and [[#img-8|Figure 8]], we analyze the factors contributing to the observed W-shaped settlement pattern. The settlement may be attributed to the following stages: 1) Before the shield arrives, the tunnel experiences overall subsidence due to the supporting pressure of the excavation face. 2) When the shield reaches the top of the left-line tunnel, the soil experiences unloading due to excavation and recovers elastic deformation, leading to a convex appearance at the top of the embankment. 3) As the shield progresses, the joint action of grouting pressure and other factors contributes to further settlement of the embankment top.
  
Fig.9 presents a comparative analysis of the lateral settling trough at the top of the levee and at depths of 10, 20, and 28 meters from the levee's top during the penetration of the left-line shield. Upon examining Fig.9, it becomes evident that the transverse settlement trough at varying depths of the levee exhibits a distinctive "V" shape, with the settling trough at the levee's summit showing a comparatively gentle slope. As we delve deeper, we observe a slight increase in soil settlement above the tunnel axis and a simultaneous reduction in the width of the settlement trough. This phenomenon can be attributed to the stress redistribution that occurs in the surrounding soil as the shield traverses the levee. This redistribution results in the formation of a lateral settlement trough. Notably, stress propagation at the top of the levee is more dispersed, leading to a less pronounced settlement effect. Moreover, with increasing depth, the soil above the tunnel axis gradually bears a greater load, contributing to a slight rise in settlement. The narrowing of the settlement trough width is a consequence of the localized nature of stress transfer.
+
====4.3.2  Deep soil settlement of levee====
  
The characteristics of deep soil settlement of the levee bear some resemblance to those induced by the φ6.34m subway shield <sup>[16]</sup>. However, it is noteworthy that the maximum soil settlement experienced by the levee and the extent to which the width of the sedimentation tank varies with depth are significantly less pronounced compared to the subway shield. In light of these observations, it is reasonable to consider that, under the specific conditions of this project, the width coefficient (i<sub>h</sub>) of the transverse sedimentation trough at varying depths closely approximates the width coefficient (i<sub>0</sub>) of the surface transverse sedimentation trough. When a large-diameter shield passes through the levee, the disparity between the settlement value at the levee's crest and the settlement value in the corresponding deep soil strata appears to be minimal. This suggests that, in practical engineering applications, the settlement value recorded at the levee's surface can effectively serve as an indirect indicator of deep soil settlement. Consequently, there may be no compelling need for the specialized arrangement of deep settlement monitoring points unless exceptionally high precision is warranted.
+
[[#img-9|Figure 9]] provides a comparative analysis of the lateral settling trough at the top of the levee and at depths of 10, 20, and 28 meters during the penetration of the left-line shield. The analysis reveals the following characteristics: 1) At varying depths, the transverse settlement trough exhibits a distinctive “V” shape, with the settling trough at the levee's summit showing a gentler slope. 2) Deeper depths show a slight increase in soil settlement above the tunnel axis and a simultaneous reduction in the width of the settlement trough. This phenomenon is attributed to stress redistribution in the surrounding soil as the shield traverses the levee, leading to the formation of a lateral settlement trough. 3) Stress propagation at the top of the levee is more dispersed, resulting in a less pronounced settlement effect. As depth increases, the soil above the tunnel axis bears a greater load, contributing to a slight rise in settlement, the narrowing of the settlement trough width is a consequence of the localized nature of stress transfer.
  
[[Image:Draft_Zhongxiang_374164140-image23.png|462px]]
+
<div id='img-9'></div>
 +
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"
 +
|-style="background:white;"
 +
|style="text-align: center;padding:10px;"| [[File:Review 537806263752-image12.png|462x462px]]
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 9'''. Lateral settlement curve of soil at different depths of levee
 +
|}
  
Fig 9 .Lateral settlement curve of soil at different depths of levee
 
  
Fig.10 presents a comparison of transverse settling troughs at various depths when the right-line shield passes through the levee. This analysis aims to shed light on the settlement patterns during this process. Upon examination, it's evident that the maximum settlement during the penetration of the right-line shield increases by 3.5 meters compared to the single-hole scenario. As we delve deeper into the analysis, several trends emerge: (a) With increasing depth, soil settlement above the tunnel axis experiences a slight increase. (b) The width of the sedimentation tank shows a minor decrease with depth. (c) Settlement in the shield area, overall, follows a "W" shape distribution.
+
The characteristics of deep soil settlement of the levee bear some resemblance to those induced by the φ6.34m subway shield [25]. However, it is noteworthy that the maximum soil settlement experienced by the levee and the extent to which the width of the sedimentation tank varies with depth are significantly less pronounced compared to the subway shield. In light of these observations, it is reasonable to consider that, under the specific conditions of this project, the width coefficient (<math display="inline">i_h</math>) of the transverse sedimentation trough at varying depths closely approximates the width coefficient (<math display="inline">i_0</math>) of the surface transverse sedimentation trough. When a large-diameter shield passes through the levee, the disparity between the settlement value at the levee's crest and the settlement value in the corresponding deep soil strata appears to be minimal. This suggests that, in practical engineering applications, the settlement value recorded at the levee’s surface can effectively serve as an indirect indicator of deep soil settlement. Consequently, there may be no compelling need for the specialized arrangement of deep settlement monitoring points unless exceptionally high precision is warranted.
  
Notably, at a depth of 28 meters below the top of the embankment (coinciding with the grouting layer), there's a pronounced "W" shape settlement distribution. This phenomenon can be explained by two key factors:
+
[[#img-10|Figure 10]] presents a comparison of transverse settling troughs at various depths when the right-line shield passes through the levee. This analysis aims to shed light on the settlement patterns during this process. Upon examination, it's evident that the maximum settlement during the penetration of the right-line shield increases by 3.5 meters compared to the single-hole scenario. As we delve deeper into the analysis, several trends emerge: 1) With increasing depth, soil settlement above the tunnel axis experiences a slight increase. 2) The width of the sedimentation tank shows a minor decrease with depth. 3) Settlement in the shield area, overall, follows a “W” shape distribution.
  
1. The depth of 28 meters coincides with the shield construction area for the right-line, resulting in significant soil disturbance near the existing left line. This leads to soil redistribution and substantial settlement.
+
<div id='img-10'></div>
 +
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"
 +
|-style="background:white;"
 +
|style="text-align: center;padding:10px;"| [[File:Review 537806263752-image13.png|408x408px]]
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 10'''. Lateral settling tank of soil in different depths of levee
 +
|}
 +
 
 +
 
 +
Notably, at a depth of 28 meters below the top of the embankment (coinciding with the grouting layer), there's a pronounced “W” shape settlement distribution. This phenomenon can be explained by two key factors:
 +
 
 +
# The depth of 28 meters coincides with the shield construction area for the right-line, resulting in significant soil disturbance near the existing left line. This leads to soil redistribution and substantial settlement.
 +
# Additionally, this depth aligns with the grouting layer location. Typically, grouting processes enhance soil strength, potentially reducing interaction between soil masses. Consequently, the soil exhibits a marked “W” shape settlement distribution at the grouting layer depth.
  
2. Additionally, this depth aligns with the grouting layer location. Typically, grouting processes enhance soil strength, potentially reducing interaction between soil masses. Consequently, the soil exhibits a marked "W" shape settlement distribution at the grouting layer depth.
 
  
 
In conclusion, this analysis reveals important settlement patterns during right-line shield penetration. These findings can have implications for future projects in similar conditions.
 
In conclusion, this analysis reveals important settlement patterns during right-line shield penetration. These findings can have implications for future projects in similar conditions.
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
====4.3.3  Horizontal displacement of deep soil mass of levee====
[[Image:Draft_Zhongxiang_374164140-image24.png|408px]] </div>
+
  
Fig. 10 Lateral settling tank of soil in different depths of levee
+
[[#img-11|Figure 11]] displays the horizontal displacement curves of the soil mass at the top of the levee and at depths of 10, 20, and 28 meters during the penetration of the left-line shield. As seen from [[#img-11|Figure 11]], the horizontal displacement curve of the levee top exhibits an approximately inverted “S” shaped distribution, symmetrically centered around the tunnel axis. The horizontal displacement directly above the axis is zero, while the horizontal displacement on both sides of the axis is non-zero and directed towards the tunnel axis. The maximum horizontal displacement is 6.6 millimeters, with symmetry appearing at positions <math display="inline">\pm i  </math> from the tunnel axis.
  
====3.3.3 Horizontal displacement of deep soil mass of levee====
+
# At a depth of 10 meters from the top of the levee, the soil is situated in a shallower position, significantly affected by surface loads and shield construction. The mechanical compression from the shield and lateral soil displacement result in substantial soil deformation, leading to a reduction in horizontal displacement by approximately 42.4% compared to the top position of the embankment.
 +
# As the depth increases to 20 meters, the soil is positioned deeper, with reduced influence from surface loads and construction, resulting in minimal soil deformation and a relatively modest reduction in horizontal displacement of around 1 millimeter.
 +
# At a depth of 28 meters from the top of the embankment, the soil is deep and in proximity to the shield grouting layer. The maximum horizontal displacement occurs at positions ±0.5i, with a reverse bend point at the centerline of the tunnel. Despite minimal impact from surface loads, shield construction induces significant deformation due to the enhanced compressive strength of the soil from the grouting layer. This results in an asymmetric displacement distribution, with some soil moving towards the tunnel direction near the shield side and away from the grouting layer.
 +
 
 +
 
 +
<div id='img-11'></div>
 +
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"
 +
|-style="background:white;"
 +
|style="text-align: center;padding:10px;"| [[File:Review 537806263752-image14.png|396x396px]]
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 11'''. Horizontal displacement of soil at different depths of the levee <br> when the left line is connected
 +
|}
 +
 
 +
 
 +
[[#img-12|Figure 12]] illustrates the horizontal displacement curves of the soil mass at the top of the levee and at depths of 10, 20, and 28 meters after shield tunneling. The horizontal displacement curve at the levee's crest exhibits an approximate “wave-like” distribution, symmetrically centered around the axis of the two tunnels. Above the axis, the horizontal displacement is zero, while on either side of the axis, it is non-zero and directed toward the tunnel axis. The maximum horizontal displacement recorded is 9.2 millimeters, symmetrically distributed at positions ±i from the tunnel axis. Distinct patterns in the variation of horizontal displacement become apparent at different depths:
 +
 
 +
# At a depth of 10 meters, the soil experiences significant influences from surface loads and shield construction, resulting in substantial deformation and a reduction in maximum horizontal displacement by approximately 42.4% compared to the top position.
 +
# At 20 meters, with reduced influence from surface loads and construction, soil deformation decreases, leading to a relatively modest reduction in horizontal displacement of approximately 1 millimeter.
 +
# At 28 meters, coinciding with the shield grouting layer, the soil displays unique behavior, with maximum horizontal displacement observed directly above the tunnel's centerline. The presence of the grouting layer enhances the soil's compressive strength, causing asymmetric displacement distribution.
 +
 
 +
 
 +
<div id='img-12'></div>
 +
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;"
 +
|-style="background:white;"
 +
|style="text-align: center;padding:10px;"| [[File:Review 537806263752-image15.png|402x402px]]
 +
|-
 +
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 12'''. Horizontal displacement of soil at different depths of levees<br> when the double lines are connected
 +
|}
  
Fig.11 displays the horizontal displacement curves of the soil mass at the top of the levee and at depths of 10, 20, and 28 meters from the top when the left-line shield penetrates the levee. As seen from Fig.11, the horizontal displacement curve of the levee top exhibits an approximately inverted "S"-shaped distribution, symmetrically centered around the tunnel axis. The horizontal displacement directly above the axis is zero, while the horizontal displacement on both sides of the axis is non-zero and directed towards the tunnel axis. The maximum horizontal displacement is 6.6 millimeters, with symmetry appearing at positions ±i from the tunnel axis. There are distinct variations in the laws of change for soils at different depths:
 
  
1. At a depth of 10 meters from the top of the levee, the soil is situated in a shallower position, making it significantly affected by surface loads and shield construction. The mechanical compression from the shield and the lateral soil displacement result in substantial soil deformation, leading to a reduction in horizontal displacement. Consequently, the maximum horizontal displacement of the soil decreases by approximately 42.4% compared to the top position of the embankment.
+
Real-time monitoring of horizontal displacement is crucial during the passage of a large-diameter shield tunnel, especially for deep soil layers within the levee. Additionally, when pile foundations or other retaining structures are present, proactive protective measures should be taken to anticipate additional shear forces and ensure levee safety.
  
2. As the depth increases to 20 meters, the soil is positioned deeper, and the influence of surface loads and construction decreases, resulting in minimal soil deformation. The reduction in horizontal displacement is also relatively modest, around 1 millimeter.
+
==5. Conclusions==
  
3. Concerning the depth of 28 meters from the top of the embankment, the soil is not only deep but also in proximity to the shield grouting layer. The maximum horizontal displacement occurs at positions ±0.5i, with a reverse bend point at the centerline of the tunnel. Despite minimal impact from surface loads, shield construction still induces significant deformation. The presence of the grouting layer likely enhances the soil's compressive strength, causing some soil to move towards the tunnel direction near the shield side and away from the grouting layer, resulting in an asymmetric displacement distribution.
+
1. When the left line shield penetrates through the levee, the maximum settlement of the top reaches 8.7mm, which meets the deformation control requirements. The surface transverse settlement trough is distributed in a “V” shape, and the width of the sedimentation trough on the top of the levee is about 5 times of the width coefficient of the sedimentation trough. The settlement distribution can be estimated by the Peck formula. When the levee is passed under the right line, the maximum settlement of the top reaches 12.2mm, which has a wider influence on the soil and a larger settlement.
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
2. When the left line shield tunneling underpasses the levee, the lateral settlement at different depths of the levee is distributed in a “V” shape, and the lateral settlement trough at the top of the levee is relatively gentle. With the increase of depth, the soil settlement above the tunnel axis slightly increases, while the width of the settlement trough slightly decreases. When the right line shield underpasses the levee, the lateral settlement at different depths of the levee shows a “W” shape distribution, and the “W” shape settlement of the soil at the depth of the grouting layer is more significant.
[[Image:Draft_Zhongxiang_374164140-image25.png|396px]] </div>
+
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
3. The horizontal displacement curves of the top of the levee and deep soil mass (except the grouting layer) are approximately inverted “S” shape when the left line shield underpasses the levee and the right line shield underpasses the levee, and the maximum horizontal displacement of the top of the levee and deep soil mass (except the grouting layer) occurs at the position of <math display="inline">\pm i</math>. And the maximum horizontal displacement of the right shield through the levee is greater than the maximum horizontal displacement of the left shield through the levee.
Fig.11 Horizontal displacement of soil at different depths of the levee when the left line is connected</div>
+
  
Fig.12 illustrates the horizontal displacement curves of the soil mass at the top of the levee and at depths of 10, 20, and 28 meters from the top after shield tunneling. As depicted in the figure, the horizontal displacement curve at the levee's crest exhibits an approximate "wave-like" distribution, symmetrically centered around the axis of the two tunnels. Above the axis, the horizontal displacement registers as zero, whereas on either side of the axis, the horizontal displacement is non-zero and directs toward the tunnel axis. The maximum horizontal displacement recorded is 9.2 millimeters, symmetrically distributed at positions ±i from the tunnel axis.
+
4. When the left line shield underpasses the levee, the maximum horizontal displacement at the grouting layer appears at the position of <math display="inline">\pm 0.5i</math>, and there is a reverse bending point at the center line of the tunnel. When the right line shield tunneling underpasses the levee, the maximum horizontal displacement at the grouting layer appears at the position of <math display="inline">\pm 0.5i</math>, and there are reverse bending points at the center line, <math display="inline">\pm i</math> and <math display="inline">\pm 0.75i</math> of the tunnel.
  
For soil at different depths, distinct patterns in the variation of horizontal displacement become apparent. At a depth of 10 meters from the levee's top, the soil, being in a shallower position, experiences significant influences from surface loads and shield construction. The mechanical compression from the shield and lateral soil displacement result in substantial soil deformation, leading to a reduction in maximum horizontal displacement of approximately 42.4% compared to the top position.
+
==References==
  
As depth increases to 20 meters, the soil is situated deeper, resulting in reduced influence from surface loads and construction activities. Consequently, soil deformation decreases, and the reduction in horizontal displacement is relatively modest, approximately 1 millimeter.
+
<div class="auto" style="text-align: left;width: auto; margin-left: auto; margin-right: auto;font-size: 85%;">
  
At a depth of 28 meters from the levee's top, coinciding with the position of the shield grouting layer, the soil displays a unique behavior. Maximum horizontal displacement is observed directly above the centerline of the tunnel. Additionally, there are inflection points in soil horizontal displacement at ±0.75i and ±i. Despite the minimal impact of surface loads, shield construction continues to induce significant soil deformation. The presence of the grouting layer appears to enhance the soil's compressive strength, causing some soil to move toward the tunnel and away from the grouting layer. This results in an asymmetric distribution of displacement.
+
[1] Xie W., Wu W. Analysis of measured settlement of the Qiantang River embankment caused by under-crossing of large-diameter slurry shield tunneling. Modern Tunneling Technology, 48(4):142-147, 2011.
  
To ensure the safety of the levee, real-time monitoring of horizontal displacement is imperative during the passage of a large-diameter shield tunnel, particularly for the deep soil layers within the levee. In cases where pile foundations or other retaining structures are present in the levee area, additional shear forces should be anticipated, necessitating proactive protective measures.
+
[2] Yang J. Finite difference analysis of impact of undercrossing super-large-diameter shield tunnel on training jetty. Chinese Journal of Geotechnical Engineering, 41(8):1569-1572, 2019.
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
[3] Zhang Y., Xia P., Zhu W., et al. Analysis of stratum deformation and stability: a case study of the Nanjing Weisan road river-crossing tunnel passing under the Yangtze River embankment. Modern Tunneling Technology, 52(4):151-157, 2015.
[[Image:Draft_Zhongxiang_374164140-image26.png|402px]] </div>
+
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
[4] Han L., Ye G., Wang J., et al. Finite element analysis of impact of under-crossing of large shallow shield tunnel on riverbank. Chinese Journal of Geotechnical Engineering, 37(S1):125-128, 2015.
Fig.12 Horizontal displacement of soil at different depths of levees when the double lines are connected</div>
+
  
===4. Conclusion===
+
[5] Li Z., Huang X. Study on settlement control for slurry shields crossing the embankment. Modern Tunneling Technology, 48(1):103-110, 2011.
  
(1) When the left line shield penetrates through the levee, the maximum settlement of the top reaches 8.7mm, which meets the deformation control requirements. The surface transverse settlement trough is distributed in a "V" shape, and the width of the sedimentation trough on the top of the levee is about 5 times of the width coefficient of the sedimentation trough. The settlement distribution can be estimated by the Peck formula. When the levee is passed under the right line, the maximum settlement of the top reaches 12.2mm, which has a wider influence on the soil and a larger settlement.
+
[6] Zhang Z., Lin C., Wu S., et al. Analysis and control of ground settlement of embankments in construction of cross-river shield tunnels. Chinese Journal of Geotechnical Engineering, 33(6):977-984, 2011.
  
(2) When the left line shield tunneling underpasses the levee, the lateral settlement at different depths of the levee is distributed in a "V" shape, and the lateral settlement trough at the top of the levee is relatively gentle. With the increase of depth, the soil settlement above the tunnel axis slightly increases, while the width of the settlement trough slightly decreases. When the right line shield underpasses the levee, the lateral settlement at different depths of the levee shows a "W" shape distribution, and the "W" shape settlement of the soil at the depth of the grouting layer is more significant.
+
[7] Wu W., Sun Y., Li L. et al. Analysis of levee deformation caused by large diameter shield tunneling in river Crossing Tunnel. China Railway Science, 37(04):78-82, 2016.
  
(3) The horizontal displacement curves of the top of the levee and deep soil mass (except the grouting layer) are approximately inverted "S" shape when the left line shield underpasses the levee and the right line shield underpasses the levee, and the maximum horizontal displacement of the top of the levee and deep soil mass (except the grouting layer) occurs at the position of ±i. And the maximum horizontal displacement of the right shield through the levee is greater than the maximum horizontal displacement of the left shield through the levee.
+
[8] Huang H. Research on settlement of Qiantang River levee caused by large-diameter shield construction. Zhejiang, Zhejiang University of Technology, 2018.
  
(4) When the left line shield underpasses the levee, the maximum horizontal displacement at the grouting layer appears at the position of ±0.5i, and there is a reverse bending point at the center line of the tunnel. When the right line shield tunneling underpasses the levee, the maximum horizontal displacement at the grouting layer appears at the position of ±0.5i, and there are reverse bending points at the center line, ±0.75i and ±i of the tunnel.
+
[9] Zhu M., Wei L., Fang Y. et al. Study on settlement analysis and control of large diameter shield tunnel under the Yellow River levee. Modern Tunnel Technology, 59(03):211-219, 2022.
  
===Reference===
+
[10] Xie X., Yang C., Wang Q. et al. Settlement analysis and control study of shield crossing levee of Yangtze River in Nanjing and Yanlu River Crossing Passage. Chinese Journal of Rock Mechanics and Engineering, 40(S2):3313-3322, 2021.
  
[1] XIE Wenbin,WU Weiyi. Analysis of measured settlement of the Qiantang River embankment caused by under-crossing of large-diameter slurry shield tunneling[J]. Modern Tunneling Technology,2011,48(4):142–147.
+
[11] Lin  Q., Cao P., Dong T. Numerical simulation of flood control levee deformation by shield tunnel excavation through the Liuyang River: A case study. Geofluids, 2023:156-163, 2023.
  
[2] YANG Jiangang. Finite difference analysis of impact of undercrossing super-large-diameter shield tunnel on training jetty[J]. Chinese Journal of Geotechnical Engineering,2019,41(8):1569–1 572.
+
[12] Wu S., Lin C., Zhang Z., et al. Risk analysis and control for slurry shield under-passing embankment. Chinese Journal of Rock Mechanics and Engineering, 30(5):1034-1042, 2011.
  
[3] ZHANG Yazhou,XIA Pengju,ZHU Wei,et al. Analysis of stratum deformation and stability: a case study of the Nanjing Weisan road river-crossing tunnel passing under the Yangtze River embankment[J]. Modern Tunneling Technology,2015,52(4):151–157.
+
[13] Epel T.,  Mooney M.A., Gutierrez N. The influence of face and shield annulus pressure on tunnel liner load development. Tunnelling and Underground Space Technology, 117, 104096, 2021.
  
[4] HAN Lei,YE Guanlin,WANG Jianhua,et al. Finite element analysis of impact of under-crossing of large shallow shield tunnel on riverbank[J]. Chinese Journal of Geotechnical Engineering,2015,37(Supp.1):125–128.
+
[14] Lee H.,  Kwak J.,  Choi J.,  Hwang B., Choi H. A lab-scale experimental approach to evaluate rheological properties of foam-conditioned soil for EPB shield tunnelling. Tunnelling and Underground Space Technology, 128, 104667, 2022.
  
[5] LI Zongliang,HUANG Xigang. Study on settlement control for slurry shields crossing the embankment[J]. odern Tunneling Technology,2011,48(1):103–110.
+
[15] Davoodi M.,  Senent S., Keshavarz A.,  Jimenez R. Three-dimensional seismic face stability of shield tunnels in undrained clay, Underground Space, 15:26-43, 2024.
  
[6] ZHANG Zhongmiao,LIN Cungang,WU Shiming,et al. Analysis and control of ground settlement of embankments in construction of cross-river shield tunnels[J]. Chinese Journal of Geotechnical Engineering,2011,33(6):977–984.
+
[16] Zizka Z.,  Schoesser B.,  Thewes M. Investigations on transient support pressure transfer at the tunnel face during slurry shield drive part 1: Case A – Tool cutting depth exceeds shallow slurry penetration depth. Tunnelling and Underground Space Technology, 118, 104168, 2021.
  
[7] Wu Weiyi, Sun Yukun, Li Liang et al. Analysis of levee deformation caused by large diameter shield tunneling in river Crossing Tunnel [J]. China Railway Science, 2016,37(04):78~82.
+
[17] Zizka Z., Schoesser B., Thewes M. Investigations on the transient support pressure transfer at the tunnel face during slurry shield drive part 2: Case B – Deep slurry penetration exceeds tool cutting depth. Tunnelling and Underground Space Technology, 118, 104169, 2021.
  
[8] Huang Hai. Research on settlement of Qiantang River levee caused by large-diameter shield construction [D] Zhejiang: Zhejiang University of Technology, 2018.
+
[18] Seol H., Won D.,  Jang J.,  Yeom Kim K.,  Sup Yun T. Ground Collapse in EPB shield TBM site: A case study of railway tunnels in the deltaic region near Nak-Dong River in Korea. Tunnelling and Underground Space Technology, 120, 104274, 2022.
  
[9] ZHU Muyuan, Wei Lifeng, Fang Yong et al. Study on settlement analysis and control of large diameter shield tunnel under the Yellow River levee [J]. Modern Tunnel Technology. 2022,59(03):211~219.
+
[19] Peck R.B. Deep excavations and tunneling in soft ground. State of the Art Report, Proceedings of 7th International Conference on Soil Mechanics and Foundation Engineering, Mexico City, 225-290, 1969.
  
[10] Xie Xiongyao, Yang Changzhi, Wang Qiang et al. Settlement analysis and control study of shield crossing levee of Yangtze River in Nanjing and Yanlu River Crossing Passage [J]. Chinese Journal of Rock Mechanics and Engineering. 201,40(S2) :3313~3322.
+
[20] O’Reilly M.P., New B.M. Settlements above tunnels in the United Kingdom-their magnitude and prediction. Proc. Tunnelling 82, Institution of Mining and Metallurgy, London,  173-182, 1982.
  
[11] Lin  Qibin ; Cao Ping ; Dong Tao. Numerical Simulation of Flood Control Levee Deformation by Shield Tunnel Excavation through the Liuyang River: A Case Study. GEOFLUIDS[J] Vol(2023):156~163.
+
[21] Jiang P.W., Zhang Z.H., Zheng H., et al. Coupling analysis method of grouting construction with deformation response of adjacent existing tunnel. Underground Space, 15:312-330, 2024.
  
[12] WU Shiming,LIN Cungang,ZHANG Zhongmiao,et al. Risk analysis and control for slurry shield under-passing embankment[J]. Chinese Journal of Rock Mechanics and Engineering,2011,30(5):1 034– 1 042.
+
[22] Abdellah W.R., Haridy A.K.A., Mohamed A.K., et al. Behaviour of horseshoe-shaped tunnel subjected to different in situ stress fields. Applied Sciences-Basel, 12(11), 5399, 2022.
  
[13] PECK R B. Deep excavations and tunneling in soft ground[C]// State of the art report. Proceedings of 7th International Conference on Soil Mechanics and Foundation Engineering. Mexico City: 1969: 225–290.
+
[23] Nguyen T.T., Do N.A., Karasev M.A., et al. Influence of tunnel shape on tunnel lining behaviour. Proceedings of the Institution of Civil Engineers-Geotechnical Engineering, 174(4):355-371, 2021.
  
[14] O’REILLY M P, NEW B M. Settlements above tunnels in the United Kingdom-their magnitude and prediction[C]// Proc. Tunnelling 82, Institution of Mining and Metallurgy, London: 173–182.
+
[24] Do N.A., Dias D., Zhang Z.X., et al. Study on the behavior of squared and sub-rectangular tunnels using the Hyperstatic Reaction Method. Transportation Geotechnics, 22, 100321, 2020.
  
[15] Tang Yiqun,Ye Weimin,Zhang Qinghe.Study of Ground Settlement Caused by Subway Tunnel Construction in Shanghai(Ⅲ) [J].Underground Space,1995,15(4):250-258.
+
[25] Sun Y., Zhou S., Gong Q. Distribution of deep displacement field during shield tunneling in soft-soil areas. Chinese Journal of Rock Mechanics and Engineering, 28(3):500-506, 2009.

Latest revision as of 13:13, 25 March 2024

Abstract

This study utilizes ABAQUS finite element software to analyze embankment deformation during shield tunneling. Results reveal that shield tunneling effects are manageable. The left tunnel line successfully tunnels beneath the embankment, causing a maximum settlement of 8.7 mm, meeting regulations. The ground surface exhibits a “V” shaped lateral settlement trough, approximately five times the tunnel centerline width. Conversely, the right tunnel line induces a more extensive impact, with a 12.2 mm maximum settlement. Different lateral settlement patterns emerge in soil at varying depths. The left tunnel creates a “V” shape, with slightly increasing settlement above the tunnel axis as depth increases, accompanied by narrower troughs. The right tunnel line results in a pronounced “W” shape, especially at the grouting layer depth. Horizontal displacement forms an “S” shape, with maximum displacement at positions ±i from the tunnel axis. The right tunnel line induces greater maximum horizontal displacement than the left. In summary, this numerical analysis provides insights into embankment deformation during shield tunneling, aiding in assessing the impact and guiding safe tunneling beneath embankments.

Keywords: Shield tunnel, finite element simulation, levee settlement, data fitting

1. Introduction

With the rapid development of urban transportation, the number and scale of urban cross-river tunnel constructions are also expanding. These tunnels, primarily built using large-diameter mud shield methods, have seen their outer diameters grow from approximately 11 metersto about 15 meters [1,2]. Notably, these river-crossing tunnels must traverse through river embankments, inevitably causing soil disturbance, particularly during the construction of super-large diameter mud shield tunnels. Given the intricacy of embankment structures and the vital role that river embankments play in safeguarding the lives and property of communities on both sides of the river, it becomes crucial to adopt rational construction parameters and measures. This approach aims to minimize the impact of shield tunneling on embankment structures.

In recent years, the engineering field has accumulated significant experience regarding the disruptive effects of shield tunneling through levees. Zhang et al. employed finite-difference software to construct a three-dimensional numerical model of the Nanjing Weisan Road super-diameter shield tunnel, which crosses the Yangtze River embankment. They calculated the settlement pattern and overall stability of the embankment [3]. Building upon the context of the Chengjiang West Road tunnel in Jiangyin, Han Lei et al. employed the displacement control finite element method, which is based on the formation loss rate, to compute the embankment settlement under various formation loss rates. They also proposed recommendations for control [4]. Utilizing monitoring data and the Peck empirical formula, Li and Huang examined the disruption and settlement patterns of the Qiantang River levee during the shield crossing construction of the Qingchun Road Tunnel in Hangzhou City [5]. They also suggested settlement control measures for the mud shield crossing the levee [5,6]. Wu et al. utilized the FLAC3D finite difference program in conjunction with field measurements to compute and analyze the deformation patterns of the levee resulting from the penetration of the Qiantang River levee by the mud-water balanced shield, which has a diameter of 11.68m [7]. Their study revealed that the transverse settling trough at the top of the levee resembled that of a subway shield with a diameter of 6.34m, displaying a Gaussian normal distribution. The Peck formula remains applicable for estimating settlement; however, the change in settlement volume and trough width with depth is not as pronounced as observed with the 6.34m diameter subway shield. Huang Hai investigated the surface settlement of the embankment caused by large-diameter shield tunneling under the Qiantang River. The study revealed that the transverse settlement trough of the levee generally follows a normal distribution pattern, with the greatest settlement occurring at the center of the tunnel axis. Settlement diminishes gradually as distance from the tunnel center increases, and settlement distribution along the tunnel axis is symmetrical, in line with the normal distribution curve proposed by Huang [8]. Relying on the Jilou Road Tunnel project, Zhu et al. analyzed the settlement patterns and characteristics of the large embankment section during the shield tunneling process through the levee using field-measured data. They employed a detailed dynamic simulation of the shield tunnel driving process and observed that the initial settlement rate of the large embankment section during the shield tail ejection stage was greater and more susceptible to disturbance compared to that of an ordinary flat section [9]. Moreover, the top of the levee is particularly sensitive to the disturbance caused by deep soil. Xie et al. utilized finite element software to construct a three-dimensional numerical model for the super-diameter shield tunneling across the Yangtze River levee [10]. They conducted a detailed simulation of the construction process, considering construction elements such as soil excavation, mud-water support, segment lining, grouting behind the wall, and formation loss. Lin et al. based their study on the Liuyang River levee engineering section of Changsha Metro Line 6, They investigated the deformation of the levee under shield tunnel excavation conditions [11]. Wu et al. analyzed the risk factors associated with mud shield tunneling through levee engineering [12]. They examined the causes, hazards, and avoidance measures of these risks and validated the effectiveness of the risk control measures by referencing the engineering example of mud shield tunneling through the levee in the Qingchun Road Tunnel in Hangzhou. Tamir et al. simplified the modeling of the tunnel boring machine shield while incorporating the hardened soil constitutive model [13]. This approach enabled them to capture crucial aspects of the shield tunneling process, such as joint segment lining, grouting pressure, time-varying grout hardening, and nonlinear soil behavior. Hyobum et al. conducted laboratory pressure blade shear tests to examine the rheological properties of sandy soils under foam conditions [14]. They investigated how varying water content, foam injection ratio (FIR), and confining pressure affected the rheological properties of soil under foam conditions. Milad et al. utilized the finite element limit analysis method to assess the limit support pressure in undrained clay [15]. Their evaluation considered horizontal pseudo-static seismic forces directed outward from the face. Zdenek et al. presented a new approach for incorporating the formation process into the analysis of pressure transfer [16,17]. This method was specifically designed for shallow slurry penetration depths smaller than the cutting depth of the tools. Hoonil et al. reported on the significant ground subsidence that ensued after the structural deformation of tunnel segments during the excavation of a cross-passage between two main tunnels [18]. These main tunnels were excavated by an earth pressure balance shielded tunnel boring machine (EPB shield TBM) and were constructed in the deltaic region near the Nak-Dong River in Korea.

Most of the existing research focuses on the analysis of levee settlement, with few studies addressing the settlement and horizontal displacement characteristics of deep soil in levees. Therefore, analyzing the three-dimensional deformation characteristics of levees after shield tunneling is beneficial for understanding the patterns and managing construction risks. In this paper, the deformation (settlement and displacement) of the levee caused by the large-diameter shield machine during the construction of the Jiangshan Road tunnel in Hangzhou is numerically calculated, and the three-dimensional deformation characteristics of the levee are studied.

2. Project overview and geological conditions

2.1 Project overview

The Genshandong Road Cross-River Tunnel is a significant infrastructure project situated between the Xisha Bridge and Jiangdong Bridge, spanning a total length of 4462.26m. This tunnel plays a vital role in the region's transportation network, and its technical aspects are noteworthy for a comprehensive understanding. This shield tunnel features a single-layer lining with segment assembly. It boasts an impressive outer diameter of 14.5m, an inner diameter of 13.3m, a segment thickness of 0.6m, and a ring width of 2m. The tunnel comprises 10 “9+1” sections, and these sections are assembled with staggered joints, tightly secured by high-strength bolts. The construction of the shield tunnel, a significant infrastructure project in the region, begins at the eastern starting well of the Big River and proceeds westward, reaching its destination at the Xisha receiving well. However, the tunnel faces a major challenge when it crosses the mighty Qiantang River due to the high deformation requirements of the levee caused by tidal influence. For a visual representation of the current status of the levee, as shown in Figure 1.

Review 537806263752-image1.png Review 537806263752-image2.png Review 537806263752-image3.png
Figure 1. Present situation of levees on the east bank of Qiantang River


2.2 Hydrogeological conditions and monitoring arrangement

The shield tunnel passes through various soil layers, including sandy silt, silt, silty clay, mucky clay, clay, Fully weathered rock, and Moderately weathered rock. Despite the levee not being strengthened before shield construction, strategic monitoring points were installed on the levee to ensure safety. These measures are crucial for maintaining project stability and safety. For detailed information about the tunnel's soil layers and monitoring setup, as shown in Figure 2.

Review 537806263752-image4.png
Figure 2. Monitoring site layout profile

3. Measured analysis of surface settlement of levees

3.1 Lateral surface settlement of levee top

In Figure 3, we observe the settlement patterns above and around the tunnel axis. The most significant settlement occurs directly above the tunnel axis. For the left-line shield tunnel, the transverse surface settlement generally forms a “V” shape distribution, with a maximum settlement of 8.7mm directly above the left line.

Conversely, after the right-line shield tunnel crosses the levee, the transverse surface settlement generally exhibits a “W” shape distribution. The maximum settlement reaches 12.2mm, which is still directly above the left line of the shield tunnel and greater than the surface settlement above the right line. Additionally, the settlement in the middle of the left and right lines is slightly less than the settlement value above the tunnel.

Review 537806263752-image5.png
Figure 3. Transverse surface settlement curve of levee top after the left and right lines cross the levee


For the prediction of ground displacement caused by single-line tunnel excavation, the horizontal distribution estimation formula of ground settlement proposed by Peck is generally adopted in current engineering practice [19]:

(1)
(2)

where is the ground subsidence caused by ground loss; the distance from the tunnel center line; the maximum settlement caused by the formation loss at the center line of the tunnel; the width coefficient of surface settling trough; the stratum loss per unit length of tunnel; the formation loss rate, is the ratio of the surface subsidence trough area to the tunnel excavation area; and is the radius of tunnel excavation.

When the Peck formula is used to predict the transverse surface settlement, the most important thing is to determine the formation loss rate and the width coefficient of sedimentation tank . O'Reilly and New [20] proposed the following values in accordance with the experience of London:

(3)

where is the sedimentation tank width parameter; and the buried depth for the tunnel axis.

In the paper, the Peck formula is utilized to fit the measured transverse surface settlement curve of the embankment top after the left line crosses the embankment. This analysis aims to compare the surface settlement of the embankment top and evaluate the applicability of the Peck formula in the project, as shown in Figure 4. The findings from this analysis are crucial for understanding the behavior of the embankment top and its response to the tunneling process.

Review 537806263752-image7.png
Figure 4. Comparison between the measured value of land surface settlement and the Peck fitting value


Upon analyzing Figure 4, we observe that the variation in the surface settlement curve closely aligns with the Peck formula. However, it's noteworthy that the measured values in the field are slightly larger than those fitted by the Peck formula. The discrepancy in measured values may arise due to the greater structural weight of the levee compared to typical strata. This difference is reflected in the higher soil and water pressure, as well as settlement, beneath the levee excavation surface at the same burial depth compared to shield tunneling in other general strata. The large diameter of the shield used influences these factors. Despite the challenges posed by large-diameter shield tunnels, the Peck formula remains applicable, allowing engineers to predict the distribution of transverse settling troughs at the top of the levee. These predictions, which consider factors such as formation loss (Vs) and the width coefficient of the settling trough (i), are invaluable for assessing and managing levee settlement and deformation during construction.

3.2 Surface settlement of shield through various stages

These stages provide a framework for understanding the settlement process during shield tunnel construction and are crucial for monitoring and managing potential impacts on surrounding structures.

  1. The shield approaches: the shield excavation face has not reached the monitoring section, and the distance is about 2 times the diameter of the shield.
  2. The shield arrives: the excavation face of the shield has not reached the monitoring section, and the distance is 0~2 times the diameter of the shield.
  3. The shield goes through: when the shield excavation face passes through the monitoring section, the distance is 1 times the diameter of the shield.
  4. The shield leaves: when the shield excavation face leaves the monitoring section about 1 times the diameter of the shield.
  5. Final stable times: about 10 times the diameter of the shield when the shield excavation face leaves the monitoring section.


As shown in Figure 5, this is the lateral surface settlement curve of the levee at each stage of left-line shield construction. By comparing these curves, it can be observed that as the shield advances, surface settlement gradually increases. The width of the transverse settlement trough in the same section remains relatively constant, and the settlement curve becomes steeper. When the shield reaches the monitoring fault front (approximately 1D~2D away from the section, D: the tunnel diameter), a noticeable transverse settlement trough with a width of about 15~20 meters is generated. The maximum land settlement occurs at the center of the tunnel, with settlement decreasing as one moves farther away from the tunnel center.

Review 537806263752-image8.png
Figure 5. Surface settlement curve of the left line shield crossing the levee at each stage


Table 1 provides the percentage distribution of settlement in different stages of shield tunnel excavation. In the early settlement stage of the monitoring point, the settlement value accounts for 7% of the total settlement, which, while relatively small, should not be overlooked. Moving on to the arrival stage of the shield, we see that it contributes 14% of the total settlement. However, the phases of shield passage and post-shield passage are particularly critical, with these two stages contributing to a larger proportion of the total settlement, specifically 43% and 23%, respectively. In fact, settlements in these stages collectively account for approximately 66% of the total settlement value, making this phase the most concerning due to its rapid and substantial settlement changes. Finally, the late sedimentation stage produces around 13% of the total settlement value. To summarize, the percentage of settlement in various stages of shield tunnel excavation varies significantly. The early and late stages contribute relatively small percentages, while the arrival stage accounts for 14% of the total settlement. However, it's during the phases of shield passage and post-shield passage that settlement percentages are most substantial, combining to make up approximately 66% of the total settlement value. This phase is of particular concern due to its rapid and pronounced settlement changes.

Table 1. Percentage of the total settlement at each stage of the monitoring point (%)
Pre-arrival settlement of shield Phase settlement of shield arrival Phase settlement of the shield as it passes through Phase settlement after shield passage Long-term subsequent stage settlement
7 14 43 23 13

4. Numerical analysis of levee deformation

To gain insights into the soil deformation characteristics of the levee during sequential passage by the left and right-line shields, we conducted an investigation using the advanced finite element software ABAQUS. This analysis aims to calculate and compare the dynamic construction processes of the shields to understand how they impact the soil and the overall construction of the tunnel.

4.1 Numerical calculation model

To gain a comprehensive understanding of the real-world engineering scenario and design parameters, we have developed a numerical calculation model using ABAQUS finite element software. This model is designed to simulate the effects of the shield tunnel's construction, which comprises two parallel lines with a centerline separation of 35m. One of the key features of the model is its representation of the grouting filler within the shield tunnel gap as a uniform and equally thick grouting layer.

In our simulation, we have incorporated the thrust generated by the shield machine and the pressure exerted by cement in the form of nodal concentrated forces. Our model encompasses various elements, including the surrounding rock, pipe segments, grouting material, and the levees, all represented using solid elements.

The numerical model’s dimensions are as follows: a longitudinal length of 56m, a transverse width of 180m, and a vertical height ranging from 76.1m to 84.6m. This extensive model is divided into 140,208 units and comprises 148,071 nodes, ensuring a comprehensive representation of the system. For a visual representation, as shown in Figure 6.

Review 537806263752-image9.png
Figure 6. Numerical calculation model


In shield construction, a ring gap will be formed between the tunnel excavation wall and the outer ring of the segment when the segment is detached from the shield tail. If the action of grouting and pressure is ignored, the soil mass will move to this gap and produce a large displacement.In the actual grouting process, the slurry fills the gap of the shield tail and penetrates into the soil body, and then gradually hardens. If the hardening of the grouting material is not considered, the settlement value will also deviate from the actual value by using the isogeneration layer simulation. Because the interaction of shield tail grouting will form a mixture of soil and cement slurry, the grouting layer can be simulated by elastic model, and the hardening of grouting material can be simulated by the variable stiffness of grouting material [21]. And an elastic model is utilized for soil. The specific material parameters for these models are detailed in Table 2.

Table 2. Material parameters of the model
Material name Weight (kN/m³) Cohesion (kPa) Internal friction Angle (°) Poisson's ratio Elastic modulus (MPa)
Sandy silt 19.8 12 32 0.2 19
Silt 19.9 13.5 33.4 0.22 20
Silty clay 17.6 19 32 0.25 8
Mucky clay 18.4 27 34 0.35 4
Clay 19.6 45 17.3 0.36 9
Fully weathered rock 22.3 50 22 0.26 320
Moderately weathered rock 22.6 200 22 0.23 550
Grouting layer 21 - - 0.2 15
Segment 25 - - 0.17 34500


The boundary conditions and deformation constraints are crucial for accurately representing the real-world scenario in our simulation. To achieve this, we have implemented free boundaries at the top and surface of the levee to allow unrestricted movement. This setup enables us to capture the natural behavior of the levee’s surface under various loading conditions. Around the levee, we have applied deformation constraints to account for the complex interactions that occur in this region. These constraints help us model the interaction between the levee and its surroundings more accurately. Additionally, at the bottom surface, we have limited deformation in the , and directions to reflect the response of the ground to the construction process. This ensures that our simulation captures the realistic behavior of the ground beneath the levee. Furthermore, we have introduced lateral constraints on the sides of the model to account for the lateral support provided by the surrounding soil. These constraints play a crucial role in simulating the lateral stability of the levee under different loading conditions.

4.2 Simulation of shield excavation process

To simulate the step excavation process of the shield accurately, we utilize the element life and death function to manage the elements representing the shield machine, the grouting layer, and the segment. This process involves several key steps:

Geostatic (Initial Step): In the initial step, We first remove the elements of the shield machine, the grouting layer and the segment, replace them with corresponding displacement boundary conditions, and carry out static analysis of the obtained model under the state of self-weight, so as to obtain a stress field result, and then assign the obtained stress field as the initial stress field to the original model for in-situ stress balance.In this way, the influence of the initial soil stress on the subsequent calculation results can be eliminated [22,23,24].

Initial Mud Support (N Analysis Step): In the N analysis step, we initiate the simulation by applying boundary conditions of mud pressure and pore pressure to the corresponding position of the palm surface, mimicking the mud support provided by the shield machine. This support pressure is continuously applied for one analysis step.

Soil Excavation (N+1 Analysis Step): In the N+1 analysis step, we replicate the excavation of soil as the shield machine advances. Here, we employ the life-death unit method to deactivate the soil elements in front of the original palm face. The shield unit in front of the machine is activated, while the shield unit in the rear is deactivated. We apply the same mud pressure and pore pressure conditions to the next palm face to simulate the shield's progress. Additionally, we activate the grouting layer and segment units at this position and apply radial grouting pressure and pore pressure boundary conditions to simulate the grouting process at the shield tail.

Grout Hardening and Dissipation (N+6 Analysis Step): In the N+6 analysis step, we conclude the process by closing the boundary conditions of grouting pressure and pore pressure at this position. Furthermore, we enhance the modulus of the isoplasmic layer material to model grout hardening and grouting pressure dissipation.

These steps are cyclically executed to achieve dynamic tunnel shield excavation.

4.3 Calculation results and analysis

4.3.1 Model verification

In Figure 7, we compare the measured settlement curve with the calculated settlement curve of the embankment top following the shield tunneling through the levee. This visual representation allows us to assess the accuracy of our calculations and their compliance with specified requirements. The comparison reveals that the maximum calculated settlement value at the top of the levee, after the completion of the left hole shield tunnel, is 8.7mm. This value meets the specified requirements outlined in the project specifications. Furthermore, the calculated settlement curve closely aligns with the measured settlement curve, indicating a high level of accuracy in both our numerical calculation model and the associated processes outlined in this paper.

Review 537806263752-image10.png
Figure 7. Measured settlement curve and calculated settlement curve of embankment top after
shield tunneling through the levee


Figure 8 presents the calculated results of roof settlement in the levee during the penetration of the right-line shield. In the figure, the maximum measured settlement of the levee after the completion of the double-hole shield is 12.2mm. Importantly, this value aligns with the requirements specified in the project's specifications. While the calculated value is slightly smaller than the measured one, the agreement in the curve further substantiates the precision of our numerical calculation model.

Review 537806263752-image11.png
Figure 8. Calculation results of the roof settlement of the levee when the right line shield
penetrates through the levee


Combined with Figure 3 and Figure 8, we analyze the factors contributing to the observed W-shaped settlement pattern. The settlement may be attributed to the following stages: 1) Before the shield arrives, the tunnel experiences overall subsidence due to the supporting pressure of the excavation face. 2) When the shield reaches the top of the left-line tunnel, the soil experiences unloading due to excavation and recovers elastic deformation, leading to a convex appearance at the top of the embankment. 3) As the shield progresses, the joint action of grouting pressure and other factors contributes to further settlement of the embankment top.

4.3.2 Deep soil settlement of levee

Figure 9 provides a comparative analysis of the lateral settling trough at the top of the levee and at depths of 10, 20, and 28 meters during the penetration of the left-line shield. The analysis reveals the following characteristics: 1) At varying depths, the transverse settlement trough exhibits a distinctive “V” shape, with the settling trough at the levee's summit showing a gentler slope. 2) Deeper depths show a slight increase in soil settlement above the tunnel axis and a simultaneous reduction in the width of the settlement trough. This phenomenon is attributed to stress redistribution in the surrounding soil as the shield traverses the levee, leading to the formation of a lateral settlement trough. 3) Stress propagation at the top of the levee is more dispersed, resulting in a less pronounced settlement effect. As depth increases, the soil above the tunnel axis bears a greater load, contributing to a slight rise in settlement, the narrowing of the settlement trough width is a consequence of the localized nature of stress transfer.

Review 537806263752-image12.png
Figure 9. Lateral settlement curve of soil at different depths of levee


The characteristics of deep soil settlement of the levee bear some resemblance to those induced by the φ6.34m subway shield [25]. However, it is noteworthy that the maximum soil settlement experienced by the levee and the extent to which the width of the sedimentation tank varies with depth are significantly less pronounced compared to the subway shield. In light of these observations, it is reasonable to consider that, under the specific conditions of this project, the width coefficient () of the transverse sedimentation trough at varying depths closely approximates the width coefficient () of the surface transverse sedimentation trough. When a large-diameter shield passes through the levee, the disparity between the settlement value at the levee's crest and the settlement value in the corresponding deep soil strata appears to be minimal. This suggests that, in practical engineering applications, the settlement value recorded at the levee’s surface can effectively serve as an indirect indicator of deep soil settlement. Consequently, there may be no compelling need for the specialized arrangement of deep settlement monitoring points unless exceptionally high precision is warranted.

Figure 10 presents a comparison of transverse settling troughs at various depths when the right-line shield passes through the levee. This analysis aims to shed light on the settlement patterns during this process. Upon examination, it's evident that the maximum settlement during the penetration of the right-line shield increases by 3.5 meters compared to the single-hole scenario. As we delve deeper into the analysis, several trends emerge: 1) With increasing depth, soil settlement above the tunnel axis experiences a slight increase. 2) The width of the sedimentation tank shows a minor decrease with depth. 3) Settlement in the shield area, overall, follows a “W” shape distribution.

Review 537806263752-image13.png
Figure 10. Lateral settling tank of soil in different depths of levee


Notably, at a depth of 28 meters below the top of the embankment (coinciding with the grouting layer), there's a pronounced “W” shape settlement distribution. This phenomenon can be explained by two key factors:

  1. The depth of 28 meters coincides with the shield construction area for the right-line, resulting in significant soil disturbance near the existing left line. This leads to soil redistribution and substantial settlement.
  2. Additionally, this depth aligns with the grouting layer location. Typically, grouting processes enhance soil strength, potentially reducing interaction between soil masses. Consequently, the soil exhibits a marked “W” shape settlement distribution at the grouting layer depth.


In conclusion, this analysis reveals important settlement patterns during right-line shield penetration. These findings can have implications for future projects in similar conditions.

4.3.3 Horizontal displacement of deep soil mass of levee

Figure 11 displays the horizontal displacement curves of the soil mass at the top of the levee and at depths of 10, 20, and 28 meters during the penetration of the left-line shield. As seen from Figure 11, the horizontal displacement curve of the levee top exhibits an approximately inverted “S” shaped distribution, symmetrically centered around the tunnel axis. The horizontal displacement directly above the axis is zero, while the horizontal displacement on both sides of the axis is non-zero and directed towards the tunnel axis. The maximum horizontal displacement is 6.6 millimeters, with symmetry appearing at positions from the tunnel axis.

  1. At a depth of 10 meters from the top of the levee, the soil is situated in a shallower position, significantly affected by surface loads and shield construction. The mechanical compression from the shield and lateral soil displacement result in substantial soil deformation, leading to a reduction in horizontal displacement by approximately 42.4% compared to the top position of the embankment.
  2. As the depth increases to 20 meters, the soil is positioned deeper, with reduced influence from surface loads and construction, resulting in minimal soil deformation and a relatively modest reduction in horizontal displacement of around 1 millimeter.
  3. At a depth of 28 meters from the top of the embankment, the soil is deep and in proximity to the shield grouting layer. The maximum horizontal displacement occurs at positions ±0.5i, with a reverse bend point at the centerline of the tunnel. Despite minimal impact from surface loads, shield construction induces significant deformation due to the enhanced compressive strength of the soil from the grouting layer. This results in an asymmetric displacement distribution, with some soil moving towards the tunnel direction near the shield side and away from the grouting layer.


Review 537806263752-image14.png
Figure 11. Horizontal displacement of soil at different depths of the levee
when the left line is connected


Figure 12 illustrates the horizontal displacement curves of the soil mass at the top of the levee and at depths of 10, 20, and 28 meters after shield tunneling. The horizontal displacement curve at the levee's crest exhibits an approximate “wave-like” distribution, symmetrically centered around the axis of the two tunnels. Above the axis, the horizontal displacement is zero, while on either side of the axis, it is non-zero and directed toward the tunnel axis. The maximum horizontal displacement recorded is 9.2 millimeters, symmetrically distributed at positions ±i from the tunnel axis. Distinct patterns in the variation of horizontal displacement become apparent at different depths:

  1. At a depth of 10 meters, the soil experiences significant influences from surface loads and shield construction, resulting in substantial deformation and a reduction in maximum horizontal displacement by approximately 42.4% compared to the top position.
  2. At 20 meters, with reduced influence from surface loads and construction, soil deformation decreases, leading to a relatively modest reduction in horizontal displacement of approximately 1 millimeter.
  3. At 28 meters, coinciding with the shield grouting layer, the soil displays unique behavior, with maximum horizontal displacement observed directly above the tunnel's centerline. The presence of the grouting layer enhances the soil's compressive strength, causing asymmetric displacement distribution.


Review 537806263752-image15.png
Figure 12. Horizontal displacement of soil at different depths of levees
when the double lines are connected


Real-time monitoring of horizontal displacement is crucial during the passage of a large-diameter shield tunnel, especially for deep soil layers within the levee. Additionally, when pile foundations or other retaining structures are present, proactive protective measures should be taken to anticipate additional shear forces and ensure levee safety.

5. Conclusions

1. When the left line shield penetrates through the levee, the maximum settlement of the top reaches 8.7mm, which meets the deformation control requirements. The surface transverse settlement trough is distributed in a “V” shape, and the width of the sedimentation trough on the top of the levee is about 5 times of the width coefficient of the sedimentation trough. The settlement distribution can be estimated by the Peck formula. When the levee is passed under the right line, the maximum settlement of the top reaches 12.2mm, which has a wider influence on the soil and a larger settlement.

2. When the left line shield tunneling underpasses the levee, the lateral settlement at different depths of the levee is distributed in a “V” shape, and the lateral settlement trough at the top of the levee is relatively gentle. With the increase of depth, the soil settlement above the tunnel axis slightly increases, while the width of the settlement trough slightly decreases. When the right line shield underpasses the levee, the lateral settlement at different depths of the levee shows a “W” shape distribution, and the “W” shape settlement of the soil at the depth of the grouting layer is more significant.

3. The horizontal displacement curves of the top of the levee and deep soil mass (except the grouting layer) are approximately inverted “S” shape when the left line shield underpasses the levee and the right line shield underpasses the levee, and the maximum horizontal displacement of the top of the levee and deep soil mass (except the grouting layer) occurs at the position of . And the maximum horizontal displacement of the right shield through the levee is greater than the maximum horizontal displacement of the left shield through the levee.

4. When the left line shield underpasses the levee, the maximum horizontal displacement at the grouting layer appears at the position of , and there is a reverse bending point at the center line of the tunnel. When the right line shield tunneling underpasses the levee, the maximum horizontal displacement at the grouting layer appears at the position of , and there are reverse bending points at the center line, and of the tunnel.

References

[1] Xie W., Wu W. Analysis of measured settlement of the Qiantang River embankment caused by under-crossing of large-diameter slurry shield tunneling. Modern Tunneling Technology, 48(4):142-147, 2011.

[2] Yang J. Finite difference analysis of impact of undercrossing super-large-diameter shield tunnel on training jetty. Chinese Journal of Geotechnical Engineering, 41(8):1569-1572, 2019.

[3] Zhang Y., Xia P., Zhu W., et al. Analysis of stratum deformation and stability: a case study of the Nanjing Weisan road river-crossing tunnel passing under the Yangtze River embankment. Modern Tunneling Technology, 52(4):151-157, 2015.

[4] Han L., Ye G., Wang J., et al. Finite element analysis of impact of under-crossing of large shallow shield tunnel on riverbank. Chinese Journal of Geotechnical Engineering, 37(S1):125-128, 2015.

[5] Li Z., Huang X. Study on settlement control for slurry shields crossing the embankment. Modern Tunneling Technology, 48(1):103-110, 2011.

[6] Zhang Z., Lin C., Wu S., et al. Analysis and control of ground settlement of embankments in construction of cross-river shield tunnels. Chinese Journal of Geotechnical Engineering, 33(6):977-984, 2011.

[7] Wu W., Sun Y., Li L. et al. Analysis of levee deformation caused by large diameter shield tunneling in river Crossing Tunnel. China Railway Science, 37(04):78-82, 2016.

[8] Huang H. Research on settlement of Qiantang River levee caused by large-diameter shield construction. Zhejiang, Zhejiang University of Technology, 2018.

[9] Zhu M., Wei L., Fang Y. et al. Study on settlement analysis and control of large diameter shield tunnel under the Yellow River levee. Modern Tunnel Technology, 59(03):211-219, 2022.

[10] Xie X., Yang C., Wang Q. et al. Settlement analysis and control study of shield crossing levee of Yangtze River in Nanjing and Yanlu River Crossing Passage. Chinese Journal of Rock Mechanics and Engineering, 40(S2):3313-3322, 2021.

[11] Lin Q., Cao P., Dong T. Numerical simulation of flood control levee deformation by shield tunnel excavation through the Liuyang River: A case study. Geofluids, 2023:156-163, 2023.

[12] Wu S., Lin C., Zhang Z., et al. Risk analysis and control for slurry shield under-passing embankment. Chinese Journal of Rock Mechanics and Engineering, 30(5):1034-1042, 2011.

[13] Epel T., Mooney M.A., Gutierrez N. The influence of face and shield annulus pressure on tunnel liner load development. Tunnelling and Underground Space Technology, 117, 104096, 2021.

[14] Lee H., Kwak J., Choi J., Hwang B., Choi H. A lab-scale experimental approach to evaluate rheological properties of foam-conditioned soil for EPB shield tunnelling. Tunnelling and Underground Space Technology, 128, 104667, 2022.

[15] Davoodi M., Senent S., Keshavarz A., Jimenez R. Three-dimensional seismic face stability of shield tunnels in undrained clay, Underground Space, 15:26-43, 2024.

[16] Zizka Z., Schoesser B., Thewes M. Investigations on transient support pressure transfer at the tunnel face during slurry shield drive part 1: Case A – Tool cutting depth exceeds shallow slurry penetration depth. Tunnelling and Underground Space Technology, 118, 104168, 2021.

[17] Zizka Z., Schoesser B., Thewes M. Investigations on the transient support pressure transfer at the tunnel face during slurry shield drive part 2: Case B – Deep slurry penetration exceeds tool cutting depth. Tunnelling and Underground Space Technology, 118, 104169, 2021.

[18] Seol H., Won D., Jang J., Yeom Kim K., Sup Yun T. Ground Collapse in EPB shield TBM site: A case study of railway tunnels in the deltaic region near Nak-Dong River in Korea. Tunnelling and Underground Space Technology, 120, 104274, 2022.

[19] Peck R.B. Deep excavations and tunneling in soft ground. State of the Art Report, Proceedings of 7th International Conference on Soil Mechanics and Foundation Engineering, Mexico City, 225-290, 1969.

[20] O’Reilly M.P., New B.M. Settlements above tunnels in the United Kingdom-their magnitude and prediction. Proc. Tunnelling 82, Institution of Mining and Metallurgy, London, 173-182, 1982.

[21] Jiang P.W., Zhang Z.H., Zheng H., et al. Coupling analysis method of grouting construction with deformation response of adjacent existing tunnel. Underground Space, 15:312-330, 2024.

[22] Abdellah W.R., Haridy A.K.A., Mohamed A.K., et al. Behaviour of horseshoe-shaped tunnel subjected to different in situ stress fields. Applied Sciences-Basel, 12(11), 5399, 2022.

[23] Nguyen T.T., Do N.A., Karasev M.A., et al. Influence of tunnel shape on tunnel lining behaviour. Proceedings of the Institution of Civil Engineers-Geotechnical Engineering, 174(4):355-371, 2021.

[24] Do N.A., Dias D., Zhang Z.X., et al. Study on the behavior of squared and sub-rectangular tunnels using the Hyperstatic Reaction Method. Transportation Geotechnics, 22, 100321, 2020.

[25] Sun Y., Zhou S., Gong Q. Distribution of deep displacement field during shield tunneling in soft-soil areas. Chinese Journal of Rock Mechanics and Engineering, 28(3):500-506, 2009.
Back to Top

Document information

Published on 25/03/24
Accepted on 03/03/24
Submitted on 06/09/23

Volume 40, Issue 1, 2024
DOI: 10.23967/j.rimni.2024.03.003
Licence: CC BY-NC-SA license

Document Score

0

Views 0
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?