基于相参雷达慢扫模式下的浪流信息提取方法
【专利摘要】本发明公开了基于相参雷达慢扫模式下的浪流信息提取方法，包括步骤为：计算选定反演区域中各散射单元的雷达回波多普勒速度，对各径向上的雷达回波多普勒速度求平均，获得各径向上的海流速度；利用各径向上的海流速度，采用最小二乘法获取海面矢量流速和流向信息；建立海浪轨道速度谱和海浪波数谱的映射关系，从而获得海浪波数谱信息。本发明的优势在于：由于浪流信息的提取是基于绝对的多普勒速度信息，相对于传统的基于回波强度的浪流提取方法，测量精度更高，而且可以避免繁琐的标定的工作及标定带来的误差。
【专利说明】
基于相参雷达慢扫模式下的浪流信息提取方法
技术领域
[0001] 本发明涉及海洋科学技术领域，特别是一种基于相参雷达慢扫模式下的浪流信息 提取方法。
【背景技术】
[0002] 在海洋观测技术领域中，微波雷达具有测量精度高、环境干扰小、测量范围大、时 空分辨率高、体积轻便等优势，在海洋科学研究中有着巨大的应用前景。
[0003] 目前，利用X波段传统导航雷达提取浪流信息的技术已经在实践中取得显著成效， 基本原理是X波段雷达向海面发射电磁波，与雷达波长相当的毛细波发生Bragg散射，并受 到长波对短波的阴影、倾斜及流体动力学调制，形成雷达海杂波图像，最后利用雷达海杂波 图像的时间序列信息来提取海浪和海面流速。
[0004]相对于传统的X波段导航雷达，相参雷达不但能够获取各分辨单元散射强度信息 而且还可以提取相位信息，因而可以获取更为丰富的海态信息。
【发明内容】
[0005] 本发明要解决的技术问题是针对上述现有技术的不足，而提供一种基于相参雷达 慢扫模式下的浪流信息提取方法，该基于相参雷达慢扫模式下的浪流信息提取方法，采用 多普勒速度来反演浪流信息，相对于传统的基于回波强度测量算法，测量精度更高。同时在 有效浪高提取算法方面，传统的基于雷达回波强度信息的反演算法，结果需要实测数据进 行标定，而本发明公布的算法是基于雷达回波多普勒速度信息，避免了繁琐的标定工作及 标定带来的误差。
[0006] 为解决上述技术问题，本发明采用的技术方案是：
[0007] 基于相参雷达慢扫模式下的浪流信息提取方法，包括以下步骤：
[0008] 步骤1，计算各散射单元的雷达回波多普勒速度:将相参雷达天线的扫描方式设置 为慢扫模式，并在扫描的雷达回波图像中选定反演区域;然后，利用互相关或FFT算法，计算 选定反演区域中各散射单元的雷达回波多普勒速度。
[0009] 步骤2,计算各径向上的海流速度:将步骤1计算出的各散射单元的雷达回波多普 勒速度，按各径向求平均，获得各径向上的海流速度。
[0010] 步骤3，计算矢量流速和流向信息:利用步骤2计算出的各径向上的海流速度，采用 最小二乘法获取海面矢量流速和流向ig息。
[0011] 步骤4，获取海浪波数谱:将步骤1计算出的各散射单元的雷达回波多普勒速度减 去步骤2计算出的各径向上的海流速度，得到各散射单元的海浪轨道速度;然后对各径向上 的海浪轨道速度分别在距离空间上做傅里叶变换得到各径向上的海浪轨道速度谱;接着， 建立海浪轨道速度谱和海浪波数谱的映射关系，从而获取海浪波数谱。
[0012] 步骤5,海波海态参数推导:结合海波频散关系，利用步骤4获取的海浪波数谱推导 出海波海态参数，海波海态参数包括有效波高、峰波波数、峰波频率。
[0013]所述步骤1中，利用FFT算法，计算选定反演区域中各散射单元的雷达回波多普勒 速度的方法为:首先对散射单元积累的脉冲数做快速傅里叶变换，得到雷达回波的频谱信 息P(f)，剔除干扰频谱信息后，谱能量最大值对应的频率即为多普勒频率fd，则各散射单元 的雷达回波多普勒速度V r,e为：
[0015] 上式中，L为雷达波长。
[0016] 所述步骤2中，雷达回波多普勒速度按各径向求平均后，得到的各径向上的海流速 度公式为：
[0018] 上式中，Ve为各径向上的海流速度；％j为各散射单元的雷达回波多普勒速度;Nr 为选定反演区域所包含的距离元数。
[0019] 所述步骤3中，米用最小二乘法获取海面矢量流速和流向彳目息时，最小二乘法的表 达式为：
[0021]上式中，V。为海面矢量流速;Dc海面流向；Ve为各径向上的海流速度;Νθ为选定反演 区域所包含的径向数;θ」为各径向波束角度。
[0022]所述步骤4中，根据线性波理论，海浪轨道速度谱和海浪波数谱的映射关系为：
[0024] 上式中，Su(k)为海浪轨道速度谱;Sn(k)为海浪波数谱;Θ为径向波束角度。
[0025] w为海浪角频率，满足海波频散关系
其中g，k，h分别代表重力 加速度、波数、水深；
[0026] 0W为海浪波向，D(9W)为海浪方向分布函数，具体为：
[0028]其中，s为能量扩散系数，一般取值为4, θ〇为主波波向；
[0029 ]通过3D-FFT算法，海浪波数谱Sn( k)的最终解析式为：
[0030] Sn(k)=XDSu(k)/〇2
[0032]所述步骤5中，有效波高的计算公式：
[0034]其中，凡表示有效波高，S(w)表示频率谱，S(w)通过以下转化关系式获得：
[0036] 其中，Sn(k)表示海浪波数谱，g，k分别表示重力加速度和海面波数。
[0037] 所述步骤5中，海波海态参数还包括峰波波数，峰波波数的计算公式为：
[0039]其中，1^表示峰波波数，Sn(k)为海浪波数谱，k表示海面波数。
[0040] 所述步骤5中，海波海态参数还包括波峰频率，波峰频率^由海波频散关系
计算获得，其中，w为海浪角频率，g，k，h各代表重力加速度、海面波数、 水深。
[0041] 本发明采用上述方法后，采用多普勒速度来反演浪流信息，相对于传统的基于回 波强度测量算法，测量精度更高。同时在有效浪高提取算法方面，传统的基于雷达回波强度 信息的反演算法，结果需要实测数据进行标定，而本发明公布的算法是基于雷达回波多普 勒速度信息，避免了繁琐的标定工作及标定带来的误差。
【附图说明】
[0042] 图1是本发明一种基于相参雷达慢扫模式下的浪流信息提取方法的流程示意图。
【具体实施方式】
[0043]下面结合附图和具体较佳实施方式对本发明作进一步详细的说明。
[0044] 如图1所示，一种基于相参雷达慢扫模式下的浪流信息提取方法，包括以下步骤。
[0045] 步骤1，计算各散射单元的雷达回波多普勒速度。
[0046] 第一步，选定反演区域。
[0047] 将相参雷达天线的扫描方式设置为慢速扫描模式，简称为慢扫模式，优选设置为 lr/min。然后在扫描的雷达回波图像中选定反演区域，假定角度分辨率为δθ，距离分辨率为 Sr，则选定反演区域所包含的径向数为Ne，包含的距离元数为Nr，各角度积累的径向脉冲数 为NP。
[0048] 第二步，计算雷达回波多普勒速度。
[0049]利用互相关或FFT算法，计算选定反演区域中各散射单元的雷达回波多普勒速度。
[0050] 下面详细介绍利用FFT算法，计算选定反演区域中各散射单元的雷达回波多普勒 速度的方法。
[0051] 首先对散射单元积累的脉冲数做快速傅里叶变换，得到雷达回波的频谱信息P (f)，剔除干扰频谱信息后，谱能量最大值对应的频率即为多普勒频率fd，多普勒频率fd的计 算公式为：
[0053]则各散射单元的雷达回波多普勒速度Vr,e为：
[0055] 上式中，L为雷达波长。
[0056] 步骤2,计算各径向上的海流速度:将步骤1计算出的各散射单元的雷达回波多普 勒速度，按各径向求平均，获得各径向上的海流速度。
[0057] 雷达多普勒速度包含海流速度和海浪轨道速度，海流周期相对于采样时间很大， 可以认为在采样过程中，海流径向速度保持不变，海浪则为随机过程，其各径向速度的平均 值近似为零，对每一径向的各距离单位的雷达回波多普勒速度做平均，即可以得到各径向 上的海流速度公式为：
[0059] 上式中，Ve为各径向上的海流速度；。y为各散射单元的雷达回波多普勒速度;Nr 为选定反演区域所包含的距离元数。
[0060] 步骤3,计算矢量流速和流向信息:利用步骤2计算出的各径向上的海流速度，采用 最小二乘法获取海面矢量流速和流向ig息。
[0061] 米用最小二乘法获取海面矢量流速和流向彳目息时，最小二乘法的表达式为：
[0063] 上式中，V。为海面矢量流速;Dc海面流向；Ve为各径向上的海流速度;Νθ为选定反演 区域所包含的径向数;θ」为各径向波束角度。
[0064] 步骤4,获取海浪波数谱:将步骤1计算出的各散射单元的雷达回波多普勒速度减 去步骤2计算出的各径向上的海流速度，得到各散射单元的海浪轨道速度;然后对海浪轨道 速度分别在距离空间上做傅里叶变换得到各径向上的海浪轨道速度谱;接着，建立海浪轨 道速度谱和海浪波数谱的映射关系，从而获取海浪波数谱。
[0065]根据线性波理论，海浪轨道速度谱和海浪波数谱的映射关系为：
[0067]上式中，Su(k)为海浪轨道速度谱;Sn(k)为海浪波数谱;Θ为雷达波束角度。
[0068] w为海浪角频率，满足海波频散关系
其中g，k，h分辨代表重力 加速度、波数、水深；
[0069] 0W为海浪波向，D(9W)为海浪方向分布函数，具体为：
[0071 ]其中，s为能量扩散系数，一般取值为4, θ〇为主波波向。
[0072]通过3D-FFT算法，海浪波数谱Sn( k)的最终解析式为：
[0073] Sn(k)=XDSu(k)/〇2
[0075]步骤5，海波海态参数推导:结合海波频散关系，利用步骤4获取的海浪波数谱推导 出海波海态参数，海波海态参数包括有效波高、峰波频率、峰波波数。
[0076]上述有效波高的计算公式：
[0078]其中，Hs表示有效波高，S(w)表示频率谱，S(w)通过以下转化关系式获得：
[0080]其中，Sn(k)表示海浪波数谱，g，k分别表示重力加速度和海面波数。
[0081 ]上述峰波波数的计算公式为：
[0083]其中，1^表示峰波波数，Sn(k)为海浪波数谱，k表示海面波数。
[0084] 上述波峰频率叫由海波频散关系
计算获得，其中，w为海浪角频 率，g，k，h各代表重力加速度、海面波数、水深。
[0085] 以上详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中 的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种等同变换，这 些等同变换均属于本发明的保护范围。
【主权项】
1. 一种基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于:包括W下步骤： 步骤1，计算各散射单元的雷达回波多普勒速度:将相参雷达天线的扫描方式设置为慢 扫模式，并在扫描的雷达回波图像中选定反演区域;然后，利用互相关或FFT算法，计算选定 反演区域中各散射单元的雷达回波多普勒速度； 步骤2,计算各径向上的海流速度:将步骤1计算出的各散射单元的雷达回波多普勒速 度，按各径向求平均，获得各径向上的海流速度； 步骤3,计算矢量流速和流向信息:利用步骤2计算出的各径向上的海流速度，采用最小 二乘法获取海面矢量流速和流向信息； 步骤4,获取海浪波数谱:将步骤1计算出的各散射单元的雷达回波多普勒速度减去步 骤2计算出的各径向上的海流速度，得到各散射单元的海浪轨道速度;然后对各径向上的海 浪轨道速度分别在距离空间上做傅里叶变换得到各径向上的海浪轨道速度谱;接着，建立 海浪轨道速度谱和海浪波数谱的映射关系，从而获取海浪波数谱； 步骤5，海态参数推导:结合海波频散关系，利用步骤4获取的海浪波数谱推导出海波海 态参数，海波海态参数包括有效波高、峰波波数和峰波频率。2. 根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于： 所述步骤1中，利用FFT算法，计算选定反演区域中各散射单元的雷达回波多普勒速度的方 法为:首先对散射单元积累的脉冲数做快速傅里叶变换，得到雷达回波的频谱信息P(f)，剔 除干扰频谱信息后，谱能量最大值对应的频率即为多普勒频率fd，则各散射单元的雷达回 波多普勒速度Vr,e为：上式中，L为雷达波长。3. 根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于： 所述步骤2中，雷达回波多普勒速度按各径向求平均后，得到的各径向上的海流速度公式 为：上式中，Ve各径向上的海流速度；为各散射单元的雷达回波多普勒速度;Nr为选定 反演区域所包含的距离元数。4. 根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于： 所述步骤3中，采用最小二乘法获取海面矢量流速和流向信息时，最小二乘法的表达式为：上式中，Vc为海面矢量流速;Dc海面流向；Ve为各径向上的海流速度;Νθ为选定反演区域 所包含的径向数;为各径向波束角度。5. 根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于： 所述步骤4中，根据线性波理论，海浪轨道速度谱和海浪波数谱的映射关系为：上式中，Su化)为海浪轨道速度谱;Sn化)为海浪波数谱;θ为径向波束角度； W为海浪角频率，满足海波频散关系，其中g，k，h分别代表重力加速 度、海面波数、水深;0W为海浪波向，D(0w)为海浪方向分布函数，具体为：其中，S为能量扩散系数，一般取值为4，θ〇为主波波向。6. 根据权利要求5所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于： 通过3D-FFT算法，海浪波数谱Sn化)的最终解析式为：.〇7. 根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于： 所述步骤5中，有效波高的计算公式：其中，出表示有效波高，S(w)表示频率谱，S(w)通过W下转化关系式获得：其中，Sn化)表示海浪波数谱，g，k分别表示重力加速度和海面波数。8. 根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于： 所述步骤5中，峰波波数的计算公式为：其中，kp表示峰波波数，Sn化)为海浪波数谱，k表示波数。9. 根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于： 所述步骤5中，波峰频率wp由海波频散关署计算获得，其中，W为海浪角频 率，其中g，k，h分别代表重力加速度、海面波数和水深。
【文档编号】G01S7/41GK106093936SQ201610766600
【公开日】2016年11月9日
【申请日】2016年8月29日
【发明人】王曙曜, 周涛, 施春荣, 陆小虎, 陈天富, 周骏
【申请人】中船重工鹏力（南京）大气海洋信息系统有限公司, 中船重工鹏力（南京）智能装备系统有限公司, 中国船舶重工集团公司第七二四研究所