面向分布式水文模型的季节性冻土区冻土层厚度模拟方法

技术领域

本发明涉及季节性冻土技术领域，具体涉及一种面向分布式水文模型的季节性冻土区冻土层厚度模拟方法。

背景技术

分布式水文模型用于季节性冻土区水文循环模拟时，往往面临着土壤冻融的影响，与暖热地区相比，寒区降雪量的比例更高，土壤冻结的时间更长，气温回升后融雪水因冻土的存在不能及时下渗补给地下水，在冻土区以上形成高含水率土壤层，甚至形成饱和土壤层，造成规模较大的春汛。水文模型如果不能较好地刻画冻土层的冻融过程，将影响地下水补给、土壤水含量、春汛径流的模拟精度。

季节性冻土的形成、发展及消融受多种因素的影响，包括气温、土温、土壤类型、土壤含水率、地表覆盖度等，在时间跨度上和空间分布上存在很大的差异。目前的分布式水文模型对冻土的动态一般都是简化处理甚至不予考虑，这在低纬度地区的水文模拟中尚能接收，但是在气候寒冷的高海拔、高纬度地区，则难以适用。

发明内容

针对现有技术中的上述不足，本发明提供的面向分布式水文模型的季节性冻土区冻土层厚度模拟方法解决了现有的分布式水文模型在气候寒冷的季节性冻土区的适用性和模拟精度都不高的问题。

为了达到上述发明目的，本发明采用的技术方案为：面向分布式水文模型的季节性冻土区冻土层厚度模拟方法，包括以下步骤：

S1、建立季节性冻土区冻土层中土壤温度与其影响因素之间的关系式；

S2、基于土壤温度与其影响因素之间的关系式，建立冻土层厚度与土壤温度之间的关系，进而实现季节性冻土区冻土层厚度模拟。

进一步地，所述步骤S1中，以日为尺度，建立的土壤温度与其影响因素之间的关系式为：

式中，T

进一步地，所述影响因子f

式中，z

其中，土壤温度稳定深度d

式中，

d

式中，ρ

进一步地，所述土表温度T

T

式中，f

其中，地表覆盖物影响土表温度的权重因子f

f

式中，f

其中，f

f

式中，veg为植被生物量和残余量，单位为kg/ha；sno为积雪当量，单位为mm；

裸土土表温度T

式中，

式中，H

进一步地，所述步骤S2具体为：

S21、将土壤剖面分为M层，基于土壤温度与其影响因素之间的关系式，通过线性插值法计算冻土层的上边界和下边界；

S22、基于冻土层的上边界和下边界，计算冻土层的厚度，进而实现季节性冻土区冻土层厚度模拟。

进一步地，所述步骤S21中，冻土层的上边界z

式中，z

冻土层的下边界z

所述步骤S22中，冻土层的厚度Δz为：

Δz＝z

本发明的有益效果为：

本发明方法通过先建立土壤温度与相关影响因素之间的关系，然后再建立冻土厚度与土壤温度之间的关系，提高了分布式水文模型在气候寒冷的季节性冻土区的适用性，同时，利用冻土区厚度监测数据率定模型中的参数变量，进一步提高模拟精度。

附图说明

图1为本发明提供的面向分布式水文模型的季节性冻土区冻土层厚度模拟方法流程图。

图2为本发明提供的季节性冻土区冻土发展过程示意图。

具体实施方式

下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。

本发明实施例中，冻土的冻结与融化直接受到土壤温度的影响，而土壤温度与气温、地表覆盖物、土壤含水率等相关，因此，如图1所示，本发明中的面向分布式水文模型的季节性冻土区冻土层厚度模拟方法，包括以下步骤：

S1、建立季节性冻土区冻土层中土壤温度与其影响因素之间的关系式；

S2、基于土壤温度与其影响因素之间的关系式，建立冻土层厚度与土壤温度之间的关系，进而实现季节性冻土区冻土层厚度模拟。

上述步骤S1中，根据水文模型开展日尺度水循环模拟的需求，以日尺度计算土壤日平均温度和冻土日平均厚度(以下无特殊说明，均以土壤温度和冻土厚度简称)，某一年第n日、深度z处的土壤温度与年平均气温、土表温度、前一天的土壤温度有关，基于此，以日为尺度，建立的土壤温度与其影响因素之间的关系式为：

式中，T

其中，初始土壤温度T

具体地，影响因子f

式中，z

从上式(1)～(3)中可以看出，土壤深度越接近地表，z

式中，

d

式中，ρ

具体地，上述土表温度T

T

式中，f

地面覆盖物一般分为两类，一是植被及其残余，二是积雪；前者不断影响土壤温度，而且影响蒸散发，因此分布式水文模型一般都具有模拟地表植被生长、成熟、休眠、收割等过程的功能，能够计算地表植被生物量及其脱落地表的残余量。模拟降雪-积雪-融雪过程是分布式水文模型的基本功能之一。基于此，地表覆盖物影响土表温度的权重因子f

f

式中，f

其中，f

f

式中，veg为植被生物量和残余量，可在水文模型中逐日计算，单位为kg/ha；sno为积雪当量，可在水文模型中逐日模拟，单位为mm；

裸土的土表温度T

式中，

式中，H

在季节性冻土区，随着气温降低，表层土壤首先冻结，然后冻土层逐日增厚，气温回升到0度以上，表层土壤开始融化，直至冻土层融通(如图2所示)；

一般而言，土壤温度沿深度垂向连续变化，这样模型可以将土壤剖面分为若干层，分别计算每层土壤中心处的土壤温度，其他深度处的温度可以插值计算，以减少模型运算的负担。如果相邻两个土壤层中心处的温度，一个大于0℃，一个小于0℃，则通过插值可以确定温度为0℃的深度，由此确定冻土层边界的深度。在土壤剖面分层方面，可等厚度平分，也可在水文模型的土壤分层基础上进一步细分，温度变化比较剧烈的表层土壤可以细分，总之，用户可根据需求灵活分层。因此，上述步骤S2具体为：

S21、将土壤剖面分为M层，基于土壤温度与其影响因素之间的关系式，通过线性插值法计算冻土层的上边界和下边界；

S22、基于冻土层的上边界和下边界，计算冻土层的厚度，进而实现季节性冻土区冻土层厚度模拟。

上述步骤S21中，冻土层的上边界z

式中，z

冻土层的下边界z

步骤S22中，冻土层的厚度Δz为：

Δz＝z

季节性冻土区冻土层在垂向上一般分为单层，局部气候异常可能导致局部区域、个别时段出现多层，可以采用上述方法分别计算各冻土层的厚度；本实施例中的季节性冻土区冻土层厚度模拟需要嵌入到分布式水文模型中，需要用户对部分参数变量赋初值，还可以利用冻土厚度监测值率定模型中的这些参数变量，提高模拟精度。