CN110222872A - 基于经验正交函数分解的海洋多要素中长期统计预测方法 - Google Patents
基于经验正交函数分解的海洋多要素中长期统计预测方法 Download PDFInfo
- Publication number
- CN110222872A CN110222872A CN201910391323.0A CN201910391323A CN110222872A CN 110222872 A CN110222872 A CN 110222872A CN 201910391323 A CN201910391323 A CN 201910391323A CN 110222872 A CN110222872 A CN 110222872A
- Authority
- CN
- China
- Prior art keywords
- marine
- ocean
- day
- formula
- sample matrix
- 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
- 230000007774 longterm Effects 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims abstract description 96
- 238000004458 analytical method Methods 0.000 claims abstract description 40
- 238000012950 reanalysis Methods 0.000 claims description 40
- 230000001932 seasonal effect Effects 0.000 claims description 14
- 238000010276 construction Methods 0.000 claims description 5
- 238000000926 separation method Methods 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 238000013277 forecasting method Methods 0.000 abstract description 3
- 230000007613 environmental effect Effects 0.000 abstract 1
- 238000010438 heat treatment Methods 0.000 abstract 1
- 230000032683 aging Effects 0.000 description 10
- 238000007619 statistical method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
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"
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于经验正交函数分解的海洋多要素中长期统计预测方法:(1)基于待分析和预测海洋的再分析资料,构建历年逐日海洋统计再分析样本矩阵;(2)将海洋各类要素的季节信号与非季节信号进行分离,得到历年逐日的海洋再分析距平样本矩阵;(3)构建海洋多要素时空四维正交模态;(4)构建海洋多要素中长期统计预测模型;(5)采用步骤4构建的海洋多要素中长期统计预测模型对海洋各类要素进行预报。本发明基于足够长时间序列的海洋再分析产品,通过统计规律,构建海洋动力热力环境中长期预测方法,使海洋预测在传统海洋预报模式的基础之上大幅度提高海洋预报的时效,为有效提高中长期海洋环境安全保障水平奠定技术基础。
Description
技术领域
本发明涉及海洋要素预报技术,特别涉及一种基于经验正交函数(Empiricalorthogonal function,EOF)分解的海洋多要素中长期统计预测方法。本发明方法主要应用于水下移动平台执行任务时的海洋要素预报工作,对海面高度、海表温度等要素做出时效为1~3个月的预报。
背景技术
目前海洋环境保障所需要的海水盐度、温度以及海流、密度等海洋环境信息,主要依赖于大量的历史资料,基于统计分析或再分析等手段,获得海洋气候态统计分析产品或长时间序列海洋再分析产品,但是这类产品是对过去海洋环境状态的描述,而非现在或未来状态的表示和预测。
水下移动平台在执行任务时,通常要在水下潜行3个月左右时间,期间的海洋环境保障主要依靠其自身携带的观测设备或海洋气候态统计分析产品,出于隐蔽性方面的考虑,需要提前预知水下移动平台周围一定范围之内的海洋环境信息,但移动平台自身所携带的观测设备只能观测到航迹上当前的海洋环境信息,而气候态产品也无法提供未来时刻的预报,故需要业务化海洋学为其提供时效为1~3个月左右海洋数值预测预报结果,作为海洋环境安全保障的支撑。
当前海洋预测技术主要依靠数值模式进行,但是模式受对应的气象驱动场最长预报时效只有10天的限制,无法提供10天以上的数值预报产品。同时,传统海洋数值预报模式在参数选择、计算平台方面有很大的要求。
综上所述,战略级海洋环境保障急需研制一种小型化、简便的,且预报时效可达1~3个月的中长期海洋环境预测产品,故本专利提出一种基于经验正交函数(Empiricalorthogonal function,EOF)分解的海洋多要素中长期统计预测方法。
发明内容
本发明针对水下移动平台在执行海洋任务时对于1~3个月左右的海洋环境数值预报结果的需求,提供一种小型化、简便的、中长期海洋环境预测方法。能够弥补传统的数值预报方法由于气象驱动时效的限制而导致的海洋要素预报时效较短的缺陷,且不需要具备高性能计算平台,对计算平台的要求较低。利用本发明中长期统计预报方法,水下移动平台在执行任务时,可以在保证自身隐蔽性的前提下对海洋各要素做出时效为1~3个月左右的环境预报,为有效提高中长期海洋环境战略级安全保障水平奠定技术基础,因此具有较大的科学意义和应用价值。
本发明所采用的技术方案是:一种基于经验正交函数分解的海洋多要素中长期统计预测方法,包括以下步骤:
步骤1,基于待分析和预测海洋的再分析资料,构建历年逐日海洋统计再分析样本矩阵;
步骤2,基于步骤1得到的历年逐日海洋统计再分析样本矩阵,将海洋各类要素的季节信号与非季节信号进行分离,得到历年逐日的海洋再分析距平样本矩阵;
步骤3,基于步骤2信号分离之后的历年逐日的海洋再分析距平样本矩阵,构建海洋多要素时空四维正交模态;
步骤4,基于步骤3得到的海洋多要素时空四维正交模态,构建海洋多要素中长期统计预测模型;
步骤5,采用步骤4构建的海洋多要素中长期统计预测模型对海洋各类要素进行预报。
进一步地,步骤1中,所述的年逐日海洋统计再分析样本矩阵构建方法为:
设海洋有K类要素,指标k代表第k类要素,取值范围(k=1,…,K);各类要素的空间点分别为(N1,…,Nk,…,NK),使用nk代表第k类要素的空间点指标,即(n1=1,…,N1),(nk=1,…,Nk),(nK=1,…,NK);要分析的时段共包含L个时刻,指标l代表第l个时刻tl,使用M年的海洋再分析产品数据集,构造M个样本,指标m代表第m个样本;
第m个样本第k类要素第l个时刻tl的Nk个空间点组成向量,如公式(1)所示:
式中,表示tl时刻、第k类要素的第m个样本,其中,该样本由Nk个空间点组成;表示其中第nk个空间点,nk=1,2,...,Nk;
第m个样本第k类要素整时段L个时刻的时空点组成向量,如公式(2)所示:
式中,表示整时段L个时刻内、第m个样本、第k类要素形成的向量;由公式(1)得到,tl=1,2,...,tL;
从而形成如公式(3)所示的样本矩阵:
式中,X即为历年逐日海洋统计再分析样本矩阵,由公式(2)得到,k=1,2,...,K,m=1,2,...,M;
历年逐日海洋统计再分析样本矩阵X是N×M的矩阵,其中总行数
进一步地,步骤2中,所述的将海洋各类要素的季节信号与非季节信号进行分离为:将步骤1得到的历年逐日海洋统计再分析样本矩阵扣除掉累年逐日海洋再分析统计分析样本矩阵;
其中,基于步骤1得到的历年逐日海洋统计再分析样本矩阵,得到累年逐日海洋再分析统计分析样本矩阵,如公式(4)所示:
式中,为累年逐日海洋再分析统计分析样本矩阵;表示第k类要素的累年逐日再分析值;
将步骤1得到的历年逐日海洋统计再分析样本矩阵扣除掉累年逐日海洋再分析统计分析样本矩阵,得到历年逐日的海洋再分析距平样本矩阵,如公式(5)所示:
式中,为历年逐日的海洋再分析距平样本矩阵;表示第1个样本、第k类要素的再分析距平值,表示整时段L个时刻内、第1个样本、第k类要素形成的向量,由公式(2)得到;σ表示标准化矩阵, 并且,表示标准化因子,表示其中第nk个空间点,表示tl时刻,k类要素中第nk个空间点的平均值。
进一步地,步骤3中,所述的海洋多要素时空四维正交模态构建方法如下:
对历年逐日的海洋再分析距平样本矩阵进行正交模态向量分解,即其中表示正交模态,A表示模态系数,并且满足公式(6):
式中,I表示全矩阵;
构建相关系数矩阵,如公式(7)所示:
式中,表示历年逐日的海洋再分析距平样本矩阵;Ckk,k=1,2,...,K代表第k类要素和第k类要素的相关系数矩阵,对角线上的元素代表同种要素的时空相关系数矩阵,非对角线元素代表不同要素的时空相关系数矩阵;
根据历年逐日的海洋再分析距平样本矩阵的展开式,相关系数矩阵变为:
式中,表示特征值矩阵;
使用雅可比分解,计算相关系数矩阵的特征向量和特征值,同时等价地得到的特征向量和特征值,如公式(9)所示:
其中,λN表示特征向量;
历年逐日海洋统计再分析样本矩阵X采用特征向量表示为:
公式(10)即为海洋多要素时空四维正交模态。
进一步地,步骤4中,所述的海洋多要素中长期统计预测模型构建为:以步骤3得到的海洋多要素时空四维正交模态为基函数,将待分析和预测海洋的当前和未来状态采用这组基函数进行展开,并通过拟合当前海洋各类要素,计算各个模态的系数,从而实现对海洋多要素中长期统计预测模型的构建。
其中,步骤4中,所述的海洋多要素中长期统计预测模型的构建方法具体包括:
设待分析和预测的海洋要素真实场向量用Xt表示,Xt对应的当前时刻之前的一段时间的海洋要素真实场为FXt,其中,F表示投影算符;当前时刻之前一段时间的海洋要素分析称为初始场,用Xb表示,初始场Xb相对于真实场Xt的误差Eb如公式(11)所示:
Eb=Xb-FXt (11)
首先采用历年逐日的海洋再分析距平样本矩阵正交模态向量分解所得的正交模态对初始场Xb进行拟合,构造目标泛函,如公式(12)所示:
式中,Jb(Γ)为目标泛函;Γ表示待优化的正交模态系数;
对目标泛函求极小值,即:
解得海洋多要素中长期统计预测模型构建为:
式中,Xa表示海洋多要素中长期统计预测模型的结果。
本发明的有益效果是:一种基于经验正交函数分解海洋多要素中长期统计预测方法,本发明针对海洋环境中的水下移动平台在执行水下任务时,提供一种全新的小型化、简便的、中长期统计预测方法和模型,与传统海洋数值模式预报相比,本方法基于足够长时间序列的海洋再分析产品,通过统计规律,构建海洋动力热力环境1~3个月的中长期预测方法,脱离了传统海洋数值预报受气象驱动时效10天左右的限制,使海洋预测在传统海洋预报模式的基础之上大幅度提高海洋预报的时效,为有效提高中长期海洋环境安全保障水平奠定技术基础,因此具有较大的科学意义和应用价值。
附图说明
图1:本发明方法流程图;
图2:本发明的海洋多要素模态分解结果图;
图3:本发明海面高度预报时效为90天的对比图;
图3a为海面高度预报时效为90天的再分析结果;
图3b为海面高度预报时效为90天的预报结果;
图4:本发明海表温度预报时效为90天的对比图;
图4a为海表温度预报时效为90天的再分析结果;
图4b为海表温度预报时效为90天的预报结果。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹例举以下实施例,并配合附图详细说明如下:
如附图1所示,一种基于经验正交函数分解的海洋多要素中长期统计预测方法,包括以下步骤:
步骤1,基于待分析和预测海洋的再分析资料,构建历年逐日海洋统计再分析样本矩阵;
步骤2,基于步骤1得到的历年逐日海洋统计再分析样本矩阵,将海洋各类要素的季节信号与非季节信号进行分离,即,将海洋最强的周期信号(季节信号)与其他信号进行分离,得到历年逐日的海洋再分析距平样本矩阵;
步骤3,基于步骤2信号分离之后的历年逐日的海洋再分析距平样本矩阵,构建海洋多要素(海面高度、海表温度等各类要素)时空四维正交模态;
步骤4,以步骤3得到的海洋多要素(海面高度、海表温度等各类要素)时空四维正交模态为基函数,将待分析和预测海洋当前和未来的状态采用这组基函数进行展开,并通过拟合当前海洋各类要素(海面高度、海表温度等各类要素),计算各个模态的系数,从而实现对海洋多要素(海面高度、海表温度等各类要素)中长期统计预测模型的构建;
同时,在海洋多要素(海面高度、海表温度等各类要素)中长期统计预测模型建立完成后,进行基于待分析和预测海洋的再分析资料的后报试验,即,对海洋多要素(海面高度、海表温度等各类要素)中长期统计预测模型进行测试,采用步骤1中的待分析和预测海洋的再分析资料进行预报,测试模型预报的结果是否准确;
步骤5,采用步骤4构建的海洋多要素中长期统计预测模型对海洋各类要素进行1~3个月的中长期预报,提高海洋预报的时效。
下面结合具体实施例和公式对本发明方法作进一步阐述。
(1)基于南海海域的长时间序列海洋再分析产品,时间序列为1958年1月1日至2016年12月31日,构建历年逐日海洋统计再分析样本矩阵。假设海洋有K类要素,指标k代表第k类要素,取值范围(k=1,…,K);各类要素的空间点分别为(N1,…,Nk,…,NK),使用nk代表第k类要素的空间点指标,即(n1=1,…,N1),(nk=1,…,Nk),(nK=1,…,NK);要分析的时段共包含L个时刻,指标l代表第l个时刻tl,使用M年的海洋再分析产品数据集,我们可以构造M个样本,指标m代表第m个样本。
第m个样本第k类要素第l个时刻tl的Nk个空间点组成向量:
式中,表示tl时刻、第k类要素的第m个样本,其中,该样本由Nk个空间点组成;表示其中第nk个空间点,nk=1,2,...,Nk;
第m个样本第k类要素整时段L个时刻的时空点组成向量:
式中,表示整时段L个时刻内、第m个样本、第k类要素形成的向量;由公式(1)得到,tl=1,2,...,tL;
从而形成如公式(3)所示的样本矩阵:
式中,X即为历年逐日海洋统计再分析样本矩阵,由公式(2)得到,k=1,2,...,K,m=1,2,...,M。
上述历年逐日海洋统计再分析样本矩阵X是N×M的矩阵,其中总行数
(2)在海洋信号中,季节信号占了总信号的很大一部分,所以需要将海洋季节信号与非季节信号进行分离,即将历年逐日海洋统计再分析样本矩阵扣除掉累年逐日海洋再分析统计分析样本矩阵,获得历年逐日的海洋再分析距平样本矩阵。
基于上述历年逐日海洋统计再分析样本矩阵X,可以获得累年逐日海洋再分析统计分析样本矩阵
式中,表示第k类要素的累年逐日再分析值;
故历年逐日的海洋再分析距平样本矩阵为:
式中,表示第1个样本、第k类要素的再分析距平值,表示整时段L个时刻内、第1个样本、第k类要素形成的向量,由公式(2)得到;σ表示标准化矩阵, 并且,表示标准化因子,表示第nk个空间点,表示tl时刻,k类要素中第nk个空间点的平均值。
至此,标准化的历年逐日的海洋再分析距平样本矩阵已构造完成。
(3)对标准化历年逐日的海洋再分析距平样本矩阵进行正交模态向量分解,即其中表示正交模态,A表示模态系数,并且满足:
式中,I表示全矩阵;
为了达到正交模态分解的目的,我们构建相关系数矩阵:
上式中每一个矩阵元素Ckk(k=1,2,...,K)都是一个分块矩阵,代表两个下标花样要素(即,第k类要素和第k类要素)的相关系数矩阵,对角线上的元素代表同种要素的时空相关系数矩阵,非对角线元素代表不同要素的时空相关系数矩阵。因此,上述相关系数矩阵不仅有效考虑了同一个海洋要素在不同时空点的相关性,还有效考虑了不同海洋要素之间在不同时空点的相关性。
根据历年逐日的海洋再分析距平样本矩阵的展开式,相关系数矩阵变为:
式中,表示特征值矩阵;
使用雅可比分解,可以计算相关系数矩阵的特征向量和特征值,同时可以等价地得到的特征向量和特征值:
其中,λN表示特征向量;
于是,原始的历年逐日海洋统计再分析样本矩阵X可采用特征向量表示为:
公式(10)即为海洋多要素时空四维正交模态。
至此,我们可以将每一年的样本看作每一个模态,并得到模态系数,如图2所示,前10模态占比达到59%,说明大部分的能量主要集中在这前10个模态中,即我们可以依据这些主要能量来代表所有模态,从而构建相应的中长期预测模型。
(4)利用上述(3)中得到的有效正交模态构造特征向量矩阵,构造中长期统计预测模型。在实际问题中,真正有效的模态只有前r个,即在样本中,大部分的能量主要集中于这前r个模态中,于是我们利用这前r个正交模态构造特征向量矩阵。
假设待分析和预测的海洋要素真实场向量用Xt表示,Xt对应的当前时刻之前的一段时间的海洋要素真实场为FXt,其中,F表示投影算符;当前时刻之前一段时间的海洋要素分析称为初始场,用Xb表示,初始场Xb相对于真实场Xt的误差Eb为:
Eb=Xb-FXt (11)
首先采用历年逐日的海洋再分析距平样本矩阵正交模态向量分解所得的正交模态对初始场Xb进行拟合,则可以构造如下目标泛函:
式中,Jb(Γ)为目标泛函;Γ表示待优化的正交模态系数;
对上述目标泛函求极小值,即:
解得据此,海洋多要素中长期统计预测模型构建为:
式中,Xa表示海洋多要素中长期统计预测模型的结果。
图3为海面高度预报时效为90天的对比图,说明依据经验正交函数分解结果,找到了主要能量。通过对比再分析结果与预报结果可以看出,预报结果把握了海面高的变化趋势,可以实现时效为90天的海面高中长期预报。
图4为海表温度预报时效为90天的对比图,通过对比再分析结果与预报结果可以看出,虽然预报结果与实际再分析结果相比有差别,但把握了总体变化趋势,故其可以实现时效为90天的海表温度中长期预报。
尽管上面结合附图对本发明的优选实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,并不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可以做出很多形式,这些均属于本发明的保护范围之内。
Claims (6)
1.一种基于经验正交函数分解的海洋多要素中长期统计预测方法,其特征在于,包括以下步骤:
步骤1,基于待分析和预测海洋的再分析资料,构建历年逐日海洋统计再分析样本矩阵;
步骤2,基于步骤1得到的历年逐日海洋统计再分析样本矩阵,将海洋各类要素的季节信号与非季节信号进行分离,得到历年逐日的海洋再分析距平样本矩阵;
步骤3,基于步骤2信号分离之后的历年逐日的海洋再分析距平样本矩阵,构建海洋多要素时空四维正交模态;
步骤4,基于步骤3得到的海洋多要素时空四维正交模态,构建海洋多要素中长期统计预测模型;
步骤5,采用步骤4构建的海洋多要素中长期统计预测模型对海洋各类要素进行预报。
2.根据权利要求1所述的一种基于经验正交函数分解的海洋多要素中长期统计预测方法,其特征在于,步骤1中,所述的年逐日海洋统计再分析样本矩阵构建方法为:
设海洋有K类要素,指标k代表第k类要素,取值范围(k=1,…,K);各类要素的空间点分别为(N1,…,Nk,…,NK),使用nk代表第k类要素的空间点指标,即(n1=1,…,N1),(nk=1,…,Nk),(nK=1,…,Nk);要分析的时段共包含L个时刻,指标l代表第l个时刻tl,使用M年的海洋再分析产品数据集,构造M个样本,指标m代表第m个样本;
第m个样本第k类要素第l个时刻tl的Nk个空间点组成向量,如公式(1)所示:
式中,表示tl时刻、第k类要素的第m个样本,其中,该样本由Nk个空间点组成;表示其中第nk个空间点,nk=1,2,…,Nk;
第m个样本第k类要素整时段L个时刻的时空点组成向量,如公式(2)所示:
式中,表示整时段L个时刻内、第m个样本、第k类要素形成的向量;由公式(1)得到,tl=1,2,…,tL;
从而形成如公式(3)所示的样本矩阵:
式中,X即为历年逐日海洋统计再分析样本矩阵,由公式(2)得到,k=1,2,…,K,m=1,2,…,M;
历年逐日海洋统计再分析样本矩阵X是N×M的矩阵,其中总行数
3.根据权利要求1所述的一种基于经验正交函数分解的海洋多要素中长期统计预测方法,其特征在于,步骤2中,所述的将海洋各类要素的季节信号与非季节信号进行分离为:将步骤1得到的历年逐日海洋统计再分析样本矩阵扣除掉累年逐日海洋再分析统计分析样本矩阵;
其中,基于步骤1得到的历年逐日海洋统计再分析样本矩阵,得到累年逐日海洋再分析统计分析样本矩阵,如公式(4)所示:
式中,为累年逐日海洋再分析统计分析样本矩阵;表示第k类要素的累年逐日再分析值,其中,表示整时段L个时刻内、第m个样本、第k类要素形成的向量,由公式(2)得到;
将步骤1得到的历年逐日海洋统计再分析样本矩阵扣除掉累年逐日海洋再分析统计分析样本矩阵,得到历年逐日的海洋再分析距平样本矩阵,如公式(5)所示:
式中,为历年逐日的海洋再分析距平样本矩阵; 表示第1个样本、第k类要素的再分析距平值,表示整时段L个时刻内、第1个样本、第k类要素形成的向量,由公式(2)得到;σ表示标准化矩阵, 并且, 表示标准化因子,表示其中第nk个空间点,表示tl时刻,k类要素中第nk个空间点的平均值。
4.根据权利要求1所述的一种基于经验正交函数分解的海洋多要素中长期统计预测方法,其特征在于,步骤3中,所述的海洋多要素时空四维正交模态构建方法如下:
对历年逐日的海洋再分析距平样本矩阵进行正交模态向量分解,即其中表示正交模态,A表示模态系数,并且满足公式(6):
式中,I表示全矩阵;
构建相关系数矩阵,如公式(7)所示:
式中,表示历年逐日的海洋再分析距平样本矩阵;Ckk,k=1,2,…,K代表第k类要素和第k类要素的相关系数矩阵,对角线上的元素代表同种要素的时空相关系数矩阵,非对角线元素代表不同要素的时空相关系数矩阵;
根据历年逐日的海洋再分析距平样本矩阵的展开式,相关系数矩阵变为:
式中,表示特征值矩阵;
使用雅可比分解,计算相关系数矩阵的特征向量和特征值,同时等价地得到的特征向量和特征值,如公式(9)所示:
其中,λN表示特征向量;
历年逐日海洋统计再分析样本矩阵X采用特征向量表示为:
公式(10)即为海洋多要素时空四维正交模态。
5.根据权利要求1所述的一种基于经验正交函数分解的海洋多要素中长期统计预测方法,其特征在于,步骤4中,所述的海洋多要素中长期统计预测模型构建为:以步骤3得到的海洋多要素时空四维正交模态为基函数,将待分析和预测海洋的当前和未来状态采用这组基函数进行展开,并通过拟合当前海洋各类要素,计算各个模态的系数,从而实现对海洋多要素中长期统计预测模型的构建。
6.根据权利要求5所述的一种基于经验正交函数分解的海洋多要素中长期统计预测方法,其特征在于,步骤4中,所述的海洋多要素中长期统计预测模型的构建方法具体包括:
设待分析和预测的海洋要素真实场向量用Xt表示,Xt对应的当前时刻之前的一段时间的海洋要素真实场为FXt,其中,F表示投影算符;当前时刻之前一段时间的海洋要素分析称为初始场,用Xb表示,初始场Xb相对于真实场Xt的误差Eb如公式(11)所示:
Eb=Xb-FXt (11)
首先采用历年逐日的海洋再分析距平样本矩阵正交模态向量分解所得的正交模态对初始场Xb进行拟合,构造目标泛函,如公式(12)所示:
式中,Jb(Γ)为目标泛函;Γ表示待优化的正交模态系数;
对目标泛函求极小值,即:
解得海洋多要素中长期统计预测模型构建为:
式中,Xa表示海洋多要素中长期统计预测模型的结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910391323.0A CN110222872B (zh) | 2019-05-12 | 2019-05-12 | 基于经验正交函数分解的海洋多要素中长期统计预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910391323.0A CN110222872B (zh) | 2019-05-12 | 2019-05-12 | 基于经验正交函数分解的海洋多要素中长期统计预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110222872A true CN110222872A (zh) | 2019-09-10 |
CN110222872B CN110222872B (zh) | 2023-04-18 |
Family
ID=67820980
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910391323.0A Active CN110222872B (zh) | 2019-05-12 | 2019-05-12 | 基于经验正交函数分解的海洋多要素中长期统计预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110222872B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110458338A (zh) * | 2019-07-23 | 2019-11-15 | 天津大学 | 一种基于再分析资料的海面高度中长期统计预测方法 |
CN110704549A (zh) * | 2019-10-09 | 2020-01-17 | 中国石油大学(华东) | 海洋环境数据服务粒度选择和构建方法、系统、介质及设备 |
CN110727750A (zh) * | 2019-10-09 | 2020-01-24 | 中国石油大学(华东) | 海洋时空过程对象提取与多尺度数据映射方法、介质及设备 |
CN111159246A (zh) * | 2019-11-15 | 2020-05-15 | 西安海的电子科技有限公司 | 一种基于大数据的海洋生物多样性统计方法 |
CN111815041A (zh) * | 2020-06-28 | 2020-10-23 | 青岛哈船海智科技有限公司 | 一种基于改进eof的海水温度场预测方法 |
CN112580899A (zh) * | 2021-02-23 | 2021-03-30 | 国能日新科技股份有限公司 | 一种融合机器学习模型的中长期风力发电预测方法及系统 |
CN113051795A (zh) * | 2021-03-15 | 2021-06-29 | 哈尔滨工程大学 | 一种面向海上平台保障的三维温盐场分析预报方法 |
CN113052370A (zh) * | 2021-03-15 | 2021-06-29 | 哈尔滨工程大学 | 一种基于时空经验正交函数的海洋环境要素统计预测方法 |
CN113297801A (zh) * | 2021-06-15 | 2021-08-24 | 哈尔滨工程大学 | 一种基于steof-lstm的海洋环境要素预测方法 |
CN114706144A (zh) * | 2022-02-24 | 2022-07-05 | 南方海洋科学与工程广东省实验室(湛江) | 海洋环境预报方法、装置、系统及存储介质 |
CN118378570A (zh) * | 2024-06-25 | 2024-07-23 | 青岛哈尔滨工程大学创新发展中心 | 一种面向长时间尺度环境信息统计分析方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102221714A (zh) * | 2011-03-11 | 2011-10-19 | 钱维宏 | 基于大气变量物理分解的低温雨雪冰冻天气的中期预报系统和方法 |
CN102289571A (zh) * | 2011-07-25 | 2011-12-21 | 浙江大学宁波理工学院 | 基于来能预测的海洋能发电系统的负载优化调度方法 |
WO2013154242A1 (ko) * | 2012-04-09 | 2013-10-17 | 주식회사 싸이트로닉 | 해양환경의 실시간 모니터링을 이용한 해양 구조물의 정적 및 동적 포지셔닝 시스템 및 방법 |
CN106446209A (zh) * | 2016-09-30 | 2017-02-22 | 浙江大学 | 一种海洋环境监测要素短期预测方法 |
CN108319772A (zh) * | 2018-01-26 | 2018-07-24 | 中国科学院海洋研究所 | 一种波浪长期数据的再分析方法 |
CN108981957A (zh) * | 2018-05-31 | 2018-12-11 | 西北工业大学 | 基于自组织神经网络及经验正交函数的水下温度场重构方法 |
-
2019
- 2019-05-12 CN CN201910391323.0A patent/CN110222872B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102221714A (zh) * | 2011-03-11 | 2011-10-19 | 钱维宏 | 基于大气变量物理分解的低温雨雪冰冻天气的中期预报系统和方法 |
CN102289571A (zh) * | 2011-07-25 | 2011-12-21 | 浙江大学宁波理工学院 | 基于来能预测的海洋能发电系统的负载优化调度方法 |
WO2013154242A1 (ko) * | 2012-04-09 | 2013-10-17 | 주식회사 싸이트로닉 | 해양환경의 실시간 모니터링을 이용한 해양 구조물의 정적 및 동적 포지셔닝 시스템 및 방법 |
CN106446209A (zh) * | 2016-09-30 | 2017-02-22 | 浙江大学 | 一种海洋环境监测要素短期预测方法 |
CN108319772A (zh) * | 2018-01-26 | 2018-07-24 | 中国科学院海洋研究所 | 一种波浪长期数据的再分析方法 |
CN108981957A (zh) * | 2018-05-31 | 2018-12-11 | 西北工业大学 | 基于自组织神经网络及经验正交函数的水下温度场重构方法 |
Non-Patent Citations (2)
Title |
---|
QI SHAO等: "Ocean Reanalysis Data-Driven Deep Learning Forecast for Sea Surface Multivariate in the South China Sea", 《AGU》 * |
白杨 等: "基于经验正交函数和机器学习的南海海面高度异常预测", 《海洋通报》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110458338B (zh) * | 2019-07-23 | 2023-01-03 | 天津大学 | 一种基于再分析资料的海面高度中长期统计预测方法 |
CN110458338A (zh) * | 2019-07-23 | 2019-11-15 | 天津大学 | 一种基于再分析资料的海面高度中长期统计预测方法 |
CN110704549A (zh) * | 2019-10-09 | 2020-01-17 | 中国石油大学(华东) | 海洋环境数据服务粒度选择和构建方法、系统、介质及设备 |
CN110727750A (zh) * | 2019-10-09 | 2020-01-24 | 中国石油大学(华东) | 海洋时空过程对象提取与多尺度数据映射方法、介质及设备 |
CN111159246A (zh) * | 2019-11-15 | 2020-05-15 | 西安海的电子科技有限公司 | 一种基于大数据的海洋生物多样性统计方法 |
CN111815041A (zh) * | 2020-06-28 | 2020-10-23 | 青岛哈船海智科技有限公司 | 一种基于改进eof的海水温度场预测方法 |
CN112580899A (zh) * | 2021-02-23 | 2021-03-30 | 国能日新科技股份有限公司 | 一种融合机器学习模型的中长期风力发电预测方法及系统 |
CN113051795A (zh) * | 2021-03-15 | 2021-06-29 | 哈尔滨工程大学 | 一种面向海上平台保障的三维温盐场分析预报方法 |
CN113052370A (zh) * | 2021-03-15 | 2021-06-29 | 哈尔滨工程大学 | 一种基于时空经验正交函数的海洋环境要素统计预测方法 |
CN113052370B (zh) * | 2021-03-15 | 2024-06-14 | 哈尔滨工程大学 | 一种基于时空经验正交函数的海洋环境要素统计预测方法 |
CN113051795B (zh) * | 2021-03-15 | 2023-04-28 | 哈尔滨工程大学 | 一种面向海上平台保障的三维温盐场分析预报方法 |
CN113297801A (zh) * | 2021-06-15 | 2021-08-24 | 哈尔滨工程大学 | 一种基于steof-lstm的海洋环境要素预测方法 |
JP2022191176A (ja) * | 2021-06-15 | 2022-12-27 | 哈爾濱工程大学 | Steof-lstmに基づく海洋環境要素予測方法 |
WO2022262500A1 (zh) * | 2021-06-15 | 2022-12-22 | 哈尔滨工程大学 | 一种基于steof-lstm的海洋环境要素预测方法 |
JP7175415B1 (ja) | 2021-06-15 | 2022-11-18 | 哈爾濱工程大学 | Steof-lstmに基づく海洋環境要素予測方法 |
CN114706144A (zh) * | 2022-02-24 | 2022-07-05 | 南方海洋科学与工程广东省实验室(湛江) | 海洋环境预报方法、装置、系统及存储介质 |
CN118378570A (zh) * | 2024-06-25 | 2024-07-23 | 青岛哈尔滨工程大学创新发展中心 | 一种面向长时间尺度环境信息统计分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110222872B (zh) | 2023-04-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110222872A (zh) | 基于经验正交函数分解的海洋多要素中长期统计预测方法 | |
WO2022194045A1 (zh) | 一种面向海上平台保障的三维温盐场分析预报方法 | |
Thompson et al. | Variability in oceanographic barriers to coral larval dispersal: Do currents shape biodiversity? | |
CN104156728B (zh) | 一种基于栈式编码和softmax的极化SAR图像分类方法 | |
CN111815041B (zh) | 一种基于改进eof的海水温度场预测方法 | |
CN113642475A (zh) | 基于卷积神经网络模型的大西洋飓风强度估算方法 | |
CN112285376A (zh) | 一种基于cnn-lstm的风速预测方法 | |
Stevens et al. | Life history demographic parameter synthesis for exploited Florida and Caribbean coral reef fishes | |
CN104021308A (zh) | 基于ERA-Interim和ERA40预测海浪有效波高的方法 | |
Raes et al. | Modeling framework to estimate and project species distributions space and time | |
CN115238937A (zh) | 一种海洋温盐预测方法及系统 | |
CN116401939A (zh) | 一种基于梯度约束神经网络的北极海冰短期预报方法 | |
CN108549961A (zh) | 一种基于cmip5预估海浪有效波高的方法 | |
CN115186959A (zh) | 一种面向海洋科考调查的海浪风险评估方法 | |
CN113946978A (zh) | 基于LightGBM模型的水下三维温盐并行预报方法 | |
CN117172091A (zh) | 一种时空多尺度的北极海冰多要素中长期智能预报方法 | |
Obakrim et al. | Statistical modeling of the space-time relation between wind and significant wave height | |
CN111273376B (zh) | 降尺度的海表净辐射确定方法、系统、设备及存储介质 | |
CN104021434B (zh) | 基于ERA‑Interim预测海浪有效波高的方法 | |
CN112330029A (zh) | 一种基于多层convLSTM的渔场预测计算方法 | |
CN117290684A (zh) | 基于Transformer的高温干旱天气预警方法、电子设备 | |
CN113052370B (zh) | 一种基于时空经验正交函数的海洋环境要素统计预测方法 | |
CN113435780B (zh) | 一种基于神经网络的应急通信感知装备体系效能评估方法 | |
CN115345207A (zh) | 一种自适应的多气象要素预测方法 | |
Becker et al. | Abundance, distribution, and seasonality of cetaceans within the US exclusive economic zone around the Hawaiian archipelago based on species distribution models |
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 |