[发明专利]一种基于空间风暴过程的设计波高计算方法

说明书

本发明公开了一种基于空间风暴过程的设计波高计算方法，涉及海洋工程领域，该方法包括确定检测点位所处的海域范围，获取海域范围内所有历史风暴过程的路径数据并确定每个独立风暴过程的发生时间段，确定每个独立风暴过程发生时间段中的有义波高的极大值，并将所有独立风暴过程中的有义波高的极大值形成样本数据集，对样本数据集进行极值分布族参数拟合，得到拟合度最高的分布并确定极值分布函数，根据分布函数计算得到特定回归周期下的设计波高。该方法以风暴过程作为锚定标准，来进行极值波高取样及计算，在物理意义上更加符合极端海况的发生物理机制原因，同时取消了先验分布函数的假设，将分布的拟合结果完全归结于样本的拟合程度。

技术领域

本发明涉及近岸工程和海洋工程技术领域，尤其是一种基于空间风暴过程的设计波高计算方法。

背景技术

在海浪统计中，每个海况的持续时间一般取三小时，统计得到的有义波高一般指的是一个三小时海况中的最大1/3波高的均值。设计海况中的有义波高大小对于海上作业的安全性至关重要，准确估算特定海域的有义波高的回归水平，是该海域近岸工程与海上建筑物设计与建造的前提。目前，国内外海岸及海洋工程设计波高的获取，大多采用概率统计分析的方法，即假设极值波高符合某种理论分布，再基于该理论分布推算出某一重现期的波高设计值，如“100年一遇”、“50年一遇”的重现期有义波高。这种重现期极值标准的实质是用有限的历史数据外延推求出更长时间段内可能出现的极值。

在推算重现期值时，常用的样本选择方法主要有三种：(1)年极值法(AnnualMaximum)，选取观测数据中的年最大值作为样本，在一定检验准则下确定理论分布曲线的最佳线型参数，外延推算至多年一遇极值，该方法由于每年只选一个最大值作为样本，数据利用率不高，所以为了保证结果的精度则要求原始数据资料年限较长。(2)全部观测值法，当关心位置没有长期固定的水文观测设备时，可通过临时观测获取短期(1～3个月)的实测数据，利用坐标转换，再外延计算至多年一遇极值；由于此方法为经验性质的外延，不同的分布选取对计算结果影响较大，且拟合的分布函数容易受到小海况波高的影响。(3)超阈法(Peak Over Threshold)，当观测序列有限时，可以确定某一观测值为下限，即阈值，凡是超过此值的峰值点均列入统计序列，而后通过GPD分布(Generalized Pareto Distribution)进行拟合，这种方法要求的序列长度不及年极值法，但如何确定阈值是其关键，对于样本较少的情况，阈值的选取十分困难。

发明内容

本发明人针对上述问题及技术需求，提出了一种基于空间风暴过程的设计波高计算方法，本发明的技术方案如下：

一种基于空间风暴过程的设计波高计算方法，包括：

确定检测点位所处的海域范围，海域范围覆盖风暴过程的发生路径；

获取海域范围内所有历史风暴过程的路径数据，并确定每个独立风暴过程的发生时间段；

确定每个独立风暴过程的发生时间段中的有义波高的极大值，并将所有独立风暴过程中的有义波高的极大值集合形成样本数据集；

对样本数据集进行极值分布族参数拟合，得到拟合度最高的分布并确定极值分布函数，极值分布族包括若干个不同的分布模型；

根据极值分布函数计算得到特定回归周期下的设计波高。

其进一步的技术方案为，任意两个独立风暴过程的间隔时间不小于3天。

其进一步的技术方案为，确定每个独立风暴过程中的有义波高的极大值，包括：

对位于检测点观测到或计算得到的有义波高按时间序列进行统计；

对统计结果进行比较，得到每个独立风暴过程的发生时间段中的有义波高的极大值。