Analysis of Ground Movement during Large-Scale Pipe Roof Installation and Artificial Ground Freezing of Gongbei Tunnel

This paper analyzes the vertical ground movement during large-scale pipe roof installation and artiﬁcial ground freezing of Gongbei tunnel of the Hong Kong-Zhuhai-Macau bridge project. The transverse ground settlement during pipe roof installation is analyzed. The ground loss volume ratio and settlement trough width coeﬃcient during pipe jacking are estimated based on the ﬁeld measurement of ground settlement. The interaction of pipes during multiple jacking is investigated. The eﬀect of frost heave control by pregrouting, limiting frozen soil thickness, and combination of the two methods is evaluated. The analysis shows that the ground settlement during pipe roof installation by jacking 37 pieces of 1620mm steel pipes is relatively small with a maximum value of 2.2 cm. The reinforcement to ground provided by the fore-jacked pipes reduces the ground loss volume ratio and, consequently, the ground settlement during the follow-up pipe jacking. The artiﬁcial ground freezing generates a relatively large ground heave with a maximum value of 7.8cm. Pregrouting plays a critical role in the frost heave control by reducing the heave by about 33%. Limiting the frozen soil thickness by heating pipes serves as an eﬀective supplement to frost heave control by reducing the heave by about 9%. The combination of the two measures reduces the ground heave by about 42%. Findings from this paper provide valuable reference to the tunnel construction using pipe roof and artiﬁcial ground freezing as presupport.


Introduction
With the advancement of tunnelling technique and increasing traffic demand, more and more large cross-section tunnels have been built in recent years. e International Tunneling Association (ITA) defines the large cross-section tunnels with an excavation area larger than 100 m 2 . Many tunnels with an excavation area larger than 200 m 2 have been built, for example, the 336.8 m 2 Gongbei tunnel in China [1], the 225 m 2 Katabira tunnel in Japan [2], and the 252 m 2 channel tunnel in Europe [3]. Ground movement control creates an engineering challenge during the excavation of large cross-section tunnels, especially for the shallowly buried ones in poor ground, such as soft clay/silt or loosely saturated sand. Previous construction experience indicates that the conventional sequential excavation supplemented with ground improvement, such as pregrouting, is difficult to excavate tunnels with the crosssectional area larger than 200 m 2 in poor ground. Strong presupport is required to excavate such tunnels. e recently built Gongbei tunnel of the Hong Kong-Zhuhai-Macau bridge project has an excavation cross-sectional area of 336.8 m 2 and burial depth of 4 to 5 m in interlayers of saturated clay/silt/sand. It applied a large-scale pipe roof in combination of artificial ground freezing as presupport. e pipe roof, composed of 36 pieces of 1620 mm diameter, 255 m length steel pipes, is designed to support the Earth load, while the artificial ground freezing is designed to seal the pipe roof to prevent water infiltration. Under such a strong presupport, the excavation induced ground movement is the minimum. However, large ground movement may occur during the installation of the pipe roof and artificial ground freezing, and the ground movement needs to be assessed. e pipe roof is installed through microtunnelling which generates ground movement. e ground settlement can be estimated by empirical, analytical, and numerical methods verified by field measurements or laboratory tests. e empirical methods provide a fast estimation of the ground movement and are widely used in engineering practice. e empirical equations assume that the transverse ground settlement follows a quasi-normal distribution, defined by the maximum settlement and the width coefficient of settlement trough [4]. e maximum settlement is characterized by the ground loss volume V s , typically estimated using the ground loss volume ratio (η). e η value is impacted by mainly factors such as the machine type, ground condition, quality of compensation grouting, and local construction experience [5][6][7][8][9][10][11] and is reported in many publications. For example, O'Reilly and New early reported that the η value ranges from 0.5-2.5% for cohesive soils and 2.0-10.0% for cohesionless soils in England [9]. Zhu and Li recently reported that the η value of cohesive and cohesionless soils is about 1.10% and 1.18%, respectively, based on summarization of the data from 20 regions in China [11]. e settlement trough width coefficient i is related to ground condition, tunnel diameter, tunnel burial depth, etc. [12][13][14]. Different equations were proposed for the characterization of the width coefficient i [9,15,16]. For example, Burland et al. reported that the normalized width coefficient k � i/z, where z is the tunnel burial depth from ground surface to tunnel centre, ranges from 0.2 to 0.3 for cohesionless soil, 0.4 to 0.5 for hard clay, and as high as 0.7 for silty clay based on the local experience in London [16]. Although the ground settlement induced by tunnelling is widely studied, the research is mostly focused on single-or double-line parallel tunnelling. Very limited study is found on the ground settlement during multiple parallel tunnelling as used for the pipe roof installation.
Frost heave control is a critical issue during artificial ground freezing. e frost heave mainly consists of two components, the in situ phase change one and the segregation one [17]. Between them, the segregation frost heave is the major one. e in situ phase change frost heave is caused by the icing of the pore water in the original soil, and the segregation frost heave is caused by the icing of pore water migrated from surrounding soils. e soil freezing is a complex hydrothermal-coupling process. e amount of segregation frost heave mainly depends on the change of temperature field and the rate of water migration. Controlling the temperature field and inhibiting the water migration are two major aspects of frost heave control [18]. Lai et al. analyzed the influence of temperature gradients on the frost heave of saturated soils. After freezing 90 h, the frost heave caused by the temperature gradient of 0.65°C/cm is 1.97 times of that caused by temperature gradient of 0.31°C/ cm [19]. Zhang et al. analyzed the frost heave of soils with high permeability in cold regions and showed that controlling the temperature field only is difficult to control the frost heave in high permeable soils. To control the segregation frost heave, it has to reduce the permeability of the soil to inhibit the water migration. Pregrouting can reduce the soil permeability effectively and is widely used to control frost heave [20]. Hu reported that when the permeability of the soil decreases from 1.7 × 10 −6 to 6 × 10 −7 cm/s, the frost heave ratio reduces from 4.5% to 2% [21]. e previous studies mainly focused on the frost method, such as pregrouting or controlling temperature field, is usually applied to control the frost heave. e combination of multiple control methods is rarely applied in engineering practice. e successful construction of the Gongbei tunnel provides valuable field measurement data.
is paper investigates the vertical ground deformation during the installation of large-scale pipe roof and artificial ground freezing using Gongbei tunnel as a case project. Following the Introduction, Section 2 of the paper provides a brief description of the Gongbei tunnel. Section 3 analyzes the ground settlement and pipe interaction during multiple parallel jacking. Section 4 investigates the effect of three frost heave control methods, including pregrouting, limiting frozen soil thickness using heating pipes, and combination of the two methods. Section 5 discusses the total vertical ground deformation during pipe roof installation and artificial ground freezing. At the end, the conclusions are made in Section 6.

General Description.
e Gongbei tunnel is one of the most technically challenging components of the Hong Kong-Zhuhai-Macau bridge project. It is located beneath the Gongbei customer port operated by the Macau special administrative region (SAR) and Zhuhai city in southeast China. It has a length of 255 m along transitional and circular curves with a curvature ranging from 886 to 906 m, as shown in Figure 1. It is a two-way, three-lane in each direction roadway tunnel with a cross-sectional area of 336.8 m 2 and burial depth of 4 to 5 m. Figure 2 shows the tunnel cross section and subsurface conditions. ere are a variety of buildings, utilities, and piles adjacent to the excavation. e closest distance between the adjacent piles and excavation boundary is about 1.6 m. e tunnel is mainly buried in saturated soft silt and clay with interlayers of fine sands. e subsurface strata consist of miscellaneous fill (stratum ①), soft organic clay and silt (stratum ③-1), fine to coarse sand (stratum ③-3), sandy clay and silt (stratum ④-3), fine sand (stratum ⑤-2), fine to coarse sand (stratum ⑥-2), and silty clay (⑦-1). e groundwater table ranges between El. 1 and El. 4.

Pipe Roof Design.
e Gongbei tunnel was excavated using the conventional method with a pipe roof in combination with artificial ground freezing as presupport. e pipe roof was designed to support the Earth loads, and the artificial ground freezing was designed to seal the pipe roof to avoid water infiltration between pipes. e pipe roof consists of 36 steel pipes with a diameter of 1620 mm and length of 255 m. Besides the 36 roof pipes, one testing pipe (0 # pipe) was designed adjacent to 5 # pipe, as shown in Figure 2.
ese two pipes were installed first to test the jacking parameters. e wall thickness of the steel pipes is 2 Advances in Civil Engineering 20 mm for the lower 19 pipes and 24 mm for the upper 18 pipes. e net spacing between adjacent pipes is 35.5-35.8 cm.
e roof pipes were installed through microtunneling using slurry shield machines. Four 1667 mm diameter AVN1200TC machines purchased from Herrenk Necht in Germany were used to jack the pipes section-bysection with a sectional length of 4 m from a vertical shaft at one end of the tunnel. e pipes were received at a vertical shaft on the other end of the tunnel. Figures 3-5 show the design of the artificial ground freezing. Figure 3 shows that the freezing pipes were installed inside the steel roof pipes. e freezing system include three different pipes which are the circular freezing pipes, irregular shape freezing pipes, and circular heating pipes. Two circular freezing pipes were installed at the sidewall of the odd number roof pipes, which were filled with lightly expansive concrete, to circulate cold brine. A circular heating pipe was installed at the crown of the odd number roof pipes to circulate hot brine with temperature of 5-10°C to limit the thickness of frozen soil to avoid the excessive frost heave. Two irregular shape freezing pipes were installed at the sidewall of the even number pipes to circulate cold brine. e even number of roof pipes were filled with concrete after the excavation was completed, and the permanent liner was installed. Figure 4 shows that, along tunnel longitudinal length, three independent freezing loops were designed with a round length of about 85 m. Figure 5 shows that the tunnel cross section was divided into 5 freezing zones, A, B1, B2, B3, and C. e frozen ring thickness of zones A, B1, and B2 was designed to be within 2 to 2.3 m and 2 to 2.6 m for zones B3 and C.

Artificial Ground Freezing Design.
To further control frost heave, the soils around the pipe roof were pregrouted to reduce the permeability. e thickness of soil grouting was designed as 2.5 to 3 m. e pregrout was conducted through the predrilled holes in the roof pipes. Figure 6 shows the arrangement of the grouting holes. Within each 4 m length pipe section, 28 grouting holes were predrilled along 4 equally spaced rolls. Each roll contains 7 holes distributed along the circumferential length of the pipe. e grouting pipes have a diameter of 32 mm and length of 1.5 m. e grout was made of microfine cement and bentonite with a water-cement ratio (W/C) of 1 : 1 and bentonite dosage of 3%. e grouting pressure of the upper, middle, and lower parts of the pipe roof was controlled at 0.05-0.15 MPa, 0.15-0.3 MPa, and 0.3-0.5 MPa, respectively. e grouting speed was controlled to 5-25 L/min at the early stage, 5-15 L/min at the middle stage, and less than 5 L/min at the end stage.

Ground Settlement Induced by Pipe Jacking
e construction of the Gongbei tunnel was started in Oct 2012 and completed in Aug 2018. e successful construction of the Gongbei tunnel provides valuable field measurement data on analysing the ground movement during the installation of the large-scale pipe roof and artificial ground freezing. To analyze the ground settlement induced by multiple parallel pipe jacking, the critical parameters, including ground loss volume ratio η and the settlement trough width coefficient i, are derived based on the field measurement. e pipe interaction during multiple parallel jacking is analyzed using analytical and numerical methods.

Field Measurement Analysis.
e jacking order of the pipes is labeled outside the pipes, as shown in Figure 2. e jacking process was divided into three stages, including stage  Advances in Civil Engineering I (0 # and 5 # testing pipes), stage II (29 # , 28 # , 9 # , and 10 # pipes), and stage III (the remaining pipes). During pipe jacking, the settlement at the ground surface was monitored. Figure 7 shows the layout of the monitoring points, including 389 points within 24 cross sections which were set along the tunnel alignment at an interval of 10 m. In each cross section, the ground settlement was measured with 2 m interval within tunnel excavation boundary and with 4, 6, 8, and 10 m interval outside the excavation boundary. To avoid the potential impact of the launching and receiving shafts, three monitoring sections CJ5, CJ15, and CJ21 in the middle of the tunnel alignment are selected for data analysis. e cumulative ground settlement at the three cross sections during the installation of the 37 pipes is shown in Figure 8. e positive direction of the x-axis represents north direction. During the stage I jacking, the ground settlement caused by jacking 0 # and 5 # pipe was about 2.7 mm and 2.2 mm, respectively, and the cumulative settlement was about 5 mm. During stage II, the installation of the 9 # , 10 # , 28 # , and 29 # pipe caused a cumulative increase of ground settlement of about 2 mm. During stage III, the maximum settlement first accumulated on the north side of the tunnel centerline and gradually approached to the tunnel centerline as the jacking process continued. At the end of jacking the 37 pipes, the ground settlement curve at each cross section was in an approximate normal distribution along horizontal direction with the maximum value occurring above tunnel centerline, and the inflection point was located at the position with a horizontal distance of about 12 m from the tunnel centerline. e maximum ground settlement at the three sections ranged from 21.4 to 22.2 mm. e ground settlement caused by single pipe jacking is calculated as the difference between the cumulative ground settlement at the end of the pipe jacking and that of the fore pipe jacking. Based on the field measurement data at the three sections, the settlement curve during each pipe jacking was fitted with Peck's equations: where S (x) is the ground settlement at a distance x from the tunnel axis, S max is the maximum ground settlement above the tunnel axis, V s � ηπR 2 is unit ground loss, in which η is the ground loss volume ratio, and i is the ground settlement trough width coefficient. i and S max obtained by curve fitting are substituted into equation (2) to back calculate the corresponding V s . en, the ground loss ratio η of single pipe jacking is derived according to V s � ηπR 2 . e width parameter k is calculated by where z is the burial depth of the roof pipe. Figure 9 presents the ground settlement trough width coefficient i, width parameter k, ground loss volume ratio η, and maximum ground settlement S max during jacking each pipe. e settlement trough width coefficient i ranges from 3.4 m to 15.4 m. e width parameter k ranges from 0.3 to 0.87 and decreases with the increase in pipe burial depth. e maximum ground settlement S max ranges from 0.2 to 2.7 mm and also decreases with increasing pipe burial depth. Figure 9(d) shows that the η value ranges from 0.3% to 1.5%. e η value of 0 # pipe jacking of three sections is about 1.5%. e average η value for the 37 pipes was about 0.8%. Based on the analysis, the ground settlement trough width coefficient i is closely related to pipe burial depth, which is consistent with the finding from literature [11,22]. In Figure 10, the relationship between the width coefficient i of each pipe and the pipe burial depth is plotted and is compared with empirical relationships. e width coefficient curve is generally in accordance with the empirical relationships. When the pipe burial depth is less than 15 m, the relationship is close to the Liu and& Hou's [23] relationship. It approaches to the O'Reilly and New's relationship when the pipe burial depth is deeper than 25 m. e regression function is i � 0.35z + 2.72, and the correlation coefficient R 2 is 0.86.

Pipe Interaction Analysis.
During the installation of pipe roof, the fore-jacked pipes have a reinforcement effect to the ground. e reinforcement increases the stiffness of the ground and reduces the ground loss during consequent Zone A

Zone B1
Zone B2 Zone B3 Zone C Figure 5: Zonation of artificial ground freezing within tunnel cross section. 4 Advances in Civil Engineering jacking [24,25]. To investigate the interaction of pipes during multiple parallel jacking, the ground settlement is estimated using analytical and numerical analysis, and the analysis results are compared with the field measurement. During analytical analysis, Peck's empirical equation is applied to single pipe jacking. e ground loss volume ratio is assumed as 1.5% which is obtained from the back calculation of settlement during jacking the first (0 # ) pipe where no pipe interaction exists, and the width coefficient i is calculated using the regression i � 0.35z + 2.72. e cumulative ground settlement is obtained from superposition of the settlement of each pipe. As such, no pipe interaction is considered in the analytical analysis. e numerical analysis is conducted using a two-dimensional plane stain model by finite element (FE) software Plaxis 2D. e ground loss volume ratio is assumed as 1.5% for all the pipes. e model was built with a width of 80 m and a depth of 38.5 m. As shown in Figure 11 e ground loss is simulated by reducing the diameter of the pipes, and the radial contraction rate is set to 1.5%. e small-strain soil hardening model (Hss) was applied to the soils. e shear strength parameters of the soil were determined based on the tri-axial consolidated drainage tests. e stiffness modulus is converted according to the Plaxis manual: E oed � E s , E 50 � E oed , E ur � 4E oed , and G 0 � 4E ur . Parameters used in numerical simulation are detailed in Tables 1-3. According to the jacking order, as shown in Figure 2, the pipe roof installation was simulated as follows: stress initialization ⟶ excavation and support of 0 # pipe ⟶ excavation and support of 5 # pipe ⟶ ······⟶ excavation and support of 20 # pipe. Figures 11(b)-11(d) show the vertical ground deformation at the end of the stage I (0 # and 5 # pipe jacking), stage II (9 # , 10 # , 28 # , and 29 # pipe jacking), and stage III (remaining pipes jacking). During stage I, the settlement is distributed symmetrically around 0 # and 5 # pipes, and the maximum settlement occurs directly above them. With the pipe jacking continuing during stage II, there is a further increase in the ground settlement above 0 # and 5 # roof pipes. e soil deformation caused by 9 # and 10 # pipe jacking is smaller than that caused by 28 # and 29 # pipe jacking in spite of the same burial depth.
at is mainly because the  Advances in Civil Engineering fore-jacked 0 # and 5 # roof pipes provide a reinforcement effect on soil. e maximum ground settlement is still above 0 # and 5 # roof pipes at the end of stage II. During stage III, the maximum ground settlement at the end of all 37 pipes jacking is above the centerline of the tunnel, and the ground settlement decreases with the increasing horizontal distance from the axis. Figure 12 shows the cumulative ground settlement at the end of three stages of pipe jacking. Since the ground settlement caused by single pipe jacking is calculated and    Advances in Civil Engineering simulated with the same ground loss volume ratio of 1.5%, the calculated values and simulated values of each stage represent the cumulative ground settlement without pipe interaction. e measured values account for the effect of pipe interaction. erefore, the pipe interaction can be analyzed by the comparison of the measured, simulated, and calculated ground settlement. Figure 12 shows that the simulation results are consistent with the calculation results.

Frost Heave Analysis
e artificial ground freezing is designed to provide water tightness between adjacent pipes. e major issue related to the artificial ground freezing is the frost heave, which, if not well controlled, may damage the adjacent structures and utilities.
e control of the frost heave consists of two aspects, limiting the thickness of the frozen soil to limit the in situ frost heave and reducing the water content and permeability of the soil by grouting to limit both the in situ and segregation frost heave.

Numerical Simulation.
e frost heave was analyzed using finite element (FE) software COMSOL Multiphysics. Based on the elastic constitutive relation of soil skeleton, the COMSOL applies a thermal-coupling model to calculate the frost heave considering the ice-water phase transition, water migration in the soil, and variation of the soil temperature during the freeze-thaw process. e basic assumptions of the model include the following: (1) the soils are isotropic elastic materials; (2) ice-water phase transition occurs at −5°C; (3) the heat conduction follows Fourier's Law; (4) the water migration takes place in the form of liquid water only and is driven by the gradient of unfrozen water following Darcy's law; (5) the soil skeleton is a low-temperature expansion material, and the expansion coefficient is defined as the  Figure 13. e model size is set as 150 × 50 m and is meshed freely by triangular elements. e minimum element size is set to 1 cm. Total of 69954 elements are established with the average mesh quality of 0.92. e subsurface strata consist of 7 layers. e groundwater table is set at 1 m below the ground surface. e mechanical properties of soils obtained from lab tests are listed in Table 4. e thermo-physical properties of soils before and after grouting are listed in Table 5. e grouted soils were obtained through field samples. In total, 72 holes were drilled which account for 8% of 900 grouting holes. e density, water content, strength, and frost heave ratio of the samples were determined by laboratory testing. e on-site hydraulic test was conducted in the drill holes to determine the soil permeability after taking samples. e physico-mechanical parameters of soils before and after grouting are detailed in Table 4. In general, the porosity of the soils decreases by approximately 13% from 0.26-0.67 to 0.22-0.60. e permeability coefficient decreases from 10 −2 -10 −7 cm/s to 10 −5 -10 −7 cm/s. e elastic modulus increases by 118-1156% from 1.6-15.0 MPa to 20.1-32.7 MPa after grouting. According to the result of freezing tests, the specific heat capacity of soils decreases by around 5%, and the frost heave ratio decreases by 12% from 3.82%-6.63% to 3.53%-5.46%. e ground surface temperature is set to 21°C. e temperature of circular freezing pipes, irregular-shaped freezing pipes, heating pipes, and empty roof pipes are set to −28°C, −24°C, 8°C, and −10°C, respectively. e soil temperature considering water migration is calculated by the coefficient differential equation PDE module. e ground deformation is calculated using the solid mechanics module.

Result Analysis.
Four scenarios were simulated. ey are without frost control measures, control scheme I: activating heating pipes, control scheme II: pregrouting, and control scheme III: the combination of pregrouting and activating heating pipes. Figure 14 shows the ground heave of four scenarios after freezing of 30, 50, and 70 days (d).
e ground heave increases with the freezing time with the development of frozen soil. e maximum heave occurs above the tunnel centerline. e ground deformation curves of scheme I, II, and III are all below that without control measures. e three schemes are beneficial to frost heave control, and scheme III is most effective. After freezing for 70 d, the maximum heave without control measures and with scheme I, II, and III control measures is 13.5, 12.3, 9.0, and 7.8 cm, respectively. e maximum ground heave decreases 9% with scheme I, 33% with scheme II, and 42% with scheme III. e decrease of maximum heave of scheme II and III is 4.5 cm and 5.8 cm respectively, which are 2.75 and 3.83 times of 1.2 cm, and the decrease of heave of scheme I indicates that pregrouting plays a critical role in the frost heave control. is is because the pregrouting forms a grout-soil composite which has a lower permeability, thus hindering the water  Note. B represents soil parameters before grouting; A represents soil parameters after grouting.   Advances in Civil Engineering migration in the soil and reducing the segregation frost heave. Limiting the frozen soil thickness using heating pipes serves a supplement to the heave control. e combination of pregrouting and heating pipes has a better control effect than any single measure.

Field Measurement Analysis.
e field temperature realtime monitoring system contains total of 32 monitoring sections along the 255 m length tunnel alignment. Figure 15 shows the layout of temperature monitoring points of each monitoring section. ere are 7 monitoring points A1-A7 at the inner wall of odd number roof pipes and 6 monitoring points B1-B6 at the inner wall of even number roof pipes in each monitoring section. Additionally, monitoring points are also installed in the soils around the pipe roof to collect the soil temperature. e monitoring section CJ21 is selected to analyze the temperature development during freezing. As shown in Figure 16, the simulation results match the field measurement well, and the frozen soil develops with the time. After freezing of 70 d, the temperature of the soil above and beneath the tunnel is around −5°C and −8°C, respectively. e frozen soil near freezing pipes develops well with the lowest soil temperature of −28°C. Near the heating pipes, the soil temperature is higher and the isotherm is denser due to the circulation of hot brine. e increase in the temperature of the soil near heating pipes breaks the balance of soil temperature field, relieves the development of frozen soil, and reduces the frost heave. Since the soil around pipes 14#-21# was not grouted, the frozen soil at the bottom of pipe roof develops rapidly and the frozen area is larger. Figure 17 shows the comparison of ground heave between numerical simulation and field measurement of section CJ21. In general, the numerical simulation results are consistent with those of field measurement. e simulation results of the maximum heave after freezing of 30, 50, and 70 d are, respectively, 7.0, 7.5, and 7.8 cm, respectively, compared with of 6.6, 7.3, and 7.6 cm from the field measurement. Figure 18 shows the total vertical ground deformation during the pipe roof installation and ground freezing. is figure shows that the maximum ground settlement and heave occur above the tunnel axis and the ground settlement and heave decrease with the distance from the axis. e maximum ground settlement and heave are around 2.2 cm and 7.8 cm, respectively. e total vertical ground deformation is calculated by the sum of ground settlement and heave. After part of the ground heave is offset by the ground settlement, the ground surface is heaved up by an average of 3.7, 5.2, and 5.6 cm, respectively, after freezing 30, 50, and 70 d.

Conclusions
is paper analyzes the vertical ground movement during large-scale pipe roof installation and artificial ground freezing of Gongbei tunnel of the Hong Kong-Zhuhai-Macau bridge project. e ground settlement during pipe roof installation and the ground heave during artificial ground freezing are analyzed analytically and numerically and are compared with the field measurements. Findings from the study are: (1) e ground settlement during pipe roof installation is relatively small. e maximum ground settlement during single pipe jacking ranges from 0.2 to 2.7 mm dependent on the pipe burial depth. e maximum cumulative ground settlement after jacking 37 pieces of 1620 mm steel pipes is 2.2 cm. (2) e ground settlement during pipe jacking is in an approximately parabola distribution in horizontal direction. e width coefficient of settlement trough during single pipe jacking ranges from 3.4 to 15.4 m. e width parameter ranges from 0.31 to 0.87 and decreases with the increase of pipe burial depth. e initial ground loss ratio during pipe jacking is 1.5%. Due to the reinforcement on the ground provided by the fore-jacked pipes, the average ground loss ratio at the end of all 37 pipe jacking is close to 0.8%.
(3) e artificial ground freezing generates a relatively large ground heave with a maximum value of 7.8 cm.
Pregrouting plays a critical role in the frost heave control by reducing the heave by about 33%. Limiting the frozen soil thickness by heating pipes serves as an effective supplement to frost heave control by reducing the heave by about 9%. e combination of the two measures reduces the ground heave by about 42%.

Data Availability
e data used to support the findings of this study are included within the article.   14 Advances in Civil Engineering

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.