[发明专利]一种海上搜救目标漂移路径的预测方法

权利要求书

1.一种海上搜救目标漂移路径的预测方法，其特征在于该方法包括以下步骤：

S₁、收集历史气象数据，建立气象数据库；

S₂、利用大气模型调用S₁数据库中的相关数据，计算获得海上遇险事故发生地点周围海区的风场；

S₃、获取遇险事故发生地点周围海区的现场风实时观测数据，更新S₁建立的气象数据库；对现场风观测数据进行数据同化，提高S₂中大气模型预报的风场的精确度；

S₄、结合S₂风场结果，采用海浪数值模型预报遇险事故发生地点周围海区的波浪场；

S₅、基于TPXO全球潮汐卫星高度计反演数据集提供的8个分潮M₂、S₂、N₂、K₂、K₁、O₁、P₁、Q₁的调和常数给出潮汐与环流的开边界条件；

S₆、以S₂大气模型提供的风场、S₄海浪数值模型提供的波浪场为强迫场，并结合S₅提供的潮汐和环流开边界条件，使用海流数值模型预报包含风生流和潮流的三维流场；

S₇、基于S₂大气模型提供的风场、S₄海浪数值模型提供的波浪场以及S₆海流数值模型提供的三维流场，构建风漂速度计算模型、海上搜救目标漂移轨迹计算模型、误差传递模型和搜救半径计算模型实现对海上遇险人员或船舶的动态漂移轨迹预测；

S₈、基于以上计算结果，利用可视化技术将搜救目标漂移路径预测结果直观显示。

2.根据权利要求1所述的海上搜救目标漂移路径的预测方法，其特征在于，S₂所述的大气模型为WRF，用该模型完成遇险事故发生地点周围海区的未来48小时的风场计算；该模型采用完全可压缩非静力欧拉方程组，水平网格采用Arakawa C网格，垂直坐标采用基于质量的地形追随η坐标。

3.根据权利要求1所述的海上搜救目标漂移路径的预测方法，其特征在于：S₃所述的获取遇险事故发生地点周围海区的现场风实时观测数据的方法为用下载工具定时自动下载观测站服务器提供的常规地面观测数据、常规探空数据、船舶观测数据、卫星海面风观测数据与卫星辐射观测数据。

4.根据权利要求1所述的海上搜救目标漂移路径的预测方法，其特征在于：S₃所述的数据同化的方法为3DVAR方法，为预报用的海流数值模型提供初始场与时变边界条件。

5.根据权利要求1所述的海上搜救目标漂移路径的预测方法，其特征在于S₆所述的三维流场计算方法为：

利用海浪数值模型SWAN来计算搜救相关海区的海浪分布情况及变化，利用大气模型每十五分钟输出一次风场并驱动海浪数值模型，海浪谱在每一个格点上取36个频率和12个方向，通过对源函数项和传播项的积分，求得未来48小时内每小时输出一次的波浪场；

基于S₂大气模型提供的风场、S₄海浪数值模型提供的波浪场以及来自TPXO全球潮汐卫星高度计反演数据集提供的8个分潮M₂、S₂、N₂、K₂、K₁、O₁、P₁、Q₁的调和常数给出潮汐与环流的开边界条件，在WRF模型和SWAN模型运算完毕后，SHELL脚本自动运行海流数值模型FVCOM，并预报获取遇险事故发生地点周围海区的未来48小时的三维流场。

6.根据权利要求1所述的海上搜救目标漂移路径的预测方法，其特征在于：S₇所述的风漂速度计算模型、海上搜救目标漂移轨迹计算模型、误差传递模型和搜救半径计算模型，其构建方法包括：

（1）风漂速度计算模型

根据海上搜救目标风漂移速度的计算方法，结合海上搜救目标的状态和几何特性，建立不同类型的海上搜救目标的风漂速度计算模型：

V_w=C_d·W

其中V_w为风漂速度，W为海表面10米风速，C_d为风拖曳系数，该系数根据搜救目标的排水状态和在水中的浸没比例等因素变化而变化。当搜救目标为船舶时，C_d在1%和5%之间。一般情况下，C_d取1%～3%，遇险人员取1%，救生筏取3%，无动力船只取2%；

（2）海上搜救目标漂移轨迹计算模型

海上搜救目标漂移轨迹计算模型为：

S=S0+∫t0t0+ΔtVtdt+S′]]>

其中，S₀为搜救目标的初始位置，V_t为t时刻搜救目标的漂移速度，为空间和时间的函数，S′为因为湍流影响发生的随机走动距离。搜救目标以速度V_t经时间步长Δt后漂移至S，该模型采用了拉格朗日追踪法，其核心就是求解V_t；这里V_t是各个海洋环境动力要素作用过程所产生的分速度的合成。本发明认为作为预测搜救目标的V_t主要由风漂速度、表面流速和波余流流速构成；

（3）误差传递模型和搜救半径计算模型

（a）误差传递模型

搜救目标的漂移轨迹求解方法如下：

S=Σi=1k(SW+SC)=Σi=1k(Cd·W[i]+V)·t]]>

这里，S为搜救目标的漂移位移，S_W和S_C分别代表风和流作用下的位移，C_d为风拖曳系数，W代表海表面10米风速，V代表流速（包括表面流速和波余流流速），t代表漂移时间；

预测的轨迹是分段线性叠加计算，最终的预测轨迹误差也是由各阶段产生的误差合成而来；每个阶段的测量可以看作独立，总误差由各个阶段产生的误差按线性方式进行合成；对于每个阶段，风、流的误差对最后预测轨迹的误差影响是风、流两者误差的叠加；

以下首先建立每个阶段的误差计算模型：

由误差理论，当间接测量量y是直接测量量x_i（i=1，2，3…n）的线性函数时，即：

y=a0+Σi=1naixi]]>

这里α₀和α_i均为非随机但是未观测到参数，y的方差D_y为：

Dy=Σi=1nai2Dxi+2Σi<jaiajσxiσxj]]>

当x₁，x₂,…x_n之间不相关时，间接测量量y的方差为：

Dy=Σi=1nai2Dxi]]>

其均方根差为：

σy=[Σi=1nai2σxi2]1/2]]>

当间接测量量y是直接测量量x_i（i=1，2,3…n）的非线性函数时，则：

y=f(x_i)（i=1，2，3…n）

其方差为：

Dy=Σi=1nai2Dxi+2Σi<jaiajσxixj]]>

这里，

ai=(∂f∂xi)x=mx]]>

ρ_ij是x_i，x_j的相关系数；

当x₁，x₂,…x_n之间不相关时，ρ_ij=0，间接测量量的方差为：

Dy=Σi=1n(∂f∂xi)x=m2]]>Dxi=Σi=1naiσxi2]]>

其均方根差为：

σy=[Σi=1nai2σxi2]1/2]]>

若y₁，y₂,…y_m均为几个间接测量量x₁，x₂,…x_n的线性函数，即：

y_r=f(x_ri)（r=1，2，3…m）

则可用向量和矩阵符号表示为：

Y=A₀+AX

其中：

Y=y1y2···ym,]]>X=x1x2···xn,]]>A0=a10a20···am0,]]>A=[ari]=a11a12···a1na21a22···a2n···am1am2···amn]]>

则其误差向量为：

σ_Y=Aσ

其中，σ=(σ₁，σ₂，…σ_n)^T

若y₁，y₂,…y_m均为几个间接测量量x₁，x₂,…x_n的非线性函数，即：

y_r=f(x_ri)（r=1，2，3…m）

则其误差向量为：

σ_Y=Aσ

其中，σ=(σ₁，σ₂，…σ_n)^T，

对于漂浮物轨迹预测问题，若不考虑风拖曳系数C_d的误差影响，直接测量量为风、流速度分量W_U，W_V，V_U，V_V，间接测量量为预测轨迹位置的北向坐标Y_N和东向坐标Y_E，则直接测量量对间接测量量的误差传递规律满足线性函数的规律。

此时，m=2，n=4，

则：Y=y1y2=YEYN,]]>X=x1x2x3x4=WUWVVUVV,]]>A₀=0，A=[ari]=a11a12a13a14a21a22a23a24]]>

对每一阶段k(k=1，2，…K)，a_ij(i=1，2；j=1，2，3，4)为：

a₁₁[k]=a₂₂[k]=C_d·t，a₁₃[k]=a₂₄[k]=t，a₁₂[k]=a₁₄[k]=a₂₁[k]=a₂₃[k]=0，即：

A=Cd·t0t00Cd·t0t]]>

则其误差向量为：

σ_Y=Aσ

其中，σ=(σWU,σWV,σVU,σVV)T]]>

若考虑风中漂移系数C_d的误差影响，则直接测量量V_U，V_V对间接测量量的误差传递规律满足线性函数的规律，而直接测量量W_U，W_V及C_d对间接测量量的误差传递规律则应按非线性函数的传递规律进行计算。根据相应的误差理论，此时的系数矩阵A的元素应为相应的偏导数，则：

Y=y1y2=YEYN,]]>X=x1x2x3x4x5=WUWVVUVVCd,]]>A₀=0，A=[ari]=a11a12a13a14a15a21a22a23a24a25]]>

对每一阶段k(k=1，2，…K)，a_ij(i=1，2；j=1，2，3，4，5)为：

a₁₁[k]=a₂₂[k]=C_d·t，a₁₃[k]=a₂₄[k]=t，a₁₂[k]=a₁₄[k]=a₂₁[k]=a₂₃[k]=0，即：

A=Cd·t0t0WU·t0Cd·t0tWV·t]]>

则其误差向量为：

σ_Y=Aσ

其中，σ=(σWU,σWV,σVV,σVV,σCd)T]]>

由此，建立了每个阶段的误差计算模型；

设每个阶段的误差为σ_i(i=1，2，…k)，则最终的预测轨迹误差应按线性方式合成，即：

σTotal=[Σi=1Kσi]1/2]]>

（b）搜救半径计算模型

通过搜救目标在指定时刻处在半径R圆内的概率P可求出R，搜救半径计算模型为：

P=1-exp[-R2ΔTr22T]]]>

其中：T为从推算开始时起的物体漂移时间；r为ΔT时间的推算均方差（根据风和流的预报误差计算）；ΔT为时段，风和流的速度和方向在这一时段内是常量；

由此，根据给定的概率P就可计算出搜寻半径R。