一种流域分布式水文模型子流域内高程带的划分方法与流程

本发明涉及水利工程领域，具体涉及一种流域分布式水文模型子流域内高程带的划分方法。

背景技术：

在流域分布式水文模型的应用过程中，能够实现气象数据在流域空间上的非均匀性处理是其重要的功能特点。气象数据包括降雨、气温、风速、日照时速、空气湿度等要素，气象数据一般来源于流域上的各个气象站点，在空间上一般是点状分布的，各个气象站点都有其空间位置，是相互独立的，各个气象站点上观测的气象数据视其所处环境，相差可能十分明显。流域水文模型在使用站点的气象数据时不是直接使用的，而是要将各个站点的气象数据通过数据处理展布为流域面上的气象数据再使用，即数据的点-面转化过程。

流域分布式水文模型对气象数据的空间展布方法有很多，最常见的是类似泰森多边形的处理方法，即首先确定模拟范围内各个子流域的形心坐标位置，再根据子流域形心坐标位置与各个气象站点坐标位置的关系，选择与子流域形心坐标位置最近的气象站点作为该子流域面上的气象数据。气象数据的空间展布效果与气象站点的数量、空间位置、模拟范围内的子流域的个数等密切相关。一般气象站点的数量越多，空间分布上越均匀，模拟范围内的子流域划分得越多，气象数据的空间展布效果越好，越能接近气象数据在整个模拟范围上的面上真实情况。

气象数据有沿海拔高程变化的特征，如降水量随高程的变化早已为人们所注意，由于水汽源的距离不同，水汽含量及输送强度不同，以及山脉的走向、高度、坡度及局部地形等因素的复杂组合，距离不远的气象站点，海拔不同时观测到的降水量有时相差很大。

气象数据是流域分布式水文模型最为基础的输入数据，若气象数据不准确，将对模型模拟结果的正确性、可靠性产生直接的影响。在模拟范围内的地表高程起伏剧烈，高程差大，气象站点个数有限的情况下，必须根据各个子流域内的高程特征对其气象数据进行修正处理，以反映气象数据沿海拔高程变化的规律，这对于山区流域的水文模拟尤为重要，因为山区地形变化的远比平原区要大。

对于一个整体流域，划分为多个子流域时，受地形起伏变化和分水岭影响，各个子流域的面积大小，以及子流域内部高程变化有可能相差很大；如果每个子流域都按照固定的高程带个数进行划分，则高程带的带宽在空间上的分布将十分不均匀，特别是一些面积较小、高程变化小的子流域；划分较多的高程带对于提高气象数据的使用精度和水文模拟精度的效果基本没有增益，反而会影响分布式水文模拟的运行效率。

因此，目前缺少一种能够根据子流域内高程分布情况自动划分高程带的方法。

技术实现要素：

本发明针对现有技术的上述不足，提供了一种能够解决现有水文模型子流域划分方法无法充分考虑地形变化等问题的流域分布式水文模型子流域内高程带的划分方法。

为解决上述技术问题，本发明采用了下列技术方案：

提供了一种流域分布式水文模型子流域内高程带的划分方法，其包括如下步骤：

s1、获取dem高程数据，计算水文模型子流域的总高程差d总；

将得到的总高程差d总与设定的基准带宽值bnd参考比较，计算子流域的高程带数量nbnd和高程带的实际带宽值bnd实际，其中高程带数量nbnd的取值范围为1至10的整数；

s3、根据高程带的实际带宽值bnd实际，依次计算每个高程带i的高程上限elehigh,i和下限值elelow,i；

s4、根据每个高程带i的高程上限elehigh,i和下限值elelow,i，计算相应高程带i占子流域的面积比例ri；

s5、根据每个高程带i占子流域的面积比例ri，计算相应高程带i的形心位置levi，完成高程带划分。

进一步地，步骤s1的具体计算方法包括：

s1-1、根据dem高程数据计算得到子流域内最低的高程特征值dmin,ini和最高的高程特征值dmax,ini；

s1-2、计算得到修正后的子流域内的最低高程值dmin和最高的高程特征值dmax，其计算公式为：

s1-3、子流域的总高程差d总为修正后高程特征值的差值，其计算公式为：

d总＝dmax-dmin。

进一步地，步骤s2的具体计算方法包括：

s2-1、当子流域的总高程差d总小于2倍基准带宽时，子流域内高程带数量nbnd为1，每个高程带的实际带宽bnd实际与总高程差d总相同，其计算公式为：

s2-2、当子流域的总高程差d总大于2倍基准带宽并且小于10倍基准带宽时，则子流域内高程带数量nbnd和高程带的实际带宽bnd实际的计算公式为：

s2-3、当子流域的总高程差d总大于10倍基准带宽时，则子流域内高程带数量nbnd和高程带的实际带宽bnd实际的计算公式为：

进一步地，步骤s3的具体计算方法包括：

进一步地，步骤s4的具体计算方法为：

在子流域中任取高程为1m的子流域高程带作为特征高程值高程带，逐块计算每一个特征高程值高程带j占目标高程带的面积比例rj，然后对rj进行求和，其计算公式为：

ri＝∑rj

其中，ri为划分的第i个高程带占子流域的面积比例。

进一步地，逐块计算每一个特征高程值高程带j占目标高程带的面积比例rj的方法为：

s4-1、当特征高程值高程j带完全位于划分的高程带范围内时，该特征高程值高程带所占的子流域面积完全计入划分的高程带内，其计算公式为：

rj＝frj

其中，frj为特征值高程带占目标子流域的面积比例；

s4-2、当特征高程值高程带j完全不位于划分的高程带范围内时，该特征高程值高程带所占的子流域面积不计入划分的高程带内，其计算公式为：

rj＝0；

s4-3、当特征高程值高程带j与目标高程带的下限部分重叠时，将重叠部分计入划分的高程带内，其计算公式为：

rj＝(elehigh,j-elelow,i)×frj

其中，elehigh,j为特征值高程带的高程上限；

s4-4、当特征高程值高程带j与目标高程带的上限部分重叠时，将重叠部分计入划分的高程带内，其计算公式为：

rj＝(elehigh,i-elelow,j)×frj

其中，elelow,j为特征值高程带的高程下限。

进一步地，步骤s5的具体计算方法为：

levi＝∑(levj×rj)/∑rj

其中，每个特征高程带j的形心高程levj的计算公式为：

levj＝(elehigh,j+elelow,j)/2。

进一步地，高程带的基准带宽的取值大于2米且小于1000米。

本发明提供的上述流域分布式水文模型子流域内高程带的划分方法的主要有益效果在于：

通过将将子流域分为1至10个高程带，有效控制水文模型中对气象数据的处理工作量，提高水文模型的模拟效率，防止计算量过大，且能有效反映年降水量在子流域内部的变化。

根据子流域高程的总带宽与用户设定的高程带宽值bnd参考之间的相对关系，以划分的带宽尽量接近用户设定的高程带宽值为原则，合理划分子流域内高程带的数量，并能使得各个子流域划分的高程带宽度尽量在整个流域范围内均匀化。

附图说明

图1为子流域内高程带的划分方法的步骤示意图。

图2为子流域内高程带划分示意图。

具体实施方式

下面结合附图对本发明作进一步说明：

如图2所示，其为子流域内高程带划分示意图。

本发明提供的流域分布式水文模型子流域内高程带的划分方法包括如下步骤，如图1所示：

s1、获取dem高程数据，计算水文模型子流域的总高程差d总。

首先将dem高程数据源中的不合理高程值，如dem的空白值清除，并清除不合理高程涉及的面积，以保证dem高程数据的可用性。

进一步地，计算水文模型子流域的总高程差d总的具体计算方法包括：

s1-1、根据dem高程数据计算得到子流域内最低的高程特征值dmin,ini和最高的高程特征值dmax,ini。

这两个值为通过dem高程数据可以直接读取得到。

s1-2、计算得到修正后的子流域内的最低高程值dmin和最高的高程特征值dmax，其计算公式为：

通过对子流域的高程特征值进行修正，以便于进一步计算分析。

s1-3、子流域的总高程差d总为修正后高程特征值的差值，其计算公式为：

d总＝dmax-dmin。

s2、将得到的总高程差d总与设定的基准带宽值bnd参考比较，计算子流域的高程带数量nbnd和高程带的实际带宽值bnd实际，其中高程带数量nbnd的取值范围为1至10的整数。

首先由用户设定子流域高程带的基准带宽值bnd参考，然后根据bnd参考与d总的大小关系，确定高程带数量nbnd和实际带宽值bnd实际。其具体计算方法包括：

s2-1、当子流域的总高程差d总小于2倍基准带宽时，子流域内高程带数量nbnd为1，每个高程带的实际带宽bnd实际与总高程差d总相同，其计算公式为：

例如，若用户设定的基准带宽为5米，但某子流域的总高程差为8米，那么由于子流域的总高程差小于基准带宽的2倍，即10米，因此该子流域整体划分为1个高程带，高程带的带宽为8米。

s2-2、当子流域的总高程差d总大于2倍基准带宽并且小于10倍基准带宽时，则子流域内高程带数量nbnd和高程带的实际带宽bnd实际的计算公式为：

其中，符号表示向下取整。例如，当子流域总高程差为28米，用户输入的基准带宽为5米，划分高程带时，子流域将被划分为5个高程带，每个高程带的实际带宽为5.60米。

s2-3、当子流域的总高程差d总大于10倍基准带宽时，则子流域内高程带数量nbnd和高程带的实际带宽bnd实际的计算公式为：

例如，当子流域总高程差为80米，用户输入的基准带宽为5米，划分高程带时，子流域将被划分为10个高程带，每个高程带的实际带宽为8米。

s3、根据高程带的实际带宽值bnd实际，依次计算每个高程带的高程上限elehigh,i和下限值elelow,i。

将划分的各个高程带，按照高程从下至上的顺序排列并标序，序号记为i。则对应每个高程带的高程上限elehigh,i和下限值elelow,i的计算公式为：

其中，i为高程带的序号。

s4、根据每个高程带i的高程上限elehigh,i和下限值elelow,i，计算每个高程带占子流域的面积比例ri。

其具体计算方式为在子流域中按照1m高程带带宽划分的子流域高程带作为特征高程值高程带，这里取各个特征高程值下的高程带，其带宽内的高程在空间上为均匀分布，这样可以根据各个特征高程值下的高程带高程上、下限与划分的各个高程带的上、下限之间的空间叠加关系，直接计算各个高程带占子流域的面积比例。

逐块计算每一个特征高程值高程带占目标高程带的面积比例rj，然后对rj进行求和，其计算公式为：

ri＝∑rj

其中，ri为划分的第i个高程带占子流域的面积比例。

进一步地，根据特征高程值高程带j与目标高程带的位置关系，其具体计算方法包括：

s4-1、当特征高程值高程带j完全位于划分的高程带范围内时，该特征高程值高程带所占的子流域面积完全计入划分的高程带内，其计算公式为：

rj＝frj

其中，frj为特征值高程带占目标子流域的面积比例；

s4-2、当特征高程值高程带j完全不位于划分的高程带范围内时，该特征高程值高程带所占的子流域面积不计入划分的高程带内，其计算公式为：

rj＝0；

s4-3、当特征高程值高程带j与目标高程带的下限部分重叠时，将重叠部分计入划分的高程带内，其计算公式为：

rj＝(elehigh,j-elelow,i)×frj

其中，elehigh,j为特征值高程带的高程上限；

s4-4、当特征高程值高程带j与目标高程带的上限部分重叠时，将重叠部分计入划分的高程带内，其计算公式为：

rj＝(elehigh,i-elelow,j)×frj

其中，elelow,j为特征值高程带的高程下限。

s5、计算每个高程带的形心位置levi，完成高程带划分。

各个高程带的形心位置高程可通过各个与之有空间重叠关系的特征值高程带的高程特征值levj及面积比例加权计算，其计算公式为：

levi＝∑(levj×rj)/∑rj

其中，每个特征高程带j的形心高程levj的计算公式为：

levj＝(elehigh,j+elelow,j)/2。

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