[发明专利]一种自动修正的多模式数值降雨集合预报方法有效
申请号: | 201610131565.2 | 申请日: | 2016-03-08 |
公开(公告)号: | CN105808948B | 公开(公告)日: | 2017-02-15 |
发明(设计)人: | 刘佳;李传哲;田济扬;于福亮;王洋 | 申请(专利权)人: | 中国水利水电科学研究院 |
主分类号: | G06F19/00 | 分类号: | G06F19/00 |
代理公司: | 成都正华专利代理事务所(普通合伙)51229 | 代理人: | 李蕊 |
地址: | 100038 *** | 国省代码: | 北京;11 |
权利要求书: | 暂无信息 | 说明书: | 暂无信息 |
摘要: | 本发明涉及一种自动修正的多模式数值降雨集合预报方法,包括以下步骤步骤1、各种数值大气预报模式的选取;步骤2、模拟预报,输出每隔T小时降雨数据;步骤3、降雨预报结果的评价;步骤4、各模式预报权重系数的确定;步骤5、发布本次预报结果。本发明在现有多模式集合降雨预报基础上,可以更加客观的评价各数值大气模式的降雨预报结果,集合降雨预报的最终结果不必过分依赖于人工决策,公布的降雨预报结果更具客观性。 | ||
搜索关键词: | 一种 自动 修正 模式 数值 降雨 集合 预报 方法 | ||
【主权项】:
一种自动修正的多模式数值降雨集合预报方法,包括以下步骤:步骤1、各种数值大气预报模式的选取;步骤2、模拟预报,输出每隔T小时降雨数据;步骤3、降雨预报结果的评价;步骤3中在各选取的数值大气模式输出6h降雨预报结果后,依据本次降雨的实测结果,对降雨预报结果进行定性与定量、时间与空间、点雨量与面雨量的综合评价,并对评价结果进行打分;步骤3中综合评价包括定性评价,所述定性评价首先将降雨的预报值与实测值对比并进行分级评估,再根据评估结果构建分类评价指标,具体来说:当使用分类指标在空间尺度进行评价时,首先对某一特定观测时刻i,处于不同观测位置上的预测值与实测值进行对比,获得降雨量等级表中的分类变量NAi、NBi和NCi,继而将观测时段内不同时刻对应的分类指标按照1‑4式进行统计平均,最终得到空间尺度上的分类评价结果,其中空间尺度评价指标为:准确率指标频率偏差指标误报率指标临界成功率指标其中,NAi、NBi和NCi分别代表第i个6h观测时段内的不同观测位置上的预报值与观测值是否在降雨量等级表对应的降雨等级中,N为观测时段的个数,而面雨量为各雨量站处雨量的平均值;在进行时间尺度评价时,首先对某一特定观测位置j上,不同观测时刻的预测值与实测值进行对比,统计降雨量等级表中的分类变量,之后将研究区域各个观测位置的分类指标按照5‑8式进行统计平均,最终得到时间尺度上的分类评价结果,其中时间尺度评价指标为:准确率指标频率偏差指标误报率指标临界成功率指标其中,NAj、NBj和NCj分别代表第j个观测位置上不同时刻的预报值与观测值是否在降雨量等级表中对应的降雨等级中,M为观测位置的个数:所述降雨量等级表为:降雨等级小雨中雨大雨暴雨大暴雨特大暴雨6h雨量(mm)0.1~2.52.6~66.1~1212.1~2525.1~60>60上述变量NAi、NBi和NCi的计算方法为:对于空间尺度评价时,某一个观测时段i内,如果观测位置降雨的预测值与观测值均在上述六个降雨等级任何之一范围内,则给NAi记1;如果降雨的观测值在上述六个降雨等级任何之一范围内,而预测值不在该范围内,但不为0,则给NBi记1;如果降雨的观测值在上述六个降雨等级任何之一范围内,而预测值为0mm,即数值大气模式未捕捉到降水,则给NCi记1;上述变量NAj、NBj和NCj的计算方法为:对于时间尺度评价时,某一特定观测位置j内,如果观测位置降雨的预测值与观测值均在上述六个降雨等级任何之一范围内,则给NAj记1;如果降雨的观测值在上述六个降雨等级任何之一范围内,而预测值不在该范围内,但不为0,则给NBj记1;如果降雨的观测值在上述六个降雨等级任何之一范围内,而预测值为0mm,即数值大气模式未捕捉到降水,则给NCj记1;步骤3中综合评价还包括定量评价,定量评价采用在误差分析中常用的4个定量评价指标,当进行时间尺度评价时,Pi和Oi分别为在观测时刻i,研究区域面平均雨量的预报值和实测值;如9‑12式所示:最大误差MEt(maximum error)=max|Pi‑Oi| (9);均方根误差平均偏差标准差其中,i为不同的观测时段,N为观测时段的个数;所述MBE为平均偏差MBEt的数值;当进行空间尺度评价时,Pj和Qj分别为在某个特定的空间位置j,整个观测时段内累积雨量的预报值和实测值,如13‑16式所示:最大误差MEs(maximum error)=max|Pj‑Oj| (13);均方根误差平均偏差标准差其中,j为不同的观测地点,M为观测点的个数;所述MBE为平均偏差MBEs的数值;步骤3中采用上述8个分类评价指标与8个定量评价指标,构建各数值大气模式降雨预报的指标体系,进而依据上述16个评价指标对各个数值大气模式的降雨预报结果进行评分;假定选用m个数值大气模式,则对每一个评价指标进行归一化处理,如对于k个数值大气模式的指标PODtk:SPODtk=(PODtk‑PODtmin)/(PODtmax‑PODtmin) (17);其中,k取1、……、m,k为数值大气模式的个数;PODtmax和PODtmin分别为m个数值大气模式对应的m个PODt的最大值和最小值;其他评价指标的归一化处理参照上述公式(17)进行计算;归一化处理后,在对每一个数值天气模式打分,综合评分用S表示,Sk表示第k个数值天气模式的综合评分:步骤4、各模式预报权重系数的确定;步骤5、发布本次预报结果。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国水利水电科学研究院,未经中国水利水电科学研究院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201610131565.2/,转载请声明来源钻瓜专利网。
- 同类专利
- 专利分类
G06 计算;推算;计数
G06F 电数字数据处理
G06F19-00 专门适用于特定应用的数字计算或数据处理的设备或方法
G06F19-10 .生物信息学,即计算分子生物学中的遗传或蛋白质相关的数据处理方法或系统
G06F19-12 ..用于系统生物学的建模或仿真,例如:概率模型或动态模型,遗传基因管理网络,蛋白质交互作用网络或新陈代谢作用网络
G06F19-14 ..用于发展或进化的,例如:进化的保存区域决定或进化树结构
G06F19-16 ..用于分子结构的,例如:结构排序,结构或功能关系,蛋白质折叠,结构域拓扑,用结构数据的药靶,涉及二维或三维结构的
G06F19-18 ..用于功能性基因组学或蛋白质组学的,例如:基因型–表型关联,不均衡连接,种群遗传学,结合位置鉴定,变异发生,基因型或染色体组的注释,蛋白质相互作用或蛋白质核酸的相互作用
G06F 电数字数据处理
G06F19-00 专门适用于特定应用的数字计算或数据处理的设备或方法
G06F19-10 .生物信息学,即计算分子生物学中的遗传或蛋白质相关的数据处理方法或系统
G06F19-12 ..用于系统生物学的建模或仿真,例如:概率模型或动态模型,遗传基因管理网络,蛋白质交互作用网络或新陈代谢作用网络
G06F19-14 ..用于发展或进化的,例如:进化的保存区域决定或进化树结构
G06F19-16 ..用于分子结构的,例如:结构排序,结构或功能关系,蛋白质折叠,结构域拓扑,用结构数据的药靶,涉及二维或三维结构的
G06F19-18 ..用于功能性基因组学或蛋白质组学的,例如:基因型–表型关联,不均衡连接,种群遗传学,结合位置鉴定,变异发生,基因型或染色体组的注释,蛋白质相互作用或蛋白质核酸的相互作用