CN111639810A - 一种基于防汛需求的降雨预报精度评估方法 - Google Patents
一种基于防汛需求的降雨预报精度评估方法 Download PDFInfo
- Publication number
- CN111639810A CN111639810A CN202010483351.8A CN202010483351A CN111639810A CN 111639810 A CN111639810 A CN 111639810A CN 202010483351 A CN202010483351 A CN 202010483351A CN 111639810 A CN111639810 A CN 111639810A
- Authority
- CN
- China
- Prior art keywords
- rainfall
- value
- forecast
- evaluation
- time
- 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.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 116
- 230000002265 prevention Effects 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 17
- 238000009825 accumulation Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000004927 fusion Effects 0.000 abstract description 4
- 230000000694 effects Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000002035 prolonged effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 241000735432 Hydrastis canadensis Species 0.000 description 1
- 241000219094 Vitaceae Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 235000005679 goldenseal Nutrition 0.000 description 1
- 235000021021 grapes Nutrition 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Environmental & Geological Engineering (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Environmental Sciences (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Ecology (AREA)
- Biodiversity & Conservation Biology (AREA)
- Atmospheric Sciences (AREA)
- Hydrology & Water Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Alarm Systems (AREA)
Abstract
本发明提供一种基于防汛需求的降雨预报精度评估方法,方法包括:第一步:网格化与面雨量计算;第二步:降雨预报量级评价指标;第三步:降雨预报空间分布评价指标;第四步:降雨预报时程分配评价指标;第五步:降雨预报综合评价指标。本发明将数值预报降雨与实测降雨值统一到标准网格下,分别进行量级、空间分布、时程分配的吻合程度评估,得到三个维度的分项指标,综合形成对预报降雨的综合评价指标。评价的结果可以为降雨预报模式的应用选取,以及多种预报模式的融合应用提供参考。
Description
技术领域
本发明属于环境预报技术领域,具体涉及一种基于防汛需求的降雨预报精度评估方法。
背景技术
降雨信息是洪水预报的重要依据,传统的洪水预报往往仅依赖实测降雨,近年来有很多文献尝试将数值降雨预报成果作为水文模型的输入,有效延长了洪水预报预见期(参见文献:[1]安莉娜,王文杰,郑敏生等.短期降雨预报信息应用在青山水库汛限水位动态控制的可行性分析[J].中国农村水利水电,2008,(6):55-56,64;[2]张卫国,范仲丽,钟伟,等.雷达回波外推方法在临近降雨预报中的应用[J].中国农村水利水电,2018,(9):69-73,120;[3]Verbunt M,Walser A,Gurtz J,et al.Probabilistic Flood Forecastingwith a Limited-Area Ensemble Prediction System:Selected Case Studies[J].Journal of Hydrometeorology,2007,8(4):897-909;[4]温娅惠,李致家,孙明坤,等.降雨输入对实时洪水预报精度与预见期的影响[J].湖泊科学,2019,31(1):39-51),可为水库、闸泵调度等抢险救灾工作赢得更多时间。当前国内外降雨数值预报模式繁多,例如美国GFS(美国国家环境预报中心研发的一种气象数值预报模式Global Forecast System)、欧洲ECWMF(欧洲中期天气预报中心研发的一种气象数值预报模式European Centre forMedium-Range Weather Forecasts)、我国中央气象台GRAPES_MESO(中国气象局数值预报中心研发的中国及周边区域数值预报产品(Global/Regional Assimilation PredictionSystem))、上海气科所SMS-WARMS(中国气象局上海气象科学研究所研发的华东区域中尺度数值预报系统)等。
以上各数值预报模式在不同预报时长、不同降雨等级、不同地区(流域)的预报精度均有所差别(参见文献:[5]刘静,叶金印,张晓红,等.淮河流域汛期面雨量多模式预报检验评估[J].暴雨灾害,2014,(1):58-64、[6]张宁娜,黄阁,吴曼丽,等.2010年国内外3种数值预报在东北地区的预报检验[J].气象与环境学报,2012,28(2):28-33),洪水预报选用的降雨数值预报成果的精度,很大程度上决定了的水文模型预报精度,也会直接影响到水利工程的调度等指挥决策工作。因此,开展降雨预报成果精度的评价研究,对洪水预报工作具有重要意义。
国内外开展了对降雨预报精度的大量研究工作,如刘硕等([7]刘硕,王国利,张琳.TlGGE降雨信息在柴河流域洪水预报中可利用性评估[J].水文,2018,38(5):17-22,84)采用TS评分和偏差Bias(TS评分:即临界成功指数,代表空报率和漏报率的综合效应。偏差Bias:即预报偏差,表示预报事件发生的次数与观测事件发生的次数的比率)来分析柴河流域定量降雨预报成果的准确率、空报率和漏报率;曹越等([8]曹越,赵琳娜,巩远发,等.ECMWF高分辨率模式降水预报能力评估与误差分析[J].暴雨灾害,2019,38(3):249-258)对ECMWF高分辨率模式与我国多个雨量测站的24h降水预报的准确度、集中度、相关度进行了评估。
现有技术大都关注流域、区域或站点的预报降雨的总量误差、场次命中率等指标,没有考虑降雨预报的空间分布、时程分配误差,但这两个要素直接影响到防汛应用场景时洪水预报成果的洪峰、峰现时间判断,而其结果也影响到水利工程调度时机选择、流域水利工程联合调度效益的发挥。实际上,支撑防汛应用场景需要的对降雨预报的综合考虑总量、空间分布、时程分配三者预报准确度的评价。将数值降雨预报成果作为水文模型的输入条件,能有效延长洪水预报预见期,为防汛工作赢得更多时间。但面对国内外多种降雨数值预报模式成果,其预报精度是选取数值预报模式的重要指标。而现有的TS评分、偏差Bias等降雨预报指标,不能适用于洪水预报等防汛应用场景需求。
为此,本发明申请提供了一种基于防汛需求的降雨预报精度评估方法,本发明综合了降雨时程分配、空间分布、总量量级三个维度的降雨预报精度评价,其评估的结果可以为防汛场景下的降雨预报模式的防汛应用选取,以及多种降雨数值预报模式的融合应用提供依据。
发明内容
针对现有技术存在的问题为此,本发明提出了一种基于防汛需求的降雨预报精度评估方法,该方法是一种综合量级、空间分布、时程分配三个维度的降雨预报精度评价方法:在量级上,建立以相对误差为基础的量级评价指标;在空间分布上,引入容差区间的概念,进行逐网格赋分,综合全流域网格后得出空间分布评价指标;在时程上,采用确定性系数评价时间过程的匹配程度,得到时程分配评价指标。其评价的结果可以为降雨预报模式的应用选取,以及多种预报模式的融合应用提供参考。所述方法包括如下步骤:
第一步:网格化与面雨量计算
将数值预报降雨值与实测降雨值统一到统一尺寸的标准网格下;按照统一设定的标准网格尺寸,逐时段采用克里金插值得到网格降雨值,将流域内各网格预报值取算术平均即得到流域面雨量预报值;实测数据主要来源于雨量测站,同样逐时段采用克里金插值得到网格降雨值,将流域内各网格降雨值取算术平均即得到流域面雨量实测值;由此,得到各网格以及全流域的逐时段实测、预报降雨量,为后续步骤提供数据基础;
第二步:降雨预报量级评价指标
对选定的流域总雨量进行量级精度评价,考虑到与时程分配、空间分布指标协调性,选取一个无量纲的值,引入相对误差值项,并按照指标的正向性考虑,取累积总降雨量的量级评价指标为:
当Esum,all<0时,令
Esum,all=0 (2)
Esum,all作为量级评价指标,值在[0,1]区间,值越大,评价指标越优;在实测降雨量很小时,会引起的表述失真;根据避免误差危害的若干原则,结合精度评价结果的防洪应用场景,当实测降雨Psum,all≤Tmm时,式(1)中除数取常数T,即:
当Esum,all<0时,令
Esum,all=0 (4)
除了上述全时段累积总降雨之外,在防汛实践中,还关注考虑最大1h、3h、6h、12h、24h等特征时段累积降雨量,为此,相应的针对最大i小时的累积雨量量级评估公式如下:
当Esum,ih<0时,令
Esum,ih=0 (6)
最后,将全时段、各个特征时段雨量量级评价指标综合,最终的量级评价指标为
式中:cih为实测最大i小时降雨量级评价值所占权重;call为全时段量级评价值所占权重,权重根据流域产汇流特性取值,或按平均权重处理;
第三步:降雨预报空间分布评价指标
(1)引入容差区间的概念
对暴雨以上等级划分做适当加密,形成为两档容差区间表:
1)评估时长≤12h时,设置容差区间一;
2)评估时长>12h时,设置容差区间二;
(2)制定判定规则
若编号为k网格的全时段累积总降雨的预报值与实测值,落在同一容差区间的,则给网格全时段累积总降雨的空间相符性指标Boolk,all赋值1,即Boolk,all=1,否则赋值为0,即Boolk,all=0;
同时关注最大1h、3h、6h、12h、24h等特征时段累积降雨量的空间分布影响:相应的针对最大i小时的累积雨量按上述规则考虑,即若编号为k的网格最大i小时降雨累计量的预报值与实测值,落在同一容差区间的,令该网格最大i小时累积降雨的空间相符性指标赋值Boolk,ih=1,否则Boolk,ih=0;
(3)得到流域空间评价指标
将第k个网格的全时段、各个特征时段的空间相符性指标加权,得到该网格的空间相符性指标为
再将流域内所有网格的空间相符性指标取算术平均,得到流域的空间分布评价指标Edis,数值在[0,1]区间,值越大,指标越优:
式中,n为流域内的总网格数;
第四步:降雨预报时程分配评价指标
定义洪水预报的常用指标“确定性系数”作为时程评价指标Epro,数值在[0,1]区间,值越大,指标越优:
第五步:降雨预报综合评价指标
将第二至第四步中得到的量级评价指标、空间分布评价指标和时程分配评价指标进行加权计算,得到降雨的综合评价指标E,其数值在[0,1]区间,值越大,指标越优;
E=Esum*Csum+Edis*Cdis+Epro*Cpro (11)
式中,Csum为量级权重系数,Cdis为空间分布权重系数,Cpro为时程分配权重系数;这三个权重系数根据对洪水洪量、洪峰等关注侧重不同而定。
优选的,在所述第二步中,当实测降雨Psum,all≤10mm时,式(1)中除数取常数T取为10,即
关注考虑最大1h、3h、6h、12h、24h等特征时段累积降雨量,式(1)中除数取常数T取为10,相应的针对最大i小时的累积雨量量级评估公式如下:
优选的,在所述第三步中,
1)评估时长≤12h时,容差区间一取为:
[0,5),[5,15),[15,30),[30,50),[50,70),[70,140),[140,+∞);
2)评估时长>12h时,容差区间二取为:
[0,10),[10,25),[25,75),[75,100),[100,150),[150,200),[200,250),[250,+∞)。
优选的,在所述第五步中,在简化处理时,量级权重系数、空间分布权重系数、时程分配权重系数取值均为1/3。
与现有方法相比,本发明具有以下优势:
(1)从防汛的场景出发,关注的是降雨引起的产汇流效应,提出的综合考虑降雨量级、空间分布、时程分配的多维度评估指标方法。
(2)该方法是一种综合量级、空间分布、时程分配三个维度的降雨预报精度评价方法:在量级上,建立以相对误差为基础的量级评价指标;在空间分布上,引入容差区间的概念,进行逐网格赋分,综合全流域网格后得出空间分布评价指标;在时程上,采用确定性系数评价时间过程的匹配程度,得到时程分配评价指标。
(3)在量级评价指标的建立方法,对相对误差项提出一种处理方法:将除数替代为某一合理常量的方法,有效避免了当降雨量级很小时,计算出的相对误差值很大,夸大降雨预报偏差对防汛的影响的恶劣情况。
(4)本发明的评估的结果可以为防汛场景下的降雨预报模式的防汛应用选取。
(5)本发明的评估结果可以为多种降雨数值预报模式的融合应用提供依据。
附图说明
图1基于防汛需求的降雨预报精度评估方法总体流程示意图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。参见附图1所示,本发明将数值预报降雨与实测降雨值统一到标准网格下,分别进行量级、空间分布、时程分配的吻合程度评估,得到三个维度的分项指标,综合形成对预报降雨的综合评价指标,具体步骤如下:
第一步:网格化与面雨量计算
由于涉及空间分布维度的评价,需要将数值预报降雨值与实测降雨值统一到统一尺寸的标准网格下。各数值预报模式水平分辨率不一,本发明按照统一设定的标准网格尺寸,逐时段采用克里金插值得到网格降雨值,将流域内各网格预报值取算术平均即得到流域面雨量预报值。实测数据主要来源于雨量测站,同样逐时段采用克里金插值得到网格降雨值,将流域内各网格降雨值取算术平均即得到流域面雨量实测值。
由此,可以得到各网格以及全流域的逐时段实测、预报降雨量,为本发明以下的评价步骤提供数据基础。
第二步:降雨预报量级评价指标
对选定的流域总雨量进行量级精度评价,考虑到与下文的时程分配、空间分布指标协调性,需选取一个无量纲的值,本文引入相对误差值项,并按照指标的正向性考虑,拟取累积总降雨量的量级评价指标为
当Esum,all<0时,令
Esum,all=0 (2)
Esum,all作为量级评价指标,值在[0,1]区间,值越大,评价指标越优。但上式在实测降雨量很小时,会引起的表述失真,例如实测降雨0.1mm,预报降雨0.2mm,按上式计算值Esumall=0,直接的表述理解为量级预报指标到恶劣下限,而从本文立意的防洪应用场景来看降雨0.2mm与0.1mm引起的水文响应差别很小,此时Esum,all为接近1的数才为合理。根据避免误差危害的若干原则,结合精度评价结果的防洪应用场景,本文规定当实测降雨Psum,all≤10mm时,上式除数取常数10,即
当Esum,all<0时,令
Esum,all=0 (4)
除了上述全时段累积总降雨之外,在防汛实践中,往往还关注考虑最大1h、3h、6h、12h、24h等特征时段累积降雨量。为此,相应的针对最大i小时的累积雨量量级评估公式如下:
当Esum,ih<0时,令
Esum,ih=0 (6)
最后,将全时段、各个特征雨量量级评价指标综合,最终的量级评价指标为
式中:cih为实测最大i小时降雨量级评价值所占权重;call为全时段量级评价值所占权重。权重可根据流域产汇流特性取值,或按平均权重处理。
第三步:降雨预报空间分布评价指标
由于流域内降雨存在空间分布不均匀的特点,为了准确描述预报降雨与实测降雨空间分布的吻合度,本发明提出了一种空间分布评价方法,:
(1)引入容差区间的概念
参考雨量等级划分的相关标准([9]国家质量监督检验检疫总局,国家标准化管理委员会.降雨量等级:GB/T 28592-2012[S],2012),同时出于区分度的考虑,对暴雨以上等级划为做了适当加密,形成为两档容差区间表:
1)评估时长≤12h时,容差区间一取为[0,5),[5,15),[15,30),[30,50),[50,70),[70,140),[140,+∞);
2)评估时长>12h时,容差区间二取为[0,10),[10,25),[25,75),[75,100),[100,150),[150,200),[200,250),[250,+∞);
(2)制定判定规则
若编号为k网格的全时段累积总降雨的预报值与实测值,落在同一容差区间的,则给网格全时段累积总降雨的空间相符性指标Boolk,all赋值1,即Boolk,all=1,否则赋值为0,即Boolk,all=0。例如,评估时长24小时,编号217号的网格,预报累积雨量167mm,实测累积雨量180mm,同属[150,200)区间,则认为该网格的空间相符性较高,Bool214,all=1.
同样地,除了全时段累积总降雨之外,在防汛实践中,也会关注最大1h、3h、6h、12h、24h等特征时段累积降雨量的空间分布影响。为此,相应的针对最大i小时的累积雨量按上述规则考虑,即若编号为k的网格最大i小时降雨累计量的预报值与实测值,落在同一容差区间的,令该网格最大i小时累积降雨的空间相符性指标赋值Boolk,ih=1,否则Boolk,ih=0.
(3)得到流域空间评价指标
将第k个网格的全时段、各个特征时段的空间相符性指标加权,得到该网格的空间相符性指标为
再将流域内所有网格的空间相符性指标取算术平均,得到流域的空间分布评价指标Edis,数值在[0,1]区间,值越大,指标越优:
式中,n为流域内的总网格数。
第四步:降雨预报时程分配评价指标
为了评价降雨预报过程与实测过程之间的吻合程度,引入洪水预报的常用指标“确定性系数”,作为时程评价指标Epro,数值在[0,1]区间,值越大,指标越优:
第五步:降雨预报综合评价指标
将以上步骤中得到的量级评价指标、空间分布评价指标和时程分配评价指标进行加权计算,得到降雨的综合评价指标E。同样地,其数值也在[0,1]区间,值越大,指标越优。
E=Esum*Csum+Edis*Cdis+Epro*Cpro (11)
式中,Csum为量级权重系数,Cdis为空间分布权重系数,Cpro为时程分配权重系数。这三个权重系数可根据对洪水洪量、洪峰等关注侧重不同而定,简化处理时取值为1/3。
本发明的技术效果如下:
(1)从防汛的场景出发,关注的是降雨引起的产汇流效应,提出的综合考虑降雨量级、空间分布、时程分配的多维度评估指标方法。
(2)空间分布指标的建立方法,在空间分布上,引入容差区间的概念,进行逐网格赋分,综合全流域网格后得出空间分布评价指标,这样就对降雨空间分布的误差得出一个较为客观的评分;
(3)在量级评价指标的建立方法,对相对误差项提出一种处理方法:将除数替代为某一合理常量的方法,有效避免了当降雨量级很小时,计算出的相对误差值很大,夸大降雨预报偏差对防汛的影响的恶劣情况。
以上所述之实施例仅为本发明的较佳实施例,并非对本发明做任何形式上的限制。任何熟悉本领域的技术人员,在不脱离本发明技术方案范围情况下,利用上述揭示的技术内容对本发明技术方案作出更多可能的变动和润饰,或修改均为本发明的等效实施例。故凡未脱离本发明技术方案的内容,依据本发明之思路所作的等同等效变化,均应涵盖于本发明的保护范围内。
Claims (4)
1.一种基于防汛需求的降雨预报精度评估方法,其特征在于:所述方法包括如下步骤:
第一步:网格化与面雨量计算
将数值预报降雨值与实测降雨值统一到统一尺寸的标准网格下;按照统一设定的标准网格尺寸,逐时段采用克里金插值得到网格降雨值,将流域内各网格预报值取算术平均即得到流域面雨量预报值;实测数据主要来源于雨量测站,同样逐时段采用克里金插值得到网格降雨值,将流域内各网格降雨值取算术平均即得到流域面雨量实测值;由此,得到各网格以及全流域的逐时段实测、预报降雨量,为后续步骤提供数据基础;
第二步:降雨预报量级评价指标
对选定的流域总雨量进行量级精度评价,考虑到与时程分配、空间分布指标协调性,选取一个无量纲的值,引入相对误差值项,并按照指标的正向性考虑,取累积总降雨量的量级评价指标为:
当Esum,all<0时,令
Esum,all=0 (2)
Esum,all作为量级评价指标,值在[0,1]区间,值越大,评价指标越优;在实测降雨量很小时,会引起的表述失真;根据避免误差危害的若干原则,结合精度评价结果的防洪应用场景,当实测降雨Psum,all≤Tmm时,式(1)中除数取常数T,即:
当Esum,all<0时,令
Esum,all=0 (4)
除了上述全时段累积总降雨之外,在防汛实践中,还关注考虑最大1h、3h、6h、12h、24h等特征时段累积降雨量,为此,相应的针对最大i小时的累积雨量量级评估公式如下:
当Esum,ih<0时,令
Esum,ih=0 (6)
最后,将全时段、各个特征时段雨量量级评价指标综合,最终的量级评价指标为
式中:cih为实测最大i小时降雨量级评价值所占权重;call为全时段量级评价值所占权重,权重根据流域产汇流特性取值,或按平均权重处理;
第三步:降雨预报空间分布评价指标
(1)引入容差区间的概念
对暴雨以上等级划分做适当加密,形成为两档容差区间表:
1)评估时长≤12h时,设置容差区间一;
2)评估时长>12h时,设置容差区间二;
(2)制定判定规则
若编号为k网格的全时段累积总降雨的预报值与实测值,落在同一容差区间的,则给网格全时段累积总降雨的空间相符性指标Boolk,all赋值1,即Boolk,all=1,否则赋值为0,即Boolk,all=0;
同时关注最大1h、3h、6h、12h、24h等特征时段累积降雨量的空间分布影响:相应的针对最大i小时的累积雨量按上述规则考虑,即若编号为k的网格最大i小时降雨累计量的预报值与实测值,落在同一容差区间的,令该网格最大i小时累积降雨的空间相符性指标赋值Boolk,ih=1,否则Boolk,ih=0;
(3)得到流域空间评价指标
将第k个网格的全时段、各个特征时段的空间相符性指标加权,得到该网格的空间相符性指标为
再将流域内所有网格的空间相符性指标取算术平均,得到流域的空间分布评价指标Edis,数值在[0,1]区间,值越大,指标越优:
式中,n为流域内的总网格数;
第四步:降雨预报时程分配评价指标
定义洪水预报的常用指标“确定性系数”作为时程评价指标Epro,数值在[0,1]区间,值越大,指标越优:
第五步:降雨预报综合评价指标
将第二至第四步中得到的量级评价指标、空间分布评价指标和时程分配评价指标进行加权计算,得到降雨的综合评价指标E,其数值在[0,1]区间,值越大,指标越优;
E=Esum*Csum+Edis*Cdis+Epro*Cpro (11)
式中,Csum为量级权重系数,Cdis为空间分布权重系数,Cpro为时程分配权重系数;这三个权重系数根据对洪水洪量、洪峰等关注侧重不同而定。
3.根据权利要求1所述的一种基于防汛需求的降雨预报精度评估方法,其特征在于:在所述第三步中,
1)评估时长≤12h时,容差区间一取为:
[0,5),[5,15),[15,30),[30,50),[50,70),[70,140),[140,+∞);
2)评估时长>12h时,容差区间二取为:
[0,10),[10,25),[25,75),[75,100),[100,150),[150,200),[200,250),[250,+∞)。
4.根据权利要求1所述的一种基于防汛需求的降雨预报精度评估方法,其特征在于:在所述第五步中,在简化处理时,量级权重系数、空间分布权重系数、时程分配权重系数取值均为1/3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010483351.8A CN111639810B (zh) | 2020-06-01 | 2020-06-01 | 一种基于防汛需求的降雨预报精度评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010483351.8A CN111639810B (zh) | 2020-06-01 | 2020-06-01 | 一种基于防汛需求的降雨预报精度评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111639810A true CN111639810A (zh) | 2020-09-08 |
CN111639810B CN111639810B (zh) | 2023-05-09 |
Family
ID=72329717
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010483351.8A Active CN111639810B (zh) | 2020-06-01 | 2020-06-01 | 一种基于防汛需求的降雨预报精度评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111639810B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112614309A (zh) * | 2020-12-18 | 2021-04-06 | 中国地质环境监测院 | 一种地质灾害气象预警准确率的计算方法 |
CN112837507A (zh) * | 2021-02-08 | 2021-05-25 | 国家海洋局东海预报中心 | 一种海浪预警报质量评估方法及装置 |
CN113009595A (zh) * | 2021-02-09 | 2021-06-22 | 海南省气象科学研究所 | 一种基于车载传感器的降雨监测方法及相关装置 |
CN114239925A (zh) * | 2021-11-29 | 2022-03-25 | 中国长江电力股份有限公司 | 一种水电站入库洪水预报偏差分割和评估方法 |
CN116068676A (zh) * | 2023-03-07 | 2023-05-05 | 南京气象科技创新研究院 | 一种降水预报跨量级通用综合评价方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006127156A (ja) * | 2004-10-28 | 2006-05-18 | Japan Weather Association | 水害リスク評価方法及びシステム、流域データーベースとその作成プログラム、水害シミュレーション方法及びプログラム |
JP2009008651A (ja) * | 2007-05-31 | 2009-01-15 | Foundation Of River & Basin Integrated Communications Japan | 全国合成レーダ雨量を用いた分布型流出予測システム |
CN105590034A (zh) * | 2016-03-15 | 2016-05-18 | 上海市水务规划设计研究院 | 防汛减灾能力评价模型 |
CN105808948A (zh) * | 2016-03-08 | 2016-07-27 | 中国水利水电科学研究院 | 一种自动修正的多模式数值降雨集合预报方法 |
CN105844709A (zh) * | 2016-03-25 | 2016-08-10 | 中国水利水电科学研究院 | 复杂河道地形流域洪水演进虚拟仿真的淹没线追踪方法 |
CN106529739A (zh) * | 2016-11-30 | 2017-03-22 | 福建四创软件有限公司 | 基于网格化的洪水预测方法 |
CN106884405A (zh) * | 2017-03-08 | 2017-06-23 | 中国水利水电科学研究院 | 一种无资料地区溃堤型山洪灾害分析评价方法 |
CN107179568A (zh) * | 2017-05-23 | 2017-09-19 | 贵州东方世纪科技股份有限公司 | 一种降水预报精度评定方法 |
CN108416468A (zh) * | 2018-02-08 | 2018-08-17 | 河海大学 | 一种洪水超前预警预报方法 |
CN108446436A (zh) * | 2018-02-08 | 2018-08-24 | 广州地理研究所 | 暴雨洪水非线性模型雨水损失参数的空间分布预警方法 |
CN109829584A (zh) * | 2019-02-02 | 2019-05-31 | 中国水利水电科学研究院 | 一种山洪风险动态评估方法 |
CN109871988A (zh) * | 2019-01-28 | 2019-06-11 | 河海大学 | 一种洪水预报预警精度分析方法 |
CN110009002A (zh) * | 2019-03-13 | 2019-07-12 | 河海大学 | 基于多维降雨特征空间最优决策的山洪快速预警预报方法 |
-
2020
- 2020-06-01 CN CN202010483351.8A patent/CN111639810B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006127156A (ja) * | 2004-10-28 | 2006-05-18 | Japan Weather Association | 水害リスク評価方法及びシステム、流域データーベースとその作成プログラム、水害シミュレーション方法及びプログラム |
JP2009008651A (ja) * | 2007-05-31 | 2009-01-15 | Foundation Of River & Basin Integrated Communications Japan | 全国合成レーダ雨量を用いた分布型流出予測システム |
CN105808948A (zh) * | 2016-03-08 | 2016-07-27 | 中国水利水电科学研究院 | 一种自动修正的多模式数值降雨集合预报方法 |
CN105590034A (zh) * | 2016-03-15 | 2016-05-18 | 上海市水务规划设计研究院 | 防汛减灾能力评价模型 |
CN105844709A (zh) * | 2016-03-25 | 2016-08-10 | 中国水利水电科学研究院 | 复杂河道地形流域洪水演进虚拟仿真的淹没线追踪方法 |
CN106529739A (zh) * | 2016-11-30 | 2017-03-22 | 福建四创软件有限公司 | 基于网格化的洪水预测方法 |
CN106884405A (zh) * | 2017-03-08 | 2017-06-23 | 中国水利水电科学研究院 | 一种无资料地区溃堤型山洪灾害分析评价方法 |
CN107179568A (zh) * | 2017-05-23 | 2017-09-19 | 贵州东方世纪科技股份有限公司 | 一种降水预报精度评定方法 |
CN108416468A (zh) * | 2018-02-08 | 2018-08-17 | 河海大学 | 一种洪水超前预警预报方法 |
CN108446436A (zh) * | 2018-02-08 | 2018-08-24 | 广州地理研究所 | 暴雨洪水非线性模型雨水损失参数的空间分布预警方法 |
CN109871988A (zh) * | 2019-01-28 | 2019-06-11 | 河海大学 | 一种洪水预报预警精度分析方法 |
CN109829584A (zh) * | 2019-02-02 | 2019-05-31 | 中国水利水电科学研究院 | 一种山洪风险动态评估方法 |
CN110009002A (zh) * | 2019-03-13 | 2019-07-12 | 河海大学 | 基于多维降雨特征空间最优决策的山洪快速预警预报方法 |
Non-Patent Citations (1)
Title |
---|
王洋;刘佳;于福亮;李传哲;田济扬;邱庆泰;: "基于数据同化的降雨数值空间分布模拟研究" * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112614309A (zh) * | 2020-12-18 | 2021-04-06 | 中国地质环境监测院 | 一种地质灾害气象预警准确率的计算方法 |
CN112614309B (zh) * | 2020-12-18 | 2021-10-01 | 中国地质环境监测院 | 一种地质灾害气象预警准确率的计算方法 |
CN112837507A (zh) * | 2021-02-08 | 2021-05-25 | 国家海洋局东海预报中心 | 一种海浪预警报质量评估方法及装置 |
CN113009595A (zh) * | 2021-02-09 | 2021-06-22 | 海南省气象科学研究所 | 一种基于车载传感器的降雨监测方法及相关装置 |
CN114239925A (zh) * | 2021-11-29 | 2022-03-25 | 中国长江电力股份有限公司 | 一种水电站入库洪水预报偏差分割和评估方法 |
CN116068676A (zh) * | 2023-03-07 | 2023-05-05 | 南京气象科技创新研究院 | 一种降水预报跨量级通用综合评价方法 |
CN116068676B (zh) * | 2023-03-07 | 2023-06-02 | 南京气象科技创新研究院 | 一种降水预报跨量级通用综合评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111639810B (zh) | 2023-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111639810A (zh) | 一种基于防汛需求的降雨预报精度评估方法 | |
Harrigan et al. | Designation and trend analysis of the updated UK Benchmark Network of river flow stations: The UKBN2 dataset | |
Keller et al. | CEH-GEAR: 1 km resolution daily and monthly areal rainfall estimates for the UK for hydrological and other applications | |
Yin et al. | Rainfall erosivity estimation based on rainfall data collected over a range of temporal resolutions | |
Rhoades et al. | Adjustment of temperature and rainfall records for site changes | |
Seckin et al. | Flood frequency analysis of Turkey using L‐moments method | |
Melching et al. | A reliability estimation in modeling watershed runoff with uncertainties | |
CN106991278A (zh) | 一种集合降水预报与实时洪水概率预报的耦合方法 | |
Katipoğlu et al. | Comparison of meteorological indices for drought monitoring and evaluating: a case study from Euphrates basin, Turkey | |
Formetta et al. | Estimating the index flood with continuous hydrological models: an application in Great Britain | |
Pruski et al. | Improved regionalization of streamflow by use of the streamflow equivalent of precipitation as an explanatory variable | |
CN110334404A (zh) | 一种流域尺度骤发干旱精准识别方法 | |
Palmen et al. | Regional flood frequency for Queensland using the quantile regression technique | |
CN111915158A (zh) | 一种基于Flood Area模型的暴雨灾害天气风险评估方法、装置及设备 | |
Otache et al. | ARMA modelling of Benue River flow dynamics: comparative study of PAR model | |
Rangsiwanichpong et al. | Assessment of sediment yield in Thailand using revised universal soil loss equation and geographic information system techniques | |
Msaddek et al. | Hydrological modeling of upper OumErRabia basin (Morocco), comparative study of the event-based and continuous-process HEC-HMS model methods | |
CN115239156A (zh) | 城镇化指标对水系结构的影响警示方法 | |
Tsakiris et al. | Regional drought identification and assessment. Case study in Crete | |
Jayadeera et al. | A Diagnostic Application of HEC–HMS Model to Evaluate the Potential for Water Management in the Ratnapura Watershed of Kalu Ganga Sri Lanka | |
Georgakakos | Probabilistic climate-model diagnostics for hydrologic and water resources impact studies | |
Coville et al. | Loss of street trees predicted to cause 6000 L/tree increase in leaf-on stormwater runoff for Great Lakes urban sewershed | |
Li et al. | Ship traffic flow prediction based on AIS data mining | |
Yasin | Mixed geographically weighted regression model (Case study: The percentage of poor households in Mojokerto 2008) | |
Mendicino et al. | Integrated drought watch system: a case study in Southern Italy |
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 |