一种基于环境轮廓线的岛礁海域极值波浪推算方法技术

本发明专利技术公开了一种基于环境轮廓线的岛礁海域极值波浪推算方法，包括选定岛礁海域目标位置，通过实地测量或数值后报获得长周期时间序列的波高和周期变化值，统计处理得到波高周期联合分布的结果，推算波高极值和周期极值的联合分布结果。本发明专利技术评估的数据结果适合于近岸工程、海洋工程和近岛礁工程领域中涉及极值波浪要素确定的工作，推算方法方便有效。推算方法方便有效。推算方法方便有效。

【技术实现步骤摘要】
一种基于环境轮廓线的岛礁海域极值波浪推算方法


[0001]本专利技术涉及海洋工程
，尤其是一种基于环境轮廓线的岛礁海域极值波浪推算方法。

技术介绍

[0002]当波浪从深水或有限水深区域传播到岛礁海域时，由于岛礁地形的影响，引起波浪破碎和底质摩擦(包括渗透)，波浪衰减特征明显，不同形态的海床底质，如泥沙、淤泥、砂石或珊瑚礁等，导致波浪衰减特征具有较大的差异性，进而影响岛礁海域波浪的极值特征。
[0003]现有的在极值波浪分析中，常用的方法有极值分布模型，韦布尔分布模型，皮尔逊III型分布模型和对数正态分布模型，这几种方法都是通过对波高概率分布进行拟合，再推算不同重现期的波高极值，然而周期极值的选择没有较为确切的模型来进行描述，特别是波高极值和周期极值的评估相对独立，如何建立波高和周期联合关系的模型以及利用该模型进行波高极值和周期极值联合分布的推算方法是一个关键问题，因此，如何定量评估岛礁海域的波浪极值条件，即波高极值和周期极值的联合分布，是近岸工程、海洋工程和近岛礁工程领域中急需解决的问题。

技术实现思路

[0004]本专利技术针对
技术介绍
中的不足，提供了一种基于环境轮廓线的岛礁海域极值波浪推算方法。
[0005]本专利技术为解决上述现象，采用以下技术方案，一种基于环境轮廓线的岛礁海域极值波浪推算方法，推算方法包括如下：
[0006]S1，选定对象：以岛礁海域选定点为对象，通过实地测量或数值后报获得到长周期时间序列的波高和周期变化值，统计处理得到波高周期联合分布的结果；<br/>[0007]S2，计算波高概率分布表达式：采用韦布尔模型进行波高分布的拟合，通过最小二乘法求解拟合经验参数，获得波高的概率分布表达式；
[0008]S3，计算给定波高下周期的条件概率分布表达式：采用对数正态分布模型，通过最小二乘法求解拟合经验参数，获得给定波高下周期的条件概率分布表达式；
[0009]S4，计算波高和周期联合分布的概率密度表达式：基于S1和S2，将波高的概率分布表达式与周期的条件概率分布表达式相乘，获得波高和周期联合分布的概率密度表达式；
[0010]S5，确定界线：给定波陡系数，确定波陡限界线；
[0011]S6，计算波高周期联合分布轮廓线：基于S4采用环境轮廓线法获得不同重现期的极值波高和周期联合分布轮廓线；基于S5的波陡限界线，识别轮廓线中波陡限界线对应的波高和周期组合值：大波高和小周期的组合值，判定这些组合值为不合理，最终获得合理的波高周期联合分布轮廓线。
[0012]作为本专利技术的进一步优选方式，步骤S2中，波高概率服从韦布尔分布，周期条件概率密度函数服从对数正态分布，波高周期联合分布概率为波高概率和周期条件概率的乘
积，波高概率表达式如下
[0013][0014]式(1)中，h表示波高，α、β和γ为分布参数，(1)式为三参数韦布尔分布，若γ＝0，则变为二参数韦布尔分布，这里采用二参数韦布尔分布。
[0015]作为本专利技术的进一步优选方式，步骤S3中，周期条件概率密度表达式如下
[0016][0017]式(2)中的μ和σ，需要通过实测数据资料或数值后报结果得到，其拟合式如下
[0018][0019][0020]在式(3)和式(4)中，a0,a1,a2,b0,b1和b2是拟合参数，采用非线性最小二乘法曲线拟合得到，因此，联合概率分布可以表示为
[0021][0022]作为本专利技术的进一步优选方式，步骤S5中，在联合概率分布表达式确定后，将波高概率分布和周期条件概率分布分别转换为标准正态U参数空间
[0023][0024][0025]接着在U形参数空间中建立规定重现期的圆。以每3小时记录的观测值为例的话，对应重现期100年等高线的圆半径为
[0026][0027]最后将圆形转换为环境参数空间对应的波高和周期轮廓线
[0028][0029][0030]波陡系数定义为
[0031]α＝(h/l)
b
·
tanhkd
ꢀꢀꢀꢀꢀꢀꢀ
(11)
[0032]式(11)中的(h/l)
b
表示为破碎波高和波长比值，k为波数，d为水深。这里波陡系数取值为0.12。
[0033]本专利技术通过基于环境轮廓线推算波浪极值条件来进行推算，以波浪实地测量结果或波浪数值后报结果为输入，获取不同重现期的波高和周期联合分布，评估的数据结果适合于近岸工程、海洋工程和近岛礁工程领域中涉及极值波浪要素确定的工作。
附图说明
[0034]图1为本专利技术的测点布局示意图；
[0035]图2为本专利技术的波高概率分布拟合结果；
[0036]图3为本专利技术的给定波高下周期条件概率拟合结果；
[0037]图4为本专利技术的波高和周期联合分布拟合结果；
[0038]图5为本专利技术的不同重现期的极值波高和周期联合分布轮廓线。
具体实施方式
[0039]下面将结合本专利技术实施例中，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0040]本专利技术提供一种技术方案：一种基于环境轮廓线的岛礁海域极值波浪推算方法，推算方法包括如下：
[0041]S1，选定对象：以岛礁海域选定点为对象，通过实地测量或数值后报获得到长周期时间序列的波高和周期变化值，统计处理得到波高周期联合分布的结果；
[0042]S2，计算波高概率分布表达式：采用韦布尔模型进行波高分布的拟合，通过最小二乘法求解拟合经验参数，获得波高的概率分布表达式；
[0043]S3，计算给定波高下周期的条件概率分布表达式：采用对数正态分布模型，通过最小二乘法求解拟合经验参数，获得给定波高下周期的条件概率分布表达式；
[0044]S4，计算波高和周期联合分布的概率密度表达式：基于S1和S2，将波高的概率分布表达式与周期的条件概率分布表达式相乘，获得波高和周期联合分布的概率密度表达式；
[0045]S5，确定界线：给定波陡系数，确定波陡限界线；
[0046]S6，计算波高周期联合分布轮廓线：基于S4采用环境轮廓线法获得不同重现期的极值波高和周期联合分布轮廓线；基于S5的波陡限界线，识别轮廓线中波陡限界线对应的波高和周期组合值：大波高和小周期的组合值，判定这些组合值为不合理，最终获得合理的波高周期联合分布轮廓线。
[0047]步骤S2中，波高概率服从韦布尔分布，周期条件概率密度函数服从对数正态分布，波高周期联合分布概率为波高概率和周期条件概率的乘积，波高概率表达式如下
[0048][0049]式(1)中，h表示波高，α、β和γ为分布参数，(1)式为三参数韦布尔分布，若γ＝0，则变为二参数韦布尔分布，这里采用二参数韦布尔分布。
[0050]步骤本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于环境轮廓线的岛礁海域极值波浪推算方法，其特征在于，推算步骤包括如下：S1，选定对象：以岛礁海域选定点为对象，通过实地测量或数值后报获得到长周期时间序列的波高和周期变化值，统计处理得到波高周期联合分布的结果；S2，计算波高概率分布表达式：采用韦布尔模型进行波高分布的拟合，通过最小二乘法求解拟合经验参数，获得波高的概率分布表达式；S3，计算给定波高下周期的条件概率分布表达式：采用对数正态分布模型，通过最小二乘法求解拟合经验参数，获得给定波高下周期的条件概率分布表达式；S4，计算波高和周期联合分布的概率密度表达式：基于S1和S2，将波高的概率分布表达式与周期的条件概率分布表达式相乘，获得波高和周期联合分布的概率密度表达式；S5，确定界线：给定波陡系数，确定波陡限界线；S6，计算波高周期联合分布轮廓线：基于S4采用环境轮廓线法获得不同重现期的极值波高和周期联合分布轮廓线；基于S5的波陡限界线，识别轮廓线中波陡限界线对应的波高和周期组合值：大波高和小周期的组合值，判定这些组合值为不合理，最终获得合理的波高周期联合分布轮廓线。2.根据权利要求1所述的一种基于环境轮廓线的岛礁海域极值波浪推算方法，其特征在于，步骤S2中，波高概率服从韦布尔分布，周期条件概率密度函数服从对数正态分布，波高周期联合分布概率为波高概率和周期条件概率的乘积，波高概率表达式如下式(1)中，h表示波高...

【专利技术属性】
技术研发人员：陈文炜，刘小龙，蔡志文，孙泽，叶永林，张海华，许凯玮，宋志勋，郭静，
申请(专利权)人：中国船舶科学研究中心，
类型：发明
国别省市：