一种考虑降雨空间分布特征的水库流域子流域划分方法
【技术领域】
[0001] 本发明涉及水文分析技术领域，具体涉及一种考虑降雨空间分布特征的水库流域 子流域划分方法。
【背景技术】
[0002] 近年来，受极端气候的影响，水库流域暴雨频发，当暴雨发生时，若没有准确可靠 的暴雨水情预测信息作为指导，水库将无法及时调整运行方式，容易造成水库发生大量弃 水，严重时甚至会影响大坝安全及水库下游防洪安全，因此，准确可靠的水库流域水雨情预 测显得尤为重要。由于流域降雨存在空间分布不均性，为提高水库流域水雨情预测准确性， 一种可行的方法是将整个流域划分为若干个小的子流域。现有的子流域划分方法主要为： 首先基于流域数字高程数据分析流域范围，然后以流域内建设的关键流量观测站点为控制 节点，划分得到各子流域。但是，该类方法没有考虑水库流域降雨空间分布特征，当流域内 流量观测站点设置较多时，子流域划分的过于细致，子流域划分的数量庞大，增大了水雨情 预测的计算量；而当流域内流量观测站点设置较少时，子流域划分的过于粗略，影响了水雨 情预测的准确性。
【发明内容】
[0003] 本发明要解决的技术问题是克服现有技术的不足，提供一种既能充分反应流域降 雨分布不均性、又能降低水雨情预测计算量的考虑降雨空间分布特征的水库流域子流域划 分方法。
[0004] 为解决上述技术问题，本发明提出的技术方案为：
[0005] -种考虑降雨空间分布特征的水库流域子流域划分方法，包括以下步骤：
[0006] (1)根据数字高程数据，采用常规方法，获取待划分的水库流域所控制的流域范 围；
[0007] (2)以所述流域范围内的流量控制站为节点，根据数字高程数据，获取各节点所控 制的子流域范围；
[0008] (3)根据所述子流域范围内的雨量站的观测数据，计算各子流域内雨量站间观测 数据的差异程度，计算公式如下：
[0009]
[0010] 式中，η为子流域内雨量站观测数据的时段数;m为子流域内雨量站个数；Xl]为子 流域第j个雨量站第i个时段的观测值3为第i个时段子流域内所有雨量站观测值的平 均；Z为子流域内雨量站间观测数据的差异程度；Z越大表示子流域内雨量站的差异程度越 大，当子流域内所有雨量站观测数据完全相同时Z的值为0 ;
[0011] (4)设子流域内雨量站间观测数据的差异程度阈值为Zt，将步骤（3)计算得到的 各子流域内雨量站间观测数据的差异程度Z与Z t进行比较，Z值小于阈值Z τ的子流域不需 再细分；Z值大于或等于阈值Zt的子流域需要进一步细化，细化的方式为根据该子流域的 水系发育特征采用常规方法将该子流域内差异程度大的雨量站划分至不同的子流域，直到 细化的子流域内雨量站差异程度Z小于阈值Z t;
[0012] (5)当各子流域均不需再分后，计算各两两相邻子流域雨量站观测数据的差异程 度S，计算公式如下：
[0013]
[0014] 式中，η为两子流域雨量观测数据的时段数；1；为子流域X第i个时段流域内所有 雨量站观测值的平均；％为与X相邻的子流域y第i个时段流域内所有雨量站观测值的平 均；S为两相邻子流域雨量站观测数据的差异程度；
[0015] (6)设两两相邻子流域雨量站观测数据的差异程度阈值为St，将步骤（5)计算得 到的各两两相邻子流域雨量站观测数据的差异程度S与S t进行比较，若两相邻子流域雨量 站观测数据的差异程度S小于阈值ST，则此两相邻子流域合并为一个；
[0016] (7)重复步骤（5)~（6)，直至各子流域均不可合并。
[0017] 上述的考虑降雨空间分布特征的水库流域子流域划分方法，优选的，取值范 围为1~5。
[0018] 上述的考虑降雨空间分布特征的水库流域子流域划分方法，优选的，取值范 围为2~5。
[0019] 与现有技术相比，本发明的优点在于：
[0020] 本发明的考虑降雨空间分布特征的水库流域子流域划分方法，充分考虑流域降雨 空间分布的不均匀性，同时将相似的子流域进行合并，尽可能减少流域水雨情的计算量，为 水库流域子流域划分提供了一种科学合理的方法。该方法思路清晰，操作方便，实用性强。
【附图说明】
[0021] 图1为本发明实施例水库流域所控制的流域范围示意图。
[0022] 图2为本发明实施例的水库流域中子流域所控制的流域范围示意图。
[0023] 图3为本发明实施例的水库流域中子流域经细化后所控制的流域范围示意图。
[0024] 图4为本发明实施例的水库流域中子流域经合并后所控制的流域范围示意图。
【具体实施方式】
[0025] 以下结合说明书附图和具体实施例对本发明作进一步描述。
[0026] -种本发明的考虑降雨空间分布特征的水库流域子流域划分方法，包括以下步 骤：
[0027] (1)根据数字高程数据，采用常规方法，通过数字高程数据填洼、栅格流向计算、栅 格累积汇流量计算、河网提取、流域范围提取等步骤获得本实施例的某水库流域所控制的 流域范围，如图1所示。
[0028] (2)以该水库流域范围内设置的4个流量控制站为节点，采用同步骤⑴的方法获 得各节点所控制的子流域范围，即将该流域划分为4个子流域（A、B、C和D)，如图2所示。
[0029] (3)根据各子流域内雨量站观测数据，计算4个子流域内雨量站间观测数据的差 异程度，计算公式如下：
[0030]
〇>
[0031] 式（1)中，η为子流域内雨量站观测数据的时段数；m为子流域内雨量站个数；Xl] 为子流域第j个雨量站第i个时段的观测值；秀为第i个时段子流域内所有雨量站观测值 的平均；Z为子流域内雨量站间观测数据的差异程度；Z越大表示子流域内雨量站的差异程 度越大，当子流域内所有雨量站观测数据完全相同时Z的值为0 ;
[0032] 按照式（1)计算（计算过程省略）得出本实施例计算结果为：子流域A内雨量站 间观测数据的差异程度Z = 1. 5 ;子流域B内雨量站间观测数据的差异程度Z = 0. 5 ;子流 域C内雨量站间观测数据的差异程度Z = 1 ;子流域D内雨量站间观测数据的差异程度Z = 6〇
[0033] (4)设子流域内雨量站间观测数据的差异程度阈值Zt，将步骤（3)计算得到的各 子流域