Next Article in Journal
Corridor Detection from Large GPS Trajectories Datasets
Previous Article in Journal
The Effectiveness of Nafion-Coated Stainless Steel Surfaces for Inhibiting Bacillus Subtilis Biofilm Formation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Variability of Grouting Layer of Shield Tunnel and Its Effect on Ground Settlement

1
School of Civil Engineering, Central South University, Changsha 410083, China
2
School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China
3
Department of Civil & Environmental Engineering, National University of Singapore, Singapore 119077, Singapore
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(14), 5002; https://doi.org/10.3390/app10145002
Submission received: 23 June 2020 / Revised: 14 July 2020 / Accepted: 15 July 2020 / Published: 21 July 2020
(This article belongs to the Section Civil Engineering)

Abstract

:
This study aims to investigate the effect of the spatial variability of grouting-layer thickness on ground-surface settlement caused by shield tunneling and to provide a rational prediction method. The spatial characteristics of grouting layers were obtained based on statistical analysis. The random finite element method was used to study the effect of spatial variability of different parameters on ground-surface settlement. Simulation results indicate that the spatial variability of the grouting layer has a negative impact on ground settlement. The surface settlement will be underestimated without considering the spatial characteristics of the grouting layer. Thus, a reliable prediction approach of the maximum ground settlement was proposed to control the construction quality.

1. Introduction

With the development of urban underground spaces, shield-tunneling methods have been widely used in urban subway engineering. However, the excavation of shield tunnels will inevitably cause ground settlement. This sometimes even leads to ground-collapse accidents. Therefore, the prediction of ground settlement caused by shield tunneling plays an important role in controlling construction quality.
During the past years, tunneling-induced ground settlement problems have attracted many researchers. The methods used to analyze this problem can be roughly classified into four categories: empirical formula method, laboratory tests, numeric simulation method and analytical method. For the (a) empirical formula method, Peck [1] derived the empirical formula for predicting the transverse settlement troughs of ground based on a large amount of field data; (b) For laboratory tests, in order to study the effect of different factors on the ground settlement, a series of small-scale models have been tested in the laboratory [2,3,4,5,6,7,8]; (c) for numeric simulation, it is common to simulate ground movement by using the finite element method (FEM) [9,10,11,12,13,14]. In practice, the mechanical properties of the soil have great spatial variability. In order to study the effect of spatial variability on the tunneling-induced ground movement, some numeric investigations with random finite element analysis (RFEA) have also been performed [15,16], providing more reliable references for tunnel construction. For the (d) analytical method, Rowe et al. [17] defined the concept of excavation gap parameter g and proposed a theoretical prediction model for predicting the ground settlement caused by shield tunneling in soft ground. On this basis, Loganathan et al. [18] further studied the ground loss rate caused by tunnel excavation and developed an analytical solution which can consider the oval-shaped deformation of the tunnel. Most of the above studies assumed that the settlement induced by tunnel boring machine (TBM) is symmetrical or in regular shape. Indeed, fine workmanship usually produces a symmetrical and normal settlement distribution.
Due to the gap between the shield tail and segmental lining, the ground settlement mainly occurs at the shield tail passing stage [19]. Backfill-grouting is critical for filling the gap and reducing ground settlement. However, due to the complexity and uncertainty of geological and environmental conditions, the shield tail gap cannot be completely refilled. Ground settlement after the shield-tail passing can be found even after sufficient grouting material is injected. This is considered to be related to the uneven distribution of the grouting layer and that some injected grouting materials are wasted [20]. Nevertheless, there is still a lack of research on the spatial variability of the grouting layer and its effect on ground settlement.
In this study, the spatial variability of the thickness of the grouting layer in shield tunnel is summarized, first based on collected test data. Based on the Mohr–Coulomb theory, a random finite element model of tunnel excavation is established in a spatial heterogeneous stratum. Two dimensional random variables were used to describe soil parameters and grouting-layer thickness. A Monte Carlo simulation was carried out to analyze the influences of spatial variability of soil parameters and grouting-layer thickness on ground settlement. Finally—to provide a more reliable prediction of the ground settlement caused by shield tunneling—a practical and realistic approach to estimate the tunneling-induced ground settlement is proposed based on the probabilistic method.

2. Spatial Variability of Grouting Layers

The shield tunneling method allows excavation while preventing tunnel collapse by assembling steel or concrete segmental blocks in a ring conformation to create a tunnel structure. In order to permit the advance of the shield machine and assemble the segments under the protection of the shield, the excavation diameter is usually larger than the external diameter of the segmental lining, thus a gap exists between the excavation boundary and segmental linings. The purpose of backfill grout injection is to immediately fill voids left by the shield machine body between the tunnel wall and the segmented blocks, preventing subsidence of surrounding ground and adjacent structures. Nevertheless, the distribution of the backfill-grouting cannot be uniformly distributed around segmental linings. As shown in Figure 1, the grouting layer in reality has spatial variability due to the geological randomness and mortar arbitrary flowing. Although the injected grouting is enough to fill the voids, a settlement can immediately occur after the passing of the shield tail. The settlement patterns and amount are considered to be closely related to the spatial variability of grouting layers [21]. The distribution of the grouting layer can be obtained by two methods, ground penetrating radar (GPR) and large visualized model test. GPR is a rapid nondestructive method that can recognize the interfaces of different objects. This method has been widely applied in the detection of tunnel defects to evaluate the quality of backfill-grouting layers [22,23]. The visualized model test is another method which can directly observe the characteristics of grouting layers, and many relevant studies can be found in the literature [24].
Based on previous studies, the statistical characteristics obtained from some monitoring results of the grouting layers are listed in Table 1. The obtained results indicate that the backfill-grouting layers have significant variability. The coefficient of variation (COV) ranges from 0.1 to 0.4. Figure 2 presents the distribution of a grouting-layer thickness; the data are collected from a monitoring section reported by Ding et al. [24]. To further test the normality hypotheses, the Anderson–Darling test within the Matlab was used. The Anderson–Darling test results indicate that the thickness of the grouting layer can be successfully fitted by the normal distribution, passing the test at 5% significance level.
It has been found that most of the ground settlement occurs at the passing of the shield tail. The distribution of the backfill-grouting layer greatly affects the convergence patterns along the excavation surface of tunnels. As shown in Figure 3, the displacement convergence patterns can be concluded as three types based on different grouting-layer distribution. In Type I, a uniform convergence pattern corresponds to a uniform distribution of the grouting layer. In Type II, the maximum convergence pattern corresponds to as that most of the grouting fluid flows to the bottom of the tunnel. The distribution of the grouting layer can be simplified as a tangent circle outside the tunnel circle. In Type III, the spatial variability convergence pattern corresponds to a grouting layer with random defects. The first two types have been widely used to estimate ground displacement and can be found in many studies [25,26]. However, the ground movement caused by the defective grouting layer is rarely discussed.

3. Modeling Aspects

3.1. FEM Model

In this study, a finite element analysis (FEA) was conducted to discuss the ground settlement induced by shield tunneling considering the backfill-grouting distribution. As shown in Figure 4, a two-dimensional nonlinear plane strain FEA established via using ABAQUS 6.14 program. The model dimension of 40 m (length) × 20 m (depth) was chosen to reduce the influence of boundary conditions on tunneling-induced stress and strain. The constitutive model of soil was based on the Mohr–Coulomb model that included soil bulk density γs, Poisson’s radio v, friction angle c, cohesion force φ, Young’s modulus E. The normal displacement of the side boundaries and the bottom boundary was fixed. The whole model adopted the four-node bilinear plane strain quadrilateral element(CPE4). The average mesh size of the stratum was 0.5 m, while the average mesh size of the tunneling area was 0.4 m. The model consisted of 3552 elements and 3675 nodes. The excavation parameters of the tunnel and soil parameters were referred to Shenzhen Metro [27]. The excavation diameter of the tunnel D was 6.0 m, the outer diameter of the lining Dl was 5.6 m, and the burial depth of the tunnel B was 9.0 m. The excavation gap g between the shield tail and lining was 200 mm. The variability of soil properties has been extensively studied in past years [28,29,30,31,32,33,34,35,36]. Phoon et al. evaluated the geotechnical variability in detail for different types of soil [34]. Based on extensive previous studies, the COVs of soil mechanical parameters c, φ and E showed significant differences, while the scale of fluctuation (SOF) of soil parameters showed a certain degree of consistency [35,36]. The reference case of this study was to simulate the clay formation. The selection of soil parameters can refer to the above research. Previous studies showed that there is low correlation between random variables of soil parameters, so it is not necessary to consider their correlation. The specific parameters of the reference case are shown in Table 2.
It is not easy to express the thickness of grouting layer in the numeric simulation software. Generally, the excavation gap parameter g is a fixed value. Therefore, this study presents an appropriate simulation method of the spatial variability of grouting layer via randomizing the convergence displacement of tunnel-excavation boundary surface. δ is equivalent to a mapping of t , in which δ exhibits the same random properties as t (the same COV). Therefore, when modeling random fields, only the randomization process of δ needs to be considered. When the grouting-layer distribution is determined, the contracted displacement around the tunnel boundary thus can be obtained,
δ = g t
where δ refers to the radial displacement at an arbitrary position around the tunnel, g refers to the voids between the excavation surface and segmental lining, t refers to the thickness of the grouting layer.
A ground contraction method was adopted in this paper to simulate the soil disturbance after the tunnel excavation. To control the convergence displacement of the excavation surface, the circumference of the tunnel was divided into 48 elements and 48 nodes. The constraint effect of tunnel supporting structure was replaced by a constraint displacement field, in which the displacement constraints were applied to 48 nodes on the circumference of the tunnel.
In ABAQUS, the value of the constraint field in the polar coordinate system was δ (negative value indicates convergence to the center of the circle, Figure 5).
In order to reflect the function of supporting structure in the tunneling process, the earth stress balance was first carried out by setting the gravity stress of the highest and lowest faces. After that, the constraint displacement field was applied at the same time as the tunnel excavation, which meant that the excavation face contract immediately after the tunnel excavation. The preliminary study showed that the mesh size and the distance from the tunnel to the side and bottom boundaries were sufficient to produce convergence under the reference case.

3.2. Random Model

It is well known that ground-surface settlement is closely related to soil properties. When the Mohr–Coulomb model was adopted, the main properties of soil are E, c, φ and v. Both practice and theory showed that the spatial variability of Young’s modulus E has a significant impact on the ground-surface settlements. When RFEM was used to simulate the spatial variation of soil shear strength, the spatial variability and interrelationship of soil properties (c and φ) also affected the surface settlement [35]. In addition, Fenton et al. pointed out that the Poisson’s ratio v of soil has smaller relative spatial variability and its importance to settlement statistics is only second–order [37]. Thus, the Poisson’s ratio in this study is fixed at 0.2 over the entire soil mass for all simulations.
The spatial variability of soil has three indices: (a) Mean of the soil parameters; (b) Coefficient of variation, which describes the variability of the soil parameters; (c) Scale of fluctuation, which describes the spatial correlation. The SOF is loosely defined as the longest distance within which a similarity between two points can be observed. Small values of SOF imply that the correlation function falls off rapidly to zero (i.e., the correlation between two points becomes rapidly smaller), which leads to rougher random fields. On the other hand, for increasing values of SOF the soil property field becomes smoother, in other words, the field shows less variability converging to a uniform field. Many scholars had done much research on the distribution range of COV and SOFs of different soils. This study intends to take the clay soil, a representative soil, as the research object. A normally distributed random field was adopted in this study to represent the soil properties and the grouting layer. The quadratic exponential autocorrelation structure of a two-dimensional random field is mathematically expressed as,
ρ = exp { [ ( x i x j s h ) 2 + ( y i y j s v ) 2 ] }
where (xi, yi) and (xj, yj) define the spatial center coordinates of two arbitrary contacts.
A standard normally distributed random field can be approximately represented by an M-term K–L expansion,
G s ( x ) i = 1 M λ i ξ i ψ i ( x )
where ξ i are independent standard normal variables, M is the truncation term of the expansion, λ i and ψ i ( x ) denote the eigenvalues and eigenfunctions of the autocorrelation function in Equation (2). x denotes a spatial vector of the random process. The convergence and the accuracy of the K–L expansion depend on the number of the term M, which relates to the domain of the random field and the scale of fluctuation. The larger the random field domain and the smaller the autocorrelation distance, the more the terms are required in Fourier series to achieve a given accuracy.
The normal random fields of the soil properties or grouting layer can be expressed as,
G s ( x ) = μ + σ G ( x )
where μ and σ are the mean and standard deviation of the grouting-layer thickness or soil properties.
The variability of soil properties has been extensively studied in past years. Phoon et al. proposed the recommended values of geotechnical variability for different types of soil [34]. Based on the suggested value of COVs, an interval of 0.1~0.3 can cover some typical soils. Some studies assumed that s h and s v are the same, which was to presume that the soil is an isotropic correlation structure. However, the variability of natural soil is usually anisotropic in reality, and the horizontal SOF s h is greater than the vertical one s v for typical soil stratum [36]. Jiang et al. summarized several previously published articles and recommended that the horizontal scale of fluctuation ranges between 10 m and 40 m, and the vertical scale of fluctuation ranges between 0.5 and 3.0 m [31]. The recommended range can be adopted in this work.
The focus of this study was to explore the influence of randomness of grouting-layer thickness on the ground-surface settlement. Under the assumption that the g is fixed, the randomness of grouting thickness would lead to the randomness of excavation boundary convergence. That is, the random distribution of radial displacement of each node reflected the spatial variability of grouting-layer thickness. The randomness of the grouting-layer thickness was simulated by controlling the displacement of each node on the excavation boundary. Regarding the randomness of grouting-layer thickness, the COV of the grouting layer ranges from 0.1 to 0.4 (see Section 2 for details). The SOF of the grouting layer is difficult to determine statistically due to the limited collected data. However, it is recognized the spatial fluctuation of the grouting layer in practice mainly relies on the grouting pipe layout. Thus, the SOF of the grouting layer is assumed as the circumferential spacing of the grouting pipes. In this study, a one-dimensional random field is employed and the SOF of the grouting layer is taken as 4.71 m. This SOF corresponds to a four grouting pipes distribution, which is the most commonly adopted in tunnel projects.
In short, parametric studies were divided into two categories. In (a) traditional soil parameters, the COV of Young’s modulus, COV of cohesion, COV of friction angle, SOF of the horizontal direction and SOF of the vertical direction, and (b) grouting-layer thickness: the mean value of grouting-layer thickness t (Corresponding to the gap parameter δ ), COV of grouting-layer thickness. In the random analysis of each parameter, Monte Carlo simulations with 300 realizations were conducted to obtain a statistically stable result, with an error of less than 0.2% (Figure 6). All random variables were considered as truncated normal distribution in order to eliminate some unrealistic data using statistical intervals of the input parameters. The details of the random parameters are shown in Table 3.

4. Results and Analysis

The maximum surface settlements in the random models were extracted and the hypothesis test was conducted. Figure 7 presents the empirical cumulative distribution function (CDF) of the maximum surface settlement. The Anderson–Darling test results indicate that the overall settlements of all the cases can be successfully fitted by the normal distribution, passing the test at the 5% significance level. After confirming that the maximum surface settlement follows the normal distribution, the probability density function (PDF) of the maximum settlements was plotted. The probability density function describes the distribution of a variable. The larger the peak point of the PDF curve, the larger the mean of the maximum surface settlement. The flatter the PDF curve, the stronger the variability of the maximum surface settlement.
Figure 8 shows the effect of the COVs of soil on the maximum surface settlement. It can be seen that the COV of the soil strength parameters, friction angle and cohesion, has little effect on the PDF of the maximum settlement. While the COV of the soil Young’s modulus has a great influence on the maximum settlement. The maximum settlement PDF becomes taller and thinner with the decrease of the COV of the soil Young’s modulus, indicating the variability of Young’s modulus is the main geological reason for the fluctuation of the surface settlement. Figure 9 presents the effect of the SOFs of soil on the maximum surface settlement. It shows that when the horizontal SOF ( s h ) increases from 10 m to 30 m, the mean value of the maximum surface settlement increases slightly. The vertical SOF ( s v ) nearly does not influence the PDF of surface settlement. This is because the vertical SOF ( s v ) is relatively smaller than the soil cover and the surface settlement is a result of stress and strain propagation in the global space.
Figure 10a shows that when the gap parameter δ increases from 10 mm to 30 mm (corresponding to the grouting-layer thickness t from 190 mm to 170 mm), the mean value and standard deviation of the maximum ground-surface settlement both increase appreciably. The apparent increase in the mean is well understood, which is mainly because the larger the radial gap, the greater the amplitude of the ground movement. It is remarkable that the standard deviation also rises observably, indicating that the increase of the gap parameter δ can also increase the variability of the ground settlement. This failure risk induced by the variability of the settlement should be fully evaluated. Figure 10b demonstrates that when COV(t) of the grouting thickness increases from 0.1 m to 0.4, the mean value changes little, whereas the standard deviation rises dramatically. It means that the random defects in the grouting layer can cause a larger maximum ground-surface settlement.

5. Engineering Applications

The construction of a tunnel needs to strictly control the ground-surface settlement, especially in urban metro projects. The above study shows that when the spatial variability of the grouting layer has a significant impact on the ground settlement, which cannot be ignored in the tunneling process. It emphasized the importance of the investigation of the grouting layer quality. The construction safety may not be guaranteed by the amount of the grouting injection volume. Additionally, the spatial variability of soil mechanics parameters is a secondary factor affecting the ground settlement. A detailed site characterization before the tunnel construction should be conducted to obtain enough information for reducing uncertainties in soil parameters.
In the design phase, a “characteristic value” is often used to account for the variability of materials [38]. Eurocode 7 recommends that the “characteristic value of a geotechnical parameter shall be selected as a cautious estimate of the value affecting the occurrence of the limit state” [38]. In the case of this study, it means selecting the 95% fractile as the “characteristic value”. Engineers can optimize the construction plan according to the predicted “characteristic ground-surface settlement”. Generally speaking, there are two methods for predicting the maximum ground settlement smax under the random field: (a) The theoretical formula derived from the classification of probability distribution; (b) The prediction curve based on the statistical value of numeric simulation. In this study, the first method is used to provide a practical and realistic approach for maximum settlement estimation.
According to the normal probability distribution function, the maximum ground settlement smax of 95% confidence degree can be obtained,
s max = μ + 1.65 σ
where μ is the mean value of a specific random case, σ is the standard deviation of a specific random case.
The definition of COV(t) is shown as below,
COV ( t ) = σ t μ t
where μt is the mean value of grouting-layer thickness, σt is the standard deviation of grouting-layer thickness.
The theoretical formula for calculating the smax can be obtained by combining Equations (5) and (6),
s max = μ t ( 1.65 × COV ( t ) + 1 )
In practice, the mean value of the grouting-layer thickness can be easily obtained using the injected grouting volume. Thus, engineers can use this approach to evaluate the maximum settlement once the grouting distribution is determined.
The characteristic value of the numeric simulations is selected as the 95% fractile of the overall simulations to draw the numeric simulation curve. Figure 11 shows the comparison between the maximum settlement obtained from the theoretical prediction approach and numeric simulations. It can be seen from this figure that the maximum settlement increase with the increasing COV of grouting-layer thickness. The slope of the predicted curve based on numeric simulation is smaller than the one based on the theoretical formula, whereas small discrepancy is noted with the increase of the COV of grouting layer. The predicted settlement is a little larger than the numeric results, indicating a conservative estimation of the surface settlement. The discrepancy is around 10% when the COV of grouting layer reaches 0.4, which can be accepted in the engineering practice.

6. Conclusions

In this study, a random finite element model of tunnel excavation was established. The effect of spatial variability of soil parameters and variability of grouting-layer thickness on the ground-surface settlement was studied, respectively. The major conclusions are summarized as follows:
The grouting layer of a shield tunnel shows a great spatial randomness along the circumferential direction. According to the statistical analysis, the thickness of the grouting layer follows a normal distribution and its coefficient of variation ranges from 0.1 to 0.4.
In terms of soil parameters, the COV of the soil strength parameters, friction Angle and cohesion, has little effect on the PDF of the maximum surface settlement, while the COV of the soil Young’s modulus shows highly sensitive to the maximum surface settlement. Both the horizontal SOF ( s h ) and the vertical SOF ( s v ) nearly have no influence on the PDF of maximum surface settlement.
The spatial variability of the grouting layer shows a significant influence on the ground settlement. The mean value and standard deviation of maximum surface settlement increase with the decrease of the grouting-layer thickness. The ground settlement will be underestimated without considering the spatial variability of the grouting layer of the shield tunnel.
In order to consider the influence of the spatial randomness of grouting layers, a practical prediction approach for ground settlement was given based on 95% fractile. It is worth noting that the value predicted by the theoretical formula is a little higher than the results from the Random Finite Element method.

Author Contributions

D.J. proposed the idea of the paper, Z.W. conducted the numerical simulations, and C.S. revised the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Project NO. 51778636). This work was also funded by the National Defense Science and Technology Innovation Special Zone.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Project NO. 51778636).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Peck, R.B. Deep excavations and tunneling in soft ground. Proc. 7th ICSMFE 1969, 2019, 225–290. [Google Scholar]
  2. Bolton, M.D.; Lu, Y.C.; Sharma, J.S. Centrifuge models of tunnel construction and compensation grouting, Int. In Proceedings of the International Symposium on Geotechnical Aspects of Underground Construction in Soft Ground, London, UK, 15 April 1996; A.A. Balkema: Brookfield, WI, USA, 1996; pp. 471–477. [Google Scholar]
  3. Nomoto, T.; Imamura, S.; Hagiwara, T.; Kusakabe, O.; Fujii, N. Shield Tunnel Construction in Centrifuge. J. Geotech. Geoenviron. Eng. 1999, 125, 289–300. [Google Scholar] [CrossRef]
  4. Chapman, D.N.; Ahn, S.K.; Hunt, D.V. Investigating ground movements caused by the construction of multiple tunnels in soft ground using laboratory model tests. Can. Geotech. J. 2007, 44, 631–643. [Google Scholar] [CrossRef]
  5. Marshall, A.M.; Mair, R.J. Tunneling beneath driven or jacked end-bearing piles in sand. Can. Geotech. J. 2011, 48, 1757–1771. [Google Scholar] [CrossRef]
  6. Mathew, G.V.; Lehane, B.M. Numerical back-analyses of greenfield settlement during tunnel boring. Can. Geotech. J. 2013, 50, 145–152. [Google Scholar] [CrossRef]
  7. Moussaei, N.; Khosravi, M.H.; Hossaini, M.F. Physical modeling of tunnel induced displacement in sandy grounds. Tunn. Undergr. Space Technol. 2019, 90, 19–27. [Google Scholar] [CrossRef]
  8. Xiang, Y.; Liu, H.; Zhang, W.; Chu, J.; Zhou, N.; Xiao, Y. Application of transparent soil model test and DEM simulation in study of tunnel failure mechanism. Tunn. Undergr. Space Technol. 2018, 74, 178–184. [Google Scholar] [CrossRef]
  9. Sharifzadeh, M.; Kolivand, F.; Ghorbani, M.; Yasrobi, S. Design of sequential excavation method for large span urban tunnels in soft ground-Niayesh tunnel. Tunn. Undergr. Space Technol. 2013, 35, 178–188. [Google Scholar] [CrossRef]
  10. Do, N.A.; Dias, D.; Oreste, P.; Djeran-Maigre, I. Three-dimensional numerical simulation of a mechanized twin tunnels in soft ground. Tunn. Undergr. Space Technol. 2014, 42, 40–51. [Google Scholar] [CrossRef]
  11. Chakeri, H.; Özçelik, Y.; Ünver, B. Investigation of ground surface settlement in twin tunnels driven with EPBM in urban area. Arab. J. Geosci. 2014, 8, 7655–7666. [Google Scholar] [CrossRef]
  12. Miro, S.; König, M.; Hartmann, D.; Schanz, T. A probabilistic analysis of subsoil parameters uncertainty impacts on tunnel-induced ground movements with a back-analysis study. Comput. Geotech. 2015, 68, 38–53. [Google Scholar] [CrossRef]
  13. Goh, A.; Xiao, Y.; Xiang, Y.; Zhang, W.; Zhang, Y. Determination of earth pressure balance tunnel-related maximum surface settlement: A multivariate adaptive regression splines approach. Bull. Int. Assoc. Eng. Geol. 2016, 77, 489–500. [Google Scholar] [CrossRef]
  14. Zheng, G.; Yang, X.; Zhou, H.; Du, Y.; Sun, J.; Yu, X. A simplified prediction method for evaluating tunnel displacement induced by laterally adjacent excavations. Comput. Geotech. 2018, 95, 119–128. [Google Scholar] [CrossRef]
  15. Xiao, L.; Huang, H.; Zhang, J. Effect of soil spatial variability on ground settlement induced by shield tunneling. Geo-Risk 2017, 2017, 330–339. [Google Scholar]
  16. Huber, M.; Hicks, M.A.; Vermeer, P.A.; Moormann, C. Probabilistic calculation of differential settlement due to tunneling. In Proceedings of the 8th International Probabilistic Workshop, Akademia Morska, Szczecin, Poland, 18–19 November 2010; pp. 1–13. [Google Scholar]
  17. Rowe, R.K.; Kack, G.J. A theoretical examination of the settlements induced by tunnelling: Four case histories. Can. Geotech. J. 1983, 20, 299–314. [Google Scholar] [CrossRef]
  18. Loganathan, N.; Poulos, H.G. Analytical Prediction for Tunneling-Induced Ground Movements in Clays. J. Geotech. Geoenviron. Eng. 1998, 124, 846–856. [Google Scholar] [CrossRef]
  19. Oh, J.Y.; Ziegler, M. Investigation on Influence of Tail Void Grouting on the Surface Settlements during Shield Tunneling Using a Stress-Pore Pressure Coupled Analysis. KSCE J. Civ. Eng. 2014, 18, 803–811. [Google Scholar] [CrossRef]
  20. Teachavorasinskun, S. Evaluation of liner responses due to non-uniform tail void grouting pressure. Geomech. Geoengin. 2018, 13, 226–232. [Google Scholar] [CrossRef]
  21. Liao, S.; Xu, J.; Jiao, Q. Identifying TBM tunneling defects through recognition of ground settlement patterns. In Proceedings of the GeoCongress 2012: State of the Art and Practice in Geotechnical Engineering, Oakland, CA, USA, 25–29 March 2012; pp. 3179–3188. [Google Scholar]
  22. Wang, Z.; Xie, X. Design and optimization of the antenna applied for detecting the voids in tunnel lining by GPR. In Proceedings of the 2012 14th International Conference on Ground Penetrating Radar, Shanghai, China, 4–8 June 2012; IEEE: New York, NY, USA, 2012. [Google Scholar]
  23. Yu, Q.M.; Zhou, H.L.; Wang, Y.H.; Duan, R.X. Quality monitoring of metro grouting behind segment using ground penetrating radar. Constr. Build. Mater. 2016, 110, 189–200. [Google Scholar] [CrossRef]
  24. Ding, W.; Duan, C.; Zhu, Y.; Zhao, T.; Huang, D.; Li, P. The behavior of synchronous grouting in a quasi-rectangular shield tunnel based on a large visualized model test. Tunn. Undergr. Space Technol. 2019, 83, 409–424. [Google Scholar] [CrossRef]
  25. Zhang, Z.X.; Liu, C.; Huang, X.; Kwok, C.Y.; Teng, L. Three-dimensional finite-element analysis on ground responses during twin-tunnel construction using the URUP method. Tunn. Undergr. Space Technol. 2016, 58, 133–146. [Google Scholar] [CrossRef] [Green Version]
  26. Wen-Jie, S.; Yan-Yong, X. A continuous upper bound limit analysis for the transverse stability of shallow shield tunnelling in undrained clays. Revue Française De Génie Civil 2020, 24, 401–422. [Google Scholar] [CrossRef]
  27. Jiang, H.; Jiang, Y.S.; Huang, M.L.; Chen, Z.Q. Study on Soil Improvement of TBM’S Commencement and Arrival in Water-Rich Sand Stratum. Appl. Mech. Mater. 2011, 90–93, 2119–2124. [Google Scholar] [CrossRef]
  28. Ching, J.; Phoon, K.-K.; Chen, Y.-C. Reducing shear strength uncertainties in clays by multivariate correlations. Can. Geotech. J. 2010, 47, 16–33. [Google Scholar] [CrossRef]
  29. Papaioannou, I.; Straub, D. Reliability updating in geotechnical engineering including spatial variability of soil. Comput. Geotech. 2012, 42, 44–51. [Google Scholar] [CrossRef] [Green Version]
  30. Li, X.-Y.; Zhang, L.-M.; Gao, L.; Zhu, H. Simplified slope reliability analysis considering spatial soil variability. Eng. Geol. 2017, 216, 90–97. [Google Scholar] [CrossRef]
  31. Jiang, S.; Li, D.; Cao, Z.; Zhou, C.; Phoon, K.K. Efficient system reliability analysis of slope stability in spatially variable soils using Monte Carlo simulation. J. Geotech. Geoenviron. Eng. 2014, 141, 04014096. [Google Scholar] [CrossRef]
  32. Phoon, K.-K.; Tang, C. Characterisation of geotechnical model uncertainty. Georisk Assess. Manag. Risk Eng. Syst. Geohazards 2019, 13, 101–130. [Google Scholar] [CrossRef]
  33. Tang, C.; Phoon, K.-K. Characterization of model uncertainty in predicting axial resistance of piles driven into clay. Can. Geotech. J. 2019, 56, 1098–1118. [Google Scholar] [CrossRef]
  34. Phoon, K.K.; Kulhawy, F.H. Characterization of geotechnical variability. Can. Geotech. J. 1999, 36, 612–624. [Google Scholar] [CrossRef]
  35. Fenton, G.A.; Griffiths, D. Bearing-capacity prediction of spatially random soils. Can. Geotech. J. 2003, 40, 54–65. [Google Scholar] [CrossRef]
  36. Lloret-Cabot, M.; Fenton, G.A.; Hicks, M.A. On the estimation of scale of fluctuation in geostatistics. Georisk Assess. Manag. Risk Eng. Syst. Geohazards 2014, 8, 129–140. [Google Scholar] [CrossRef] [Green Version]
  37. Fenton, G.A.; Griffiths, D.V. Probabilistic foundation settlement on spatially random soil. J. Geotech. Geoenviron. Eng. 2002, 128, 381–390. [Google Scholar] [CrossRef] [Green Version]
  38. BSI. BS EN 1997–2:2007: Eurocode 7. Geotechnical Design. Ground Investigation and Testing; BSI: London, UK, 2007. [Google Scholar]
Figure 1. Schematic diagram of backfill-grouting-layer distribution.
Figure 1. Schematic diagram of backfill-grouting-layer distribution.
Applsci 10 05002 g001
Figure 2. Distribution of the grouting-layer thickness.
Figure 2. Distribution of the grouting-layer thickness.
Applsci 10 05002 g002
Figure 3. Displacement convergence patterns for different grouting-layer distributions.
Figure 3. Displacement convergence patterns for different grouting-layer distributions.
Applsci 10 05002 g003
Figure 4. Model geometry and mesh gird.
Figure 4. Model geometry and mesh gird.
Applsci 10 05002 g004
Figure 5. Displacement-constrained field.
Figure 5. Displacement-constrained field.
Applsci 10 05002 g005
Figure 6. Validation of mean and standard deviation.
Figure 6. Validation of mean and standard deviation.
Applsci 10 05002 g006
Figure 7. Anderson–Darling goodness-of-fit test of normal distribution.
Figure 7. Anderson–Darling goodness-of-fit test of normal distribution.
Applsci 10 05002 g007
Figure 8. Effect of the soil COVs on settlement. (a) COV(c); (b) COV(φ); (c) COV(E).
Figure 8. Effect of the soil COVs on settlement. (a) COV(c); (b) COV(φ); (c) COV(E).
Applsci 10 05002 g008aApplsci 10 05002 g008b
Figure 9. Effect of the soil SOFs on settlement. (a) s v ; (b) s h .
Figure 9. Effect of the soil SOFs on settlement. (a) s v ; (b) s h .
Applsci 10 05002 g009
Figure 10. Effect of the grouting-layer parameters on settlement. (a) Grouting-layer thickness; (b) COV(t).
Figure 10. Effect of the grouting-layer parameters on settlement. (a) Grouting-layer thickness; (b) COV(t).
Applsci 10 05002 g010
Figure 11. Predicted curve for ground-surface settlement.
Figure 11. Predicted curve for ground-surface settlement.
Applsci 10 05002 g011
Table 1. Statistic parameters of grouting-layer thickness.
Table 1. Statistic parameters of grouting-layer thickness.
Source.MeanSpanCOVTest Method
Yu et al. (2016) [23]300 mm220 mm–380 mm0.1–0.20Ground penetrating radar (GPR)
Ding et al. (2019) [24]169 mm105 mm–240 mm0.15–0.40Large visualized model test
Note: COV refers to the coefficient of variation.
Table 2. Input parameters for the reference case.
Table 2. Input parameters for the reference case.
Deterministic VariablesValueInput VariablesMean Value (μ)COV (%)SOF (m)
s h s v
Dl (m)5.6c (kPa)100.20102.0
D (m)6.0φ (°)250.10
g (mm)200E (MPa)200.15
γs (kN/m3)20 δ (mm)200.20
v0.2 t (mm)180
Note: SOF refers to the scale of fluctuation.
Table 3. Reference case and parametric study.
Table 3. Reference case and parametric study.
RemarksCOV (c)COV (φ)COV (E)SOF (soil)/m t (mm) δ (mm)COV (t)
svsh
Reference case0.20.10.152.010180200.2
Effect of COV of soil cohesion0.1 *0.10.152.010180200.2
0.2 *0.10.152.010180200.2
0.3 *0.10.152.010180200.2
Effect of COV of soil friction angel0.20.1 *0.152.010180200.2
0.20.2 *0.152.010180200.2
0.20.3 *0.152.010180200.2
Effect of COV of soil Young’s modulus0.20.10.05 *2.010180200.2
0.20.10.15 *2.010180200.2
0.20.10.25 *2.010180200.2
Effect of soil SOF0.20.10.152.0 *10180200.2
0.20.10.152.0 *10180200.2
0.20.10.153.0 *10180200.2
0.20.10.152.010 *180200.2
0.20.10.152.020 *180200.2
0.20.10.152.030 *180200.2
Effect of grouting-layer thickness0.20.10.152.010190 *10 *0.2
0.20.10.152.010180 *20 *0.2
0.20.10.152.010170 *30 *0.2
Effect of COV of grouting layer0.20.10.152.010180200.1 *
0.20.10.152.010180200.2 *
0.20.10.152.010180200.3 *
0.20.10.152.010180200.4 *
Note: “*” refer to association of variables.

Share and Cite

MDPI and ACS Style

Wang, Z.; Jin, D.; Shi, C. Spatial Variability of Grouting Layer of Shield Tunnel and Its Effect on Ground Settlement. Appl. Sci. 2020, 10, 5002. https://doi.org/10.3390/app10145002

AMA Style

Wang Z, Jin D, Shi C. Spatial Variability of Grouting Layer of Shield Tunnel and Its Effect on Ground Settlement. Applied Sciences. 2020; 10(14):5002. https://doi.org/10.3390/app10145002

Chicago/Turabian Style

Wang, Zhongzheng, Dalong Jin, and Chenghua Shi. 2020. "Spatial Variability of Grouting Layer of Shield Tunnel and Its Effect on Ground Settlement" Applied Sciences 10, no. 14: 5002. https://doi.org/10.3390/app10145002

APA Style

Wang, Z., Jin, D., & Shi, C. (2020). Spatial Variability of Grouting Layer of Shield Tunnel and Its Effect on Ground Settlement. Applied Sciences, 10(14), 5002. https://doi.org/10.3390/app10145002

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