本发明专利技术涉及一种基于再分析数据的海浪有效波高的长期趋势预测方法，其特征在于步骤包括：一、收集ERA-Interim各时次气象预报数据；二、获取各格点坐标；三、计算SLP距平值及标准偏差；四、SLP距平值主成分分析；五、对海区数据进行Box-Cox变换；六、计算海浪有效波高的预测因子；七、计算有效波高和预测因子的标准偏差；八、预测因子带入预测模型；九、有效波高滞后值带入模型；十、在EOF基础上的SLP场预测；十一、优化选择预测因子；十二、模型预测海浪有效波高；十三、评估预测水平；十四、计算海浪有效波高长期趋势；十五、绘制有效波高长期趋势图。本发明专利技术可预报多时次的有效波高的长期趋势，且准确率高。
【技术实现步骤摘要】
本专利技术属于海浪参数预报
，特别是涉及一种基于再分析数据的海浪有效波高的长期趋势预测方法。
技术介绍
海浪对人们的生产生活有着不可忽视的影响，如海上航行、沿海港口建设、航道工程、渔业生产等都与海浪有密切关系。另外，海上石油平台的安全也与海浪息息相关。海浪有效波高就是反映海浪特征的一个重要参数，因此分析预测海浪有效波高的趋势具有重要的现实意义。传统的观测手段如浮标等，虽然能够精确的获得海浪波高的变化信息，但它们只能获得海浪在固定点的变化，而且覆盖面也非常有限，目前很难在中国海域得到超过20年的连续的海面波浪的浮标观测数据。随着卫星遥感技术的成熟，卫星数据逐渐被应用，有关海浪波高的卫星数据虽有较广覆盖范围，但最多也只是近20年的资料，这就严重制约了对海浪有效波高长期趋势研究的可靠性。如何克服现有技术的不足已成为当今海浪参数预报
中亟待解决的重点难题之一。
技术实现思路
本专利技术的目的是为克服现有技术的不足而提供一种基于再分析数据的海浪有效波高的长期趋势预测方法，本专利技术利用全球先进稳定的再分析数据源，采用Box-Cox变换对原始数据进行修正，再依据修正后的海平面气压、海浪有效波高等气象数据，采用主成分分析方法和长期波高趋势公式，计算和预测各时次的海浪有效波高的长期趋势，具有很强的可操作性。根据本专利技术提出的一种基于再分析数据的海浪有效波高的长期趋势预测方法，其特征在于包括如下具体步骤：步骤一，收集基于格点模式的欧洲中尺度天气预测中心的ERA-Interim再分析数据集的20～30年时间段的各时次气象预报数据，其中各时次气象预报数据是指包括4～8小时一次的海平面气压SLP和海浪有效波高Hs；步骤二，获取所收集的各时次气象预报数据所标格点的坐标，以该坐标为依据，提取与所述各时次气象预报数据所标格点的坐标相对应的海平面气压矩阵S，如(1)式所示，海浪有效波高矩阵H，如(2)式所示，其中包括m个空间点，每个空间点含有n次观测数据：步骤三，计算基于格点模式的ERA-Interim各时次的海平面气压SLP的均值M，再用原始值S减去均值M，得到基于格点模式的各时次的SLP的距平值P，并计算出SLP距平值P的标准偏差S，如(3)式所示：上述(3)式中：mi=Σj-1nSijn;]]>步骤四，对SLP距平值P做EOF分析，得到不同成分及各成分对总方差的贡献率，保留前30个EOF和主成分；其中：对P进行协方差计算，得到实对称矩阵Lm×m·，其中：Lm×m=1nPPT,]]>T表示矩阵的转置；然后求协方差矩阵Lm×m的特征向量V和特征值Λ，如(4)式所示，以满足LV＝ΛV,其中矩阵V是正交矩阵，矩阵V的第j列元素就是特征值λj对应的特征向量；根据实对称矩阵Lm×m的特征向量V和特征值FL计算每个特征向量的方差贡献率和前几个特征向量的累计方差贡献率；按照特征值从大到小的顺序对L进行排序，排在第一位的为EOF1，以此类推；步骤五，对根据步骤一收集的基于格点的各时次的原始海平面气压SLP和海浪有效波高Hs进行Box-Cox变换，得到变换后的海平面气压trGt和海浪有效波高trHt；步骤六，对每个格点上对应的trHt，用PCk，t和PCk，t-4计算其相关系数，并取相关系数最高时的28个PCk,t或PCk,t-4作为海浪有效波高的预测因子；步骤七，计算海浪有效波高的标准偏差SHl和30个预测因子Xk,t的标准偏差SXk，保存备用；步骤八，将预测因子带入预测模型，用F统计量比较第i个模型和第i+1个模型的预测结果，从而选出最优的预测因子；步骤九，将滞后一步的海浪有效波高也带入模型，作为预测因子之一，综合预测下一时次的各格点的海浪有效波高，优化模型参数，得到最终模型；其中模型如(5)式所示：Ht=a+Σk-1KbkXk,t+Σp-1PcpHt-p+ut---(5),]]>上述(5)式中Ht是每个网格点上的经过变换的海浪有效波高，Ht-p是滞后p的海浪有效波高，P是跟预报量相关的参变量的滞后系数，Xk，t是第k个基于SLP的预报因子，ut可以用M阶自回归模型来表示，如果M＝0，ut就是白噪声；步骤十，在步骤四得到的前30个EOF的基础上对各时次的SLP场进行预测，得到PCk,t；步骤十一，用步骤七保存备用的SXk衡量选择30个预测因子Xk,t；步骤十二，将所有预测因子带入步骤九的最终模型，预测目标时期内各时次的海浪有效波高，将预测出的有效波高值还原到Box-Cox变换前的值，保存为格点模式文件；步骤十三，采用RMSE等评估指标评估预测水平；步骤十四，以步骤十二预测的海浪有效波高为依据，用趋势计算公式计算海浪有效波高的长期趋势，最终得到海浪有效波高的长期趋势结果；步骤十五，根据步骤十二的结果，对应到相应的格点坐标，绘制出海浪有效波高长期趋势图。本专利技术与与现有技术相比其显著优点：一是本专利技术利用全球先进稳定的再分析数据源，将数据建立在具有几十年甚至跨百年的海浪有效波高资料的再分析数据的基础上，从而解决了海浪有效波高长期趋势分析的可靠性；二是本专利技术采用Box-Cox变换对原始数据进行修正，再依据修正后的海平面气压、海浪有效波高等气象数据，采用主成分分析方法和短期波高趋势公式，计算和预测各时次的有效波高的长期趋势，预报海浪有效波高长期趋势的准确率高；三是本专利技术可有效指导各沿海区域的海浪防护工作，对于维护沿海地带的安全稳定、减轻海浪灾害起到至关重要和不可或缺的作用，具有很强的可操作性。附图说明图1为本专利技术提出的一种基于再分析数据的海浪有效波高的长期趋势预测方法的流程方框示意图。图2为应用本专利技术提出的一种基于再分析数据的海浪有效波高的长期趋势预测方法绘制的中国某海域冬季最大海浪有效波高的长期趋势结果示意图。具体实施方式下面结合附图和实施例对本专利技术的具体实施方式作进一步的详细说明。现以中国某海区域为例，应用本专利技术提出的一种基于再分析数据的海浪有效波高的长期趋势预测方法来预报海浪有效波高的长期趋势，结合图1，其具体步骤包括如下：步骤一，收集基于格点模式的中国某海区域的ERA-Interim再分析数据集的1981-2000年期间的各时次海平面气压SLP和有效波高Hs数据，该数据间隔为每6小时一次；但不限于此，可收集基于格点模式的欧洲中尺度天气预测中心的ERA-Interim再分析数据集的20本文档来自技高网...
【技术保护点】
一种基于再分析数据的海浪有效波高的长期趋势预测方法，其特征在于包括以下具体步骤：步骤一，收集基于格点模式的欧洲中尺度天气预测中心的ERA‑Interim再分析数据集的20～30年时间段的各时次气象预报数据，其中各时次气象预报数据是指包括4～8小时一次的海平面气压SLP和海浪有效波高Hs；步骤二，获取所收集的各时次气象预报数据所标格点的坐标，以该坐标为依据，提取与所述各时次气象预报数据所标格点的坐标相对应的海平面气压矩阵S，如(1)式所示，海浪有效波高矩阵H，如(2)式所示，其中包括m个空间点，每个空间点含有n次观测数据：步骤三，计算基于格点模式的ERA‑Interim各时次的海平面气压SLP的均值M，再用原始值S减去均值M，得到基于格点模式的各时次的SLP的距平值P，并计算出SLP距平值P的标准偏差S，如(3)式所示：上述(3)式中：mi=Σj-1nSijn;]]>步骤四，对SLP距平值P做EOF分析，得到不同成分及各成分对总方差的贡献率，保留前30个EOF和主成分；其中：对P进行协方差计算，得到实对称矩阵Lm×m其中：Lm×m=1nPPT,]]>T表示矩阵的转置；然后求协方差矩阵Lm×m的特征向量V和特征值Λ，如(4)式所示，以满足LV＝ΛV,其中矩阵V是正交矩阵，矩阵V的第j列元素就是特征值λj对应的特征向量；根据实对称矩阵Lm×m的特征向量V和特征值‑I,计算每个特征向量的方差贡献率和前几个特征向量的累计方差贡献率；按照特征值从大到小的顺序对L进行排序，排在第一位的为EOF1，以此类推；步骤五，对根据步骤一收集的基于格点的各时次的原始海平面气压SLP和海浪有效波高Hs进行Box‑Cox变换，得到变换后的海平面气压trGt和海浪有效波高trHt；步骤六，对每个格点上对应的trHt，用PCk，t和PCk，t‑4计算其相关系数，并取相关系数最高时的28个PCk,t或PCk,t‑4作为海浪有效波高的预测因子；步骤七，计算海浪有效波高的标准偏差SHl和30个预测因子Xk,t的标准偏差SXk，保存备用；步骤八，将预测因子带入预测模型，用F统计量比较第i个模型和第i+1个模型的预测结果，从而选出最优的预测因子；步骤九，将滞后一步的海浪有效波高也带入模型，作为预测因子之一，综合预测下一时次的各格点的海浪有效波高，优化模型参数，得到最终模型；其中模型如(5)式所示：Ht=a+Σk-1KbkXk,t+Σp-1PcpHt-p+ut---(5),]]>上述(5)式中Ht是每个网格点上的经过变换的海浪有效波高，Ht‑p是滞后p的海浪有效波高，P是跟预报量相关的参变量的滞后系数，Xk，t是第k个基于SLP的预报因子，ut可以用M阶自回归模型来表示，如果M＝0，ut就是白噪声；步骤十，在步骤四得到的前30个EOF的基础上对各时次的SLP场进行预测，得到PCk,t；步骤十一，用步骤七保存备用的SXk衡量选择30个预测因子Xk,t；步骤十二，将所有预测因子带入步骤九的最终模型，预测目标时期内各时次的海浪有效波高，将预测出的有效波高值还原到Box‑Cox变换前的值，保存为格点模式文件；步骤十三，采用RMSE等评估指标评估预测水平；步骤十四，以步骤十二预测的海浪有效波高为依据，用趋势计算公式计算海浪有效波高的长期趋势，最终得到海浪有效波高的长期趋势结果；步骤十五，根据步骤十二的结果，对应到相应的格点坐标，绘制出海浪有效波高长期趋势图。...
【技术特征摘要】
1.一种基于再分析数据的海浪有效波高的长期趋势预测方法，其特征在于
包括以下具体步骤：
步骤一，收集基于格点模式的欧洲中尺度天气预测中心的ERA-Interim再分
析数据集的20～30年时间段的各时次气象预报数据，其中各时次气象预报数据
是指包括4～8小时一次的海平面气压SLP和海浪有效波高Hs；
步骤二，获取所收集的各时次气象预报数据所标格点的坐标，以该坐标为依
据，提取与所述各时次气象预报数据所标格点的坐标相对应的海平面气压矩阵
S，如(1)式所示，海浪有效波高矩阵H，如(2)式所示，其中包括m个空间
点，每个空间点含有n次观测数据：
步骤三，计算基于格点模式的ERA-Interim各时次的海平面气压SLP的均值
M，再用原始值S减去均值M，得到基于格点模式的各时次的SLP的距平值P，并
计算出SLP距平值P的标准偏差S，如(3)式所示：
上述(3)式中：mi=Σj-1nSijn;]]>步骤四，对SLP距平值P做EOF分析，得到不同成分及各成分对总方差的贡
献率，保留前30个EOF和主成分；其中：
对P进行协方差计算，得到实对称矩阵Lm×m其中：
Lm×m=1nPPT,]]>T表示矩阵的转置；
然后求协方差矩阵Lm×m的特征向量V和特征值Λ，如(4)式所示，以满足
LV＝ΛV,其中
矩阵V是正交矩阵，矩阵V的第j列元素就是特征值λj对应的特征向量；
根据实对称矩阵Lm×m的特征向量V和特征值-I,计算每个特征向量的方差贡
献率和前几个特征向量的累计方差贡献率；按照特征值从大到小的顺序对L进行
排序，排在第一位的为EOF1，以此类推；
步骤五，对根据步骤一收集的基于格点的各时次的原始海平面气压SLP和海
浪有效波高Hs进行Box-Cox变换，得到变换后的海平面气压trGt和海浪有效波
高trHt；
步骤六，对每个格点上对应的trHt，用PCk，t和PCk，t-4计算其相关系数，并
取相关系数最高时的28个PCk,t或PCk,t-4作为海浪有效波高的预测因子；
步骤七，计算海浪有...
【专利技术属性】
技术研发人员：吴玲莉，张玮，吴腾，焦楚杰，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏;32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。