Citation (APA) Vu, M. N., Broere, W., & Bosch, J. (2016). Volume loss in shallow tunnelling. Tunnelling and Underground Space Technology, 59,

Transcription

1 Delft University of Technology Volume loss in shallow tunnelling Vu, M.N.; Broere, Wout; Bosch, Johan Publication date 216 Document Version Accepted author manuscript Published in Tunnelling and Underground Space Technology Citation (APA) Vu, M. N., Broere, W., & Bosch, J. (216). Volume loss in shallow tunnelling. Tunnelling and Underground Space Technology, 59, Important note To cite this publication, please use the final published version (if applicable). Please check the document version above. Copyright Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons. Takedown policy Please contact us and provide details if you believe this document breaches copyrights. We will remove access to the work immediately and investigate your claim. This work is downloaded from Delft University of Technology. For technical reasons the number of authors shown on this cover page is limited to a maximum of 1.

2 218 Manuscript version made available under CC-BY-NC-ND 4. license Volume loss in shallow tunnelling Minh Ngan Vu a,b,, Wout Broere a, Johan Bosch a a Delft University of Technology, Geo-Engineering Section, Stevinweg 1, 2628 CN Delft, the Netherlands b Hanoi University of Mining and Geology, Faculty of Civil Engineering, Vietnam Abstract Although volume loss has an important effect in estimating the ground movements due to tunnelling in the design stage, this parameter is often determined by experience. This makes it difficult to estimate the impact on volume loss when changing project parameters like soil conditions, depth of the tunnel or sensitivity of the surroundings. This paper investigates the relationship between volume loss and cover-to-diameter ratio in shallow tunnelling. Based on a number of (empirical) relations from literature, such as the stability number method and an analysis of the bentonite and grout flows, volume loss at the face, along the shield and at the tail is determined. Long-term volume loss behind the shield is also estimated by means of consolidation. In this way a band width of achievable volume loss for future projects is derived. Keywords: cover-to-diameter ratio; volume loss; TBM; tail void; shallow bored tunnel 1. Introduction Tunnelling often leads to settlements of the soil surface due to overexcavation, soil relaxation and inefficient tail void filling. The magnitude of volume loss is influenced by tunnelling management, characteristics of the tunnel boring machines (TBM), and the geotechnical conditions. In predictions of surface settlement (Peck, 1969) and subsurface settlement (Mair et al., 1993), the volume loss is often determined by engineering experience Corresponding author: Tel.: ; Fax: address: N.Vuminh@tudelft.nl/vuminhngan@gmail.com (Minh Ngan Vu) Preprint submitted to Tunnelling and Underground Space Technology March 25, 216

3 and data from previous cases. This makes it difficult to correctly assess the volume loss for a future project under radically different conditions like a shallow depth of the tunnel and/or very different soil parameters. A ground movement analysis in Vu et al. (215a) shows the important role of volume loss for settlement calculations and in predicting the effects on existing buildings induced by tunnelling. Especially for (very) shallow tunnels near building foundations, the impact of changes in volume loss is large. Most previous studies on volume loss start from a given volume loss and establish deformation patterns from that or correlate surface observations to volume loss at the tunnel for specific projects. Mair et al. (1982); Attewell et al. (1986); Macklin (1999); Dimmock and Mair (27) studied the volume loss with a summary of projects in overconsolidated clay relating to the volume loss at the tunnelling face. Verruijt and Booker (1996); Verruijt (1997); Strack (22) applied analytical methods for predicting the ground loss around the tunnel. Loganathan (211) proposed volume loss calculations but only approximated volume loss along the shield with the worst case, and does not take the consolidation into account. Meanwhile, Bezuijen and Talmon (28) showed the effect of grouting pressure on the volume loss around the TBM but none of these includes a detailed method to estimate volume loss along the TBM. This paper aims to estimate the volume loss when tunnelling with limited ratios (i.e. less than 1) in various soils with a focus on slurry shield tunnelling. On the basis of the studies by Attewell and Farmer (1974), Cording and Hansmire (1975) and Mair and Taylor (1999), the volume loss in the tunnelling progress can be estimated by the sum of the following components as shown in Figure 1: - Volume loss at the tunnelling face: soil movement towards the excavation chamber as a result of movement and relaxation ahead of the face, depending on the applied support pressures at the tunnelling face; - Volume loss along the shield: the radial ground loss around the tunnel shield due to the moving soil into the gap between the shield and surrounding soil, which can be caused by overcutting and shield shape. The bentonite used in the tunnelling face flows into the gap, while the grout used in the shield tail also flows in the opposite direction. Due to the drop of bentonite and grout flow pressures in a constrained gap, soil can still move into the cavity when the soil pressure is larger than the bentonite pressure or grout pressure; - Volume loss at the tail: when precast segments are placed, the advance of the shield results an annular cavity between the segments and surrounding 2

4 Figure 1: Volume loss components soil. Grout is used in order to prevent surrounding soil moving into the gap. Volume loss at the tail depends on applied grouting pressure at the tail and proper volume control, where high grout volume and pressure may lead to local heave and low volume to increase settlements as indicated in Figure 1; - Volume loss behind the shield tail due to consolidation: in this void along the tunnel lining, grout consolidates and forms a grout cake, and the stress changes induced in the soil may lead to long-term consolidation settlements in soil volume above the tunnel. Other causes of volume loss are shrinkage of grout and long-term lining deformations. However, their contributions to the total volume loss are small comparing to the above factors. The total volume loss V L in tunnelling progress can be given as: V L = V L,f +V L,s +V L,t +V L,c (1) where V L,f is volume loss at the tunnelling face, V L,s is volume loss along the shield, V L,t is volume loss at the tail, and V L,c is volume loss due to consolidation. To illustrate the impact of the different contributions in different soil conditions, estimates are made for a number of ideal soil profiles which are derived from Amsterdam North-South metro line project (Gemeente-Amsterdam, 3

5 29), consisting of a single soil type with most important properties as defined in Table 1, where γ is volumetric weight, ϕ is the friction angle, K is the initial coefficient of lateral earth pressure, c is cohesion, C s is compression constant, C swel is swelling constant, ν is Poisson s ratio and E s is the stiffness modulus of the ground. Table 1: Soil parameters used in design of Amsterdam North-South metro line project (Bosch and Broere, 29; Gemeente-Amsterdam, 29) Soil type γ(kn/m 3 ) ϕ( o ) K( ) c(kn/m 2 ) C s ( ) C swel ( ) ν( ) E s (kn/m 2 ) Sand Clayey sand Clay Organic clay Peat Volume loss at the tunnelling face When tunnelling, the soil ahead of the excavation chamber generally has the trend to move into the cavity which is created by the tunnelling machine. The soil volume moving towards the face depends on applied support pressures and can be controlled by adjusting the support pressures. In stability analysis for tunnelling, the stability number N proposed by Broms and Bennermark (1967) is widely used. By studying the relationship between this stability number and volume loss at tunnelling face, Attewell et al. (1986), Mair et al. (1982), Mair (1989), Macklin (1999) and Dimmock and Mair (27) presented a method to determine the expected volume loss based on observed data. The stability number N is given by: N = γ(c +D/2) s c u (2) where s is the support pressure at the tunnelling face and c u is undrained shear strength of the soil. In shallow tunnelling, the support pressure at the tunnelling face should be high enough to avoid the collapse to the excavation chamber but also limited to prevent blow-out and fracturing. Firstly, the required support 4

6 pressure must be higher than or at least equal to the total of water pressure and horizontal effective soil pressure taking into account three dimensional arching effects. The wedge model, which was studied by Anagnostou and Kovári (1994), Jancsecz and Steiner (1994) and Broere (21), is commonly applied to determine the minimum support pressure s min. In the case of shallow tunnelling, the minimum support pressure s min can be derived from the wedge model, as follows: s min = σ, h +p = σ, vk A3 +p = γ, zk A3 +p (3) where p is pore pressure and K A3 is the three dimensional earth pressure coefficient determined in Jancsecz and Steiner (1994). Secondly, the maximum support pressures are often estimated as to avoid blow-out and fracturing. According to Vu et al. (215b), the maximum support pressures in the case of blow-out are given by: s,t,max = γ (H π ) 8 D +2 H ( c+ 1 ) D 2 HK yγ tanϕ ad (4) 4 s,b,max = γ (H π ) 8 D +2 H D ( c+ 1 ) 2 HK yγ tanϕ +γ T πd+ ad 4 (5) where s,t,max and s,b,max are the maximum support pressures at the top and bottom of the tunnel. In normally consolidated soil, according to Mori et al. (1991), the maximum pressure in the case of fracturing is presented as: s f = σ, vk +p+c u (6) However, field data show that the higher allowable support pressures are often applied in the tunnelling face, according to NEN-365 (212) and reports by BTL (Boren van Tunnels en Leidingen), in the Netherlands. Therefore, the support pressures boundaries are determined with the minimum support pressure and the maximum support pressure in the case of blow-out as indicated in Equations 4 and 5 above. Figure 2 shows the relationship between the required support pressures and the ratio with the tunnel diameter D = 6m in clay. We will elaborate the calculation method for these conditions and present overall results for different diameters and soil conditions later in Figure 1. According to Vu et al. (215b), only ratios larger than.4 are studied, as less cover 5

7 s min s,t,max s,b,max s fracturing support pressure(kn/m 2 ) Figure 2: The range of support pressures at the tunnelling face of a tunnel with a diameter D = 6m in clay would lead to unreasonable large volume loss, and the upper 3 to 4 meter of soil in urban areas are often taken up by various utilities and therefore would not be available for tunnelling. The support pressures calculated here are the minimum support pressure from a wedge model and the maximum support pressures for fracturing and blow-out at the top and the bottom of the tunnel. Figure 3 shows the calculated stability number N for these support pressures. Since the applied support pressures are derived from the wedge stability model, fracturing and blow-out conditions, N values in this figure are smaller than 2. It means that the tunnelling face is stable with these support pressures. O Reilly (1988) indicated that a relation exists between the volume loss at tunnelling face V L,f and the load factor LF, which is estimated by the ratio of working stability number N and the stability number at collapse N TC, as follows: LF = N N TC (7) 6

8 with s min with s,t,max with s fracturing N Figure 3: The range of stability number N for a tunnel with a diameter D = 6m in clay where N TC is estimated according to Mair and Taylor (1999) as: For 1: ( ) 2C N TC = 2+2ln D +1 For 1 1.8: ( ) 2C N TC = 4ln D +1 Figure 4 shows the relationship between the load factor LF and the ratio. The load factor is less than.6 for the minimum support pressure and has the trend of reducing when the ratio increases. This means that the tunnel becomes safer with regards to estimating the support pressures when the ratio becomes larger. From the analysis of case history data of the load factor LF and the volume loss at the tunnelling face (Figure 5), Macklin (1999) presented a formula to calculate the volume loss at the tunnelling face V L,f as: V L,f (%) =.23e 4.4LF (8) Equation 8 can be used to convert the load factor LF to the volume loss V L,f 7

9 .6 with s min.5 with s,t,max with s fracturing.4 LF Figure 4: Relationship between load factor LF and ratio for a tunnel with diameter D=6m in clay estimates, which leads to Figure 6. This shows the range of volume loss at the tunnelling face V L,f with various ratios for a tunnel with D = 6m in clay. In shallow tunnels with.4 1 the range of possible volume loss V L,f is large, ranging from.12% to 3.1%. This means that if tunnelling uses the minimum pressure in the excavation chambers, the volume loss V L,f will increase significantly. Meanwhile, the volume loss V L,f in the case of 1 2 ranges from.27% to 1.5%. The difference in volume loss V L,f between the minimum pressure and maximum pressures due to blow-out and fracturing is clearly reduced. Therefore, in the case of very shallow tunnels ( 1) the support pressures applied at the tunnelling face should be kept near to the maximum pressure in order to avoid increasing the volume loss. 3. Volume loss along the shield The diameter of the cutting wheel in front of the TBM is often larger than the diameter of the shield. This leads to an overcut when tunnelling (Figure 7). Also, the TBM is often tapered, which creates a gap between 8

10 Figure 5: Volume loss V L,f and load factor LF (Macklin, 1999) 9

11 3.5 with s min with s,t,max with s,b,max with s fracturing Vl(%) Figure 6: Volume loss at tunnelling face for a tunnel with diameter D = 6m in clay the shield skin and the surrounding soil. Additional gapping can also occur when the TBM moves in curves as indicated in Festa et al. (215). In this study, the effect of curves is not included. This gap is often filled by bentonite, which flows from the tunnelling face and/or grout which comes from the shield tail. In practice, the grout and bentonite pressures are often larger than the vertical soil pressure at the tunnelling face and tail. From the observation of Bezuijen (27), there are three possible bentonite and grout flows that can occur along the shield when tunnelling. Firstly, the bentonite flows from the tunnelling face to the tail and pushes the grout at the joint between the tail and the TBM. Secondly, the grout flows from the tail to the tunnelling face and pushes the bentonite away. Thirdly, the grout flows from the tail to the tunnelling face and the bentonite also flows in the opposite direction. The flows of bentonite and grout were also simulated in Nagel and Meschke (211). In shallow tunnelling, due to the possibility of blow-out and fracturing, there is a limitation of applied grout and bentonite pressures at the tunnelling face and the tail. According to Bezuijen (27), both liquids, the grout in the shield tail and the bentonite applied at the tunnelling face are assumed to behave as the 1

12 Figure 7: Bentonite and grout flows along the shield and lining segments Table 2: Input parameters of tunnel boring machine Diameter of shield D 6m Length-to-diameter P/D ratio of the shield 1 Reduction of shield diameter a.2% Overcutting h overcut.15m Shear strength of grout τy grout 1.6kP a Shear strength of bentonite τy bentonite.8kp a Bingham liquids, such that the yield stress is governing in the flow behaviour. The flow pressures in grout and bentonite reduce along the shield as in Figure 7. The reduction of grout pressure along the shield is given by: p = x w j τ y (9) where p is the change of the pressure due to flow, x is a length increment alongthetbm,w j isthejointwidthbetweenthetunnelandthesurrounding soil and τ y is a shear strength of the grout around the TBM. In this study, the volume loss along the shield is calculated with input parameters as indicated in Table 2 with the following approach. As an example, the calculation is carried out with a case of tunnel with 11

13 14 12 s,t,max grout pressure (kn/m 2 ) VL=.7% 4 VL=.56% VL=.64% smin sfracturing 2 grout pressure bentonite pressure soil pressure distance from tail (m) Figure 8: Bentonite and grout pressures along a shield with D = 6m in clay D = 6m and =.75 in clay. Figure 8 shows the change of grout pressure and bentonite pressure along the shield. It is assumed that when the grout pressure and bentonite pressure are less than the vertical soil pressure, the soil is moving into the cavity. The volume loss is estimated as the void volume that is filled by soil. The volume loss will not occur if the grout pressure and the bentonite pressure are larger than the vertical soil pressure. In that case, the gap along the shield is assumed to be completely filled by grout and bentonite. From this figure, the volume loss along the shield depends on the bentonite pressure, which is applied at the tunnelling face and the grout pressure at the tail. When the bentonite and grout pressures are equal to the minimum required pressure as calculated in previous section, the volume loss will be maximal. On the other hand, when the maximum allowable pressures are applied, there is no volume loss along the shield. In order to investigate the effect of grout pressure on the shield in different soils, we assume the bentonite pressure applied at the tunnelling face is the average of the minimum support pressure and the maximum support pressure for fracturing. Figure 9 shows the change of volume loss with different tunnel diameters in various soils. When the ratio increases, the range of 12

14 volume loss(%) D=6 D=8 D= (a) in sand volume loss (%) D=6 D=8 D= (b) in clayey sand volume loss(%).4.3 volume loss(%) D=6 D=8 D= D=6 D=8 D= (c) in clay (d) in organic clay volume loss(%) D=6 D=8 D= (e) in peat Figure 9: Volume loss along the shield in various soils 13

15 the volume loss along the shield V L,s is larger. With a particular ratio of the tunnel, the larger the tunnel diameter is, the smaller the volume loss V L,s is with the constant overcut parameter of the shield. In the case of tunnelling in sand, the upper boundary of the volume loss V L,s increases nearly linearly with the ratio from.4 to.8, then becomes almost constant when the ratio increases whereas the lower boundary reduces linearly when the ratio increases. This also appears in the cases of tunnelling in clay and organic clay. At this point, basically, the entire annulus is filled by the surrounding soil, leading to a maximum attainable volume loss along the tail. In the case of tunnelling in peat, in the range of this analysis with.4 2, the upper boundary of the volume loss V L,s lightly rises and the lower boundary linearly decreases. It is noted that the maximum upper boundary volume loss along the shield V L,s is the same for a given tunnel diameter. Regardless of soil conditions, for a tunnel with D = 6m, it follows that V L,s,max =.7%, with D = 8m, V L,s,max =.57% and with D = 1m, V L,s,max =.5%. Figure 1 shows the boundary of the volume loss along the shield V L,s in relationship with ratios for different tunnel diameters in different soils. The upper boundary for.4.6 corresponds to the case of tunnelling in peatinallthreetunneldiameters. Whentunnellingwith.6 1, the upper boundary is given by tunnelling in organic clay and when tunnelling with the ratio larger than 1, the upper boundary becomes constant and depends on the tunnel diameter D. The maximum volume loss along the shield V L,s is about.7% for D = 6m, V L,s,max =.57% for D = 8m and V L,s,max =.5% for D = 1m. For the lower boundary, there is a decreasing trend of the minimum volume loss along the shield V L,s when the ratio increases. In the case of D = 6m the maximum V L,s of the lower boundary is about.47% when =.5. The maximum volume loss along the shield V L,s of the lower boundary is about.38% with D = 8m and.32% with D = 1m when =.5. When = 2, V L,s,max =.2% for D = 6m, V L,s,max =.13% for D = 8m and V L,s,max =.9% for D = 1m. 4. Volume loss behind the shield When precast segments are placed, the advance of the shield results in an annular cavity between the segments and the surrounding soil due to the shape of the TBM and the overcut as discussed above. Grout is injected 14

16 .7.65 upper boundary.6.55 volume loss(%) lower boundary.35.3 in sand in clay.25 in peat (a) D = 6m.6.55 upper boundary.5.45 volume loss(%) lower boundary.25.2 in sand in clay.15 in peat (b) D = 8m.5.45 upper boundary.4.35 volume loss(%) lower boundary.15 in sand.1 in clay in peat (c) D = 1m Figure 1: Volume loss along the shield with different tunnel diameters 15

17 rapidly in order to prevent the surrounding soil to move into the gap. It is assumed that the void is filled by the grout. The injected grout pressure induces the loading on the soil around the tunnel lining. This might lead to immediate displacements and long-term consolidation of the soil. These are two components of the volume loss behind the shield: the volume loss at the tail and the volume loss due to consolidation Volume loss at the tail When the grout is injected with high pressures at the tail, the ground around the tunnel will be deformed. In order to estimate the surface settlement induced by tunnelling, there are some analytical solutions proposed by Sagaseta (1988); Verruijt (1997); Strack (22) based on cavity expansion and taking the influence of a free surface into account. However, the effect of the range of support pressures has not taken into account in these methods and resulting solutions, for instance expressed as a Laurent series expansion in the case of Verruijt (1997), require an increasing number of terms for a stable numerical integration if the distance between free surface and tunnel reduces. On the other hand, the cavity expansion developed for the case of a cavity in infinite medium has been implemented in tunnelling studies by Taylor (1993); Yu (213) and results in far more elegant and practical solution for a first estimate of the effect of grout pressures on soil stresses and deformations around the TBM. To determine the effect of grouting at the tail on volume loss at the tail and consolidation, in this study, the cavity expansion method for tunnelling, which is proposed by Yu (213), is therefore applied as a simplified method. In this cavity-expansion theory, it is assumed that the soil around the tunnel is a Tresca medium. The stresses in the soil and the settlement at the surface can be calculated by the cavity-expansion theory. According to Yu (213), the plastic zone will deform around the tunnel wall, as can be seen in Figure 11, with the radius R p of the plastic zone estimated from the following equation: R p = D 2 exp ( p s Y k 1+k k ) (1) where p is the pre-tunnelling pressure; k = 1 or 2 corresponding to cylindrical or spherical cavity models; Y = 2c u or 2c u corresponding to the case of contraction or expansion of the tunnel. 16

18 Similar to Yu (213) and Taylor (1993), the pre-tunnelling pressure p can be estimated as: ( p = γ C + D ) (11) 2 The soil displacement u s in the elastic zone is given by: Figure 11: Deformations around a shallow tunnel at the tail Yr u s = 2(k +1)G ( Rp r ) 1+k (12) where r is the distance from the calculated point to the tunnel centre and G = E/2(1+ν) is the shear modulus of soil. The soil displacement u s in the plastic zone is given by: ( ) k ( ) Y D D u s = 2(k +1)G 2r 2 exp (1+k)(p s) 1 (13) ky In this case, the effect of grouting pressures at the tail is analysed with a cylindrical cavity model and is calculated with the minimum and maximum 17

19 support pressures. Thus, in Equation 1, k equals 1. It can be assumed that the volume loss around the tunnel due to grouting at the tail equals the volume of ground settlement at the surface. In order to identify the contribution of soil deformation at the tail on the total volume loss, the displacement of the ground surface is estimated. According to assumptions in Yu (213), the tunnel will collapse when the plastic zone expands to the ground surface. It means that when the tunnel is stable, the radial displacement of ground u s at the surface is in the elastic zone and can be calculated with Equation 12. The surface settlement at the tail can be estimated from: u t = u s sinθ (14) where θ is the angle between the calculated point to the tunnel centre and the horizontal axis (see Figure 11). The volume loss at the tail V L,t can be estimated as: V L,t = V s,t π(d/2) 2 (15) where V s,t is the volume of the surface settlement due to grouting pressures at the tail (see Figure 11). Figure 12 shows the boundaries of the volume loss at the tail V L,t in various soils for tunnels with D = 6,8 and 1m with the range of the support pressures from the vertical soil stress to the maximum support pressure at the top of the tunnel derived from Equation 4. The figure shows that the larger the tunnel diameter is, the larger the range of volume loss V L,t is. When the support pressure equals the vertical soil stress at the top of the tunnel lining, there is a contraction in the cavity and this leads to positive values of the lower boundary of volume loss at the tail. When a high support pressure is used, the cavity will expand. The negative volume loss V L,t values indicate that the soil above the tunnel lining is pushed upward and there might be heave at the ground surface. In practice, this heave might not be observed becausethesettlementduetovolumelossatthetunnellingfaceandalongthe shield could be larger. When a high support pressure is applied at the tail, a heave can occur in order to compensate the volume loss at the tunnelling face and along the TBM as can be seen in Figures 12a, 12b, 12c and 12d. However, in the case of very shallow tunnelling, there is no heave due to the small margin in the range of allowable support pressures as indicated in Vu 18

20 .4 lower boundary for D=1m.5 lower boundary for D=1m.2 volume loss at the tail V L, t (%) lower boundary for D=6m lower boundary for D=8m upper boundary for D=6m upper boundary for D=8m upper boundary for D=1m volume loss at the tail V L, t (%) -.5 lower boundary for D=6m lower boundary for D=8m upper boundary for D=6m upper boundary for D=8m upper boundary for D=1m (a) in sand (b) in clayey sand 1 4 lower boundary for D=1m.5 3 lower boundary for D=1m volume loss at the tail V L, t (%) lower boundary for D=6m lower boundary for D=8m upper boundary for D=6m upper boundary for D=8m upper boundary for D=1m volume loss at the tail V L, t (%) lower boundary for D=8m lower boundary for D=6m upper boundary for D=6m upper boundary for D=8m upper boundary for D=1m (c) in clay (d) in organic clay 18 volume loss at the tail V L, t (%) lower boundary for D=1m upper boundary for D=1m lower boundary for D=8m upper boundary for D=8m lower boundary for D=6m upper boundary for D=6m (e) in peat Figure 12: Volume loss at the tail with different tunnel diameters in various soils 19

21 et al. (215b). In Figure 12e, when tunnelling in peat, the volume loss at the tail is positive with a high value, especially in the case of a tunnel diameter D = 1m. It means that shallow tunnelling with a large diameter in peat might be difficult due to the large expected volume loss. This conclusion coincides with the conclusion indicated in Vu et al. (215b) for the range of support pressure for shallow tunnelling in peat with low ratios. Figure 13 shows the dependence of V L,t values on soils in various tunnel diameters. When tunnelling in peat, the range of V L,t values is significantly large compared to tunnelling in sand, clay and organic clay, especially in the case of tunnels with large diameters as indicated above Volume loss due to consolidation For the volume loss due to consolidation, in the cavity behind the tail, two consolidation processes occur along the tunnel lining. Firstly, the newly injected grout is consolidating and forms a consolidated grout cake in the cavity along the tunnel lining (Talmon and Bezuijen, 29). In the case of tunnelling in clay, the consolidation in grout might not occur (Bezuijen and Talmon, 213) and the length of liquid grout on the lining is much longer. Although the grout pressure decreases along the lining, the injected grout may flow along 2 to 3 following segments and the appearance of the grout cake will prevent the movement of the soil above. It is often assumed that there is no volume loss in the grout consolidating. The other volume loss is due to the subsequent shrinkage of grout, which is estimated at about 7 to 1 percent of total tail gap according to Loganathan (211). However, the contribution of this volume loss to the total volume loss is small comparing to the other volume losses. This volume loss, therefore, is not taken into account in this study. The second process is the consolidation of the soil volume above the tunnel behind the tail. When grout is applied at the tail, the soil stress in the above soil volume will change. This will induce consolidation in the long term behind the tail. The volume loss due to consolidation V L,c is derived from the consolidation settlement of the soil volume above tunnel lining. In the case of tunnelling in sand, consolidation of soil will probably not occur or be minimal. For tunnels in clay or peat, this may be a notable contribution. The consolidation settlement of the soil can be estimated from Terzaghi s formula, as follows: u c = 1 ( ) σsoil ln (16) C s 2 σ

22 14 volume loss at the tail V L, t (%) in sand in clay in peat (a) D = 6m 6 volume loss at the tail V L, t (%) in sand in clay in peat (b) D = 8m in sand volume loss at the tail V L, t (%) in clay in peat (c) D = 1m Figure 13: Volume loss at the tail with different tunnel diameters 21

23 Figure 14: Soil stresses at the tail where C s is the compression constant depending on soil type (as can be seen in Table 1), σ soil is the vertical stress in the soil and σ is the initial vertical stress in the soil. In case the vertical stress is lower than the initial vertical stress, unloading occurs and Equation 16 would be modified to: u c = 1 C swel ln ( σsoil σ ) (17) where C swel is the swelling constant depending on soil type (as can be seen in Table 1). The stress in the soil σ soil is estimated from the radial and tangential stresses derivedbythecavityexpansiontheoryascanbeseeninfigure14. According to Yu (213), the stresses in the elastic zone are given by: σ r = p + ky 1+k σ θ = p Y 1+k 22 ( Rp r ( Rp r ) (1+k) (18) ) (1+k) (19)

24 where σ r and σ θ are the radial and tangential stresses as shown in Figure 14. In the plastic zone, the stresses are given by: σ r = p + ky 1+k +ky ln R p r (2) σ θ = p Y 1+k +ky ln R p (21) r In order to estimate the consolidation settlement, the soil volume above the tunnel lining is divided into n layers. The final consolidation settlement is derived by summing deformations of these layers, which are calculated by Equations 16 and 17. The final consolidation settlement is given by: u j c = n i=1 u (j,i) c z (22) where u c (j,i) and z (i) are the deformation due to consolidation and the depth of the i th layer at the j th location along the surface. Figure 15 shows settlement troughs in the case of tunnelling in clay with a diameter D = 1m and the ratio of = 1, as an example. It can be seen that a heave and a settlement can occur depending on what particular support pressure is applied. By integrating the final consolidation settlements over the surface, the volume of consolidation settlement at the surface V cons can be estimated as: V cons = m j=1 u (j) c x (23) where x is a length increment along the surface consolidation settlement and m is the increment number. The volume loss due to consolidation settlement is then estimated as: V L,c = V cons π(d/2) 2 (24) Figure 16 shows the relationship between the consolidation volume loss V L,c and the ratio for tunnels with different diameters in clay, organic clay and peat. With in the case of tunnelling in clay and 23

25 .5 consolidation settlement(m) lower boundary support pressure upper boundary support pressure distance to the tunnel axis(m) Figure 15: Consolidation settlement troughs when tunnelling in clay with D = 1m and = in the case of tunnelling in organic clay, it can be seen that the maximum support pressure applied at the tail can lead to a heave on the surface. The volume loss due to consolidation V L,c when maximum support pressure is applied becomes smaller than when minimum support pressure is applied. An example shown in Figure 15 shows that the volume of consolidation settlement V cons when applying lower boundary of support pressure is smaller than the value of V cons when applying the maximum support pressure. When the tunnel is located at a deeper level, the volume loss V L,c when applying the maximum support pressure is higher than the volume loss V L,c in the case of applying minimum support pressure. The dependence of the volume loss due to consolidation V L,c on soil type is shown in Figure 17 for tunnels with diameters D = 6, 8, and 1m. It can be seen that the volume loss V L,c in the case of tunnelling in peat is much higher compared to tunnelling in clay and organic clay. 24

26 1.4 volume loss at the tail V L,c (%) upper boundary for D=1m lower boundary for D=1m upper boundary for D=8m lower boundary for D=8m upper boundary for D=6m lower boundary for D=6m (a) in clay 2.5 volume loss at the tail V L,c (%) lower boundary for D=1m lower boundary for D=8m upper boundary for D=8m lower boundary for D=6m upper boundary for D=6m upper boundary for D=1m (b) in organic clay 35 3 lower boundary for D=1m 25 lower boundary for D=8m volume loss at the tail V L,c (%) upper boundary for D=8m 5 lower boundary for D=6m upper boundary for D=1m upper boundary for D=6m (c) in peat Figure 16: Volume loss due to consolidation V L,c with different tunnel diameters in various 25 soils

27 16 5 volume loss at the tail V L,c (%) in clay in organic clay in peat in clay in organic clay in peat (a) D = 6m (b) D = 6m (in detailed) 25 5 volume loss at the tail V L,c (%) in clay in organic clay in peat in clay in organic clay in peat (c) D = 8m (d) D = 8m (in detailed) 35 5 volume loss at the tail V L,c (%) in clay in organic clay in peat in clay in organic clay in peat (e) D = 1m (f) D = 1m (in detailed) Figure 17: Volume loss due to consolidation V L,c with various tunnel diameters 26

28 5. Total volume loss and case studies 5.1. Total volume loss From Equation 1, the total volume loss is derived by summing the volume loss of tunnelling face, along the shield, at the tail and due to consolidation. Figures 18 and 19 show the total volume loss in the case of shallow tunnelling in sand and clayey sand. It can be seen that the range of the total volume loss decreases with the increase of the ratio and the tunnel diameter D. In the case of a ratio from.4 to 1, a volume loss in shallow tunnelling of less than.5% can be achieved with the condition of careful monitoring. The highest expected volume loss in this range of the ratio is about 3.7% for tunnelling in sand and 5% for tunnelling in clayey sand when less optimal but still stable support and grout pressures are applied. When the ratio larger than 1, the maximum volume loss is less than 1.5% with the range of support pressures in this study. These figures also show that a result of no volume loss can be achieved when tunnelling with 2. Figure 2 shows the relationship between the total volume loss V L and total volume loss(%) lower boundary for D=6m lower boundary for D=8m lower boundary for D=1m 1.5 upper boundary for D=6m upper boundary for D=8m upper boundary for D=1m Figure 18: Total volume loss for tunnelling in sand with various diameter D 27

29 5.5 5 total volume loss(%) lower boundary for D=1m lower boundary for D=8m lower boundary for D=6m upper boundary for D=6m upper boundary for D=8m upper boundary for D=1m Figure 19: Total volume loss for tunnelling in clayey sand with various diameter D ratios when tunnelling in clay. The total volume loss V L when tunnelling has just finished (not taking into account the consolidation) is shown in Figure 2a. It can be seen that a total volume loss V L less than.5% after tunnelling is feasible even with 1. This figure also shows that for very shallow tunnelling with.6, a tunnel with a large diameter has a larger range of expected volume loss. With deeper tunnelling when 1 2, the maximum value of the total volume loss reduces and becomes less than 2%. Figure 2b shows the total volume loss V L for tunnelling with various diameters D = 6,8 and 1m in clay including consolidation of soil layers above the tunnel. It also follows that the lower the ratio is, the larger the range of volume loss is. The total volume loss of tunnelling in clay would be at maximum about 6% with D = 1m, 5.5% with D = 8m and 5% with D = 6m when =.4. The lower boundary corresponding with the minimum support pressure applied has a reducing trend when the ratio increases. This means there might be a larger volume loss when the tunnel becomes shallower. At the upper boundary of the total volume loss, corresponding with the maximum support pressure applied, the final volume loss 28

30 of tunnelling with D = 6m can reach just over % after consolidation has been taken into account Case studies Tables 3 and 4 show a summary of volume loss in case histories for tunnelling in sand and clay with the ratio less than 2. In Table 3, there are two case studies the Ayshire Joint Drainage Scheme and WNTDC Lumb Brook Sewer, derived from the study of O Reilly and New (1982) and three case studies Second Heinenoord Tunnel, Botlek Railway Tunnel and Sophia Railway Tunnel, derived from Netzel (29). Table 4 shows volume loss data from various projects all over the world including Madrid Metro Extension, Heathrow Express Trail Tunnel, Waterloo, Garrison Dam test tunnel, Baulos 25, Barcelona Subway and London Transport Experimental Tunnel. The data in Table 3 are plotted in Figure 21a in order to compare the volume loss derived from this analysis to field data in the case of shallow to medium deep tunnelling in sand. It can be seen that most of field data falls in the boundaries of volume loss in Figure 21a. Figure21bshowsthevalidationofthecalculatedvolumelosstothefielddata in Table 4 in the case of shallow tunnelling in clay. It also shows that all the field data is in agreement with the boundaries of volume loss derived in this study. Only the Madrid Metro Extension which is known as a successful tunnelling project has one data point below the lower boundary of volume loss for a low ratio. The agreement between derived boundaries of volume loss in sand and clay and field data shows that the approach of estimating volume loss in this study can successfully predict volume loss in shallow tunnelling. 29

31 total volume loss(%) lower boundary for D=1m lower boundary for D=8m lower boundary for D=6m upper boundary for D=6m upper boundary for D=8m upper boundary for D=1m (a) not including consolidation of soil layers above the tunnel 6 total volume loss(%) lower boundary for D=1m lower boundary for D=8m lower boundary for D=6m 1 upper boundary for D=6m upper boundary for D=8m upper boundary for D=1m (b) including consolidation of soil layers above the tunnel Figure 2: Total volume loss for tunnelling in clay with various diameter D 3

32 Ayrshire Joint Drainage Scheme WNTDC Lumb Brook Sewer Second Heinenoord Tunnel Botlek Railway Tunnel Sophia Railway Tunnel total volume loss(%) upper boundary for D=6m lower boundary for D=6m lower boundary for D=8m lower boundary for D=1m 1 upper boundary for D=8m.5 upper boundary for D=1m (a) in sand total volume loss(%) upper boundary for D=6m lower boundary for D=1m upper boundary for D=8m Madrid Metro Extension Heathrow Express Trail Tunnel Waterloo Garrison Dam test tunnel Baulos 25 Barcelona Subway London Transport experimental tunnel lower boundary for D=8m lower boundary for D=6m 1.5 upper boundary for D=1m (b) in clay Figure 21: Validation for volume loss in shallow tunnelling 31

33 32 Case H(m) C(m) D(m) V L (%)Soil conditions Reference Ayrshire Joint Drainage Scheme fine to medium slightly silty sand; loose and medium density O Reilly and New (1982) WNTDC Lumb Brook Sewer medium/dense sands and O Reilly and New gravel with a little clay (1982) Second Heinenoord Tunnel dense sand Netzel (29) Botlek Railway Tunnel Holoceen and Pleistoceen sand Netzel (29) Sophia Railway Tunnel Pleistoceen sand Netzel (29) Table 3: Volume loss of tunnelling in sand projects

34 33 Case H(m) C(m) D(m) V L (%)Soil conditions Reference Madrid Metro Extension stiff clay Melis et al. (22) Heathrow Express Trail Tunnel London clay Bowers et al. (1996) Waterloo London clay Harris et al. (1994) Garrison Dam Test Tunnel clay-shale and lignite Peck (1969) Baulos Frankfurt clay Macklin (1999) Barcelona Subway red and brown clay with some Ledesma and Romero London Transport Experimental Tunnel gravel dense sandy gravel overlain with made ground of soft clay with sand and gravel (1997) O Reilly and New (1982) Table 4: Volume loss of tunnelling in clay projects

35 6. Conclusion Volume loss is a major parameter in the calculation of ground movement by tunnelling. The range of attainable volume loss can be estimated by combining stability analysis at tunnelling face, along and behind the shield. In thistheoreticalstudy, itisfoundthatinthecaseoftunnellingwith 1, the volume loss at the tunnelling face has a major impact in total volume loss. The volume loss along the shield can be optimized by selecting optimal bentonite and grout pressures applied at tunnelling face and tail. The proposed calculation method estimates attainable upper and lower boundaries of volume loss along the shield for a particular tunnel. This paper also presents methods to identify the volume loss behind the shield. The volume loss at the tail when tunnelling in peat has a large impact, especially in the case of shallow tunnels ( 1). The volume loss behind the tail was estimated by the volume loss due to shrinkage of grout and consolidation of above soil volume. The volume loss due to consolidation depends on the surrounding soil and the ratio. The total volume losses for tunnelling in sand, clayey sand and clay are derived and have a good agreement with case studies. Overall, the range of volume loss increases when tunnelling with shallower overburden. By controlling the applied support pressure at the tunnelling face and tail, the volume loss can be minimized. Still, a direct volume loss around 1% is a reasonable minimum for very shallow tunnels ( =.4) where for deeper tunnels no volume loss should be attainable. If pressure control is less optimal but still controlled, a direct volume loss up to 5.5% is not unreasonable to expect for very shallow tunnels. Analysis also shows that consolidation after the TBM has passed can contribute considerably to the final surface settlements and can be of this same order as direct volume loss effects in clay and even larger in very soft soils like peat. This effect, however, is more pronounced in deeper tunnels, where it could easily double the direct volume loss. References Anagnostou, G., Kovári, K., The face stability of slurry-shield-driven tunnels. Tunnelling and Underground Space Technology 9 (2),

36 Attewell, P., Farmer, I., Ground disturbance caused by shield tunnelling in a stiff, overconsolidated clay. Engineering Geology 8 (4), Attewell, P. B., Yeates, J., Selby, A. R., Soil movements induced by tunnelling and their effects on pipelines and structures. Methuen, Inc., New York, NY. Bezuijen, A., 27. Bentonite and grout flow around a tbm. In: Underground Space The 4th Dimension of Metropolises, Three Volume Set+ CD-ROM: Proceedings of the World Tunnel Congress 27 and 33rd ITA/AITES Annual General Assembly, Prague, May 27. CRC Press, p Bezuijen, A., Talmon, A., 28. Processes around a tbm. In: Proceedings of the 6th International Symposium on Geotechnical Aspects of Underground Construction in Soft Ground (IS-Shanghai 28). pp Bezuijen, A., Talmon, A., 213. Grout properties and their influence on back fill grouting. In: Proceedings of the 5th International Symposium on Geotechnical Aspects of Underground Construction in Soft Ground. p Bosch, J. W., Broere, W., 29. Small incidents, big consequences. leakage of a building pit causes major settlement of adjacent historical houses. amsterdam north-south metro line project. In: Safe Tunnelling for the City and Environment: ITA-AITES World Tunnel Congress, May 29, Budapest, Hungary. Hungarian Tunneling Association. Bowers, K., Hiller, D., New, B., Ground movement over three years at the heathrow express trial tunnel. In: Proceedings of the international symposium on geotechnical aspects of underground construction in soft ground, London. pp Broere, W., 21. Tunnel face stability & new cpt applications. Ph.D. thesis, Delft University of Technology. Broms, B. B., Bennermark, H., Stability of clay at vertical openings. Journal of Soil Mechanics & Foundations Div. 35

37 Cording, E. J., Hansmire, W., Displacements around soft ground tunnels. In: 5th Pan American Congress on Soil Mechanics and Foundation Engineering. Vol. 4. Buenos Aires, pp Dimmock, P. S., Mair, R. J., 27. Estimating volume loss for open-face tunnels in london clay. Proceedings of the ICE-Geotechnical Engineering 16 (1), Festa, D., Broere, W., Bosch, J., 215. Kinematic behaviour of a tunnel boring machine in soft soil: Theory and observations. Tunnelling and Underground Space Technology 49, Gemeente-Amsterdam, 29. Geotechnical base report alternative design deep stations and risk assessment bored tunnel. Tech. rep., Projectbureau Noord/Zuidlijn. Harris, D., Mair, R., Love, J., Taylor, R., Henderson, T., Observations of ground and structure movements for compensation grouting during tunnel construction at waterloo station. Geotechnique 44 (4), Jancsecz, S., Steiner, W., Face support for a large mix-shield in heterogeneous ground conditions. In: Tunnelling94. Springer, pp Ledesma, A., Romero, E., Systematic backanalysis in tunnel excavation problems as a monitoring technique. In: Proceedings of the International Conference on Soil Mechanics and Foundation Engineering- International Society for Soil Mechanics and Foundation Engineering. Vol. 3. AA BALKEMA, pp Loganathan, N., 211. An innovative method for assessing tunnelling-induced risks to adjacent structures. New York: Parsons Brinckerhoff Inc. One Penn Plaza New York. Macklin, S., The prediction of volume loss due to tunnelling in overconsolidated clay based on heading geometry and stability number. Ground engineering 32 (4), Mair, R., Discussion leaders report on session 9: Selection of design parameters for underground construction. In: Proceedings of the 12th International Conference on Soil Mechanics and Foundation Engineering, Rio de. Vol. 5. pp