[发明专利]一种基于数值集合预报极值修正的中长期径流预报方法有效
申请号: | 201810036068.3 | 申请日: | 2018-01-15 |
公开(公告)号: | CN108171007B | 公开(公告)日: | 2019-02-05 |
发明(设计)人: | 杨明祥;雷晓辉;蒋云钟;王浩;张云辉;张梦婕;刘珂;闻昕;权锦 | 申请(专利权)人: | 中国水利水电科学研究院 |
主分类号: | G06F17/18 | 分类号: | G06F17/18;G01W1/10;G01W1/14 |
代理公司: | 北京市盛峰律师事务所 11337 | 代理人: | 于国富 |
地址: | 100038 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 预报 数值集合 径流 预报区域 预报时段 预报数据 数据处理领域 修正 趋势预测 统计预报 水文 尺度 预测 | ||
1.一种基于数值集合预报极值修正的中长期径流预报方法,其特征在于,所述方法包括:
S1,确定预报区域和预报时段;
S2,采用数值集合预报方法对预报时段的预报区域进行预报,得到平均径流预报数据结论α,记为Q_Model_Forecast;
S3,采用海温遥相关方法对预报时段的预报区域进行预报,得到平均径流预报数据结论β,记为Q_SST_Forecast;
根据公式(1)计算Q_SST_Forecast的距平值:
Ano为距平值,为预报时段的同期多年平均值;
S4,预报区域在预报时段的预报结论
当|Ano|≥50%,则基于海温遥相关方法得到的预报时段的径流出现极端情况的可能性高,用Q_SST_Forecast代替Q_Model_Forecast作为最终预报结论;
当|Ano|<50%,基于海温遥相关方法得到的预报时段的径流出现极端情况的可能性低,用Q_Model_Forecast作为最终预报结论;
S2,具体为:
S21,从全球气候预报系统中获取最近5天发布的降水速率预报数据;
S22,对所述降水速率预报数据解析获取全球格网式的降水速率数据,对相同时刻的格网式降水速率数据取算数平均,即得每个栅格的逐6h平均降水速率;
S23,计算获取每个网格的逐6h累积降水量;
S24,根据公式(2)获取预报区域逐24h面平均雨量P;
Pi是预报区域内第i个网格点的6h累积降水量,n为预报区域网格点的数量;
依据预报区域逐24h面平均雨量,估计预报区域在预报时段中所指定的各个自然月的月累积降水量预报数据;
S25,设预报流域在预报时段中指定的自然月t的月累积降水量预报数据Pt,获取历史上该自然月t与Pt最接近的多个实测降水量即各个实测降水量对应的径流量,将各个径流量的算术平均值作为预报流域在自然月t的径流预报意见,重复S25,直至得到所有指定自然月的径流预报意见,将得到的所有径流预报意见组成平均径流预报数据结论α,记为Q_Model_Forecast;
S3,具体为:
S31,获取全球格网型海温数据;
S32,获取自1981年以来,预报区域的待预报断面在预报时段中指定月x的各年历史平均径流数据,记为序列一;
获取全球格网型海温数据中不同网格提前m月份的各年的温度序列,记为序列二,m≤12;
将序列一与序列二进行相关性分析,得到与序列一相关性系数最大的海面位置与提前月份;将与序列一相关性系数最高的序列二记为序列a;
S33,设全球格网型海温数据中第i_r行、第j_r列提前leadMonth的海温数据序列为序列a,绘制序列a与序列一的对照图,序列a的曲线处于时间最末端处的海温数据A对应的径流数据为待预报断面在预报时段中指定月x的待预报径流;
S34,从海温数据历史数据中选取与海温数据A最接近的多个海温数据点,并获取每个海温数据点所对应的历史流量数据,然后将得到的多个历史流量数据的算数平均值,即平均径流预报数据作为待预报断面在预报时段中指定月x的待预报径流,将得到的多个历史流量数据的最大值和最小值作为平均径流预报数据的区间,将平均径流预报数据及其区间作为平均径流预报数据结论β,记为Q_SST_Forecast,完成采用海温遥相关方法对预报时段的预报区域的预报。
2.根据权利要求1所述基于数值集合预报极值修正的中长期径流预报方法,其特征在于,S1,具体为:
接收预报需求,并提取所述预报需求中的待预报断面和预报时段,将待预报断面的上流域作为预报流域。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国水利水电科学研究院,未经中国水利水电科学研究院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810036068.3/1.html,转载请声明来源钻瓜专利网。