Abstract

In order to accurately and efficiently assess the impact of tunnel excavation on overlying existing pipeline, an analytical method is proposed to solve this problem. First, the vertical free displacement of the surrounding soil due to tunnel excavation can be derived by the Loganathan formula. Next, the overlying existing pipeline can be treated as a Timoshenko beam resting on the Vlasov foundation model, and the influence of the surrounding soil on the both sides of the existing pipeline is taken into consideration. Finally, an analytical solution for the longitudinal deformation of the existing pipeline can be obtained by using the integral method. Case analysis results demonstrate that the calculated results of this method closely in line with measured data. Compared to the degenerate analytical solution given by this method, the result from this method is more consistent with the measured data. Further parameter studies show that the volume loss rate, diameter of new tunnel, skew angle, and vertical distance between tunnel and pipeline are significant factors affecting the existing pipeline response due to tunneling underlying.

1. Introduction

In recent years, with the flourishing development of urban rail transit, the pressure on urban surface transportation has been significantly alleviated, making the development of urban subways a consensus among urban planners [19]. However, in the process of underground space development, existing subway lines and pipeline are inevitably threatened by nearby excavations, which can affect their safety and comfort, potentially leading to issues like tunnel or pipeline lining cracking, water seepage, and even structural damage. In many cases, shield tunnels crossing adjacent existing pipelines can pose challenges to ensuring the stability and safety of nearby pipelines, making it a hot topic in the field of geotechnical engineering. This issue has garnered increasing interest from scholars.

Compared to the labor-intensive and time-consuming methods of finite element analysis [1013] and model testing [1416], analytical methods [1728] can be served as a fast and lower cost means to estimate the underlying excavation-triggered pipeline deformation. Especially those based on the two-stage approach, offer a quicker and more accurate means of assessing the impact of tunnel excavations on existing tunnels or pipelines. These analytical methods have two stages: the first stage relies on classical formulas, often based on the work of Loganathan and others, to estimate the influence of tunnel excavations on surrounding soil. The second stage involves applying soil displacements in reverse to the existing tunnel structure, based on various beam and foundation models, to obtain the stress and deformation responses of the existing tunnel due to nearby tunnel excavations [2430].

However, Winkler foundation models [31], only consider the elastic behavior of the soil, and their calculated data often deviates significantly from field monitoring data. To address this limitation, scholars have proposed dual-parameter Pasternak [32] and Vlasov and Leontev [33] and three-parameters Kerr foundation models [4, 5, 3436] that consider soil shear deformations. However, the control equation of tunnel responses is complicated by the multiplicity of parameters in the Kerr foundation model. Its inherent complexity in parameter determination and it cannot be applied in actually engineering projects. At the same time, the Pasternak foundation model does not always offer an accurate calculation of subgrade parameters. It may often overlook disturbances induced by tunneling in the surrounding soil. Determining the shear modulus from the Pasternak and Kerr foundation models in real-world scenarios remains challenging. Consequently, many researchers have sought alternative means to equivalently substitute the second parameter, avoiding the difficulty of accurate determination. Instead, the Vlasov foundation model can be used to overcome many of these limitations [37]. Meanwhile, the metro tunnel and pipeline in soft areas are usually constructed by shield tunneling, and the tunnel is lined with prefabricated concrete segment linked by bolted joint. The Euler–Bernoulli beam, without taking the shear behavior into consideration, is not appropriate to estimate the mechanical of tunnel and pipeline response. To overcome this shortcoming, Timoshenko [38] beam was proposed and applied in an amount of investigations [28, 29].

In this paper, a Vlasov foundation model is selected with explicit parameters, which is less commonly used in practical engineering. Additionally, the lateral soil effects on pipeline structural deformation are taken into consideration, which always been overlooked. Then, the existing pipeline is treated as a Timoshenko beam placed on the Vlasov foundation model, accounting for the influence of lateral soil forces on the pipeline structure. By utilizing an integral method, an analytical solution is derived for the stress and deformation of the existing pipeline caused by tunnel excavations and a field verification is taking into discussion. Finally, an exhaustive parameter analysis, including volume loss rate, diameter of new tunnel, skew angle, and vertical distance between tunnel and pipeline, is conducted.

2. Analysis Method

To evaluate the interaction between tunnel excavation underneath, soil, and pipeline, the two-stage method is introduced here in this section. Some basic assumptions are as follows:(1)The soil movement along the existing pipeline can be computed by using the Loganathan and Poulos [39] equation.(2)The existing pipeline is assumed to be an infinitely long Timoshenko beam resting on an elastic foundation model.(3)The elastic foundation model employed is the Vlasov foundation model.(4)Deformations of the pipeline and soil are assumed to be mutually consistent, and cases involving the separation of the pipeline from the soil are not considered.(5)The spatial variability of soil properties and the effect of water pressure in this study are not taking into consideration.

2.1. Adjacent Tunnel Excavation-Triggered Unloading Stress

The free displacement of the surrounding soil due to tunnel excavation (as shown in Figure 1), as suggested by Loganathan and Poulos [39], is given by the following equation:where R represents the excavation radius, H is the depth of the tunnel axis, x is the horizontal distance of the centerline between the new tunnel and existing pipeline, z is the vertical distance from the ground surface, ε is the equivalent ground loss ratio, and υ is the Poisson’s ratio of the soil. Noted that when the new tunnel and existing pipeline have a nonzero inclination angle θ0 in the horizontal plane, Equation (1) is replaced as follows:

Consequently, the additional stress at the pipeline axis due to tunnel excavation is given by the following equation:

Merging Equations (1)–(3), q(x) can be easily solved.

2.2. Pipeline–Soil Interaction

According to the Valsov model (Figure 2), the foundation reaction beneath the pipeline p(x) can be expressed as follows:where w(x) denotes the pipeline deformation, k is the spring stiffness of the soil, and 2t is the shear stiffness of the foundation layer. The values of k and 2t can be determined as follows:where Es is the elastic modulus of the soil, υ is the Poisson’s ratio of the soil, h = h(y) is a function varying in the y-direction, and Hs is taken as 2.5 times the diameter of the pipeline.

Most previous research has ignored the lateral soil effects on pipeline–soil interaction. Therefore, the paper introduces the lateral soil effects besides the pipeline. The equilibrium of the pipeline element under these lateral forces is analyzed, as shown in Figure 3. According to existing literature [40], the lateral forces from the soil on both sides of the pipeline satisfy the following equation:

Based on Timoshenko’s [38] beam, the equilibrium equation is as follows:where κ is equivalent cross-section coefficient, pipeline is circular cross-section, which can be taken as 0.5, G is the tunnel shear modulus, A is the annular section area of the tunnel, and EI is the stiffness of the segmental rings.

Then, considering the vertical equilibrium of the existing pipeline:

Combing Equations (4)–(8), the vertical equilibrium of the existing pipeline can be expressed as follows:where

The homogeneous equation of Equation (9) is as follows:

Then, the solution of Equation (11) can be obtained as follows:where C1, C2, C3, and C4 can be expressed as follows:

Considering Equation (9), which is a fourth-order differential equation, the analytical solution is obtained using an approach that involves assuming an arbitrary force P acting at the x = 0 point of the infinitely long existing pipeline. The boundary of the pipeline in this case is defined as follows:

Substituting Equation (14) into Equation (12) yields the analytical solution for the deformation of the infinite existing pipeline under the action of any force P at x = 0:

Assuming that tunnel excavation underneath causes an additional stress q(x) at any point on the axis of the existing pipeline, the additional stress acting on the existing pipeline at any pipeline point can be calculated from Equation (15). This additional stress causes vertical displacement dw (x) at any point on the pipeline axis:

Integrating Equation (16) over the range of the distribution of the additional stress on the tunnel results in the longitudinal deformation of the existing pipeline caused by tunneling underlying:

Additionally, the analytical solutions for the flexural moments and shear forces in the existing pipeline can be derived as follows:

Therefore, the analytical solutions for the longitudinal deformation and internal forces of existing pipeline induced by tunneling underneath can be obtained. It is important to note that when there are no lateral soil effects on both sides (T1 = T2 = 0), the analytical results degenerate into the TV model (in which the pipeline is treated as Timoshenko beam resting on Vlasov foundation model), and when the lateral soil effects and soil shear stiffness are not considered, the results degenerate into the TW model (in which the pipeline is treated as Timoshenko beam resting on Winkler foundation model).

3. Verification

To validate the correctness of the methodology presented in this paper, the author collected measured data from Shenzhen metro and have a discussion between the results given by this study and its degenerate solutions. Therefore, a typical case, where the horizontal location between tunnel and pipeline is almost perpendicularity, is introduced here to verify the rationality of this study, as illustrated in Figure 4. The engineering parameters are as follows: the diameter and axis depth of existing pipeline are D = 3.0 and z0 = 8.7 m, respectively. The diameter and axis depth of new tunnel are 2R = 6.0 and H = 14.4 m, respectively. According to existing literature [41], the soil’s elastic modulus and Poisson’s ratio can be regarded as Es = 8.2 MPa and v = 0.3, the compression stiffness and shear modulus of this model can be calculated as 2,873 and 14,741 kN/m3, respectively. The excavation of the tunnel resulted in an average volume loss rate ε = 0.84%, and the existing pipeline’s flexural stiffness and shearing stiffness can be treated as 5.87 × 107 kN·m2 and 1.16 × 107 kN/m [42, 43].

In order to obtain a more accurate calculation result and eliminate influence of the size of length of existing pipeline, the calculated deformation of the pipeline is within 200 m from both sides of the neutral axis of the pipeline, and the calculation results of interception are within the range of 100 m in the middle. Therefore, the results obtained by using the proposed method are compared with engineering monitoring data in Figure 5, where the monitoring data are given by Zhang and Zhang [41]. Figure 5 indicates that the results obtained using our method and its degenerate solution are consistent with the trend of the measured data. The significant deformation of the existing pipeline is concentrated within a 30-m range on either side of the pipeline’s central axis, while the measured data indicates deformations greater than 2 mm distributed within a 20-m range on either side of the neutral axis of the pipeline. The data obtained from this paper and its degenerate solution are in better agreement with the distribution of the measured data. By using our method, the maximum displacement of the pipeline is calculated to be 7.6 mm, while the degenerate TV model yields a smaller displacement with a peak of 5.9 mm. The TW model produces significantly larger results with a maximum pipeline displacement of 10.9 mm. This phenomenon is because the TV model neglects the influence of lateral soil on the pipeline–soil interaction, leading to an underestimation of the tunnel’s impact on the existing pipeline. On the other hand, the Winkler foundation model does not consider the effect of soil shearing on soil–structure interaction, resulting in a significant deviation from the measured values. Compared to the degenerate solution of our method, the maximum pipeline deformation from the measured data [41] is 8.0 mm, indicating that the results obtained in this paper are more in line with the monitoring data.

4. Parameter Analysis

Herein, a series of parameter studies are conducted to examine influencing factors on the existing pipeline’s response to tunneling underneath, including the volume loss rate, diameter of new tunnel, the skew angle, and vertical distance between tunnel and pipeline. The typical engineering parameters are assumed as in Table 1.

4.1. Volume Loss Rate

Figure 6 shows the calculated results of the existing pipeline response triggered by tunneling underneath by using our method with different volume loss rate. The volume loss rate can be regarded as 0.5%, 1.0%, 1.5%, 2.0%, and 2.5%.

As shown in Figure 6, the longitudinal displacement and bending moment of the existing pipeline are symmetrically distributed along the central axis of the pipeline, and the maximum longitudinal displacement and maximum bending moment values occur at the pipeline’s central axis. Moreover, both the longitudinal displacement and bending moment of the existing pipeline linearly increase with an increase in volume loss rate. This is because the additional stress on the pipeline linearly increases with the volume loss rate, leading to linear variations in the longitudinal displacement at each position of the existing pipeline. Therefore, in practical engineering, it is advisable to minimize the volume loss rate caused by shield tunneling to reduce the adverse effects on nearby structures.

4.2. Diameter of New Tunnel

Figure 7 illustrates the existing pipeline’s response to a new tunnel excavation underlying by using our method with different diameters of new tunnel. Here, the diameter of new tunnel can be assumed as 3.0, 3.5, 4.0, 4.5, and 5.0 m.

Figure 7 demonstrates that the longitudinal displacement and bending moment of the existing pipeline are symmetrically distributed along the central axis of the pipeline. Meanwhile, as the diameter of new tunnel increases from 3.0 to 5.0 m, the maximum longitudinal displacement of the existing pipeline gradually decreases from 4.4 to 12.7 mm, with an increase rate of up to 1.9 times, and the growth rate tends to increase gradually. Moreover, with an increase in the diameter of new tunnel, the maximum bending moment of the pipeline increases from 1.4 to 3.9 MN·m2, with an increase rate of ∼1.8 times. This indicates that increasing the diameter of new tunnel has a significant impact on the overlying existing pipeline, and the stress and strain peaks on the pipeline increase rapidly. Therefore, in practical engineering, efforts should be made to minimize the diameter of new tunnel to reduce the impact of ground stress changes on the existing pipeline.

4.3. Skew Angle between Tunnel and Pipeline

Figure 8 shows the calculated results of deformations and internal forces given by the proposed method with different angles between tunnel and pipeline, which can be taken as θ = 15°, 30°, 45°, 60°, and 75°.

Figure 8 indicates that the deformation and bending moment of the existing pipeline are symmetrically distributed along the central axis of the pipeline. Meanwhile, as the angle between the new tunnel and the existing pipeline increases from 15° to 75°, the maximum longitudinal displacement of the existing pipeline gradually decreases from 8.1 to 4.6 mm, a decrease of ∼43%, and the rate of decrease becomes progressively smoother. However, it also shows that with an increase in the angle between the new tunnel and the existing pipeline, the maximum bending moment of the pipeline increases from 0.52 to 1.41 MN·m2, an increase of ∼1.7 times, and the growth rate becomes smoother as well. This is because, as the relative position of the new tunnel and the existing pipeline changes from a more “coincidental” state to a “vertical” state, the impact of ground settlement caused by the new tunnel on the existing pipeline axis decreases. Meanwhile, the additional stress on the existing pipeline changes from “uniform load” to “concentrated load” mode, resulting in an increase in the maximum additional stress on the pipeline. This indicates that the angle between the new tunnel and the existing pipeline is a sensitive parameter affecting the deformation and internal forces of the overlying existing pipeline. In practical engineering, attention should be paid to the influence of changes in the skew angle between the new tunnel and the existing pipeline on the existing pipeline response.

4.4. Vertical Distance between Tunnel and Pipeline

Figure 9 shows the calculated results of the existing pipeline response triggered by tunneling underneath by using our method with different vertical distances between tunnel and pipeline. The relationship d = nD is adhered to, where D represents the diameter of the pipeline. For this study, values of n such as 0.5, 1.0, 1.5, 2.0, and 2.5 are implemented.

As shown in Figure 9, the longitudinal displacement and bending moment of the existing pipeline are symmetrically distributed along the central axis of the pipeline, and the maximum longitudinal displacement and maximum bending moment values occur at the pipeline’s central axis. Meanwhile, as the vertical distance between tunnel and pipeline increases from 0.5D meters to 2.5D, the maximum longitudinal displacement of the existing pipeline gradually decreases from 6.0 to 2.8 mm, and the rate of decrease becomes progressively smoother. Moreover, with an increase in the vertical distance between tunnel and pipeline, the maximum bending moment of the pipeline decreases from 1.6 to 0.87 MN·m2, with an decrease rate of ∼46%. This indicates that increasing the vertical distance between tunnel and pipeline has a significant impact on the overlying existing pipeline, and it is useful for enlarging the vertical distance between tunnel and pipeline to ensure the existing pipeline safety.

5. Conclusions

Based on a two-stage approach, this paper presents a theoretical method for predicting the deformation response of overlying existing pipelines when shield tunnels are constructed beneath them. The conclusions obtained are as follows:(1)The pipeline is simplified as an infinitely long Timoshenko beam placed on Vlasov foundation model, and the influence of lateral soil on both sides of the pipeline is introduced. The integral method is used to analyze the stress and deformation response of the pipeline.(2)Compared to measured data from a specific project (Shenzhen Project), the results obtained by using our method are in good agreement with the measured data. Compared with the degenerate solutions (TV and TW methods), the suggested method is more in line with monitoring data.(3)Increasing the ground loss rate effectively increases the displacement and internal forces of the existing pipeline, but the displacement and internal forces of the existing pipeline decrease with the increase the vertical distance between tunnel and pipeline. The displacement and bending moment values of the pipeline gradually increase with an increase in the diameters of the tunnel passing beneath. Increasing the skew angle between the new tunnel and the existing pipeline leads to a reduction in the pipeline’s displacement but an increase in its internal forces, with the rate of change becoming smoother. Therefore, the skew angle between tunnel and the existing pipeline is a sensitive parameter that should be carefully considered in practical engineering to assess the deformation and stress response of the existing pipeline.(4)The suggested analytical solution is appropriate for beams that are supported by linear foundation soils. The absence of consideration for the nonlinearity of pipeline–ground interaction is evident in this context. Future investigations should pay closer attention to these deficiencies.

Data Availability

All data, models, or codes that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Authors’ Contributions

Yao Rong contributed to methodology, formal analysis, writing–original draft, and software. Guohui Feng contributed to conceptualization, writing–review and editing, and funding acquisition. Yang Sun contributed to formal analysis and writing–review and editing. Yujie Li contributed to resources and data curation. Guanyu Chen contributed to writing–review and editing and visualization. Haibin Ding contributed to writing–review and editing. Changjie Xu contributed to supervision and funding acquisition.

Acknowledgments

This research was funded by the National Key R&D Program of China, grant number 2023YFC3009400; the National Natural Science Foundation of China, grant numbers 52068033, 52238009; the National Science Fund of Jiangxi Province, grant numbers 20223BBG71018, 20232BCJ23069; and the foundation of the Department of Transportation of Jiangxi Province, grant numbers 2021C0006, 2022Z0001, 2022Z0002.