Experimental Investigation of Ground and Air Temperature Fields of a Cold-Region Road Tunnel in NW China

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. /e conclusion is as follows: the temperature outside the tunnel and the ambient temperature are affected by wind speed and light. /e 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.8m/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. /e 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.


Introduction
Permafrost regions account for approximately 23% of Earth's surface, mainly in Russia, Canada, China, and the US state of Alaska [1][2][3]. 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 centrenorth. 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 [4][5][6]. 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. e 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.
e Yuximolgai tunnel built in Xinjiang have ice jam due to the impact of low temperature and frost damage. e 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]. us, 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 [10][11][12], 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 [13][14][15]. e 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. rough 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 coldregion 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. e 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.

Project Overview
e Qingshashan tunnel on the Tibetan Plateau, China ( Figure 2), was constructed as part of a highway improvement program from Pin'an to Adai. e tunnel is a separated tunnel with traffic lanes that run in opposite directions. e Qingshashan tunnel is the longest highway tunnel in Qinghai Province. e tunnel's two traffic pathways are 3350 m and 3340 m long, respectively. e 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. e 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. e climate is that of a typical continental semiarid plateau. e Qingshashan tunnel is in an area in which the seasonal frost depths are between moderate and deep and receives moderate frost damage [30]. e 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. e 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  Advances in Civil Engineering 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. e tunnel excavation area is rich in water, and there is a risk of frost damage in the tunnel.

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. at 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. e 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. ermometric 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. eir locations are shown in Figure 5. e 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. e temperature was monitored every six hours: 2:00 am, 8:00 am, 2:00 pm, and 8: 00 pm.

Air Temperature Monitoring.
e air temperature inand 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. e 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]. e longitudinal distribution of the thermohygrometers along the tunnel is shown in Figure 7       year. e 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. e 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. e third stage is from the middle of August to November. e rate of temperature decrease is faster than that of temperature rise, so the curve of temperature decrease is steeper than that of temperature rise.

Discussion of the Ground Temperature
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.
e 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. e 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. e 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. e third stage is from the middle of July to November. e 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. e lowest temperature of the testing sections and the corresponding time is listed in Tables 2 and 3. e lowest temperature of section K31 + 668 is lower than that of section K33 + 218, which is contrary to what is generally assumed. e 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.  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. e 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.

Longitudinal Distribution of Ground
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. e water in the fault originates in a surface water source. e 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. us, the temperature curves tend to be more intensive from the exit to the entrance. (2) e 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. erefore, the wind direction   Advances in Civil Engineering in the tunnel is the key parameter that affects the ground temperature distribution.

Radial Distribution of Ground
Temperature. e 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. e 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. e 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. e 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. e two measured temperature curves are highly regular. us, they can be used to predict the average temperature of the rock masses. Figure 11, the air temperature near the tunnel entrance is   Figure 11: e air temperature curves outside the tunnel vs. time. 8 Advances in Civil Engineering lower than that near the tunnel exit throughout the year. e 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. e 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. roughout 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. 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. e 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. erefore, the temperature of the liner and the rock mass at a certain depth is generally below 0°C for a long period.

Air Temperature inside the Tunnel. As shown in
roughout the year, the temperature, which is generally below 0°C, is sustained for 210 days in the section 150 m from the exit. us, 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. Figure 13, the air temperature distribution in the longitudinal direction presents an inverted V-shaped pattern from January to March.

Longitudinal Air Temperature Distribution. As shown in
at 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. e exit air temperature behaves similarly. e 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). e 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.
e air temperature outside the tunnel can be expressed as follows: (1) e air temperature in the section that is 3 m from the tunnel portal can be expressed as follows: e air temperature in the section that is 500 m from the tunnel portal can be expressed as follows: T 500 � 2.9 + 8.9 × sin 2πt 365 e 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 inor outside the tunnel (°C).

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. e 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 Advances in Civil Engineering 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. e 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. e 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. e 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. ose can be expressed by correlation function relations. e National Highway 317 Line Lushan tunnel and the Dabanshan tunnel also showed similar function relations, both of them belonged to long tunnels. e law of the temperature field in the cold-region tunnel with similar altitude, length, wind direction, and wind speed is similar. (1) When the temperature in the tunnel is −12°C, the surface of the concrete upper lining is 0.2°C-0.5°C.
(2) e temperature in the tunnel in the cold region changes regularly along the depth, which is approximately a parabolic distribution.
(1) ere is no pressure cold-proof drain tunnel under the tunnel's main hole.
(2) e tunnel's surface insulation layer is rigid foam profile and rare earth insulation material.
(3) e tunnel outlet insulation door, the tunnel outside the circular end buried thermal insulation outlet.
e analysis shows that the maximum frozen depth of surrounding rock changes little along the whole length of the tunnel with a length of about 600 m, and the frozen depth outside the tunnel is larger than that inside the tunnel.
(1) e antifreeze scheme is designed with drainage as the core, using highquality waterproof materials, as well as wave-shaped longitudinal drainage pipes. (2)  (1) e temperature change in-and outside the tunnel is periodic, which can be fitted by sine function regression method. (2) After the tunnel is through, the average temperature inside the tunnel decreases by 3-4°C, after the tunnel is through, the air flow has an impact on the temperature distribution inside the tunnel.
(1) In the length range where frost damage may occur, the heat insulation layer with surface heat insulation treatment shall be adopted. (2) e proper drainage system shall be adopted, the deep buried water or cold proof ditch shall be adopted to prevent the water discharged from the tunnel from overflowing and road freezing in the portal section due to the inability to discharge from the tunnel. e temperature inside the tunnel lining presents an inverted "U" distribution along the longitudinal direction of the tunnel. e temperature rise at the entrance of the tunnel is small, while the temperature rise at the exit is large. e main influencing factor is that the wind speed in the tunnel increases with the depth of the surrounding rock.
(1) Cold proof ditch shall be set at the two tunnel openings, respectively, and the distance from the ditch to the outside of the tunnel shall be more than 20 meters. (2) In addition, multiple sets of electric heaters are set in the cold proof ditch to relieve the freezing of drainage channel. (3) Within the length of the hole where freezing damage may occur, the surface heat insulation layer shall be used.
More similarities and differences are shown in Table 5. e 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]. ere 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. rough 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
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.