寒区涵洞建设是路基工程中必不可少的一部分,它的设置改善了冻土地区路基横向排水条件,优化了交通路线,但同时也改变了路基原有的水、热平衡状态,对路基的热稳定性有一定的影响。到目前为止,针对多年冻土区涵洞对路基热稳定性的影响,很多学者展开了研究,也取得了很多成果:俄罗斯专家Passek等综合运用数值模拟和现场观测方法,研究了大口径波纹管涵在道路工程中的应用;Kunitake等[1]采用三维温度场的有限元法对一混凝土箱涵进行了数值模拟;王凯等[2]针对青藏铁路自身所处的特殊自然条件,对涵节预制、涵节运输、涵洞基础施工和涵节现场拼装等技术作了技术改进,克服了现浇涵洞中所遇到的一些困难缺点;吴少海等[3]根据青藏公路涵洞的经验,结合青藏铁路的特点,提出了适合多年冻土区的桥涵结构、支座类型、墩台结构形式、基础类型、混凝土材料及其施工方案;张学富等[4]通过对青藏铁路某处涵洞现浇混凝土基础水化热进行数值分析,发现现浇混凝土水化热对涵洞周围冻土的热状况有较大影响,并计算分析了在涵底基础下铺设保温层时冻土热状况;张旭芝[5]通过对青藏铁路沱沱河试验段两座试验涵洞进行地温现场监测和对观测数据的分析,研究了适合青藏高原特殊施工环境的涵洞施工工艺、最佳施工季节、施工对多年冻土的影响以及沿涵洞轴向多年冻土上限的变化特征;孙文等[6]运用渗流理论和传热学理论,计算并讨论了渗流对涵洞地基温度场的影响以及涵洞对路基热稳定性的影响。由于早期季节冻土区工程经验不足,因此有关季节冻土区涵洞对路基热稳定性的影响研究不多。
哈尔滨-大连高速铁路(简称“哈大高铁”)是我国在严寒地区设计、建造的第一条铁路客运专线,该线穿越东北地区中部山前平原,处于中深季节冻土区,随着它的建成通车,为铁路建设提供了大量宝贵的资料和经验,基于此,许健等[7]运用数值方法模拟分析了季节冻土区保温路基设计参数;Liu等[8]分析了季节性冻融路堤和地基不同的温度变化过程,牛富俊等[9]对哈大高铁某涵洞的路涵过渡段路基土体的冻结特征展开监测,并分析了路涵过渡段路基及其下地基土的冻结特征和时空变化过程及其对路基稳定性的影响,取得了一定的成果,但对于季节冻土区涵洞的研究还可以更加深入。
本文以哈大高铁典型涵洞为原型,以热学理论为基础,建立了季节冻土区涵洞的温度场计算模型,并与监测数据进行对比分析该模型的可靠性,通过改变参数和尺寸,对不同模型进行数值计算,分析了冷季涵底基础的对流效果、地基土含水量和涵洞尺寸对最大冻结深度发展的影响,对深季节冻土区涵洞工程建设和维护具有实际意义。汪江红,等:深季节冻土区涵洞热状况分析
哈大客运专线涵洞防冻胀试验研究试验场地位于吉林省中西部的公主岭市西北郊范家屯附近,海拔204 m,属于温带半湿润季风气候,其特点为:冬冷夏热,温度、雨量、光照等季节性变化显著,春季干旱多大风,回暖迅速;夏季热而多雨;秋季温暖多晴朗天气;冬季漫长而寒冷。年平均气温5.6 ℃。据气象资料及现场监测的天然场地冻结深度变化数据显示,最大冻结深度为1.56 m。
现场试验涵洞设计里程为DK673+820,是1孔2.0 m的钢筋混凝土箱型涵。该涵为排地表水而设,与线路正交。上下游均采用标准铺砌,翼墙采用直径ϕ12 mm、间距15 cm的护面钢筋,基础底地基处理方法同相邻路基处理。上游侧路基边坡上设置简易检查梯一座。
试验涵洞处的底层岩性特征为:
黏质黄土:厚5~16 m,分布于地表,软塑为主,局部硬塑,Ⅲ级硬土。
粗砂:厚1.0~3.0 m,透镜状分布于黏质黄土之下,中密,饱和,Ⅰ级松土。
泥岩夹砂岩:成岩作用差,全风化泥岩岩芯呈泥柱状、土状,强风化泥岩岩芯呈块状、碎石状、短柱状,锤击易碎,风化层厚度大于10 m,Ⅲ级硬土,弱风化,泥岩为主,泥质胶结,Ⅳ级软石,为极软岩,泥岩具有弱膨胀性。
地下水位埋深1~4 m,地下水对混凝土不具有侵蚀性,监测资料显示,地下水位在3到4月份保持在最低值,之后随降水量的增加逐渐上升,到次年7~8月降水量最大时达到最高位,然后随降水量的减少而缓慢下降,并保持水位在1 m以内到11月底,之后,降水量很小,地下水位迅速下降,并于3月初达到最低值。
试验涵洞段的路堤表层上是厚0.1 m的强化基床表层混凝土、其下是0.65 m的级配碎石砂砾土、1.00 m的非冻胀性A/B组填料、1.30 m的普通A/B组填料、基床以下路堤是1.66 m的普通A/B/C组填料,再下面是0.50 m的褥垫层和黏质黄土地基。路涵过渡段路基采用的是一段式过渡,涵洞高2.50 m,涵侧混凝土壁厚0.22 m,涵顶和涵底厚0.28 m,涵身以下设置了厚1.02 m的涵底地基,总厚度达1.30 m,各构件组成图如图 1所示。
目前有关冻土区温度场数值模型的研究比较多[10-14],本文采用了一种较成熟的数值模拟方法,旨在研究深季节冻土区涵洞模型参数的选取以及分析各参数对其温度场的影响。
控制微分方程及其有限元公式
只考虑热传导和相变的情况下,温度场的数学模型方程为[15]
式中:C、λ、ρd分别为模型中各材料的容积热容量、导热系数和干容重,Wi为含冰量;L为冰、水的相变潜热。
在冻结过程中,未冻水含量是温度的函数[16],不考虑水分迁移,则
将式(2)带入式(1),简化为
其中C为模型材料的视热容量。
假设剧烈相变发生在(Tm-ΔT, Tm)范围内,且在这范围内,体积热容量C及导热系数λ均随着温度线性变化[9-10]。则
冻融过程的热传导是一个二维变系数非线性问题,采用具有迭代的有限元法,并通过减少时间步长提高精度,计算模型的有限元公式导出为
式中:Δt为时间步长;T0、T1分别表示每个计算步长开始之前和结束时温度列向量;K表示导热系数矩阵;M表示总变温矩阵;Q1与边界条件有关的热流量矢量列阵。
计算时,首先根据上个时间步的结果Tn-1做为初始温度条件,确定各单元的初始状态,材料视热容量C1n和导热系数λ1n,由此计算T1n,再根据T1n计算C2n, λ2n, T2n,不断重复,直到迭代i次以后使得|Tin-Ti-1n|小于给定的精度就停止该时间步,根据这一步计算下一步,周而复始,直到所要计算的时间点为止[15]。
采用Geo-Studio软件中TEMP/W(地下热传递分析软件)进行模拟计算,Geo-Studio是一套适用于岩土工程和岩土环境模拟计算的仿真软件。在TEMP/W软件中,不再需要编制繁杂的求解器,它含有全面的地热模型,用户可以根据自己的需要设定,能对水转化成冰和冰转化成水的过程中释放和吸收的相变潜热问题进行计算分析,并且将其表示为未冻水含量的函数。TEMP/W是一种对各种不同地面条件进行分析的强大软件。
试验涵洞在几何形状上是关于横截面和中截面对称的,忽略阴阳坡效应等一些不对称因素,将涵洞及其过渡段简化为一半来计算。根据《高速铁路设计规范(试行)》(2009年),综合考虑现场实际情况(图 1),将计算区域延伸到天然地表以下14 m,沿路基方向取25 m,模型示意简图如图 2。
由于涵顶在2009年底做了堆载预压处理,于2010年6月底卸除完毕,并从2010年7月开始以每周一次的频率进行正式的地温监测,所以以2010年7月10日这一天的实测温度场数值作为初始温度场。
涵洞及其过渡段温度场初始温度条件为:T(t)|t=0=T0。
监测结果显示涵洞处深部地温为8 ℃左右,且年变幅只有0.1 ℃,因此取模型的下部边界条件为:T=8 ℃。
左边界是在对称轴上,两侧无热量交换,当轴向距离足够时,可以认为右边界两边路基温度没有变化,取模型的左右边界条件为:$\frac{{\partial T}}{{\partial n}} = 0$。
对于上边界和涵洞内壁温度边界条件,为了使得计算模型更加符合实际情况,减小由于边界条件的选择而产生的误差,拟采取实测温度(如图 3)进行计算。
涵洞内壁没有太阳直射,取室外大气温度作为计算温度,室外温度理论上应该是一天内各个时间段气温的加权平均值,采用平均温度之半加上最低温度之半所得的值作为涵洞内壁温度边界条件的计算数值,对于上边界混凝土路面的温度边界条件,采用地表的实测温度值作为上边界的温度边界条件。
为了方便模型的建立,引入如下基本假定:
1) 假定在每次的模拟过程中涵洞及路基填料无水分补给和排水作用,即忽略水分迁移,这样可以将含水量固定为一个常量,同时也是考虑到试验涵洞现行的排水措施的有效性。
2) 假定填料及地基土介质为刚性体,在冻结过程中土体不发生变形。
3) 各区域的材料均匀且各向同性。
4) 考虑热传导和相变作用,忽略填料内部热对流的影响。
5) 土层热物理参数的选择,由于温度对填料的体积热容量影响比较小,为了简便起见,计算过程中,土体在未相变区和冻结相变区的体积热容量只取其融、冻两种状态的数值,即融土和冻土,同样,忽略温度对导热系数的影响,在未相变区和冻结相变区土体的导热系数只考虑冻融状态[17-18]。
2011年春季以后,该处涵洞进行了修补操作,在涵底加铺了20 cm的混凝土层等,考虑到水化热等影响较大,本文所要研究的主要内容是围绕最大冻结深度等内容展开,而4月份以后气温已经升到0 ℃以上,所以所模拟的时间区间为2010-07-10到2011-04-02。
根据现场实测资料和有关参考资料的研究[7, 18],土层的热物理参数见表 1。
为了验证计算模型的可行性,取涵底下0.3、0.8、1.3、1.8和2.3 m各位置上的温度变化情况(如图 4中的模拟值1),与实测结果分别进行对比。
对涵底以下各探头处实测值与计算值进行对比,可以看到,按表 1所示的热物理参数计算,暖季涵下温度一致性比较好,但是冷季温度远高于监测数据,从而由此得出的冻结深度也远小于监测数据。考虑到涵底基础是由片石混凝土或浆砌片石铺砌,涵底在暖季热传导起主流作用,在冷季热对流会对其产生明显影响,所以在模拟过程中忽略冷季热对流效果会导致很大误差,假定4)并不适用。因此假设涵底基础冷季当量导热系数为暖季的2倍(具体热物理参数如表 2所示),并进行计算对比,对于不同涵底基础热对流效果对涵底最大冻深的影响在后面进行讨论。
图 4中模拟值2是假设涵底基础冷季当量导热系数为暖季的2倍时涵底下不同深度探头的温度计算值,对比实测值可以看到,在模拟时间段内,各探头点的模拟计算结果与实测数据吻合的相当好,误差都比较小,该模拟结果不仅在温度上升下降的整体趋势是与实测值相符合,在对涵洞进行定量计算都有一定的可靠性。
为了进一步分析模拟结果与实测结果之间的差值随深度的变化规律,取2010-12-04和2011-02-26涵底不同位置温度示意图(图 5),可以看出进行参数修正以后,模拟精度得到提高,尤其在冻结层附近的温度拟合精度较高。
根据第二次模拟结果,对涵洞的热状况进行分析。
如图 6所示,取路基下0.65和3.25 m距涵洞中心不同距离处的温度进行对比,可以看出涵洞对涵周路基温度分布有显著影响,加大了路基温度随气温变化而改变的幅度,对气温改变的响应更为迅速,而且越接近涵壁,这种影响更为剧烈。
如图 7所示的涵顶填料中不同位置温度时程曲线图,可以看到在高温段和低温段,涵顶温度发展较均衡,温差较小,而且由于太阳辐射的作用,从上往下温度降低,另外在降温和升温过程中,分别表现出上下低于中心和上下高于中心的特点,双向冻结和双向融化明显。
为便于对涵洞下冻结深度的发展过程有个直观了解,根据有限元模型的计算结果,绘制了涵下基础中不同深度温度时程曲线(如图 8)。其中横坐标表示时间历程,纵坐标表示从涵洞底向下的深度。
计算模型温度边界条件采用的是实测温度,在冻结初期,由于温度的扰动产生了一段波动,但随着温度整体趋势趋于下降,涵洞下土体正式进入冻结阶段,正式冻结开始后,涵底1 m范围内温度下降迅速,在1周内即整体下降到0 ℃以下,然后趋于平缓,以一定的速率向下发展,经历大概100 d,达到最大冻结深度;之后进入融化阶段,冻深慢慢减小,仔细观察包络线,会发现融化结束前0 ℃线几乎是垂直的,表明涵洞中心下温度梯度很小。
如图 9所示的涵周路基冻结圈发展历程,在冻结期内,涵洞周围路基内形成的冻结圈呈“倒钟”形。涵周填土首先形成冻结圈,然后涵顶双向冻结,涵底单向冻结,第4 d涵顶填料就整体冻结,形成“倒钟”形冻结圈,再经过大约90 d,于2011-02-22形成最大冻结圈;再经过40 d左右,涵顶和大部分涵下冻层消失,一个月以后,于2011-05-06冻层彻底消失,这些与现场实测数据体现的冻结圈发展历程基本吻合。涵底中心冻结层最大厚度为2.26 m。
综上对比分析计算结果和实测数据,可以看出,这种根据试验涵洞实际铺砌情况及实测温度建立的涵洞段温度场模型,不仅在定性分析涵洞热状态变化趋势上与现场实际情况有较好的吻合,而且在定量计算上都有一定的可靠性。这说明在温度场模型建立过程中,应该根据填料热对流实际情况,确定其冷季当量导热系数,由此建立模型计算结果才能够准确的反应温度状况,对其它涵洞工程的设计具有一定的参考价值。
涵底基础是由片石混凝土或浆砌片石铺砌的,在前面的分析讨论中看到,涵底基础的冷季对流效果明显,加深了涵底最大冻结深度。
在多年冻土区,碎石通风路基是一种积极保护多年冻土的工程措施,工作原理是:在寒冷季节,冷空气有较大的密度,在自重和风的作用下使片石间隙中的热空气上升,冷空气下降并进入地基;而在温暖的季节,热空气密度小,很难进入地基,类似于热开关效应。碎石路基在多年冻土区对于保护地基土冻结状态具有优势,被国内外广泛应用于冻土区房屋建筑,相关资料指出,某试验中块石层的有效导热系数在冬季是夏季的12.2倍[19]。
但在季节性冻土区,由于防止冻胀灾害关键在于减小冻结深度,减小冻胀量,而非多年冻土区的保护冻土层,防止冻土层的退化。在上面建立的涵洞温度场计算模型中,考虑到对流作用,对涵底基础取2倍于暖季导热系数的当量导热系数为其冷季热参数值,与实测数据的比较也能看出,这种处理办法是符合客观规律的。
为了讨论季节性冻土区涵底碎石对流效果对涵底最大冻结深度的影响,采用上面建立的温度场模型,取冷季当量导热系数分别为暖季的1~12倍,其它参数不变进行计算。
如图 10所示最大冻结深度与冷、暖季涵底基础当量导热系数比值之间的关系,可以用下面的函数公式进行较好的拟合:
随着涵底碎石冷、暖季当量导热系数比值的增加,涵底冻结深度呈增函数形式变化,而且增长先快后慢,当λc/λw=2时,最大冻结深度较没有热对流时增加了0.38 m,从而看出,对于深季节冻土区,碎石层的对流传热对防止冻胀灾害是极其有害的,所以在季节性寒区涵洞设计工作中,应该避免一般涵洞设计中,涵底基础统一用片石混凝土或浆砌片石,应该采用新工艺,减少碎石中的孔隙,或者减少孔隙连通性,减少对流传热对涵底冻深发展的影响,使得涵底碎石基础在稳固路基的同时,还能够起到很好的保温效果。
为了研究地基土含水量对最大冻深的影响,分别取不同含水量的地基土所对应的热学参数(如图 11、12所示),带到上面有限元模型中进行计算。
不同体积含水量地基土涵底最大冻结深度变化趋势如图 13所示,随着含水量的增加,最大冻结深度的变化规律大致可分为3段:1)当含水量W < 17.5%,冻深随着含水量的增大而减少;2)当17.5% < W < 27.5%时,涵下最大冻结深度减小速度减缓,冻深减小速度减缓,接着稍稍有些加深;3)当W > 27.5%时冻深随着含水量的增加而减小。
分析图 11、12中融、冻地基土的导热系数、体积热容量与含水量之间的关系,可以看到,融土导热系数随含水量增加的变化规律是前面增加的快,后面增加的慢,冻土的变化曲线分3个阶段,两头缓中间陡;融土的容积热容量随含水量呈直线增大,冻土的容积热容量随含水量呈折线增大,整体上冻土容积热容量随含水量增大的斜率较融土的缓[18]。据此可以分析最大冻结深度随含水量变化规律。
第1阶段,干燥即结合水未达到最大值时,这段时间内,冻、融土的导热系数与体积热容量基本上都呈线性变化趋势,同时,随着含水量的增加,潜热变大,从而导致涵底最大冻结深度不断减小;第2阶段,冻、融土的体积热容量依然呈线性变化趋势,融土的导热系数随含水量的变化依然保持第1阶段的增速,但冻土的导热系数随着含水量的变化快速增大,从而导致在冷季,地基土中能量大量流失,从而导致涵下最大冻结深度减小速度减缓,随着含水量的持续增加,冻深稍稍有些加深;第3阶段,随着含水量的增加,冻、融土的导热系数和体积热容量的变化规律又趋于一致,由于潜热的增大,涵底最大冻结深度又开始不断减少。这与冻土导热系数随含水量变化规律分3个阶段是一致的。整体上看,涵底最大冻结深度随着含水量的增加呈减小趋势。
在实测中,涵底最大冻深为天然地层的两倍左右,究其原因,通过上面的计算讨论,最有可能是有两个方面:1)冷季涵底基础热对流效果明显,当量导热系数远大于天然地层,散热严重;2)由于天然地层的含水量较涵底碎石基础的要高,根据最大冻结深度随含水量变化规律,天然地层的冻深较浅显而易见。涵底冻深较大,一方面体现出涵洞及路基内采取的隔水措施是有效的;另一方面,也表明当区域地基土是冻胀敏感性土,且隔、排水措施完善,人为冻深也需要考虑的相对较深。
试验涵洞为钢筋混凝土箱型涵,为了排地表水而设,根据《铁路桥涵设计基本规范》[20],泄水涵洞的孔径不宜小于2.0 m,涵洞顶至轨底的填方厚度不小于1.2 m,并且对涵洞净高与涵洞长度做出了规定。为了研究涵洞几何尺寸对季节冻土区涵底最大冻结深度的影响,利用上面的模型,针对不同尺寸的涵洞计算其涵底最大冻结深度,为了满足排水需求,而且在规范要求内,将涵洞孔径和净高均控制在2.00~2.95 m之间,同时计算过程中不考虑强度、稳定性等要求,仅仅讨论净高、内径的改变对涵底冻深的影响,求得最合适的涵洞尺寸。计算结果如图 14所示。
从图 14可以看出,箱涵的净高对涵底最大冻结深度的影响很小,当内径恒定为2 m,箱涵的净高从2 m增加到2.95 m时,最大冻结深度仅增加了4.5 cm,增加率仅为2.22%;相较于净高,内径对涵底最大冻结深度的影响相当明显,随着内径的增大,最大冻结深度呈线性加深,当净高恒定为2.95 m,内径从2 m增加到2.95 m时,最大冻结深度增加了0.251 m,增加率为12.23%。所以在寒区涵洞的设计过程中,在满足强度、稳定性、功能性等一系列要求外,尺寸方面优先考虑内径,取允许的最小内径,由于净高对涵底冻结深度方面的影响很小,可以不考虑。
1) 根据试验涵洞实际铺砌情况及现场气象资料,建立了考虑相变潜热的涵洞段温度场模型,对比分析计算结果和实测数据,表明考虑冷季热对流效果后,计算结果不仅在定性分析涵洞热状态变化趋势上与现场实际情况有较好的吻合,而且在定量计算上都有一定的可靠性,能够准确的反应温度状况,对季节冻土区其它涵洞工程设计有一定的参考价值。
2) 在多年冻土区,忽略对流传热的计算是基于一种不利的情况进行计算,对于工程安全稳定趋于保守的,但在季节冻土区,忽略对流传热的影响相当于是基于一种有利的情况进行计算,得到的人为冻结深度要浅于实际人为冻结深度,由此指导设计的涵洞工程就不能达到预定的安全系数。
3) 在深季节性冻土区涵洞设计工作中,冷季涵底基础对流传热的影响不容忽视,应该避免一般的涵洞设计中,涵底基础统一用片石混凝土或浆砌片石;应该采用新工艺,减少碎石中的孔隙,或者减少孔隙连通性,从而减小对流传热对涵底冻深发展的影响,使得涵底碎石基础在稳固路基的同时,还能够起到很好的保温效果。
4) 涵底最大冻结深度随着含水量变化规律呈三段:当含水量W < 17.5%,冻深随着含水量的增大而减少;当17.5% < W < 27.5%时,冻深减小速度减缓,接着稍稍有些加深;当W > 27.5%时冻深随着含水量的增加而减小。涵底最大冻结深度整体上随着含水量的增加呈减小趋势。据此,在隔、排水措施完善的区域,人为冻深相较更深,需要给予更多注意。
5) 研究了寒区涵洞尺寸对涵底冻深的影响,分别计算了不同净高和不同内径下涵底最大冻结深度值,通过对比可以看出其内径对涵底冻结深度方面的影响很大,在满足强度、稳定性和功能性等一系列要求后应该优先考虑,而净高对涵底冻结深度方面的影响很小,可以不考虑。