CN111914432A - 一种基于大数据的水文预报方法 - Google Patents

Abstract

本发明公开了一种基于大数据的水文预报方法，它包括：步骤1、将整个流域按范围划分为各个基本计算单元，步骤2、对每个计算单元进行产流计算；步骤3、根据产流计算结果对每个计算单元进行汇流计算；步骤4、对每个计算单元输出流量过程进行河网演进叠加计算，最终得到整个流域断面出口流量；根据流域断面出口流量水文预报；解决了现有技术中水文预报主要采用新安江模型存在的由于模型的参数得不到准确的率定，最终导致新安江模型计算的预报结果存在准确度较差等问题；参数率定的难度和应用地区的限制，导致新安江模型很难向无水文资料地区进行大力的技术推广和应用等技术问题。

Description

一种基于大数据的水文预报方法

技术领域

本发明属于水文预报技术，尤其涉及一种基于大数据的水文预报方法。

背景技术

水文预报对防洪、抗旱和水资源合理利用等领域中有着重要意义。在现有技术中，我国一般常用的水文预报模型是新安江模型。目前，新安江模型针对无资料地区和干旱地区的水文预报存在一定的限制。在新安江模型中，土壤蒸发机制采用的是三层蒸发模式，对于无水文资料的地区，模型的蒸发参数较难率定，很难通过流域下垫面因子对蒸发参数进行率定估值；另外，新安江模型的产流计算机制只适用于湿润地区，而不适应干旱地区。参数率定的难度和应用地区的限制，导致新安江模型很难大范围普及，尤其是向无资料地区进行技术推广和应用；同时，由于模型的参数得不到准确的率定，最终导致新安江模型计算的预报结果存在准确度较差等问题。

发明内容

本发明要解决的技术问题是：提供一种基于大数据的水文预报方法，以解决现有技术中水文预报主要采用新安江模型存在的由于模型的参数得不到准确的率定，最终导致新安江模型计算的预报结果存在准确度较差等问题；参数率定的难度和应用地区的限制，导致新安江模型很难向无水文资料地区进行大力的技术推广和应用等技术问题。

本发明的技术方案是：

一种基于大数据的水文预报方法，它包括：

步骤1、将整个流域按范围划分为各个基本计算单元，

步骤2、对每个计算单元进行产流计算；

步骤3、根据产流计算结果对每个计算单元进行汇流计算；

步骤4、对每个计算单元输出流量过程进行河网演进叠加计算，最终得到整个流域断面出口流量；根据流域断面出口流量水文预报。

它还包括：步骤5、对最大下渗率F₀、稳定下渗率F_c、霍顿下渗公式系数Beta、流域土壤田间持水能力Wm和参数Em进行率定。

步骤1所述基本计算单元的划分方法包括：将流域内互不嵌套的分水岭作为基本计算单元；一个计算单元对应一个计算点，计算点分为叶子和区间两种类型，上游计算点通过马斯京根算法与下游计算点进行连接。

每个计算单元的降水输入采用气象部门的5平方公里降雨栅格数据；按流域内各个分水岭范围进行分别提取。

步骤2所述对每个计算单元进行产流计算的方法包括：

步骤2.1、当降水量或降雨量h到达地面时，按流域蒸发能力Em进行降水扣损,得到扣损后的降雨量PE，损失的降雨量为记为E_S，即PE+E_S＝h；PE<0时激活土壤蒸发计算，并置PE为0；激活土壤蒸发计算后，执行步骤2.4完成一次产流迭代计算；

步骤2.2、当降水扣损后，如果PE>0，则按水文标准算法，通过流域蓄水容量曲线进行蓄满产流计算；输出净雨R和蓄满吸收量DW，将输出净雨R通过二水源算法得到地上净雨Rs和地下净雨hg；得到地上净雨Rs后，根据霍顿下渗公式结合土壤含水量W，推算下渗率f_v；

步骤2.3、根据当前下渗率f_v和产流面积比a，进行第一次壤中净雨分配，得出壤中净雨的第一部分dhss1和地表净雨hs；再进行第二次壤中净雨分配，从蓄满吸收量DW得到第二部分壤中净雨dhss2，并修正DW；最后合并dhss1和dhss2得到最终的壤中净雨hss；

步骤2.4、通过渗漏机制，根据壤中自由水转换率SSCR和地下水渗漏转换率GCR，实现壤中净雨h_ss向地下净雨h_g，以及地下净雨h_g向深层地下水h_dg的转换。

对计算单元进行汇流计算的内容包括：地表汇流Ds、壤中汇流Dss和地下汇流Dg的计算；地表汇流采用瞬时单位线方法实现，壤中汇流Dss和地下汇流Dg采用线性水库方法实现；

将地表净雨通过瞬时单位线转换为流量过程的算法：

其中，U_i为瞬时单位线第i时刻对应的值；D_ui为地表净雨h_s对应的第i时刻的流量值，F为流域面积。

所述进行河网演进叠加的方法为马斯京根河道演进算法。

对于一个区间结点，出流量Q等于上游入流Q_uI演进到结点的出流Q_uO加结点区间的总汇流量D，即：

D_i＝D_si+D_ssi+D_gi

Q_i＝Q_uOi+D_i

其中，Q_uOi-1表示区间河道第i-1时刻演进出流量，Q_uIi表示区间河道第i时刻入流量,

为平均流速，X为流量比重系数，L为区间河道长度，

Q_uOi表示区间河道第i时刻的出流量。

最大下渗率F₀、稳定下渗率F_c、霍顿下渗公式系数Beta、流域土壤田间持水能力Wm和参数Em的率定方法包括：

根据土壤质地和概化土壤间隙度，推导最大下渗率F₀、稳定下渗率F_c和霍顿下渗公式系数Beta；

根据霍顿下渗公式试算Wm，并引用土地利用类型的影响系数进行修正试算；依据土壤含水量分区和实测样本进行验证，并进行相应的反馈优化；一直循环进行参数优化，直到获得准确的成果，最后输出Wm、F₀、F_c和Beta栅格成果；

所述参数Em率定方法为：

基于30年雨站点统计计算干旱指数和汛期权重系数；根据经验关系拟定干旱指数曲线；

根据干旱指数公式，反推年蒸发能力，并结合汛期权重系数计算汛期蒸发能力；

对比实测样本，然后反馈优化干旱指数曲线；

得到结果后，再加入土地利用影响系数进行进一步优化；直到获得最终成果；最后输出Em栅格。

所述土壤蒸发的计算方法为：

建立土壤含水量w与土壤蒸发量E_w的相关曲线，简称土壤蒸发曲线；

通过土壤蒸发曲线中的Em、Wm和Wb来反应不同流域的蒸发特性；Em为流域日蒸发能力(mm/d)，Wm为流域土壤田间持水能力(mm)，Wb为流域土壤毛管断裂含水量(mm)。土壤蒸发曲线的数学公式如下：

式中，K₁和K₂是全局的归一化形状系数，用于控制曲线在Wm和Wb两个区域的弯曲度；

根据流域某个时刻的土壤含水量w计算这个时刻的土壤蒸发量Ew；土壤蒸发量计算必须限制Ew+Es<＝Em。

当前下渗率f_v推导逻辑表达式为：

f_v＝F(W,F₀,F_c,Beta)

其中，F₀为最大下渗率；F_c为稳定下渗率；Beta为霍顿下渗公式系数。

第一次壤中净雨分配公式为：

hs＝Rs-dhss1

第二次壤中净雨计算公式为：

hss＝dhss1+dhss2

以上步骤在计算hss时，需要将hhs累加到壤中自由水总量Wss，同时Wss应该受壤中自由水总量上限MaxWss的控制，如果超出上限，则需要将hss重新修正，并将多余部分归为地表净雨hs；MaxWss的值等于WM和壤中自由水容量系数WSSC的乘积；

dhss＝F(W,WM,F_c,SSCR,h_ss,h_g)

h_dg＝F(F_c,GCR,h_g)

上式中，dhss为壤中净雨到地下净雨的转换量；h_dg为地下净雨到深层地下水的转换量；计算中设定深层地下水不再产生汇流量。

对于每次产流迭代计算，必须根据水量平衡原理保证下述公式成立：

PE＝h_s+h_ss+h_g+h_dg。

本发明有益效果：

本发明基于传统水文计算的基本原理，并参考了新安江模型的部分算法思想重新构建而成的水文计算方法；本发明将土壤、植被、土地利用和地质等下垫面的物理因素与EC模型的参数体系进行有效关联，从而解决无资料地区水文预报的技术难题；本发明水源划分采用了分层模式，并结合流域下渗曲线和动态渗漏算法实现了水源的划分机制；其中，稳定下渗率F_c用于控制地下水分配，下渗率f_v用于控制壤中净雨分配量，并通过WSSC对壤中流线性水库进行总量控制；SSCR和GCR两个参数是与稳定下渗F_c严格相关的，在一定程度上可以用于模拟干旱地区的产流现象；本发明基于瞬时单位线，根据净雨强度动态计算，提升了流域计算单元的划分密度。

本发明充分利用现代GIS(遥感遥测和地理信息)系统技术，利用GIS提取流域地形、土壤、植被的强大功能，建立具有物理概念基础的大数据的水文预报方法，使计算参数客观化，尽量避免人为调试或选取参数的主观性，从而加强算法的通用性，以解决无资料地区中小流域洪水预报的困难；采用大数据分析方法对参数进行反馈优化率定，提高了算法参数适应度，有助于向无资料地区进行技术推广与应用。

本发明优点：

提升了模型降雨输入的精度。

通过土壤蒸发曲线简化了相关蒸发参数的估值和率定工作。

充分利用土壤的下渗曲线，通过额外的渗漏机制实现三水源划分，提升了模型对干旱地区产流计算的适应度。

通过大数据分析方法，建立了模型参数与土壤、植被、土地利用和地质等下垫面物理因素的映射关系；为无资料地区的洪水预报，提供了技术支撑。

解决了现有技术中水文预报主要采用新安江模型存在的由于模型的参数得不到准确的率定，最终导致新安江模型计算的预报结果存在准确度较差等问题；参数率定的难度和应用地区的限制，导致新安江模型很难向无水文资料地区进行大力的技术推广和应用等技术问题。

具体实施方式

一种基于大数据的水文预报方法，它包括：

步骤1、计算单元划分：将整个流域按范围划分为各个基本计算单元；

步骤1所述将整个流域按范围划分为各个基本计算单元的方法包括：将互不嵌套的分水岭作为基本计算单元；一个计算单元对应一个计算点，计算点分为叶子和区间两种类型，上游计算点通过马斯京根算法与下游计算点进行连接。

本发明每个计算单元的降水输入采用气象部门的5平方公里降雨栅格数据；按流域内各个分水岭范围进行分别提取；而传统新安江模型应用一般是根据雨量站点分布，按泰森多边形进行降雨提取；因此本发明的水文预报精度相对现有技术有很大的提高。

本发明利用30米DEM数据，建立了全流域的水系和30km^2以下小流域的分水岭的数据集成果，计算单元划分比传统新安江模型更加合理和精细。

步骤2、对每个计算单元进行产流计算；

步骤2所述对每个计算单元进行产流计算的方法包括：

步骤2.1、当降水量或降雨量h到达地面时，按流域蒸发能力Em进行降水扣损,得到扣损后的降雨量PE，损失的降雨量为记为E_S,即PE+E_S＝h；

PE<0时激活土壤蒸发计算，并置PE为0；激活土壤蒸发计算后，执行步骤2.4完成一次产流迭代计算。

所述土壤蒸发的计算方法为：

建立土壤含水量w与土壤蒸发量E_w的相关曲线，简称土壤蒸发曲线；

通过土壤蒸发曲线中的Em、Wm和Wb来反应不同流域的蒸发特性；Em为流域日蒸发能力(mm/d)，Wm为流域土壤田间持水能力(mm)，Wb为流域土壤毛管断裂含水量(mm)。土壤蒸发曲线的数学公式如下：

式中，K₁和K₂是全局的归一化形状系数，用于控制曲线在Wm和Wb两个区域的弯曲度。

在产流计算中，可以根据流域某个时刻的土壤含水量w计算这个时刻的土壤蒸发量Ew；土壤蒸发量计算必须限制Ew+Es<＝Em，否则要修正Ew的值来保证Ew+Es＝Em。

步骤2.2、当降水扣损后，如果PE>0，则按水文标准算法，通过流域蓄水容量曲线进行蓄满产流计算；输出净雨R和蓄满吸收量DW，将输出净雨R通过二水源算法得到地上净雨Rs和地下净雨hg；得到地上净雨Rs后，根据霍顿下渗公式结合土壤含水量W，推算下渗率f_v。

步骤2.3、根据当前下渗率f_v和产流面积比a，进行第一次壤中净雨分配，得出壤中净雨的第一部分dhss1和地表净雨hs；再进行第二次壤中净雨分配，从蓄满吸收量DW得到第二部分壤中净雨dhss2，并修正DW；最后合并dhss1和dhss2得到最终的壤中净雨hss。

当前下渗率f_v推导逻辑表达式为：

f_v＝F(W,F₀,F_c,Beta)

其中，F₀为最大下渗率；F_c为稳定下渗率；Beta为霍顿下渗公式系数。

第一次壤中净雨分配公式为：

hs＝Rs-dhss1

第二次壤中净雨计算公式为：

hss＝dhss1+dhss2

以上步骤在计算hss时，需要将hhs累加到壤中自由水总量Wss，同时Wss应该受壤中自由水总量上限MaxWss的控制，如果超出上限，则需要将hss重新修正，并将多余部分归为地表净雨hs。MaxWss的值等于WM和壤中自由水容量系数WSSC的乘积。

步骤2.4、通过渗漏机制，根据壤中自由水转换率SSCR和地下水渗漏转换率GCR，实现壤中净雨h_ss向地下净雨h_g，及地下净雨h_g向深层地下水h_dg的转换。

dhss＝F(W,WM,F_c,SSCR,h_ss,h_g)

h_dg＝F(F_c,GCR,h_g)

上式中，dhss为壤中净雨到地下净雨的转换量；h_dg为地下净雨到深层地下水的转换量。模型中，设定深层地下水不再产生汇流量。

对于每次产流迭代计算，必须根据水量平衡原理保证下述公式成立：

PE＝h_s+h_ss+h_g+h_dg

步骤3、根据产流计算结果对每个计算单元进行汇流计算；

对计算单元进行汇流计算的内容包括：地表汇流Ds、壤中汇流Dss和地下汇流Dg的计算。地表汇流采用瞬时单位线方法实现，壤中汇流Dss和地下汇流Dg采用线性水库方法实现。

将地表净雨通过瞬时单位线转换为流量过程的算法：

其中，U_i为瞬时单位线第i时刻对应的值；D_ui为地表净雨h_s对应的第i时刻的流量值，F为流域面积。

瞬时单位线确定的关键问题是确定流域洪水滞时TL；

本发明采用公式：

T_L＝T(H,F,C,J,Ke,LLR)

进行计算；其中T(X)为滞时概化率定函数，H为地表净雨强度，F为流域面积，C为坡面流速系数，J为流域平均坡度比，Ke为形状系数，LLR为蓄滞系数。

步骤4、对每个计算单元输出流量过程进行河网演进叠加计算，最终得到整个流域断面出口流量。

所述进行河网演进叠加的方法为马斯京根河道演进算法。

对于一个区间结点，出流量Q等于上游入流Q_uI演进到结点的出流Q_uO加结点区间的总汇流量D。即：

D_i＝D_si+D_ssi+D_gi

Q_i＝Q_uOi+D_i

其中，Q_uOi-1表示区间河道第i-1时刻演进出流量，Q_uIi表示区间河道第i时刻入流量,

为平均流速，X为流量比重系数，L为区间河道长度，Q_uOi表示区间河道第i时刻的出流量。

最大下渗率F₀、稳定下渗率F_c、霍顿下渗公式系数Beta和流域土壤田间持水能力Wm的率定：

1、根据土壤质地，概化土壤间隙度，推导最大下渗率F₀、稳定下渗率F_c和霍顿下渗公式系数Beta。

2、根据霍顿下渗公式试算Wm，并引用土地利用类型的影响系数进行修正试算。

3、依据中国土壤含水量分区和实测样本进行验证，并进行相应的反馈优化。

4、多次循环，进行参数优化。直到获得较满意的成果。最后输出Wm、F₀、F_c和Beta栅格成果。

参数Em率定

1、基于30年雨站点统计计算干旱指数和汛期权重系数；根据经验关系拟定干旱指数曲线。

2、根据干旱指数公式，反推年蒸发能力，并结合汛期权重系数计算汛期蒸发能力。

3、加入模型演算，并对比实测样本，然后反馈优化干旱指数曲线。

4、得到较理想结果后，再加入土地利用影响系数进行进一步优化。

5、执行多次优化步骤，直到获得满意成果。最后，输出Em栅格。

参数的率定作为模型应用的一个重要环节，大部分的流域水文模型特别是中小流域的一些参数不能直接通过观测试验确定。它们的取值却与流域的下垫面特征有着一定的关系，但现实中却不能直接与流域的下垫面特征建立起关系，所以对于流域水文模型来说参数的确定仍然是一个比较困难的问题。

本发明借助大数据分析的方法，基于一些理论数学公式，通过不断的迭代和反馈试算，获取一个相对合理的参数估值结果；以解决现有技术的不足。

Claims (10)

1.一种基于大数据的水文预报方法，它包括：

步骤1、将整个流域按范围划分为各个基本计算单元，

步骤2、对每个计算单元进行产流计算；

步骤3、根据产流计算结果对每个计算单元进行汇流计算；

步骤4、对每个计算单元输出流量过程进行河网演进叠加计算，最终得到整个流域断面出口流量；根据流域断面出口流量水文预报。

2.根据权利要求1所述的一种基于大数据的水文预报方法，其特征在于：它还包括：步骤5、对最大下渗率F₀、稳定下渗率F_c、霍顿下渗公式系数Beta、流域土壤田间持水能力Wm和参数Em进行率定。

3.根据权利要求1所述的一种基于大数据的水文预报方法，其特征在于：步骤1所述基本计算单元的划分方法包括：将流域内互不嵌套的分水岭作为基本计算单元；一个计算单元对应一个计算点，计算点分为叶子和区间两种类型，上游计算点通过马斯京根算法与下游计算点进行连接。

4.根据权利要求1所述的一种基于大数据的水文预报方法，其特征在于：每个计算单元的降水输入采用气象部门的5平方公里降雨栅格数据；按流域内各个分水岭范围进行分别提取。

5.根据权利要求1所述的一种基于大数据的水文预报方法，其特征在于：步骤2所述对每个计算单元进行产流计算的方法包括：

步骤2.1、当降水量或降雨量h到达地面时，按流域蒸发能力Em进行降水扣损,得到扣损后的降雨量PE，损失的降雨量为记为E_S，即PE+E_S＝h；PE<0时激活土壤蒸发计算，并置PE为0；激活土壤蒸发计算后，执行步骤2.4完成一次产流迭代计算；

步骤2.2、当降水扣损后，如果PE>0，则按水文标准算法，通过流域蓄水容量曲线进行蓄满产流计算；输出净雨R和蓄满吸收量DW，将输出净雨R通过二水源算法得到地上净雨Rs和地下净雨hg；得到地上净雨Rs后，根据霍顿下渗公式结合土壤含水量W，推算下渗率f_v；

步骤2.3、根据当前下渗率f_v和产流面积比a，进行第一次壤中净雨分配，得出壤中净雨的第一部分dhss1和地表净雨hs；再进行第二次壤中净雨分配，从蓄满吸收量DW得到第二部分壤中净雨dhss2，并修正DW；最后合并dhss1和dhss2得到最终的壤中净雨hss；

步骤2.4、通过渗漏机制，根据壤中自由水转换率SSCR和地下水渗漏转换率GCR，实现壤中净雨h_ss向地下净雨h_g，以及地下净雨h_g向深层地下水h_dg的转换。

6.根据权利要求1所述的一种基于大数据的水文预报方法，其特征在于：对计算单元进行汇流计算的内容包括：地表汇流Ds、壤中汇流Dss和地下汇流Dg的计算；地表汇流采用瞬时单位线方法实现，壤中汇流Dss和地下汇流Dg采用线性水库方法实现；

将地表净雨通过瞬时单位线转换为流量过程的算法：

其中，U_i为瞬时单位线第i时刻对应的值；D_ui为地表净雨h_s对应的第i时刻的流量值，F为流域面积。

7.根据权利要求1所述的一种基于大数据的水文预报方法，其特征在于：所述进行河网演进叠加的方法为马斯京根河道演进算法。

对于一个区间结点，出流量Q等于上游入流Q_uI演进到结点的出流Q_uO加结点区间的总汇流量D，即：

D_i＝D_si+D_ssi+D_gi

Q_i＝Q_uOi+D_i

其中，Q_uOi-1表示区间河道第i-1时刻演进出流量，Q_uIi表示区间河道第i时刻入流量,

为平均流速，X为流量比重系数，L为区间河道长度，Q_uOi表示区间河道第i时刻的出流量。

8.根据权利要求2所述的一种基于大数据的水文预报方法，其特征在于：最大下渗率F₀、稳定下渗率F_c、霍顿下渗公式系数Beta、流域土壤田间持水能力Wm和参数Em的率定方法包括：

根据土壤质地和概化土壤间隙度，推导最大下渗率F₀、稳定下渗率F_c和霍顿下渗公式系数Beta；

根据霍顿下渗公式试算Wm，并引用土地利用类型的影响系数进行修正试算；依据土壤含水量分区和实测样本进行验证，并进行相应的反馈优化；一直循环进行参数优化，直到获得准确的成果，最后输出Wm、F₀、F_c和Beta栅格成果；

所述参数Em率定方法为：

基于30年雨站点统计计算干旱指数和汛期权重系数；根据经验关系拟定干旱指数曲线；

根据干旱指数公式，反推年蒸发能力，并结合汛期权重系数计算汛期蒸发能力；

对比实测样本，然后反馈优化干旱指数曲线；

得到结果后，再加入土地利用影响系数进行进一步优化；直到获得最终成果；最后输出Em栅格。

9.根据权利要求5所述的一种基于大数据的水文预报方法，其特征在于：所述土壤蒸发的计算方法为：

建立土壤含水量w与土壤蒸发量E_w的相关曲线，简称土壤蒸发曲线；

通过土壤蒸发曲线中的Em、Wm和Wb来反应不同流域的蒸发特性；Em为流域日蒸发能力(mm/d)，Wm为流域土壤田间持水能力(mm)，Wb为流域土壤毛管断裂含水量(mm)。土壤蒸发曲线的数学公式如下：

式中，K₁和K₂是全局的归一化形状系数，用于控制曲线在Wm和Wb两个区域的弯曲度；

根据流域某个时刻的土壤含水量w计算这个时刻的土壤蒸发量Ew；土壤蒸发量计算必须限制Ew+Es<＝Em。

10.根据权利要求5所述的一种基于大数据的水文预报方法，其特征在于：下渗率f_v推导逻辑表达式为：

f_v＝F(W,F₀,F_c,Beta)

其中，F₀为最大下渗率；F_c为稳定下渗率；Beta为霍顿下渗公式系数。

第一次壤中净雨分配公式为：

hs＝Rs-dhss1

其中，f_v为当前下渗率，通过下渗曲线和土壤含水量推导取值，a为产流面积比；

第二次壤中净雨计算公式为：

hss＝dhss1+dhss2

以上步骤在计算hss时，需要将hhs累加到壤中自由水总量Wss，同时Wss应该受壤中自由水总量上限MaxWss的控制，如果超出上限，则需要将hss重新修正，并将多余部分归为地表净雨hs；MaxWss的值等于WM和壤中自由水容量系数WSSC的乘积；

dhss＝F(W,WM,F_c,SSCR,h_ss,h_g)

h_dg＝F(F_c,GCR,h_g)

上式中，dhss为壤中净雨到地下净雨的转换量；h_dg为地下净雨到深层地下水的转换量；计算中设定深层地下水不再产生汇流量。

对于每次产流迭代计算，根据水量平衡原理保证下述公式成立：

PE＝h_s+h_ss+h_g+h_dg。