基于海用相参雷达测量海浪参数的方法技术

基于海用相参雷达测量海浪参数的方法，本发明专利技术涉及海浪参数测量方法。本发明专利技术的目的为了解决基于非相参雷达计算复杂度大，获得的参数精度低的问题。过程为：一：利用海浪径向速度重构海面波高场；二：对重构海面波高场序列进行子图像选取及标准化；三：对重构海面波高场子图像序列进行三维傅里叶变换获取三维海浪谱；四：对三维海浪谱进行高通滤波，去除噪声；五：对去除噪声后的三维海浪功率谱进行带通滤波，再次去除噪声，获得海浪相关三维海浪谱；六：对获得的海浪相关三维海浪谱进行频率维度积分，获得二维海浪方向谱；七：根据二维海浪方向谱获得海浪参数。本发明专利技术用于海用相参雷达测量领域。

【技术实现步骤摘要】
基于海用相参雷达测量海浪参数的方法
本专利技术涉及海浪参数测量方法。
技术介绍
海浪参数如表面流速，波向及有效波高的测量对海岸保护，航运以及海上工业运营都有着至关重要的影响。准确的海洋环境监测数据可以为上述工作提供有效指导和参考，从而确保海洋资源的合理开采及海上航运的安全性。X波段海用雷达的工作频率为8-12GHz，标准工作波长为3cm。与传统浮标和测量计等单点测量设备相比，X波段雷达具有监控范围广、维护成本低、成像分辨率高的特点。鉴于上述优点，X波段海用雷达在海洋遥感领域的发展尤为迅速，正逐步取代过去几十年间使用的浮标成为目前测量海浪参数的主要工具之一。在利用X波段海用雷达进行海浪测量方向，针对于非相参雷达提出的傅里叶变换算法已经相对成熟，且已经成功应用于一些商业产品。该算法可以从雷达图像序列中提取出海浪谱、表面流速以及有效波高等可以实时观测海况信息的参数。然而由于重力波对雷达回波的调制作用使得海浪雷达图像的强度与海浪波高呈非线性关系，传统三维傅里叶变换算法不可避免地需要使用调制传递函数对雷达图像谱进行校正以获得海浪谱。然而，调制传递函数的形式复杂多样，求取复杂，且受众多因素的影响，导致计算复杂度大，获得的海浪谱参数精确度低。虽然目前已经提出一些调制传递函数的经验公式，但是都存在不同程度的局限性。
技术实现思路
本专利技术的目的为了解决基于非相参雷达计算复杂度大，获得的参数精度低的问题，而提出基于海用相参雷达测量海浪参数的方法。基于海用相参雷达测量海浪参数的方法具体过程为：步骤一：利用海浪径向速度重构海面波高场；步骤二：对重构海面波高场序列进行子图像选取及标准化，得到重构海面波高场子图像序列；步骤三：对重构海面波高场子图像序列进行三维傅里叶变换获取三维海浪谱；步骤四：结合经验关系，对三维海浪谱进行高通滤波，去除噪声，获得去除噪声后的三维海浪功率谱；步骤五：根据频散关系，对步骤四获得的去除噪声后的三维海浪功率谱进行带通滤波，再次去除噪声，获得海浪相关三维海浪谱；步骤六：对步骤五获得的海浪相关三维海浪谱进行频率维度积分，获得二维海浪方向谱；步骤七：根据二维海浪方向谱进一步计算获得海浪参数：海浪频谱、平均角频率、平均周期、峰值角频率、主波方向及有效波高。本专利技术的有益效果为：本专利技术X波段相参雷达相比于非相参雷达，不仅多获得一维速度信息，且可以提供更高的距离分辨率以检测更小尺寸的海浪。由于水质点径向速度与海浪波高之间存在直接联系，因此可以直接使用海浪径向速度重构海面波高场。重构的海面波高场未经过复杂的成像调制，可以取代原始海浪雷达图像，重构海面波高场与原始波高具有更高的相似度，利用重构海面波高场计算的海浪参数具有更高的精确度，这不仅避免了传统基于海用雷达获取海浪参数的频谱分析方法必须要使用调制传递函数校正图像谱的问题，获得海浪参数的整个过程无需使用调制传递函数，降低了计算复杂度，而且可以避免调制传递函数对结果精度的影响以获得更高精度的海浪参数。如图6所示，本专利技术提出的改进算法所获得的海浪频谱比传统算法获得的海浪频普更接近原始曲线。与此同时获得的海浪参数包括平均角频率为0.84弧度/米，偏差0.06弧度/米；平均周期7.51秒，偏差0.51秒；峰值角频率0.74弧度/秒，偏差0.05；主波方向-3.46度，偏差3.46度；有效波高3.38米，偏差0.12米。因此可以证明，本专利技术提出的基于海用相参雷达测量海浪参数的三维傅里叶变换算法具有可行性，且可以获得高精度海浪参数。附图内容图1为基于海用相参雷达测量海浪参数的方法流程图；图2为P-M谱派生的ITTC双参数谱示意图；图3为仿真的t＝0秒时的海浪径向速度场示意图；图4a为t＝0秒时的重构海面波高场示意图；图4b为t＝0秒时的原始海面波高场示意图；图5a为随机方向的重构波高与原始波高比较1图，横坐标为距离，单位为米，纵坐标为波高；图5b为随机方向的重构波高与原始波高比较2图，横坐标为距离，单位为米，纵坐标为波高；图5c为随机方向的重构波高与原始波高比较3图，横坐标为距离，单位为米，纵坐标为波高；图6为传统方法与改进方法获得海浪谱的比较图。具体实施方式具体实施方式一：结合图1说明本实施方式，本实施方式基于海用相参雷达测量海浪参数的方法具体过程为：步骤一：利用海浪径向速度重构海面波高场；步骤二：对重构海面波高场序列进行子图像选取及标准化，得到重构海面波高场子图像序列；步骤三：对重构海面波高场子图像序列进行三维傅里叶变换获取三维海浪谱；步骤四：结合经验关系，对三维海浪谱进行高通滤波，去除噪声，获得去除噪声后的三维海浪功率谱；步骤五：根据频散关系，对步骤四获得的去除噪声后的三维海浪功率谱进行带通滤波，再次去除噪声，获得海浪相关三维海浪谱；步骤六：对步骤五获得的海浪相关三维海浪谱进行频率维度积分，获得二维海浪方向谱；步骤七：根据二维海浪方向谱进一步计算获得海浪参数：海浪频谱、平均角频率、平均周期、峰值角频率、主波方向及有效波高。具体实施方式二：本实施方式与具体实施方式一不同的是：所述步骤一中利用海浪径向速度重构海面波高场；具体过程为：步骤一一、对某一时刻的整个雷达观察范围内的1副海浪径向速度进行距离维度积分获得速度势函数其中，rmin是雷达盲区，u(r,θ,t)是海浪径向速度，r为径向距离，θ为方位角，t为时间；步骤一二、根据小振幅波理论及自由海面的运动学边界条件，确定重构海面波高场为其中，g是重力加速度，ωη是自由表面水质点的垂直运动速度；假设观测海面满足线性波理论且不受岸基浪流作用，即存在时间平稳性和空间均匀性，具有高精度的重构海面波高场可近似为步骤一三、对随时间变化的64幅海浪径向速度场重复步骤一一、步骤一二获得重构海面波高场序列。其它步骤及参数与具体实施方式一相同。具体实施方式三：本实施方式与具体实施方式一或二不同的是：所述步骤二中对重构海面波高场序列进行子图像选取及标准化，得到重构海面波高场子图像序列；具体过程为：由于雷达回波受高的海浪的遮挡，回波会产生一定程度的丢失，这种现象被称为阴影调制。阴影调制使得距离雷达较远处的海浪径向速度丢失或不准确，进而引起重构海面波高的不准确。为了尽量避免阴影调制以及边界区域对海浪参数估计带来的影响，要对重构海面波高场序列做如下操作。步骤二一、对重构海面波高场进行子图像选取；步骤二二、将子图像由极坐标系变换到笛卡尔坐标系下；步骤二三、对子图像进行归一化处理，即将子图像每个位置的波高减去64幅图中该位置波高值构成的时间波高序列的平均值；例如获得64幅图像都是1000×1000的矩阵，每一幅图每一个位置有一个值，第一幅图中比如(1,1)位置有一个值a1，第二幅图同一位置(1,1)位置有一个值a2，第三幅图同一位置(1,1)位置有一个值a3，a1、a2、a3……a64这些序列构成时间波高序列；步骤二四、对64幅重构海面波高场重复步骤二一、步骤二二、步骤二三获得子图像序列。其它步骤及参数与具体实施方式一或二相同。具体实施方式四：本实施方式与具体实施方式一至三之一不同的是：所述对重构海面波高场进行子图像选取，过程为：选取海面原始波高与重构海面波高相关系数大于0.8的矩形区域进行后续处理，确保选取的子图像区域具有良好的重构效果。其它步骤及参数与本文档来自技高网...

【技术保护点】
1.基于海用相参雷达测量海浪参数的方法，其特征在于：所述方法具体过程为：步骤一：利用海浪径向速度重构海面波高场；步骤二：对重构海面波高场序列进行子图像选取及标准化，得到重构海面波高场子图像序列；步骤三：对重构海面波高场子图像序列进行三维傅里叶变换获取三维海浪谱；步骤四：对三维海浪谱进行高通滤波，去除噪声，获得去除噪声后的三维海浪功率谱；步骤五：对步骤四获得的去除噪声后的三维海浪功率谱进行带通滤波，再次去除噪声，获得海浪相关三维海浪谱；步骤六：对步骤五获得的海浪相关三维海浪谱进行频率维度积分，获得二维海浪方向谱；步骤七：根据二维海浪方向谱获得海浪参数：海浪频谱、平均角频率、平均周期、峰值角频率、主波方向和有效波高。

【技术特征摘要】
1.基于海用相参雷达测量海浪参数的方法，其特征在于：所述方法具体过程为：步骤一：利用海浪径向速度重构海面波高场；步骤二：对重构海面波高场序列进行子图像选取及标准化，得到重构海面波高场子图像序列；步骤三：对重构海面波高场子图像序列进行三维傅里叶变换获取三维海浪谱；步骤四：对三维海浪谱进行高通滤波，去除噪声，获得去除噪声后的三维海浪功率谱；步骤五：对步骤四获得的去除噪声后的三维海浪功率谱进行带通滤波，再次去除噪声，获得海浪相关三维海浪谱；步骤六：对步骤五获得的海浪相关三维海浪谱进行频率维度积分，获得二维海浪方向谱；步骤七：根据二维海浪方向谱获得海浪参数：海浪频谱、平均角频率、平均周期、峰值角频率、主波方向和有效波高。2.根据权利要求1所述基于海用相参雷达测量海浪参数的方法，其特征在于：所述步骤一中利用海浪径向速度重构海面波高场；具体过程为：步骤一一、对某一时刻的整个雷达观察范围内的1副海浪径向速度进行距离维度积分获得速度势函数其中，rmin是雷达盲区，u(r,θ,t)是海浪径向速度，r为径向距离，θ为方位角，t为时间；步骤一二、根据小振幅波理论及自由海面的运动学边界条件，确定重构海面波高场为其中，g是重力加速度，ωη是自由表面水质点的垂直运动速度；重构海面波高场近似为步骤一三、对随时间变化的64幅海浪径向速度场重复步骤一一、步骤一二获得重构海面波高场序列。3.根据权利要求1或2所述基于海用相参雷达测量海浪参数的方法，其特征在于：所述步骤二中对重构海面波高场序列进行子图像选取及标准化，得到重构海面波高场子图像序列；具体过程为：步骤二一、对重构海面波高场进行子图像选取；步骤二二、将子图像由极坐标系变换到笛卡尔坐标系下；步骤二三、对子图像进行归一化处理，即将子图像每个位置的波高减去64幅图中该位置波高值构成的时间波高序列的平均值；步骤二四、对64幅重构海面波高场重复步骤步骤二一、步骤二二、步骤二三获得子图像序列。4.根据权利要求3所述基于海用相参雷达测量海浪参数的方法，其特征在于：所述对重构海面波高场进行子图像选取，过程为：选取海面原始波高与重构海面波高相关系数大于0.8的矩形区域。5.根据权利要求4所述基于海用相参雷达测量海浪参数的方法，其特征在于：所述步骤三中对重构海面波高场子图像序列进行三维傅...

【专利技术属性】
技术研发人员：谢俊好，李美薇，姚国伟，席堃，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江,23