Factor’s Influence Analysis of Frost Heave during the Twin-Tunnel Construction Using Artificial Horizontal Ground Freezing Method

Hong, Rong-bao; Cai, Hai-bing; Li, Meng-kai; Li, Ya-ru

doi:https://doi.org/10.1155/2020/6650182

Advances in Civil Engineering

On this page

Abstract Introduction Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Advancements in Tunnel Construction

View this Special Issue

Research Article | Open Access

Volume 2020 | Article ID 6650182 | https://doi.org/10.1155/2020/6650182

Factor’s Influence Analysis of Frost Heave during the Twin-Tunnel Construction Using Artificial Horizontal Ground Freezing Method

Rong-bao Hong,¹Hai-bing Cai,¹Meng-kai Li,¹and Ya-ru Li¹

Academic Editor: Ma Jianjun

Received07 Oct 2020

Revised05 Nov 2020

Accepted18 Nov 2020

Published11 Dec 2020

Abstract

In order to analyze the influence of different twin-tunnel parameters on the frost heave of the ground, six tunnel clear distances (0.25 D, 1.00 D, 1.50 D, 2.00 D, 2.50 D, and 3.00 d), three tunnel buried depths (6 m, 12 m, and 18 m), and two freezing orders (simultaneous freezing and sequential freezing) are selected to establish the two-dimensional numerical calculation models, using ABAQUS finite-element program and the user subroutine of frost heaving deformation, and considering the orthotropic deformation characteristics of frozen soil. Numerical simulation results show that the interaction of twin-tunnel freezing is weakened with the increase in tunnel clear distance. Correspondingly, the heaving displacement of the ground surface also decreases. Besides, the heaving displacement curve of the ground surface gradually changes from the unimodal N-type to bimodal M-type as the tunnel clear distance increases. However, the trough of the bimodal M-type curve gradually disappears as the tunnel buried depth increases. SPSS mathematical analysis results show that tunnel clear distance has the highest significance, tunnel buried depth ranks the second, and freezing order ranks the last. It is important to select the appropriate tunnel clear distance during the design of the twin-tunnel construction using the artificial horizontal ground freezing method.

1. Introduction

The artificial horizontal ground freezing method has been widely applied to subway tunnel construction due to its temporary support capacity and sealing capacity. At present, many successful cases of artificial freezing tunnel construction are available for research. In China, the geological survey report of Tianhe passenger station of Guangzhou metro line 3 shows that the tunnel construction section passes through the sandy clay layer, which has poor compactness and rich water that often lead to quicksand and landslides [1]. To this end, the solution that combines the full section horizontal freezing method with CRD excavation method [2] was proposed and successfully applied to solve this engineering problem. The same situation also occurred in the tunnel project from Minzu University station to convention center station of Nanning rail transit line 1. The stratum of tunnel crossing is round gravel layer, which has large permeability that often leads to water flush accident, and the artificial horizontal freezing method works perfectly [3]. Hong Kong-Zhuhai-Macau Bridge, which is the key construction project in China, had attracted the attention of domestic and foreign engineers. Gongbei tunnel, an important part of its connecting line, also faced severe geological conditions. Many researchers [4–7] put forward different design schemes from the aspects of controlling the ground deformation, limiting the flow direction of groundwater and protecting the surrounding buildings. At last, the solution that combines the pipe curtain method with artificial freezing method [8] had been adopted with various efforts and made this tunnel engineering a good quality to resist natural disasters [9].

In view of this, the numerical simulation prediction accuracy of the ground frost heave has been further improved in the subway tunnel construction using artificial freezing technology based on the above successful cases. Konrad [10] first proposed the frost heave of segregation potential model and developed a two-dimensional numerical calculation program (SP-2D). This model revealed the frost heave law of the ground during the tunnel freezing construction of Kobe in Japan. Cheng and Zhang [11] used ADINA finite-element program to study the coupling law of the ground temperature and displacement based on a tunnel freezing project. Qiao and Tao [12] used FLAC finite-element program to study the ground displacement in the process of tunnel freezing construction. Yang et al. [13] used MARC finite-element program to establish a two-dimensional moisture-heat-stress model, which simulated the ground frost heaving displacement law in the connecting passage construction of Nanjing metro line 1. Yan and Fang [14] used GEO-SLOPE program to simulate the freezing and excavation process of connecting passage, based on the coupled thermo-solid theory considering the phase transition. Gao et al. [15] used COMSOL finite-element program to analyze the distribution of the temperature field, seepage field, and displacement field in the horizontal freezing engineering of connecting passage. Cai et al. [16] used ABAQUS finite-element program combining the subroutine of the frost heaving deformation to establish the thermal mechanical coupling numerical analysis method of the ground displacement during horizontal freezing of tunnel. The results of the above numerical calculation models are all consistent with the actual project, which also shows that the numerical simulation is an effective means to predict the ground change law in the tunnel freezing construction.

Based on the above point of view, it is meaningful to research the frost heaving phenomenon caused by the different twin-tunnel parameters during the horizontal freezing process of the twin tunnels through ABAQUS finite-element program, and the frost heaving law obtained by the numerical simulation will be of great guiding significance to the twin-tunnel design scheme using the artificial freezing technology.

2. Numerical Simulation of Frost Heave

2.1. Mathematical Model of Freezing Temperature Field

Freezing temperature field of tunnel construction is a transient heat conduction problem with soil-water phase change. According to heat transfer theory [17, 18], the governing partial differential equation for 2D freezing temperature can be written aswhere is called the equivalent volumetric heat capacity (kcal/(m³· C)), T is the soil temperature (°C), t is time (day), and is called the equivalent thermal conductivity (kcal/(m·d·°C)).

and are also given bywhere C_f and C_u denote the volumetric heat capacity (kcal/(m³·°C)) of frozen soil and unfrozen soil, and C_f = ρ_f·c_f and C_u = ρ_u·c_u, where ρ_f and ρ_u are the density of the frozen soil and unfrozen soil (kg/m³), c_f and c_u are the mass specific heat of the frozen soil and unfrozen soil (kcal/(kg·°C)). k_f and k_u are, respectively, the thermal conductivity (kcal/(m·d·°C)) of the frozen soil and unfrozen soil. T_d is the freezing temperature (°C) of soil; T_r is the thawing temperature (°C) of soil. L is the latent heat of phase change for one cubic meter of soil (kcal/m³).

2.2. Initial and Boundary Conditions of Freezing Temperature Field

The initial condition of freezing temperature field iswhere T₀ is the initial temperature (°C) of soil, assuming that the formation temperature distribution is uniform.

In the calculation model, each freeze pipe could be simplified into one point, so the boundary condition for each freeze pipe is expressed aswhere (x_p, y_p) is the Cartesian coordinate of each point (freeze pipe) and T_c (t) is the surface temperature (°C) of the freezing pipes.

The boundary condition of soil in far field is given by

On the surface of the ground, the boundary condition of convective heat transfer between atmosphere and soil is given bywhere α₁ is the convective heat transfer coefficient between the atmosphere and soil (kcal/(m²·d·°C)), T_a is the atmospheric temperature (°C), and n₁ is the normal direction vector of the surface.

2.3. Orthotropic Frost Heave Deformation of Frozen Soil

Existing studies show that [19–22] there is a big difference in the mechanical properties between the freezing front and the heat flow direction during the freezing of soil, and the frost heave deformation of frozen soil exhibits orthotropic characteristic. This means that the heat flow direction is the main principle direction, while the other two main directions are perpendicular to the heat flow direction, and the main frost heave deformation occurs along the heat flow direction.

When the temperature of soil reaches the freezing temperature (freezing point), the frost heave ratio of soil, which can be obtained from the indoor frozen soil test, is defined as the ratio of a volume increase to the initial volume [23]. Therefore, the frost heave ratio can be regarded as the temperature-dependent instantaneous volumetric strain, which is symbolized by here, and the theoretical calculation method is shown as follows:where is the frost heave ratio of soil, is the longitudinal frost heave of sample, and is the original height of sample.

As shown in Figure 1, in the Cartesian coordinate system (x, y, z), a local coordinate system (1, 2, 3) is established along the heat flow direction, and 1, 2, and 3 represent the three principal axes of the local coordinate system, respectively.

When the soil temperature is colder than the freezing point, the total strain of the frozen soil consists of two parts, namely, the strain caused by load and the instantaneous volumetric strain caused by temperature. In the local coordinate system, the instantaneous thermal strain components can be obtained by the following equation:where η is an orthotropic parameter. When η = 1/3, the frost heave deformation of the frozen soil exhibits isotropic characteristic; when η = 1, the frost heave deformation of frozen soil exhibits unidirectional characteristic, and its direction is along the temperature gradient.

In the local coordinate system, the total strain component of frozen soil is calculated according to the elastic Hooke theorem, as follows:where E is elastic modulus of frozen soil and μ is Poisson’s ratio of frozen soil.

For the plane strain problem, there is , so that

By substituting equation (10) into (9), the total strain component of the frozen soil in the local coordinate system can be expressed as

The instantaneous thermal strain components can be written as

As shown in Figure 1, the heat flow angle θ is the angle between the local coordinate system and the Cartesian coordinate system. The instantaneous thermal strain components in the Cartesian coordinate system can be calculated by

Substituting equation (12) into (13), we can get

Therefore, considering orthotropic frost heave deformation of frozen soil, the frost heave strain components (i.e., the instantaneous thermal strain components) can be obtained by equation (14) when the soil temperature is colder than the freezing point.

2.4. User Subroutine of Freeze-Thaw Deformation of Frozen Soil

The developed frost heave model considering orthotropic frost heave deformation of frozen soil is implemented by ABAQUS software. ABAQUS employs user-defined subroutines to extend the capabilities of the applications. The developed frost heave model is programmed by three user subroutine interfaces (USDFLD, UEXPAN, and GETVRM). The flowchart of the user-defined subroutine is shown in Figure 2.

At each incremental step, GETVRM is called by USDFLD and the heat flow components of the integration point are used as state variables and transmitted to UEXPAN. In the user subroutine interface UEXPAN, the heat flow angle in Figure 1 is calculated according to the heat flow components, and then the thermal strain components of each integration point are calculated according to equation (14). According to temperature determination of integral points, the frost heave deformation can be determined whether it occurs or not at the current incremental step.

In this paper, ABAQUS finite-element program and user subroutine are used to simulate the frost heave of the ground during the tunnel freezing construction. The analysis type is coupled temperature/displacement analysis. The main program of numerical simulation can simulate the whole construction process of tunnel freezing, excavation, and thawing.

3. Numerical Calculation Model

In order to verify the reliability of the numerical simulation method proposed in this paper, Cai [1] applied the method to the tunnel freezing project of Guangzhou Metro Line 3. The soil parameters input in ABAQUS finite-element software were obtained from the test report of Guangzhou artificial frozen soil. Eventually, the frost heaving displacement of the ground surface acquired from the numerical simulation is in good agreement with that of the field measurement. In view of this, using the same model and soil parameters to analyze the influence of tunnel parameters on the frost heaving displacement of the ground surface will produce certain engineering application value.

3.1. Implementation of the Frost Heave Model

The length (x-axis) of the model is set as 80 m, and the width (y-axis) of the model is set as 50 m for eliminating the boundary effect according to modeling rules proposed by Cai et al. [24]. In addition, the cross section of each tunnel is circular with the outside diameter (D) of 6 m, and the centers of the two tunnels are on the same horizontal line. Correspondingly, the radius of horizontal freezing pipes layout circle (R_d) is set as 4m in order to form a frozen wall with an effective thickness of 2 m after freezing 95 days.

In this paper, three tunnel buried depths (6 m, 12 m, and 18 m), six tunnel clear distances (0.25 D, 1.00 D, 1.50 D, 2.00 D, 2.50 D, and 3.00 D), and two freezing orders (simultaneous freezing and sequential freezing) are selected to establish two-dimensional models, and as a result, 18 models are established through a random combination between the tunnel clear distances and buried depths, which are shown in Table 1. In model 6-1, model 12-1, and model 18-1, the tunnel clear distance is 0.25 D, namely, 1.5 m, which is less than the effective thickness of frozen wall (2 m); therefore, three freezing pipes are shared in the middle of the two tunnels, and the total freezing pipes are 37, as shown in Figure 3(a). Meanwhile, the total freezing pipes are 40 in the remaining models, as shown in Figure 3(b).

(a)

(b)

Finally, the ground frost heaving phenomenon is analyzed under the 18 conditions by ABAQUS finite-element program and the user subroutine of the frost heaving deformation based on the above two-dimensional geometric models. Due to the large number of models, only two representative meshing diagrams are shown here. Model 6-1 consists of 5384 elements shown in Figure 4(a), while model 12-3 consists of 6005 elements shown in Figure 4(b). The quadrilateral element (CPE4RT) with good stability is selected as the element type.

(a)

(b)

3.2. Calculation Method of Thermodynamic Parameters

Unfrozen water content is not only an important index to evaluate the characteristics of water migration in frozen soil, but also a commonly used parameter in the thermal calculation of frozen soil. Xu et al. [25] put forward a rapid method for determining unfrozen water content, who thinks that there is a dynamic equilibrium relation between the initial water content of soil and the freezing temperature, as shown in the following equation:where is the initial water content of soil (%), is the freezing temperature of soil (°C), and a and b are the empirical constants related to the soil properties.

Research [26] shows that the empirical constant b can be determined when the soil properties and the salt content of the soil are determined. Ulteriorly, Xu and Deng [27] put forward that the unfrozen water content can be calculated as follows:where is the water content of unfrozen soil (%), and T is the temperature of frozen soil (°C).

The specific heat calculation equation of frozen soil is as follows:where is the specific heat of soil particles, it is taken as 0.18 kcal/(kg·°C) in this paper, is the specific heat of ice, it is taken as 0.5 kcal/(kg·°C) in this paper, and is the specific heat of water, it is taken as 1.0 kcal/(kg·°C) in this paper.

The calculation equation of thermal conductivity of frozen soil is as follows:where Q is the heat quantity, is the temperature gradient (°C/m), is the area (m²), and t is the time.

The calculation equation of latent heat of frozen soil is as follows:where is the latent heat of phase change per unit mass of water, it is taken as 79.6 kcal/kg in this paper, and is the dry density of soil (kg/m³).

The calculation results of thermophysical parameters of soil are shown in Table 2 based on the above calculation equations.

The freezing temperature of soil is often colder than 0°C as the soil contains a certain amount of salt. Hou et al. [28] had effectively identified the variation characteristics of soil temperature, moisture, and salinity during freezing through a laboratory unidirectional freezing experiment. The salt in soil inhibited water migration, which obviously reduced the freezing point of soil. In this paper, the phase transition temperature range is considered to reflect this process. The freezing point of the soil selected is −0.5°C referring to the frozen soil test report [29]. The values of “Solidus Temp” and “Liquidus Temp” [30] are needed to be input in ABAQUS calculation program. Therefore, the Solidus Temp is set as −0.5°C, and the Liquidus Temp is set as 0°C in the simulations.

The mechanics parameters [27] of soil are shown in Table 3. T is the negative temperature of soil, and the frost heave ratio of soil is 1%. Besides, the freeze-thaw deformation characteristics of orthotropic soil are considered during the process of numerical calculation, where the deformation characteristic coefficient is 0.9 [31].

3.3. Initial and Boundary Conditions

The initial temperature of soil is set as 20°C according to the field measured ground temperature report [29]. The “ammonia-brine” cycle refrigeration system is widely adopted in practical projects, which can stabilize the temperature of brine at −25°C [32]. Therefore, the surface temperature of the freezing pipes is assumed to be −25°C. Besides, the temperature of frozen soil is much colder than that of air, which will inevitably lead to heat exchange between soil and air. Therefore, the convective heat transfer coefficient between soil and air is assumed to be 175 kcal/(m²·d·°C) [13].

For the bottom boundary of the calculation model, vertical and horizontal displacements are constrained. Horizontal displacement is not allowed on the left and right boundaries of calculation model. Only the upper boundary of the calculation model is free [24].

3.4. Analysis of Calculation Results

3.4.1. Effect of Tunnel Clear Distance

(1) Heaving Displacement Field of the Ground. After simultaneous freezing for 95 days, the distributions of the vertical displacement field are shown in Figure 5 with the tunnel buried depth of 6 m.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5 shows that the vertical displacement field is symmetrically distributed along the center line of twin tunnels after simultaneous freezing for 95 days. In addition, the soil at the bottom of tunnel produces the downward frost heave displacement, while the soil at the top of tunnel produces the upward frost heave displacement. The vertical displacement isopleth of soil inside tunnel is distributed horizontally. It is caused by the fact that the soil inside tunnel is squeezed by the uniform frost heaving force in the movement of the frozen internal front. With the tunnel clear distances increasing from 0.25 D to 3.00 D, the maximum vertical displacements of the tunnel top soil are 60.5, 61.5, 56.3, 54.2, 52.2, and 50.9 mm, while the maximum vertical displacements of the tunnel bottom soil are −11.8, −11.2, −11.9, −12.9, −14.0, and −14.4 mm, respectively.

(2) Heaving Displacement of the Ground Surface. The heaving displacements of the ground surface with three tunnel buried depths are shown in Figure 6, which directly reflect the frost heaving phenomenon caused by the change of tunnel clear distance during freezing. Figure 6 shows that there are two types of the frost heave curves of the ground surface, which are unimodal N-type and bimodal M-type. The number of bimodal M-type curves decreases gradually with the increase of the tunnel buried depth under six tunnel clear distances; correspondingly, the number of the unimodal N-type curves increases gradually. Besides, the frost heave curves of the ground surface are all unimodal N-type under three tunnel buried depths as the tunnel clear distance is 0.25 D. Finally, the maximum vertical displacement of the ground surface decreases with the increase of the tunnel clear distance from 1.00 D to 3.00 D, which all appear at the center line of the tunnel.

(a)

(b)

(c)

3.4.2. Effect of Tunnel Buried Depth

(1) Heaving Displacement Field of the Ground. After simultaneous freezing for 95 days, the distributions of the vertical displacement field are shown in Figure 7 with the tunnel clear distance of 1.00 D.

(a)

(b)

(c)

Figure 7 shows that they have the same laws that the vertical displacement field is symmetrically distributed along the center line of the twin tunnels after simultaneous freezing for 95 days. In addition, the soil at the bottom of the tunnel produced the downward frost heave displacement, while the soil at the top of the tunnel produced the upward frost heaving displacement. The vertical displacement isopleth of the soil inside the tunnel is distributed horizontally. It is caused by the fact that the soil inside the tunnel is squeezed by the uniform frost heaving force in the movement of the frozen internal front. With the tunnel buried depths increasing from 6 m to 18 m, the maximum vertical displacements of the tunnel top soil are 61.5, 64.2, and 71.0 mm, while the maximum vertical displacements of the tunnel bottom soil are −11.2, −19.8, and −26.2 mm, respectively.

(2) Heaving Displacement of the Ground Surface. The heaving displacements of the ground surface with six tunnel clear distances are shown in Figure 8, which directly reflect the frost heaving phenomenon caused by the change of tunnel buried depth during freezing. Figure 8(a) shows that the frost heave curves of the ground surface are all unimodal N-type under the three tunnel buried depths as the tunnel clear distance is 0.25 D. Meanwhile, the vertical displacement of the ground surface decreases with the increasing distance from the center line of twin tunnels. Besides, it can be seen from the other figures that the slope of the rising section of the vertical displacement curve of the ground surface decreases as the tunnel buried depth increases from 6 m to 18 m. It is easy to be found that the vertical displacement curve of the twin-tunnel freezing is composed of two normal curves [33].

(a)

(b)

(c)

(d)

(e)

(f)

3.4.3. Effect of Freezing Order

(1) Heaving Displacement Field of the Ground. The distributions of the vertical displacement field are shown in Figure 9 after sequential freezing for 190 days.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 9 shows that the vertical displacement field is no longer symmetrically distributed along the center line of twin tunnels after sequential freezing for 190 days. It turns to be that the displacement nephogram area of the right half of the model is larger than the left half of the model. It needs to be declared here that the freezing pipes of the left line tunnel have been first opened for 95 days of freezing, and immediately, the freezing pipes of the right line tunnel continued to be opened for 95 days of freezing in the process of numerical simulation. With the tunnel clear distances increasing from 0.25 D to 3.00 D, the maximum vertical displacements of the tunnel top soil are 63.34, 63.79, 56.39, 53.09, 51.88, and 51.05 mm, while the maximum vertical displacements of the tunnel bottom soil are −13.10, −11.46, −11.92, −12.91, −13.84, and −14.65 mm, respectively.

(2) Heaving Displacement of the Ground Surface. The heaving displacement curves of the ground surface are shown in Figure 10 after simultaneous freezing for 95 days and sequential freezing for 190 days, respectively. Figure 10 shows that the vertical displacement curve of the ground surface after simultaneous freezing for 95 days is included in the vertical displacement curve of the ground surface after sequential freezing for 190 days when the tunnel clear distance is 0.25 D and 1.00 D. However, two curves are basically coincident as the tunnel clear distance increases from 1.50 D to 3.00 D. Besides, the vertical displacement of the ground surface at the center line of the left line tunnel is slightly less than that of the right line tunnel. Therefore, regardless the perspective of the construction period or the vertical displacement of the ground surface, simultaneous freezing is more economical and safer than sequential freezing [34].

(a)

(b)

(c)

(d)

(e)

(f)

3.4.4. Factor Significance Analysis

(1) Factor Influence Level. The vertical displacement of the ground surface is accumulated by the frost heave of the ground [35], which is one of the items that must be tested in the field. To this end, the vertical displacements of the ground surface corresponding to six tunnel clear distances, three tunnel buried depths, and two freezing orders have been imported to SPSS mathematical analysis program for the factor influence levels analysis, and the analysis results can be seen in Table 4.

It is assumed that the tunnel clear distance, tunnel buried depth, and freezing order are all significantly impacting the vertical displacement of the ground surface; namely, the value of significance is less than 0.05 [36]. Table 4 shows that the significance values of the tunnel clear distance, tunnel buried depth, and freezing order are 0, 0.069, and 0.410. Only the significance value of the tunnel clear distance is less than 0.05, which means that only tunnel clear distance has accepted the original hypothesis. Finally, the influence level of the three factors ranks as follows: tunnel clear distance > tunnel buried depth > freezing order.

(2) Post Hoc Test. Six tunnel clear distances (0.25 D, 1.00 D, 1.50 D, 2.00 D, 2.50 D, and 3.00 D) and three tunnel buried depths (6 m, 12 m, and 18 m) are further analyzed by a post hoc test, which can be seen in Tables 5 and 6. It needs to be explained here that the number of the subfactors of the freezing order is less than 3, which is not satisfied to the condition of the post hoc test. Cai et al. [24] analyzed the freezing sequence of the uplink tunnel and downlink tunnel by model test and numerical simulation based on the actual twin-tunnel freezing project of the Shanghai Metro and found that the difference of the maximum ground surface frost heaving displacement in simultaneous freezing mode and sequential freezing mode is only 15%. Compared with the simultaneous freezing mode, the sequential freezing method extends the freezing time and increases the cost, but it can prevent frost heave damage and ensure project safety. So, the freezing order is not discussed below again.

Table 5 shows that the influence level of the tunnel clear distance ranks as follows: 3.00 D < 2.50 D < 2.00 D < 1.50 D < 0.25 D < 1.00 D, which shows that the tunnel clear distance of less than or equal to 1.00 D should be avoided as far as possible during the design process of horizontal freezing engineering of twin tunnels. Table 6 shows that the influence level of the tunnel buried depth ranks as follows: 18 m < 12 m < 6 m, which shows that the frost heave of the ground surface caused by the artificial horizontal freezing method in the construction process of shallow buried twin tunnels cannot be ignored.

4. Conclusions

In this paper, six tunnel clear distances, three tunnel buried depths, and two freezing orders have been selected to study the ground deformation laws during the construction process of twin tunnels using the artificial horizontal freezing method, and the conclusions are as follows:(1)The ground displacement field is symmetrically distributed along the center line of the twin tunnels in the process of simultaneous freezing, while it is not symmetrically distributed in the process of sequential freezing, using ABAQUS finite-element program and the user subroutine of frost heave deformation.(2)The vertical displacement curve of the ground surface gradually changes from unimodal N-type to bimodal M-type as the tunnel clear distance increases; however, the trough of the bimodal M-type curve gradually disappears as the tunnel buried depth increases.(3)The influence level of the three factors ranks as follows: tunnel clear distance > tunnel buried depth > freezing order. It is also recommended that the tunnel clear distance of less than or equal to 1.00 D should be avoided as far as possible during the design process of the horizontal freezing engineering of twin tunnels, and it is necessary to strengthen the monitoring of the frost heave of the ground surface in the shallow buried twin tunnels.(4)Simultaneous freezing is more economical and safer than sequential freezing in the construction of twin tunnels.

Data Availability

The datasets generated and analyzed during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant no. 51778004), Research Activities Funding for Reserve Candidate of Academic and Technical Leaders of Anhui Province, China (Grant No. 2018H170), and Academic Funding for Top-notch Talents in University Disciplines (Majors) of Anhui Province, China (Grant No. gxbjZD10).

References

H. Cai, Prediction Methods and Their Engineering Application of Ground Frost Heave and Thawing Settlement in Subway Tunnel Horizontal Freezing Construction, Central South University, Changsha, China, 2012.
Z. Yao, H. Cai, H. Cheng et al., “Construction technique for the shallow-buried subway tunnel with large-section by long-distance horizontal ground freezing method and undermining method,” China Railway Science, vol. 32, no. 1, pp. 75–80, 2011.
View at: Google Scholar
Z. Jin, Research and Application on Horizontal Freezing Method in Connected Aisle of Nanning Rail Transit Line 1, Guangxi University, Guangxi, China, 2017.
P. Zhang, B. Ma, and W. Zhao, “The largest curved pipe roofing tunnel project in the world,” in Proceedings of the Pipelines 2013, pp. 953–963, Fort Worth, TX, USA, June 2013.
View at: Google Scholar
P. Zhang, B. Ma, C. Zeng, H. Xie, X. Li, and D. Wang, “Key techniques for the largest curved pipe jacking roof to date: a case study of Gongbei tunnel,” Tunnelling and Underground Space Technology, vol. 59, pp. 134–145, 2016.
View at: Publisher Site | Google Scholar
Y. Kang, Q. Liu, Y. Cheng, and X. Liu, “Combined freeze-sealing and new tubular roof construction methods for seaside urban tunnel in soft ground,” Tunnelling and Underground Space Technology, vol. 58, pp. 1–10, 2016.
View at: Publisher Site | Google Scholar
H. Zhang, P. Zhang, W. Zhou, S. Dong, and B. Ma, “A new model to predict soil pressure acting on deep burial jacked pipes,” Tunnelling and Underground Space Technology, vol. 60, pp. 183–196, 2016.
View at: Publisher Site | Google Scholar
X. Hu, S. Deng, and H. Ren, “In situ test study on freezing scheme of freeze-sealing pipe roof applied to the Gongbei tunnel in the Hong Kong-Zhuhai-Macau bridge,” Applied Sciences, vol. 7, no. 1, p. 27, 2017.
View at: Google Scholar
J. Li, Z. Li, and X. Hu, “Analysis on water sealing effect of freezing-sealing pipe roof method,” Chinese Journal of Underground Space and Engineering, vol. 11, no. 3, pp. 751–758, 2015.
View at: Google Scholar
J.-M. Konrad, “Prediction of freezing-induced movements for an underground construction project in Japan,” Canadian Geotechnical Journal, vol. 39, no. 6, pp. 1231–1242, 2002.
View at: Publisher Site | Google Scholar
H. Cheng and H. Zang, “Quasi-coupling numeric analysis of frost heave effect of horizontal freezing in artificial ground,” Chinese Journal of Geotechnical Engineering, vol. 25, no. 1, pp. 87–90, 2003.
View at: Google Scholar
J. Qiao and L. Tao, “Numerical simulation analysis of stress and displacement field in horizontal local freezing construction of subway tunnel,” Railway Engineering, vol. 2, pp. 39–41, 2004.
View at: Google Scholar
P. Yang, J. Ke, J. Wang et al., “Numerical simulation of frost heave with coupled water freezing, temperature and stress fields in tunnel excavation,” Computers and Geotechnics, vol. 33, no. 6, pp. 330–340, 2006.
View at: Publisher Site | Google Scholar
Q. Yan and X. Fang, “Thermal-solid coupling analysis of the connected cross aisle in metro constructed by horizontal freezing in ground,” China Railway Science, vol. 33, no. 1, pp. 54–59, 2012.
View at: Google Scholar
J. Gao, M. Feng, and Q. Gao, “THM coupling analysis of connected aisle in metro construction by artificial freezing method,” Journal of Glaciology and Geocryology, vol. 35, no. 4, pp. 904–912, 2013.
View at: Google Scholar
H. Cai, H. Cheng, and Z. Yao, “Numerical analysis of ground displacement due to orthotropic frost heave of frozen soil in freezing period of tunnel,” Chinese Journal of Rock Mechanics and Engineering, vol. 34, no. 8, pp. 1667–1676, 2015.
View at: Google Scholar
Y. Lai, S. Liu, Z. Wu, and W. Yu, “Approximate analytical solution for temperature fields in cold regions circular tunnels,” Cold Regions Science and Technology, vol. 34, no. 1, pp. 43–49, 2002.
View at: Google Scholar
H. Song, H. Cai, Z. Yao, C. Rong, and X. Wang, “Finite element analysis on 3D freezing temperature field in metro cross passage construction,” Procedia Engineering, vol. 165, pp. 528–539, 2016.
View at: Publisher Site | Google Scholar
J.-M. Konrad and M. Shen, “2-D frost action modeling using the segregation potential of soils,” Cold Regions Science and Technology, vol. 24, no. 3, pp. 263–278, 1996.
View at: Publisher Site | Google Scholar
K. Kim, W. Zhou, and S. L. Huang, “Frost heave predictions of buried chilled gas pipelines with the effect of permafrost,” Cold Regions Science and Technology, vol. 53, no. 3, pp. 382–396, 2008.
View at: Publisher Site | Google Scholar
H. Zheng, S. Kanie, F. Niu, S. Akagawa, and A. Li, “Application of practical one-dimensional frost heave estimation method in two-dimensional situation,” Soils and Foundations, vol. 56, no. 5, pp. 904–914, 2016.
View at: Publisher Site | Google Scholar
K. Kim, Multi-Dimensional Frost Heave Modeling with SP Porosity Growth Function, Gradworks, Atlanta, GA, USA, 2011.
China Coal Industry Press, Artificial Frozen Soil Physics Mechanics Performance Test: MT/T 593-2011, China Coal Industry Press, Beijing, China, 2011.
H. Cai, S. Li, Y. Liang, Z. Yao, and H. Cheng, “Model test and numerical simulation of frost heave during twin-tunnel construction using artificial ground-freezing technique,” Computers and Geotechnics, vol. 115, Article ID 103155, 2019.
View at: Publisher Site | Google Scholar
X. Xu, J. Wang, and L. Zhang, Frozen Soil Physics, Science Press, Beijing, China, 2001.
H. A. Цытович, Frozen Soil Mechanics, Science Press, Beijing, China, 1985.
X. Xu and Y. Deng, Experimental Study on Water Migration in Frozen Soil, Science Press, Beijing, China, 1991.
R.-J. Hou, T.-X. Li, Q. Fu et al., “Research on the distribution of soil water, heat, salt and their response mechanisms under freezing conditions,” Soil and Tillage Research, vol. 196, Article ID 104486, 2020.
View at: Publisher Site | Google Scholar
Anhui Institute of Architecture and Technology, Test Report on Physical and Mechanical Properties of Conventional Soil and Frozen Soil in Shanghai Metro Line 4 Restoration Project, Anhui University of Science and Technology, Huainan, China, 2006.
H. Karlsson, Sorenson, ABAQUS, Vélizy-Villacoublay, France, 2002.
J. Qiao, L. Tao, and S. Mi, “Simulation study on deformation characteristic of ground subway in horizontal freezing construction of metro tunnel,” Chinese Journal of Rock Mechanics and Engineering, vol. 23, no. 15, pp. 2643–2646, 2004.
View at: Google Scholar
Y. Cui, Brief Handbook of Well Construction Engineering, Coal Industry Press, Beijing, China, 2003.
F. T. She, S. J. Shao, J. Fang, and Y. Y. Jiao, “Calculation model for ground settlement process of Xi’an subway tunnel non-synchronous excavation with double-line,” Applied Mechanics and Materials, vol. 90–93, pp. 2221–2225, 2011.
View at: Publisher Site | Google Scholar
Y. Li, Study on Frost Heave Law of Stratum during Horizontal Freezing Construction of Metro Twin-Track Tunnel, Anhui University Science and Technology, Huainan, China, 2019.
H. Cai, L. Peng, and T. Zheng, “A duration prediction model of surface frost heave induced by tunnelling with horizontal freezing method,” Rock and Soil Mechanics, vol. 33, no. 6, pp. 1761–1768, 2012.
View at: Google Scholar
C. Wu, Probability Theory of Economic Mathematics and Mathematical Statistics, Higher Education Press, Beijing, China, 2015.

Copyright

Copyright © 2020 Rong-bao Hong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

542

Downloads

1026

Citations