CN113589404B - 一种场次暴雨径流量预测方法 - Google Patents
一种场次暴雨径流量预测方法 Download PDFInfo
- Publication number
- CN113589404B CN113589404B CN202110868146.8A CN202110868146A CN113589404B CN 113589404 B CN113589404 B CN 113589404B CN 202110868146 A CN202110868146 A CN 202110868146A CN 113589404 B CN113589404 B CN 113589404B
- Authority
- CN
- China
- Prior art keywords
- rainfall
- data
- runoff
- field
- prediction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/14—Rainfall or precipitation gauges
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
Landscapes
- Environmental & Geological Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Hydrology & Water Resources (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种场次暴雨径流量预测方法,通过收集目标流域内各个水文气象站点的降雨径流数据;其次,根据曼‑肯德尔法秩次相关检验法和R/S分析法分别对流域内场次极端暴雨量和次洪量进行趋势检验和持续性检验;最后,利用次降雨量拟合法、降雨因子组合拟合法、次降雨分时段雨量拟合法、上包络线法4种方法对目标流域内场次暴雨径流量进行预测,本发明能够较为简单快速的对单次暴雨洪水进行预测,为政府决策部门防汛减灾、预测预警、调洪调度工作的开展以及流域水资源合理开发利用提供理论支撑和技术支持。
Description
技术领域
本发明属于水文预测技术领域,具体涉及一种场次暴雨径流量预测方法。
背景技术
降雨不仅是气候变化的主导因子,且是产流产沙的动力因子,使得流域径流量对降雨变化极为敏感,近年来,全国各地汛期时极端降雨也即暴雨频次和总量都呈现增加态势,暴雨危害和破坏性比普通降水更为明显,黄土高原的年径流量和年侵蚀量基本由少数几场暴雨引起,且暴雨事件在流域内易引起大水和大沙,易导致水库、淤地坝溃坝等风险,城镇地区的极端降水会导致城市防洪压力进一步加大,同时,由于城市的排水系统的设计标准偏低,导致我国城市内涝现象日益严重。
而现有技术中普遍将暴雨情形下的径流量预测与普通降雨特性认为一致,没有考虑到暴雨情形和普通降雨情形两者间性质的不同。
因此,如何准确地对场次暴雨产生的径流量进行预测,是本领域技术人员有待解决的技术问题。
发明内容
本发明的目的是为了解决现有技术无法对场次暴雨所产生的径流量进行准确地预测,提出了一种场次暴雨径流量预测方法。
本发明的技术方案为:一种场次暴雨径流量预测方法,包括以下步骤:
S1、获取目标流域内的历史数据,所述历史数据具体为所述目标流域内所有雨量站和所有水文站中记载的降雨数据和洪水摘录表数据,并对所述历史数据进行预处理得到预处理数据;
S2、基于曼肯德尔法划分出所述预处理数据的现状期并得到现状期数据;
S3、通过R/S分析法在所述现状期数据中确定出预测数据;
S4、基于所述预测数据对场次暴雨径流量进行预测。
进一步地,所述步骤S1中的预处理具体包括以下步骤:
S11、将所述降雨数据和洪水摘录表数据统一时间步长;
S12、根据统一时间步长后的各站点对应的降雨数据确定所述目标流域中场次降雨对应的场次降雨量;
S13、根据统一时间步长后的各站点对应的洪水摘录表确定所述目标流域中场次降雨对应的场次洪量;
S14、将大于50mm的场次降雨量对应的场次降雨和场次洪量分别组成场次降雨序列和场次洪量序列,并将所述场次降雨序列和所述场次洪量序列组合为所述预处理数据。
进一步地,所述步骤S12中确定所述目标流域中场次降雨对应的场次降雨量具体通过泰森多边形法进行确定,所述泰森多边形法如下公式所示:
式中,P为场次降雨量,Pi为流域内第i个雨量站雨量,αi为流域内第i个雨量站的面积权重系数,n为流域内雨量站个数。
进一步地,所述步骤S2具体包括以下分步骤:
S21、通过曼-肯德尔法判断所述预处理数据是否通过显著性检验,若是,则执行步骤S22,若否,则执行步骤S23;
S22、将所述预处理数据进行突变点检验,并以突变点检验结果为依据划分所述现状期以及所述现状期数据;
S23、将所述预处理数据作为所述现状期,并将所述预处理数据作为现状期数据。
进一步地,所述步骤S3具体为通过R/S分析法确定出所述现状期数据的赫斯特指数,并根据所述赫斯特指数确定出预测数据。
进一步地,所述步骤S4具体包括以下分步骤:
S41、基于所述预测数据,通过场次降雨量拟合法、降雨因子组合拟合法、场次降雨量分时段雨量拟合法和上包络法建立四种径流量预测计算式;
S42、基于所述场次降雨量拟合法、降雨因子组合拟合法、场次降雨量分时段雨量拟合法和上包络法建立的四种径流量预测计算式分别对场次暴雨径流量进行预测,并得到四个预测结果;
S43、基于所述四个预测结果确定所述场次暴雨径流量的预测区间,并将所述预测区间作为最终预测结果。
进一步地,所述场次降雨量分时段雨量拟合法建立的径流量预测计算式具体如下:
W=K1P1+K2P2-1+K3P3-2+K4P4-3+K5P5-4+K6P6-5+K7PP-6+c;
式中,W为径流量预测值,P1、P2-1、P3-2、P4-3、P5-4、P6-5、PP-6为各分时段雨量,K1、K2、K3、K4、K5、K6、K7、c为方程拟合系数。
与现有技术相比,本发明具备以下有益效果:
(1)本发明通过收集目标流域内各个水文气象站点的降雨径流数据;其次,根据曼-肯德尔法次相关检验法和R/S分析法分别对流域内场次极端暴雨量和次洪量进行趋势检验和持续性检验;最后,利用次降雨量拟合法、降雨因子组合拟合法、次降雨分时段雨量拟合法、上包络线法4种方法对目标流域内场次暴雨径流量进行预测,本发明能够较为简单快速的对单次暴雨洪水进行预测,为政府决策部门防汛减灾、预测预警、调洪调度工作的开展以及流域水资源合理开发利用提供理论支撑和技术支持。
(2)本发明得到四种预测公式,不同降雨因子和不同的降雨预测效果不同,通过不同公式预测径流量的区间,达到预测结果更加精准。
附图说明
图1所示为本发明实施例提供的一种场次暴雨径流量预测方法的流程示意图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
本申请提出了一种场次暴雨径流量预测方法,如图1所示为一种场次暴雨径流量预测方法的流程示意图,该方法包括以下步骤:
步骤S1、获取目标流域内的历史数据,所述历史数据具体为所述目标流域内所有雨量站和所有水文站中记载的降雨数据和洪水摘录表数据,并对所述历史数据进行预处理得到预处理数据。
在本申请实施例中,所述步骤S1中的预处理具体包括以下步骤:
S11、将所述降雨数据和洪水摘录表数据统一时间步长;
S12、根据统一时间步长后的各站点对应的降雨数据确定所述目标流域中场次降雨对应的场次降雨量;
S13、根据统一时间步长后的各站点对应的洪水摘录表确定所述目标流域中场次降雨对应的场次洪量;
S14、将大于50mm的场次降雨量对应的场次降雨和场次洪量分别组成场次降雨序列和场次洪量序列,并将所述场次降雨序列和所述场次洪量序列组合为所述预处理数据。
在具体的应用场景中,统一时间步长具体是等时间间隔数据,也即一小时,场次洪量的计算公式如下:
式中,W为次洪水量,亿m3;t为时段,h;i为次洪过程时段次序,n为次洪过程时段划分数量;Q为流量,m3/s。
所述步骤S12中确定所述目标流域中场次降雨对应的场次降雨量具体通过泰森多边形法进行确定,所述泰森多边形法如下公式所示:
式中,P为场次降雨量,Pi为流域内第i个雨量站雨量,αi为流域内第i个雨量站的面积权重系数,n为流域内雨量站个数。
步骤S2、基于曼-肯德尔法划分出所述预处理数据的现状期并得到现状期数据。
在本申请实施例中,所述步骤S2具体包括以下分步骤:
S21、通过曼-肯德尔法判断所述预处理数据是否通过显著性检验,若是,则执行步骤S22,若否,则执行步骤S23;
S22、将所述预处理数据进行突变点检验,并以突变点检验结果为依据划分所述现状期以及所述现状期数据;
S23、将所述预处理数据作为所述现状期,并将所述预处理数据作为现状期数据。
在具体应用场景中,曼-肯德尔法具体实施步骤如下:
对于一个具有n个样本量的时间序列X,首先构造一秩序列:
式中:Sk是第i时刻数值大于j时刻数值的累计个数,Xi,Xj分别是第i时刻和j时刻的样本值。
定义统计量:
以α=0.05为给定的显著性水平值,即U0.05=±1.96。当UFi<0时,表明序列呈下降趋势;当UFi>0时,表明序列呈上升趋势。一般情况下,若UFi>U0.05时,认为序列存在明显的上升或下降趋势变化。
对通过显著性检验的序列进行突变点检验,以突变点检验结果为依据进行时期的划分。若未通过趋势的显著性检验,则整个序列皆为现状期;认为序列受气候变化和人类活动的影响小,整个研究时期内未发生显著变化。
若UFi>U0.05,则表明序列存在明显的趋势变化,将时间序列X按逆序排列,再按曼-肯德尔法计算所得的UFi值定义为UBi,如果UBi和UFi两条曲线出现交点且交点在临界直线之间,那么交点对应的时刻就是突变时刻。
步骤S3、通过R/S分析法在所述现状期数据中确定出预测数据。
在本申请实施例中,所述步骤S3具体为通过R/S分析法确定出所述现状期数据的赫斯特指数,并根据所述赫斯特指数确定出预测数据。
在具体应用场景中,对现状期下的场次暴雨量序列和径流量时间序{ζ(t)}t=1,2,…,n,对于任意的正整数τ≥1,定义如下的统计量:
均值序列:
累计离差:
极差:
R(τ)=maxX(t,τ)-minX(t,τ)τ=1,2…
标准差:
如果R/S=R(τ)/S(τ)存在以下关系:
R/S∝τH
则序列存在Hurst现象也即赫斯特指数,两边取对数,即Log(R/S)=Logc+Hlogτ,c为某一常数值,Log(R/S)~Logτ图形中,直线斜率H即为Hurst指数。0.5<Hurst<1,表明未来趋势与现状下一致;若0<Hurst<0.5,表明未来趋势与现状趋势相反。
以R/S分析法检测结果判断数据是否具有预测性。若现状期数据通过R/S持续性检验,即0.5<Hurst<1,表明未来暴雨径流趋势与现状期趋势一致,序列具有预测性。若0<Hurst<0.5,表明未来暴雨径流与现状期趋势相反,序列不具有预测性。
步骤S4、基于所述预测数据对场次暴雨径流量进行预测。
在本申请实施例中,所述步骤S4具体包括以下分步骤:
S41、基于所述预测数据,通过场次降雨量拟合法、降雨因子组合拟合法、场次降雨量分时段雨量拟合法和上包络法建立四种径流量预测计算式;
S42、基于所述场次降雨量拟合法、降雨因子组合拟合法、场次降雨量分时段雨量拟合法和上包络法建立的四种径流量预测计算式分别对场次暴雨径流量进行预测,并得到四个预测结果;
S43、基于所述四个预测结果确定所述场次暴雨径流量的预测区间,并将所述预测区间作为最终预测结果。
在本申请实施例中,所述场次降雨量分时段雨量拟合法建立的径流量预测计算式具体如下:
W=K1P1+K2P2-1+K3P3-2+K4P4-3+K5P5-4+K6P6-5+K7PP-6+c;
式中,W为径流量预测值,P1、P2-1、P3-2、P4-3、P5-4、P6-5、PP-6为各分时段雨量,K1、K2、K3、K4、K5、K6、K7、c为方程拟合系数。
上式中,每场降雨过程中7个分时段最大雨量,即最大1h降雨量(P1)、最大2h降雨量减去最大1h降雨量(P2-1)、…、最大6h降雨量减去最大5h降雨量(P6-5)、最大6h降雨量之外的降雨量(Pp-6)。
在具体应用场景中,次降雨量拟合法为以场次降雨量为横坐标,场次径流量为纵坐标,拟合现状期下的径流量预测计算式。
降雨因子组合拟合法为影响场次降雨产流的因子包括降雨量及雨强,为同时考虑两者对径流量的影响,将场次降雨量与平均雨强的乘积PI作为自变量,径流量作为因变量,拟合现状期下降雨因子与径流量的计算式。
上包络线法为尽可能考虑复杂不利条件的影响,选取场次降雨量与径流量关系图的上边缘点据,点绘现状期的上包络线并建立径流量预测计算式。
需要说明的是,场次降雨量拟合法、降雨因子组合拟合法和上包络法建立的径流量预测计算式,是根据实测点数据,定义自变量和因变量后,放到统计软件中进行非线性拟合,最终以决定性系数为评判标准择优选择径流量预测计算式,例:次降雨量拟合法为以场次暴雨降雨量为横坐标,场次径流量为纵坐标,输入序列后,显示二次项计算式拟合较好,所以最后计算式为W1=8×10-5P2-0.003P+0.13,不同的区域计算式不同,也有可能是呈指数关系,三次关系等,主要抓住不同的降雨因子作为自变量,这是由本领域技术人员根据实际情况进行选择实施的。
为验证本发明所提出的方法能够对流域内暴雨情景下的径流量进行预测,选取汾河静乐站控制流域作为验证流域,流域内包括怀道、宁化堡和海子背等雨量站,场次径流量利用汾河静乐站记录数据。
数据选取的时间跨度为1971~2014年汛期场次产流降雨的相关数据,共89场数据,将其整理为时间间隔1h的降雨、径流过程数据,统计计算每场洪水的降雨量、洪量、降雨历时、平均雨强、各分时段雨量。
利用Mann-Kendall非参数秩次相关检验法对场次暴雨量和径流量序列进行了趋势性检验和突变性检验。场次极端暴雨量和对应的次洪量的UFi分别为0.61和1.74(表1),序列有上升趋势,但都没有通过显著性检验,表明其上升趋势不明显。无突变点,整个序列作为现状期。
表1流域场次暴雨洪水特征变化趋势检验
利用R/S分析法对场次暴雨量和径流量序列径流持续性分析。场次暴雨量和径流量序列的Hurst指数分别为0.61和0.71,均大于0.5,表明流域内的降雨径流序列在未来一段时间内仍持续为不显著增长趋势。可以用现状降雨径流序列预测未来场次暴雨下的可能径流量。
按照步骤S4建立次降雨量拟合法、降雨因子组合拟合法、次降雨分时段雨量拟合法、上包络线法对应的四种径流量预测计算式,见表2。并利用水文频率分析法和历史极端降雨组合法设计极端降雨情景,频率达到百年一遇和三百年一遇并带入计算式,与实际的同等量级场次降雨洪量对比发现,预测效果较好。
表2径流量预测计算式
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (6)
1.一种场次暴雨径流量预测方法,其特征在于,包括以下步骤:
S1、获取目标流域内的历史数据,所述历史数据具体为所述目标流域内所有雨量站和所有水文站中记载的降雨数据和洪水摘录表数据,并对所述历史数据进行预处理得到预处理数据;
S1中的预处理具体包括以下步骤:
S11、将所述降雨数据和洪水摘录表数据统一时间步长;
S12、根据统一时间步长后的各站点对应的降雨数据确定所述目标流域中场次降雨对应的场次降雨量;
S13、根据统一时间步长后的各站点对应的洪水摘录表确定所述目标流域中场次降雨对应的场次洪量;
S14、将大于50mm的场次降雨量对应的场次降雨和场次洪量分别组成场次降雨序列和场次洪量序列,并将所述场次降雨序列和所述场次洪量序列组合为所述预处理数据;
S2、基于曼肯德尔法划分出所述预处理数据的现状期并得到现状期数据;
S3、通过R/S分析法在所述现状期数据中确定出预测数据;
S4、基于所述预测数据对场次暴雨径流量进行预测。
3.如权利要求1所述的场次暴雨径流量预测方法,其特征在于,所述步骤S2具体包括以下分步骤:
S21、通过曼-肯德尔法判断所述预处理数据是否通过显著性检验,若是,则执行步骤S22,若否,则执行步骤S23;
S22、将所述预处理数据进行突变点检验,并以突变点检验结果为依据划分所述现状期以及所述现状期数据;
S23、将所述预处理数据作为所述现状期,并将所述预处理数据作为现状期数据。
4.如权利要求1所述的场次暴雨径流量预测方法,其特征在于,所述步骤S3具体为通过R/S分析法确定出所述现状期数据的赫斯特指数,并根据所述赫斯特指数确定出预测数据。
5.如权利要求1所述的场次暴雨径流量预测方法,其特征在于,所述步骤S4具体包括以下分步骤:
S41、基于所述预测数据,通过场次降雨量拟合法、降雨因子组合拟合法、场次降雨量分时段雨量拟合法和上包络法建立四种径流量预测计算式;
S42、基于所述场次降雨量拟合法、降雨因子组合拟合法、场次降雨量分时段雨量拟合法和上包络法建立的四种径流量预测计算式分别对场次暴雨径流量进行预测,并得到四个预测结果;
S43、基于所述四个预测结果确定所述场次暴雨径流量的预测区间,并将所述预测区间作为最终预测结果。
6.如权利要求5所述的场次暴雨径流量预测方法,其特征在于,所述场次降雨量分时段雨量拟合法建立的径流量预测计算式具体如下:
W=K1P1+K2P2-1+K3P3-2+K4P4-3+K5P5-4+K6P6-5+K7PP-6+c;
式中,W为径流量预测值,P1、P2-1、P3-2、P4-3、P5-4、P6-5、PP-6为各分时段雨量,K1、K2、K3、K4、K5、K6、K7、c为方程拟合系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110868146.8A CN113589404B (zh) | 2021-07-30 | 2021-07-30 | 一种场次暴雨径流量预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110868146.8A CN113589404B (zh) | 2021-07-30 | 2021-07-30 | 一种场次暴雨径流量预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113589404A CN113589404A (zh) | 2021-11-02 |
CN113589404B true CN113589404B (zh) | 2023-02-03 |
Family
ID=78252359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110868146.8A Active CN113589404B (zh) | 2021-07-30 | 2021-07-30 | 一种场次暴雨径流量预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113589404B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114023049B (zh) * | 2021-11-25 | 2022-08-16 | 西安理工大学 | 一种山洪灾害预警指标检验复核方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288229A (zh) * | 2011-05-11 | 2011-12-21 | 中国水利水电科学研究院 | 径流量模拟预测方法 |
CN102495937A (zh) * | 2011-10-18 | 2012-06-13 | 南京信息工程大学 | 一种基于时间序列的预测方法 |
CN103713336A (zh) * | 2013-12-24 | 2014-04-09 | 广西壮族自治区气象服务中心 | 基于gis分区的水电站流域面雨量气象预报方法 |
CN105046051A (zh) * | 2015-06-19 | 2015-11-11 | 长安大学 | 一种实际环境冻融作用次数预估方法 |
CN107220496A (zh) * | 2017-05-26 | 2017-09-29 | 上海市气象灾害防御技术中心 | 一种城市暴雨内涝评估建模方法 |
CN108491682A (zh) * | 2014-07-06 | 2018-09-04 | 韦雪 | 降雨径流预报系统 |
CN109886461A (zh) * | 2019-01-18 | 2019-06-14 | 昆仑(重庆)河湖生态研究院(有限合伙) | 一种径流预报方法及装置 |
CN112784479A (zh) * | 2021-01-12 | 2021-05-11 | 河海大学 | 一种洪水流量预测的方法 |
KR20210077228A (ko) * | 2019-12-17 | 2021-06-25 | (주)지오시스템리서치 | 침수 범람 예측 시스템 |
-
2021
- 2021-07-30 CN CN202110868146.8A patent/CN113589404B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288229A (zh) * | 2011-05-11 | 2011-12-21 | 中国水利水电科学研究院 | 径流量模拟预测方法 |
CN102495937A (zh) * | 2011-10-18 | 2012-06-13 | 南京信息工程大学 | 一种基于时间序列的预测方法 |
CN103713336A (zh) * | 2013-12-24 | 2014-04-09 | 广西壮族自治区气象服务中心 | 基于gis分区的水电站流域面雨量气象预报方法 |
CN108491682A (zh) * | 2014-07-06 | 2018-09-04 | 韦雪 | 降雨径流预报系统 |
CN105046051A (zh) * | 2015-06-19 | 2015-11-11 | 长安大学 | 一种实际环境冻融作用次数预估方法 |
CN107220496A (zh) * | 2017-05-26 | 2017-09-29 | 上海市气象灾害防御技术中心 | 一种城市暴雨内涝评估建模方法 |
CN109886461A (zh) * | 2019-01-18 | 2019-06-14 | 昆仑(重庆)河湖生态研究院(有限合伙) | 一种径流预报方法及装置 |
KR20210077228A (ko) * | 2019-12-17 | 2021-06-25 | (주)지오시스템리서치 | 침수 범람 예측 시스템 |
CN112784479A (zh) * | 2021-01-12 | 2021-05-11 | 河海大学 | 一种洪水流量预测的方法 |
Non-Patent Citations (7)
Title |
---|
R/S和Mann—Kendall法综合分析水文时间序列未来的趋势特征;于延胜;《水资源与水工程学报》;20080630;全文 * |
城市场次降雨径流污染负荷快速估算方法;何胜男;《湖泊科学》;20210120;全文 * |
基于Mann-Kendall和R/S法的水文序列变化趋势分析——苏州市为例;杨金艳;《水利水电技术》;20170220;摘要,正文第28-30页 * |
基于R_S和M-K分析的高良涧闸年径流量变化趋势研究;黄璐;《江苏水利》;20191130;全文 * |
广西极端降水时空变化规律研究;孙桂凯;《广西大学学报( 自然科学版)》;20210430;全文 * |
汉江流域安康段降雨径流特征分析及预测;刘易文;《人民珠江》;20210624;摘要,正文第68页 * |
镇汀市降雨场次划分与暴雨特性分析;皋云;《水资源保护》;20210531;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113589404A (zh) | 2021-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111310968B (zh) | 一种基于互信息的lstm神经网络循环水文预报方法 | |
Sedki et al. | Simulation-optimization modeling for sustainable groundwater development: a Moroccan coastal aquifer case study | |
KR20170005553A (ko) | 지능형 정보제공을 위한 홍수, 가뭄정보의 평가 및 예측 기법 개발 | |
CN112686426B (zh) | 一种基于水电站流域关键点的来水量预警方法及系统 | |
CN112528563B (zh) | 一种基于svm算法的城市内涝预警方法 | |
CN112036683B (zh) | 一种适用于未来气候变化情景下的水库防洪风险预估方法 | |
Li et al. | Real-time flood forecast using the coupling support vector machine and data assimilation method | |
CN113589404B (zh) | 一种场次暴雨径流量预测方法 | |
Rashid et al. | Statistical characteristics of rainfall in the Onkaparinga catchment in South Australia | |
CN114756817A (zh) | 一种基于Copula函数的复合洪涝灾害联合概率分析方法 | |
Jiang et al. | Reservoir flood season segmentation and optimal operation of flood-limiting water levels | |
CN109902266B (zh) | 一种基于Copula函数的河道流量演算方法 | |
Fang et al. | Non-identical models for seasonal flood frequency analysis | |
CN116167513A (zh) | 基于单变量优化dmca模型的流域洪水响应时间计算方法 | |
CN113506000B (zh) | 一种计算城市汇水区域暴雨溢流风险的方法 | |
CN115619045A (zh) | 一种河口潮差预报方法和系统 | |
CN114528761A (zh) | 平原水网圩区-圩外系统滞蓄关系优化方法及系统 | |
CN110807475B (zh) | 一种基于确定性系数的洪水分类识别预报方法 | |
Suljić et al. | Analysis of Time Oscillations of Water on Lake Modric as a Multi-Purpose Reservoir | |
Bacalhau et al. | Water supply reservoir operation in relation to climate variability: Pirapama river basin (Pernambuco-Brazil) | |
Zhang et al. | Streamflow pattern variations resulting from future climate change in middle Tianshan Mountains region in China | |
CN114970082B (zh) | 非一致性设计洪水估计方法 | |
CN117216468B (zh) | 基于深度学习及联合概率分布的湖泊生态缺水预警方法 | |
CN117150600B (zh) | 一种用于城市防内涝的屋顶蓄排水控制方法 | |
CN117494002A (zh) | 一种流域气象-水文旱涝遭遇特征评估方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |