冻结层上水的分布及工程影响研究现状与展望

　　摘要：多年冻土约占地球陆地面积的1/4，冻结层上水作为一种多年冻土区的特殊水文地质现象，其分布受局地因素的控制，且随活动层的季节性冻融而变化，影响着地表水和地下水循环以及多年冻土环境水热平衡。多年冻土将冻结层上水限制在一个狭窄的空间内，在暖季，冻结层上水横向和竖向的渗流传热将加剧多年冻土的退化，也会对上覆工程构筑物的稳定运营造成极大威胁。目前关于冻结层上水的研究主要集中在分布特征、变化规律、流量计算、渗流模拟、水热耦合等方面，以区域分布特征为出发点，对冻结层上水的研究现状进行了归纳和总结，并对其工程影响有关的渗流携热理论研究成果进行了梳理，最终对今后需要进行深入研究的方向进行了展望，这有助于全面理解冻结层上水在冻土区水文过程中的功能，为掌握冻结层上水的研究现状提供了进一步的理论参考。

　　关键词：冻结层上水;水热耦合;水文地质;水量平衡;活动层

　　0引言

　　冻结层上水是指分布于多年冻土层上、活动层中的潜水，存在于平坦的分水岭、平原和河流阶地，随着气温的变化而产生季节性冻融[1]。多年冻土层作为隔水层或弱透水层，很大程度上阻碍了冻结层上水的下渗和排泄，同时限制了地表水与冻土层下部融土层之间的水分交换和水力联系，决定了冻结层上水的分布格局和循环方式[2]。

　　冻结层上水的形成、演化、运移和动力过程在水温、渗流量等方面都具有不同于常规地区潜水的特殊性，不仅是寒区能水循环中的一个关键组成部分，而且与寒区生态环境变化密切相关，在多年冻土区的陆地水循环中有着特殊作用[3]。冻结层上水在含水层厚度、水量、水温、相态、水化学特征及水动力特征等方面均具有季节性变化特性，且受局地因素如地形、地貌、植被、降雨、岩性等的控制而产生区域性变化。

　　冻结层上水和多年冻土相互制约、相互作用，使寒区冻土水文地质条件具有特殊性：随着活动层的冻结，冻结层上水的排泄占整个径流消退的60%~80%[4-5];随着活动层的融化，冻结层上水是多年冻土与大气系统、植被系统、地表水进行水热交换的重要纽带，其空间分布控制着包气带含水量的变化特征和地表植被的分布类型[3]。全球变暖导致的冻土退化将增大冻结层上含水层厚度，势必影响到冻结层上水的分布特征和水量变化，但因其动态过程的复杂性和观测研究的诸多困难，尚缺乏对其运动规律、驱动因素和机制的系统认知[6-7]。

　　本文对冻结层上水的研究现状进行了综述，有助于系统地认知和掌握冻结层上水在多年冻土环境中的功能。从冻结层上水的赋存特征切入，分析了各类冻结层上水的分布、运移、补给以及排泄特征。以此为基础，根据全球多年冻土的分布范围，分别归纳了低纬度多年冻土区(青藏高原、祁连山)和高纬度多年冻土区(大兴安岭、俄罗斯以及北美)典型区域冻结层上水的分布特征。此外，还概括了冻结层上水对于多年冻土区既有工程构筑物稳定性的影响研究，并对冻结层上水在水循环系统中水量平衡和运移过程中渗流携热效应的理论研究和数值计算进行了总结和整理。

　　1冻结层上水的赋存特征

　　冻结层上水的定义是根据含水层与多年冻土层之间的空间关系而确定的[8-9]。1932年，俄罗斯学者将多年冻土区地下水分为冻结层上水、冻结层间水和冻结层下水，将冻结层上水又细分为冻结层上层滞水、季节性融化层水和冻结层上地下水3种[10]。第一种冻结层上水即为活动层中的水;第二种中，部分冻结层上水被季节性低温层覆盖;第三种冻结层上水为冻结层上多年融区的水。

　　随着对其认知的发展，寒区地下水的分类问题仍然没有完全解决，但对冻结层上水的形成和发展过程已有明确解释。冻结层上在冷季发生冻结，当地表开始融化后，随着融化厚度的增加，从土体内排出的水分将渐增并聚集到上限之上，便形成了一定厚度的潜流层(即为冻结层上水)[11]。冻结层上水的冻结温度在-0.5~0℃之间，其补给主要来源于大气降水和融雪融水，其排泄一部分以下降泉的方式溢出地表并汇流成溪，一部分通过沼泽湿地垂直蒸发排泄，另一部分则通过融区补给冻结层下水。

　　由于冻结层上水的产流和补给都发生在活动层，因此其补给面积和分布面积基本重合[1]。多年冻土对于冻结层上水具有阻隔、蓄水调节和蒸发抑制的作用[12-13]，目前主要通过定性研究来分析二者的相互关系。多年冻土层产生的水文效应影响到了水循环动态规律和地下水的补给排泄关系，使得降水、地表水、地下水失去了循环通道，对其三者之间的转换关系造成了制约[14]。

　　多年冻土的阻隔作用导致了土体导热导水性能的改变，直接影响着冻结层上水的下渗和蒸发[15]，多年冻土的退化将增加地表水和冻结层上水流量，其中包括冻结层上水水量、地下水和地表水的交换量、基流流量的增大。相反地，冻结层上水随着地形的变化在多年冻土层上部产生侧向流动，流动过程所引起的传热效应将造成下伏多年冻土的退化[16-17]。总而言之，冻结层上水与多年冻土在热交换、地表水-地下水动态循环及区域水文效应等方面有着密切联系。

　　2低纬度多年冻土区

　　青藏高原被称为“世界屋脊”、“亚洲水塔”、“世界第三极”，是我国重要的生态安全屏障，也是全球气候变暖的“先兆区”和“放大器”[18-19]。青藏高原地区是世界上中低纬度海拔最高、面积最大的多年冻土区，在目前气候背景下，多年冻土面积约为120×104km2，随着高原平均气温的升高，青藏高原多年冻土将发生退化，在气温年增加0.02℃和0.052℃两种不同气候背景下50a后多年冻土面积分别减少8.8%和13.4%[20-21]。

　　此外，受地质构造、地形条件、水文特征和工程影响，青藏高原地区自然灾害分布呈高强度、高频率、突发性、季节性、准周期性特征[22]。而冻结层上水作为地下水的一种，成为了多年冻土区典型的水文地质现象，在青藏高原冻土退化及灾害发生过程中有着不可忽视的影响。青藏高原冻结层上水分布范围广，含水层以第四纪松散层、山区基岩风化破碎带为主[23]。地下水的补给、流动和排泄均受到多年冻土的控制，同时，也受局地因素如地形、地貌、构造、地表水文的影响，因此，每个水文地质单元中的冻结层上水具有差异性。此外，由于多年冻土层厚度随高程的增大而加厚，这进一步改变了地下水动态变化规律，使得冻结层上水在含水层厚度、流速、水温、水文化学、动力学等方面都具有季节变化[24-26]。

　　冻结层上水含水层厚度受与冻土上限有关，在青藏高原地区一般为1~3m，尤其是河流和小湖塘附近，含水层厚度较大，储水量较高。含水层的空间分布特征明显受到局部地形的控制，地势低洼地带储水量较大，在地势较高地区受径流及蒸发的影响，因储水量低而无法形成统一的含水层。青藏高原地区含水层的形成从每年4月土层中的冰体融化开始，一直到9月份厚度达到最大[27]，冻土层上水的循环周期一般为1~4a[3]。

　　目前对于青藏高原地区冻结层上水的研究主要集中在三个区域：青藏工程走廊，黄河源区和祁连山地区。青藏工程走廊内由于集中了5条重大线性工程而备受关注，冻结层上水的存在不但对工程稳定性有着不可忽视的影响，还对走廊沿线冻土环境特征的演化过程产生着不同程度的影响[28-29];黄河源区作为黄河流域的水源地，分布有许多小盆地、径流和湖泊，丰富的地表水和地下水使其在多年冻土区具有特殊的水文地质特征[30-31]，作为流域的重要补给来源，冻结层上水在黄河源区的水循环过程中起着调控作用;祁连山地区分布有大量内陆水系，下垫面特征复杂、冰川众多、储水量较大，冻结层上水广泛发育，在水分疏导、地下水交互方面具有不可替代的作用。

　　2.1青藏工程走廊

　　与其他区域类似，青藏工程走廊内的冻结层上水埋藏和分布状况也受地形、地貌、地层、冻土上限深度等因素的制约，补给来源为大气降水、融雪、冰川融水[32]。廊内山区地带如昆仑山区、可可西里山区、风火山区、开心岭山区、唐古拉山区、头二九山区基本无冻结层上水，相反，在一些河流盆地地区如楚玛尔河断陷盆地、通天河盆地，由于表层为粉砂层、卵砾石层，地形平缓，冻土上限深，冻结层上水较为发育[33]。以沱沱河为例，冻结层上水在其流域内主要赋存于基岩类山区的活动层中，一般存在于每年5-9月，11月至次年4月全部冻结。

　　含水层厚度也受多年冻土层的控制，一般小于3 m。这类冻结层上水通过风化裂隙及构造裂隙接受大气降水的补给后，产生的径流沿着冻结面进行流动，在山脚和山坡转折处以线状、片状的方式泄于地表或顺构造裂隙带直接补给冻结层下水[34-35]。在沱沱河北部及南部的高山与高海拔丘陵地区，由于降水较少、寒冻风化较弱、节理裂隙不发育等原因，导致补给水源不足，加之所处地区储水空间有限，因此活动层内土体的富水性较差[36]。

　　在沱沱河东部的极高山地区，由于降水充沛、寒冻风化强烈、节理裂隙发育，加之山体宽厚、岩屑坡广泛分布，因此整个地区有利于大气降水入渗补给，富水性较好。除冻结层上水的空间分布特征之外，还有对其时间变化特征的研究，通过研究发现，风火山地区典型高寒草甸坡面冻结层上水位的季节动态变化与活动层土壤温度和水分变化类似，活动层土壤温度控制了冻结层上水季节性动态变化格局。研究表明，距地表60cm深度以下土壤含水率和地下水补给来源决定了冻结层上水水位动态变化的位相和幅度[7]。

　　在青藏工程走廊内，热融湖塘作为富冰冻土区的一种常见地貌形态，其演化过程对于冻土退化、地表水和地下水动态平衡均具有重要影响，湖塘的存在及积水范围的扩大成为了工程安全运营的潜在隐患，因此湖塘的演化特征在多年冻土区备受关注[37-38]，而冻结层上水在湖塘周边广泛发育，其排泄和补给过程与湖塘的变化密切相关。青藏工程走廊内青藏公路两侧从西大滩到可可西里山之间分布着数百个热融湖塘，且随着多年冻土的退化，湖塘数量呈增加趋势[39-40]。

　　冻结层上水的流入和流出影响着热融湖塘的水文过程，为了评估冻结层上水的热效应对于热融湖塘的发展的影响，可以通过水位、水温、冻结层上水流量的观测，结合探地雷达、电阻率层析成像、钻孔温度监测等方式对活动层厚度的时空变化以及湖塘周围冻土分布进行研究。研究发现，廊内典型湖塘周围的冻结层上水每年融化时长持续5个月，湖岸边温度较低的冻结层上水将补给湖水，在这一过程中伴随着与下伏冻土热交换的产生，将引起活动层厚度和湖塘周围冻土特征的定向差异[41]。

　　2.2黄河源区

　　黄河源区位于青藏高原东北部、青海省玛多县境内，承担着35%以上的黄河流域径流量，是整个黄河流域重要的水源地，极易受到人为和自然因素的影响而发生冻土环境退化[42-43]。

　　黄河源区冻融和风化作用强烈，地表岩石发生破碎，产生大量的裂隙，这有利于冻结层上水的储存。源区内地下水水位近几十年来呈现明显下降的趋势，甚至有些区域出现了冻结层上水消失的现象。多年冻土的退化使得冷生隔水层下移，导致了包气带土壤层的含水量减少，该区域出现了植被草场退化、植物多样性减少、沼泽湿地萎缩、荒漠化加剧及黄河断流等一系列生态环境问题[30,44]。整个源区内冻结层上水的发育与冻土类型密切相关，主要发育在含冰量较高的冻土层上部的季节融化带内，埋藏深度受岩性、地层结构、坡向、植被覆盖率等诸多因素综合作用控制，埋藏深度一般为0.5~2m[45]。

　　冻结层上水占水资源总量的88.3%，受冻结层上水供给的洼地和沼泽面积在逐年缩小。夏季和秋季入流量大，而春季和冬季由于活动层的冻结作用而入流量明显降低[24]。黄河源区的冻结层上水主要分布在地形相对低洼的地表水和地下水汇集区，如巴颜喀拉山主峰-尕拉拉错周围的大型古冰川雪粒盆、皱玛曲源头的冰蚀洼地、以及冲积扇前缘及沟谷中[46]。

　　黄河源区的水文地质结构、地下水位的变化以及过度放牧和开采沙金等人为作用均会对冻结层上水的空间分布产生影响，这些作用不仅影响着冻土类型的退化，还会引发一系列生态问题[46]。源区内冻结层上水的主要补给来源是大气降水和地表水。在降水期内，大气降水和地表水给予地下水充足的水量补给，同时也伴随着水热交换，从而引起冻结层上界的消融和退化。冻结层上水主要通过河流、泉和蒸发排泄，有时通过融区补给冻结层下水[45]。从沼泽草甸排泄出的冻结层上水汇聚成溪流，成为源区主要河流的发源地，是河流的主要补给来源之一。

　　2.3祁连山

　　祁连山地区多年冻土广泛发育，且温度较高，是我国第二大内陆河黑河的上游地区[47]。此地区冻地下水的补给来源主要是降水，在海拔3800m以上的多年冻土区，地下水以冻结层上水和冻结层下水形式存在;在海拔3800m以下的非多年冻土区，地下水主要以潜水形式存在[48]。基岩山区冻结层上水储存于风化裂隙或孔隙中，含水层厚度随岩性的差异而不同，厚度一般为1~2m，山区降水充沛，补给条件良好，冻结层上水较发育且水量丰富。

　　各盆地内冻土平原冻结层上水则储存于第四系冰碛-冰水相含泥砂砾石层中，含水层厚度0.5~2m[49]。冻结层上水与地表水、其他类型的地下水之间具有水力联系。冻结层上水成为了冰川和积雪融水从高山往平原地区迁移的通道，作为蓄水层，冻结层上水对于水分的疏导能力产生的季节性变化主要体现在两方面：在暖季沿着地势从高高低流向产生溪流的沟谷地带;在冷季流向因河冰和季节性冻结造成的临时承压含水层，这也说明冷季平原地区储存的地下水和河槽中的基流进行了缓慢释放[50]。

　　为了研究冻土和地下水的交互影响，通过对祁连山葫芦沟流域三维水文地质模型的计算结果表明，此流域地下水在基流的形成中占43%，地表温度升高会使多年冻土的分布范围减小，增加地下水向河流的排泄。整个流域内冻结层上水距有明显的流动性特征，在渗流过程中有少部分向下穿透了多年冻土层。冻结层下水水头较高，而冻结层上水水头较低，向上的水力梯度驱动着温度更高的地下水流向地表，促进暗渠和融区的形成，增加地下水与地表水的交互作用[51]。

　　3高纬度多年冻土区

　　3.1大兴安岭

　　大兴安岭位于欧亚大陆多年冻土区的南部边缘，海拔1100~1500m，冻结层上水在山间沟谷和河谷谷地最为发育[52]。区域内植被茂密、森林众多、降雨量大，水系有黑龙江和嫩江及其支流，如额木尔河、盘古河、呼玛河、讷谟尔河、乌裕尔河等，湿地沼泽较为发育，地下水的赋存特征受地质构造、岩性、风化作用和地形地貌等因素的影响。冻结层上水主要分布于大兴安岭北部较大的河谷平原区，含水层厚度、富水性随季节融化深度的增加而增大，此地区冻结层上水主要接受大气降水入渗补给，丰水期接受地表水补给，枯水期排泄于地表水[53-54]。此区域冻结层上含水层的分布、厚度、演变程度既与地形地貌有关，又受植被覆盖度的影响，还与地表粗糙度等地形因子有关。

　　植物根系及堆积的枯枝败叶对降水起到良好的输导入渗和涵养功能，完整基岩就构成了稳定的隔水底板，地下水在基岩风化带里形成了冻结层上水。冻结层上水分布空间狭窄，富水条件普遍不佳，再加上接受降水补给时间短，水量贫乏，动态波动特征明显[55]。基岩风化带冻结层上水按照流量大小分为水量中等和水量贫乏两类：中等水量的冻结层上水分布于伊勒呼里山分水岭北部大片中低山区，含水层岩性较为复杂，主要有侏罗系火山角砾岩;水量贫乏的冻结层上水主要分布于伊勒呼里山分水岭南部中山区以及分水岭北部零星山区，含水层岩性主要为侏罗系火山岩[56]。

　　3.2俄罗斯

　　俄罗斯多年冻土面积占国土面积的一半以上，是北半球多年冻土区的重要组成部分[57]。俄罗斯典型地区的冻结层上水动态特征与气温变化密切相关，且含水层具有冻结层上水和冻结层间水的特点。在东西伯利亚勒拿河盆地，冻结层上水分布广泛，占多年冻结层面积的40%~60%，冻结层上水和冻结层间水混合体的滞留时间为1~55a，多年冻土退化将增大地下水的补给[58]。

　　在季曼-伯朝拉地区，地下水水位动态变化与气温之间联系密切。该地区冻结层上水水位一年最小值出现在冬季末期，一年之中，近地面大气平均昼夜温度超过0℃后，包气带水分交换过程和冻结层上水水位开始出现季节性变化[59]。俄罗斯科学院西伯利亚分院冻土研究所和俄罗斯地质部的东北水文地质分部共同研究了不同地貌和岩性条件下冻结层上水的季节性变化，总结了三种典型地区的水位变化特征。

　　当冻结层上水水位距地表较浅时，依靠积雪融化和大气降水进行补给，水位在衰减时期也能上升;当水位距离地表较深时，基本不受大气降水的影响，因此水位在夏季呈总体下降的趋势;当水位与季节性融化层水资源枯竭无关，与含水层自上而下和自下而上的冻结有关[60]。除此之外，一些特殊地区如河床下融区产生的冻结层上水不受季节性冻融过程的直接影响，影响程度取决于河谷的地貌、冻土特征、水文地质条件等[61]。高纬度多年冻土区人类活动较多，除了道路等线性工程之外，还分布有市政建筑等众多基础设施，城市的发展和扩张也对冻结层上水产生影响。

　　通过已有文献可知，俄罗斯高纬度地区人类活动对冻结层上水补给的影响主要包括耕作区部分灌溉入渗和来自城市管道、蓄水系统水的渗漏两方面，通过对俄罗斯北部城市雅库茨克市区开展的综合性水文监测工作的研究表明，虽然人类活动因素产生的影响呈局部性特征，但却使冻结层上水的水文情势发生了本质性变化[60]。

　　3.3北美

　　北美地区的多年冻土主要分布于阿拉斯加和加拿大北部，对于冻结层上水的研究主要集中在几条主要河流的流域盆地如育空河、马更些河以及环北极地区的康沃利斯岛等区域[62]。

　　育空河全长约3300km，穿越了加拿大及阿拉斯加地区，流域内多年冻土的分布和厚度差异性较大。随着多年冻土的融化，近30a以来地下水占河流总流量的比例正以每年0.7%~0.9%的速度持续增大[63]。马更些河是加拿大第一长河，马更些河盆地占整个加拿大国土面积的近四分之一，盆地中有75%的区域分布有连续或不连续多年冻土[64]，在春季融化过程中，冻结层上水占整个排泄量的90%，这表明在此盆地中，冻结层上水是极地水文的重要组成部分[65]。

　　对于育空河盆地地区冻结层上水的研究主要集中在活动层冻融过程、气候-湖泊-地下水之间关系、侧向流动特征三个方面，研究方法以实地观测和模拟计算为主，将二者结合，能够透过现象揭示机理，为查明冻结层上水的变化特征奠定基础。育空河盆地多年冻土区地下水向河流的排泄量逐年增大，地下水总量将有增大的趋势，这首先表现为地下水在低洼地带的排泄范围增大，除此之外还发现，此地区冻结层上水对于基流的贡献率有所降低[66]。

　　在气候变化背景下，湖泊与地下水之间在融区形成过程中有着密切关系，气温变暖会加剧冻土融化、增加冻结层上水流量。在寒区，除了气候是冻土环境变化的主导因素外，湖面面积和地下水流动在多年冻土融化和融区形成过程中有着重要影响。湖边缘季节性的冻融循环影响着冻结层上水的运移和冻结过程，在湖面产流区域，多年冻土的分布影响着湖中水体在冻结层上下之间的交换，与此同时，冻土分布特征也受湖水及冻结层的影响[67]。

　　冻结层上水在坡面地区将产生侧向流动，与垂向流动一并形成了其特殊的水文地质现象，狼溪(WolfCreek)盆地坡面垂向的渗透和蒸发过程在南坡较为明显，而北坡的下伏多年冻土阻止了深层渗漏，使得冻结层上水得以存储[68]。加拿大的多年冻土主要分布在北部的极地地区，对于冻结层上水的研究主要集中在坡面、河流流域以及极地岛屿等典型区域。在极地地区，产生坡面地表径流的原因在于冻结层上水水位的上升和活动层土壤的饱和，坡面径流取决于冻结层上水水位，并与地表形态有关。

　　从径流源头来看，融雪所占比例最大，降雨次之，排泄方式依次为蒸发、侧向径流。侧向径流包括地表流和地下流，地下流则跟导水系数、活动层厚度、水位等因素有关[69]。在加拿大北部的砾石地貌中，有湿地、积水和径流，所有的这些地貌都与冻结层上水的局地渗流有关。

　　砾石具有很高的单位产水量和导水系数，能够促进水分流动，但是，浅层冻结锋面影响着深层的渗透，因此，当有融雪和降雨的充足补给时，在地表与水位相交的地方便有渗流，且具有季节性变化特征。通过水量平衡计算可知，上坡砾石区域的渗流是湿地和积水的主要来源，连续性的渗流改变了植被特征和土壤含水率[70]。此外，冻结层上水对河流流量有一定程度的贡献，在河道下方的活动层厚度通常更大，甚至出现融区，原因在于冻结层上水的携热效应[71]。

　　4冻结层上水对工程稳定性的影响

　　由上述可知，多年冻土对于冻结层上水的影响特征主要有以下四个方面：①减缓水循环过程;②阻隔水分运动的垂向通道，影响地表水和地下水的转换;③使含水层的变化具有季节性特点和动态规律;④冻胀作用改变了土壤结构并增强了入渗能力[14]。随着多年冻土区人类活动的加剧和各类工程的修筑，冻结层上水的存在及其渗流携热效应严重影响着各类工程基础的热收支特征，直接对工程稳定性产生威胁，冻结层上水渗流、冻土退化及工程稳定性之间的相互影响，但三者的定量关系有待研究。从已有文献来看，冻结层上水对于工程稳定性的影响研究主要集中在公路、铁路以及机场。对于多年冻土区的道路工程而言，路基修筑后冻土上限上升，将有一部分冻结层上水在上限抬升处进入并横穿路基。

　　随着运营时间的增加，活动层和路基中冻结层上水(2~4℃)的渗流将造成多年冻土温度升高，造成多年冻土退化，同时还会加速融区的出现，加剧路基不均匀沉降，导致路基病害的出现[76]。通过室内模型试验，结合现场监测资料，能够从冻结层上水夏季渗流导致路基下伏多年冻土融化的角度，进一步证明路面横向裂缝、纵向裂缝以及路基不均匀沉降的成因[77]。针对冻结层上水对公路路基的影响，可采用防水材料，阻隔水分迁移和渗透通道，对路基下部过水进行隔水治理可有效降低路基沉降，同时降低因为活动层集水而造成的冻胀[78]。

　　5理论研究

　　对于冻结层上水的理论研究主要可以归纳为三个方面：一是水文模型的建立与改进;二是水量平衡计算，三是渗流携热效应。含水层土体的导水系数、液态水储存能力、冻结层上水水位随着土体的冻融与温度密切相关，而且活动层的双向冻结控制和改变了冻结层上水的渗流过程，向下冻结控制着地表径流的产生和冻结层上水的补充，向上冻结影响着地下水尤其是冻结层上水的排出[87]。由此可见，冻结层上水的渗流过程和水位变化特征受局地因素的控制，极为复杂，因此截至目前，仍然没有有效的数学模型和理论来预测多年冻土区地下水的变化特征[88]。

　　5.1水文模型

　　国际上第一代水文模型很好地描绘了流域产流的基本规律，如斯坦福模型、水箱模型、新安江模型，但这类模型忽略了地形、植被等的空间变化对流域产流过程的影响。随后发展出了一批分布式水文模型，如TOPMODEL、SHE、VIC、DHSVM、GBHM等[89]。近十多年来，由于气象学、生态学、冰冻圈等学科的兴起和发展，分布式水文模型与多学科的交叉研究越来越多，这使得建模难度增加，分布式水文模型开始描述不同圈层之间的关系，即“多圈层水文模型”[90-91]。目前中国寒区水文模型的主要限制因素是缺乏全面、同步的观测数据，可通观测、试验、模拟的方法，获取寒区水文过程的基本规律、经验公式和关键参数[92]。基于物理学基础的水文数值模型能够定性和定量地分析含水层冻融过程中的关键影响因素，且能够反映出水文要素的变化。

　　在已有水文模型的基础上，可以采用热传导方程模拟冻土上限在水文模型中的演变过程[93]，结合能量平衡计算地表热状况的空间分布特征，以此为基础可以得到土壤剖面的融化深度[94]，也可以通过水分迁移模型和热传导模型的耦合，可以研究多年冻土形成和演化过程中的水力特性，同时研究地下水和地表水的季节性变化特征[95]，还可以采用达西方程、流体连续性方程、溶质连续性方程、能量迁移方程来计算冰川作用下地下水流动[96]。此外，也可以通过在已有的简易地下水模型(SIMGM)、陆面模式(GLM)中增加未冻水含量、相变、土层性质等参数等方式进行改进[97]，建立面向对象的模型系统，来模拟和预测中小尺度的区域中水循环特征[98]及对气候的反馈[99-100]。

　　6小结与展望

　　冻结层上水的研究正从定性的特征分析向定量的机理研究发展，已有的观测和模拟研究已经揭示了冻结层上水的影响因素、赋存特征、补给排泄以及动态变化。冻结层上水的存在使多年冻土区的水文地质条件变得更为复杂，冻结层上水与多年冻土层的相互影响、相互制约的交互作用将会对整个冻土区的生态环境产生深远影响。在今后有必要在以下几个方面进行深入研究：

　　(1)对局地因素条件下冻结层上水的水文要素时空变化特征进行研究。如在坡面地区，需要继续深入研究冻结层上水对于坡面径流的调控机理，可针对典型地区不同的水文地理和水文生态单元开展比较研究，用以获得春季融化和秋季冻结期间冻结层上水对于坡面径流的补给占比以及年度变化特征。

　　(2)对多年冻土区浅表层地下水与地表水的相互关系及转化特征进行研究。在地势缓和且冻结层上水发育地区，往往伴随着一定的地表积水，应该对地表积水入渗至活动层以及冻结层上水在冻土退化过程中迁移至地表的机理、变化特征进行深入研究，揭示地表积水与冻结层上水之间的转换规律。

　　(3)对多年冻土区路基稳定性受冻结层上水的影响机理进行研究。多年冻土区路基/地基下部冻结层上水的存在成为其稳定性的潜在隐患，有必要对路基沉降变化特征与冻结层上水之间的关系进行研究，并对冻结层上水的渗流携热效应对路基变形的影响程度进行定量研究，以此来修正有冻结层上水存在条件下的路基变形计算模型。

　　(4)对气候变暖背景下冻结层上水与多年冻土交互作用及变化特征进行预测。需要对多年冻土与冻结层上水的相互作用进行研究，揭示活动层的冻融过程对多年冻土退化的影响，以及多年冻土对冻结层上水流动特征的限制，并结合气候变化模式预测二者的动态发展规律，获得更具有普适意义的冻结层上水动态规律和驱动机制。

　　(5)对冻结层上水水位、流速、构成等基本参数在测试过程中的新方法和新手段进行探索。目前对于冻结层上水的监测主要集中在水位和水温方面，受客观条件制约，冻结层上水的流量和流速很难被精确测量，针对这一问题，应采用新的监测方法和传感器，提高测量精度，从整体上掌握冻结层上水的存在范围和流动特征，为理论研究提供更为精确的参数。

　　参考文献(References):

　　[1]DingYongjian,ZhangShiqiang,ChenRensheng,etal.Introductiontohydrologyincoldregions[M].Beijing:SciencePress,2017.[丁永建,张世强,陈仁升,等.寒区水文导论[M].北京:科学出版社,2017.]

　　[2]MetcalfeRA,ButtleJM.SoilpartitioningandsurfacestorecontrolsonspringrunofffromaborealforestpeatlandbasininnorthcentralManitoba,Canada[J].HydrologicalProcesses,2001,15(12):2305-2324.

　　[3]ChengGuodong,JinHuijun.GroundwaterinthepermafrostregionsontheQinghai-TibetPlateauanditchanges[J].HydrologyandEngineeringGeology,2013,40(1):1-11.[程国栋,金会军.青藏高原多年冻土区地下水及其变化[J].水文地质工程地质,2013,40(1):1-11.]

　　[4]LiuGuangsheng,WangGenxu,SunXiangyang,etal.VariationcharacteristicsofstableisotopesinprecipitationandriverwaterinFenghuoshanpermafrostwatershed[J].AdvancesinWaterScience,2012,23(5):621-627.[刘光生,王根绪,孙向阳,等.多年冻土区风火山流域降水河水稳定同位素特征分析[J].水科学进展,2012,23(5):621-627.]

　　[5]WangGX,HuHC,LiTB.Theinfluenceoffreeze-thawcyclesofactivesoillayeronsurfacerunoffinapermafrostwatershed[J].JournalofHydrology,2009,375(3/4):438-449.

　　作者：柴明堂，马巍，穆彦虎

转载请注明来自发表学术论文网：http://www.fbxslw.com/nylw/28701.html