基于海浪谱的三维非线性海浪模拟方法

1.本发明涉及海洋遥感模拟场景构建技术领域，尤其是涉及一种基于海浪谱的三维非线性海浪模拟方法。


背景技术：

2.在海洋遥感领域，能够获取的实测数据有限，因此需要用模拟数值作补充，为海浪反演等研究提供数据。目前海浪模拟方法主要有基于物理和基于构造的两类方法。基于物理的方法通过求解n-s方程组得到流体质点在各时刻的状态从而生成海浪，仿真结果精确，但计算量大不适合用于产生几百甚至上千平方米范围的模拟海浪。基于构造的方法通过数学模型构造海浪，其中基于海浪谱的方法通过波高的统计特征确定波形，仿真效果逼真且计算量适中，适合用于仿真大范围的模拟海面。
3.基于海浪谱的方法将海浪看作随机过程，把海洋波动近似为由多个振幅不同、频率不同、相位杂乱的线性波形叠加的结果，同时要求海浪能量分布与海浪谱模型吻合。基于海浪谱的方法以谱模型作为海浪波形构建的约束条件，海洋遥感是从图像或回波信号的统计参数中提取海浪信息，二者都关注海浪的统计特性，因此基于谱的海浪模拟方法更适合应用于海洋遥感领域模拟数据的生成。
4.在基于海浪谱的模拟中，以正弦波为组成波的线性叠加模型直观简洁且易于实现，适用于大范围海浪模拟。但实际海浪受多种动力学因素影响，组成波的形状偏离正弦波，因此有多位学者开展了非线性波浪的相关研究。丁平兴等在精确导出任意均匀水域上三维随机波动一阶解、二阶解的基础上，更严密、更合理地重新推导了波面高度分布的偏度与二阶功率谱的理论表达式,并首次给出了有限水深波面二阶谱的理论形式。nin等利用源造波技术对无限水深的完全非线性波浪进行了数值模拟研究。bai和taylor采用高阶边界元方法在数值波浪水槽中模拟造波板实时运动产生波浪，分别对完全非线性规则波、聚焦波与垂直圆柱相互作用问题及完全非线性波浪对固定和漂浮结构物作用问题进行了模拟。波高在统计特性方面的表现为，线性叠加模型的波高服从正态分布，而实际海浪波高则偏离正态分布。kinsman最先使用 gram-charlier级数拟合海浪波面高度的统计分布；longue-higgins解释了这种拟合的依据，从理论上给出了edgeworth形式的gram-charlier级数；孙孚依据线性海浪模式和波动射线理论导出了波高分布；ochi等使用变换关系ζ＝lnh将波面高度ζ所服从的高斯分布转换为有效波高h的对数—正态分布；侯一筠等在变换ζ＝lnh的基础上引入非线性形状参数β，采用变换关系ζ＝lnh/β将其推得的波面高度分布转换为一种非线性有效波高的分布。
5.非线性波的统计特性被证实更接近于真实海浪，因此用非线性的有限振幅波代替微幅波作为组成波模拟波浪，产生了很多研究成果。王岚等通过采用基于jonswap谱的线性滤波并依据β分布特性进行非线性修正的方法对深水无破碎二维非线性海浪进行数值模拟；焦甲龙等在正弦波的表达式中引入一些修正系数,得到一种非线性水波自由面模型，通过确定其修正系数可以快速高效地模拟非线性水波自由面。日本学者tsuchiya，yasud提出
用孤立波作为子波模拟非线性不规则波，并取得了良好的结果，肖波等采用椭圆余弦波作为组成波建立了浅水不规则波的随机模型。在这个研究分支中，采用斯托克斯波作为组成波的模型，因其精准的模拟效果而得到持续性研究。斯托克斯波是非线性重力波的一种波动解，这种波动是斯托克斯于1847年提出的，同时给出了两个关键性的结论：第一，在非线性系统中，上下不对称的周期波列是可能的；第二，弥散关系与振幅有关。瑞利自1876年以后，采用将流场化成定常流场的复势函数方法，通过逐级近似，求得了各阶stokes波弥散关系和波面函数，形成了经典斯托克斯波理论。沈正等1993年给出了深水三阶托斯托克斯波的一种新的近似解，得到了与经典三阶以上斯克斯波不同弥散关系和波面函数，并且证明了与经典的复势解法相比逼近程度更优。薛亚东等2017年提出基于三阶斯托克斯波的海浪数值模拟方法，并通过实测数据验证了模型的有效性。
6.现有的非线性模型虽然在线性波模型的基础上进行了改进，能够满足波高偏离正态分布的特征，但采用的非线性波形状与具体海况无关，导致模拟波面的非线性程度不会随仿真场景设定的改变而变化。而现实中不同场景中，由于海况、水深等因素的影响，不同场景海浪的非线性程度不同，模拟海浪的波形需要与具体的仿真场景设定参数(如风速等)建立联系，因此需对现有模型进行改进。


技术实现要素：

7.本发明的目的是提供一种基于海浪谱的三维非线性海浪模拟方法，模型中含有非线性控制参数，该参数与仿真场景设定的具体参数有关(如风速等)，可以根据仿真场景对波形进行调节，形成与输入参数对应的三维非线性海浪模型，使不同设定场景的模拟海浪具有不同的非线性程度，符合真实海浪的特征。
8.为实现上述目的，本发明提供了一种基于海浪谱的三维非线性海浪模拟方法，具体步骤如下：
9.步骤s1：设置仿真场景，确定海浪频率谱s(ω)和方向分布函数g(ω,θ)类型，确定输入参数u；
10.步骤s2：计算非线性控制参数δ,计算公式如下：其中，α和β均为相对应的海浪频率谱s(ω)的有关常数；
11.步骤s3：计算海浪频率方向谱s(ω,θ)；
12.步骤s4：计算频率ω和波向角θ的分段节点数值和间隔；
13.步骤s5：确定三维非线性海面波高模型η；
14.步骤s6：将场景参数代入三维非线性海面波高模型η中产生模拟三维波高序列。
15.优选的，在步骤s1中，设置仿真场景，步骤s1具体为：
16.步骤s11：根据模拟需求，选择确定类型的海浪频率谱s(ω)和方向分布函数g(ω,θ)，其中，ω表示圆频率，θ表示波向角；
17.步骤s12：根据s11中选择的海浪频率谱s(ω)和方向分布函数g(ω,θ)，设置模型所需的输入参数，输入参数与海浪频率谱s(ω)相对应。
18.优选的，海浪频率谱
19.其中，g＝9.8m/s2,u为输入参数风速,α为表征波高及平均周期有关的常数，β为与平均周期和u有关的常数。
20.优选的，在步骤s2中，计算非线性控制参数δ，步骤s2具体为：
21.步骤s21：根据s12中确定的海浪频率谱s(ω)，计算谱的零阶矩m0，计算公式如下：
22.即
23.步骤s22：计算步骤s12中海浪频率谱s(ω)对应的波数谱s(k)形式，计算公式如下：
[0024][0025]
其中，km＝363rad/m，g＝9.8m/s2；
[0026]
步骤s23：计算波数谱峰对应的波数k
p
，计算方法如下：
[0027]kp
为时k的取值，
[0028]
步骤s24：计算非线性控制参数δ，计算公式如下：
[0029]
δ＝m0·kp
，即
[0030]
优选的，在步骤s3中，根据选择的频率谱s(ω)和方向分布函数g(ω,θ)，构成海浪频率方向谱s(ω,θ)，海浪频率方向谱s(ω,θ)的计算公式如下：
[0031]
s(ω,θ)＝s(ω)
·
g(ω,θ)。
[0032]
优选的，在步骤s4中，计算频率ω和波向角θ的分段节点数值和间隔，步骤s4具体为：
[0033]
步骤s41：根据所选海浪频率谱的能量分布情况，确定海浪能量集中频率范围和波向角范围；
[0034]
步骤s42：根据模拟场景的精细度需求，设置频率分割数n和波向角分割数m；
[0035]
步骤s43：采用等频率分割方法，计算各分割频段频点ωn和频率间隔δω；
[0036]
步骤s44：根据方向分布函数中波向角的分布类型，计算各段波向角θm和波向角间隔δθ。
[0037]
优选的，步骤s5具体为：
[0038]
采用三阶随机斯托克斯波ζ(t)为组成波，结合斜向余弦波函数，构建三维非线性海面波高模型：
[0039][0040]
其中，η表示波浪相对于平静水面的瞬时波高，x和y表示水质点的二维坐标，δ为非
线性控制参数，a
n,m
表示组成波的振幅，ωn表示组成波的频率，θm表示组成波的波向角，ε
n,m
表示组成波的随机初始相位且在0-2π间服从均匀分布，g为表示重力加速度。
[0041]
优选的，步骤s6具体为：
[0042]
步骤s61：计算各组成波振幅，计算公式为：
[0043][0044]
步骤s62：设定观测距离范围和时间序列，计算公式为：
[0045]
xi＝x1 (i-1)
·
dx,i∈n
[0046]
yj＝y1 (j-1)
·
dy,j∈n
[0047]
t
p
＝t1 (p-1)
·
dt,p∈n
[0048]
其中，x1表示模拟海面x轴方向起点坐标，i表示x轴方向节点序号，dx 表示x轴方向节点间距，xi表示x轴方向第i个节点坐标；
[0049]
y1表示模拟海面y轴方向起点坐标，j表示y轴方向节点序号，dy表示y 轴方向节点间距，yj表示y轴方向第j个节点坐标；
[0050]
t1表示模拟起始时刻，p表示时间节点序号，dt表示时间节点间距，t
p
表示第p个时刻；
[0051]
步骤s63：根据模型计算各节点在各时刻的波高值，构成模拟海面，计算公式如下：
[0052][0053]
其中
[0054]
因此，本发明采用上述一种基于海浪谱的三维非线性海浪模拟方法，构建含有非线性控制参数的三维非线性海面波高模型，非线性控制参数与场景输入参数相关，可以根据仿真场景调节组成波波形，模拟不同非线性程度的三维海浪，模拟波高的统计分布偏离正态分布，不同非线性控制参数的模拟结果偏离程度不同，符合实际海浪波高分布规律。
[0055]
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
[0056]
图1为本发明一种基于海浪谱的三维非线性海浪模拟方法流程；
[0057]
图2为本发明不同仿真场景对应的不同非线性控制参数仿真结果对比图。
具体实施方式
[0058]
实施例1
[0059]
图1为本发明一种基于海浪谱的三维非线性海浪模拟方法流程，如图1 所示，一种基于海浪谱的三维非线性海浪模拟方法，具体步骤如下：
[0060]
步骤s1：设置仿真场景。
[0061]
步骤s11：根据模拟需求，选择确定类型的海浪频率谱s(ω)和方向分布函数g(ω,θ)，其中，ω表示圆频率，θ表示波向角。
[0062]
海浪频率谱
[0063]
其中，g＝9.8m/s2,u为输入参数风速,α为表征波高及平均周期有关的常数，β为与平均周期和u有关的常数。
[0064]
仿真场景为深水海域，风区均匀，海浪类型为充分成长的海浪。根据此场景特点，选择海浪频率谱s(ω)为pierson-moskowitz谱(以下称p-m谱)，选择方向分布函数g(ω,θ)为stereo wave observation project(以下称swop) 方向分布函数。其中p-m谱表示为：
[0065][0066]
其中，g＝9.8m/s2，u为风速。
[0067]
swop方向分布函数表示为：
[0068][0069]
其中，ω
p
为海浪频率谱s(ω)取最大值时的频率。
[0070]
步骤s12：根据s11中选择的海浪频率谱s(ω)和方向分布函数g(ω,θ)，设置模型所需的输入参数，选取的海浪频率谱s(ω)不同，输入参数不同，如风速和周期等；
[0071]
本实施例选取p-m谱，根据p-m谱形式，需要设置海面高度19.5米处风速u作为输入参数，本实施例中设置u为10m/s
[0072]
步骤s2中，计算非线性控制参数δ。
[0073]
步骤s21：根据s12中确定的海浪频率谱s(ω)，计算谱的零阶矩m0，计算公式如下；
[0074]
即
[0075]
步骤s22：计算步骤p-m谱s(ω)对应的波数谱s(k)形式，计算公式如下：
[0076][0077]
步骤s23：计算波数谱峰对应的波数k
p
，计算方法如下：
[0078]
时k的取值，在本实施例中计算得：
[0079][0080]
步骤s24：计算非线性控制参数δ，计算公式如下：
[0081]
δ＝m0·kp
，在本实施例中计算得：
[0082][0083]
步骤s3：计算海浪频率方向谱s(ω,θ)。
[0084]
根据模拟需求选择海浪频率谱s(ω)和方向分布函数g(ω,θ)，构成海浪频率方向谱s(ω,θ)，计算公式如下：
[0085]
s(ω,θ)＝s(ω)
·
g(ω,θ)
[0086]
由得谱峰频率ω
p
＝8.5/u
[0087]
将u＝10m/s可得：
[0088][0089]
步骤s4：计算频率ω和波向角θ的分段节点数值和间隔。
[0090]
步骤s41：根据所选海浪频率谱的能量分布情况，确定海浪能量集中频率范围、波向角范围。根据p-m海浪谱特性可知，海浪能量主要集中在圆频率ω∈[0.3,1.8]rad/s范围内，swop方向分布函数的波向角θ的实际变化范围在θ∈[-π/2,π/2]内。
[0091]
步骤s42：根据模拟场景的精细度需求，设置频率分割数n和波向角分割数m，本实施例n＝45，m＝25。
[0092]
步骤s43：采用等频率分割方法，计算各分割频段频点ωn和频率间隔δω，计算过程如下：
[0093]
ω1＝0.3rad/s；
[0094][0095]
ωn＝0.033rad/s
·
n；
[0096]
其中n＝1,2,
…
,45。
[0097]
步骤s44：根据方向分布函数中波向角的分布类型，计算各段波向角θm和波向角间隔δθ。
[0098]
swop方向分布函数的波向角服从均匀分布，即：
[0099][0100]
θm＝0.13rad
·m[0101]
其中m＝1,2,
…
,25。
[0102]
步骤s5：确定三维非线性海面波高模型η。
[0103]
采用三阶随机斯托克斯波ζ(t)为组成波，结合斜向余弦波函数，构建三维非线性
海面波高模型：
[0104][0105]
其中，η表示波浪相对于平静水面的瞬时波高，x和y表示水质点的二维坐标，δ为非线性控制参数，a
n,m
表示组
[0106]
成波的振幅，ωn表示组成波的频率，θm表示组成波的波向角，ε
n,m
表示组成波的随机初始相位且在0-2π间服从均匀分布，g为表示重力加速度。
[0107]
步骤s6：将场景参数代入三维非线性海面波高模型η中产生模拟三维波高序列。
[0108]
步骤s61：计算各组成波振幅，计算公式为：
[0109][0110]
步骤s62：设定观测距离范围和时间序列，计算公式为：
[0111]
xi＝x1 (i-1)
·
dx,i∈n；
[0112]
yj＝y1 (j-1)
·
dy,j∈n；
[0113]
t
p
＝t1 (p-1)
·
dt,p∈n
[0114]
其中，其中，x1表示模拟海面x轴方向起点坐标，i表示x轴方向节点序号，dx表示x轴方向节点间距，xi表示x轴方向第i个节点坐标；
[0115]
y1表示模拟海面y轴方向起点坐标，j表示y轴方向节点序号，dy表示y 轴方向节点间距，yj表示y轴方向第j个节点坐标；
[0116]
t1表示模拟起始时刻，p表示时间节点序号，dt表示时间节点间距，t
p
表示第p个时刻。本实施例中：
[0117]
x＝0:5:1500m
[0118]
y＝0:5:1500m
[0119]
t＝0:0.5:40s
[0120]
步骤s63：根据模型计算各节点在各时刻的波高值，构成模拟海面，计算公式如下：
[0121][0122]
其中
[0123]
每个时刻t对应一个模拟三维海表面，81个时刻构成的时间序列最终将产生81个三维波高表面构成的序列，每个波高表面有301*301个波高值，统计81*301*301个波高值的统计分布，结果如图2(a)所示。
[0124]
实施例2
[0125]
本实施例中设置u为50m/s
[0126]
步骤s2中，计算非线性控制参数δ。
[0127]
步骤s21：根据s12中确定的海浪频率谱s(ω)，计算谱的零阶矩m0，计算公式如下；
[0128][0129]
步骤s22：计算步骤p-m谱s(ω)对应的波数谱s(k)形式，计算公式如下：
[0130][0131]
步骤s23：计算波数谱峰对应的波数kp，计算方法如下：
[0132]
时k的取值，在本实施例中计算得：
[0133][0134]
步骤s24：计算非线性控制参数δ，计算公式如下：
[0135]
δ＝m0
·
kp，在本实施例中计算得：
[0136][0137]
模拟结果的波高统计分布如图2(b)所示。
[0138]
实施例2与实施例1其他步骤一致。
[0139]
因此，本发明采用上述一种基于海浪谱的三维非线性海浪模拟方法，构建含有非线性控制参数的三维非线性海面波高模型，可以根据仿真场景调节非线性控制参数，模拟不同非线性程度的三维海浪，模拟波高的统计分布偏离正态分布，不同控制参数的模拟结果偏离程度不同，符合实际海浪波高分布规律。
[0140]
最后应说明的是：以上实施例仅用以说明本发明的技术方案而非对其进行限制，尽管参照较佳实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对本发明的技术方案进行修改或者等同替换，而这些修改或者等同替换亦不能使修改后的技术方案脱离本发明技术方案的精神和范围。