基于多源卫星降水信息和产流约束校正的无资料地区洪水实时预报方法

技术领域

本发明属于洪水实时预报技术领域，涉及一种基于多源卫星降水信息和产流约束校正的无资料地区洪水实时预报方法。

背景技术

洪水实时预报可有效减轻洪水造成的损失。由于降雨是驱动水文模型的首要信息，因此精确的降雨信息在洪水实时预报中起着至关重要的作用。目前，基于地面站观测的降雨信息被广泛应用于洪水预报。但是，我国河流众多，受投资资金、地形气候以及设备维护等条件限制，目前仍有很多地区无降雨监测站点。近年来，随着卫星遥感技术的迅猛发展，定量降雨分析技术为解决无资料地区洪水预报提供了新的思路。目前，全球多机构通过融合多源卫星传感器发布了多种卫星降水产品，且每种产品包含多个版本，因此卫星降水产品版本繁多且各有优劣。

作为一种间接估测结果，卫星遥感估测降雨产品往往存在着系统误差和不确定性，其中近实时产品的误差尤为显著。即便如此，这仅有的近实时卫星降水产品对于无降雨资料地区是尤为可贵的，因为对于洪水实时预报而言，卫星降水产品的滞时长短直接影响洪水预报的预见期和时效性。大量卫星降水产品评价的研究表明，尽管近实时产品的精度明显低于后处理产品，但它们都能一定程度上描述降雨的特征和趋势。因此，亟需发明一项针对近实时卫星降水产品的洪水实时预报方法，从而充分利用近实时降水产品的有效信息。

许多研究在提高卫星降雨产品驱动的洪水预报精度方面做出尝试，其中最为常见的方法就是以卫星降雨产品为输入重率定水文模型，但由此得到的模型参数会较大程度地偏离参数的物理意义范围，导致模型参数不具备代表性，洪水预报结果存在较大不确定性。此外，进一步提高卫星降雨产品驱动的洪水预报精度的方法主要包含三类。第一类，在驱动水文模型之前，根据少量的或流域周边的地面站观测降雨与多源卫星降雨产品的误差统计关系，校正产品的误差，再以校正后的产品驱动水文模型进行洪水预报；第二类，利用多源卫星降雨产品分别直接驱动水文模型，再基于权重或贝叶斯平均等方法融合水文模型的模拟结果；第三类，在水文模型的模拟过程中利用实时观测的流量或土壤含水量更新校正洪水预报结果。多研究对比表明，由于流量可以被直接观测，利用实时流量校正预报结果在上述三类方法中表现最好。但利用观测流量直接校正预报流量大大缩短了洪水预报的预见期，因此利用观测流量间接校正产流或降雨成为了更优选择。此外，由于校正过程中缺乏约束，往往导致校正后的洪水过程跳动现象严重，洪水预报结果不稳定。但目前仍没有任何一种方法是基于多源卫星降水产品和产流约束校正来实时预报无资料地区的洪水。

发明内容

针对现有技术存在的问题，本发明提供一种基于多源卫星降水信息和产流约束校正的无资料地区洪水实时预报方法。

本发明采用的技术方案为：

一种基于多源卫星降水信息和产流约束校正的无资料地区洪水实时预报方法，包括以下步骤：

第一步，拆分流域出口流量组成，获得本发明的校正对象。

流域尺度的水文模型可以被表达为：

Q＝f(P，E，X，θ) (1)

其中，Q＝[Q

产流量作为一个可以综合反映函数(1)中自变量的中间变量，数值上等于净雨，表达式为：

R

其中，R

Q＝Q(R，θ) (3)

其中，R＝[R

当任一时段的降雨、蒸散发和流域状态变量驱动水文模型时，可以得到：

q

其中，q

q

其中，q

式(6)的矩阵形式为：

Q＝LA (7)

其中，A＝[1，1，…，1]

第二步，采用多源卫星降水产品驱动水文模型后得到第一步获得的校正对象的多个序列，将多序列的最值形成的范围作为约束，进行产流约束校正和流量预报。

(1)约束

对于任一时段，有h种卫星降雨产品可供驱动水文模型，由此可得h种产流量，即

(2)校正

Q

Q

A′＝(ΔL

此外，α

为了约束校正产流量，A′中的元素也被约束，即当α

(3)预报

对于预报而言，不仅A′中前(m′-T)个元素很重要，后T个元素也很重要，因为未来的流量过程受前m′个时段的所有的降雨的影响。考虑到A′中的后T个元素也是流量预报的关键，本发明令其都等于α

Q

其中，

第三步，基于本校正和预报方法的参数重率定

(1)评价指标

本发明选用以下4个指标评价每场洪水的预报表现：纳什系数(NS)、洪峰相对误差(RPF)、峰现时间误差(EPT)和产流相对误差(RRD)。表达式如下：

EPT＝T

其中，NS由场次洪水计算得到；Q

(2)参数率定

θ的组成元素根据水文模型的选取而不同，并且需要优化算法进行率定。

采用NS

本发明的效果和益处为：本发明通过将多源卫星降水信息转化为产流的上下限，从而对产流进行有约束的校正，可大幅提高卫星降雨产品驱动的水文模型的洪水预报精度，以支持无降雨资料地区洪水预报，并为多源卫星降雨产品的应用提供了新思路。

附图说明

图1是本发明涉及的流量过程分解，图中字母含义与公式中涉及字母的含义相同。

图2是本发明以T＝1和m’＝3为例的校正和预报运算。

图3是基于本发明的参数率定框架。

具体实施方式

本发明在以多源卫星降水信息驱动水文模型的基础上，提出一种基于多源卫星降水信息和产流约束校正的无资料地区洪水实时预报方法。

下面通过实施例，对本发明做进一步说明。

小二沟流域位于嫩江西岸，发源于内蒙古自治区大兴安岭地区，控制流域面积16761km

表1选用的降水产品的基本信息

具体步骤如下：

(1)本次实例研究中选用新安江模型作为集总式水文模型，需首先根据小二沟流域的特性设置部分模型参数，如流域面积、计算步长等。本次实例研究中选用粒子群优化算法，需设置粒子群优化算法的运行参数，学习因子、惯性权重、迭代次数、搜索空间纬度等。

(2)参数率定：参数率定过程如图3所示，优化算法生成初始的模型参数粒子种群后，以4种近实时卫星降水产品驱动本发明提出的方法。当进行第m′次更新时，以T＝1，m′＝3为例的校正和预报过程如图2所示，涉及的各向量和矩阵的前三阶用于校正过程，利用最小二乘法解得校正系数向量A′，其中仅α

(3)模型验证：基于率定得到的模型参数最优参数，以4种近实时卫星降水产品驱动本发明提出的方法。当进行第m′次更新时，以T＝1，m′＝3为例的校正和预报过程如图2所示，涉及的各向量和矩阵的前三阶用于校正过程，利用最小二乘法解得校正系数向量A′，其中仅α

小二沟流域基于上述4种近实时卫星降水信息和产流约束校正的洪水预报方法，在率定期和验证期的预报效果评价指标如表2所示。可见除率定期的峰现时间误差和验证期的产流相对误差本发明方法略劣于未校正方法，其余所有基于本发明的评价指标在率定期和验证期都显著优于未校正方法的评价指标。此外，基于本发明的预报效果评价指标在率定期和验证期的表现稳定，说明基于本发明方法率定得到的参数具有代表性，能够有效描述流域的雨洪转化传递规律。

表2本发明在小二沟流域的率定期和验证期的预报效果评价指标

以上所述实施例仅表达本发明的实施方式，但并不能因此而理解为对本发明专利的范围的限制，应当指出，对于本领域的技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些均属于本发明的保护范围。