Abstract

To fully understand the temperature distribution of cold regions and the variation law of temperature fields in cold-region tunnels, this paper presents a case-history study on a tunnel located on the eastern Qinghai-Tibet Plateau, China. The conclusion is as follows: the temperature outside the tunnel and the ambient temperature are affected by wind speed and light. The law of the temperature field in the tunnel is greatly affected by wind speed and wind direction. According to the field test, the wind speed in the tunnel is about 2.8 m/s in winter, and the daily average temperature at the exit of the tunnel is basically lower than that at the entrance. From the central to the entrance, the temperature in the tunnel decreases by 0.11°C every 10 meters along the longitudinal direction; from the central to the exit, the temperature in the tunnel increases by 0.07°C every 10 meters. In this regard, for the problems of lining frost damage and central drainage pipe freezing, it is suggested to adopt the way of heating and drainage, but heating the freezing area outside the drainage pipe should be avoided. The test results can provide references for the design, construction, and research of the temperature field of the tunnel antifreezing system in the cold region. It is hoped that the test results can be useful in the design and construction of frost damage prevention systems and the investigation of temperature fields in cold-region tunnels.

1. Introduction

Permafrost regions account for approximately 23% of Earth’s surface, mainly in Russia, Canada, China, and the US state of Alaska [13]. Of these regions, the permafrost regions in China, which account for approximately 10% of Earth’s surface, represent the world’s third-largest permafrost area. China’s permafrost regions are primarily distributed in the Tibetan Plateau, which is known as the world’s “third pole.” Additionally, seasonal frozen-ground regions in China, which account for approximately 53.5% of the nation’s total land area, are distributed in the centre-north. Tunnelling in cold regions changes the thermal balance of rock masses. However, the temperature range in which cold-region engineering is practiced is significantly influenced by the thermal conductivity of the soil-rock media, which is related to engineering problems caused by frost heave [46]. Generally, when frost heave occurs because of low temperatures, the properties of the rock mass change, and frost deformation significantly affects tunnel stability. Frost damage, such as peeling, lining surface fractures, seepage, icing, slippery collapse, and tunnel freeing port freezing, is primarily caused by freeze-thaw cycles. The substantial difference in air temperature between morning and evening significantly weakens service capacity and life.

In Northeast China, the Baikal tunnel, Xiluo tunnel, and Cuiling tunnel. The Yuximolgai tunnel built in Xinjiang have ice jam due to the impact of low temperature and frost damage. The thickest ice accumulated in the tunnel is more than 1.6 meters, which results in the impassable of the tunnel every winter [7]. Cold-region tunnels that suffer cold damage are widely distributed in China (Figure 1). To ensure the normal operations of such tunnels, a substantial quantity of manpower, material, and financial resources are required. Cold damage results in substantial monetary losses for the government [8, 9]. Thus, research on cold-region tunnels is highly important.

In a subzero temperature environment, tunnels generate a few problems. Although numerous studies on cold-region tunnels have been conducted in recent decades [1012], the frost-heaving force in rock masses and the ground and air temperature fields have been insufficiently analysed. To understand the temperature distribution of cold regions and the variation law of temperature fields in cold-region tunnels, the meteorological conditions in cold regions and the ground temperature in the tunnels must be measured systematically [1315]. The following results for different aspects of cold-region tunnels have been obtained. For instance, in the US, the Cold Regions Research and Engineering Laboratory (CRREL) [16, 17] has obtained a series of valuable research results based on a long-term (since 1963) temperature observation of permafrost tunnels. Using 1974 meteorological data for New Delhi, Bansal et al. [18] studied the annual variation of the ground temperature distribution for different surface conditions in warm regions. Using a test tunnel in the cold region of the Heilongjiang Province, China, Li [19] determined a temperature distribution curve for in- and outside the tunnel. He et al. [20] studied the freezing and thawing of the rock that surrounds the tunnels in cold regions, employing a combined convection-conduction model of turbulent airflow in tunnels. Lai et al. [21] performed nonlinear analyses for the coupled problem of temperature and seepage fields in cold-region tunnels using Galerkin’s method. Through numerical investigation, field testing, and local data observation, Goreing and Kumar [22] investigated the protective effect of a riprap embankment on permafrost. According to the measured temperature fields in a cold-region tunnel, Yan et al. [23] analysed the law that regulates the annual and daily temperature in- and outside the tunnel. Considering the temperature variations according to coordinates, Qin et al. [24, 25] and Prashantk [26, 27] employed the superposition principle and the variable separation approach to obtain analytical solutions for a transient temperature field in a multilayer cylinder under the convective boundary condition. Based on the extensive testing of inner and outer tunnel temperature fields, Chen and Zan [28, 29] determined the temperature change regularity and the change characteristics of the maximum frozen soil depth.

To the best of our knowledge, the construction of tunnels in cold regions is a difficult task. It is of positive significance for the design of the tunnel structure to master the variation law of the ground temperature and air temperature field. However, the conclusion of the law of ground temperature and air temperature field cannot be based on a single tunnel in the cold region. We need to monitor and summarize the temperature field of many tunnels in cold regions. To understand the freezing and thawing of the rock mass and the distribution of the air temperature in- and outside such tunnels and to provide reference data for practical engineering designs, this paper presents a case study on a cold-region tunnel. Based on the long-term field testing at the Qingshashan tunnel, China, 14 test sections were selected. In addition, a regressive analysis of the temperature field data was conducted. The analysis included the inside and surface temperature fields of the liner, the rock mass temperature field, and the air temperature field. Specifically, the relationship between the temperature in the tunnel and the outer atmospheric temperature is discussed. Additionally, the temperature’s longitudinal variation curve without mechanical ventilation in the tunnel is presented.

2. Project Overview

The Qingshashan tunnel on the Tibetan Plateau, China (Figure 2), was constructed as part of a highway improvement program from Pin’an to Adai. The tunnel is a separated tunnel with traffic lanes that run in opposite directions. The Qingshashan tunnel is the longest highway tunnel in Qinghai Province. The tunnel’s two traffic pathways are 3350 m and 3340 m long, respectively. The Qingshashan tunnel section looks like horseshoe, and the clear width and height of tunnel construction clearance are 9.25 m and 5.0 m, respectively. The tunnel is in the Qingshashan area belonging to the denudation of mountain topography, where the elevation ranges from 3004.59 m to 3398.81 m, and the highest height difference is approximately 354 m. The climate is that of a typical continental semiarid plateau. The Qingshashan tunnel is in an area in which the seasonal frost depths are between moderate and deep and receives moderate frost damage [30]. The ground of the Qingshashan area is a typical seasonal frozen ground with alpine meadows as the dominant vegetation type and with complicated topography, steep slope, water erosion, U-shaped gully development, and a substantial temperature difference between morning and evening. The winters are bitterly cold and long, whereas the summers are cool. Based on long-term historical meteorological data, the mean annual air temperature is 5.90°C, the extreme high temperature is 38.70°C, the extreme low temperature is −26.60°C, and the maximum frost depth is approximately 1.48 m. Rainfall in this area increases with the gradually rising terrain, with a mean annual rainfall of approximately 500 mm and a mean annual evaporation of 1,636 mm, which primarily occurs from June to September. Southeast winds are typical, and the average wind speed is 1.8∼3.9 m/s. Additionally, there are large numbers of faults in the tunnel location, and the fissure phreatic water is primarily distributed in the fractured fault zone, whose water yield property is medium. The tunnel excavation area is rich in water, and there is a risk of frost damage in the tunnel.

3. Methods and Instrumentation

3.1. Ground Temperature Monitoring

After a tunnel is opened through a mountain, the convection of air in the tunnel changes the initial thermal balance of the rock that surrounds the tunnel. That is, freezing and thawing occurs in the tunnel rock mass because of the change of air temperature in the tunnel. Because the ground temperature field can provide a beneficial reference for the design of tunnel antifreezing insulation, we installed submerged platinum metal thermistors (JMT-36X) (Figure 3) at equal radial intervals along the tunnel to monitor the ground temperature. As shown in Figure 4, the JMZX200X comprehensive tester was used with the temperature sensors. The JMT-36X temperature sensor is a new type of sensor with high accuracy, high stability, high reliability, moisture-proof, and good insulation. It is often used in large-scale projects such as railway and highway. Its sensitivity is 0.1°C, accuracy is within ±0.5°C, and its range is −30°C ∼ +120°C. Thermometric elements, which were longitudinally installed along the tunnel in symmetrical distribution, were buried in the tunnel’s rock masses to measure the temperature of the rock. Their locations are shown in Figure 5. The distance between the thermometric element (number) and the liner surface is detailed in Table 1. Additionally, considering that the temperature data of the left tube were seriously affected by tunnelling, only the temperature data of the Qingshashan tunnel’s right tube were analysed. The temperature was monitored every six hours: 2:00 am, 8:00 am, 2:00 pm, and 8:00 pm.

3.2. Air Temperature Monitoring

The air temperature in- and outside the tunnel was recorded automatically using a temperature and humidity automeasurement system (ZJI-2B) (Figure 6) and was corrected strictly in accordance with the requirements. The accuracy of temperature range of the ZJI-2B weekly temperature hygrometer is ±1°C, and its accuracy of humidity range is ±5% RH. To measure the significant temperature change that occurs outside the tunnel, thermohygrometers were installed at the tunnel portals and at sites in the tunnel 500 m from the entrance. Additionally, a thermohygrometer was installed in the middle section of the tunnel, and outdoor thermohygrometers were installed outside the tunnel according to the meteorological requirements, i.e., at sites outside the tunnel and 50 m from the entrance and exit. Moreover, the outdoor thermohygrometers were placed in shelters that were 1.5 m high and distant from obstacles to obtain free-moving air standard values, whereas the thermohygrometers in the tunnel were suspended from the liner surface at the height of 1.5 m to reduce sensitivity to temperature and interference from dust. According to international meteorological test practices, the temperature was monitored every six hours: 2:00 am, 8:00 am, 2:00 pm, and 8:00 pm [31]. The longitudinal distribution of the thermohygrometers along the tunnel is shown in Figure 7.

4. Results

4.1. Discussion of the Ground Temperature Field
4.1.1. Ground Temperature of All Testing Sections

The ground temperature distribution at the different testing sections over time is shown in Figures 8(a)8(f). It can be observed that the ground temperature distribution, which consists of reading from sites 3 m, 8 m, and 23 m from the entrance and sites 5 m and 15 m from the exit, can be divided into three stages based on time. The first stage is from November of the previous year to early April of the second year. The deeper the surrounding rock is buried, the higher the measured temperature value is. In this period, the temperature after the tunnel lining and within the range of 0.5 m surrounding rock has been negative. In the coldest February, the negative temperature also occurs for a long time in the depth of 1 m surrounding rock. The second section is from the beginning of April to the middle of August, during which the temperature in the surrounding rock can be regarded as from negative to positive and then reaches the highest temperature in the whole year. The third stage is from the middle of August to November. The rate of temperature decrease is faster than that of temperature rise, so the curve of temperature decrease is steeper than that of temperature rise.

Similarly, the ground temperature distribution along the tunnel (the sites in the tunnel that are 876 m and 1000 m from the entrance and the site that is 800 m from the exit) can also be divided into three stages based on time.

The first stage is from November of the first year to the end of March of the next year. In this stage, the temperature in the tunnel is the lowest, but the temperature in the tunnel is quite different from that measured in the surrounding rock. The second stage is from the beginning of April to the middle of July. During this period, the surrounding rock temperature rises sharply, the temperature value is unstable, and the beating is frequent. The temperature curve of each measuring point in the surrounding rock is readjusted according to the position. In winter, the curve of the low temperature rises to the highest temperature, and the curve of the highest temperature reaches the lowest temperature. The third stage is from the middle of July to November. The temperature drop curve in the surrounding rock at the entrance of the tunnel is synchronous with the change of the atmospheric temperature in the tunnel area, while the temperature in the surrounding rock at the exit of the tunnel lags behind the temperature in the tunnel.

The lowest temperature of the testing sections and the corresponding time is listed in Tables 2 and 3. The lowest temperature of section K31 + 668 is lower than that of section K33 + 218, which is contrary to what is generally assumed. The primary reason for this lower temperature is that the natural wind blows from the exit to the entrance during winter. In addition, the distance between section K31 + 668 and the entrance is longer than that between section K33 + 218 and the exit, which can easily lead one to assume that the temperature of section K33 + 218 should be significantly lower than that of section K31 + 668. From our analysis, the reason for this phenomenon is that the water storage in the rock mass near section K33 + 218, which is likely groundwater, is abundant.

4.1.2. Longitudinal Distribution of Ground Temperature

The line length of the thermometric elements with burial depths of 0.5 m, 1.0 m, 1.5 m, 2.0 m, and 2.5 m is 1.0 m, 1.5 m, 2.0 m, 2.5 m, and 3.0 m, respectively. In addition, the test results of the thermometric elements with burial depths of 1.5 m, 2 m, 2.5 m, and 3 m are shown in Figures 9(a)9(d). A comparison of Figure 9(a) with Figure 9(b) reveals that the ground temperature distribution on the interior of the liner always appears in an inverted U shape from November of the previous year to May of the following year. Additionally, overall, the temperature of the rock masses gradually increases from the tunnel portal to the middle of the tunnel. However, the ground temperature of the section that is 10 m from the exit is lower than that at the tunnel exit. Moreover, the temperature increase curves present an overall trend of a small increase at the entrance and a large increase at the exit. Furthermore, the ground temperature at the entrance is slightly higher than at the exit when the air temperature is very low, and the ground temperature at the entrance is lower than at the exit when the air temperature is very high. The main reasons for this phenomenon are that the tunnel exit only has one air inlet and that the distance between section K31 + 536 and the entrance is longer than that between section K33 + 218 and the entrance.

An exception to this pattern appears in the temperature curve of section K31 + 536 with the line length that ranges from 2 m to 3 m: the ground temperature of section K31 + 536 is lower than that of the same depth at the tunnel entrance and at a shallow depth in the same section when the air temperature is very low. However, the ground temperature of this section rapidly rises with the increasing air temperature after April. From our analysis, this phenomenon occurs because there is a substantial quantity of underground water at this area around the depth of 2.0∼3.0 m which may have infiltrated from the land surface. In addition, based on design drawings, on-site geological data, and construction records, the thermometric elements with the line length that ranges from 2 m to 3 m pass precisely through a fault. The water in the fault originates in a surface water source. The thermometric elements with the line length that ranges from 1 m to 1.5 m are installed outside the fault.

From the longitudinal ground temperature curves, we can see that the curves are relatively dense in all of the testing sections near the tunnel entrance, except section K31 + 536, because the temperature centralization of the measuring lines is acceptable and the monthly measured temperature difference is low. However, at the tunnel exit, the temperature curves tend to become more intensive with the increase in the burial depth of the thermometric elements and the increase in the distance between testing sections and the tunnel exit, which is consistent with theory and practice. In addition, there are obvious differences between the exit and the entrance, and the primary reasons for which are as follows: (1) the natural wind continually blows from the exit to the entrance. Thus, the temperature curves tend to be more intensive from the exit to the entrance. (2) The entrance of the tunnel’s right tube, which is used as a construction passage for the left tube, has a tremendous effect on the temperature in the tunnel. Therefore, the wind direction in the tunnel is the key parameter that affects the ground temperature distribution.

4.1.3. Radial Distribution of Ground Temperature

The radial temperature curves of the rock masses in February and July are shown in Figures 10(a)10(d) since those months are the coldest and the hottest ones among the entire year, respectively, It is important to note that the temperature curves in the other months are omitted for brevity which have similar change trends and ranges. The test results indicate that with the increase in the depth of the rock mass, the temperature gradually increases, and that the temperature increment in the shallow of the rock mass is slightly larger than in the deep site.

As shown in Table 4, the maximum negative air temperature and the maximum negative ground temperature are not in the same section. Although the entrance is located on the right side of the mountain, its maximum depth of the negative temperature is less than that at the exit. The primary reason is that the exit is the tunnel’s air inlet, and the temperature distribution is influenced by the tunnelling.

When the air temperature outside the tunnel is high, the thermometric elements close to the tunnel portal are nearly at the same temperature only with a small fluctuation. The temperature in the tunnel’s middle section also remains consistent. However, the temperature fluctuation of the elements at the entrance and the exit is large, which indicates that the air temperature at the entrance and the exit is unstable. Section K33 + 218 is located at the tunnel exit and experiences minor construction interference. The two measured temperature curves are highly regular. Thus, they can be used to predict the average temperature of the rock masses.

4.2. Discussion of the Air Temperature Field
4.2.1. Air Temperature outside the Tunnel

As shown in Figure 11, the air temperature near the tunnel entrance is lower than that near the tunnel exit throughout the year. The average temperature difference between the entrance and the exit reaches its maximum value when the air temperature decreases to its lowest, and the maximum difference is approximately 4°C. According to the annual wind direction monitored in the tunnel, the natural wind blows from the tunnel exit to the entrance during winter, and the entrance average wind speed is approximately 2.8 m/s. In summer, the wind direction in the tunnel is uncertain with a relatively low average wind speed. The tunnel entrance and exit are located on the right side and the southern slope of the mountain, respectively, which is most likely a key parameter that influences the air motion in the tunnel. Throughout the year, the daily average temperature, which is generally below 0°C, is sustained for 180 days at the exit and 165 days at the entrance.

4.2.2. Air Temperature inside the Tunnel

As shown in Figures 12(a) and 12(b), the air temperature difference between the entrance and exit is large during the coldest month. On the one hand, the air temperature at the entrance is influenced by the construction interference. The entrance of the tunnel’s right tube, which is the construction passage of the left tube, substantially affects the temperature in the tunnel. On the other hand, the natural wind continually blows from the exit to the entrance, which is the key parameter that influences the air temperature in the tunnel. At the tunnel exit, because of the stable wind direction in winter, the low-temperature air outside the tunnel fully transfers the heat energy to the secondary liner and rock masses. Therefore, the temperature of the liner and the rock mass at a certain depth is generally below 0°C for a long period. Throughout the year, the temperature, which is generally below 0°C, is sustained for 210 days in the section 150 m from the exit. Thus, the temperature fields of the liner and the rock masses at the exit of the tunnel’s right tube can be adopted as the evaluation basis for the freeze proofing and heat insulation in this cold-region tunnel. Because of ventilation after the holing-through of the tunnel, heat exchange occurs between the air in the tunnel and the rock mass. With a faster wind speed, more cold air enters the tunnel, and the air temperature in the tunnel decreases accordingly.

4.2.3. Longitudinal Air Temperature Distribution

As shown in Figure 13, the air temperature distribution in the longitudinal direction presents an inverted V-shaped pattern from January to March. That is, the highest temperature occurs in the middle section of the tunnel. In addition, the monthly average temperature difference between the tunnel portal and the middle section can be as much as 11°C. Moreover, the air temperature distribution in the longitudinal direction presents an approximately W-shaped pattern from April to May, and the temperature in the middle section of the tunnel remains high. During this period, the temperature at the exit is rapidly increasing and the temperature difference between the exit and the outdoors is small. Meanwhile, with the minimum temperature occurring in the section that is 150 m from the entrance, the air temperature distribution at the entrance changes from one way increasing into a V-shaped pattern.

In addition, the air temperature at the entrance increases sharply from January to March, which indicates that the air temperature in the tunnel is significantly higher than that outside the tunnel during the coldest month. The exit air temperature behaves similarly. The primary reason for this phenomenon is that (considering the smaller cross section and the greater length) the heat exchange between the air in the tunnel and the atmosphere outside the tunnel is not significant enough. When the monthly average air temperature is slightly above 0°C, the air temperature in the tunnel’s first 500 m is lower than at the entrance. Between the section that is 500 m from the tunnel portal and the middle section of the tunnel, a 0.11°C increase in the air temperature per 10 m is observed in a longitudinal direction from the entrance. Similarly, a 0.07°C increase in the air temperature per 10 m is observed in a longitudinal direction from the exit.

According to previous studies, the climate temperature in one year can be considered to change with time as a sine function (or a cosine function). The extreme temperature of the climate around the Qingshashan tunnel differs according to altitude and location, and the design target time is November 20. Based on the geological investigation and the observational data collected in the field experiments, the air temperature functions in- and outside the tunnel can be obtained as follows.

The air temperature outside the tunnel can be expressed as follows:

The air temperature in the section that is 3 m from the tunnel portal can be expressed as follows:

The air temperature in the section that is 500 m from the tunnel portal can be expressed as follows:

The air temperature in the section that is 1700 m from the tunnel portal can be expressed as follows:where t is the time (s) and T is the climate temperature in- or outside the tunnel (°C).

5. Discussion and Conclusion

A regressive analysis of temperature field data was conducted based on the long-term field tests at the Qingshashan tunnel, which was located on the eastern Qinghai-Tibet Plateau. The conclusion is as follows: the air temperature outside the tunnel area is lower than that at the tunnel entrance in the whole year, and the average temperature difference outside the tunnel area reaches the maximum when the temperature is the lowest, and the maximum difference can reach 4°C. According to the annual wind direction measurement in the tunnel, the wind direction blows from the exit to the entrance in winter, and the wind speed at the exit of the tunnel can reach 2.8 m/s; the days when the daily average temperature at the exit of the tunnel is lower than 0°C each year are about 180 days, and the days at the entrance of the tunnel are about 165 days. The minimum value of the air temperature curve at 500 m away from the tunnel is 30 days longer than the minimum value of the air temperature curve outside the tunnel. The negative temperature days in the tunnel at 500 m away from the tunnel exit can reach more than 210 days in the whole year, 30 days longer than the total negative temperature days outside the tunnel entrance. In addition to some extremely low temperatures in the middle of the tunnel, the temperature remained above 0°C throughout the year.

During the tunnel operation, because of the increase in the wind speed due to the “piston” effect, the invasion of cold energy is intensified, the temperature in the tunnel environment is decreased, and the ground temperature in the tunnel is decreased. The temperature of the left line is lower than that of the right line, and the inlet of the left line (one side of the air inlet, south side) is significantly lower than that of the right line (one side of the air inlet, south side), while the difference between the outlet of the left line and the inlet of the right line is small.

From the central to the entrance, the temperature in the tunnel decreases by 0.11°C every 10 meters along the longitudinal direction; from the central to the exit, the temperature in the tunnel increases by 0.07°C every 10 meters. Those can be expressed by correlation function relations. The National Highway 317 Line Lushan tunnel and the Dabanshan tunnel also showed similar function relations, both of them belonged to long tunnels. The law of the temperature field in the cold-region tunnel with similar altitude, length, wind direction, and wind speed is similar. More similarities and differences are shown in Table 5. The research results provided fundamental dates for the investigation of temperature fields and antifreeze insulation measures and references that can be used in the construction of cold-region tunnels [32, 33].

There is a cold proof ditch at the entrance of this long tunnel. Considering the extremely low temperature, there are many electric heaters in the cold proof ditch. For tunnels that may suffer from frost damage, the surface insulation treatment shall be adopted, and then the surrounding rock shall be grouted at low temperature to block the fissure water in the bedrock, and the insulation ditch shall be set in the tunnel to drain the water after lining.

When the central drainage pipe is frozen and thawed, it shall be heated in time to ensure the smooth drainage of the tunnel. Through the temperature field test, the tunnel may be frozen in the range of heat preservation treatment, and the melted water on the back of the lining shall be removed in time, and in principle, the frozen area shall not be heated actively.

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 work was financially supported by the Fundamental Research Funds for Ministry of Transport of the People’s Republic of China (no. 2015319812140), Project on Social Development of Shaanxi Provincial Science and Technology Department (nos. 2018SF-382 and 2016SF-412), and Open Fund for Shaanxi Key Laboratory of Geotechnical and Underground Space Engineering (no. YT201905).