【摘要】 通过对严寒地区隧道现场基本气象条件的分析 ,建立了隧道内空气与围岩对流换热及固体导热的综合模型 ;在隧道内空气以紊流为主的情况下 ,分析预报了祁连山区大坂山隧道开通运营后洞内温度及围岩冻结、融化状况 ,并与在隧道内空气以层流为主的情况下的已有结论进行了比较 .更多还原

【Abstract】 Based on the analyses of fundamental meteorological and hydrogeological conditions at the site of the Dabanshan Tunnel in the Qilian Mountains, a combined convection-conduction model was constructed for turbulent air flow in the tunnel and temperature field in the surrounding rock. Then, with the in situ conditions of air temperature, atmosphere pressure, wind force, as well as the hydro-thermal conditions, the relationship between the temperature on the surface of the tunnel wall and the air temperature at the entry and exit of the tunnel has been obtained, the freeze-thaw conditions in the surrounding rock wall of the tunnel is predicted, and the simulated result is compared with that obtained in the case of laminar air flow in the tunnel. Many tunnels, constructed in the cold regions of the Tibetan Plateau, are at an elevation above 4 000 m and nearby the ridge of mountain. Since there is cold air flowing in the cold regions almost all the time and there is no wind protective screen, the wind-velocity in situ of the tunnels is more than 5 m/s, and because of the difference of atmospheric pressure, air temperature between the entry and exit of the tunnel, the air flow in the tunnel would be turbulent. In order to deal with the complex random unsteady nature of the turbulent air flow, the Reynolds time-average equations were used, that is, based on the classical Navier-Stokes equations, introduced the pulsate kinetic energy equation (K-equation) and dissipativity equation (ε-equation), then, by the Boussinesq assumption, the algebraic relationship of the turbulent kinematic viscosity v t , pulsate kinetic energy K and dissipativity ε are obtained. Combining these equations about the turbulent air flow in tunnel with the convective and conductive equations about the air temperature in the tunnel and about the temperature field with phase-change in surrounding rock wall, the whole mathematical model with a system of equations was constructed. In order to predict the freeze-thaw conditions for the Dabanshan Tunnel, the parameters about the turbulent air flow in the model are chosen by the routine methods in air fluid dynamics, and the thermal parameters and initial and boundary conditions in the model are defined as follows. The air density ρ=0.774 kg/m 3 , the thermal capacity of air C p =1.874 4 kJ/kg·K, heat conductivity λ=2.0×10 -2 W/m·K and the dynamic viscosity μ=9.218×10 -6 kg/m·s, the thermal diffusivity α=1.378 8×10 -5 m 2 /s and the kinematic viscosity v=1.19×10 -5 m 2 /s. In the surrounding rock wall, the dry volumetric weight γ d =2 400 kg/m 3 , the content of water and unfrozen water in rock are 3% and 1%, respectively, and the thermal conductivity λ u =1.9 W/m·K, λ f =2.0 W/m·K, heat capacity C v =0.8 kJ/kg·K, and C f =0.8+2.1(V-V u )+4.182V u 1+V×r d , C u =(0.8+4.182V)1+V×γ d .  The wind speed at the entry and exit is approximated as (t)=[0.028×(t-7) 2 +2.5] (m/s), where tis the t-th month in a year. The initial wind speed in the tunnel are set to be  U(0,x,r)=U a (1-(rR) 2 ), V(0,x,r)=0 The effective air pressure p=0 at the entry of wind andp=(1+βL2R)×U 2 2 at the exit,where β is the coefficient of resistence along the tunnel wall, andpr｜ r=R =0. The kinetic energy K=0.01×(U 2 +V 2 )/2 and the dissipativity ε=C 3/4 v K 3/2 0.4(R-r) at the entry, and K and ε are calculated by the local one direction method at the exit, and K=0, ε=0 on the surface of surrounding rock wall. The initial and boundary values of temperature Tare set to be  T(t,0,r)=T(t,L,r) =A+Bsin(2π8760t-π2)℃, T(0,x,R 0 )=[f(x)-R 0 ]×0.03 ℃; T(0,x,r)=(f(x)-r)×0.03,R<r≤R 0 A-B, r≤R. where f(x) is the distance from the vault to the permafrost base(f(x) is about 190 m in central part of the tunnel and about 40～50 m near the entry and exit), and R 0 =25 m is the radius of controll更多还原