一種改進的地波雷達鄰近距離單元格一階海雜波對消方法＊

錢文振，紀永剛，王祎鳴，許樂達，戴永壽，于長軍

（1.中國石油大學（華東），山東 青島266555；2.國家海洋局 第一海洋研究所，山東 青島266061；3.哈爾濱工業大學（威海），山東 威海264209）

高頻地波雷達利用垂直極化高頻電磁波沿海面繞射的原理，能夠對海面艦船及低空飛行目標進行超視距探測［1］。在地波雷達回波中會摻雜許多雜波，如海雜波、地雜波、電離層雜波等，對目標檢測形成干擾。地波雷達時域信號經過脈沖壓縮、波束合成及相干積累等信號處理后，可得到距離－多普勒波譜數據。雖然通過長時間的相干積累可以獲得所需的多普勒頻率分辨力，但海雜波、地雜波和電離層等雜波干擾依然存在。特別是海洋回波信號所產生的一階雜波幅度很高，在進行目標檢測時就很容易將這些雜波當做目標被檢測出來，以致于虛警率較高，有時目標還會被海雜波一階譜、二階譜等給淹沒了，造成檢測盲區［2］。所以，在利用波譜數據開展目標檢測前，需要對海雜波、地雜波進行處理，來有效抑制海雜波、地雜波，以降低虛警率、提高探測率。

目前，對這種基于波譜數據的海雜波、地雜波抑制方面的研究相對較少，西安電子科技大學尚海燕利用待檢測“距離和”單元的時頻分布減去鄰近“距離和”單元的平均時頻分布的方法來削弱強雜波［3］，但此方法是假設艦船目標在相干積累時間內沒有距離模糊。本研究針對實驗獲取的高頻地波雷達多普勒數據特性特點，分析了特定相干積累時間內艦船目標的運動特性，特別是定量分析了目標距離模糊特性，通過采用去目標距離模糊的鄰近單元格波譜數據海雜波對消的方法來抑制海雜波，并將本研究海雜波抑制方法應用于目標檢測來驗證其有效性。

1 相干積累時間對目標回波的影響分析

本研究利用某地波雷達站獲取的并經過波束合成及相干積累之后形成的R－D回波譜數據。圖1為發射頻率為4.7MHz時波束1方向的R－D回波譜數據，其相干積累時間為291.23s（重疊率為79.4%），距離單元格數為170，每個距離單元格對應1.5km距離分辨。從圖中可以看出存在目標、海雜波、地物雜波、電離層雜波、噪聲干擾等。圖中標出了不同雜波區域。

在進行多普勒濾波處理時需要對每個波束同一距離單元的回波信號在時域上進行相干積累，在多普勒頻譜上區分運動目標和固定雜波，對目標進行分類和初步識別；對于積累時間的選擇則影響到信噪比的改善、目標回波譜是否展寬以及目標在積累時間內是否運動出距離單元格。在滿足海浪信號平穩性條件下，相干積累時間越長，獲得的信號功率越強，信噪比也越大；然而相干積累時間較長時，目標由于航行方向和速度的變化，其多普勒頻率往往是時變的，信號功率譜有較大的擴展而不能得到有效積累。另外，在進行海上艦船目標探測時，高頻地波雷達為保證足夠的多普勒頻率分辨率其相干積累時間不能太短［4］。針對海上艦船目標采用短相干積累時間更有利于目標探測，一般情況下對艦船目標的積累時間為3～5min。

假定艦船的最大速度為30kn，距離單元格為1.5km，那么此艦船可以在97.2s后跨越此距離單元。積累時間大于97.2s時，進行相干積累處理后，艦船會在這兩個甚至更多的距離單元格內同時出現。針對獲取的波譜數據來看，其積累時間為291.23s，在這段時間內艦船目標（速度30kn）可以航行4.5km，跨越了3個距離單元格，會產生距離模糊。這樣在一個積累時間內，進行數據處理時會在幾個距離單元格內同時出現目標。實測數據處理如圖2所示，標示的目標在多普勒向頻率保持一致，距離向遠離雷達。圖3為在頻率0.145 9，0.149 4，0.163 1Hz位置處的譜線，距離向范圍為37.5～52km，其中可以看出45km 和46.5 km處的回波幅度相差不大，48km處的回波幅度有所下降，分析其原因可能是目標在積累時間內剛駛入該距離單元格，49.5km以后該目標幅度下降且目標信息逐漸消失。由此可以看出該艦船目標在積累時間內跨越了3個距離單元格，且在頻譜范圍上有所展寬。圖4分別是45，46.5，48km處的波譜圖，該目標多普勒頻率均為0.149 4Hz，幅度分別為80.4，80.52，69.6dB。在45km和46.5km 處的回波幅度相差不大，在48km處幅度下降11dB左右。通過對距離向和多普勒向的分析，可得出在相干積累時間內，目標跨越了多個距離單元格，即在R－D波譜圖中目標可以同時在多個距離單元格出現。

圖4 3處距離的多普勒回波譜圖Fig.4 The Doppler echo spectra at 3distance

2 去距離模糊的一階海雜波對消原理

相鄰距離單元格海雜波空時特性有相似性［5］，且在圖1中可以看出不同距離單元上海雜波一階峰空時特性差異較小。若同一波束方位內、相鄰距離單元海雜波多普勒頻率差異小于1～2個多普勒分辨單元（多普勒分辨率取決于計算空時譜所采用的調頻周期數），可采用相鄰距離單元回波數據抑制海雜波［3］。在一個發射綜合波束內，鄰近距離單元海雜波的幅頻特征具有一定的相關性，因此西安電子科技大學尚海燕選擇利用高頻地波雷達海洋回波時頻圖減去平均海雜波時頻分布圖來抑制海雜波［3］。這樣可以削弱強海雜波，使海雜波區以外的目標得以凸顯。由于相干積累時間較長，艦船目標在相干積累時間內可能跨越幾個距離單元格，在選取距離單元時要避開自身的模糊單元。不同距離單元格內多普勒回波數據進行對消，兩對消距離單元格之間相隔3個距離單元。假設船的速度為30kn，這樣在291.23s的積累時間內，艦船能航行4.5 km，所以艦船不可能到達選取的用來對消的單元格內（至少相隔4.5km），從而不會對消掉目標。假如兩對消單元格內都有目標，目標徑向速度、船體大小不同也不會對消掉。對相鄰距離單元回波平均數據進行對消處理，得到一階海雜波抑制掉的幅頻圖。對消處理可以將大部分一階海雜波、地雜波抑制掉，其對電離層雜波抑制也起到一定得作用。

方法流程：首先，選擇某一波束第N個距離單元回波數據作為待檢測距離單元數據；其次，在該波束下，選擇待檢測距離單元鄰近單元格的平均回波數據作為參考對消回波數據；最后，采用待檢測距離單元回波數據減去鄰近距離單元的平均回波數據的方法來抑制海雜波。具體公式為

式中，FN0為第N距離單元格原始回波譜數據；FN＋4，FN＋5，FN－4，FN－5分別為第N＋4，N＋5，N－4，N－5距離單元格回波譜數據；FN為第N距離單元格海雜波對消后結果。

由于進行對消處理后，目標單元格幅度有所下降，為不削弱目標RCS值，需要對其進行幅度補償。將對消前該單元格幅度平均值及對消后幅度平均值做差，求出所需要的幅度補償值，將對消后幅頻圖進行幅度補償，把幅度還原到原始水平。

3 實測數據處理結果與比對

3.1 去距離模糊的一階海雜波對消處理

下面對第28距離單元格的實測回波數據進行海雜波抑制處理，并選取與其相隔3～5個距離單元的平均數據作為參考對消單元數據。將28距離單元格數據減去其對應的參考對消單元數據后再進行幅度補償即可得到海雜波抑制后回波數據。圖5為選取波束1方向的做對消處理的第28距離單元格（42km處）幅頻圖，圖中實線為第42km處的幅頻圖，虛線為鄰近平均距離單元格的幅頻圖。從圖中可以看出，相鄰單元格之間的海雜波有相關性［6］。做對消處理可以使海雜波有較大抑制，而目標的多普勒頻率與海雜波多普勒頻率不同，對消時不會對消掉。利用待檢測距離單元格的原始數據除以鄰近距離單元格的平均數據再對其取對數，即兩單元格原始數據取對數后相減。將第28距離單元格數據減去鄰近平均距離單元格數據，采用波束1、距離42km處，對消結果如下圖所示：

從圖6可以看出，因對消后幅度有所下降，為使目標RCS值恢復原始水平，需要對其進行幅度補償。將對消前該單元格幅度平均值及對消后幅度平均值做差，求出所需要的幅度補償值，將對消后幅頻圖進行幅度補償，把幅度還原到原始水平。進行38.1dB幅度補償后的幅頻圖結果如圖7所示。圖7中實線為對消前第42km處的幅頻圖，虛線為幅度補償后的幅頻圖。根據海雜波的特性分析，海雜波正負一階Bragg峰分別在頻率0.221 5Hz和－0.218 1Hz處，地雜波大體在零頻附近。從圖7可以看出，正負一階峰Bragg得到了較大的抑制，并且地雜波也得到了較大的抑制。而在零頻與正負一階Bragg峰之間頻率為0.180 3Hz，－0.152 8Hz處的目標沒有被抑制掉，且經幅度補償后幅度仍在60dB左右。從圖中可以推斷出，此單元格內可能有兩個目標，一個目標正朝向雷達站運動，另一個目標在遠離雷達站。圖8為對消后局部放大R－D波譜圖，圖9和圖10分別為忽略距離模糊對消結果及局部放大R－D波譜圖，從圖9、圖10中可以看出改進前方法將目標一起對消掉了。通過綜合對比可以得出去距離模糊的鄰近距離單元格海雜波對消法的有效性。

3.2 整體抑制效果比對分析

圖11為波束2、第一時相數據在所有距離單元格上的原始二維圖及對消后二維圖。

圖11 二維多普勒回波圖及對消圖Fig.11 2DDoppler echo spectrum and cancellation spectrum

從圖11中可以看出海雜波對消的整體效果，海雜波一階峰、二階峰以及地雜波得到了較大抑制，而目標卻較明顯的顯現出來。利用海雜波的相關性進行不同距離單元格海雜波對消效果較好。

下面利用常規CFAR檢測算法［7］進行驗證，采用波束2，第一時相數據進行處理。本研究主要針對海雜波、地雜波抑制方法進行研究，目標檢測算法在這不做詳細描述，直接引用其做驗證處理。選取雷達波束2方向，所有距離單元格，對海雜波抑制后的二維數據進行目標檢測處理。圖12為原始二維圖和對消前目標檢測結果的合成圖，圖13為原始二維圖和對消后目標檢測結果的合成圖，圖中標出的小圓圈為目標檢測結果，也就是目標。可以看出，進行鄰近距離單元格雜波對消處理后，海雜波、地雜波得到了有效抑制，然而目標卻都還保留著。從圖中可以明顯的看出，在波束2方向，0～250km范圍內的目標分布情況。當然也不排除有特殊情況，兩對消格內都有目標且速度和大小相同時，此時會對消掉目標。這種情況可以通過在距離單元格選取上的改進而避免，或者用航跡關聯來判斷是否存在目標。

4 結語

本研究針對強海雜波背景特點，開展了基于波譜數據的海雜波特性分析?？紤]到地波雷達相干積累時間相對較長而造成目標距離單元模糊，提出了改進的鄰近距離單元格雜波對消的方法對地波雷達實測波譜數據進行預處理，并用常規目標檢測算法進行了驗證。從目標檢測結果可以看出，檢測前對海雜波進行對消預處理效果比較好，而且是很有必要的。這種預處理有效的避免了海雜波對目標檢測的影響，降低了虛警率。

需要指出的是，在進行對消處理時，淹沒在海雜波一階譜中的目標會被對消掉，從而出現漏警的問題。因此在以后的工作中需要繼續對海雜波一階譜進行特性分析，結合不同海雜波對消方法進行分析，并通過時頻分析來確定淹沒在海雜波一階譜中的目標。

（References）：

［1］GAO X B，ZONG C G，YUAN Y S.Sea clutter canceling for HF ground wave shipborne OTH radar［J］.Acta Electronica Sinica，2000，28（3）：5－6.高興斌，宗成閣，袁業術.高頻地波艦載超視距雷達的海雜波對消［J］.電子學報，2000，28（3）：5－6.

［2］LEI Z Y，WEN B Y，CHENG F.Target detection in first order Bragg peaks with adaptive cancellation［J］.Chinses Journal of Radio Science，2006，21（3）：365－366.雷志勇，文必洋，程豐.基于自適應對消法檢測一階Bragg峰內目標的研究［J］.電波科學學報，2006，21（3）：365－366.

［3］SHANG H Y.The maneuvering target detection of coast－ship bistatic/multistatic surface wave over the horizon radar［D］.Xi＇an：Xi Dian University，2008：42－66.尚海燕.岸艦雙/多基地地波超視距雷達機動目標的檢測［D］.西安：西安電子科技大學，2008：42－66.

［4］ZHOU H.Study on some problems in HF radar signal processing［D］.Wuhan：Wuhan Univ ersity，2004：36－37.周浩.高頻雷達信號處理中的若干問題研究［D］.武漢：武漢大學，2004：36－37.

［5］CHEN D F.Study on some issues of coast－ship bistatic surface wave over－the－horizon radar［D］.Xi＇an：Xi Dian University，2008：95－110.陳多芳.岸－艦雙基地波超視距雷達若干問題研究［D］.西安：西安電子科技大學，2008：95－110.

［6］YU J S，MA H X，ZHU Z Y.Correlation character of radar echo－based sea clutter suppress－ion algorithm［J］.Ship Electronic Engineering.2009，180（6）：121－122.於建生，馬紅星，朱征宇.基于雷達回波相關特征的海雜波抑制方法［J］.艦船電子工程，2009，180（6）：121－123.

［7］LU X L.Enhanced detection of small targets in ocean clutter for high frequency surface wave radar［D］.Canada：University of Victoria，2009：65－91.