(中文)严寒地区隧道围岩冻融状况分析的导热与对流换热模型.doc

严寒地区隧道围岩冻融状况分析的导热与对流换热模型 何春雄 吴紫汪 朱林楠 （中国科学院寒区旱区环境与工程研究所冻土工程国家重点实验室） （华南理工大学应用数学系） 摘 要 通过对严寒地区隧道现场基本气象条件的分析，建立了隧道内空气与围岩对流换热及固体导热的综合模型;用此模型对大兴安岭西罗奇2号隧道的洞内气温分布进行了模拟计算，结果与实测值基本一致;分析预报了正在开凿的祁连山区大坂山隧道开通运营后洞内温度及围岩冻结、融化状况. 关键词 严寒地区隧道 导热与对流换热 冻结与融化 在我国多年冻土分布及邻近地区，修筑了公路和铁路隧道几十座.由于隧道开通后洞内水热条件的变化;，普遍引起洞内围岩冻结，造成对衬砌层的冻胀破坏以及洞内渗水冻结成冰凌等，严重影响了正常交通.类似隧道冻害问题同样出现在其他国家（苏联、挪威、日本等)的寒冷地区.如何预测分析隧道开挖后围岩的冻结状况，为严寒地区隧道建设的设计、施工及维护提供依据，这是一个亟待解决的重要课题. 在多年冻土及其临近地区修筑的隧道，多数除进出口部分外从多年冻土下限以下岩层穿 过.隧道贯通后，围岩内原有的稳定热力学条件遭到破坏，代之以阻断热辐射、开放通风对流为特征的新的热力系统.隧道开通运营后，围岩的冻融特性将主要由流经洞内的气流的温度、速度、气—固交界面的换热以及地热梯度所确定.为分析预测隧道开通后围岩的冻融特性，Lu-nardini借用Shamsundar研究圆形制冷管周围土体冻融特性时所得的近似公式，讨论过围岩的冻融特性.我们也曾就壁面温度随气温周期性变化的情况，分析计算了隧道围岩的温度场[3].但实际情况下，围岩与气体的温度场相互作用，隧道内气体温度的变化规律无法预先知道，加之洞壁表面的换热系数在技术上很难测定，从而由气温的变化确定壁面温度的变化难以实现.本文通过气一固祸合的办法，把气体、固体的换热和导热作为整体来处理，从洞口气温、风速和空气湿度、压力及围岩的水热物理参数等基本数据出发，计算出围岩的温度场. 1数学模型 为确定合适的数学模型，须以现场的基本情况为依据.这里我们以青海祁连山区大坂山公路隧道的基本情况为背景来加以说明.大坂山隧道位于西宁一张业公路大河以南，海拔3754.78~3801.23 m，全长1530 m ,隧道近西南—东北走向. 由于大坂山地区隧道施工现场平均气温为负温的时间每年约长8个月，加之施工时间持续数年，围岩在施土过程中己经预冷，所以隧道开通运营后，洞内气体流动的形态主要由进出口的主导风速所确定，而受洞内围岩地温与洞外气温的温度压差的影响较小；冬季祁连山区盛行西北风，气流将从隧道出曰流向进口端，夏季虽然祁连山区盛行东偏南风，但考虑到洞口两端气压差、温度压差以及进出口地形等因素，洞内气流仍将由出口北端流向进口端.另外，由于现场年平均风速不大，可以认为洞内气体将以层流为主 基于以上基本情况，我们将隧道简化成圆筒，并认为气流、温度等关十隧道中心线轴对称，忽略气体温度的变化对其流速的影响，可有如下的方程: 其中t为时间，x为轴向坐标，r为径向坐标;U, V分别为轴向和径向速度，T为温度，P为有效压力(即空气压力与空气密度之比少，V为空气运动粘性系数，a为空气的导温系数，L为隧道长度，R为隧道的当量半径，D为时间长度, 分别为围岩的冻、融区域. ,分别为冻、融状态下的热传导系数，,分别为冻、融状态下的体积热容量， X=(x,r) , 为冻、融相变界面,To为岩石冻结临界温度(这里具体计算时取To=-0.10)，为水的相变潜热. 2 求解过程 由方程(1)知，围岩的温度的高低不影响气体的流动速度，所以我们可先解出速度，再解温度. 2.1 连续性方程和动量方程的求解 由于方程((1)的前3个方程不是相互独立的，通过将动量方程分别对x和r求导，经整理 化简，我们得到关于压力P的如下椭圆型方程: 于是，对方程(1)中的连续性方程和动量方程的求解，我们按如下步骤进行: (1)设定速度,; ( 2)将,代入方程并求解，得 (3)联立方程(1)的第一个和第二个方程，解得一组解,; (4)联立方程((1)的第一个和第三个方程，解得一组解,; (5)对((3) ,(4)得到的速度进行动量平均,得新的,返回(2) ; (6)按上述方法进行迭代，直到前后两次的速度值之差足够小.以,,作为本时段 的解,下一时段求解时以此作为迭代初值. 2. 2 能量方程的整体解法 如前所述，围岩与空气的温度场相互作用，壁面既是气体温度场的边界，又是固体温度场的边界，壁面的温度值难以确定，我们无法分别独立地求解隧道内的气体温度场和围岩温度场.为克服这一困难，我们利用在洞壁表面上，固体温度等于气