青藏铁路高温多年冻土区典型路基的长期热稳定性研究
2022-02-12徐岳震申明德周志伟李国玉
徐岳震, 申明德, 周志伟, 马 巍, 李国玉
(1.中国科学院 西北生态环境资源研究院 冻土工程国家重点实验室,甘肃 兰州 730000;2.中国科学院大学,北京 100049)
0 引言
青藏铁路是全球海拔最高和里程最长的高原高寒地区铁路,全线运营里程约1 925 km,格尔木-拉萨段约1 118 km,其中在海拔4 000 m 以上约965 km。青藏铁路穿越高原多年冻土区达到632 km,包括:岛状、片状不连续多年冻土区约82 km,大区域连续多年冻土区约550 km,其中高温高含冰量多年冻土区约为355 km[1-4]。青藏铁路在多年冻土区的安全运行首要解决的问题是保证高温高含冰量冻土路基的融沉变形在可控的安全范围内,相关“主动冷却路基”工程措施的目的就是确保多年冻土热稳定性,使得多年冻土上限不发生明显下移,控制冻土地温的升温速率,避免多年冻土区路基产生过大的融化沉降变形和高温冻土的蠕变沉降变形[5-9]。
青藏铁路采用“主动冷却路基”思路积极保护多年冻土,利用块碎石材料、通风管、遮阳棚(板)、热棒等工程措施,通过调整和控制热对流、太阳辐射或热传导的思路,降低传入路基内部的热量,避免路基下多年冻土温度的升高[10-14]。其中,通过调控热传导换热机制的块石基底路基、块石护坡路基、U型块石路基在不同的路段被大量采用,三种类型的块石路基结构在青藏铁路多年冻土区运营里程中占到了60%以上[15-17]:(1)块石基底路基结构是不同高度的路基填料层下铺设一定厚度的块石层基底;(2)块石护坡路基结构是指路基边坡铺设一定厚度的块石层护坡;(3)U型路基结构是综合上述两种路基结构,在路基填料层下铺设块石层,在边坡再铺设一定厚度的块石层护坡。块石路基主动降温的原理主要是由于块石层内存在的强制对流和自然对流,基于冬季的空气对流换热作用和夏季的隔热作用,可以有效地主动冷却块石路基下的多年冻土层[18-20]。但也有研究表明,块石路基在低温冻土区的降温效果较好,而在高温多年冻土区降温效果不太理想[21-23]。此外,青藏高原多年冻土区对近年来全球气候变暖的响应极其敏感,青藏高原冻土工程走廊内线性工程的修筑和运营,又进一步叠加了人类活动对多年冻土地区温度场的影响,并引起了广大多年冻土区的渐进退化,威胁到多年冻土区内的青藏铁路路基热稳定,尤其是位于高温、富冰区的路基地段[4,24-28]。到目前为止,多年冻土区内的青藏铁路已建成并安全运行了16 年以上。在长期的运营过程中,多年冻土区内未采取冷却措施的某些普通路基路段,已出现路基下多年冻土加速退化的趋势,特别是在高温及富冰的多年冻土区[29-31]。对普通路基的维修补强措施主要以块石护坡、护道为主,并在个别地段增设了热棒,对路基下的多年冻土取得了一定的主动降温效果[32-34]。在青藏铁路穿越的个别沙漠化较为严重的多年冻土区段,风积沙堵塞了块石路基的孔隙,降低了块石结构的孔隙率,影响了块石结构调节多年冻土温度的能力[35-37]。此前已有大量文献通过现场监测资料,分析了不同路基结构的热状况[9,12,14,21-23],但这些研究大多基于短期的现场监测数据,长期运营后路基热稳定状态的相关研究仍然较少,考虑到高温冻土区内的多年冻土对温度变化的敏感性更高,因此,对于青藏铁路高温多年冻土区内的不同路基结构,其在长周期条件下的多年冻土路基热稳定性是迫切值得关注的问题。
本文基于青藏铁路长期监测系统所获取的长达15年完整周期的地温监测资料(2006—2020年),分别选取了高温多年冻土区内三个典型路基结构断面(普通路基、U型块石路基、块石护坡路基)的地温监测数据,分析和讨论了15年间三个监测断面的温度场、年平均地层温度和年最高地层温度的长期演化规律以及引起相应变化的主要原因,旨在加深对高温多年冻土区内这几类路基长期热稳定性的认识,为青藏铁路多年冻土区路基的维护与补强提供科学依据及指导意见。
1 监测断面介绍
选用青藏铁路沿线多年冻土区内的三个长期监测断面进行分析,三个监测断面(图1 中P17、P30、P43)均位于高温多年冻土区,年平均地温高于 -1 ℃,且含冰丰富,为饱冰、富冰冻土(以体积含冰量iv划分,20%<iv≤30%,富冰冻土;30%<iv≤50%,饱冰冻土[38])。其中P17 断面位于二道沟西段谷地坡脚,P30 断面位于布曲河西岸阶地,P43 断面位于多年冻土南界附近的托居谷地西侧斜坡,监测断面具体信息如表1。监测断面的路基结构是青藏铁路最常见的三种典型路基结构形式,即普通路基(P17 断面)、U 型块石路基(P30 断面)以及块石护坡路基(P43断面)。从3个监测断面的路基走向来看,路基左侧对应阳坡侧,右侧对应阴坡侧(路基左右侧是以格尔木到拉萨的方向定义)。U型块石路基断面,基底为厚度1.0 m、粒径25~30 cm的块石层,左右两侧护坡厚度分别为1.2 m、0.7 m,碎石粒径10~15 cm。块石护坡路基断面,左右两侧护坡厚度分别为1.0 m、0.8 m,碎石粒径10~15 cm。路基边坡坡率为1∶1.5,在历年维修中,U 型块石路基及块石护坡路基的左右侧补强了宽度2~3 m、粒径25~30 cm 的块石护道,具体维修时间不详,由于未发现补强块石护道导致的左、右路肩地温变化规律的异常,故在分析这三个监测断面时,仍按原路基结构分析各孔的地温变化规律。
表1 监测断面信息Table 1 The information of the observation profiles
图1 青藏铁路及三个监测断面位置Fig.1 Qinghai-Tibet Railway and the selected observation profiles
三个监测断面均布设了天然场地、左路肩和右路肩3 个测温孔,其测温孔布设位置如图2 所示,其中天然场地测温孔距路基左坡脚的垂直距离为20 m,天然场地测温孔深度为18 m,左、右路肩测温孔深度均为20 m。三个测温孔的温度探头自地表到10 m 深度每隔0.5 m 布设一个,10 m 深度以下每隔1.0 m 布设一个,温度探头采用冻土工程国家重点实验室研制的热敏电阻,测量精度为±0.05 ℃,温度数据的记录采用DT500 型自动数据采集仪,测温频率为每日1 次。目前已积累了2005 年10 月—2021 年1 月总计15 年完整周期的长期连续地温监测资料。
图2 监测断面地温孔布设示意图Fig.2 Layout of the ground temperature boreholes at observation profile
2 监测结果与分析
2.1 地温整体变化
图3分别为三个监测断面左右路肩及天然孔的地温等温线图,路基左右路肩地温孔深度自路基表面算起,天然地温孔深度自天然地表算起,以下各图的孔深计算点无特殊说明的均依照此规定。
如图3(a)所示,普通路基左路肩孔的年最大融化深度在2014 年之前增速较为缓慢,自4.8 m 增加至5.6 m,在2015 年后由于年最大冻结深度不能到达年最大融化深度位置,该深度土层出现融化夹层[39-42],随着融化夹层的出现,导致其下的多年冻土加速退化,截至2020年,6年时间内最大融化深度自5.6 m 增加至6.8 m。右路肩孔的年最大融化深度在2014 年之前变化与左路肩基本一致,自4.4 m 增加至5.6 m,在2014 年之后融化速率增大,至2020年,6 年时间内最大融化深度自5.6 m 增加至6.4 m,增大速率略低于左路肩孔。天然孔的年最大融化深度在2016 年前基本保持在4.0 m 深,在2016 年之后逐渐增加,至2020 年约4.8 m 深,同时4.0~6.0 m 深度的地温相对之前的年份有升高的趋势。分析认为2014 年之后左右路肩孔融化速率增大的原因与融化夹层的出现有关。考虑到该路基断面高度仅为0.9 m,这表明由于路基工程及运营的热扰动影响,普通路基断面下的多年冻土人为上限已略低于原天然冻土上限,这不利于路基的长期热稳定。值得注意的是,左路肩的融化深度始终大于右路肩的融化深度,其差值除了铁路运营第一年(2006 年)达到1.0 m 以外,其他年份均小于0.5 m,阴阳坡效应[43-47]对左右路肩融深差异的影响不显著。
图3 三个监测断面的等温线图Fig.3 Isothermals of the three observation profiles:traditional embankment(a),U-shaped crushed rock embankment(b),crushed rock revetment embankment(c)
如图3(b)所示,对于U 型块石路基,左路肩在2015年之前其最大融化深度基本保持在5.6~6.0 m之间,2015 年之后最大融化深度略微减小,保持在4.4~5.5 m 之间。右路肩各年份的最大融化深度基本保持在约4.2~4.5 m 之间,相同年份下右路肩的最大融化深度始终小于左路肩。天然孔2006—2013 年的最大融化深度保持在2.4~4.0 m 之间,2012年后融化速度变快,2014年后的数据缺失。考虑到该断面路基高度为5.7 m,这表明U 型块石路基对多年冻土温度场的主动降温效果明显,人为多年冻土上限已提升至路基本体内,但是阴阳坡效应使得左右路肩下部地温出现较为明显的差值。
如图3(c)所示,块石护坡路基左路肩在2005—2012 年以及2015—2018 年两个连续时间段出现了融化夹层,最大融化深度在7.2~8.2 m,最大冻结深度始终在多年冻土上限之上,且保持了1.0~2.0 m的深度差距。由于块石护坡以及后期维修增加的块石护道的存在,导致左路肩下多年冻土的融化速度没有因为融化夹层的出现而变快,这与普通路基略有区别。而右路肩的最大融化深度保持在4.7~5.3 m 之间,且没有出现融化夹层。天然孔的最大融化深度有逐年增加的趋势,且这种趋势在2016年之后较为明显,最大融化深度维持在3.2~4.6 m。考虑到该断面路基高度为7.0 m,这表明块石护坡路基下的人为多年冻土上限有一定程度提升,右路肩下的人为多年冻土上限已提升至路基本体内。综合对比左右路肩和天然孔的地温等温线表明,块石护坡路基对控制多年冻土上限的退化具有一定作用,但是阴阳坡效应降低了块石护坡路基的降温效果。
2.2 年平均地层温度变化
图4 分别为15 年间三个监测断面左右路肩及天然孔不同深度的年平均地层温度曲线,整体可以看出,浅部地层年平均温度变化幅度较大,这里我们重点关注0 ℃年平均地层温度在地层中的位置及不同地温孔浅部地层的地温状态。
图4 三个监测断面年平均地层温度Fig.4 Mean annual ground temperature of the three observation profiles:traditional embankment(a),U-shaped crushed rock embankment(b),crushed rock revetment embankment(c)
如图4(a)所示,普通路基断面左右路肩和天然孔的0 ℃年平均地层温度分别在6.5 m、6.0 m 和4.0 m处,且三处0 ℃年平均地层温度所对应的地层深度在15 年中基本没有变化。左路肩和天然孔的年平均地层温度最大值出现在0~1.5 m 深的近地表,该深度地层的年平均温度在2006 年时最高,在2006—2008 年下降迅速,2008 年之后基本稳定。右路肩近地表层的年平均温度在2006—2008 年下降迅速,在2008—2012 年趋于稳定,2012 年后出现波动上升,由于普通路基断面高度仅0.9 m,右路肩年平均地层温度最大值对应的地层深度约为1.0~1.5 m,再结合此路基的工程实际,右路肩近地表层的年平均温度的异常波动可能与该侧的侧向水热扰动有关,未来可通过获取现场调查数据等方式进一步开展相关研究。对比来看,排除近地表3 m 内年平均地层温度异常的情况,3 m 以下年平均地层温度变化速率随深度增加而减小,以3 m 深地层的年平均温度为例,左路肩和天然孔年平均地层温度差值最高不超过0.93 ℃,经历了一个先减小后增大的过程,说明普通路基下的土层经历了长时间的温度调整后,并不能保持地温的稳定性。左右路肩年平均地层温度差值最高不超过0.8 ℃,经历了一个先增加后稳定的过程,说明路基边坡热效应导致的阴、阳坡路肩浅层地温差异不明显,这或与该截面较低的路基高度有关。
如图4(b)所示,在2006 年初,U 型块石路基断面左右路肩及天然孔的0 ℃年平均地层温度分别在6.0 m、2.5 m 和1.5 m处,15年间各孔0 ℃年平均地层温度变化差异较大。左路肩0 ℃年平均地层温度所对应的地层深度,近15 年以来没有明显变化,基本保持在6 m 深度。左路肩年平均地层温度大于0 ℃的深度范围内,除2018 年和2020 年数据变化波动较大,其他年份6 m 深度范围内的年平均地层温度变化波动较小,均在0.5~2.5 ℃范围内变化。对于右路肩0 ℃年平均地层温度,对应的深度在0~2.5 m 之内,值得注意的是2006 年地表层年平均温度在6.3 ℃,而其他年份的地表层平均温度均在0~-4 ℃范围内,这是由于青藏铁路运营前的工程施工热扰动对路基地表层的影响还较为显著所致。对于天然孔,由于2014 年之后的数据缺失,图中只给出了2006—2014 年的年平均地层温度数据。2006 年的0 ℃年平均地层温度对应深度为1.5 m,除了2007 年、2009 年和2014 年出现某一深度地层年平均温度超过0 ℃以上,其他年份所有深度的年平均地层温度均小于0 ℃。对比来看,尽管不同年份4 m 以上深度年平均地层温度数据规律性不强,但是在4 m 以下深度的年平均地层温度数据变化较有规律。对于4~8 m 深度土层,同等深度的年平均地层温度,左路肩高于天然孔,天然孔高于右路肩,温差随深度增加而减小。考虑该断面的路基高度为5.7 m,以6 m 深度为例,左路肩和天然孔年平均地层温度差值自2006 年的0.50 ℃逐渐降低至2014年的0.34 ℃,左右路肩年平均地层温度差值自2006年的0.60 ℃逐渐增加至2020 年的0.91 ℃,这表明U型块石路基结构对路基下的浅部土层主动降温效果明显,但对路基阴、阳坡路肩下的土体温差调节效果较差,这可能与两侧块碎石护坡的碎石粒径及铺设厚度等有关系。
如图4(c)所示,在2006 年初,块石护坡路基断面左右路肩及天然孔的0 ℃年平均地层温度分别在7.5 m、3.0 m和1.5 m处,15年间左路肩7.5 m以上深度的年平均地层温度在-0.75~1.6 ℃之间波动,而7.5 m 以下深度的年平均地层温度均在-0.75~0 ℃小范围波动。右路肩在近地表1 m 深度范围内不同年份的年平均地层温度波动较大,特别是地表层从2006 年的-2.0 ℃升温至2020 年的9.2 ℃。天然孔3.5 m 深度以下所有年份的年平均地层温度均在0 ℃以下,而在3.5 m 深度以上2016—2020 年之间的年平均地层温度明显高于其他年份的年平均地层温度,且在同一深度(约0.5 m)出现年平均地层温度峰值,结合文献[48]中的气象数据,这或与近年来的气候变暖有关。考虑该断面的路基高度为7.0 m,以7.0 m 深度地层的年平均温度为例,左路肩和天然孔年平均地层温度差值自2006 年的0.36 ℃逐渐降低至2020年的0.11 ℃,左右路肩年平均地层温度差值自2006 年的0.12 ℃逐渐增加至2020 年的0.46 ℃,这表明块石护坡路基结构对路基下浅部土层具有一定的降温效果,但对阴阳坡效应造成的浅部土层温差增大的情况同样不能有效控制。
为进一步对比和分析三个监测断面深部土层年平均地层温度的发展变化趋势,图5 给出了三个监测断面左右路肩以及天然孔15 m 深度的年平均地温变化曲线。如图5(a)所示,普通路基左路肩年平均地温在2006 年略低于右路肩0.025 ℃,随运营时间的增加,左路肩年平均地温增长速率高于右路肩,至2020 年左路肩年平均地温略高于右路肩0.027 ℃。值得注意的是,到了2020 年左右路肩以及天然孔的年平均地温值基本接近吻合,最大差值为0.028 ℃。一般而言,地温梯度在冻土层上部为负,下部为正(规定地温随地层深度增加时为正,反之为负)[49],即冻土地温在某一深度的年平均地层温度最小,在此深度以下的年平均地层温度随深度增加而增加,这里选取的路基左右路肩孔15.0 m 深度(深度起算点为路基表面)对应天然孔14.1 m深度(天然孔深度起算点为天然地表,路基高度0.9 m),因此,在2006 年初期,天然孔15.0 m 深度的年平均地温要高于路基左右路肩孔15.0 m 深度的年平均地温。铁路的施工及运营给路基下的冻土带来了新增热扰动,随运营时间的增加,这种热扰动由路基浅部土层传递到深层冻土,表现为路基左右路肩孔15.0 m 深度年平均地温的增长速率高于天然孔的增长速率,且左路肩增长速率高于右路肩。
图5 三个监测断面年平均地温变化Fig.5 The variations of mean annual ground temperature in the three observation profiles:traditional embankment(a),U-shaped crushed rock embankment(b),crushed rock revetment embankment(c)
如图5(b)所示,U 型块石路基左右路肩的年平地温均保持在-0.43~-0.67 ℃之间,而天然孔的年平均地温始终保持较高水平,从2006 年的-0.33 ℃到2012年的-0.29 ℃。左路肩年平均地温增速与天然孔的年平均地温增速基本持平,而右路肩年平均地温增速略小于左路肩,这表明地表阴阳坡效应对路基温度场的影响在15 m 深处仍然有体现。尽管此断面2012年之后的天然孔年平均地温数据缺失,但是对比图中三个地温孔的年平均地温数据,仍然可以看出U型块石路基对于路基下多年冻土温度场的控温效果明显,在青藏铁路运营达15年后,U型块石路基对路基下土层的降温效果依然有效。
如图5(c)所示,块石护坡路基在运营的15 年中,左右路肩的年平均地温差始终保持在0.03 ℃以上,其阴阳坡效应较为明显。值得注意的是,从2007年到2016年,左右路肩年平均地温均出现了先增大(2007—2011 年)后减小(2011—2016 年)的特征,两阶段左右路肩的年平均地温差值规律基本与此保持相同趋势,这表明深部土层受阴阳坡效应等外部环境的影响也是同步的。此外天然孔年平均地温始终与左路肩的年平均地温大体一致,表明在铁路施工和运营的初期,块石护坡路基就对冻土起到了保温作用,且随着运营时间的增加,块石护坡路基与天然场地15 m 深度处地层的年平均地温始终保持着大体一致,对路基下的多年冻土具有一定的保温效果。
2.3 多年冻土上限变化
图6 分别为三个监测断面各地温孔15 年来的年最高地层温度曲线图,选取地层每一深度每年出现的年最高温度,作为反映多年冻土上限变化的参考。如图6(a)所示,普通路基人为多年冻土上限深度变化情况为,左路肩从2006 年的5.7 m 下降到2020 年的6.5 m;右路肩从2006 年的4.7 m 下降到2020 年的6.3 m。与左路肩相比,尽管青藏铁路运营的第一年(2006 年)人为多年冻土上限深度有1.0 m 的差距,但是在运营15 年以后人为多年冻土上限深度只有0.2 m 的差值。对于天然孔,在2006—2020 年之间,天然多年冻土上限由4.0 m 下降到4.7 m。对比来看,该断面的左右路肩人为多年冻土上限深度已低于原天然上限0.7~0.9 m(路基高度0.9 m),且左路肩人为多年冻土上限以下土层的年最高地层温度持续高于天然场地相同深度土层的年最高地层温度,这再一次说明高温多年冻土区内的普通路基不能保持路基下多年冻土的热稳定性。
图6 三个监测断面年最高地层温度Fig.6 Maximum ground temperature of the three observation profiles:traditional embankment(a),U-shaped crushed rock embankment(b),crushed rock revetment embankment(c)
如图6(b)所示,U 型块石路基人为多年冻土上限深度变化情况为,左路肩由2006 年的6.3 m 下降到2020 年的6.5 m,15 年间仅下降了0.2 m;右路肩由2006 年的4.5 m 抬升至2020 年的4.3 m,15 年间上限保持稳定,甚至抬升了0.2 m,右路肩的人为多年冻土上限深度高于左路肩约2 m。天然孔多年冻土上限深度由2006 年的2.5 m 下降到2011 年的2.8 m,下降速度较为缓慢,但可以发现,2012—2013 年这两年时间内,多年冻土上限由2.8 m 下降到3.9 m,结合文献[48]中的气象数据,这或与部分年份的外界气温偏高有关。对比来看,高温多年冻土区内U 型块石路基断面对路基下多年冻土的主动降温效果明显,但左右路肩下的人为多年冻土上限差异过大的问题不容忽视。
如图6(c)所示,块石护坡路基人为多年冻土上限深度变化情况为,左路肩由2006 年的8.6 m 抬升至2020 年的6.8 m,这表明块石护坡路基对多年冻土上限的抬升起到了积极作用;右路肩2006 年为5.5 m,2020 年抬升至5.2 m,值得注意的是该侧地表年最高地层温度从2006 年的10.5 ℃增长到2020年的27.5 ℃,但这17 ℃的温度增加对该侧人为多年冻土上限深度的影响并不大,同样说明了块石护坡路基对保护多年冻土的有效性。天然孔2006 年的多年冻土上限深度为3.4 m,至2020 年下降到4.4 m,这与左右路肩人为多年冻土上限深度的上升成相反的变化趋势。对比来看,15 年间块石护坡路基左路肩的人为多年冻土上限抬升显著,右路肩略有上升,但天然孔却有明显下降,同时左右路肩的人为多年冻土上限深度差值减小。以上分析表明,块石护坡路基对路基下冻土层起到了一定的主动控温作用,且缩小了路基左右侧人为多年冻土上限差值。
3 结论
本文对2006—2020 年青藏铁路运营15 年中,高温多年冻土区内的三种典型路基结构(普通路基结构、U 型块石路基结构、块石护坡路基结构)的温度场数据进行分析和讨论,得到以下主要结论:
(1)随着运营年份的增加,普通路基的人为多年冻土上限略低于原天然冻土上限,左右路肩的融深差值不明显。U型块石路基和块石护坡路基左右路肩的人为多年冻土上限均已抬升至路基本体内,左右路肩的融深差值在1.0~1.5 m。
(2)综合三个监测断面年平均地层温度的长期变化情况可以看出,高温多年冻土区内的普通路基左右路肩下冻土层的年平均地层温度增长速率高于天然场地同深度冻土层的增长速率,普通路基下的多年冻土不能保持长期热稳定性;U 型块石路基对其下的多年冻土主动降温效果明显,但左右路肩下的地温差异不可忽视;块石护坡路基对其下的多年冻土具有一定的主动降温效果,左右路肩下的地温差异要小于U型块石路基。
(3)整体来看,在全球气候变暖的背景下,青藏高原高温多年冻土区内的天然场地经历了多年冻土上限下降和不同深度地层温度上升的冻土退化过程,叠加可能出现的路基侧向水热扰动,高温多年冻土区内的青藏铁路普通路基需要采取增设块石护坡等补强措施,才能保持路基下多年冻土的热稳定性。同时,U 型块石路基及块石护坡路基虽然具有主动降温效果,但其左右路肩下的地温差异仍需引起重视,可在路基左右两侧增设不对称的热棒,以消除路基下地温场左右不均匀分布的问题。