CN111639810B - 一种基于防汛需求的降雨预报精度评估方法 - Google Patents

一种基于防汛需求的降雨预报精度评估方法 Download PDF

Info

Publication number
CN111639810B
CN111639810B CN202010483351.8A CN202010483351A CN111639810B CN 111639810 B CN111639810 B CN 111639810B CN 202010483351 A CN202010483351 A CN 202010483351A CN 111639810 B CN111639810 B CN 111639810B
Authority
CN
China
Prior art keywords
rainfall
forecast
evaluation
index
value
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
Application number
CN202010483351.8A
Other languages
English (en)
Other versions
CN111639810A (zh
Inventor
张卫国
顾巍巍
江雨田
赵思远
张焱
朱从飞
范仲丽
钟伟
肖涵
王新龙
林宇
张玉兰
王晓晓
孙飞飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ningbo Water Conservancy And Hydropower Planning And Design Institute Co ltd
Original Assignee
Ningbo Water Conservancy And Hydropower Planning And Design Institute Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ningbo Water Conservancy And Hydropower Planning And Design Institute Co ltd filed Critical Ningbo Water Conservancy And Hydropower Planning And Design Institute Co ltd
Priority to CN202010483351.8A priority Critical patent/CN111639810B/zh
Publication of CN111639810A publication Critical patent/CN111639810A/zh
Application granted granted Critical
Publication of CN111639810B publication Critical patent/CN111639810B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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等降雨预报指标，不能适用于洪水预报等防汛应用场景需求。
为此，本发明申请提供了一种基于防汛需求的降雨预报精度评估方法，本发明综合了降雨时程分配、空间分布、总量量级三个维度的降雨预报精度评价，其评估的结果可以为防汛场景下的降雨预报模式的防汛应用选取，以及多种降雨数值预报模式的融合应用提供依据。
发明内容
针对现有技术存在的问题为此，本发明提出了一种基于防汛需求的降雨预报精度评估方法，该方法是一种综合量级、空间分布、时程分配三个维度的降雨预报精度评价方法：在量级上，建立以相对误差为基础的量级评价指标；在空间分布上，引入容差区间的概念，进行逐网格赋分，综合全流域网格后得出空间分布评价指标；在时程上，采用确定性系数评价时间过程的匹配程度，得到时程分配评价指标。其评价的结果可以为降雨预报模式的应用选取，以及多种预报模式的融合应用提供参考。所述方法包括如下步骤：
第一步：网格化与面雨量计算
将数值预报降雨值与实测降雨值统一到统一尺寸的标准网格下；按照统一设定的标准网格尺寸，逐时段采用克里金插值得到网格降雨值，将流域内各网格预报值取算术平均即得到流域面雨量预报值；实测数据主要来源于雨量测站，同样逐时段采用克里金插值得到网格降雨值，将流域内各网格降雨值取算术平均即得到流域面雨量实测值；由此，得到各网格以及全流域的逐时段实测、预报降雨量，为后续步骤提供数据基础；
第二步：降雨预报量级评价指标
对选定的流域总雨量进行量级精度评价，考虑到与时程分配、空间分布指标协调性，选取一个无量纲的值，引入相对误差值项，并按照指标的正向性考虑，取累积总降雨量的量级评价指标为：
Figure BDA0002518173650000031
当Esum,all＜0时，令
Esum,all＝0 (2)
式中：Psum,all为实测全时段累积总降雨量；
Figure BDA0002518173650000032
为预报全时段累积总降雨量；
Esum,all作为量级评价指标，值在[0,1]区间，值越大，评价指标越优；在实测降雨量很小时，会引起的表述失真；根据避免误差危害的若干原则，结合精度评价结果的防洪应用场景，当实测降雨Psum,all≤Tmm时，式(1)中除数取常数T，即：
Figure BDA0002518173650000033
当Esum,all＜0时，令
Esum,all＝0 (4)
除了上述全时段累积总降雨之外，在防汛实践中，还关注考虑最大1h、3h、6h、12h、24h等特征时段累积降雨量，为此，相应的针对最大i小时的累积雨量量级评估公式如下：
Figure BDA0002518173650000034
当Esum,ih＜0时，令
Esum,ih＝0 (6)
式中：Psum,ih为实测最大i小时累积降雨量；
Figure BDA0002518173650000035
为实测最大i小时降雨对应时段的预报累积降雨量；i＝1,3,6,12,24；
最后，将全时段、各个特征时段雨量量级评价指标综合，最终的量级评价指标为
Figure BDA0002518173650000036
式中：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个网格的全时段、各个特征时段的空间相符性指标加权，得到该网格的空间相符性指标为
Figure BDA0002518173650000041
再将流域内所有网格的空间相符性指标取算术平均，得到流域的空间分布评价指标Edis，数值在[0,1]区间，值越大，指标越优：
Figure BDA0002518173650000042
式中，n为流域内的总网格数；
第四步：降雨预报时程分配评价指标
定义洪水预报的常用指标“确定性系数”作为时程评价指标Epro，数值在[0,1]区间，值越大，指标越优：
Figure BDA0002518173650000051
式中，Pt为第t时段流域实测面雨量，
Figure BDA0002518173650000052
为第t时段流域预报面雨量，
Figure BDA0002518173650000053
为流域实测面雨量系列的均值；
第五步：降雨预报综合评价指标
将第二至第四步中得到的量级评价指标、空间分布评价指标和时程分配评价指标进行加权计算，得到降雨的综合评价指标E，其数值在[0,1]区间，值越大，指标越优；
E＝Esum*Csum+Edis*Cdis+Epro*Cpro (11)
式中，Csum为量级权重系数，Cdis为空间分布权重系数，Cpro为时程分配权重系数；这三个权重系数根据对洪水洪量、洪峰等关注侧重不同而定。
优选的，在所述第二步中，当实测降雨Psum,all≤10mm时，式(1)中除数取常数T取为10，即
Figure BDA0002518173650000054
关注考虑最大1h、3h、6h、12h、24h等特征时段累积降雨量，式(1)中除数取常数T取为10，相应的针对最大i小时的累积雨量量级评估公式如下：
Figure BDA0002518173650000055
优选的，在所述第三步中，
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所示，本发明将数值预报降雨与实测降雨值统一到标准网格下，分别进行量级、空间分布、时程分配的吻合程度评估，得到三个维度的分项指标，综合形成对预报降雨的综合评价指标，具体步骤如下：
第一步：网格化与面雨量计算
由于涉及空间分布维度的评价，需要将数值预报降雨值与实测降雨值统一到统一尺寸的标准网格下。各数值预报模式水平分辨率不一，本发明按照统一设定的标准网格尺寸，逐时段采用克里金插值得到网格降雨值，将流域内各网格预报值取算术平均即得到流域面雨量预报值。实测数据主要来源于雨量测站，同样逐时段采用克里金插值得到网格降雨值，将流域内各网格降雨值取算术平均即得到流域面雨量实测值。
由此，可以得到各网格以及全流域的逐时段实测、预报降雨量，为本发明以下的评价步骤提供数据基础。
第二步：降雨预报量级评价指标
对选定的流域总雨量进行量级精度评价，考虑到与下文的时程分配、空间分布指标协调性，需选取一个无量纲的值，本文引入相对误差值项，并按照指标的正向性考虑，拟取累积总降雨量的量级评价指标为
Figure BDA0002518173650000071
当Esum,all＜0时，令
Esum,all＝0 (2)
式中：Psum,all为实测全时段累积总降雨量；
Figure BDA0002518173650000072
为预报全时段累积总降雨量。
Esum,all作为量级评价指标，值在[0,1]区间，值越大，评价指标越优。但上式在实测降雨量很小时，会引起的表述失真，例如实测降雨0.1mm，预报降雨0.2mm，按上式计算值Esumall＝0，直接的表述理解为量级预报指标到恶劣下限，而从本文立意的防洪应用场景来看降雨0.2mm与0.1mm引起的水文响应差别很小，此时Esum,all为接近1的数才为合理。根据避免误差危害的若干原则，结合精度评价结果的防洪应用场景，本文规定当实测降雨Psum,all≤10mm时，上式除数取常数10，即
Figure BDA0002518173650000073
当Esum,all＜0时，令
Esum,all＝0 (4)
除了上述全时段累积总降雨之外，在防汛实践中，往往还关注考虑最大1h、3h、6h、12h、24h等特征时段累积降雨量。为此，相应的针对最大i小时的累积雨量量级评估公式如下：
Figure BDA0002518173650000074
当Esum,ih＜0时，令
Esum,ih＝0 (6)
式中：Psum,ih为实测最大i小时累积降雨量；
Figure BDA0002518173650000075
为实测最大i小时降雨对应时段的预报累积降雨量；i＝1,3,6,12,24。
最后，将全时段、各个特征雨量量级评价指标综合，最终的量级评价指标为
Figure BDA0002518173650000081
式中：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个网格的全时段、各个特征时段的空间相符性指标加权，得到该网格的空间相符性指标为
Figure BDA0002518173650000082
再将流域内所有网格的空间相符性指标取算术平均，得到流域的空间分布评价指标Edis，数值在[0,1]区间，值越大，指标越优：
Figure BDA0002518173650000091
式中，n为流域内的总网格数。
第四步：降雨预报时程分配评价指标
为了评价降雨预报过程与实测过程之间的吻合程度，引入洪水预报的常用指标“确定性系数”，作为时程评价指标Epro，数值在[0,1]区间，值越大，指标越优：
Figure BDA0002518173650000092
式中，Pt为第t时段流域实测面雨量，
Figure BDA0002518173650000093
为第t时段流域预报面雨量，Pt为流域实测面雨量系列的均值。
第五步：降雨预报综合评价指标
将以上步骤中得到的量级评价指标、空间分布评价指标和时程分配评价指标进行加权计算，得到降雨的综合评价指标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)
式中：Psum,all为实测全时段累积总降雨量；为预报全时段累积总降雨量；
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)
式中：Psum,ih为实测最大i小时累积降雨量；为实测最大i小时降雨对应时段的预报累积降雨量；i＝1,3,6,12,24；
最后，将全时段、各个特征时段雨量量级评价指标综合，最终的量级评价指标为
式中：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]区间，值越大，指标越优：
式中，Pt为第t时段流域实测面雨量，为第t时段流域预报面雨量，为流域实测面雨量系列的均值；
第五步：降雨预报综合评价指标
将第二至第四步中得到的量级评价指标、空间分布评价指标和时程分配评价指标进行加权计算，得到降雨的综合评价指标E，其数值在[0,1]区间，值越大，指标越优；
E＝Esum*Csum+Edis*Cdis+Epro*Cpro; (11)
式中，Csum为量级权重系数，Cdis为空间分布权重系数，Cpro为时程分配权重系数；这三个权重系数根据对洪水洪量、洪峰关注侧重不同而定。
2.根据权利要求1所述的一种基于防汛需求的降雨预报精度评估方法，其特征在于：在所述第二步中，
当实测降雨Psum,all≤10mm时，式(1)中除数取常数T取为10，即
关注考虑最大1h、3h、6h、12h、24h特征时段累积降雨量，式(1)中除数取常数T取为10，相应的针对最大i小时的累积雨量量级评估公式如下：
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。
CN202010483351.8A 2020-06-01 2020-06-01 一种基于防汛需求的降雨预报精度评估方法 Active CN111639810B (zh)

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 CN111639810A (zh) 2020-09-08
CN111639810B true 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)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112614309B (zh) * 2020-12-18 2021-10-01 中国地质环境监测院 一种地质灾害气象预警准确率的计算方法
CN112837507A (zh) * 2021-02-08 2021-05-25 国家海洋局东海预报中心 一种海浪预警报质量评估方法及装置
CN113009595B (zh) * 2021-02-09 2023-05-02 海南省气象科学研究所 一种基于车载传感器的降雨监测方法及相关装置
CN114239925A (zh) * 2021-11-29 2022-03-25 中国长江电力股份有限公司 一种水电站入库洪水预报偏差分割和评估方法
CN116068676B (zh) * 2023-03-07 2023-06-02 南京气象科技创新研究院 一种降水预报跨量级通用综合评价方法

Citations (12)

* Cited by examiner, † Cited by third party
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 中国水利水电科学研究院 复杂河道地形流域洪水演进虚拟仿真的淹没线追踪方法
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 河海大学 基于多维降雨特征空间最优决策的山洪快速预警预报方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106529739A (zh) * 2016-11-30 2017-03-22 福建四创软件有限公司 基于网格化的洪水预测方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
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 中国水利水电科学研究院 复杂河道地形流域洪水演进虚拟仿真的淹没线追踪方法
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)

* Cited by examiner, † Cited by third party
Title
王洋 ; 刘佳 ; 于福亮 ; 李传哲 ; 田济扬 ; 邱庆泰 ; .基于数据同化的降雨数值空间分布模拟研究.中国水利水电科学研究院学报.2018,(第03期),全文. *

Also Published As

Publication number Publication date
CN111639810A (zh) 2020-09-08

Similar Documents

Publication Publication Date Title
CN111639810B (zh) 一种基于防汛需求的降雨预报精度评估方法
CN112070286B (zh) 复杂地形流域的降水预报预警系统
Xu et al. Assessment of runoff and sediment yield in the Miyun Reservoir catchment by using SWAT model
CN106991278A (zh) 一种集合降水预报与实时洪水概率预报的耦合方法
CN111507375B (zh) 一种城市内涝风险快速评估方法及系统
KR20170005553A (ko) 지능형 정보제공을 위한 홍수, 가뭄정보의 평가 및 예측 기법 개발
Pruski et al. Improved regionalization of streamflow by use of the streamflow equivalent of precipitation as an explanatory variable
CN112381285B (zh) 一种基于遥感的洪涝淹没预测方法
Li et al. Analysis of the relationship between soil erosion risk and surplus floodwater during flood season
Msadala et al. Sediment yield prediction for South Africa: 2010 edition
CN110334404A (zh) 一种流域尺度骤发干旱精准识别方法
Chang et al. Spatial and temporal variations of precipitation during 1979–2015 in Jinan City, China
Rangsiwanichpong et al. Assessment of sediment yield in Thailand using revised universal soil loss equation and geographic information system techniques
CN111915158A (zh) 一种基于Flood Area模型的暴雨灾害天气风险评估方法、装置及设备
Aksoy et al. SPI-based drought severity-duration-frequency analysis
Gharib et al. Evaluation of ModClark model for simulating rainfall-runoff in Tangrah watershed, Iran.
Yan et al. Study on the dynamic characteristics of shallow groundwater level under the influence of climate change and human activities in Cangzhou, China
Zhao et al. Effect of soil and water conservation measures on the reduction of runoff and sediment load in a loess hilly-gully region
Herma Data processing and model choice for flood prediction
Deby et al. Analysis of Wanggu River Flood Inundation Kendari City Southeast Sulawesi Province Using HEC RAS 5.0. 6
CN115358587A (zh) 区域多部门协同的基础设施规划方法及系统
Li et al. Runoff characteristics of the Nen River Basin and its cause
CHEN et al. Evaluation of land-use efficiency based on regional scale:—a case study in Zhanjiang, guangdong province
Zhang et al. Sustainable effect of coal mining on hydrological drought risk: a case study on Kuye river basin, China
CN110751398A (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