Influence of Underground Tunnel Excavation on Upper-level Existing Pipelines

Winkler elastic foundation beam theory is often used to observe the influence of underground engineering construction on pipelines. Excavation under existing pipelines disturbs the lower strata to varying extent, which creates significant differences in the subgrade coefficients. Ignoring these differences can severely impact the accuracy of calculation results. A theoretical model of the effects of underground excavation on upper-level existing pipelines was developed in this study based on Winkler elastic foundation beam theory and considering the difference of subgrade coefficients under the existing pipeline. Methods for determining the subgrade coefficient under the influence of construction disturbance and other relevant parameters are proposed. ,e theoretical model can be used to calculate the settlement of an existing pipeline under the influence of underground excavation. ,e calculated settlement was compared with settlement measurements to find that the subgrade coefficient of the excavation section is the most important parameter in the calculation process.,e suitable value range of the excavation subgrade coefficient is relatively large, as it has a relatively minor influence on the calculation result. When the value of subgrade coefficient of each section is appropriate, the calculated settlement curve is in close accordance with the measured settlement curve.


Introduction
e gradual development and utilization of urban underground spaces have produced a number of underground excavation projects affecting existing pipelines [1][2][3]. In the process of an underground engineering project, disturbance to the stratum inevitably deforms the adjacent underground pipeline. When the deformation exceeds the allowable deformation of the existing pipeline and its joints, destruction of the existing pipeline may result in even serious secondary disasters. It is difficult to directly monitor the deformation of existing pipelines due to the distribution of points above ground throughout a ground excavation project. e points instead can only be monitored indirectly [4], which does not readily allow for effective accident prediction or prevention.
ere is an urgent demand for a straightforward and easily applied method to evaluate the safety of existing pipelines and classify construction risks. Such evaluations would facilitate the proper selection of preconstruction reinforcement or relocation treatment for existing pipelines to minimize the failure risk of existing pipelines, save costs, and shorten the construction period.
Previous studies regarding the impact of tunnel construction on existing pipelines have mainly centered on numerical simulations and theoretical analyses. eoretical analysis methods are mainly based on the Winkler elastic foundation beam theory [5,6], two-parameter Pasternak elastic foundation beam theory [7,8], or energy method [9] to solve parameters such as the deformation, bending moment, and shear force of an existing pipeline as affected by tunnel construction and to select corresponding standards to judge its safety. Unlike the two-parameter elastic foundation beam theory, the Winkler elastic foundation beam theory does not reflect the continuity of foundation deformation. However, a method based on the Winkler elastic foundation beam is simpler, includes more accumulated experiences in parameter values, and is most widely used in engineering. e two-parameter elastic foundation beam is more comprehensive in theory but includes less experience in the parameter value selection process; whether the parameter value is appropriate is directly related to the accuracy of calculation results. e Winkler elastic foundation beam theory was utilized in the present study to develop a calculation model of the influence of excavation on the upper-level pipeline.
e stratum under an existing upper-level pipeline is disturbed as an excavation project passes under it. e disturbance is most intense above the excavation section though there is also disturbance in a certain area near the excavation area. e mechanical properties of the soils in these two areas change significantly at this point compared with the original stratum, so there is also a substantial difference in the subgrade coefficient. Previous theoretical calculations have rarely considered the differences in subgrade coefficient in various areas under the existing pipeline, though ignoring such differences drives down the accuracy of calculation results and subsequent pipeline safety assessments.
Based on Winkler elastic foundation beam theory and considering the different subgrade coefficients in each area under the existing pipeline, a theoretical calculation model of the influence of concealed excavation construction on upper existing pipelines was established in this study. Important parameters in the model were acquired, and the influence of different subgrade coefficient values on the predicted deformation of the existing pipeline was observed. A case analysis was conducted to test the proposed method.
is work may provide an important reference for the safety assessment of existing pipelines in similar projects.

Theoretical Calculation Model
Generally, underground pipelines are laid along a roadway. Most underground tunnels cross the road vertically and also intersect with the existing pipelines vertically. e impact of the underground excavation project intersecting with the existing pipeline vertically is under investigation here ( Figure 1). Flexible pipes are allowed to rotate after loading, while rigid pipes are not allowed to rotate after loading. e pipe considered in this paper allows rotation angle, so it is suitable for the calculation of flexible pipe.
It is also assumed here that the excavation channel is rectangular and that the foundation is homogeneous. e existing pipeline is regarded as an elastic foundation beam on a Winkle foundation. It is assumed that the settlement of any point in the foundation is only proportional to the pressure on the unit area of the point, regardless of the pressure at other points. ere is a relatively long longitudinal extension of the buried pipeline, which is assumed to be an infinite beam. During the underground excavation project, the stress area of the existing pipeline at the upper level is divided into five sections. As shown in Figure 2, the CD section is the excavation area, BC and DE sections are soil-loosening sections caused by excavation, and AB and EF sections are areas not affected by excavation. Construction disturbance creates different subgrade coefficients for each section. It is assumed that the subgrade coefficient of the AB and EF section is k 1 , that of the BC and DE section is k 2 , and that of the CD section is k 3 . e maximum load of the buried pipeline mainly originates in the gravity of the soil above it. e Terzaghi loose Earth pressure calculation model is generally used to determine the vertical Earth pressure on a buried pipeline [10]. In the Terzaghi model, the dislocation of soil particles produces stress transfer, which makes the surrounding soil restrict the downward moving soil so that the vertical Earth pressure borne by the underground structure is less than the gravity of the soil column. As shown in Figure 2, if other external loads on the upper part of the existing pipeline are not considered, the uniformly distributed load q 1 of sections AB and EF of the existing pipe jacking is not affected by excavation and can be calculated using the Terzaghi model. e soil above the excavation section and the excavation-affected section is impacted by construction disturbance. e mechanical properties of the soil and the holding capacity of the downward part are reduced, which increases the Earth pressure borne by the underground pipeline. It is assumed that the Earth pressure borne by the pipeline above the excavation section CD is q 3 and that of the excavation-affected sections BC and DE is q 2 .
As shown in Figure 2, the deflection differential equation of each section can be obtained as follows:

Road extension direction
Existing pipeline Underground excavation passageway Figure 1: Locations of the underground excavation project and existing pipeline. where EI is the bending stiffness of the existing pipeline, N•m 2 ; ω i is the deflection of each section of the existing pipeline, m, i � 1, 2, 3; x is the distance between each point in the right half of the existing pipeline and point O, m; k i is the subgrade coefficient of each section under the existing pipeline, N/m 3 , i � 1, 2, 3; q i is the Earth pressure borne by each section of the existing pipeline, Pa, i � 1, 2, 3; and b is the calculated width of the foundation reaction, m. e characteristic value of flexibility is determined by the calculated width of the foundation reaction.
where i � 1, 2, 3. e general solution for each section of the deflection differential equation and the first, second, and third derivative of each section can be obtained as follows: If the deflection of the beam is known, the rotation angle, bending moment, and shear force of any section of the beam can be obtained according to the knowledge of Material Mechanics.
where θ denotes the corners of existing pipelines; M is the bending moment of the existing pipeline, N•m; and Q is the shear force of the existing pipeline, N.
In an actual construction project, the pipeline is quite long. Due to the constraints of the soil around the pipeline on the existing pipeline, the excavation under the pipeline has a limited range of influence. At the end far away from the excavation range (x, ∞), the boundary condition can be assumed to be fixed support. erefore, the following boundary conditions and deformation coordination conditions hold.
Point A (x, ∞) is the fixed end: Point O is the symmetrical point: According to the coordination among deflection, rotation angle, bending moment, and shear force at the left and right ends of points B and C, the following supplementary conditions hold.
Point B: Point C: A system of 12 equations and 12 unknowns can be established by substituting formulas (5)-(8) into formula (3). e analytical solutions of the undetermined coefficients A i , B i , C i , and D i were obtained here in MATLAB and then substituted into formulas (5)- (8) to obtain the deflection, rotation angle, bending moment, and shear force generated at each point in the existing pipeline AE during excavation of the BC section. e safety of the existing pipeline can be evaluated based on the results by selecting appropriate judgment standards.

Calculation of Vertical Earth Pressure on the Existing
Pipeline. In this study, the Terzaghi loose Earth pressure calculation model was selected to calculate the vertical Earth Advances in Civil Engineering pressure q 1 on the existing pipeline not affected by excavation. e soil above the excavation section and the excavation-affected section is impacted by the construction disturbance. e mechanical properties of the soil and holding capacity of the downward moving part are reduced, thus increasing the Earth pressure on the pipeline. Based on the simplified calculation and partial safety considerations, the holding capacity of the surrounding soil to the downward moving part can be ignored. e Earth column theory can be directly used to calculate the Earth pressure q 2 and q 3 on the existing pipe jacking above the excavation section and the excavation-affected section.

Determination of the Excavation-Affected Section Length.
e angle between the fracture surface and the tunnel excavation surface can be assumed as (45°−φ/2) in the soilloosening area caused by tunnel excavation [11] (Figure 3). e length of the loose section can be calculated as follows: where a is the length of the excavation-affected section, m; H is the distance between the top of the excavation area and the existing pipeline, m; h is the height of the excavation area, m; and φ is the internal friction angle of soil,°.

Determination of Foundation Reaction Width.
Assuming that the diameter of the existing pipeline is D, the calculated width of soil reaction b is as follows [12]: When b > 1 m, where d is the outer diameter of the existing pipeline, m.

Determination of the Subgrade Coefficient under Construction Disturbance.
e subgrade coefficient is the most important parameter in this theoretical model. ere have been few previous studies on the influence of construction disturbance on subgrade coefficients, but there are many research results regarding other soil mechanical properties under construction disturbance. Xu Yongfu [13] pointed out that the soil disturbance caused by shield tunneling decreases the deformation modulus by 30-70%, and that the cohesion and internal friction angle are related to the soil strain rate: when the soil strain rate increases from 0 to 6%, the cohesion is completely lost, and the internal friction angle decreases linearly to 50% of the original value. Liu Zhe [14] divided tunnel excavation disturbance into three areas according to volumetric strain rate ( Figure 4): area A, where the volumetric strain exceeds 3%; area B, where the volumetric strain is 1-3%; and area C, where volumetric strain is less than 1%. e change range of deformation modulus in each region was found to be greater than that of cohesion and internal friction angle, while the overall change in internal friction angle was minimal. At different points in area A, the deformation modulus decreased by 37-91.7%, the cohesion decreased by 31.4-63.7%, and the internal friction angle decreased by 3.6-12.3%. At different points in area B, the deformation modulus decreased by 8.1-24.1%, the cohesion decreased by 2.5-23.1%, and the internal friction angle decreased by 1.1-2.3%.
For specific projects, the subgrade coefficients of undisturbed soil and disturbed soil should be determined by an in-situ plate load test or indoor consolidation test. If no test data are available, the method proposed by Vesic [15] can be used to determine the subgrade coefficient (formula (8)), where it is not only related to the mechanical properties of the foundation soil itself but also to the action width and stiffness of the beam [16,17]. e subgrade coefficient can be also determined by the known deformation modulus and Poisson's ratio. According to the spatial relationship between the existing pipeline and the underground excavation space, and based on the existing literature, the underground excavation disturbance zoning can be defined as shown in Figure 4. e average deformation modulus and average internal friction angle of the soil in each area can then be estimated. e Poisson's ratio can be obtained by formulas (11)- (12) according to the average internal friction angle, and then, it can be substituted into formula (8) together with the estimated deformation modulus to calculate the subgrade coefficient under the influence of construction disturbance.
e lateral pressure coefficient can be measured with a tri-axial shear meter or lateral pressure meter. When test data are lacking, it can also be determined empirically based on the internal friction angle of each section of soil φ i .

Case Calculation and Analysis
e case analysis data used in this study were taken from the literature [18]. e Wuhan Zhuodaoquan interchange underground pedestrian passage project is located on the east side of the intersection of Luoyu Road, Zhuodaoquan South Road, and Zhuodaoquan North Road. Zhuodaoquan Park is on the south side of the passage, and Zhuodaoquan Middle School is on the north side. It crosses the Zhuodaoquan interchange and Luoyu Road, connecting the sidewalks on the north and south sides. e buried depth of the underground pedestrian passage is 3 m, the width is 8 m, and the height is 3.55 m. e overburdened soil of the underpass is mainly silty clay with cohesion of 15 kPa, internal friction angle of 12°, gravity of 19.3 kN/m 3 , and subgrade coefficient of 1.4 × 10 4 kN/m 3 . ( is subgrade coefficient was calculated by the method in Section 3.3 according to the deformation modulus given in the literature.) ere is a DN100 cast iron pipe above the underpass with an outer diameter of 114 mm and a buried depth of 1.5 m. e existing pipeline settlement calculated by the proposed method was compared against the field-measured settlement under different subgrade coefficients (Table 1) in the excavation section and excavation-affected section (Figures 5-6).
As shown in the figures above, when the subgrade coefficient of the excavation section is 10% of the original subgrade coefficient and the subgrade coefficient of the excavation-affected section is 30% of the original subgrade coefficient, the calculated settlement curve is consistent with the measured settlement curve (calculated value 1 and measured value). e maximum settlement occurs at the midpoint of the concealed excavation section, and the calculated value is slightly greater than the measured value. e widths of both the left and right sides of the settlement curve calculated and measured are about 12 m. ere is only a slight difference between the measured and calculated values of each measuring point. e calculated value slightly exceeds the measurement for the area with a large settlement (0-3 m from the midpoint of the excavation section) but still satisfies the engineering requirements for a safe prediction. With an increase in the distance from the center, the measured curve shows a continuous decreasing trend while the calculated curve shows a sudden change. e reason for this difference is the inherent deficiency of the Winkle elastic foundation beam theory, which does not reflect the continuity of foundation deformation. e results of calculations show that the subgrade coefficient of the excavation section is the most critical parameter in the calculation, which directly affects the settlement calculation results of the maximum settlement and the area with maximum settlement (the interval 0-3m from the midpoint of the excavation section). In this case, when the subgrade coefficient of the excavation section is 25% less than the applicable value (0.1k 1 ), the calculated maximum settlement value increases by 35%. When the subgrade coefficient of the excavation section is 25% greater than the applicable value (0.1k 1 ), the calculated maximum settlement value increases by 20%. It can be seen that the value of this parameter has a great impact on the calculation results and is directly related to the accuracy of prediction and evaluation. erefore, if possible, indoor and in-site tests shall be used to determine the subgrade coefficient of the excavation section. e subgrade coefficient of the excavation-affected section appears to have negligible impact on the maximum settlement, little impact on the settlement calculation results of the maximum settlement area (0-3 m from the midpoint of the concealed excavation section), and a relatively large impact on the second-largest settlement area (3-9 m from the midpoint of the concealed excavation section). If the subgrade coefficient value in the excavation section is accurate and its value in the excavation-affected section (relative to 0.3k 1 ) increases or decreases by 30%, the difference in the calculated maximum settlement and maximum settlement area values is relatively small, and all values are conservative. e difference among calculated values in

Conclusion
Based on the Winkler elastic foundation beam theory and considering the difference of the subgrade coefficient in each area under the existing pipeline, a theoretical calculation model for the effects of underground excavation on an upper-level upper existing pipeline was established in this study. e deflection, rotation angle, bending moment, and shear force of each point of the existing pipeline can be calculated using this model, and then, the safety of the existing pipeline can be evaluated by selecting appropriate judgment criteria. Further, model parameters such as the vertical Earth pressure of each section on the existing pipeline, the length of the excavation-affected section, the calculated width of the foundation reaction, and the subgrade coefficient under the influence of construction disturbance were empirically determined. When there is no test data for reference, the estimated deformation modulus and Poisson's ratio values can be determined based on the disturbance zoning. e subgrade coefficient under construction disturbance can be calculated according to the method proposed by Vesic. e proposed method was used to calculate the settlement of an existing pipeline under the influence of excavation, and the results were compared against measured settlement values. When the subgrade coefficient is appropriate, the calculated settlement curve is in close agreement with the measured settlement curve. With an increase in the distance from the center, the measured curve shows a continuous decreasing trend while the calculated curve shows a sudden change. e reason for this difference is the inherent deficiency of the Winkle elastic foundation beam theory, which does not reflect the continuity of foundation deformation.
e foundation bed coefficient of the excavation section is the most important parameter in the calculation process. Its value significantly affects the final calculation results and is directly related to the accuracy of predictions and evaluations. If possible, both indoor and on-site tests should be conducted to determine the subgrade coefficient of the excavation section. e appropriate value range of the subgrade coefficient of the excavation-affected section is relatively large, and the value deviation is within 30%, which indicates a relatively slight impact on the pipeline settlement calculation results.

Data Availability
All data included in this study are available upon request to the corresponding author.

Conflicts of Interest
e authors declare that they have no conflicts of interest.