一种基于典型设计洪水过程线的小流域短期洪水预报方法与流程

本发明涉及洪水预报领域，特别是涉及一种基于典型设计洪水过程线的小流域短期洪水预报方法。

背景技术：

小流域通常是指集水面积不超过数百平方公里的小河小溪，但并无明确限制，一般认为流域面积在300～500km²以下可认为是小流域。从水文角度看小流域具有流域汇流以坡面汇流为主、水文资料缺乏、集水面积小等特性。

小流域设计洪水计算广泛应用于中、小型水利工程中，如修建农田水利工程的小水库、撇洪沟，渠系上交叉建筑物如涵洞、泄洪闸等，铁路、公路上的小桥涵设计，与大、中流域相比，小流域设计洪水具有以下四方面的特点：

(1)绝大多数小流域都没有水文站，即缺乏实测径流资料，甚至降雨资料也没有。因此小流域设计洪水计算一般为无资料情况下的计算。

(2)小流域面积小，自然地理条件趋于单一，拟定计算方法时，允许作适当的简化，即允许作出一些概化的假定。例如假定短历时的设计暴雨时空分布均匀。

(3)小流域分布广、数量多。因此，所拟定的计算方法，在保持一定精度的前提下，将力求简便，一般借助水文手册即可完成。

(4)小型工程一般对洪水的调节能力较小，工程规模主要受洪峰流量控制，因此对设计洪峰流量的要求高于对洪水过程线的要求。

在小流域中计算典型设计洪水过程线的方法，是由推理公式法或地区经验公式法计算小流域的设计洪峰流量，然后根据综合瞬时单位线法或多边形过程线法(三角形、五点形、多峰形等)，结合小流域的设计洪峰流量，计算小流域的典型设计洪水过程线。

由于我国小流域众多，小流域上的洪灾发生频率高、危害范围广、对国民经济影响严重。正因为小流域洪灾产生的灾害损失如此巨大，所以小流域短期洪水预报减少甚至避免洪灾的灾害损失尤为重要，合理的预报结果对于兴利防洪具有重大意义。

目前，常有的河道洪水预报方法有流量演算法和相应水位法。流量演算法的实质是成因分析法，主要分为两种，即水文学方法以及水力学方法。其中，水文学方法为马斯京根法和特征河长法；而水力学方法则采用解析法和数值法。相应水位法的实质是数理统计法，是研究相应水位(或相应流量)间的关系(统计规律)，依此为依据进行河段洪水预报的方法，主要包括洪峰水位预报法。

流量演算法中的水文学方法由于其参数具有较强经验型，在水流情况复杂的流域会面临参数难以正确率定的问题；而水力学方法存在所需不同类型资料数据较多，计算时间长等缺点。相应水位法则需要对洪水波进行大量的数理统计，数据精度要求高，统计计算费时，且主要预报洪峰处的水位。

技术实现要素：

本发明的目的在于克服现有技术的不足，提供一种基于典型设计洪水过程线的小流域短期洪水预报方法，对于小流域中某一河道，可以通过综合瞬时单位线法或多边形过程线法，结合不同年限的设计洪峰流量，获取相应的典型设计洪水过程线，利用上一时间点的实测流量确定预报时间，并根据典型设计洪水过程线，计算确定该河道逐时段的洪水预报结果，具有预报简单、快速的优势。

本发明的目的是通过以下技术方案来实现的：一种基于典型设计洪水过程线的小流域短期洪水预报方法，其特征在于：包括以下步骤：

s1.获取小流域河道的不同年限的设计洪峰流量和典型设计洪水过程线：

在小流域中，确定河道不同年限t1＜t2＜t3＜……＜tn-1＜tn的设计洪峰流量q1＜q2＜q3＜……＜qn-1＜qn；其中第n个年限tn对应的设计洪峰流量为qn，n＝1,2,3,…,n；

然后根据综合瞬时单位线法或多边形过程线法(三角形、五点形、多峰形等)，结合该河道不同年限的设计洪峰流量，计算确定该河道不同年限的典型设计洪水过程线；

s2.根据实测流量判定河道是否开始进行短期洪水预报：

获取河道的实测流量q，当实测流量q不低于某一年限tn的设计洪峰流量qn时，即满足条件q≥qn，就判定该河道开始进行短期洪水预报，进入步骤s3，此时的实测流量q为实测洪水流量；当实测流量q低于某一年限tn的设计洪峰流量qn时，即满足条件q＜qn，重复执行步骤s2，直至该河道开始进行短期洪水预报为止；

s3.根据河道的实测洪水流量确定短期洪水预报的预报时间：

每一个年限的设计洪峰流量都对应一个洪水预报时间间隔，当实测洪水流量q与设计洪峰流量qn、qn+1满足关系：qn≤q＜qn+1，n＝1,2,3,…,n-1，确定此时洪水的预报时间间隔为qn对应的δtn；根据实测时间tn和预报时间间隔δtn，结合公式tn+1＝tn+δtn，确定下一次的预报时间tn+1；

s4.根据相应年限的典型设计洪水过程线确定斜率kn：

在tn年限的设计洪水过程线上确定出tn时刻的洪水流量tn+1时刻的洪水流量根据斜率计算公式进行计算，当tn时刻实测流量相对于tn-1时刻实测流量降低，即qn-1＞qn，则kn＝-k；当tn时刻实测流量相对于tn-1时刻实测流量增加，qn-1＜qn时，则kn＝k；

s5.根据tn时刻的实测洪水流量qn和斜率kn，确定tn+1时刻预测洪水流量q′n+1：

根据直线方程公式q′n+1＝qn+kn(tn+1-tn)或q′n+1＝qn+knδtn,确定出tn+1时刻预测洪水流量q′n+1，这个预测值q′n+1是向社会公布的tn+1时刻的预测洪水流量数值；

s6.根据tn+1时刻的实测洪水流量判定是否可以结束短期洪水预报：

当实测洪水流量qn+1低于某一年限tn的设计洪峰流量qn时，即满足条件qn+1＜qn，判定该河道结束短期洪水预报，进入步骤s8；当实测洪水流量qn+1不低于某一年限tn的设计洪峰流量qn时，即满足条件qn+1≥qn，进入步骤s7；

s7.根据tn+1时刻的实测洪水流量qn+1和tn+1时刻的预测洪水流量q′n+1，对计算斜率kn的设计洪水过程线的年限tn校正；

根据比值关系式计算确定出预测洪水流量和实测洪水流量的比值a；

若90％≤a≤110％，预测洪水流量在实测洪水流量的合理范围内，则预测洪水流量满足精度要求，说明此时的河道洪水满足tn年限的设计洪水，则tn年限不需要校正，返回步骤s3；

若a＞110％，预测洪水流量高于实测洪水流量合理范围上限，预测洪水流量不满足精度要求；若此时的实测洪水流量比上一时刻增加，说明此时应满足低年限tn-1的设计洪水，则tn年限需要校正为tn-1，其中tn-1＜tn，并返回步骤s3；若此时的实测洪水流量比上一时刻减少，说明此时应满足高年限tn+1的设计洪水，则tn年限需要校正为tn+1，其中tn+1＞tn，返回步骤s3；

若a＜90％，预测洪水流量低于实测洪水流量合理范围下限，预测洪水流量不满足精度要求；若此时的实测洪水流量比上一时刻增加，说明此时应满足高年限tn+1的设计洪水，则tn年限需要校正为tn+1，其中tn+1＞tn，返回步骤s3；若此时的实测洪水流量比上一时刻减少，说明此时应满足低年限tn-1的设计洪水，则tn年限需要校正为tn-1，其中tn-1＜tn，返回第3步；

s8.结束短期洪水预报。

进一步地，步骤s2中作为判定条件的设计洪峰流量的相应年限tn要根据河道特征及保护对象重要性预先进行选择确定。

进一步地，设计洪峰流量qn越大，或者tn年限越大，洪水预报的时间间隔δtn就越小(δt1＞δt2＞δt3＞……＞δtn-1＞δtn)。

进一步地，步骤s4中tn时刻对应的实测洪水流量为qn，在tn年限设计洪水过程线上确定的tn时刻洪水流量为tn+1时刻对应的实测洪水流量为qn+1，在tn年限设计洪水过程线上确定的tn+1时刻洪水流量为

进一步地，步骤s4中第一次确定计算预测洪水流量的斜率kn的年限tn，可根据步骤s2中判定短期洪水预报的年限tn确定。

进一步地，步骤s5中tn时刻对应的实测流量qn可以预先向相关部门获取，在本申请中为已知参数；步骤s6中该河道tn+1时刻的实测洪水流量qn+1可以预先向相关部门获取，在本申请中为已知参数。

进一步地，步骤s6中作为判定条件的设计洪峰流量的年限tn要根据河道特征及保护对象重要性进行预先选择确定。

进一步地，步骤s7中设计洪水过程线的年限tn一般是用相邻年限tn-1或tn+1进行校正，若预测洪水流量计算误差太大(超过设定误差范围)，可考虑对tn用tn-2、tn+2或其它年限进行校正。

本发明的有益效果是：本发明提供的基于典型设计洪水过程线的小流域短期洪水预报方法，对小流域某河道通过水文方法获取设计洪峰流量和典型设计洪水过程线，通过实测流量确定是否开始和结束洪水预报以及预报时间，并根据典型设计洪水线确定预测洪水流量，最后对典型设计洪水过程线的年限进行修正，获得整场洪水的预报结果。该方法比较现有的流量演算法和相应水位法具有简单、快速等特点。

附图说明

图1为本发明的方法流程图。

图2为实施例中某河道不同年限的设计洪水过程线。

图3为实施例中小流域某河道，某场洪水的预测洪水流量和实测洪水流量的对比图，其中(a)(b)分别对应于洪水预报的两种方案，方案的预报时间间隔不同。

具体实施方式

下面结合附图进一步详细描述本发明的技术方案，但本发明的保护范围不局限于以下所述。

如图1所示，基于典型设计洪水过程线的小流域短期洪水预报方法，对小流域中某一河段进行两场洪水预报，主要包括以下几步：

s1.获取小流域河道的不同年限的设计洪峰流量和典型设计洪水过程线：

对这一河道，由推理公式法或地区经验公式法计算该河道6种年限(2、5、10、20、30、50年)的设计洪峰流量(q1＜q2＜q3＜……＜q6)，然后根据综合瞬时单位线法或多边形过程线法(三角形、五点形、多峰形等)，结合该河道6种年限的设计洪峰流量，计算确定该河道相应的典型设计洪水过程线。图2为该河道的6种典型设计洪水过程线。

s2.根据实测流量判定河道可以进行短期洪水预报：

向相关部门获取该河道的实测流量q，当实测流量q不低于2年一遇的设计洪峰流量238m³/s时，即满足条件q≥238m³/s，就判定该河道开始进行短期洪水预报，进行步骤s3，此时的实测流量q为实测洪水流量；当实测流量q低于2年一遇的设计洪峰流量238m³/s时，即满足条件q＜238m³/s，重复步骤s2，直至该河道可以进行短期洪水预报为止。

s3.根据河道的实测洪水流量确定短期洪水预报的预报时间：

每一个年限的设计洪峰流量都对应一个洪水预报时间间隔，当实测洪水流量q与设计洪峰流量qn、qn+1满足关系：qn≤q＜qn+1，可以确定此时洪水的预报时间间隔为δtn。根据实测时间tn和预报时间间隔δtn，结合公式tn+1＝tn+δtn，确定下一次的预报时间tn+1。

表1为不同年限的设计洪峰流量对应的预报时间间隔的两种方案。

s4.根据相应年限的典型设计洪水过程线确定斜率kn：

在tn年限的设计洪水过程线上可以确定出tn时刻的洪水流量tn+1时刻的洪水流量根据斜率计算公式进行计算，当tn时刻实测流量相对于tn-1时刻实测流量降低，即qn-1＞qn，则kn＝-k；当tn时刻实测流量相对于tn-1时刻实测流量增加，qn-1＜qn时，则kn＝k。计算初始斜率kn的设计洪水过程线的年限取为2年。

s5.根据tn时刻的实测洪水流量qn和斜率kn，确定tn+1时刻预测洪水流量q′n+1：

根据直线方程公式q′n+1＝qn+kn(tn+1-tn)或q′n+1＝qn+knδtn,确定出tn+1时刻预测洪水流量q′n+1，这个预测值q′n+1是向社会公布的tn+1时刻的预测洪水流量数值。

s6.根据tn+1时刻的实测洪水流量判定是否可以结束短期洪水预报：

当实测洪水流量qn+1低于2年一遇的设计洪峰流量238m³/s时，即满足条件qn+1＜238m³/s，就可以判定该河道可以结束短期洪水预报，进行步骤s8；当实测洪水流量qn+1不低于2年一遇的设计洪峰流量238m³/s时，即满足条件qn+1≥238m³/s，进行步骤s7。

s7.根据tn+1时刻的实测洪水流量qn+1和tn+1时刻的预测洪水流量q′n+1，对计算斜率kn的设计洪水过程线的年限tn校正。

根据比值关系式计算确定出预测洪水流量和实测洪水流量的比值a。

若90％≤a≤110％，预测洪水流量在实测洪水流量的合理范围内，则预测洪水流量满足精度要求，说明此时的河道洪水满足tn年限的设计洪水，则tn年限不需要校正，返回步骤s3；

若a＞110％，预测洪水流量高于实测洪水流量合理范围上限，预测洪水流量不满足精度要求。若此时的实测洪水流量比上一时刻增加，说明此时应满足低年限tn-1的设计洪水，则tn年限需要校正为tn-1(tn-1＜tn)，返回步骤s3；若此时的实测洪水流量比上一时刻减少，说明此时应满足高年限tn+1的设计洪水，则tn年限需要校正为tn+1(tn+1＞tn)，返回步骤s3；

若a＜90％，预测洪水流量低于实测洪水流量合理范围下限，预测洪水流量不满足精度要求。若此时的实测洪水流量比上一时刻增加，说明此时应满足高年限tn+1的设计洪水，则tn年限需要校正为tn+1(tn+1＞tn)，返回步骤s3；若此时的实测洪水流量比上一时刻减少，说明此时应满足低年限tn-1的设计洪水，则tn年限需要校正为tn-1(tn-1＜tn)，返回步骤s3；

s8.结束短期洪水预报。

如图3所示，在本申请的实施例中，图3(a)位方案一的预测洪水过程图，在方案一中，相应设计年限2、5、10、20、30、50年预报的时间间隔分别为60、50、40、30、20、10分钟；图3(b)位方案二的预测洪水过程图，在方案二中，相应设计年限2、5、10、20、30、50年预报的时间间隔分别为30、25、20、15、10、5分钟，两种方案的设计洪峰流量相等，方案二比方案一的预报时间间隔减小。结果表明：方案一中的预测和实测洪量的比值为0.999，预测和实测洪水流量的比值基本满足90％到110％，但预测和实测最大洪水流量的比值为0.807；方案二中的预测和实测洪量的比值为1.002，预测和实测洪水流量的比值基本满足90％到110％，且预测和实测最大洪水流量的比值为0.954。通过两种方案对比可知，将预报时间间隔减小，洪水预报具有良好的预测性，预测洪水流量和洪量的计算精度基本满足要求，该方法的预报结果已能较好地描述整场洪水的发展过程,且计算方法省时简单，相比于流量演算法和相应水位法具有简单、快速等特点。

以上所述是本发明的优选实施方式，应当理解本发明并非局限于本文所披露的形式，不应该看作是对其他实施例的排除，而可用于其他组合、修改和环境，并能够在本文所述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。