一种通过高频地波雷达反演涌浪波高的方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN109782287A(43)申请公布日2019.05.21(21)申请号201910184143.5(22)申请日2019.03.12(71)申请人国家海洋局北海预报中心地址266061山东省青岛市崂山区云岭路27号(72)发明人王炜荔(74)专利代理机构济南知来知识产权代理事务所(普通合伙)37276代理人于志波(51)Int.Cl.G01S13/95(2006.01)权利要求书2页说明书7页附图1页(54)发明名称一种通过高频地波雷达反演涌浪波高的方法(57)摘要本发明公开一种通过高频地波雷达反演涌浪波高的方法，包括以下步骤：选取用于反演的高频地波雷达回波多普勒谱，对每一个雷达回波多普勒谱进行单独计算确定其当时当地的噪声水平；以能量较大的Bragg单峰作为参照，对雷达回波多普勒谱进行去海流化；筛选出可以进行后续谱形分析的谱；识别谱中涌浪谱峰；确定谱峰的实测位置；对所获谱峰的实测位置进行质量检验；计算获得涌浪角度；计算涌浪波高。本发明利用高频地波雷达观测涌浪不仅是以一个全新的手段提供涌浪观测数据，而且可以有效地辅助该手段提高海面波浪的观测精度，具有重要的实用和科学价值。CN109782287ACN109782287A权利要求书1/2页1.一种通过高频地波雷达反演涌浪波高的方法，其特征在于包括以下步骤：步骤1，选取用于反演的高频地波雷达回波谱，对每一个雷达回波谱进行单独计算确定其当时当地的噪声水平；计算基于白噪声假设，根据白噪声平均功率谱密度与标准差相等的特性有：式中x为噪声的谱值，横杠代表取均值，Nsum为样本数量；计算过程中将每一个回波谱的数据点由小到大进行排列，逐步列入噪声数据进行计算，当达到上述白噪声条件时即得到回波谱中的噪声水平；根据噪声水平检验回波谱，必须具有两个Bragg谱峰；步骤2，以能量较大的Bragg单峰作为参照，对雷达回波谱进行去海流化；回波谱逆多普勒频移移动相同的频率差，使一阶Bragg谱峰频率fB回归理论无海流位置式中k0为雷达电磁波波数，g为重力加速度；步骤3，筛选出可以进行后续谱形分析的回波谱；使用至少10个互相独立的回波谱进行平滑处理，该平均要求所跨时间长度须大于20分钟风浪稳定期，采取对当前回波谱数据于时间上向前向后各一小时取数据做谱平均，以抑制风浪信号，突出涌浪信号；经过叠加后有可能出现奇异值，应检验平滑后Bragg双峰间距的对称性，通常给予两个分辨率(δf)的宽容度；步骤4，识别谱中涌浪谱峰；涌浪谱峰频率(fj)的最大值是在布拉格峰两侧的四个特定的区间寻找；这些区间的选定对应着涌浪频率的设定和回波谱本身特性的限制；首先，涌浪谱峰的寻找必须避开一阶谱的位置；阈值是根据统计的Bragg峰半功率谱宽和涌浪峰半功率谱宽的和来确定的，取6个分辨率长度；||fj|-fB|＞6δf；其次，限定查找范围步骤5，确定谱峰的实测位置；对每一个谱峰的认定都采用五个数据点的方法，然后对这些数据点进行高斯拟合，从而确定谱峰的实测位置；步骤6，对所获谱峰的实测位置进行质量检验；谱峰宽度的阈值取统计数量中的前80％作为有效值，其他作为特异性展宽剔除；步骤7，计算获得涌浪角度；根据涌浪谱峰位置计算得到涌浪角度其中，Δf-和Δf+分-+别是负、正回波谱部分两个涌浪谱峰的多普勒频率差值，即：Δf＝f2-f1；Δf＝f4-f3；2CN109782287A权利要求书2/2页步骤8，计算涌浪波高；采用以下公式计算涌浪有效波高其中，r是实测标准化涌浪谱峰能量为其自身回波谱对频率的能量积分与相邻一阶谱峰能量积分之比；2-2系数C经简化后可以表达为C≈(1+(Ks/k0)/4+m1Kscosθs/k0)；耦合系数Γ可以表达为Γ＝ΓEM-iΓH，其中i是复数流体力学耦合系数电磁耦合系数其中，k0为雷达电磁波波数；K1和K2为共同产生电磁波后向散射共振的两组海浪的波数；箭头代表矢量；m1和m2分别取值±1；ω和ωB为分别对应f和fB的圆频率；海面阻抗Δ在未经特别测定的情况下，可以取平均值在浅水情形下，电磁耦合系数保持不变，但流体力学耦合系数需要进行水深修正：式中d为水深，K′1和K′2是波数修正项：K′1＝K1tanh(K1d)K′2＝K2tanh(K2d)水深修正仅在水深20米-50米的区域进行；将两个波数中的任意一个设想为目标涌浪，并通过数值迭代的方式求得其解，得到涌浪波高。2.根据权利要求1所述的一种通过高频地波雷达反演涌浪波高的方法，其特征在于，步骤1中：设定一阶峰值幅度为高于噪声电平6dB。3.根据权利要求1所述的一种通过高频地波雷达反演涌浪波高的方法，其特征在于，步骤6中：要求在回波谱中同时存在4个涌浪谱峰，二阶谱振幅高于噪声电平3dB。4.根据权利要