一种场次暴雨径流量预测方法与流程

1.本发明属于水文预测技术领域，具体涉及一种场次暴雨径流量预测方法。


背景技术：

2.降雨不仅是气候变化的主导因子，且是产流产沙的动力因子，使得流域径流量对降雨变化极为敏感，近年来，全国各地汛期时极端降雨也即暴雨频次和总量都呈现增加态势，暴雨危害和破坏性比普通降水更为明显，黄土高原的年径流量和年侵蚀量基本由少数几场暴雨引起，且暴雨事件在流域内易引起大水和大沙，易导致水库、淤地坝溃坝等风险，城镇地区的极端降水会导致城市防洪压力进一步加大，同时，由于城市的排水系统的设计标准偏低，导致我国城市内涝现象日益严重。
3.而现有技术中普遍将暴雨情形下的径流量预测与普通降雨特性认为一致，没有考虑到暴雨情形和普通降雨情形两者间性质的不同。
4.因此，如何准确地对场次暴雨产生的径流量进行预测，是本领域技术人员有待解决的技术问题。


技术实现要素：

5.本发明的目的是为了解决现有技术无法对场次暴雨所产生的径流量进行准确地预测，提出了一种场次暴雨径流量预测方法。
6.本发明的技术方案为：一种场次暴雨径流量预测方法，包括以下步骤：
7.s1、获取目标流域内的历史数据，所述历史数据具体为所述目标流域内所有雨量站和所有水文站中记载的降雨数据和洪水摘录表数据，并对所述历史数据进行预处理得到预处理数据；
8.s2、基于曼肯德尔法划分出所述预处理数据的现状期并得到现状期数据；
9.s3、通过r/s分析法在所述现状期数据中确定出预测数据；
10.s4、基于所述预测数据对场次暴雨径流量进行预测。
11.进一步地，所述步骤s1中的预处理具体包括以下步骤：
12.s11、将所述降雨数据和洪水摘录表数据统一时间步长；
13.s12、根据统一时间步长后的各站点对应的降雨数据确定所述目标流域中场次降雨对应的场次降雨量；
14.s13、根据统一时间步长后的各站点对应的洪水摘录表确定所述目标流域中场次降雨对应的场次洪量；
15.s14、将大于50mm的场次降雨量对应的场次降雨和场次洪量分别组成场次降雨序列和场次洪量序列，并将所述场次降雨序列和所述场次洪量序列组合为所述预处理数据。
16.进一步地，所述步骤s12中确定所述目标流域中场次降雨对应的场次降雨量具体通过泰森多边形法进行确定，所述泰森多边形法如下公式所示：
[0017][0018]
式中，p为场次降雨量，p
i
为流域内第i个雨量站雨量，α
i
为流域内第i个雨量站的面积权重系数，n为流域内雨量站个数。
[0019]
进一步地，所述步骤s2具体包括以下分步骤：
[0020]
s21、通过曼
‑
肯德尔法判断所述预处理数据是否通过显著性检验，若是，则执行步骤s22，若否，则执行步骤s23；
[0021]
s22、将所述预处理数据进行突变点检验，并以突变点检验结果为依据划分所述现状期以及所述现状期数据；
[0022]
s23、将所述预处理数据作为所述现状期，并将所述预处理数据作为现状期数据。
[0023]
进一步地，所述步骤s3具体为通过r/s分析法确定出所述现状期数据的赫斯特指数，并根据所述赫斯特指数确定出预测数据。
[0024]
进一步地，所述步骤s4具体包括以下分步骤：
[0025]
s41、基于所述预测数据，通过场次降雨量拟合法、降雨因子组合拟合法、场次降雨量分时段雨量拟合法和上包络法建立四种径流量预测计算式；
[0026]
s42、基于所述场次降雨量拟合法、降雨因子组合拟合法、场次降雨量分时段雨量拟合法和上包络法建立的四种径流量预测计算式分别对场次暴雨径流量进行预测，并得到四个预测结果；
[0027]
s43、基于所述四个预测结果确定所述场次暴雨径流量的预测区间，并将所述预测区间作为最终预测结果。
[0028]
进一步地，所述场次降雨量分时段雨量拟合法建立的径流量预测计算式具体如下：
[0029]
w＝k1p1 k2p2‑1 k3p3‑2 k4p4‑3 k5p5‑4 k6p6‑5 k7p
p
‑6 c；
[0030]
式中，w为径流量预测值，p1、p2‑1、p3‑2、p4‑3、p5‑4、p6‑5、p
p
‑6为各分时段雨量，k1、k2、k3、k4、k5、k6、k7、c为方程拟合系数。
[0031]
与现有技术相比，本发明具备以下有益效果：
[0032]
(1)本发明通过收集目标流域内各个水文气象站点的降雨径流数据；其次，根据曼
‑
肯德尔法次相关检验法和r/s分析法分别对流域内场次极端暴雨量和次洪量进行趋势检验和持续性检验；最后，利用次降雨量拟合法、降雨因子组合拟合法、次降雨分时段雨量拟合法、上包络线法4种方法对目标流域内场次暴雨径流量进行预测，本发明能够较为简单快速的对单次暴雨洪水进行预测，为政府决策部门防汛减灾、预测预警、调洪调度工作的开展以及流域水资源合理开发利用提供理论支撑和技术支持。
[0033]
(2)本发明得到四种预测公式，不同降雨因子和不同的降雨预测效果不同，通过不同公式预测径流量的区间，达到预测结果更加精准。
附图说明
[0034]
图1所示为本发明实施例提供的一种场次暴雨径流量预测方法的流程示意图。
具体实施方式
[0035]
下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
[0036]
本技术提出了一种场次暴雨径流量预测方法，如图1所示为一种场次暴雨径流量预测方法的流程示意图，该方法包括以下步骤：
[0037]
步骤s1、获取目标流域内的历史数据，所述历史数据具体为所述目标流域内所有雨量站和所有水文站中记载的降雨数据和洪水摘录表数据，并对所述历史数据进行预处理得到预处理数据。
[0038]
在本技术实施例中，所述步骤s1中的预处理具体包括以下步骤：
[0039]
s11、将所述降雨数据和洪水摘录表数据统一时间步长；
[0040]
s12、根据统一时间步长后的各站点对应的降雨数据确定所述目标流域中场次降雨对应的场次降雨量；
[0041]
s13、根据统一时间步长后的各站点对应的洪水摘录表确定所述目标流域中场次降雨对应的场次洪量；
[0042]
s14、将大于50mm的场次降雨量对应的场次降雨和场次洪量分别组成场次降雨序列和场次洪量序列，并将所述场次降雨序列和所述场次洪量序列组合为所述预处理数据。
[0043]
在具体的应用场景中，统一时间步长具体是等时间间隔数据，也即一小时，场次洪量的计算公式如下：
[0044][0045]
式中，w为次洪水量，亿m3；t为时段，h；i为次洪过程时段次序，n为次洪过程时段划分数量；q为流量，m3/s。
[0046]
所述步骤s12中确定所述目标流域中场次降雨对应的场次降雨量具体通过泰森多边形法进行确定，所述泰森多边形法如下公式所示：
[0047][0048]
式中，p为场次降雨量，p
i
为流域内第i个雨量站雨量，α
i
为流域内第i个雨量站的面积权重系数，n为流域内雨量站个数。
[0049]
步骤s2、基于曼
‑
肯德尔法划分出所述预处理数据的现状期并得到现状期数据。
[0050]
在本技术实施例中，所述步骤s2具体包括以下分步骤：
[0051]
s21、通过曼
‑
肯德尔法判断所述预处理数据是否通过显著性检验，若是，则执行步骤s22，若否，则执行步骤s23；
[0052]
s22、将所述预处理数据进行突变点检验，并以突变点检验结果为依据划分所述现状期以及所述现状期数据；
[0053]
s23、将所述预处理数据作为所述现状期，并将所述预处理数据作为现状期数据。
[0054]
在具体应用场景中，曼
‑
肯德尔法具体实施步骤如下：
[0055]
对于一个具有n个样本量的时间序列x，首先构造一秩序列：
[0056][0057]
式中：s
k
是第i时刻数值大于j时刻数值的累计个数，x
i
,x
j
分别是第i时刻和j时刻的样本值。
[0058]
定义统计量：
[0059][0060]
式中：uf1＝0；n为样本数；e(s
k
)是累计数s
k
的均值；var(s
k
)是累计数s
k
的方差；
[0061]
以α＝0.05为给定的显著性水平值，即u
0.05
＝
±
1.96。当uf
i
<0时，表明序列呈下降趋势；当uf
i
>0时，表明序列呈上升趋势。一般情况下，若uf
i
>u
0.05
时，认为序列存在明显的上升或下降趋势变化。
[0062]
对通过显著性检验的序列进行突变点检验，以突变点检验结果为依据进行时期的划分。若未通过趋势的显著性检验，则整个序列皆为现状期；认为序列受气候变化和人类活动的影响小，整个研究时期内未发生显著变化。
[0063]
若uf
i
>u
0.05
，则表明序列存在明显的趋势变化，将时间序列x按逆序排列，再按曼
‑
肯德尔法计算所得的uf
i
值定义为ub
i
，如果ub
i
和uf
i
两条曲线出现交点且交点在临界直线之间，那么交点对应的时刻就是突变时刻。
[0064]
步骤s3、通过r/s分析法在所述现状期数据中确定出预测数据。
[0065]
在本技术实施例中，所述步骤s3具体为通过r/s分析法确定出所述现状期数据的赫斯特指数，并根据所述赫斯特指数确定出预测数据。
[0066]
在具体应用场景中，对现状期下的场次暴雨量序列和径流量时间序{ζ(t)}t＝1,2,
…
,n，对于任意的正整数τ≥1，定义如下的统计量：
[0067]
均值序列：
[0068][0069]
累计离差：
[0070][0071]
极差：
[0072]
r(τ)＝maxx(t,τ)
‑
minx(t,τ)τ＝1,2
…
[0073]
标准差：
[0074][0075]
如果r/s＝r(τ)/s(τ)存在以下关系：
[0076]
r/s
∝
τ
h
[0077]
则序列存在hurst现象也即赫斯特指数，两边取对数，即log(r/s)＝logc hlogτ，c为某一常数值，log(r/s)～logτ图形中，直线斜率h即为hurst指数。0.5<hurst<1，表明未来趋势与现状下一致；若0<hurst<0.5，表明未来趋势与现状趋势相反。
[0078]
以r/s分析法检测结果判断数据是否具有预测性。若现状期数据通过r/s持续性检验，即0.5<hurst<1，表明未来暴雨径流趋势与现状期趋势一致，序列具有预测性。若0<hurst<0.5，表明未来暴雨径流与现状期趋势相反，序列不具有预测性。
[0079]
步骤s4、基于所述预测数据对场次暴雨径流量进行预测。
[0080]
在本技术实施例中，所述步骤s4具体包括以下分步骤：
[0081]
s41、基于所述预测数据，通过场次降雨量拟合法、降雨因子组合拟合法、场次降雨量分时段雨量拟合法和上包络法建立四种径流量预测计算式；
[0082]
s42、基于所述场次降雨量拟合法、降雨因子组合拟合法、场次降雨量分时段雨量拟合法和上包络法建立的四种径流量预测计算式分别对场次暴雨径流量进行预测，并得到四个预测结果；
[0083]
s43、基于所述四个预测结果确定所述场次暴雨径流量的预测区间，并将所述预测区间作为最终预测结果。
[0084]
在本技术实施例中，所述场次降雨量分时段雨量拟合法建立的径流量预测计算式具体如下：
[0085]
w＝k1p1 k2p2‑1 k3p3‑2 k4p4‑3 k5p5‑4 k6p6‑5 k7p
p
‑6 c；
[0086]
式中，w为径流量预测值，p1、p2‑1、p3‑2、p4‑3、p5‑4、p6‑5、p
p
‑6为各分时段雨量，k1、k2、k3、k4、k5、k6、k7、c为方程拟合系数。
[0087]
上式中，每场降雨过程中7个分时段最大雨量，即最大1h降雨量(p1)、最大2h降雨量减去最大1h降雨量(p2‑1)、
…
、最大6h降雨量减去最大5h降雨量(p6‑5)、最大6h降雨量之外的降雨量(p
p
‑6)。
[0088]
在具体应用场景中，次降雨量拟合法为以场次降雨量为横坐标，场次径流量为纵坐标，拟合现状期下的径流量预测计算式。
[0089]
降雨因子组合拟合法为影响场次降雨产流的因子包括降雨量及雨强，为同时考虑两者对径流量的影响，将场次降雨量与平均雨强的乘积pi作为自变量，径流量作为因变量，拟合现状期下降雨因子与径流量的计算式。
[0090]
上包络线法为尽可能考虑复杂不利条件的影响，选取场次降雨量与径流量关系图的上边缘点据，点绘现状期的上包络线并建立径流量预测计算式。
[0091]
需要说明的是，场次降雨量拟合法、降雨因子组合拟合法和上包络法建立的径流量预测计算式，是根据实测点数据，定义自变量和因变量后，放到统计软件中进行非线性拟合，最终以决定性系数为评判标准择优选择径流量预测计算式，例：次降雨量拟合法为以场次暴雨降雨量为横坐标，场次径流量为纵坐标，输入序列后，显示二次项计算式拟合较好，所以最后计算式为w1＝8
×
10
‑5p2‑
0.003p 0.13，不同的区域计算式不同，也有可能是呈指数关系，三次关系等，主要抓住不同的降雨因子作为自变量，这是由本领域技术人员根据实际情况进行选择实施的。
[0092]
为验证本发明所提出的方法能够对流域内暴雨情景下的径流量进行预测，选取汾
河静乐站控制流域作为验证流域，流域内包括怀道、宁化堡和海子背等雨量站，场次径流量利用汾河静乐站记录数据。
[0093]
数据选取的时间跨度为1971～2014年汛期场次产流降雨的相关数据，共89场数据，将其整理为时间间隔1h的降雨、径流过程数据，统计计算每场洪水的降雨量、洪量、降雨历时、平均雨强、各分时段雨量。
[0094]
利用mann
‑
kendall非参数秩次相关检验法对场次暴雨量和径流量序列进行了趋势性检验和突变性检验。场次极端暴雨量和对应的次洪量的uf
i
分别为0.61和1.74(表1)，序列有上升趋势，但都没有通过显著性检验，表明其上升趋势不明显。无突变点，整个序列作为现状期。
[0095]
表1流域场次暴雨洪水特征变化趋势检验
[0096][0097]
利用r/s分析法对场次暴雨量和径流量序列径流持续性分析。场次暴雨量和径流量序列的hurst指数分别为0.61和0.71，均大于0.5，表明流域内的降雨径流序列在未来一段时间内仍持续为不显著增长趋势。可以用现状降雨径流序列预测未来场次暴雨下的可能径流量。
[0098]
按照步骤s4建立次降雨量拟合法、降雨因子组合拟合法、次降雨分时段雨量拟合法、上包络线法对应的四种径流量预测计算式，见表2。并利用水文频率分析法和历史极端降雨组合法设计极端降雨情景，频率达到百年一遇和三百年一遇并带入计算式，与实际的同等量级场次降雨洪量对比发现，预测效果较好。
[0099]
表2径流量预测计算式
[0100][0101][0102]
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。