Next Article in Journal
Modified Accuracy of RANS Modeling of Urban Pollutant Flow within Generic Building Clusters Using a High-Quality Full-Scale Dispersion Dataset
Next Article in Special Issue
Laboratory Experimental Study on the Pressure Relief Effect of Boreholes in Sandstone under High-Stress Conditions
Previous Article in Journal
Study on Spatial and Temporal Differences of Water Resource Sustainable Development and Its Influencing Factors in the Yellow River Basin, China
Previous Article in Special Issue
Study on the Stress Distribution and Stability Control of Surrounding Rock of Reserved Roadway with Hard Roof
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Surface Settlement Induced by Shield Tunnelling: Grey Relational Analysis and Numerical Simulation Study on Critical Construction Parameters

1
School of Civil Engineering, Qingdao University of Technology, Qingdao 266520, China
2
School of Mechanical and Electrical Engineering, China University of Mining and Technology, Beijing 100083, China
3
State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Beijing 100083, China
*
Authors to whom correspondence should be addressed.
Sustainability 2023, 15(19), 14315; https://doi.org/10.3390/su151914315
Submission received: 28 August 2023 / Revised: 25 September 2023 / Accepted: 26 September 2023 / Published: 28 September 2023
(This article belongs to the Special Issue Advancing Sustainability in Geotechnical Engineering)

Abstract

:
Excessive surface settlement poses significant challenges to shield tunnelling construction, resulting in damage to adjacent buildings, infrastructure, and underground pipelines. This study focused on investigating the surface settlement induced by shield tunnelling during the construction of Qingdao Metro Line 6 between Haigang Road Station and Chaoyang Road Station. Firstly, the settlement data from the left line of the shield tunnel were evaluated by grey relational analysis. The relational coefficients were calculated to assess the correlation degrees of each influential parameter. Subsequently, the four critical influential parameters with the highest relational degrees were chosen to investigate their effects on surface settlement through numerical simulations under different scenarios. The results show that the four parameters with the highest relational degrees were thrust, grouting pressure, earth pressure, and strata elastic modulus. It should be noted that the strata elastic modulus significantly affects surface settlement, while the grouting pressure influences the settlement trough width in weak strata. Moreover, improper thrust magnitude can lead to an increase in surface settlement. Based on these findings, recommendations are proposed for the right-line tunnel construction and practical countermeasures for surface settlement during shield tunnelling construction are provided.

1. Introduction

Rapidly escalating urbanization calls for efficient transport systems and has made metro construction crucial to meet growing urban transportation demand. The shield tunnelling method has become the primary construction method for underground metro projects, owing to the rapid construction speed, automated operation, and minimal impact on surface traffic. However, shield tunnelling inevitably disturbs the surrounding rocks, leading to stress redistribution within the soil [1] and inducing surface settlement [2,3]. The excessive surface settlement can damage adjacent buildings and infrastructure [4,5], jeopardizing the safety of pedestrians, vehicles, and roads [6]. This in turn can lead to work stoppages and result in additional costs [7]. Moreover, it may trigger soil liquefaction and landslides, leading to severe construction accidents. Therefore, it is essential to understand and predict the behavior of surface settlement induced by shield tunnelling to ensure the safety and stability of urban areas.
An approach using Gaussian normal distribution to represent the lateral surface settlement induced by shield tunnelling was initially proposed by Peck [8]. The prediction formulas to describe surface settlement and its influential parameters were introduced and improved by subsequent researchers [9,10,11,12,13,14]. Notably, these formulas typically rely on empirical parameters and require local measurements for theoretical adjustments. In addition, some researchers have also studied surface settlement through model tests [15,16,17], but these tests have limitations in accurately imitating field conditions and actual construction scenarios. For greater authenticity and credibility, researchers utilized the measured data to analyze the correlation between the surface settlement and the influential construction parameters. For instance, Kannangara et al. [18] investigated the impact of shield tunnelling parameters on surface settlement. They identified the correlation between earth pressure and grouting pressure on the peak of surface settlement. Based on the construction data from EPB-TBM in the Moscow metro, Ter-Martirosyan et al. [19] studied the effect of metro tunnel construction parameters on surface settlement. The results indicated a close mutual relationship between earth pressure and surface settlement. Dongku Kim et al. [20] analyzed data from the construction section of the Hong Kong Metro and identified parameters such as boring speed, thrust, and cutterhead torque presented a modest correlation with the surface settlement. However, it cannot offer specific functional relationships or in-depth analysis of the complex interactions between the construction parameters and the surface settlement.
In recent years, numerical simulation methods have been widely employed to solve the interdisciplinary problems encountered in underground engineering, involving the surrounding rock deformation affected by novel metallic rock support materials [21,22] and the surface settlement issues induced by shield tunnelling [23,24,25,26,27,28,29,30]. Wang et al. [31] employed particle flow to simulate surface settlement in sandy soils, revealing its dynamic nature involving soil arch formation and destruction. Alzabeebee et al. [32] studied traffic load, microtunnel diameter, and backfill height using a finite element model. They found traffic loads significantly impact microtunnel-induced pavement settlement. Through finite element analysis, Yuan et al. [33] investigated the effect of parameters like deformation modulus, burial depth, and tunnel spacing on surface settlement during shield construction in weathered mudstone areas. Numerical simulations can theoretically model the entire shield construction process. However, existing simulations frequently suffer from inaccurate selection and validation of essential simulation parameters due to the multitude of construction variables. In addition, these simulations typically oversimplify the surface settlement process by focusing on a single and homogeneous stratum, leaving room for improvement to consider complex geological conditions and various influential parameters.
Existing data analysis methods have limitations in quantifying settlement impact patterns, and numerical simulations encounter challenges in selecting critical construction parameters. Hence, this study proposes a novel hybrid numerical approach combined with grey relational analysis to analyze the critical construction parameters influencing surface settlement comprehensively. The grey relational analysis serves to quantify the correlation degree of various parameters with surface settlement by calculating grey relational coefficients [34]. When coupled with numerical simulations, it facilitates optimized parameter selection and enhances the assessment of specific parameters’ impact on surface settlement. The engineering background was based on the shield construction section of Qingdao Metro Line 6 between Haigang Road Station and Chaoyang Road Station. This tunnel section presents complex geological conditions with diverse soil properties and various stratigraphic characteristics, where severe settlements occurred during the construction of the left line. In this study, the settlement monitoring data from the left line were evaluated by grey relational analysis to assess the relational degrees of each influential parameter. Subsequently, numerical simulations were conducted to investigate the effect of critical construction parameters on surface settlement, and recommendations are proposed for the right-line tunnelling. Finally, practical countermeasures for surface settlement during shield tunnelling construction are provided based on the research findings to ensure the safety and stability of future urban subway projects.

2. Project Overview

The shield construction section of Qingdao Metro Line 6, between Haigang Road Station and Chaoyang Road Station, comprises two single-line single-tunnel sections, starting and ending at Y(Z)DK 31 + 248.393~32 + 039.183. The right-line section is 790.790 m in length, while the left is 791.732 m. The depth of vault ranges from 12.9 m to 25.3 m. CR Equip. Machine 179, produced by China Railway Engineering Equipment Group Co., Ltd., headquartered in Zhengzhou, China, was employed for this tunnel section. The shield machine excavated and constructed the tunnel from Chaoyang Road Station and advanced toward Haigang Road Station, featuring an excavation diameter of ϕ = 6303 mm. The cutterhead comprised six spoke-type cutters and six face plates, totaling 34 face cutters, 12 gauge cutters, and six center cutters. The main parameters of the shield machines are detailed in Table 1.
The topography of the construction section slopes from northwest to southeast. The YDK31 +248.393~+339.682 section exhibits an erosional slope landform, while the YDK31 +339.682~32 + 039.183 section displays an erosional accumulation and gently sloping landform. Within the section, there are six standard layers and 19 sublayers. The left line of the construction section traverses several main strata, including powdery clay, clayey gravelly sand, strongly weathered granite porphyry, moderately weathered granite porphyry, slightly weathered granite porphyry, blocky granite (fractured), strongly weathered tuff, moderately weathered tuff, and sandy tuff (fractured). The groundwater in the section mainly consists of quaternary pore water and bedrock fissure water. Additionally, the construction section is influenced by the Lingshanwei fracture and its secondary fracture. The fault’s exposure along the Lingshanwei area shows a northeast strike of 40° to 55°, predominantly northwest tendency, and a dip of 70° to 88°. This section reveals three fracture zones and one severely affected zone of the Lingshanwei fault. Table 2 provides more detailed stratigraphic physical and mechanical parameters, sourced from the “Geotechnical Investigation Report of Qingdao Metro Line 6 Phase I Project, Haigang Road Station—Chaoyang Road Station,” supplied by Qingdao Municipal Research Institute of Surveying and Mapping. The geological profile of the left-line construction section is illustrated in Figure 1. Figure 2 and Figure 3 depict the layout of monitoring points and an aerial view in the vicinity of Haigang Road Station, respectively.

3. Grey Relational Analysis

Grey relational analysis is a quantitative method used to assess the correlation of trends among different parameters in a complex system. It is particularly suitable for systems with multiple independent variables and unclear interconnections. The method involves the following key steps.
  • Determine the evaluation index sequence and the influential parameter sequence.
  • Transform the obtained sequences into dimensionless form.
  • Calculate the difference sequences and identify the maximum and minimum differences.
  • Compute the relational coefficients and determine the grey relational degree.
The above four steps enable researchers to evaluate the influence of different factors on the evaluation index and facilitate the subsequent optimization analysis.

3.1. Model Index Determination

When calculating the relational degree in the grey relational model, it is necessary to determine the evaluation index sequence (dependent variables) and the influential parameter sequence (independent variables) that characterize the system behavior. The selection of the evaluation index sequence directly affects the rationality of the relational degree calculation results:
X 0 = x 0 1 x 0 2 x 0 n
X i = x i 1 x i 2 x i n
where X 0 is the evaluation index sequence and X i is the influential parameter sequence.
Previous studies [35,36,37,38,39] have identified several influential factors for the soil movement and surface settlement induced by shield tunnelling:
  • Changes in the stress state of the excavated soil.
  • Squeezing and forward movement of the soil in front of the excavation face and the surrounding soil in contact with the shield.
  • Over-excavation during the tunnelling process.
  • Formation of gaps around the shield due to the shield shell radius being smaller than the cutterhead radius, causing the surrounding soil to move towards the void and result in surface settlement.
  • Improper grouting at the shield tail leads to soil compression and settlement.
  • Interaction between the soil and the lining, as well as between the lining and the segment.
  • Creep or secondary consolidation of the soft clay soil following construction disturbance around the shield tunnel, leading to continuous secondary consolidation settlement.
Among the above influential factors, factors 1, 2, and 3 affect the stability of the excavation face and can be controlled by adjusting tunnelling parameters during construction. Factors 4, 5, and 6 are related to grouting quality and can be controlled by adjusting grouting parameters. Factor 7 can be considered stratum parameters.
Based on the above correlations between the influential factors and the surface settlement, this study selected the surface settlement as the evaluation index sequence in the relational degree calculation. To ensure the measuring data’s coherence, surface settlement data were obtained from the monitoring points above the centerline of the left-line tunnel section. The influential parameter sequences selected were earth pressure, tunnelling speed, and thrust as tunnelling parameters; grouting pressure and grouting volume as grouting parameters; and the weighted average of cohesion, internal friction angle, elastic modulus, and Poisson’s ratio as stratum parameters. The impact on surface settlement was evaluated by analyzing the relational degrees of these parameters. Table 3 summarizes the influential construction parameters and corresponding surface settlement data collected and collated for each ring during the left-line tunnelling.

3.2. Grey Relational Degree Calculation

The grey relational degree calculation quantifies the relative influences and trends of selected influential parameter sequences compared to the evaluation index sequence. To account for the varying dimensions and magnitudes of the system’s influential parameters and minimize errors in the calculation, mean normalization is performed on the original data columns for dimensionless standardization:
X 0 k = X 0 k / X ¯ 0 k = 1,2 , , n
X i k = X i k / X ¯ i k = 1,2 , , n
where X ¯ 0 and X ¯ i represent the average values of the evaluation index sequence and influential parameter sequences, respectively.
The absolute differences between all index sequences and the influential parameter sequence are then calculated using the following formula:
Δ 0 i k = X 0 k X i k i = 1,2 , , n
Subsequently, the maximum and minimum absolute differences are obtained with the following formulas:
Δ m a x = M a x i M a x k Δ 0 i k
Δ m i n = M i n i M i n k Δ 0 i k
where Δ 0 i k is the absolute difference between the ith influential parameter sequence and the evaluation index sequence at the kth data point, and Δ m a x and Δ m i n represent the maximum and minimum values in the set, respectively.
Next, the relational coefficient and grey relational degree are calculated using the following expressions:
l 0 i = Δ m i n + ρ Δ m a x Δ 0 i ( k ) + ρ Δ m a x
r 0 i = 1 n k = 1 n l 0 i
where l 0 i is the relational coefficient between the evaluation index sequence and the ith influential parameter sequence at the kth data point. The resolution coefficient ρ is introduced, with smaller ρ values indicating greater resolution ( ρ ∈ (0,1), and in this paper, ρ = 0.5). r 0 i stands for the relational degree between the evaluation index sequence and the kth influential parameter sequence, with a value domain of (0,1).

3.3. Results and Analysis

The calculation results of relational coefficients for each data point are shown in Figure 4 and Table 4. The overall relational degrees of influential construction parameters are shown in Table 5.
Analysis of Figure 4, Table 4 and Table 5 yield the following observations:
  • Grey relational coefficients above 0.7 indicate significant correlations between certain parameters and surface settlement. Figure 4a shows that the relation of cohesion and internal friction angle with the surface settlement is relatively low within the Data ID range of 9–15 (corresponding to tunnelling rings 280–247). Similarly, in the interval of Data ID 34–38 (corresponding to tunnelling rings 65–41), the relation between Poisson’s ratio and surface settlement is also notably low, which corresponds to fractured zone areas. This indicates a weak correlation between stratum parameters (excluding elastic modulus) and surface settlement within the fractured zone region.
  • Earth pressure exhibits the highest grey relational degree among the construction parameters. Combined with Table 4 and Table 5, in Data ID 38–34, which exhibits high relational coefficients, surface settlement decreases as earth pressure gradually increases, consistent with previous studies [18]. Notably, slightly different from other findings [20], the grey correlation degree of tunnelling speed does not reach a significant level. Tunnelling speed briefly demonstrates a strong relational coefficient only when maintained near its maximum value (Data ID 25–30, corresponding to tunnelling rings 183–116). Therefore, flexibly adjusting the tunnelling speed based on the coefficients during right-line tunnelling helped mitigate surface settlement.
  • While the grouting parameters did not exhibit the highest relational degrees, they still held prominent positions [40]. Notably, grouting pressure contributes slightly more to surface settlement than grouting volume. In addition, as shown in Figure 4c, the relational coefficients between the two are higher in the fractured zone than in other intervals. Therefore, ensuring the stability of grouting pressure and controlling the grouting volume are both critical for achieving the optimal grouting effect.

4. Numerical Simulation

Grey relational analysis reveals the relational degrees of different parameters with the surface settlement. However, it cannot offer specific functional relationships or in-depth analysis of the complex interactions between the construction parameters and the surface settlement. To explore the variation patterns governing the surface settlement influenced by the critical construction parameters, numerical simulations were conducted to analyze the effect of the four parameters (thrust, grouting pressure, earth pressure, and strata elastic modulus) with the highest relational degrees. In this study, we utilized Flac3D as a numerical simulation tool, which is based on the finite difference method and designed for addressing underground geotechnical problems. It employs an implicit solution algorithm, seamlessly combining the finite difference method with the analysis of two-body dynamics to accurately simulate complex underground phenomenon.

4.1. Simulation Setup

The following simplifying assumptions were made to simulate the deformation behaviors of rock and soil zones in the model.
  • The tunnel structure and surrounding rock were assumed uniform and isotropic, exhibiting typical elastoplastic behavior.
  • The initial stress field was induced only by the self-weight of the surrounding rock, regardless of the potential influence of groundwater on the stress distribution.
  • The deformation of the rock and soil zones follows the Mohr–Coulomb strength criterion. Throughout the simulation, no progressive damage, nonlinear bending behavior, or loading–unloading cyclic loading occurred. Therefore, the lining and grouting zones were assumed as elastic material, and the shield shell zone was treated as a rigid material.
Besides the above simplifying assumptions, the following adjustments were also made to enhance the relevance to the actual conditions.
  • Normal constraints were set on the surroundings and bottom of the model to simulate actual constraints on the surrounding rocks. No constraints were applied on the top.
  • A fixed time-step value was implemented in the simulation, ensuring that the parameters of the grouting layer within the simulation varied in sync with the time step (time).
  • The simulation considered the standard gravitational acceleration (9.8 m/s2), and the convergence threshold of the computational iterations was set at 10−5 to maintain numerical simulation accuracy.

4.2. Simulation Modelling

Based on the geological profile of Haigang Road Station to Chaoyang Road Station, the tunnelling rings from 140 to 180 (corresponding to Data ID 27–23) were selected as the representative interval to construct the numerical model. The stratum of the selected interval consists of powdery clay, clayey gravelly sand, strongly weathered tuff, moderately weathered tuff, slightly weathered tuff, and blocky granite (fractured). This interval lies within the transition zone between granite and weathered tuff, adjacent to fractured zones, encompassing both weak and stable strata. Such transitional zones frequently exhibit nonuniform physical properties like particle size distribution, porosity, and crack distribution, presenting intricate rock traits and mechanical behavior. Therefore, selecting this interval for simulation can provide comprehensive influence patterns of parameter variations under different conditions.
According to Saint-Venant’s principle, the width of the model along the X direction was set to 3–5 times the tunnel diameter to avoid boundary effects, resulting in a total width of 30 m. Similarly, the height of the model along the Z direction was 30 m, and the length of the model along the Y direction was 60 m. Additionally, to make the simulation more representative of the actual tunnelling process, grouting, lining, and segment zones were introduced to consider the influence of support processes and grouting parameters on the strata. The thickness of the shield shell zones was set to 60 mm, and the thickness of the lining zones was set to 300 mm.
To determine the optimal meshing scheme, we conducted simulations with different grid densities. The scheme that showed no significant result changes with further refinement was chosen. The final mesh scheme divided the model into grids with a 0.375 m interval, with additional refinement at the geological transitions and tunnel excavation sections. The three-dimensional numerical model is shown in Figure 5 and consists of 373,722 zones, encompassing 274,797 grid points.

4.3. Simulation Scheme

Based on the maximum and minimum values of actual construction parameters during shield tunnelling, the range for the four critical parameters were set as follows.
  • The earth pressure ranges in stress from 20 to 138 kPa.
  • The grouting pressure ranges in stress from 134 to 334 kPa.
  • The thrust ranges in force from 6886 to 11,460 kN.
  • The elastic modulus changes from 50% to 150% (compared with the exploration data).
Each parameter was set at five levels, resulting in 17 construction simulation schemes, which shown in Table 6. Each one of critical parameters varied within the specified range, while the remaining parameters were kept at the mean values. In each simulation scheme, the entire numerical simulation process was conducted using a ring (1.5 m) as the smallest circulation unit. The procedure involved the following steps.
  • Initially, the shield shell zones of the shield machine were advanced by 1.5 m to simulate the initiation of a new shield tunnelling loop.
  • Subsequently, the tunnel zones, segment zones, and grouting zones were removed to mimic the excavation process of shield construction.
  • After the release of surrounding rock stresses, the earth pressure, the thrust, and the grouting pressure were applied to replicate the grouting process of shield construction.
  • Following the grouting process, the grouting pressure was canceled, and material parameters were assigned to grouting and the segment zones to simulate the assembly process of shield construction.
  • The parameters of the grouting zones were defined as a function of time, gradually transitioning from pre-solidification parameters to post-solidification parameters with each calculation step.
  • The calculation was then allowed to converge, completing one cycle of excavation.
  • Subsequently, the shield shell zones of the shield machine were advanced by another 1.5 m, initiating the next tunnelling cycle.
This cyclic tunnelling process continued until the model was fully constructed. The geotechnical zone’s parameters in the model were adopted from the geological survey report (refer to Table 2), while the material parameters for the tunnel and shield shell are presented in Table 7.

4.4. Model Validation and Analysis

Simulation scheme 3, where all critical construction parameters were given the mean values, was as per the initial choice numerical simulation. Subsequently, the simulation results were compared with the actual monitoring data. The regression coefficient (R2) and root mean square error (RMSE) were used as evaluation metrics to assess the accuracy of the simulation, which is calculated as follows:
R 2 = 1 i = 1 N Y i Y i 2 i = 1 N Y i Y 2
R M S E = 1 N i = 1 N Y i Y i 2
where N denotes the number of data samples, Y i represents the observed surface settlement values, Y i represents the predicted surface settlement values, and Y is the mean of the observed surface settlement values.
Figure 6 displays the contour map of surface settlement, while Table 8 presents the assessment results of the simulation compared to the actual. In Table 8, the R2 between the simulation data and the actual monitoring data is 0.911, and the RMSE is 0.235. This high consistency between the simulation results and the monitoring settlement data not only validates the reliability of numerical simulation but also confirms the effectiveness of grey relational analysis in capturing the critical construction parameters.
Further analysis of Figure 6 yields:
  • Lateral settlement varies with strata zones, reaching up to 24.2 mm with a broad trough (approximately 20 m wide) in weak strata (140–170 rings). In more stable strata (170–180 rings), maximum surface settlement decreases to 2.1 mm with a narrower trough (about 12 m).
  • Longitudinal settlement peaks at 24 mm in the 140–160 rings, with local minima at 140–145 and 152–155 rings. In stable strata, maximum surface settlement rapidly decreases from 24.2 mm to 10.3 mm.

4.5. Simulation of Influence Patterns

To further investigate the influence patterns of critical construction parameters on surface settlement, systematic numerical simulations were conducted following the simulation scheme described in Section 4.3. Figure 7 and Figure 8 present the results of maximum lateral and longitudinal settlements under different schemes of varying parameters.
From Figure 7 and Figure 8, it can be seen that:
  • Elastic modulus significantly impacts surface settlement, showing a nonlinear relationship (Figure 7a). Strata with smaller elastic modulus demonstrate increased sensitivity to shield tunnelling, resulting in larger settlement and wider troughs. In Figure 8a, reducing the elastic modulus from 150% to 50% leads to a 248% increase in maximum settlement (from 6.2 mm to 21.6 mm) in intervals with weak strata. Meanwhile, the impact on stable strata intervals is less than 50%. The actual settlement data across fractured zones (Table 3, Data ID 73–78) also coincide with this observation.
  • Earth pressure’s influence on surface settlement parallels elastic modulus trends (Figure 7b and Figure 8b). In stable strata, it has a moderate impact, but in poorer geological conditions, reducing earth pressure from 138 kPa to 20 kPa increases maximum settlement from 8.6 mm to 18.3 mm, consistent with previous findings [41,42].
  • As indicated in previous studies [18,36], the influence of grouting pressure on surface settlement is evident. Moreover, it greatly influences the settlement trough width in weak strata (Figure 7c). Increasing grouting pressure from 134 kPa to 334 kPa decreases the settlement trough from 10 m to 6 m. Somewhat differently from a previous study [40], no surface uplift occurs, likely due to reasonable grouting pressure within the soil’s bearing capacity.
  • Thrust magnitude directly affects the final surface settlement [33] (Figure 8d). Excessive thrust that is beyond the strata’s resistance during tunnelling leads to a significant surface settlement with surface uplift and subsequent sinking after shield tunnelling. Similarly, insufficient thrust results in incomplete compensation of strata displacement, leading to an increase in surface settlement.

5. Discussion

This study investigated the surface settlement in the construction section of Qingdao Metro Line 6 between Haigang Road Station and Chaoyang Road Station through a hybrid numerical approach comprised of grey relational analysis and numerical simulation. The research findings provide scientific evidence and engineering guidance for mitigating surface settlement during shield tunnel construction. Based on the results of grey relational analysis and numerical simulations, the following recommendations are proposed for the right-line tunnelling and future urban subway project.
As illustrated in Figure 7 and Figure 8, it is important to note that the elastic modulus of the stratum has a significant influence on the surface settlement. Given its inherent property as a critical stratum parameter, precise characterization of the elastic modulus can provide a direct link to the potential surface settlement and offer a scientific basis for developing effective mitigation strategies. Geotechnical investigations and predictive modelling before shield tunnelling are required. In addition, it is advisable to conduct geological predictions in advance and promptly adjust construction parameters based on real-time stratum conditions during tunnel evacuation. In regions with small elastic moduli, suitable reinforcement measures like pre-excavation grouting or enhanced ground stabilization are also recommended.
Similarly, the influence of earth pressure on the surface settlement should also be considered to mitigate the settlement fluctuations during shield tunnelling. Careful adjustments in earth pressure in sections with weak geological conditions are vital to prevent excessive surface settlement, and utilizing air pressure can enhance tunnelling efficiency in strata with stable geological conditions. In light of the notable influence of grouting pressure and thrust magnitude on surface settlement, it becomes crucial to determine suitable magnitude levels of these parameters through numerical modelling and field validation. According to the simulation results, maintaining the grouting pressure between 250 and 300 kPa and keeping the thrust force matched with the initial earth pressure can effectively mitigate surface settlement. Findings (Table 5) highlight that maintaining stable and reasonable tunnelling speed throughout the process can mitigate its impact on surface settlement. For the right-line tunnelling, maintaining a tunnelling speed similar to the left line can help control the surface settlement.
The research findings offer scientific insights and practical guidance in mitigating the surface settlement during shield tunnel construction. However, there are also some limitations and constraints in this study. One constraint lies in the assumption of the initial stress field, which considers only the self-weight of the surrounding rock and overlooked parameters like rock displacement and external influences such as tectonic shifts, geological history, and seismic activity. In addition, there are also limitations when employing the Mohr–Coulomb strength criterion for the deformation of rock and soil. It cannot fully capture the complexities of various geological conditions and mechanical phenomena, especially the deformation and fracture behavior of soft rocks. Furthermore, due to limitations in computational resources, simulation time, and data availability, consideration of influential parameters like the groundwater level shifts and soil permeability’s impact on surface settlement was omitted in this study. There are also certain insufficiencies in model dimensions created through numerical simulation.
To address these limitations, future study will concentrate on three main areas. Firstly, more accurate data of material strength will be acquired through model testing and field monitoring. Consequently, a more realistic initial in situ stress field model will be constructed based on these measuring data, thus enhancing simulation accuracy. Secondly, advanced mechanical models will be developed to comprehensively represent the variations in material strength, especially in describing the deformation and fracture behavior of soft rocks. Finally, efforts will be made to extend the numerical model size to better fit the engineering context. Additional factors will be taken into account, including groundwater levels, infiltration coefficients, fracture development, and the stratigraphic properties of other intervals. By incorporating diverse engineering scenarios and geological contexts, more precise engineering recommendations and effective control strategies for surface settlement will be provided.

6. Conclusions

This study investigated surface settlement in the urban subway construction section of Qingdao Metro Line 6 between Haigang Road Station and Chaoyang Road Station. By utilizing a hybrid numerical approach, the main conclusions are as follows.
  • Earth pressure emerges as the most influential tunnelling parameter, while tunnelling speed has minimal impact. Grouting pressure has a greater effect on surface settlement than grouting volume, and strata elastic modulus significantly influences settlement.
  • Numerical simulations reveal that strata with lower elastic moduli are more sensitive to shield tunnelling, resulting in larger settlement. Earth pressure and grouting pressure show similar trends, with pronounced effects in weak strata.
  • For stable strata, it is feasible to consider a judicious increment in tunnelling speed and employ compressed air pressure for excavation. Keeping grouting pressure between 250 and 300 kPa and adjusting thrust with initial earth pressure can mitigate settlement.

Author Contributions

Conceptualization, M.L. and X.W.; methodology, M.L. and D.W.; field testing, Z.L.; data curation, M.L. and Z.L.; data analysis, M.L.; writing—original draft preparation, M.L.; writing—review and editing, D.W. and X.W.; supervision, X.W. and D.W.; project administration, X.W. and D.W.; funding acquisition, X.W. and D.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the National Natural Science Foundation of China (grant 52204115), the Foundation of Research Institute for Deep Underground Science and Engineering (grant XD2021022), and the Shield/TBM Construction Risk Consultation Project in Qingdao Metro Line 6.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data included in this study are available upon request by contact with the corresponding author.

Acknowledgments

Thanks are due to Guangzhao Zhang (Qingdao University of Technology, Qingdao) for assistance with the numerical simulation, and Zelin Lu (Qingdao University of Technology, Qingdao) for valuable advice.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ji, M.; Wang, X.; Luo, M.; Wang, D.; Teng, H.; Du, M. Stability Analysis of Tunnel Surrounding Rock When TBM Passes through Fracture Zones with Different Deterioration Levels and Dip Angles. Sustainability 2023, 15, 5243. [Google Scholar] [CrossRef]
  2. Tian, Y.; Motalleb Qaytmas, A.; Lu, D.; Du, X. Stress Path of the Surrounding Soil during Tunnel Excavation: An Experimental Study. Transp. Geotech. 2023, 38, 100917. [Google Scholar] [CrossRef]
  3. Liu, C.; Yang, S.; Liu, W.; Wang, Z.; Jiang, Y.; Yang, Z.; Jiang, H. Three-Dimensional Numerical Simulation of Soil Deformation during Shield Tunnel Construction. Math. Probl. Eng. 2022, 2022, 5029165. [Google Scholar] [CrossRef]
  4. Peng, S.; Huang, W.; Luo, G.; Cao, H.; Pan, H.; Mo, N. Failure Mechanisms of Ground Collapse Caused by Shield Tunnelling in Water-Rich Composite Sandy Stratum: A Case Study. Eng. Fail. Anal. 2023, 146, 107100. [Google Scholar] [CrossRef]
  5. Seol, H.; Won, D.; Jang, J.; Kim, K.Y.; Yun, T.S. Ground Collapse in EPB Shield TBM Site: A Case Study of Railway Tunnels in the Deltaic Region near Nak-Dong River in Korea. Tunn. Undergr. Space Technol. 2022, 120, 104274. [Google Scholar] [CrossRef]
  6. Yao, Q.; Di, H.; Ji, C.; Zhou, S. Ground Collapse Caused by Shield Tunneling in Sandy Cobble Stratum and Its Control Measures. Bull. Eng. Geol. Environ. 2020, 79, 5599–5614. [Google Scholar] [CrossRef]
  7. Bayati, M.; Hamidi, J.K. A Case Study on TBM Tunnelling in Fault Zones and Lessons Learned from Ground Improvement. Tunn. Undergr. Space Technol. 2017, 63, 162–170. [Google Scholar] [CrossRef]
  8. Peck, B.B. Deep excavation and tunnelling in soft ground, State of the art volume. In Proceedings of the 7th ICSMFE, Mexico City, Mexico, 25–29 August 1969; pp. 225–290. Available online: https://www.issmge.org/publications/publication/deep-excavations-and-tunneling-in-soft-ground (accessed on 27 August 2023).
  9. Ren, Y.; Zhang, C.; Zhu, M.; Chen, R.; Wang, J. Significance and Formulation of Ground Loss in Tunneling-Induced Settlement Prediction: A Data-Driven Study. Acta Geotech. 2023, 18, 4941–4956. [Google Scholar] [CrossRef]
  10. Jin, H.; Yuan, D.; Jin, D.; Wu, J.; Wang, X.; Han, B.; Mao, J. Shield Kinematics and Its Influence on Ground Settlement in Ultra-Soft Soil: A Case Study in Suzhou. Can. Geotech. J. 2022, 59, 1887–1900. [Google Scholar] [CrossRef]
  11. Lou, P.; Li, Y.; Tang, X.; Lu, S.; Xiao, H.; Zhang, Z. Influence of Double-Line Large-Slope Shield Tunneling on Settlement of Ground Surface and Mechanical Properties of Surrounding Rock and Segment. Alex. Eng. J. 2023, 63, 645–659. [Google Scholar] [CrossRef]
  12. Sun, F.; Jin, Z.; Wang, C.; Gou, C.; Li, X.; Liu, C.; Yu, Z. Case Study on Tunnel Settlement Calculations during Construction Considering Shield Disturbance. KSCE J. Civ. Eng. 2023, 27, 2202–2216. [Google Scholar] [CrossRef]
  13. Ahn, C.-Y.; Park, D.; Moon, S.-W. Analysis of Surface Settlement Troughs Induced by Twin Shield Tunnels in Soil: A Case Study. Geomech. Eng. 2022, 30, 325–336. [Google Scholar] [CrossRef]
  14. Wang, F.; Du, X.; Li, P. Predictions of Ground Surface Settlement for Shield Tunnels in Sandy Cobble Stratum Based on Stochastic Medium Theory and Empirical Formulas. Undergr. Space 2023, 11, 189–203. [Google Scholar] [CrossRef]
  15. Sohaei, H.; Hajihassani, M.; Namazi, E.; Marto, A. Experimental Study of Surface Failure Induced by Tunnel Construction in Sand. Eng. Fail. Anal. 2020, 118, 104897. [Google Scholar] [CrossRef]
  16. Bel, J.; Branque, D.; Camus, T. Physical Modelling of EPB TBM in Dry Sand and Greenfield Conditions. Eur. J. Environ. Civ. Eng. 2023, 27, 3236–3259. [Google Scholar] [CrossRef]
  17. Berthoz, N.; Branque, D.; Wong, H.; Subrin, D. TBM Soft Ground Interaction: Experimental Study on a 1 g Reduced-Scale EPBS Model. Tunn. Undergr. Space Technol. 2018, 72, 189–209. [Google Scholar] [CrossRef]
  18. Kannangara, K.K.P.M.; Ding, Z.; Zhou, W.-H. Surface Settlements Induced by Twin Tunneling in Silty Sand. Undergr. Space 2022, 7, 58–75. [Google Scholar] [CrossRef]
  19. Ter-Martirosyan, A.Z.; Cherkesov, R.H.; Isaev, I.O.; Shishkina, V.V. Surface Settlement during Tunneling: Field Observation Analysis. Appl. Sci. 2022, 12, 9963. [Google Scholar] [CrossRef]
  20. Kim, D.; Pham, K.; Park, S.; Oh, J.Y.; Choi, H. Determination of Effective Parameters on Surface Settlement during Shield TBM. Geomech. Eng. 2020, 21, 153–164. [Google Scholar] [CrossRef]
  21. Wang, D.; He, M.; Tao, Z.; Guo, A.; Wang, X. Deformation-Softening Behaviors of High-Strength and High-Toughness Steels Used for Rock Bolts. J. Rock Mech. Geotech. Eng. 2022, 14, 1872–1884. [Google Scholar] [CrossRef]
  22. Wang, D.; He, M.; Jia, L.; Sun, X.; Xia, M.; Wang, X. Energy Absorption Characteristics of Novel High-Strength and High-Toughness Steels Used for Rock Support. J. Rock Mech. Geotech. Eng. 2023, 15, 1441–1456. [Google Scholar] [CrossRef]
  23. Chen, R.; Meng, F.; Ye, Y.; Liu, Y. Numerical Simulation of the Uplift Behavior of Shield Tunnel during Construction Stage. Soils Found. 2018, 58, 370–381. [Google Scholar] [CrossRef]
  24. Jia, J.; Gao, R.; Wang, D.; Li, J.; Song, Z.; Tan, J. Settlement Behaviours and Control Measures of Twin-Tube Curved Buildings-Crossing Shield Tunnel. Struct. Eng. Mech. 2022, 84, 699–706. [Google Scholar] [CrossRef]
  25. Jallow, A.; Ou, C.-Y.; Lim, A. Three-Dimensional Numerical Study of Long-Term Settlement Induced in Shield Tunneling. Tunn. Undergr. Space Technol. 2019, 88, 221–236. [Google Scholar] [CrossRef]
  26. Eskandari, F.; Goharrizi, K.G.; Hooti, A. The Impact of EPB Pressure on Surface Settlement and Face Displacement in Intersection of Triple Tunnels at Mashhad Metro. Geomech. Eng. 2018, 15, 769–774. [Google Scholar] [CrossRef]
  27. Miliziano, S.; De Lillis, A. Predicted and Observed Settlements Induced by the Mechanized Tunnel Excavation of Metro Line C near S. Giovanni Station in Rome. Tunn. Undergr. Space Technol. 2019, 86, 236–246. [Google Scholar] [CrossRef]
  28. Shahin, H.M.; Nakai, T.; Okuno, T. Numerical Study on 3D Effect and Practical Design in Shield Tunneling. Undergr. Space 2019, 4, 201–209. [Google Scholar] [CrossRef]
  29. Bahri, M.; Mascort-Albea, E.J.; Romero-Hernández, R.; Koopialipoor, M.; Soriano-Cuesta, C.; Jaramillo-Morilla, A. Numerical Model Validation for Detection of Surface Displacements over Twin Tunnels from Metro Line 1 in the Historical Area of Seville (Spain). Symmetry 2022, 14, 1263. [Google Scholar] [CrossRef]
  30. Anato, N.J.; Chen, J.; Tang, A.; Assogba, O.C. Numerical Investigation of Ground Settlements Induced by the Construction of Nanjing WeiSanLu Tunnel and Parametric Analysis. Arab. J. Sci. Eng. 2021, 46, 11223–11239. [Google Scholar] [CrossRef]
  31. Wang, Z.; Li, G.; Wang, A.; Pan, K. Numerical Simulation Study of Stratum Subsidence Induced by Sand Leakage in Tunnel Lining Based on Particle Flow Software. Geotech. Geol. Eng. 2020, 38, 3955–3965. [Google Scholar] [CrossRef]
  32. Alzabeebee, S.; Keawsawasvong, S. Numerical Assessment of Microtunnelling Induced Pavement Settlement. Geotech. Geol. Eng. 2023, 41, 2173–2184. [Google Scholar] [CrossRef]
  33. Yuan, C.; Zhang, M.; Ji, S.; Li, J.; Jin, L. Analysis of Factors Influencing Surface Settlement during Shield Construction of a Double-Line Tunnel in a Mudstone Area. Sci. Rep. 2022, 12, 22606. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, Y.; Jiang, W.; Wang, M.; Li, Y. Risk Assessment and Implementation of Deformation Disaster for Operation Tunnel Based on Entropy Weight-Grey Relational Analysis. Geomat. Nat. Hazards Risk 2022, 13, 1831–1848. [Google Scholar] [CrossRef]
  35. Ma, J.; He, S.; Cui, G.; He, J.; Liu, X. Construction Stability and Reinforcement Technology for the Super-Large Rectangular Pipe-Jacking Tunnel Passing beneath the Operational High-Speed Railway in Composite Stratum. Geomat. Nat. Hazards Risk 2023, 14, 2208720. [Google Scholar] [CrossRef]
  36. Broere, W.; Festa, D. Correlation between the Kinematics of a Tunnel Boring Machine and the Observed Soil Displacements. Tunn. Undergr. Space Technol. 2017, 70, 125–147. [Google Scholar] [CrossRef]
  37. Goh, A.T.C.; Zhang, W.; Zhang, Y.; Xiao, Y.; Xiang, Y. Determination of Earth Pressure Balance Tunnel-Related Maximum Surface Settlement: A Multivariate Adaptive Regression Splines Approach. Bull. Eng. Geol. Env. 2018, 77, 489–500. [Google Scholar] [CrossRef]
  38. Shi, Y.-J.; Xiao, X.; Li, M.-G. Long-Term Longitudinal Deformation Characteristics of Metro Lines in Soft Soil Area. J. Aerosp. Eng. 2018, 31, 04018080. [Google Scholar] [CrossRef]
  39. Pan, Y.; Ou, S.; Zhang, L.; Zhang, W.; Wu, X.; Li, H. Modeling Risks in Dependent Systems: A Copula-Bayesian Approach. Reliab. Eng. Syst. Saf. 2019, 188, 416–431. [Google Scholar] [CrossRef]
  40. Moghtader, T.; Sharafati, A.; Naderpour, H.; Gharouni Nik, M. Estimating Maximum Surface Settlement Caused by EPB Shield Tunneling Utilizing an Intelligent Approach. Buildings 2023, 13, 1051. [Google Scholar] [CrossRef]
  41. Vinoth, M.; Aswathy, M.S. 3D Evolution of Soil Arching during Shield Tunnelling in Silty and Sandy Soils: A Comparative Study. Int. J. Numer. Anal. Methods Geomech. 2023, 47, 299–322. [Google Scholar] [CrossRef]
  42. Wang, J.; Feng, K.; Wang, Y.; Lin, G.; He, C. Soil Disturbance Induced by EPB Shield Tunnelling in Multilayered Ground with Soft Sand Lying on Hard Rock: A Model Test and DEM Study. Tunn. Undergr. Space Technol. 2022, 130, 104738. [Google Scholar] [CrossRef]
Figure 1. Geological profile of the left-line construction section between Haigang Road Station and Chaoyang Road Station.
Figure 1. Geological profile of the left-line construction section between Haigang Road Station and Chaoyang Road Station.
Sustainability 15 14315 g001
Figure 2. Layout plan of monitoring sites at Haigang Road Station.
Figure 2. Layout plan of monitoring sites at Haigang Road Station.
Sustainability 15 14315 g002
Figure 3. Aerial photo of the proposed tunnel construction section between Haigang Road Station and Chaoyang Road Station.
Figure 3. Aerial photo of the proposed tunnel construction section between Haigang Road Station and Chaoyang Road Station.
Sustainability 15 14315 g003
Figure 4. Calculation results of grey relational coefficients. A relational coefficient of 0.7 is placed in each subfigure to represent the correlation threshold between the influential parameters and the surface settlement: (a) relational coefficients of stratum parameters; (b) relational coefficients of tunnelling parameters; (c) relational coefficients of grouting parameters.
Figure 4. Calculation results of grey relational coefficients. A relational coefficient of 0.7 is placed in each subfigure to represent the correlation threshold between the influential parameters and the surface settlement: (a) relational coefficients of stratum parameters; (b) relational coefficients of tunnelling parameters; (c) relational coefficients of grouting parameters.
Sustainability 15 14315 g004aSustainability 15 14315 g004b
Figure 5. Three-dimensional numerical model.
Figure 5. Three-dimensional numerical model.
Sustainability 15 14315 g005
Figure 6. Contour map of surface settlement.
Figure 6. Contour map of surface settlement.
Sustainability 15 14315 g006
Figure 7. Maximum lateral surface settlement under different variations of critical construction parameters: (a) maximum lateral surface settlement under different elastic moduli; (b) maximum lateral surface settlement under different earth pressure; (c) maximum lateral surface settlement under different grouting pressure; (d) maximum lateral surface settlement under different thrust.
Figure 7. Maximum lateral surface settlement under different variations of critical construction parameters: (a) maximum lateral surface settlement under different elastic moduli; (b) maximum lateral surface settlement under different earth pressure; (c) maximum lateral surface settlement under different grouting pressure; (d) maximum lateral surface settlement under different thrust.
Sustainability 15 14315 g007aSustainability 15 14315 g007b
Figure 8. Maximum longitudinal surface settlement under different variations of critical construction parameters: (a) maximum longitudinal surface settlement under different elastic modulis; (b) maximum longitudinal surface settlement under different earth pressure; (c) maximum longitudinal surface settlement under different grouting pressure; (d) maximum longitudinal surface settlement under different thrust.
Figure 8. Maximum longitudinal surface settlement under different variations of critical construction parameters: (a) maximum longitudinal surface settlement under different elastic modulis; (b) maximum longitudinal surface settlement under different earth pressure; (c) maximum longitudinal surface settlement under different grouting pressure; (d) maximum longitudinal surface settlement under different thrust.
Sustainability 15 14315 g008aSustainability 15 14315 g008b
Table 1. Main parameters of the shield machine.
Table 1. Main parameters of the shield machine.
ParameterCR Equip. Machine 179
Tunnel boring typeEarth pressure balance
Excavation diameter (mm)6303
Shield diameter (mm)6250
Cutterhead rotation speed (rpm)0~3.7
Maximum thrust (kN)39,000
Rated torque (kN/m)6000
Unblocking torque (kN/m)7200
Horizontal turning capability (m)250
Maximum climbing capability (‰)±50
Cutterhead opening ratio (%)34
Table 2. Stratigraphic physical and mechanical parameters of the construction section between Haigang Road Station and Chaoyang Road Station.
Table 2. Stratigraphic physical and mechanical parameters of the construction section between Haigang Road Station and Chaoyang Road Station.
StratumDensity (kg/m3)Cohesion (kPa)Internal Friction (°)Elastic Modulus 1 (MPa)Poisson’s Ratio
Plain fill soil165081480.3
Powdery clay19801816160.4
Clayey gravelly sand20501428120.3
Sandy tuff (fractured)24502533400.3
Blocky granite (fractured)250048351000.3
Slightly weathered tuff26107104512500.2
Moderately weathered tuff2600400405500.3
Strongly weathered tuff23002230200.3
Slightly weathered granite porphyry26805304510800.2
Moderately weathered granite porphyry2600300404000.3
Strongly weathered granite porphyry23001530180.3
1 The elastic modulus in the table is the tangential modulus.
Table 3. Influential construction parameters and corresponding surface settlement.
Table 3. Influential construction parameters and corresponding surface settlement.
Data IDSettlement (cm)Tunnelling ParametersGrouting ParametersStratum Parameters
Thrust (kN)Earth
Pressure (kPa)
Tunnelling Speed (mm/min)Grouting Pressure (kPa)Grout
Volume (m3)
Cohesion (kPa)Internal Friction (°)Poisson’s RatioElastic Modulus (MPa)
0−3.1010,42483152723.423739.70.23471
1−2.34981794332443.118037.60.22477
2−4.29765669262073.123738.90.21518
3−1.5411,460121132985.022339.90.22542
4−0.6110,138113152444.816438.20.22489
5−6.2710,43818122603.821337.60.21448
6−2.44985590162233.223042.90.20533
7−3.12944081123003.823742.90.20585
8−5.70688617221945.317240.90.21391
9−7.09744137331855.714538.70.21323
10−9.28779121202035.616038.40.20238
11−9.90741921132035.319135.40.20158
12−8.79796569372145.916139.50.21150
13−8.71948982152205.617540.10.21207
14−7.53704573281885.97735.90.26220
15−7.56764582292705.813737.90.24298
16−4.50892176273345.314240.10.23383
17−3.92780962272235.022542.90.21499
18−3.36937954391905.219840.60.22475
19−1.619723108142956.119940.40.22516
20−0.9210,344118322776.123742.90.22524
21−1.6110,07466242325.417740.40.23435
22−1.3410,217112272324.415539.30.24408
230.0310,678126262924.814740.10.24429
24−0.2110,735125282835.610735.70.28586
25−2.09996164442646.215939.60.27447
26−1.689572104252744.14032.60.28495
27−2.429374125342746.24033.30.28461
28−1.1311,109105262476.24031.00.28520
29−1.259547109342374.44029.90.28510
30−2.308677142362333.64027.00.28454
31−1.0910,429145302683.64029.50.28515
32−1.1110,260126412433.24031.20.28518
33−3.119504126381943.54033.50.27427
34−6.2310,70374342213.14025.90.28160
35−7.14704164332003.24025.70.28147
36−7.76757455371813.14023.50.26240
37−10.39756326291354.110223.30.2391
38−10.37717025281503.314628.70.24168
39−3.089289107272283.916328.10.26470
40−2.829150101122373.723725.40.28654
Table 4. Calculation results of grey relational coefficients.
Table 4. Calculation results of grey relational coefficients.
Data IDRelational Coefficients of
Tunnelling Parameters
Relational Coefficients of
Grouting Parameters
Relational Coefficients of
Stratum Parameters
ThrustEarth
Pressure
Tunnelling SpeedGrouting
Pressure
Grouting
Volume
Cohesion Internal Friction Poisson’s RatioElastic
Modulus
00.870.720.440.980.440.620.780.580.95
10.790.740.970.680.400.730.920.590.85
20.540.730.880.690.500.540.700.710.74
30.760.920.360.940.800.760.990.560.91
40.680.720.340.560.660.630.730.520.68
50.560.550.650.670.960.440.590.950.67
60.810.710.400.600.500.670.670.510.96
70.780.710.410.780.640.620.620.450.73
80.520.520.760.761.000.470.520.810.85
90.710.750.430.880.830.420.510.970.83
100.840.870.610.670.740.480.420.800.76
110.870.970.560.620.670.400.460.910.87
120.860.660.380.670.590.500.420.760.91
130.540.580.610.650.650.480.410.840.91
140.670.750.470.980.640.500.570.570.91
150.820.670.460.540.660.410.510.700.83
160.800.830.780.530.790.580.620.990.91
170.540.640.920.730.910.560.560.680.82
180.790.560.600.550.910.600.700.720.98
190.690.790.370.930.830.780.940.560.85
200.760.810.600.710.900.840.840.500.78
210.770.510.480.580.950.950.940.670.68
220.780.790.540.560.630.730.900.720.61
230.740.760.400.690.610.480.770.590.55
240.780.790.460.680.830.330.580.930.83
250.800.530.700.760.700.960.940.870.75
260.660.760.520.780.600.370.570.760.80
270.690.860.890.880.670.390.660.750.82
280.930.710.490.600.810.350.490.920.79
290.620.760.780.570.630.360.470.800.78
300.560.710.830.630.540.380.450.820.79
310.800.820.570.680.420.350.460.820.78
320.760.920.840.580.380.350.500.840.79
330.790.760.630.550.590.410.730.790.83
340.530.910.460.930.640.740.640.500.64
350.640.890.440.970.870.930.720.460.70
360.830.910.430.960.720.610.670.500.98
370.770.870.411.001.000.551.000.521.00
380.890.890.420.871.000.460.640.460.78
390.731.000.870.680.660.680.520.890.95
400.680.870.400.690.600.640.440.730.64
Table 5. Grey relational degrees of influential construction parameters.
Table 5. Grey relational degrees of influential construction parameters.
ParameterRelational DegreeRanking
Tunnelling ParametersEarth Pressure (kPa)0.7622
Tunnelling Speed (mm/min)0.5758
Thrust (kN)0.7303
Grouting ParametersGrouting Pressure (kPa)0.7264
Grout Volume (m3)0.7056
Stratum ParametersCohesion (kPa)0.5629
Internal Friction (°)0.6487
Elastic Modulus (MPa)0.8091
Poisson’s Ratio0.7085
Table 6. Simulation scheme sets for different levels of construction parameters.
Table 6. Simulation scheme sets for different levels of construction parameters.
SchemesThrust
(kN)
Earth Pressure
(kPa)
Elastic Modulus
(%)
Grouting Pressure
(kPa)
1688679100234
2803079100234
3917379100234
410,31779100234
511,46079100234
6917320100234
7917350100234
89173109100234
99173138100234
1091737950234
1191737975234
12917379125234
13917379150234
14917379100134
15917379100184
16917379100284
17917379100334
Table 7. Material parameters of the tunnel and shield zones.
Table 7. Material parameters of the tunnel and shield zones.
Zone NameDensity (N/m3)Elastic Modulus (MPa)Poisson’s Ratio
Segment rings325030000.2
Shield shell9500200,0000.3
Grouting layer (pre-solidification)21002400.3
Grouting layer (post-solidification)21007200.3
Table 8. Validation results of the simulation compared to the actual.
Table 8. Validation results of the simulation compared to the actual.
Data NameActual DataSimulation DataR2RMSE
ZDBC65−2.42−2.370.9110.235
(Data ID 27)
ZDBC64−1.68−2.04
(Data ID 26)
ZDBC63−2.09−2.23
(Data ID 25)
ZDBC62−0.21−0.56
(Data ID 24)
ZDBC610.03−0.36
(Data ID 23)
ZDBC60−1.34−1.09
(Data ID 22)
DBC10-01−0.79−0.92
DBC10-01−1.86−2.01
DBC10-01−1.95−2.02
DBC10-04−0.88−1.01
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luo, M.; Wang, D.; Wang, X.; Lu, Z. Analysis of Surface Settlement Induced by Shield Tunnelling: Grey Relational Analysis and Numerical Simulation Study on Critical Construction Parameters. Sustainability 2023, 15, 14315. https://doi.org/10.3390/su151914315

AMA Style

Luo M, Wang D, Wang X, Lu Z. Analysis of Surface Settlement Induced by Shield Tunnelling: Grey Relational Analysis and Numerical Simulation Study on Critical Construction Parameters. Sustainability. 2023; 15(19):14315. https://doi.org/10.3390/su151914315

Chicago/Turabian Style

Luo, Minhe, Ding Wang, Xuchun Wang, and Zelin Lu. 2023. "Analysis of Surface Settlement Induced by Shield Tunnelling: Grey Relational Analysis and Numerical Simulation Study on Critical Construction Parameters" Sustainability 15, no. 19: 14315. https://doi.org/10.3390/su151914315

APA Style

Luo, M., Wang, D., Wang, X., & Lu, Z. (2023). Analysis of Surface Settlement Induced by Shield Tunnelling: Grey Relational Analysis and Numerical Simulation Study on Critical Construction Parameters. Sustainability, 15(19), 14315. https://doi.org/10.3390/su151914315

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop