Abstract

Groundwater seepage significantly affects the temperature field of a cold region tunnel. Laboratory model tests are carried out to evaluate its effects, yielding four main results. First, groundwater seepage can increase tunnel air temperature and decrease the thickness and length of the tunnel insulation layer. Second, groundwater seepage and tunnel ventilation exert a coupling effect on the surrounding rock temperature. This effect is related to the surrounding rock depth. Third, the influence of the groundwater seepage velocity on the temperature of the interface between the lining and surrounding rock demonstrates a spatial difference, and the groundwater seepage leads to an uneven temperature distribution at the interface between the lining and surrounding rock. Furthermore, under groundwater seepage, the shape and size of the tunnel cross section have significant effects on the interface temperature. Fourth, the cold region tunnel has an antifreezing capability that is mainly related to the frost heaving of the surrounding rock and the groundwater seepage velocity. This capability should be fully utilized in the design of cold region tunnels. The experimental data presented can be used to verify the reliability of the theoretical calculation model for tunnel temperatures in cold regions.

1. Introduction

Half of the world’s land area is covered by transient permafrost, seasonal frozen soil, or permafrost [1, 2]. Cold region tunnels have been increasingly built in recent years but freezing damage to such tunnels threatens their structural and operational safety [36]. Aiming to solve the problem of freezing damage to cold region tunnels, investigations on their temperature field have increased. For example, Lai et al. [7] used dimensionless parameters and perturbation methods to approximate the temperature distribution of the surrounding rock of cold region tunnels. Zhang et al. [8] proposed an analytical solution for heat conduction in tunnels, considering both composite medium and time-dependent boundary conditions. Lai et al. [9] derived the governing equations for the coupled problems of seepage, temperature, and stress fields, with consideration of the ice-water phase change. Zhang et al. [1012] derived an element calculation formula for the thermo-hydro coupling model to explore how different construction seasons, initial temperatures, and thermal insulation thicknesses affect the temperature field of cold region tunnels. Tan et al. [1315] used the “three-zone” model to derive the governing equation for the temperature field of surrounding rock. The airflow inside the tunnel and the heat convection between the air and lining were simulated using the k-s turbulence model and the wall function. Li et al. [16] established a water-thermal coupling mathematical model for permafrost to determine the optimal thickness of the tunnel insulation layer. Yan et al. [17] performed a numerical investigation on the influences of different ventilation velocities on the design parameters of the thermal insulation layer.

The above theoretical studies demonstrate the complexity of tunnel heat transfer in cold regions. Model tests can also be used to study the temperature field of the cold region tunnel. Feng et al. [18] developed a cold region tunnel model consisting of a refrigeration system, a circulation system, and a temperature control system with a 1 : 25 geometric reduced scale to study the temperature field and frost heave force of cold region tunnels and the reliability of thermal insulation layers. Zhang et al. [19] investigated the effects of the heat released during construction and of boundary temperatures on the temperature field of surrounding rock in permafrost tunnels using a model experiment with a 1 : 26.83 geometric reduced scale. Speeding up construction and installing insulation layers effectively weakened the negative effects of thermal disturbance during construction and of boundary temperatures. Zeng et al. [20] explored the influences of inlet airflow temperatures and mechanical ventilation on the temperature distribution of surrounding rock via a model test with a 1 : 30 geometric reduced scale under ventilation conditions. They found that mechanical ventilation in the positive ventilation direction adversely affects the freezing damage to cold region tunnels. Liu et al. [21] built a tunnel model with a 1 : 37 geometric reduced scale to explore the distribution of the temperature field and the expansion law of frost front under different inlet airflow temperatures and wind speeds.

Groundwater seepage has significant effects on the temperature field of the surrounding rock in cold region tunnels [9, 13, 2224]. The frozen depth of cold region tunnels is significantly affected by the seepage velocity [9]. However, the above model tests did not consider its effects on the temperature field of cold region tunnels. In this paper, heat transfer model tests of cold region tunnels experiencing groundwater seepage are carried out. The model test results are used to analyze the distribution characteristics of the temperature field and to reveal the gradual freezing process of cold region tunnels under groundwater seepage. The results are also used to analyze the causes of freezing damage to cold region tunnels and to thereby provide a basis for its prevention and control.

2. Model Test

2.1. Investigating the Freezing Damage to a Prototype Tunnel

The freezing damage (shown in Figure 1) to two highway tunnels, one in the Inner Mongolia Autonomous Region and the other in Hebei Province, is investigated.

The highway tunnel from the Inner Mongolia Autonomous Region is a separated tunnel that is 3,960 m long on the left and 3,915 m long on the right, and the width of each tunnel is 12 m. The average wind speed is 1.2 m/s, and the mean annual temperature is −2.6°C. The yearly frost-free period lasts for approximately 95 days. The frost penetration depth is 3 m. The tunnel was opened to traffic in November 2012. In January 2013, the average temperature in was below −25°C, and a large number of leakages and extensive road ice occurred in the tunnel.

The highway tunnel from Hebei Province is a separated tunnel that is 2,946 m long on the left and 2,992 m long on the right, and the width of each tunnel is 12 m. The average wind speed is 2.4 m/s, and the mean annual temperature is 3.2°C. The yearly frost-free period lasts for approximately 120 days. The frost penetration depths of the entrance and exit of the tunnel are 1.4 m and 2 m, respectively. In January 2016, the tunnel experienced extreme cold weather. The average temperature in January 2016 was close to −25°C. Freezing damage, such as leakage and water freezing, occurred in the tunnel.

Based on the freezing damage to two highway tunnels in the Inner Mongolia Autonomous Region and Hebei Province, a model test is conducted to investigate the distribution of the tunnel temperature field and the mechanism of freezing damage under various groundwater seepage velocities.

2.2. Similarity Criteria for the Prototype and Model Tunnels

The transient heat transfer in the concrete lining and surrounding rock of the tunnel prototype is governed by Fourier’s law in the polar coordinates [25], as shown in the following equation:where subscript p indicates the tunnel prototype, is the temperature, r is the distance from the center of the tunnel, t is the time, is the thermal diffusivity, and is the inner radius of the tunnel lining.

The transient heat transfer in the concrete lining and surrounding rock of the tunnel model is governed by Fourier’s law in the polar coordinates, as shown in the following equation:where subscript m indicates the tunnel model.

According to similarity theory [18, 21], the heat transfer similarity criteria can be defined as follows:where and are the heat transfer similarity criteria, is the latent heat of the phase change due to groundwater freezing, and is determined as follows:

Based on equation (3), the following similarity relationships of temperature and time can be obtained:

The groundwater seepage in the rock surrounding the tunnel can be modeled as a homogeneous flow in a medium with an effective porosity [26]. The groundwater seepage velocities of the tunnel prototype and tunnel model can be obtained using Darcy’s law, as shown in the following equation, respectively:where is the groundwater seepage velocity, is the water head, and is the hydraulic conductivity.

The water head of the tunnel model is the same as that of the tunnel prototype. Based on equation (6), the following similarity relationship for groundwater seepage velocity can be obtained:

Considering the research results of previous tunnel heat transfer model tests [16, 1821] and previous test conditions [22, 23], a geometric similarity ratio of 1 : 30 between the tunnel model and tunnel prototype is selected. According to equations (5) and (7), similarity ratios of 1 : 1, 1 : 900, and 30 : 1 for temperature, time, and velocity, respectively, are obtained.

2.3. Model Test Design

Figure 2 describes the design of the cold region tunnel model test under the conditions of tunnel ventilation and groundwater seepage. In the model test, the fractured rock mass is considered to be a homogeneous medium with an effective porosity [26], the groundwater seepage in the tunnel surrounding rock is modeled as a homogeneous flow in a homogeneous medium with an effective porosity, and the sand soil is used to simulate the tunnel surrounding rock. Constant water head tests are used to simulate the groundwater seepage field of the tunnel surrounding rock, and water with a constant inlet temperature is used to simulate the underground temperature field. All of the tests have the same inlet temperature which is 13°C. To ensure a uniform groundwater seepage field, groundwater seepage occurs from bottom to top. The groundwater seepage velocity is controlled by adjusting the water head difference. Cold airflow is blown into the tunnel through the closed circulation ventilation ducts to freeze the groundwater in the tunnel surrounding rock, and the wind speed variation is controlled by changing the voltage applied to the air blower. The temperatures of the surrounding rock, tunnel lining, and air in the tunnel are measured by temperature sensors. The wind speed is measured by using an ultrasonic anemometer, and the groundwater seepage velocity is measured by using a liquid mass flow meter.

2.4. Model Test Apparatus and Materials

Figure 3 presents the apparatus of the laboratory model test, which consists of the tunnel model box, the tunnel ventilation system, the groundwater seepage system, and the measuring system.

Table 1 summarizes the model test equipment and devices.

2.4.1. Tunnel Model

The tunnel model consists of a 1.4 m × 1.2 m × 1.2 m stainless steel box and a 1.8 m concrete tube with an outer diameter of 40 cm and a thickness of 2 cm. The concrete tube, which has a circular cross section, is made of cement, sand, gravel, and water in a ratio of 1 : 2.4 : 3.6 : 0.65. Sand is used in the test. The sand-raining method is used to fill the tunnel model with sand from the Yangtze River, the density and thermal conductivities of which are 2.05 g/cm3 and 1.58 W/(m·°C), respectively. To minimize the heat transfer boundary effect, a 5-cm thick thermal insulation layer is installed around the outer surface of the stainless steel box. The plastic permeable boards are installed at the top and bottom of the stainless steel box to ensure that the groundwater seepage spreads over the surrounding rock. The permeable plastic boards are wrapped with geotextile to prevent the groundwater seepage from carrying away the sand.

2.4.2. Tunnel Ventilation System

The ventilation system consists of a refrigerator, an air blower, a thermal insulation box, and a thermal insulation pipe. The thermal insulation pipe connects the air blower, the thermal insulation box, and the concrete tube to form closed circulation ventilation ducts. The air blower is installed on the top of the thermal insulation box, and the refrigerator is installed in the thermal insulation box. The air inside the refrigerator can be cooled to −30°C, and the airflow velocity in the ducts can reach 4 m/s.

2.4.3. Groundwater Seepage System

The groundwater seepage system consists of a constant temperature water tank, a pump, two thermal insulation water tanks, two lifting platforms, a water supply pipe, and a return pipe. The water supply pipe and return pipe connect the constant temperature water tank, the pump, the thermal insulation water tanks, and the stainless steel box to form a closed circulation system. Both thermal insulation water tanks are placed on the lifting platforms. The groundwater seepage velocity is controlled by adjusting the height difference between the two thermal insulation water tanks.

2.4.4. Measuring System

To monitor the temperature variations in the air, surrounding rock, and the interface between the tunnel lining and surrounding rock, a temperature monitoring cross section is installed in the middle of the model box. The platinum resistance temperature sensor TA is installed in the center of the concrete tube to monitor the temperature of the airflow. The platinum resistance temperature sensors TS7 to TS11 are installed in the sand to monitor the temperature variation of the surrounding rock temperature field at different positions and depths. The platinum resistance temperature sensors TS1 to TS6 are installed on the concrete tube surface, equally spaced along the cross section, to monitor the temperature variation of the interface between the tunnel lining and surrounding rock. The wind speed is monitored using an ultrasonic anemometer that is placed at the tunnel exit. Finally, the liquid mass flow meter is connected to the water return pipe of the groundwater seepage system.

2.5. Model Test Strategy

Two factors are considered in the model test, namely, wind speed and groundwater seepage velocity. Zhang et al. [22] presented a systematic study of the heat transfer in tunnels under groundwater seepage. Groundwater seepage velocities of less than 0.864 m/d had negligible effects on the temperature field of the surrounding rock. Thus, groundwater seepage velocities of 0, 0.923 and 2.743 m/d are chosen for this model. The average wind speed of the tunnel in the Inner Mongolia Autonomous Region and in Hebei Province is 2.2 m/s and 2.4 m/s, respectively. Thus, wind speeds of 0.606, 1.002, and 2.398 m/s are chosen for this model. Table 2 presents the model test strategy. The experiment time of each test is different, as groundwater seepage can prevent the surrounding rock from freezing. Hence, the time that the temperature of surrounding rock down to 0°C varies with groundwater seepage velocities. The experiment times of 11.5, 23, and 32.5 h for groundwater seepage velocities of 0, 0.923, and 2.743 m/d, respectively, are chosen for this model. These experiment times can guarantee a part of the surrounding rock around the lining will freeze.

3. Results and Discussion

3.1. Effect of Groundwater Seepage on the Tunnel Air Temperature

Figure 4 presents the variations in tunnel air temperature over time under the effects of groundwater seepage and ventilation in the model test. Specifically, it shows that the air temperature in the tunnel decreases. Furthermore, the air temperature reduction rate decreases over time, and the temperature curve fluctuates every 12 h due to the deicing of the refrigeration equipment. Figure 4 also shows that groundwater seepage velocity and wind speed have significant influences on tunnel air temperature. Figure 5 presents the variations in tunnel air temperature for various groundwater seepage velocities at different wind speeds after cooling for 11.5 h. In Figure 5, the air temperature in the tunnel increases almost linearly as the groundwater seepage velocity increases. Furthermore, the effect of groundwater seepage velocity on air temperature in the tunnel is related to wind speed. The lowest air temperature in the tunnel is −22.55°C under a groundwater seepage velocity of 0 m/s and a wind speed of 0.606 m/s. The highest air temperature in the tunnel is −14.58°C under a groundwater seepage velocity of 2.743 m/d and a wind speed of 2.398 m/s. Thus, the higher the groundwater seepage velocity and wind speed are, the higher the air temperature is in the tunnel.

The above analysis shows that groundwater seepage has a significant effect on the air temperature in a tunnel. The higher the groundwater seepage velocity is, the higher the air temperature is in the tunnel. Low air temperature in a tunnel may induce tunnel freezing damage. The thickness and length of the tunnel insulation layer are determined by the air temperature in the tunnel. Specifically, the lower the air temperature is in the tunnel, the thicker and longer the tunnel insulation layer needs to be to avoid freezing damage. It can be speculated that groundwater seepage has a significant influence on the thickness and length of the tunnel insulation layer in cold regions, causing it to thin and short. Thus, the influence of groundwater seepage should be considered when designing the insulation layers of cold region tunnels.

3.2. Effect of Groundwater Seepage on the Surrounding Rock Temperature

Figure 6 presents the variations in the surrounding rock temperature over time under the influence of groundwater seepage and ventilation in the model test. Specifically, it shows that the surrounding rock temperature at TS9, TS10, and TS11 decreases. The cooling rates and amplitudes of the surrounding rock temperature differ significantly between TS9, TS10, and TS11. The cooling rate and amplitude of the surrounding rock temperature are related to the depth of the tunnel surrounding rock, the groundwater seepage velocity, and the wind speed. The cooling rate and amplitude of the surrounding rock temperature at the TS9 monitoring point are the largest and those at the TS11 monitoring point are the smallest. Thus, the deeper the tunnel surrounding rock is, the lower the cooling rate and amplitude of the tunnel surrounding rock temperature are. Lower groundwater seepage velocities and higher wind speeds lead to greater cooling rates and amplitudes of the surrounding rock temperature. The influence of groundwater seepage and tunnel ventilation on the surrounding rock temperature field is related to the depth of the surrounding rock. The increase in wind speed has a strong influence on the temperature field at TS9 and TS10 but a weak influence on that at TS11. Figure 6 also shows that the wind speed has the most significant influence on the surrounding rock temperature field when the groundwater seepage velocity is 0.923 m/d but has the least significant influence on the surrounding rock temperature field when the groundwater seepage velocity is 0 m/d. This can be explained by the coupling effect of groundwater seepage and tunnel ventilation.

Figure 7 presents the variations in the surrounding rock temperature drops under groundwater seepage velocities at different wind speeds after 11.5 h of cooling. It also shows that the temperature drops in the tunnel surrounding rock decrease as the groundwater seepage velocity increases. The temperature drops differ at TS9, TS10, and TS11 under different wind speeds. The temperature drop at TS11 at different wind speeds decreases linearly as the groundwater seepage velocity increases. At wind speeds of 0.606 m/s and 1.002 m/s, the temperature drops at TS9 and TS10 decrease linearly as the groundwater seepage velocity increases. At a wind speed of 2.398 m/s and groundwater seepage velocity under 0.923 m/d, the temperature drops at TS9 and TS10 slowly decrease as the groundwater seepage velocity increases. At a wind speed of 2.398 m/s and groundwater seepage velocity over 0.923 m/d, the temperature drops at TS9 and TS10 sharply decrease as the groundwater seepage velocity increases.

The above analysis shows that groundwater seepage and tunnel ventilation have a coupling effect on the surrounding rock temperature. Specifically, higher groundwater seepage velocities and lower wind speeds lead to higher surrounding rock temperatures. Thus, the coupling influence of groundwater seepage and tunnel ventilation should be considered when designing the insulation layers of cold region tunnels.

3.3. Effect of Groundwater Seepage on the Temperature at the Interface between the Lining and Surrounding Rock

Figure 8 presents the temperature variation over time at the interface between the lining and surrounding rock under the influence of groundwater seepage and ventilation in the model test. It shows that the interface temperatures of TS1, TS3, TS4, TS5, and TS6 decrease over time when groundwater seepage and ventilation are considered in the model test. Groundwater seepage has a significant effect on the interface temperature. The effect is related to the location of the monitoring points. Therefore, the cooling rates and amplitudes of the interface temperature for TS1, TS3, TS4, TS5, and TS6 differ significantly. The groundwater seepage has the most significant effects on measuring points TS3 and TS4 and has the least significant effects on measuring points TS1 and TS6. Wind speed also has a significant effect on the interface temperature. The effect is related to the groundwater seepage velocity. These results indicate that groundwater seepage and tunnel ventilation have a coupling effect on the interface temperature and that this coupling effect is related to the interface location.

Figure 9 shows the variations in the temperature drops at TS1, TS3, TS4, TS5, and TS6 with groundwater seepage velocity under tunnel ventilation after 11.5 h of cooling. It shows that the temperature drops at TS3 and TS4 upstream of the groundwater seepage field decrease as the groundwater seepage velocity increases. However, the temperature drops at TS1 and TS6 downstream of the groundwater seepage field increase as the groundwater seepage velocity increases from 0 to 0.923 m/d. Furthermore, the temperature drops at TS1 and TS6 decrease as the groundwater seepage velocity increases from 0.923 to 2.743 m/d. The smallest temperature drop occurs at TS3 (upstream of the groundwater seepage field), and the largest temperature drop occurs at TS1 (downstream of the groundwater seepage field) under groundwater seepage. These results indicate that groundwater seepage leads to an uneven temperature distribution at the interface between the lining and surrounding rock and that tunnel ventilation enhances this uneven distribution.

Figure 10 presents the temperature differences between TS3, TS5, and TS6. The temperature at TS3 is higher under groundwater seepage than those at TS5 and TS6. The interface temperature upstream of the groundwater field is higher than that downstream. Figure 10 also shows that the temperature difference between TS3 and TS6 increases as the groundwater seepage velocity increases. Thus, the higher the groundwater seepage velocity is, the higher the temperature difference is between TS3 and TS6. The temperature difference between TS3 and TS5 increases as the groundwater seepage velocity increases from 0 to 0.923 m/d and decreases as the groundwater seepage velocity increases from 0.923 to 2.743 m/d. Groundwater seepage has a great effect on the interface temperature. This effect has a spatial difference. Thus, both the shape and size of the cross section of the tunnel significantly affect the interface temperature under groundwater seepage.

3.4. Effect of Groundwater Seepage on the Freezing Damage to Cold Region Tunnels

The tunnel drainage system is located between the lining and the surrounding rock. Monitoring points TS3, TS5, and TS6 are located at the top, middle, and foot of the tunnel drainage system, respectively. The temperatures at TS3, TS5, and TS6 determine the antifreezing ability of cold region tunnels. Therefore, the higher the temperatures at TS3, TS5, and TS6 are, the stronger the antifreezing capability of the tunnel in cold regions is . Figure 8 shows that when the groundwater seepage velocity is 0 m/d, the measuring points TS3, TS5, and TS6 freeze synchronously. If the tunnel surrounding rock is free of frost heaving, cold region tunnels may not need a thermal insulation layer. When the groundwater seepage velocity is 0.923 m/d, the monitoring point TS6 freezes first. The interface then gradually freezes from TS6 to TS5 along the circumferential direction of the tunnel section. This explains why cold region tunnels suffer tunnel lining leakage and road icing. The results of the model test reveal the cause of tunnel lining leakage and road icing. To avoid tunnel leakage in cold region tunnels, a thermal insulation layer should be properly designed and a heating device should be installed in its longitudinal drainage pipe when the groundwater seepage velocity is low. When the groundwater seepage velocity is 2.743 m/d, the monitoring point at TS6 freezes but the monitoring points at TS3 and TS5 do not. Furthermore, their temperatures do not decrease significantly as cooling time increases. Therefore, the tunnel has good antifreezing capability when the groundwater seepage velocity is high. This good antifreezing capability should be fully utilized. Furthermore, the foot of the tunnel drainage system should be strengthened against freezing via an insulation layer and heating device.

The above analysis demonstrates that cold region tunnels have an antifreezing capability that is mainly related to the frost heaving of the surrounding rock and the groundwater seepage velocity. This antifreezing capability should be fully utilized in the antifreezing design of cold region tunnels.

4. Conclusions

A systematic study on the temperature field of cold region tunnels is presented in which a model test that considers the coupling effects of groundwater seepage and tunnel ventilation is conducted. The effects of groundwater seepage on the air temperature in the tunnel, the surrounding rock temperature, the temperature at the interface between the lining and surrounding rock, and the freezing damage of the tunnel in cold regions are studied. The following four main conclusions are drawn.

First, groundwater seepage has a significant effect on the air temperature in the tunnel. The higher the groundwater seepage velocity, the higher the air temperature in the tunnel. Groundwater seepage has a significant influence on the thickness and length of the tunnel insulation layer in cold regions. The thickness and length of the tunnel insulation layer can decrease under groundwater seepage.

Second, groundwater seepage and tunnel ventilation have a coupling effect on the surrounding rock temperature. Higher groundwater seepage velocity and lower wind speed lead to higher surrounding rock temperature. The influence of groundwater seepage and tunnel ventilation on the surrounding rock temperature field is related to the depth of the surrounding rock.

Third, groundwater seepage leads to an uneven temperature distribution at the interface between the lining and surrounding rock. The tunnel ventilation further enhances this uneven temperature distribution. The groundwater seepage and tunnel ventilation have a coupling effect on the interface temperature. This coupling effect is related to the interface location. The influence of the groundwater seepage velocity on the interface temperature has a spatial difference. The shape and size of the cross section of the tunnel have a significant influence on the interface temperature under groundwater seepage.

Fourth, cold region tunnels have an antifreezing capability that is mainly related to the frost heaving of the surrounding rock and the groundwater seepage velocity. The antifreezing capability of the tunnel should be fully utilized in the antifreezing design of cold region tunnels.

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.

Acknowledgments

This study was supported by the Key Laboratory of Geotechnical and Underground Engineering of Ministry of Education, Tongji University (No. KLE-TJGE-B1605), the Key Laboratory of Rock Mechanics and Geohazards of Zhejiang Province (No. ZJRM-2018-Y-01), and the Sichuan Science and Technology Program (No. 2018GZ0359).