第49卷第4期2021年4月
华中科技大学学报（自然科学版）
J. Huazhong Univ. of Sci. &Tech. (Natural Science Edition)
Vol.49 No.4
Apr. 2021
D O I：10.13245/j.hust.210420
基于BP神经网络的局频雷达有效浪局场反演光固化打印
文必洋汤文成田应伟
(武汉大学电子信息学院，湖北武汉430072)
摘要针对传统利用一、二阶峰功率，二阶与一阶峰功率比(R S B)等雷达回波参数勾浪高的简化关系模型反演有 效浪高，存在较大误差的问题，提出一种基于B P(反向传播)神经网络的有效浪高反演方法，获取雷达回波参数与 浪高间更准确的关系，以提高反演精度.首先将上述三种雷达回波参数作为输入，海洋数
值模型提供的有效浪高 作为输出，训练回波参数与浪高的关系模型；然后将雷达回波参数代入训练好的模型中，即可实现有效浪高场的 反演.现场实验表明：雷达在不同方位角/距离元上所反演的浪高时间序列与数值模型有效浪高间的相关系数为微电解填料
0.78〜0.98，均方根误差为0.21〜0.80 m，二者统计均值分别为0.93和0.38 m.与假设R S B和浪高的对应关系是线
性模型的方法相比，该方法具有更高的反演精度.
关键词高频雷达；有效浪高；B P神经网络；一阶峰：二阶峰
中图分类号T N95文献标志码A文章编号1671-4512(2021)04-0114-06
Significant wave height field inversion of high frequency radar based on BP
neural network
WEN Biyang TANG Wencheng TIAN Yingwei
(Electronic and Information School, W u h a n University, W u h a n430072, China)
Abstract Aiming at the problem that there i s comparatively large error in the traditional inversion of the significant w a v e height by using the simplified relationship model of radar echo parameters and w a v e heightsuch as the first-order (Bragg) and second-order peak power, the radio of the second-order to the Bragg (RSB) peak power etc, a significant w a v e height inversion method was proposed based on back propagation (BP) neural network to accurately m a p the relationship between radar echo parameters and wav e height so as to improve the inversion accuracy. First, three parameters of radar echo were used as input and ocean numerical model's significant w a v e height as output to train the relational model between radar echo parameters and w ave height. T h e n，the radar echo parameters were substituted into the trained model, the significant w a v e height field could be retrieved. Field experiment shows that the retrieved w a v e height time series of radar and numerical m odel’s significant wav e height at different azimuth/distance bin s h o w a correlation coefficient of 0.78〜0.98 (mean value 0.93) and a root-mean-square error of 0.21 〜0.80 m (mean value 0.38 m). C o m p a r e d with the method assuming that the relationship between R S B and w a v e height is a linear model* this method has higher inversion accuracy.
K e y w o r d s high frequency radar；significant w a v e height：back propagation (BP) neural network；first-order peak：second- order peak
海洋环境状况的变化影响着人类海上生产生活 的安全，对海况进行实时监测十分必要m.海浪场 作为海洋动力学参数三要素之一，是人类进行海洋遥感的重要观测对象.高频表面波雷达(HFSWR)具有超视距、全天候、低成本[:1等特点，是近岸浪 场遥感的理想设备，而采用单极子交叉环接收天线
收稿日期2020-09-01.
作者简介文必洋(1963-)，男，教授，E-mail: ************.基金项目国家自然科学基金资助项目(41706200).
第4期文必洋，等：基于B P神经网络的高频雷达有效浪高场反演•115.
的紧凑型131HFSWR是其重要发展方向[41.
基于HFSWR的浪高反演理论框架最早是由 Barrick于1972年提出的一、二阶散射截面方程导 出，其主要思想是从海洋回波一阶峰和二阶连续谱 中获取海浪谱信息，进而反演浪高[5_7].但这些方法 都建立在窄波束雷达的前提下，紧凑型HFSWR波束较宽的特性，使上述方法受到限制.为此一些学 者改变思路，利用一、二阶谱峰提取海浪信息.文 献[8]提出利用一阶布拉格峰值功率估计浪高的经 验算法，但由于高海况下布拉格峰饱和与诺顿衰减 因子的影响，限制了算法的浪高探测上界.文献 [9]提出二阶峰与一阶峰功率比(RSB)与浪高的线性 回归模型，提升了浪高探测的上界，减轻了衰减因 子的影响.
但由于二阶峰信噪比(SNR)远低于一阶 峰SNR，因此该算法在低海况和探测距离上受到了 限制.文献[10]理论推导传播距离与方向扩展对 RSB的影响，将改进的线性模型运用于大范围的浪 场反演，使得浪场反演精度进一步提高.
事实上，上述雷达回波参数与浪高之间的关系 往往不是简单的解析关系，采用简单的解析公式量 化该关系往往会带来浪高反演误差因此，本研 究利用BP神经网络较强的学习能力和非线性拟合 能力™，将雷达回波参数中提取的一阶峰功率、二阶峰功率、RSB作为网络输入，欧洲中期气象预报 4^ (European Centre for Medium-Range Weather Fore-casts，ECMWF)数据集[13]提供的有效浪高数据 作为网络输出，训练出回波参数与浪高之间更准确 的关系模型.然后将雷达回波参数代入模型，实现 紧凑型HFSWR浪场的高精度反演.对同一区域不 同时段的浪高场反演实验表明该算法精度较高，优于采用线性模型拟合RSB与浪高关系的方法. 1紧凑型HFSWR浪高反演原理
由雷达方程和海洋表面一二阶散射截面可知：雷达回波多普勒(Doppler)谱中，一阶峰功率、二阶 峰功率及RSB，均会随着海态的增加而增加.因 此，可以通过对三者功率的估计实现浪高反演.三 种功率与浪高的模型可统一表示为
尸=/("s)，
式中：P为一阶峰功率P,、二阶峰功率/>2或1158 强度；Z/s为有效浪高；/(•)为凡和户的对应关系•在已有的研宄中，通常采用线性模型表示RSB强度和//s的关系，采用对数模型表示弋或尺与/^的 关系.一般地，采用雷达估计的回波参数与原位浮标等提供的有效浪高进行拟合，获取对应的拟合模 型，从而实现浪高反演.
在实际中，这些雷达回波参数与浪高之间往往 不是简单的映射关系，难以通过固定的拟合模型得 到111]，而简化的拟合模型会使得反演的浪高产生较 大误差.B P神经网络作为解决数据回归分析问题 的常用工具，因其精度高、适应性强的特性，在业 界己有广泛应用〜.本研宄将通过雷达回波参数和 有效浪高数据训练神经网络，得到回波参数与浪高 间更准确的关系模型，以提高浪场反演精度.
2 BP神经网络拟合模型
2.1BP神经网络构建
浪高场是由一定距离和角度范围的单点浪高共 同组成的.实验中浪场的时间分辨率为0.5 h.BP 神经网络由输入层、隐藏层和输出层三部分构成，隐藏层根据网络复杂程度可选择一层或多层.各层 含有一至多个神经元，相邻层间神经元通过权值和 阈值连接.图1展示了所搭建的浪高反演拟合网络 基本结构.
图1有效浪高反演的BP神经网络基本结构硬质合金锯片铣刀
与现有紧凑型HFSWR浪高反演算法一样，本研究仍然选取回波多普勒谱中的一阶峰功率、二阶 峰功率和RSB作为网络模型输入，由此确定输入 层神经元个数为3.选取有效浪高作为网络模型输 出，确定输出层神经元个数为1.依据经验，BP神经网络隐藏层层数选择为1层，即可有效训练回波 参数与浪高的关系，节点选择依据经验公式有M2)
p = \/ m + n + q,(1)式中：p，m和《分别为隐藏层、输入层和输出层 神经元个数；^为1〜10的常数.通过实验比对，选择隐藏层祌经元个数为6.另外，隐藏层激励函 数为S型函数tansig，输出层激励函数为线性传递 函数 purelin
.
l6562.116.华中科技大学学报（自然科学版）第49卷
2.2网络训练数据集
网络训练的输入参数，采用2013年2月份(共 28 d)于福建下垵采集的高频雷达探海实验数据.实验设备为武汉大学雷达与信号处理实验室研制的 OSMAR-S114%雷达，工作频率约13 MHz，接收天 线采用单极子/交叉环天线|161，天线法向100°(正北 方向起始，顺时针方向).根据回波参数提取情况，选定雷达覆盖区域：角度范围为60°〜140°，其中 角度分辨率为5°；距离范围为2.5〜50 km，距离分 辨率为2.5 km.由此雷达覆盖区域包含的网格单元 数量为17 x 19 = 323.图2给出了实验海域的地理 位置及周边浮标布置情况，雷达站名用XIAN表 示，雷达覆盖区域用红线圈出，图中：儿5, C, D，G为浮标位置；为纬度；£■为经度.将2月份 的雷达实验数据通过MUSIC估角的方式划分到各 网格单元，对空间和时间维度上的坏值进行剔除，缺失值通过线性插值的方法补全，并根据浪场反演 的时间分辨率0.5 h，将各网格单元的时间序列分 成48 x 28 = 1344组作为网络输入训练集.另外，由于实验期间，正多普勒频率区域的回波强度总体 高于负频率，因此训练样本均是从正多普勒频率区 域选择，这利于提高浪高反演的性能.
117.5 117.8 118.1 118.4
E/D
图2实验地理位置及周边浮标分布
网络训练的输出参数是与输入训练数据对应时 刻的有效浪高值.雷达覆盖区域包含323个网格单 元，但该区域内的原位浮标仅有3站(儿C，G).为此，实验利用 ECMWF(wf.int)数据集的历史再分析有效浪高数据作为参考值.它可 以提供原始空间分辨率为0.125°x〇.125°、时间分辨 率为6 h的浪高数据.实验中通过双线性插值的方 法将数据处理成与输入数据对应的时空分辨率，并利用3浮标数据校准ECMWF数据获取校准因子. 校准后的ECMWF数据与浮标C,G对应的浪 高时间序列的相关系数和均方根误差精度可以达到 很高1131，可以作为网络输出训练集.2.3 BP神经网络训练
网络训练采用trainlm作为训练函数.通过信 号在网络中前向传播、误差反向传播的方式，不断 修正网络各层间的权值和阈值，直到达到目标最小 误差.另外，在确定隐藏层神经元数量方面，依据 式(1)，优先选取最大数量的神经元.将输入训练 数据预先进行归一化处理，并随机分成80%和20% 两部分，分别作为训练集和测试集，前者用于训练 网络参数，后者测试网络性能.实验前期由于隐藏 层神经元数量选取过大，出现训练集拟合过好而测 试集表现较差的过拟合现象.通过逐一减少隐藏层 神经元数量，使得训练集和测试集性能表现逐渐接 近；当数量过少时，测试集性能又呈现下降趋势. 综合考虑，选用训练集和测试集性能最接近时的隐 藏层神经元数量6作为最佳值.最后，在确定最终 网络结构
的情况下，将所有输入训练数据作为训练 集，训练神经网络，当达到目标最小误差并保持相 对稳定时结束训练，生成浪高反演网络模型.
3浪场反演结果
将雷达回波参数输入网络模型，即可反演对应 时刻的有效浪高场.利用2013年福建下垵高频雷 达探海实验3月1〜14日长达14d的数据，验证浪 高反演算法的精度和所训练网络的泛化能力.
3.1浪场反演结果和单点浪高的精度
图3展示了 3月2日6:00, 3月4日23:30，3月8日17:30三个时刻雷达反演的浪场(分别对应高、中、低海况)，用于验证本文算法在不同海况条件 下浪场反演的稳定性.明显地，浪高在空间上保持 平稳，三种海况下本文算法反演浪场与ECMWF所提供浪场间的均方根误差分别为0.36, 0.25, 0.28 m，表明浪场反演结果具有较好的稳定性.
将雷达反演的浪高与ECMWF所提供的浪高进 行比对，以进一步分析本文算法反演浪高的性能.从边缘方位角近距离元、中心方位角中间距离元、边缘方位角远距离元，依次选取[方位角，距离元] 为[140。，4]，[100。，10], [70。，16]的3 个网格单 元(分别对应图3中的位置I,II,III)，与ECMWF 浪高时间序列进行比对.另外，基于RSB算法[91反演的浪场被用于反演性能的对比，其中RSB与浪 高之间
的关系采用线性模型表示，该线性模型参数 通过ECMWF浪高与雷达RSB最小二乘拟合获取. 图4〜6分别给出了 3月份雷达反演浪高和ECMWF
浪髙间的时间序列、散点图、误差分布统计结果.
第4期文必洋，等：基于B P神经网络的高频雷达有效浪高场反演•117.
24.2 -
曲嘉瑞
c 24.0 | 23.8 -
-2 麵1
-23.6
■023.4 -
B
117.4 117.8 118.2
"/(〇)
(a)高海况
图3不同海况下本文算法反演浪场(标单位:m)
4
$
(a)位置I
2
"2/m
(a)位置I
(b)位置II
■E C M W F；线性拟合；…神经网络.
图4有效浪高时间序列
(c)位置III
「f夕4
•.A i l1b
-2
111__1
0 2 4 0
"2/m
(b)位置II
线性拟合；•神经网络.
图5有效浪高散点图
2
H2/m
(c)位置m
(a)位置I(b)位置II
一线性拟合；一神经网络
(c)位置m
图6浪高误差分布统计结果
图5中：//，为反演浪高；//2SE C M W F浪高.图6 中：《为统计数量；//,2为反演浪高与ECMWF浪高 的差值.
由图4浪高的时间序列结果分析可以看出：三 个位置处本文算法反演的浪高与ECMWF所提供的 浪高时间序列间具有很高的一致性，其相关系数分 别达到了 0.95, 0.97，0.92，均方根误差分别为 0.32，0.29, 0.41m；采用RSB方法反演的浪高与ECMWF所提供的浪高间的相关系数分别为0.89，0.88，0.87，均方根误差分别为0.46, 0.55，0.53 m.即本文算法在相关系数和均方根误差的指标方 面优于RSB方法.另外，两种方法随着距离增加 和方位角偏离法向，反演数据波动性都随之增大，而RSB方法的波动性要稍大于本文算法.此外，在海况较低的情况下，本文算法反演浪高的结果出
现了明显的下限.主要原因是用于训练网络模型的
•118*华中科技大学学报（自然科学版）第49卷2月份浪高数据下限的影响，通过比对发现训练数
据的下限值普遍高于验证阶段的最低浪高值.
从图5散点图的结果可以发现：在边缘方位角 远距离元位置(位置III处)，两种反演方法散点图的 离散程度都增大，主要原因在于雷达回波数据一、二阶峰受距离衰减和方向扩展等因素的影响，相同 海况下回波功率起伏性大.图中RSB方法反演结 果的离散程度要大于本文算法.
从图6误差统计结果发现：本文算法反演的浪 高与ECMWF所提供的浪高时间序列间的误差主要 集中在±0.4 m内，而RSB方法误差波动范围相对 较大，这也说明本文算法可靠性更高.
3.2有效浪高场的精度
为了反映实验中浪场反演的精度，通过将各网 格单元的反演浪高时间序列一一同ECMWF所提供 的浪高进行比较，计算了二者之间的相关系数和均 方根误差.
图7和8分别绘制了本文方法反演结果的相关 系数与均方根误差的空间分布，另外，RSB算法的 结果也被展示用于对比.本文算法在距雷达2.5〜35 km的范围内(2〜14距离元)反演效果很好，与ECMWF浪高数据的相关系数均在0.9以上，均方 根误差在0.5 m以内；在35〜50 km距离范围（15〜20距离元)且靠近天线法向，相关系数和均方根 误差也保持着较好的水平.统计了反演浪高与ECMWF浪高时间序列的相关系数和均方根误差. 本文算法的总体相关系数分布在0.78〜0.98之间 (86%在0.9以上)，均方根误差分布在0.21〜0.80 m 之间(66%在0.4 m以内)，二者统计均值分别为0.93 和0.38 m;RSB算法总体相关系数分布在0.54〜0.92之间(仅有11%在0.9以上)，均方根误差分布在 0.37〜丨.44 m之间(仅有5
%在0.4 m以内），二者统计 均值分别为0.81和0.67 m.
值得注意的是：当远距离偏离法向时，精度会 略微下降.分析产生上述现象的原因：一方面，较远距离的精度低于较近距离，主要原因是受到距离 衰减因素的影响，信噪比起伏较大，信号估计准确 性降低；另一方面，远距离偏离法向方向的精度低 于法向，主要原因是信号在方位角偏离法向的情况 下，总存在某一个环天线的天线增益在该方位角较 小，信号估角准确性降低.图9给出了较远距离元 RSB时间序列的标准差随角度变化的情况(浪高范 围为l.l±0.25m)，可以发现：在一定浪高范围内，不同距离元的数据波动性都呈现随方位角从法向到 两端而逐渐增大的趋势.
24.2
24.0
Sr 23.8
23.6
23.4
(a)本文算法
N/n卡环弯制
(b)RSB算法
图7浪高时间相关系数空间分布
0.0 0.4 0.8 1.2
24.2
24.0
23.8
23.6
23.4
(a)本文算法
24.2
24.0
23.8
23.6
23.4
117.4 117.8 118.2
N/n
(b)RSB算法
图8浪高均方根误差空间分布(标单位:m)
总体而言，本文算法相对于RSB方法拥有更 高的反演精度和性能；但在远距离偏离法向的反演
性能上，还有待进一步提升.
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议