Abstract

Peck method and stochastic medium method are the two most commonly used methods to estimate surface settlement caused by tunnel excavation. However, the Peck method was not suitable for a shallow-buried tunnel, and the calculation process of the stochastic medium theory was complicated. To solve this problem, in this paper, a simple and accurate prediction approach for surface settlement was obtained by improving the Peck method based on the basic idea of stochastic medium theory. In detail, the over-excavation area of the tunnel was divided into n independent units, and the surface settlement caused by the collapse of each unit was calculated, respectively. Then, the total surface settlement can be obtained by superimposing surface settlement induced by each unit. Taking the shallow-buried section of Mulingguan tunnel entrance as a case, the surface settlement calculated by the modified Peck formula and original Peck formula was compared with the observed data, respectively. The comparison results indicated that the surface settlement calculated by the modified Peck formula is closer to the observed data than that calculated by the original Peck formula in the calculation process of surface settlement of shallow-buried tunnel. The table of recommendation for the number of units can be obtained by a discussion of reasonable n values. Finally, the difference between the original Peck formula and the modified Peck formula was analysed, and the results showed that the change rule of the surface settlement is consistent under the tunnel depth, internal friction angle, and ground loss of the tunnel. However, the calculation error of the surface settlements calculated by the original Peck formula is more significant than that calculated by modified Peck formula under the tunnel diameter ratio being less than 1.75. The modified Peck formula is more suitable for calculating the surface settlement under internal angle friction being less than 20° or greater than 40°. The research results expand the scope of application of the original Peck formula and enrich the calculation approach of surface settlement induced by underground excavation in tunnel construction.

1. Introduction

Subway tunnels and other public facility tunnels are often constructed by the subsurface excavation method due to the problem of traffic congestion and the environmental quality requirements. The underground excavation of urban subway tunnels is bound to cause surface settlement and even cause surface collapse, thus affecting the safety of driving and pedestrians [18]. Therefore, it is of great practical significance to predict surface settlement induced by the underground excavation of urban subway tunnels to prevent excessive surface settlement or surface collapse.

At present, there are many methods for predicting surface settlement induced by tunnel construction, for example, the empirical formula method [912], elastic strain method [1316], the airy stress function method [1721], stochastic medium theory [2224], numerical simulation method [25, 26], and centrifuge test method [27, 28]. According to the above literature review, although there are many methods to predict surface settlement, the Peck method and stochastic medium theory are still the two most widely applied approaches to predict surface settlement. Therefore, many scholars modified and extended the Peck formula and stochastic medium theory. For example, according to the research of Jacobsz et al. [29] and Vorster et al. [30], the width of the settlement trough is related to the nature of soil. It indicated that measured data conform to the normal distribution law and its correlation coefficient was above 0.8 by using the Peck formula to fit the results of up to 51 settlement curves in different parts of the world [31]. Han et al. evaluated the applicability of the Peck formula in different regions by discussed and analysed surface settlement observed data collected from more than 30 groups of subway construction in China; the results showed that the Peck method was not suitable for shallow-buried tunnels [32].

Another method to predict surface settlement and deformation induced by tunnel construction is the stochastic medium theory. The stochastic medium theory was proposed by the Polish scholar Litwiniszyn to study the influence of coal seam mining on surface displacement. The rock and soil are regarded as a stochastic medium, and surface displacement caused by soil excavation is considered to be a stochastic process. Subsequently, the stochastic medium theory was introduced by Chinese scholars into the study of predicting surface settlement caused by urban subway construction [33, 34]. Compared with the Peck method, stochastic medium theory can be applied to a shallow-buried tunnel and various excavation sections, which has more advantages [35, 36]. However, its calculation method is complicated, so it is difficult to be widely used in engineering.

With respect to the problems as mentioned earlier, a simple and accurate prediction method for surface settlement was presented by improving the Peck formula based on the basic idea of stochastic medium theory in this paper. The shallow-buried section of the Mulingguan tunnel entrance is taken as an actual case to verify the rationality of the modified Peck formula. Furthermore, the calculation precision of the various number of units is analysed to obtain a reasonable number of division units. Finally, the difference between the modified Peck formula and original Peck formula is discussed, which further verified the rationality of the modified Peck formula.

2. Stochastic Medium Theory and Peck Formula

2.1. Stochastic Medium Theory

Litwiniszyn proposed the stochastic medium theory in 1957. He considered that surface settlement caused by the complete excavation of a subsurface unit is a stochastic event. In the plane, the excavated unit dξ dη at the subsurface depth of H completely collapses, and the settlement is [22]where ξ and η are the horizontal and vertical coordinates of the collapsed unit in the plane, respectively; is the influence radius; and , β is the influence angle of overlying rock and soil layer.

In tunnel engineering, the actual volume of the excavated soil is usually more massive than the design limit of the tunnel due to the limitation of construction technology and geological conditions. This part of the over-excavation continues to sink with the progress of construction until the lining of the tunnel is closed and formed. In the stochastic medium theory, the soil mass of the over-excavation part of the tunnel can be regarded as complete collapse, as shown in Figure 1. If to integrate part over-excavation, then we can get the surface settlement expressions:

It was difficult for the engineers to calculate the surface settlement due to the above integrand being nonintegrable.

2.2. Peck Formula

In 1969, Peck proposed that the ground settlement conforms to the law of normal distribution by summarising the measured data of ground settlement caused by a large number of tunnel construction at that time [9]. The following formula can be used to calculate surface settlement:where is the surface settlement value at a distance of x from the central axis of the tunnel; Smax is the maximum value of surface settlement, generally located in the axis of the tunnel; and i is the width of the surface settlement trough and is the horizontal distance from the axis of the tunnel to the inflexion point of the settlement curve (as shown in Figure 1). The Peck method was not suitable for the shallow-buried tunnel according to the analysis of the applicability of the Peck formula in different regions [32].

3. Modified Peck Formula Based on Stochastic Medium Theory

3.1. Correlation between Stochastic Medium Theory and Peck Formula

Assume that the surface settlement calculated by the Peck method for the same project is equal to that calculated by stochastic medium theory. According to the concept of ground loss in the Peck method, the lost volume of soil (Vl) is equal to the volume enclosed by the surface settlement curve and the original surface level:

The mathematical meaning of Vl is the area enclosed by the normal distribution curve and the horizontal coordinate axis. According to equation (5), this area is equal to the triangle area with height Smax and base (as shown in Figure 1). Therefore, the following correlation can be established between and i:

Equation (3) can be rewritten as

It can be seen that the Peck formula regards the overall over-excavation part (Vl) of the tunnel as a unit by comparing equations (6) and (1). The surface settlement caused by excavation conforms to the normal distribution law. Therefore, it can be considered that the Peck formula is an approximate calculation of stochastic medium theory, and the difference lies in that stochastic medium theory regards the over-excavation area as countless tiny units, while the Peck formula regards the whole over-excavation part as a tiny unit. It is assumed that the over-excavation area of the tunnel can be considered as a tiny unit in comparison with the tunnel depth and the area affected by the surface settlement caused by construction. Then the Peck formula is an approximate method of calculating the surface settlement of stochastic medium theory under the condition of large depth. The calculation accuracy of the Peck formula will increase with the increase of the tunnel depth of the tunnel and the influence area of surface settlement. With the decrease of the tunnel depth and the affected area, the accuracy of the surface settlement value calculated by the Peck formula becomes worse.

3.2. Modified Peck Formula

It was found that the Peck formula refers to treating the over-excavation part of the tunnel as a tiny unit by comparing equations (6) and (1), and the surface settlement caused by its complete collapse conforms to the normal distribution law, so there is an error in calculating the surface settlement. Equation (2) is the theoretical solution of surface settlement based on the stochastic medium theory. However, the solution is extremely complex and difficult to be widely applied in the engineering field due to the nonintegrability of its original function.

It is assumed that the over-excavation part of the tunnel is centred on the centre of the tunnel circle and evenly distributed on the ring around the outer diameter of the tunnel [34], as shown in Figure 2. The over-excavation part of the tunnel is divided into n tiny units with equal area. It is assumed that the collapse between tiny units is independent of each other, and surface settlement caused by the collapse of each tiny unit is calculated, respectively. Based on the principle of displacement superposition, the total settlement induced by tunnel excavation is equal to the sum of surface settlement caused by n independent units. This method for calculating surface settlement induced is called the modified Peck method in this paper. The more the number of the unit, the more accurate the surface settlement obtained by this method will be.

According to equations (1) and (3), the surface settlement induced by the i-th unit can be expressed aswhere R is the outer diameter of the tunnel, θi is the angle between i unit and the vertical axis of the tunnel, and ii is the width of the settlement trough caused by i units.

The width of the settlement trough caused by tunnel construction can be expressed aswhere zi is the depth of the unit and φ is the internal friction angle of the soil.

The expression of the total surface settlement curve is

In the actual calculation, the number n of over-excavation sections of the tunnel can be determined according to the tunnel depth and section size of the tunnel. The larger the n is, the more accurate the calculation result will be.

4. A Case Study

In this section, the engineering object of the shallow-buried section at the entrance of the Mulingguan tunnel [37] was used to analyse the variation law of the surface settlement curve under different n values to verify the modified Peck formula.

The depth of at the entrance of the Mulingguan tunnel is 5.1 m, the grade of surrounding rock is V grade, the ground loss rate is 0.45% [38], and the friction angle of soil is 30° (Code for Design of Railway Tunnel, China [39]). For the convenience of calculation, the original tunnel section is equivalent to a circular section according to the principle of equal area and buried depth. After #1 tunnel has been excavated for 50 m, #2 tunnel was excavated, and the three-step and seven-step method was adopted for construction. For the convenience of calculation, the equivalent radius of the tunnel section is taken to be 7.5 m, as shown in Figures 3(a) and 3(b). A total of 18 observing points are arranged above the tunnel at an interval of 5 m to study the deformation law of surface settlement, as shown in Figure 3(c). The final surface settlement curve of #1 tunnel after the excavation is shown in Figure 4(a). The maximum surface settlement value is 32.8 mm, and the actual width of the settlement trough obtained by fitting observing data is 9.0 m. The final surface settlement curve of #2 tunnel after the excavation is shown in Figure 4(b). The maximum surface settlement value of #1 tunnel is 34.2 mm, and the maximum surface settlement value of #2 tunnel is 37.9 mm.

According to equation (9), the settlement curve expressions of the tunnel surface at the time when the division number of the over-excavation part is 24, 20, 16, 12, 8, and 4 are calculated, respectively. The calculation results are shown in Table 1 and Figure 4.

According to Table 1, when the number of units increases from 4 to 8, the relative error of the width of the settlement trough decreased from 18.22% to 3.0%. If the number n continues to increase, the variation of relative error is less than 1%. When the number of partition units is greater than 16, the relative error is unchanged with the increase of n.

The shape of the surface settlement curve can be seen according to the number of different units given in Figure 4(a). When the number of units is n = 4, the settlement curve calculated by using equation (9) and the settlement curve when n = 8 differ greatly in the value of the maximum settlement, but the difference in the width of the settlement trough is not obvious. When the number of units n increases from 8, there is no significant change in either the width of the settlement trough or the value of the maximum settlement, and the settlement curve is nearly coincident.

Based on the analysis results of Table 1 and Figure 4(a), when the number n = 4, the calculated maximum settlement value is smaller than the measured maximum settlement value 9.75%, and when n is greater than or equal to 8, the calculated maximum settlement is 3.54%∼5.21% larger than the observed maximum settlement. The maximum surface settlement calculated with the Peck formula was 9.15% larger than the observed one, which indicates that the surface settlement calculated by the modified Peck formula is closer to the measured data and that the surface settlement calculated by the Peck formula has a significant error in the shallow-buried tunnel.

The surface settlement curve does not change when the n value continues to increase. Therefore, for this project, when n = 8 for the division of the over-excavation area, the calculation precision of the surface settlement curve obtained by using equation (9) can meet the accuracy of engineering estimation. Therefore, the engineering example also verified the rationality of the modified Peck formula based on the idea of stochastic medium theory, and there is a reasonable number of n when the modified Peck formula is used to calculate surface settlement.

As can be seen from Figure 4(b), the surface settlement curve calculated by the modified Peck formula is consistent with the observed surface settlement curve after the excavation of #2 tunnel. Affected by the construction of #2 tunnel, the maximum observed settlement value of #1 tunnel does not appear above the tunnel crown, but on the right side of the tunnel crown. However, the maximum calculated value of #1 tunnel is equal to the observed settlement. The observed surface settlement curve of #2 tunnel is consistent with the calculated settlement curve, but the maximum observed value of #2 tunnel is 9.8% larger than the calculated value. It shows that the modified Peck formula can better predict the surface settlement of single-line tunnel, but there will be some errors for the double-line tunnel.

5. Analysis and Discussion

5.1. Reasonable n Values

According to the above case study, the higher the n value is, the higher the calculation accuracy will be. However, the amount of calculation and the difficulty of calculation increase with the increase of n value. The reasonable unit number n is further discussed to meet the needs of engineering evaluation and to achieve the purpose of simple calculation. In this paper, the reasonable value of n is presented in the clay and sand soil layers commonly encountered in the construction process of urban subway tunnels. The variation range of ground loss rate of clay is 0.19%∼6.9%, that of sand is 0.08%∼3.01% [34], and the friction angle of soil is 30°. If the relative error = 1%, the control condition of reasonably divided number n iswhere Sn·max and are the maximum settlement values when the number of units is divided into n and n + 1, respectively.

The variation trend of the relative error with the number of n is shown in Figure 5.

As can be seen from Figure 5, when the tunnel diameter ratio (the ratio of tunnel depth to the equivalent diameter of the tunnel Z/D) is constant, the relative error does not change with the diameter but decreases with the increase of the tunnel diameter ratio. It is also found that the relative error does not change with the formation loss. When the tunnel diameter ratio is greater than or equal to 2, the relative error is reduced to less than 1%; only the over-excavation part needs to be divided into 4 units, which indicates the reasonable units n = 4 under this condition. When the tunnel diameter ratio is less than 2, to reduce the relative error to less than 1%, the over-excavation part needs to be divided into 8 units, indicating the reasonable division unit number n = 8 under this condition.

The urban subway tunnel passes through various soil layers in the construction process. Different φ values represent various soil layers; the obtained relative error δ with the n change trend is shown in Figure 6.

5.2. Comparison of the Modified Peck Formula and Peck Formula

To further analyse and verify the rationality of the modified Peck formula, we took urban subway construction as the engineering background and compared the differences of the surface settlement curves calculated by the modified Peck formula and Peck formula under different tunnel depth, soil internal friction angle, and ground loss rate conditions.

5.2.1. Comparison of Surface Settlement Curves with Different Tunnel Depths

To further analyse and verify the rationality and feasibility of the modified Peck formula, the paper took the urban subway construction in the soft and clay soil layer as the engineering background and compared and analysed the differences of the surface settlement curves calculated by the modified Peck formula and Peck formula under different tunnel buried depth conditions. Based on the comparative analysis and discussion, the applicable range of the Peck formula for engineering allowable error range is discussed.

Tunnel diameter is 6.0 m, the ground loss rate is 3.0%, and the internal friction angle of soil is 12.5°. The over-excavation part of the tunnel is divided into 8 units, and the surface settlement curves at the tunnel depth of 6.0 m, 7.5 m, 9.0 m, 10.5 m, and 12.0 m are calculated, respectively. Table 2 shows the modified Peck formula and Peck formula to calculate the parameters of the surface settlement curve at different tunnel depths. Figures 7(a)7(e), respectively, show the settlement curves obtained by two calculation formulas at different tunnel depths.

According to Table 3 and Figure 7, the maximum settlement calculated by the Peck method is larger than that obtained by the modified Peck formula, but the width of the settlement trough calculated by the modified Peck formula was smaller than that obtained by the Peck formula. The settlement curve calculated by the Peck formula is broader and shallower than that achieved by the modified formula. With the increase of the tunnel depth, the surface settlement curve tends to be broad and shallow. As for the width of the settlement trough, the relative error between Peck formula and the modified formula is 18.0% when the tunnel depth was 6.0 m. When the tunnel depth is 9.0 m, the relative error is reduced to 6.16%. When the tunnel depth increases to 12 m, the relative error decreases to 3.35%.

Therefore, the curve calculated by the Peck formula is more and more similar to the settlement curve calculated by the modified formula, and the value of the maximum settlement and the width of the settlement trough are also more and more similar. However, the surface settlement calculated by the Peck formula has a large error in the tunnel diameter ratio less than 1.75, which also shows that the Peck method is not suitable for the shallow-buried tunnel. These results are similar to those of Han et al. [32].

5.2.2. Comparison of Surface Settlement Curves with Different Internal Friction Angles

The tunnel diameter is 6.0 m, the ground loss rate is set at 3.0%, and the tunnel depth is 6.0 m. The over-excavation part of the tunnel is divided into 8 units, and the surface settlement curves of the soils with friction angles of 0°, 10°, 2°, 30°, 40°, 50°, and 60° are calculated, respectively. Table 4 shows the modified formula and the Peck formula to estimate the parameters of the surface settlement curve at different internal friction angles. Figures 8(a)8(f), respectively, shows the comparison of settlement curves obtained by the two formulas at different internal friction angles.

According to Table 3 and Figure 8, the maximum settlement value and settlement trough width calculated by the Peck formula and the modified Peck formula both increase with the increase of internal friction angle. When the internal friction angle 0°∼30°, the maximum settlement and the width of the settlement trough calculated by the modified Peck formula are smaller than the Peck formula, and the relative error of the two decreases with the increase of the internal friction angle. When the internal friction angle is 30°, the relative error between them is the smallest. When the internal friction angle is 30°∼50°, the maximum settlement and the width of the settlement trough calculated by the modified Peck formula are larger than the Peck formula, and the relative error between them increases with the increase of the internal friction angle.

According to the above analysis, when the internal angle friction is about 30°, the settlement curve obtained by the Peck formula and the modified Peck formula is equal. It also shows that the modified Peck formula is more sensitive to the friction angle. Therefore, when the internal angle friction is about 30°, this paper suggests using the Peck formula or the modified Peck formula to calculate the surface settlement. When internal angle friction is less than 20° or greater than 40°, the modified Peck formula is adopted to calculate the surface settlement.

5.2.3. Comparison of Surface Settlement Curves with a Different Ground Loss Rates

The tunnel diameter is 6.0 m, the internal friction angle is 12.5°, and the tunnel depth is 6.0 m. The over-excavation part of the tunnel is divided into 8 units, and the surface settlement curves when the ground loss rate is 1.0%, 2.0%, 3.0%, 4.0%, 5.0%, and 6.0% are calculated, respectively. Table 5 shows the modified formula and the Peck formula to estimate the parameters of the surface settlement curve under different bottom loss rates. Figures 9(a)9(f), respectively, shows the comparison of settlement curves under different ground loss rates under two calculation formulas.

According to Table 5 and Figure 9, the maximum settlement calculated by the modified Peck formula was smaller than that of the Peck formula. The maximum settlement calculated by the Peck formula and the modified Peck formula increases with the increase of ground loss rate, but the width of the settlement trough hardly changed. The relative error between the maximum settlement and the width of the settlement trough calculated by the two methods hardly changed with the increase of formation loss rate. It shows that when the tunnel diameter ratio and the internal friction angle of the soil are constant, the relative errors of the two methods are also constant, and the width of the settlement trough will not change the settlement value.

6. Conclusions

The Peck method was not suitable for shallow-buried tunnels, and the calculation process of the stochastic medium theory was complicated. To address this problem, in this study, a simple and accurate prediction method for surface settlement was obtained by improving Peck method based on the basic idea of stochastic medium theory. The modified Peck formula was verified by an engineering example. The obtained conclusions are summarised as follows:(1)The Peck formula is an approximate algorithm of the stochastic medium theory under the condition of large depth and a small section of the tunnel by comparing the calculation principle of surface settlement with the stochastic medium theory and the Peck formula.(2)The case study indicated that the surface settlement calculated by the modified Peck formula is closer to the measured data than that calculated by the original Peck formula. It also showed that the modified Peck formula could better predict the surface settlement of single-line tunnel, but there will be some errors for the double-line tunnel.(3)The reasonable number of n does not change with the change of tunnel diameter when the tunnel diameter ratio is constant. With the control condition of the relative error is equal to 1%, the reasonable number of n with the tunnel diameter greater than or equal to 2 for the soil layer of whatever nature is equal to 4. When the tunnel diameter ratio is less than 2, the internal friction angle is greater than or equal to 0° and less than or equal to 10°; the reasonable unit number n equals 8. Meanwhile, internal friction angle φ is more than 10°; the reasonable unit number is equal to 10.(4)The change rule of surface settlement is consistent under the tunnel depth, internal friction angle, and ground loss of the tunnel through comparing the surface settlement curve calculated by the modified Peck formula and Peck formula. However, the calculation error of the surface settlements calculated by the original Peck formula is larger than the surface settlements calculated by the modified Peck formula under the tunnel diameter ratio being less than 1.75. The modified Peck formula is more suitable for calculating the surface settlement under internal angle friction being less than 20° or greater than 40°.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was subsidized and supported by the National Natural Science Foundation of China (No. 51578447), the Science and Technology Project of Ministry of Housing Urban-Rural Construction (No. 2017-K4-032), and the project funded by China Postdoctoral Science Foundation (No. 2018M643809XB). The financial supports are gratefully acknowledged by the authors.