Abstract
There are many advantages of reasonably determining the affected zone and its construction sequence for the multiline parallel pipe-jacking construction, such as convenient to reasonably arrange the pipeline section distribution mode, reducing the construction safety risk, and construction difficulty. Combined with the engineering construction practice of the multiline parallel pipe-jacking project for the north city drainage and flood control project through Chu River in Hefei, the strength reduction method has been utilized to study the surrounding stratum FOS for the pipe-jacking construction unloading under different overburden thickness Hs, different clear distance D, and different river water depth in this paper. The formula of the minimum critical overburden thickness Hsmin and the minimum critical clear distance Dmin of the surrounding stratum self-stability during pipe-jacking construction unloading have been derived. The division method of the affected zone for the multiline parallel river-crossing pipe-jacking construction which behaves as the mutual influence nonself-stability zone, the mutual influence self-stability zone, and the no mutual influence self-stability zone has been proposed. For further discussion, the construction sequence of the multiline parallel river-crossing pipe jacking has been studied.
1. Introduction
Because of many advantages such as the little impact on the environment, the small construction area, the fast construction speed, and the high degree of mechanization, slurry balance pipe jacking has become an important construction method of underground pipeline construction through rivers, roads, buildings, etc. [1–4]. Due to the factors of the topographical conditions, functional requirements, and others, there are more and more multiline parallel pipe-jacking projects with small clear distance D and thin overburden thickness Hs, such as four parallel pipe jacking under the Guan River [5, 6] and three-hole parallel adjacent jacking in the drainage project of Meilan International Airport [7]. The construction process of multiline parallel pipe jacking is the mechanical process in which an incomplete structure with a gradual change in geometric shape and material properties is subjected to the change in construction loads in time and space. Early construction pipe jacking will impact on the surrounding environment of postconstruction pipe jacking by influencing the displacement and stress of its surrounding strata, while postconstruction pipe jacking will impact on early construction pipe jacking that has been completed during the construction process [8–12]. The repeated cross influence of multiple construction for multiline parallel pipe jacking will make the stress repeatedly adjusted and form the stress field closely related to the spatial distribution of pipe jacking with a significant space-time effect and group-hole effect [13, 14]. They have important theoretical significance and engineering application value to explore the self-stability ability of the surrounding strata during the construction of multiline parallel pipe jacking and define the reasonable construction sequence.
The main methods for studying on the surrounding stratum self-stabilizing capacity and the deformation characteristics of the multiline parallel pipe-jacking construction include the theoretical analysis method [15–19], model test method [20, 21], in situ test method [22, 23], and numerical computational method [24–29]. At present, the numerical calculation method has become an increasingly important tool in urban underground engineering research due to its advantages of the high speed, high accuracy, low cost, time-saving and labor-saving, and strong adaptability to complex conditions and processes. The finite element strength reduction method (FESRM) is one of the numerical calculation methods. Because of many advantages of the strict mechanical basis, such as being quantifiable, visible, dynamic, and intuitive, FESRM has been widely used in the self-stability evaluation of surrounding rock stability of tunnels and underground projects in recent years [30–36].
At present, the research on the multiline parallel pipe-jacking construction mainly focuses on the deformation of surrounding strata and pipe-jacking segments, construction parameters, and other aspects. The research on the subject of the affected zone and construction sequence is mainly carried out in combination with specific projects, which is not systematic. Combined with engineering construction practice of the multiline parallel pipe-jacking project for the north city drainage and flood control project through Chu River in Hefei China, the strength reduction method has been utilized to study the surrounding stratum FOS for pipe-jacking construction unloading under different overburden thickness Hs, different clear distance D, and different river water depth . The affected zone division idea and the construction sequence determination method of multiline parallel pipe-jacking construction are summarized and improved in this paper. The research results provide a case support and the theoretical basis for the multiline parallel river-crossing pipe-jacking construction with the thin overburden thickness and the narrow clear distance.
2. Mechanics and Engineering Profile
2.1. Mechanical Mechanism of Multiline Parallel Pipe-Jacking Construction
The construction method of slurry balance pipe jacking is carried out by using the pipe-jacking machine located in the starting well through the main pushing system. The cutter head of the pipe-jacking machine cuts the stratum in front of the pipe-jacking face and then discharges the stratum soil through the mud treatment system and transportation system. At the same time, under the continuous pushing of the main pushing system, the subsequent segments are pushed following the pipe-jacking machine one by one, until reaching the receiving well. In addition, the drag reduction slurry shall be properly injected behind the pipe-jacking segment according to the actual situation, and the whole pipe-jacking construction process is completed. The acting force acted on surrounding strata during pipe-jacking construction mainly includes the jacking force acting on the position of the pipe-jacking face and the frictional force between the pipe-jacking segment and the surrounding strata, as shown in Figure 1. When the pipe-jacking project is under normal construction, the stratum stress state near the pipe-jacking head in front of the machine is extremely complex due to the effect of the pipe-jacking force. The stratum stress state in front of the pipe-jacking machine increases due to the extrusion effect of the front additional thrust of the pipe-jacking face. After the pipe-jacking machine passes, as the outer diameter of the pipe-jacking machine is slightly larger than the outer diameter of the pipe segment, the stratum loss is formed between the subsequent pipe segment and the surrounding strata. The stratum moves towards the gap, and then, stress diffusion occurs. Simultaneous grouting behind the pipe segment causes the surrounding strata to be encapsulated. When pipe-jacking construction is completed, the surrounding stratum consolidation settlement will produce under the action of weight. During the whole process of pipe-jacking construction, the stress states of the pipe-jacking strata around are in a continuous dynamic change, which causes the continuous change of stratum displacement. When the overburden thickness is thin, the deformation of the strata around the pipe extends to the surface and causes surface deformation.

When the clear distance between the early and the later stages of pipe-jacking construction is large, the stratum disturbance during the early stage has no impact on the later stage; meanwhile, the disturbance of the strata during the later stage has no impact on the early stage. When the clear distance between the early and the later stages of pipe-jacking construction is small, the later stage of pipe-jacking construction is completed in the strata affected by the disturbance of the early stage. At the same time, the stratum disturbance generated by the later stage has an impact on the early stage. The mutual influence degree becomes stronger as the clear distance decreases. When the clear distance between multiple parallel pipe-jacking is small, the change in the stratum stress state has been caused by the second pipe-jacking construction generating additional stress on the surrounding stratum of the first pipe jacking, which will change the stress and displacement. This effect is greater on the side close to each other and less on the side far away from each other, as shown in Figure 2(a). The third pipe-jacking construction has the same phenomenon, which generates the new equilibrium state and reacts on the first pipe jacking, the second pipe jacking, and the third pipe jacking through the stress state. Among them, the impact between two jacking pipes is greater when they are close to each other and less when they are far away from one side, as shown in Figure 2(b). In this case, three construction sequence problems arise (see Figure 3). When there are 4 pipe-jacking operations, the relative position relationship between pipes increases to 6 forms. When there are npipe-jacking operations, the relative position relationships between pipes increase to C2n forms.

(a)

(b)

2.2. Project Overview
The north city drainage and flood control project is located in Hefei, China. The total length of the project line is 7015 m. Y25-Y26 of this project pipeline network crosses the main channel of Chu River. In this section, the width of the river channel is 55 m, the width is 25 m, the embankment height is 9.46 m, the slope is 45°, and the river depth is in the range of 2.0∼5.0 m, as shown in Figure 4(a). Slurry balance mechanical pipe-jacking construction is adopted. Three circular reinforced concrete pipes with an inner diameter of 3.0 m and an outer diameter of 3.6 m are used for pipe jacking. The clear distance between two adjacent pipe jacking is 2.8 m. The stratum where pipe jacking crosses is silty clay. The thickness of the pipe-jacking covering layer in the riverbed section is 3.0 m, as shown in Figure 4(b). The vertical horizontal distance from the riverside of the pipe-jacking working shaft to the edge of the valley is about 10 m.

(a)

(b)
The project site belongs to the undulating plain landform of Jianghuai, and the microgeomorphology is the hillock with a concave groove. The regional geological structure belongs to the southern edge of the North China Plateau, and the secondary tectonic unit belongs to the Hefei Basin. The stratum distribution from top to bottom is as follows: fill soil, 0.5–1.5 m; silty chalky clay, 1.2–2.7 m; clay, 5.9–7.5 m; chalky clay, 4.8–6.0 m; and strongly weathered mudstone, 1.7–2.6 m, and the following are medium and slightly weathered mudstones. The physical and mechanical parameters of each stratum are shown in Table 1.
3. Study on the Division of the Affected Zone under Construction Unloading of Multiline Parallel River-Crossing Pipe Jacking
3.1. Research Methodology
3.1.1. Finite Element Strength Reduction Method
Due to the high aspect ratio, the problem of multiline parallel pipe-jacking construction unloading can be considered as a plain-strain one. The failure complies with the Mohr–Coulomb criterion which can be expressed as shown in the following formula:where I1 is the first invariant of the stress tensor, J2 is the second invariant of the partial stress tensor, c is cohesion, φ is the angle of internal friction, and θσ is Lodder’s angle.
The finite element strength reduction method (FESRM) is a method that combines the strength reduction technique, the principle of ultimate equilibrium, and the principle of elastic-plastic finite element calculations. First, the state of force and deformation under the original parametric working condition of the stratum was calculated. Then, the stratum strength parameters c and φ were simultaneously discounted according to equation (2) to obtain a new set of strength parameters c′ and φ′, and they were used as the new material strength parameters for calculation. Finally, the calculation was carried out by continuously adjusting the discount factor k until the surrounding strata were in ultimate equilibrium, and the critical rupture surface was obtained; at the same time, the discount factor k of the material was FOS. The calculation principle of FESRM [23] is shown in Figure 5:where k is the discount factor.

3.1.2. Computational Model
Finite element numerical calculation software of Midas GTS NX was utilized to study the surrounding stratum self-stability characteristics of the pipe-jacking construction. The numerical calculation was based on the plane strain problem. Pipe jacking took a circle with an outer diameter of 3.6 m. The distance between the left and the right boundaries of the calculation model was more than 3 times of the pipe-jacking outer diameter, and the horizontal displacement constraint was applied. The distance between the lower boundary and the bottom of the calculation model was more than 3 times of the pipe-jacking outer diameter, and the vertical displacement constraint was applied to the lower boundary, with the upper surface free. The model boundary calculation cell grid was set to 1.0 m longitudinally and horizontally, and the pipe perimeter grid was 0.5 m, as shown in Figure 6. The calculation used the DP4 equivalence Mohr–Coulomb yield criterion. Jacking construction unloading took place in the one-time full-section excavation. The initial stress took into account the self-weight of the ground soil and river water pressure, and no other construction process factors were considered. Numerical calculations of the physical and mechanical parameters of the strata, the dimensions of pipe jacking, and the relative position relationships were based on dependent engineering parameters as basic data and were extended to the general case on this basis. The method of controlling single variable analysis was adopted while discussing the self-stability characteristics of the unloading stratum during pipe-jacking construction; that is, when one calculation parameter changed, other calculation parameters remained unchanged.

3.2. Study on the Surrounding Stratum Self-Stability Characteristics of the Single-Line Pipe-Jacking Construction
In this study, was taken as 2.0 m, 5.0 m, and 8.0 m, respectively, and the stratum of pipe-jacking crossing was silty clay stratum. The FOS calculation results of the pipe-jacking construction surrounding strata under different Hs are shown in Table 2.
From Table 2, it was shown that the overall change trend of FOS increased first and then decreased with an increase in Hs under the same conditions of other factors (see Figure 7). It illustrated the problem of the surrounding stratum self-stabilizing ability for pipe-jacking construction increased first and then decreased with an increase in Hs. The fundamental reason for the above phenomenon was that the self-stability characteristics of the pipe-jacking construction surrounding strata were mainly determined by the relative relationship between the stratum strength itself and the surrounding stratum redistribution stress caused by pipeline construction unloading. The surrounding stratum redistribution stress by pipeline construction unloading was positively related to the stratum initial stress itself under the same conditions of other factors. The stratum initial stress was composed of the effect for stratum self-weight and for the river water self-weight. With an increase in stratum depth, the effect of the former on the stratum initial stress gradually increased, while the latter gradually decreased. Thus, the stratum initial stress decreased first and then increased with an increase in stratum depth. Furthermore, the surrounding stratum redistribution stresses the same rule. The stratum strength itself was basically unchanged. Therefore, the ratio of stratum strength to the surrounding stratum redistribution stress increased first and then decreased.

FOS could quantitatively evaluate the degree of surrounding stratum self-stability by pipe-jacking construction unloading. Taking FOS = 1.25 as limitation, the calculation results of the surrounding stratum minimum critical overburden thickness d, the maximum critical overburden thickness r, and the single-line pipe-jacking construction under different river water depth are shown in Table 3. The mathematical fitting equations between and are shown in equations (3) and (4), respectively:
The basic goal of an excellent underground space development plan was to make full use of the stratum self-stability ability, minimize the engineering auxiliary measures, and pay attention to the economy while ensuring safety. The realization of this goal mainly depends on the geological conditions of the site, the understanding of these geological conditions, and the ability to use them as design. Taking as the upper and the lower limitation, respectively, the surrounding stratum self-stability distribution zone of single-line pipe-jacking construction unloading is shown in Figure 8. For the specific environmental conditions of pipe-jacking construction, there are many advantages for the pipe-jacking vertical section that was designed in the surrounding stratum self-stability distribution zone, such as reducing the construction safety risk and construction difficulty.

3.3. Study on the Surrounding Stratum Self-Stability Characteristics for the Two-Line Parallel Pipe-Jacking Construction
In this study, was taken as 2.0 m, 5.0 m, and 8.0 m, respectively, Hs was taken as 3.0 m, and the stratum of pipe-jacking crossing was the silty clay stratum. The FOS calculation results of the pipe-jacking construction surrounding strata under different D are shown in Table 4.
From Table 4, it was shown that the overall change trend of FOS increased first and then decreased with an increase in D under the same conditions of other factors (see Figure 9). With an increase in D, the potential failure pattern of the surrounding strata for double-line parallel pipe-jacking construction unloading transited from the double pipe-jacking overall collapse to the gradual separation collapse, and the final collapse pattern was consistent with that of single pipe-jacking (see Figure 10). It was demonstrated that when D was small, the mutual influence was strong of double-line parallel pipe-jacking construction unloading. With D increased, the mutual influence was gradually weakened, and when D increased to a certain extent, it would be no longer affecting each other.


(a)

(b)

(c)
Taking FOS = 1.25 as the limitation, the calculation results of the surrounding stratum minimum critical clear distance Dmin for the double-line pipe-jacking construction under different buried depths Hs are shown in Table 5. The mathematical fitting equation among the minimum critical clear distance Dmin, the overburden thickness Hs, and the river water depth of double-line parallel pipe-jacking construction unloading is shown in equation (5). The spatial feature distribution map among them is drawn in Figure 11, which provided the theoretical basis for the cross-sectional design of the multiline parallel pipe-jacking construction:

3.4. Study on the Surrounding Stratum Self-Stability Characteristics for the Three-Line Parallel Pipe-Jacking Construction
In this study, was taken as 2.0 m, 5.0 m, and 8.0 m, respectively, Hs was taken as 3.0 m, and the stratum of pipe-jacking crossing was the silty clay stratum. Three working conditions were taken for the construction sequence, which were 1# ⟶ 2# ⟶ 3#, 2# ⟶ 1# ⟶ 3#, and 1# ⟶ 3# ⟶ 2# (see Figure 9). The FOS calculation results of the pipe-jacking construction surrounding strata under different D are shown in Table 6.
From Table 4, it was shown that with the pipe-jacking number increased, FOS tends to decrease gradually on the whole, and the smaller the D was, the more significant the trend was. When the pipe-jacking construction was completed, FOS was not affected by the construction sequence. However, when the sequence was different, FOS in the construction process was different (see Figure 12). The overall change trend of FOS increased first and then decreased with an increase in D under the same conditions of other factors (see Figure 13). With an increase in D, the potential failure pattern of the surrounding strata for three-line parallel pipe-jacking construction unloading transited from the double pipe-jacking overall collapse to the gradual separation collapse, and the final collapse pattern was consistent with that of single pipe jacking (see Figure 14). It was demonstrated that when D was small, the mutual influence was strong of three-line parallel pipe-jacking construction unloading; with D increased, the mutual influence was gradually weakened, and when D increased to a certain extent, it would be no longer affecting each other.



(a)

(b)

(c)
Taking FOS = 1.25 as the limitation, the calculation results of the surrounding stratum minimum critical clear distance for the three-line pipe-jacking construction under different buried depths Hs are shown in Table 7. The mathematical fitting equation among the minimum critical clear distance , the overburden thickness Hs, and the river water depth of double-line parallel pipe-jacking construction unloading is shown in equation (6). The spatial feature distribution map among them is drawn in Figure 15:

3.5. Dividing Mutual Influence Zones of Multiline Parallel Pipe-Jacking Construction Unloading
The surrounding stratum self-stability degree of multiline parallel pipe-jacking construction unloading could be determined quantitatively by FOS, and the mutual influence degree of them could be determined quantitatively by the change in FOS. For the convenience of expression, when the multiline parallel pipe-jacking construction has been completed, FOS was uniformly called FOS-M. Meanwhile, the single-line pipe-jacking construction has been completed, and FOS was uniformly called FOS-S. Taking FOS = 1.25 as the limitation, based on the surrounding stratum self-stability degree and the mutual influence degree, the surrounding stratum zone of the multiline parallel pipe-jacking construction unloading could be divided.
When FOS-M ≤ FOS-S < 1.25, it meant that the mutual influence would occur during multiline parallel pipe-jacking construction unloading, and the surrounding stratum zone did not meet the self-stability requirement. Thus, the surrounding stratum zone was called the mutual influence nonself-stability zone.
When 1.25 < FOS-M < FOS-S, it meant that the mutual influence would occur during multiline parallel pipe-jacking construction unloading, and the surrounding stratum zone met the self-stability requirement. Thus, the surrounding stratum zone was called the mutual influence self-stability zone.
When FOS-M = FOS-S > 1.25, it meant that the mutual influence would not occur during multiline parallel pipe-jacking construction unloading, while the surrounding stratum zone met the self-stability requirement. Thus, the surrounding stratum zone was called the nonmutual influence self-stability zone.
4. Study on the Construction Sequence of the Multiline Parallel River-Crossing Pipe Jacking
4.1. Basic Assumptions
The stratum initial stress had been changed by pipe-jacking construction, which led to a series of complex physical and mechanical effects. It was not only related to the physical properties of the stratum but also closely related to the construction method, construction process, and relevant construction parameters. When using numerical software for calculation, it was difficult to take all factors into account to completely reflect the construction process. Therefore, proper simplification should be carried out during numerical calculation. It could not only meet the requirements of calculation software but also make the numerical calculation result reflect the construction process well.
According to the actual situation of the supporting project, the numerical calculation was based on the following assumptions in this paper:(1)All strata were homogeneous, continuous, and isotropic ideal elastic-plastic material. The influence of groundwater is ignored in the calculation process. The initial stress only considers the stratum self-weight stress. Stratum settlement was only considered due to pipe-jacking construction, and the consolidation settlement was ignored.(2)Jacking pressure of the pipe-jacking excavation surface was applied to the whole excavation surface in the form of a circular uniformly distributed load. The pressure was taken as the lateral static earth pressure at the center of the pipe-jacking excavation surface. According to the project actual situation, the pressure was 88.5 kpa in this paper, as shown in Figure 16(a).(3)The influence of grouting pressure on the construction process was not considered. The grouting unit around pipe jacking was distributed along the radial and equal thickness of the segment. The elastic modulus was taken as 1/50 of the original formation unit, and the thickness of the grouting equivalent layer was 2 cm, as shown in Figure 16(b).(4)The pipe joint material was an isotropic linear elastic body, and the indirect head effect of a pipe joint was ignored. The friction resistance between pipeline and stratum acted on the outer surface of the pipe casing and the inner surface of the soil around the pipe, with the same size and opposite direction. The frictional resistance was a certain value and evenly distributed along the pipe-jacking direction. The frictional resistance was achieved by setting a friction coefficient on the contact surface, and the value was 3.5 kpa in this paper, as shown in Figure 16(c).

(a)

(b)

(c)
4.2. Calculation Model and Implementation Process
The length, the width, and the height of the numerical calculation model for the multiline parallel pipe-jacking construction were 80 m, 80 m, and 30 m, respectively. The left and right boundaries of the calculation model were subject to horizontal displacement constraints, the lower boundary was subject to vertical displacement constraints, and the upper surface was free. The pipe-jacking segment was regarded as an elastic material, and the thickness was 0.30 m. The shell of the pipe-jacking machine adopted the linear elastic model. The material property adopted a 2D plate element, and the thickness was 0.06 m, as shown in Table 8. The material properties of strata and segments adopted 3D unit entities, which were divided into meshes by geometry and then expanded by 2D meshes. The calculation unit grid around pipe jacking was set to 0.5 m in both vertical and horizontal directions, while the model boundary was 2.0 m. The values of the stratum physical and mechanical parameters are shown in Table 1, and the calculation model is shown in Figure 17.

The numerical calculation process of the multiline parallel pipe-jacking construction was divided into the following stages: pipe-jacking excavation, pipe-jacking advance, and segment application. The specific steps were as follows: (1) The initial boundary conditions of the model and the stratum self-weight were applied, the initial displacement was cleared, and the analysis of the initial stress field was started. (2) The jacking pressure of the pipe-jacking face was applied, the function of the stratum excavated by pipe-jacking construction was deactivated, and the friction resistance was activated; at the same time, the casing of the pipe jacking machine was activated. (3) The shell of the pipe-jacking machine was deactivated, the pipe-jacking machine was jacked, and the pipe-jacking segment was applied. (4) The stratum of the next pipe segment was excavated, and the above steps were repeated.
4.3. Calculation and Analysis
There were differences in the degree of disturbance to the surrounding strata caused by the different jacking sequences of the multiline parallel pipe-jacking construction, which made the stratum displacement different. It was very important to control the stratum displacement by relying on the pipe-jacking project passing through the river channel. The cumulative vertical displacement of the stratum under various working conditions is shown in Figure 17. Displacement monitoring points were set at the river bottom surface section at the river center, which was perpendicular to the pipe-jacking axis. The surface displacement monitoring points directly above 1#, 2#, and 3# pipe jacking were marked as S1, S2, and S3 in turn (see Figure 18).

(a)

(b)

(c)
The calculation results of the stratum deformation at S1, S2, and S3 points in the construction steps of the construction sequences for 1# ⟶ 2# ⟶ 3#, 2# ⟶ 1# ⟶ 3#, and 1# ⟶ 3# ⟶ 2# were extracted, respectively, as shown in Figure 19.

(a)

(b)

(c)
From Figure 19, it was shown that the overall law of stratum displacement at S1, S2, and S3 points in the construction steps under different construction sequences was basically consistent. On the whole, the stratum uplifted in the front of the pipe-jacking cutter head and subsided in the rear. The former was caused by the jacking force of the pipe-jacking face, and the latter was caused by the formation loss caused by the cutter head of the pipe-jacking machine slightly larger than the outer diameter of the segment. The mutual influence of stratum displacement was very strong caused by the three-line parallel pipe-jacking construction.
The calculation results of the pipe jacking in the construction sequence of “1# ⟶ 2# ⟶ 3#” were taken as examples for detailed explanation. During the construction of 1# pipe jacking, with the distance between the cutter head and the river center continuously approaching, the stratum displacement at S1, S2, and S3 was continuously uplifted and the value gradually increased. When the cutter head was pushed directly below the river center, the uplift of the stratum displacement at S1, S2, and S3 reached the maximum. When the cutter head passed through the river center, the stratum displacement at S1 rapidly sank to the negative value (settlement), and the displacement at S2 slightly sank. With the distance between the cutter head and the river center getting further away, the stratum displacement at S1, S2, and S3 was becoming slower and tending to be stable. During the construction of 2# pipe jacking, with the distance between the cutter head and the river center continuously approaching, the stratum displacement at S2 and S3 was continuously uplifted and the value gradually increased, while S1 showed the rising trend, but the value was small. When the cutter head was pushed directly below the river center, the uplift of the stratum displacement at S2 and S3 reached the maximum, while S1 rose to the maximum value. When the cutter head passed through the river center, the stratum displacement at S2 rapidly sank to a negative value (settlement), while S1 and S3 slightly sank. With the distance between the cutter head and the river center getting further away, the stratum displacement at S1, S2, and S3 was becoming slower and tending to be stable. 3# pipe-jacking construction had little impact on the stratum displacement at S1. During the construction of 3# pipe-jacking, with the distance between the cutter head and the river center continuously approaching, the stratum displacement at S3 was continuously uplifted and the value gradually increased, while S2 showed the rising trend, but the value was small. When the cutter head was pushed directly below the river center, the uplift of the stratum displacement at S3 reached the maximum, while S2 rose to the maximum value. When the cutter head passed through the river center, the stratum displacement at S3 rapidly sank to a negative value (settlement), while S2 slightly sank. With the distance between the cutter head and the river center getting further away, the stratum displacement at S2 and S3 was becoming slower and tending to be stable.
When the first, the second, and the third pipe-jacking has been completed, respectively, in the sequence of 1# ⟶ 2# ⟶ 3#, 2# ⟶ 1# ⟶ 3#, and 1# ⟶ 3# ⟶ 2#, the calculation results of river bottom stratum deformation at the river bottom surface section in the river center, which was perpendicular to the pipe jacking axis, are shown in Figure 20.

(a)

(b)

(c)
From Figure 20, it is shown that the overall law of stratum displacement at the river center section under different construction sequences was basically consistent. The mutual influence of stratum displacement was very strong caused by the three-line parallel pipe-jacking construction. When the single-line pipe-jacking construction was completed, the horizontal surface displacement curve was approximately in a normal distribution. The maximum settlement was located at the pipe-jacking axis, and the maximum settlement value was 18.73 mm and gradually decreased from the axis to both sides. When the double-line pipe-jacking construction was completed and the clear distance D was narrow (in the sequence of 1# ⟶ 2# or 2# ⟶ 1#), the horizontal surface displacement curve showed the partial v-shaped distribution, and the surface settlement maximum value was 25.10 mm, which was located between two pipe jacking and in the pipe-jacking side that was constructed first. When the clear distance D was wide (in the sequence of 1# ⟶ 3# or 3# ⟶ 1#), the horizontal surface displacement curve showed the partial w-shaped distribution. When the three-line parallel pipe-jacking construction was completed in the sequence of 1# ⟶ 2# ⟶ 3#, 2# ⟶ 1# ⟶ 3#, and 1# ⟶ 3# ⟶ 2#, the maximum values of surface settlement were 26.46 mm, 26.64 mm, and 23.82 mm, respectively.
4.4. Construction Sequence
There was no construction sequence problem of the multiline parallel pipe-jacking in the mutual influence non-self-stability zone and the nonmutual influence self-stability zone. The impact of the multiline parallel pipe-jacking construction in the mutual influence self-stability zone on the surrounding stratum was not that the simple superposition of the individual single-line pipe-jacking construction. Early construction pipe jacking would impact on the surrounding strata of postconstruction pipe jacking by influencing the displacement and stress, while postconstruction pipe jacking would impact on early construction pipe jacking that has been completed during the construction process. Different construction sequence had different impacts on the surrounding strata. Pipe-jacking construction should avoid or reduce the mutual influence of the adjacent pipe-jacking construction as much as possible. Therefore, the reasonable selection of the pipe-jacking construction sequence was related to the safety, efficiency, and even the success or failure of the project construction to a certain extent.
The construction sequence scheme of the multiline parallel pipe-jacking construction for the north city drainage and flood control project through Chu River was established according to the following steps. First, we ensure that technology was feasible, safe, and reliable. Then, the characteristics of stratum deformation caused by the multiline parallel pipe-jacking construction were comparatively analyzed, and the working condition with the minimum stratum deformation was taken. Then, comparative analyses were carried out for construction organization, construction period, project economy, and other factors. Finally, the jacking sequence of “1# ⟶ 3# ⟶ 2#” as the optimal sequence was confirmed. The project was officially launched on June 10 in 2022, and the construction was successfully completed on June 30 in 2022. The project implementation effect was good.
5. Conclusions
(1)It was revealed that the variation law of the surrounding stratum self-stability characteristics for single-line pipe-jacking construction unloading varies with the overburden thickness Hs. The surrounding stratum minimum critical overburden thickness Hsmin and the maximum critical overburden thickness Hsmax of single-line pipe-jacking construction unloading were obtained. The surrounding stratum self-stability distribution zone of single-line pipe-jacking construction unloading was drawn, which provided the theoretical basis for the vertical-sectional design of river-crossing pipe jacking with overburden thickness Hs drastic change.(2)It was revealed that the variation law of the surrounding stratum self-stability characteristics for multiline parallel pipe-jacking construction unloading varies with the clear distance D. The mathematical fitting equations among the minimum critical clear distance Dmin, the overburden thickness Hs, and the river water depth of multiline parallel pipe-jacking construction unloading were obtained. The spatial feature distribution map among the minimum critical clear distance Dmin, the overburden thickness Hs, and the river water depth was drawn, which provided the theoretical basis for the cross-sectional design of the multiline parallel pipe-jacking construction.(3)Based on the surrounding stratum self-stability degree and the mutual influence degree of multiline parallel pipe-jacking construction unloading, the surrounding stratum zone was divided into the mutual influence nonself-stability zone, the mutual influence self-stability zone, and the nonmutual influence self-stability zone.(4)Combined with the relying project, the construction sequence of “1# ⟶ 3# ⟶ 2#” was concluded to be the optimal construction sequence through the comparative analysis of the ground deformation characteristics caused by the multiline parallel pipe-jacking construction and the combination of various factors such as engineering safety, economy, and convenience.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.
Authors’ Contributions
Ziguang Zhang was the person in charge of the paper and completed the construction of thesis framework system and thesis writing. Jiesheng Zhang completed the collection, sorting, and analysis of relevant project case data in the paper. Ruijin Mao, Xuefeng Wang, and Mengqing Zhang assisted in completing the related work of numerical calculation in the paper.
Acknowledgments
This study was funded by the Provincial Natural Science Research Project of Colleges and Universities in Anhui Province-Key projects (KJ2021A0611), the Scientific Research Project of Anhui Jianzhu University (no. XJ2019000502), the Science and Technology Plan of Housing and Urban-Rural Construction in Anhui Province (2022-YF096 and 2020-YF38), and the Science and Technology Development Project (HYB20190152, HYB20220092, and HYB20220162).