天地波雷达回波的直达波抑制

朱道建，皮晓珊，岳显昌，徐心毅

(1.武汉大学电子信息学院,武汉 430072；2.中国电子科学研究院,北京 100041)

天地波雷达回波的直达波抑制

朱道建1，皮晓珊1，岳显昌1，徐心毅2

(1.武汉大学电子信息学院,武汉 430072；2.中国电子科学研究院,北京 100041)

高频天地波雷达是一种天发地收的新体制雷达，其信号传播信道主要包括电离层反射传播信道和海面绕射传播信道。电离层反射造成信号相位污染以及多径效应，接收机所接收到的天波直达波会发生频移和展宽，其强度一般远高于海洋回波，受污染后的直达波将会覆盖部分回波信息，影响海洋回波/海面目标的提取。本文根据天地波雷达直达波特性，提出了基于距离域自适应滤波抑制天地波雷达中的直达波强干扰的方法。该方法选择只含有直达波或直达波旁瓣信号的邻近距离元，作为自适应对消器的输入信号，实现天地波雷达回波信号中直达波及其旁瓣信号的抑制。实测结果表明，该方法能够很好的对直达波进行抑制，有利于海洋回波/海面目标信息的提取。

高频天地波雷达；直达波抑制；自适应滤波；多普勒展宽

0 引 言

高频天地波雷达是一种新体制雷达，与传统的地波雷达和天波雷达相比，有着其独特的优势[1-4]。系统中天波发射的高功率信号经由低层大气进入电离层，被电离层反射后再经由低层大气传播。部分信号照射到海面，经海面散射后以绕射的方式到达地波雷接收机，而也有部分信号直接到达接收天线，成为直达波信号。直达波具有信号能量强，方向性强，其旁瓣影响多个距离元，多普勒频率反映电离层的运动状态以及存在多径现象等特点。基于直达波能量强、方位性稳定的特点，可以将其作为参考信号，用于天地波雷达接收阵列的校正。但同时，直达波作为强干扰，其强大的旁瓣可能淹没附近的慢速目标和海洋信息[5]；受多径效应的影响，不同层反射的直达波及其旁瓣将相互叠加，严重影响回波信号中有用信息的提取。因此，天地波雷达在进行信息提取之前，必须对直达波进行抑制，降低它对回波信号提取的影响。

利用自适应天线进行干扰置零、通过设置参考通道实现旁瓣对消[6-8]是传统的对直达波进行抑制的方法，根据直达波的特性，主要从时域、空域和多普勒域对其加以抑制。时域中直达波及其多径信号抑制的方法有通道均衡和自适应方法[9-11]，譬如利用时域置零法对双基地地波雷达中的直达波进行抑制[12]，该方法简单易行，但当目标回波与直达波有部分重叠时，在抑制了直达波的同时也会造成目标回波能量的损失。在空域通常采用自适应波束形成的方法抑制直达波[13-14]。另外，根据目标子空间和直达波子空间的正交性[15]，利用空间相关矩阵来构造直达波子空间并通过空间投影实现直达波抑制。尽管这些方法都是有效的，但并不能完全解决直达波抑制问题。针对具体情况，需要综合运用这些方法，或者发展新的方法。

在天地波雷达回波多普勒谱中，直达波类似于高频雷达系统的射频干扰，在距离域有很强的相关性。高频雷达中的射频干扰在时域表现为短时性[16]，根据干扰信号的带宽可分为窄带射频干扰(narrowband RFIs，NBRFI)和宽带射频干扰 (wideband RFIs，WBRFI)。其中，窄带射频干扰在时域具有较强的相关性，相关时间长，经过多普勒处理后，其所有距离元数据都有很强的相关性；而宽带射频干扰的时间相关性相对较差，且仅与相邻距离元数据有强的相关性[16-19]。针对窄带射频干扰和宽带射频干扰在距离域的不同特性，C. Pan等人[19]利用自适应滤波技术，通过选择不同的输入信号分别对它们进行处理，结果显示，射频干扰得到有效的抑制，海洋回波得以显现。因此考虑在距离域采用自适应滤波算法对电离层反射直达波及其旁瓣信号进行抑制。

本文取只含有直达波或直达波旁瓣信号的邻近距离元，作为自适应对消器的输入信号，在距离域进行自适应滤波来抑制天地波雷达中的直达波。

1 自适应滤波算法

如图1为自适应对消器的结构框图[20]，对消器的输入为期望信号和输入信号，期望信号为某一待处理距离元的回波信号：

(1)

输入信号为其它距离元的回波数

(2)

其中

(3)

T为转置符号。一般，期望信号中包括目标信号和干扰，而输入信号只由干扰信号组成；并且要求期望信号和输入信号中的干扰信号成分是相关的。

图1 自适应对消器结构框图

定义自适应系统的权矢量为W，其中W=[w1,w2,…,wM]T，通过寻找最优权矢量，从而估计出期望信号中的干扰成分，进而将其消除以实现干扰的抑制。系统的输出为：

(4)

其中WHYin为估计出的期望信号中的干扰成分，H表示共轭转置。

为保证杂波的抑制程度，W须使滤波器的输出功率最小，即：

(5)

由最小梯度法得到最佳权矢量为：

(6)

其中RYY为输入信号Yin的自相关矩阵，PYY为输入信号与期望信号的互相关矩阵。将最优权矢量代入式(4)即可得带去除干扰的信号。

(7)

假设待处理距离元的数据为x(n)，则期望信号矩阵为X=[x(0)，x(1)，…，x(N-1)]，其与输入信号矩阵的互相关矩阵由以下公式估计：

(8)

将式(7)、式(8)代入式(6)即可得到最优权矢量。

2 实测数据分析

2015年1月，武汉大学海态实验室进行了天地波混合组网高频超视距雷达海洋环境观测试验。雷达系统中的天波发射站位于湖北咸宁，包含5台发射系统和电离层环境监测设备，雷达信号波形采用调频连续波FMCW体制；地波接收站位于福建龙海、赤湖和东山，采用小型相控阵接收阵列，每站包括8根鞭天线，单通道接收波束宽度为 120°。天波发射和地波接收系统都配有GPS模块以实现同步控制。具体实验过程及参数设置请参考文献[18]。

图2(a)为2015年1月19日龙海站接收的一个采样周期的回波距离-多普勒谱，此时天波站发射一路天波，频率为12.47 MHz，信号扫频周期为0.125 s，距离分辨率为10 km，采样周期为5 min。采样启动时间为12:31:56，有效的探测距离元数目为80。从探测结果图可以看到，直达波及其旁瓣占据了多个距离元，而且由于存在E层和F层两层反射，两种回波叠加在了一起，由于电离层造成的多普勒频偏的不同，使两层的直达波都对另一层的海洋回波信号造成了遮挡。

利用自适应滤波算法对数据中的直达波进行抑制。具体做法如下：对于第k(3

对比图2(a)和图2(b)可以看到，直达波及其旁瓣信号均得到了很大程度的抑制，被E层直达波旁瓣信号遮蔽的F层反射海洋回波负一阶谱以及E层反射海洋回波的正一阶谱清晰的显现出来。

图2 直达波抑制结果

进一步地，给出典型距离元数据的抑制结果，如图3所示，采用的相干积累时间为64 s。其中图3(a)～(c)是E层反射回波的直达波所在距离元及其邻近距离元数据的多普勒谱，可以看到抑制后噪声基底都有所降低，直达波被抑制了将近30 dB；图3(d)为E层反射的非直达波距离元数据的处理结果，可以看到直达波被抑制了将近40dB，原本被F层反射回波的直达波旁瓣遮蔽的正一阶Bragg峰正常显现，而负一阶Bragg峰的强度以及相对位置在处理前后基本没发生变化。图3(e)～(h)是F层反射回波相应距离元数据的多普勒谱，同样可以看到处理后噪声基底降低，直达波及其旁瓣信号的能量将近30 dB，图3(h)中负一阶谱区附近的直达波旁瓣信号被有效抑制，而原来被噪声基底掩盖的正一阶峰明显显现。

图3 直达波及其旁瓣信号抑制前后对比 (数据采样时间12:31)

下面给出另一组采样数据的处理结果，数据同样来自于2015年1月19日龙海站，天波发射频率也为12.47 MHz，采样起始时间为13:56，图4为直达波及其旁瓣信号抑制前后的距离多普勒谱图。从图4(a)可以看到，F层反射回波的谱展宽较严重，直达波旁瓣呈现宽条带状，而E层回波直达波旁瓣所占频谱宽度较窄；图4(b)是应用本文方法进行直达波抑制后的结果，可以看到抑制效果明显，一阶海洋回波均显现出来。

图4 直达波抑制结果

图5 直达波及其旁瓣信号抑制前后对比(数据采样时间13:56)

图5给出的是相应的E层和F层回波中直达波以及非直达波所在距离元的多普勒谱。图5(a) 、5(c)是直达波抑制结果，可以看到噪声基底降低，直达波能量极大的抑制；图5(b) 、5(d)是直达波旁瓣的抑制结果，可以看到图5(b)中被展宽严重的F层直达波旁瓣淹没的负一阶Bragg峰重新显现，而正一阶Bragg峰的信噪比亦得到提高，而图5(d)在能量极强的直达波旁瓣被有效抑制，抑制能量达30 dB。

分析图3和图5可知，直达波抑制前后海洋回波一阶Bragg峰的相对位置没发生改变，表明本文处理算法未造成回波中的海流信息的损失。

3 结 语

本文提出了基于距离域自适应滤波算法抑制天地波雷达回波中的直达波信号的方法。以只含有直达波或直达波旁瓣信号的邻近距离元，作为自适应对消器的输入信号，在距离域进行自适应滤波来抑制天地波雷达中的直达波及其多径效应。实测数据的处理结果表明，直达波及其旁瓣信号得到了明显有效的抑制，抑制能量在30 dB左右，表现出很好的信噪比，并且抑制前后海洋有回波一阶谱的位置没有变化，证明本文算法是有效的。

另外，实测数据处理表明，对于沿距离域呈条形分布、具有强的距离相关性的直达波及其旁瓣信号，本文方法都有较好的抑制效果。在处理的过程中，需要注意如下几个问题：1) 由于海洋回波仅与相邻距离元回波的距离相关性较高，根据自适应滤波算法对输入信号和期望信号的要求，为避免在抑制直达波的同时损失海洋回波，在选择输入信号时需避开与待处理距离元数据紧邻的距离元回波信号；2)由于直达波仅与相邻一到两个距离元数据有强的相关性，在进行距离域自适应滤波时，可供选用的输入信号有限，从而影响了抑制效果；3) 对处于边界的距离元数据，采用本文的抑制算法时需谨慎选择输入数据；4) 对于电离层污染严重回波信号，在进行直达波抑制后海洋回波仍然很模糊，因而在此之前应进行电离层解污染处理。

[1] R.J.Riddolls.Limits on the Dectec-tion of Low-Doppler Targets by a High Frequency Hybrid Sky-Surface Wave Radar System[C]. IEEE Radar Conference. Rome Italy, 2008: 1-4.

[2] 焦培南, 杨龙泉,凡俊梅.短波天波反射/地波绕射组合新传播模式及其可能应用[J]电波科学学报,2007, 22(5):746-773.

[3] 吴雄斌, 张兰, 柳剑飞. 海洋雷达探测技术综述[J]. 海洋技术学报, 2015, 34(3).

[4] 纪永刚, 张杰, 王祎鸣,等. 浅谈天发地收超视距雷达在南海监测中的应用[J]. 海洋开发与管理, 2013, 30(B12):102-104.

[5] 宋孝果. 天地波高频雷达阵列校准和直达波抑制[D]. 哈尔滨工业大学, 2011.

[6] 冯明月, 何明浩, 郁春来,等. 自适应旁瓣对消抗灵巧噪声干扰效果研究[J]. 现代防御技术, 2014, 42(3).

[7] R.T.Compton,Jr.Adaptivearrayscon-cepts and performance [M]. Englew-ood Cliff, NJ: Prentice-Hall, 1988: 75-153.

[8] 王俊,赵洪立,张守宏,保铮.非合作连续波雷达中存在强直达波和多径杂波的运动目标检测方法[J].电子学报, 2005, 33(3): 419-422.

[9] 单涛, 陶然, 卓智海,等. 一种直达波及其多径干扰的抑制方法:, CN 103344945 B[P]. 2015.

[10]朱家兵，洪一，陶亮.基于自适应分数延迟估计的无源雷达直达波干扰对消方法[J].火控雷达技术，2006, 35(2): 46-50.

[11]王俊，水鹏朗，保铮等.基于分数延迟的外辐射源雷达杂波相消算法[J].西安电子科技大学学报(自然科学版), 2005, 32(3): 378-382.

[12]陈多芳.岸-舰双基地波超视距雷达若干问题研究[D]. 西安：西安电子科技大学博士学位论文，2008, 12.

[13]X.B. Wu, L. Li, Y.X. Shao, Y. Li, T. Guo. Experimental determi-nation of significant waveheight byOSMAR071: Comparison with resu-lts from buoy[J]. Wuhan University Journal of NaturalSciences, 2009, 14(6): 499-504.

[14]G. Fabrizio, F. Colone, L. Lombardo, etal. AdaptiveBeamformingforHigh-FrequencyOver-the-HorizonPassiveRadar[J].IETRadar,Sonarand Navigation, 2009, 3(4): 384-405.

[15]陈多芳，陈伯孝，刘春波等. 基于子空间投影的双基地地波超视距雷达直达波抑制方法[J]. 电子与信息学报, 2008, 30(11): 2702-2705.

[16]万显荣. 高频地波雷达数字接收机设计与抗干扰研究[D]. 武汉大学博士学位论文,2005.

[17]X.R. Wan, F. Chen, H.Y. Ke. Mainbeam cochannel interference suppression by range adaptive processing for HFSWR [J]. IEEE Signal Processing Letter, 2006,13(1):29-32.

[18]H.Zhou B.Y. Wen, S.C. Wu. Dense radio frequency interference suppression in HF radars [J]. IEEE Signal Processing Letter, 2004, 12(5):361-364.

[19]C. Pan, H.Zhou B.Y. Wen. Radio frequency interference suppression in high frequency surface wave radar based on range-domain correlation [J]. Journal of Electromagnetic Waves& Applications, 2013, 27(4): 448-457.

[20]姚天任, 孙洪. 现代数字信号处理[M]. 武汉: 华中科技大学出版社, 1999.

[21]张兰. 天地波混合组网体制海态雷达系统设计与实现[D]. 武汉大学博士学位论文,2015.

A Method on Direct Arrival Wave Suppression for High Frequency Hybrid Sky-surface Wave Radar

ZHU Dao-jian1, PI Xiao-shan1, YUE Xian-chang1, XU Xin-yi2

(1. Electronic Information School, Wuhan University,Wuhan 430072; 2. China Academy of Electronic and Information Technology, Beijing 100041)

High Frequency hybrid sky-surface wave radar is a new mode of radar system with sky wave transmitting and surface wave receiving, whose signal path is mainly composed of the ionosphere reflection propagation and the sea surface diffraction propagation. The modulation by perturbations of ionosphere will result in the frequency shift and band width in the electromagnetic wave and induce multi-path effect in radar signal. The direct arrival wave of the radar is usually decades dB stronger than the sea echoes in intensity, therefore the sea echoes are partly covered by the contaminated direct arrival wave which bring more difficulty to extract the sea echoes from radar signal. Considering the high correlation of the direct arrival wave in range domain in hybrid sky-wave radar, a method based on adaptive filtering in range dimension to suppress the strong direct arrival wave interference is developed in this article. Data of adjacent range cells containing the direct arrival wave or its side-lobes are chosen as the input of the adaptive canceler to achieve the suppression of the direct arrival wave or its side-lobes in radar echoes at each range bin. Results of the processing of radar measurements indicate that this method works to suppress the direct arrival wave echoes and improve significantly the signal-to-noise-ratio of sea echoes.

high-frequency hybrid sky-surface wave radar; direct wave suppression; adaptive filtering; doppler broadening

10.3969/j.issn.1673-5692.2017.01.006

2016-12-30

2017-01-15

“十二五”国家863课题(2012AA091701,2012AA091702)、国家自然科学(61401316)

朱道建(1991—)，男，湖北人，硕士，主要研究方向为海洋动力学及海洋遥感；

E-mail：zhudaojian@whu.edu.cn

皮晓珊(1989—)，女，湖北人，硕士，主要研究方向为无线电海洋遥感；

岳显昌(1976—)，男，辽宁人，副教授，主要研究方向为中高层大气动力学及无线电海洋遥感；

徐心毅(1981—),女，吉林人，博士，主要研究方向为网络空间安全。

P733.6

A

1673-5692(2017)01-031-06