劉生鋒，嚴(yán) 勇，陸建兵

(1.中國兵器工業(yè)試驗測試研究院， 陜西華陰714200;2.南京電子技術(shù)研究所， 南京210039)

0 引言

雷達(dá)的最大不模糊距離與雷達(dá)的最大不模糊速度是一個矛盾，要增加測速的范圍必然要減少測距的范圍。當(dāng)最大不模糊距離外存在目標(biāo)物時，就會在雷達(dá)回波圖上出現(xiàn)因距離折疊而形成的二(多)次回波［1］。在多普勒天氣雷達(dá)中，尤其是對中小尺度強(qiáng)對流天氣系統(tǒng)進(jìn)行觀測時，為了避免產(chǎn)生速度模糊，往往需要采用比較高的重復(fù)頻率，這時遠(yuǎn)距離回波就會產(chǎn)生折疊、與正?；夭ɑ殳B在一起，使數(shù)據(jù)質(zhì)量大為下降［2-3］。迄今為止，已提出了多種方法解決距離折疊，但各種方法均有不足的地方，最好的方法也只能是減輕它的影響［4］。

在多普勒天氣雷達(dá)中，常用的退距離折疊方法有:批處理方法［5］、分離掃描法［5］、系統(tǒng)相位法［6］和隨機(jī)相位法［7］等。本文詳細(xì)闡述使用相位編碼退距離折疊的原理，并將該技術(shù)應(yīng)用于某X波段多普勒天氣雷達(dá)。實(shí)際探測結(jié)果表明，使用偽隨機(jī)相位編碼可以有效地解決多普勒天氣雷達(dá)的距離折疊問題，但也會在一定程度上使真實(shí)回波有所減弱，并消去邊緣的弱回波。

1 距離折疊的產(chǎn)生機(jī)制

1.1 最大不模糊距離

多普勒天氣雷達(dá)的測距原理與一般雷達(dá)的測距原理相同，即它也是通過測量目標(biāo)物的后向散射波回到雷達(dá)時所需的時間來計算目標(biāo)距離的。因此，多普勒天氣雷達(dá)的最大不模糊距離可以描述為:當(dāng)雷達(dá)發(fā)出一個脈沖遇到該距離處的目標(biāo)物，產(chǎn)生的后向散射波回到雷達(dá)時，下一個雷達(dá)脈沖剛好發(fā)出。該距離可用下式表示

式中:c為光速;PRF為雷達(dá)的重復(fù)頻率。

由于光速c是一個常數(shù)，最大不模糊距離Rmax與脈沖重復(fù)頻率PRF成反比。

1.2 最大不模糊速度

多普勒天氣雷達(dá)測量速度時，通常不是直接測量多普勒頻移，而是通過測量相繼返回的脈沖對之間的相位差來確定目標(biāo)物的徑向速度Vr，見下式

式中:T為脈沖重復(fù)周期;ΔR為相繼脈沖的時間間隔內(nèi)T目標(biāo)物沿徑向變化的距離;λ為雷達(dá)的波長;Δφ為相繼返回的兩個脈沖之間的相位差;PRF為雷達(dá)的重復(fù)頻率。

由于多普勒天氣雷達(dá)的采樣是離散的，即每個脈沖采樣一次，這就造成了Δφ的取值范圍只能在(-π，π)之間。這里，相移π所對應(yīng)的目標(biāo)物徑向速度值稱為最大不模糊速度Vmax，它與脈沖重復(fù)頻率PRF和波長λ的關(guān)系可以通過將π帶入式(2)中得到

對于給定的雷達(dá)，波長λ為一定值，因此最大不模糊速度Vmax與脈沖的重復(fù)頻率PRF成正比。

1.3 最大不模糊距離與最大不模糊速度之間的關(guān)系

將式(1)和式(3)相乘后，可以得到

式(4)就是最大不模糊距離與最大不模糊速度之間的關(guān)系。由此可見，當(dāng)雷達(dá)波長選定后，二者的乘積是一個常數(shù)，也就是說，最大不模糊距離Rmax與最大不模糊速度Vmax是相互制約的，當(dāng)增大其中的一個時，必然會導(dǎo)致另外一個減小;而且波長越短，約束條件越苛刻。

1.4 距離折疊

當(dāng)目標(biāo)物離雷達(dá)的距離大于Rmax且位于kRmax+Rs(k表示任意整數(shù)，Rs為小于Rmax的正數(shù))位置處時，在雷達(dá)回波上該目標(biāo)將出現(xiàn)在距離Rs上，這就是“距離折疊”。距離在Rmax內(nèi)的回波稱為一次回波，距離在Rmax和2Rmax內(nèi)(k=1)的回波稱為二次回波，依次類推，甚至可能存在3次(k=2)或更高次(k＞2)的回波。當(dāng)使用多普勒天氣雷達(dá)觀測中小尺度的強(qiáng)對流天氣［8］時，為了減少速度模糊，常選用較高的PRF，這時Rmax就會減小，導(dǎo)致距離折疊。出現(xiàn)距離折疊時，二次或更高次回波將折疊在一次回波中，導(dǎo)致我們對回波出現(xiàn)的位置產(chǎn)生錯誤判斷并影響產(chǎn)品的數(shù)據(jù)質(zhì)量，圖1給出了一個比較明顯的例子。

圖1a)為高重頻下的強(qiáng)度回波，最大探測距離為75 km;圖1b)為較低重頻下的強(qiáng)度回波，最大探測距離為150 km。由圖1b)可知，在距離雷達(dá)75 km范圍內(nèi)幾乎沒有回波，而在75 km以外，有一片回波;而在圖1a)中，75 km范圍內(nèi)卻出現(xiàn)了回波，這片回波是虛假回波，是由于距離折疊產(chǎn)生的二次回波，嚴(yán)重影響了回波數(shù)據(jù)的質(zhì)量。

圖1 高重頻下的距離折疊

2 隨機(jī)相位編碼退距離模糊的原理

隨機(jī)相位法的基本原理是對不同的脈沖串，在發(fā)射之前對其先隨機(jī)的移一個相角，由于不同距離段所移的相角不同，故能把不同距離段的信號區(qū)分開來。當(dāng)?shù)谝痪嚯x段和第二距離段都存在氣象目標(biāo)時，第n個發(fā)射脈沖的回波信號可表示為

多普勒氣象雷達(dá)是一個使用速調(diào)管發(fā)射機(jī)的全相參雷達(dá)。雷達(dá)由接收機(jī)產(chǎn)生一個優(yōu)化的隨機(jī)相位碼去控制發(fā)射機(jī)的射頻移相器，從而使每一個發(fā)射脈沖的起始相位都是隨機(jī)的，該相位相當(dāng)于式(5)中的φn，當(dāng)我們將接收回波序列進(jìn)行相干時，有

式中:［N］表示所有白噪聲總功率。

由于φn為(-π，π)均勻分布的隨機(jī)相位，因此上述兩式中的和成了白噪聲，從而達(dá)到了第一距離段和第二距離段真實(shí)回波信號的分離。在上述分離方法中，由于其中一個距離段的白化會導(dǎo)致另一距離段的信噪比下降，特別是對于第二距離段而言，由于第一距離段的回波信號較強(qiáng)，白化后使第二距離段信噪比變得很差。因此，在實(shí)際估計時通常采用自適應(yīng)濾波方法。假設(shè)我們要處理第一距離段回波信號，先對第二距離段信號進(jìn)行相干，恢復(fù)其功率譜(由于發(fā)射脈沖具有不同的隨機(jī)初始相位，該回波剛到達(dá)接收機(jī)時是被白化的)后，構(gòu)造一個濾波器，將該相干功率濾除，使其成為一白色譜，然后用IFFT反變換后再與第一距離段進(jìn)行重相干，這樣就可以恢復(fù)第一距離段信號譜，并且該信號譜的信噪比也得到了改善，接下來再按照FFT方法的各種處理步驟計算雷達(dá)回波參數(shù)。

3 應(yīng)用實(shí)例

我們對某X型波段多普勒天氣雷達(dá)中進(jìn)行了改造，根據(jù)隨機(jī)相位編碼退距離模糊的原理，在發(fā)射每個脈沖時，由接收機(jī)產(chǎn)生一個隨機(jī)相位碼(0，π/2，π，3π/2)去控制發(fā)射機(jī)的射頻移相器，從而使每一個發(fā)射脈沖的起始相位都是隨機(jī)的，在接收時再將不同距離段的回波信號進(jìn)行分離。圖2給出了一次實(shí)際探測的結(jié)果。

圖2 使用隨機(jī)相位編碼退距離模糊的對比

圖2a)為量程300 km的強(qiáng)度回波圖，在接近300 km、方位210°附近有一片回波(圖中的A區(qū))，圖2b)為量程150 km、未使用隨機(jī)相位編碼退距離折疊的強(qiáng)度回波圖，因此在B區(qū)產(chǎn)生了圖2a)中A區(qū)回波的二次回波，圖2c)的量程同樣為150 km，但由于使用了隨機(jī)相位編碼，因此在B區(qū)沒有回波，說明該區(qū)域因距離折疊產(chǎn)生的二次回波已經(jīng)被消去。同時，仔細(xì)對比圖2三個圖中的C區(qū)，因為使用了隨機(jī)相位編碼技術(shù)，真實(shí)回波有所減弱，邊緣弱信號更是由于信噪比下降，在低于雷達(dá)檢測門限后無法檢測得到。

4 結(jié)束語

距離折疊是多普勒天氣雷達(dá)中經(jīng)常會出現(xiàn)的現(xiàn)象，嚴(yán)重影響雷達(dá)的數(shù)據(jù)質(zhì)量，甚至?xí)?dǎo)致對天氣形勢的誤判。本文在分析了距離折疊的產(chǎn)生機(jī)制后，詳細(xì)闡述了利用隨機(jī)相位編碼進(jìn)行退距離模糊的原理，并結(jié)合某型X波段多普勒天氣雷達(dá)對該技術(shù)進(jìn)行了成功的應(yīng)用。探測結(jié)果表明，采用隨機(jī)相位編碼技術(shù)可以有效消去由距離折疊而產(chǎn)生的二(多)次回波;但與此同時，也會使真實(shí)回波有所減弱，邊緣弱信號甚至無法檢測到。這是使用隨機(jī)相位編碼退距離折疊需要考慮的問題，也是我們下一步要解決的問題。

［1］ 丁鷺飛，耿富錄.雷達(dá)原理［M］.西安:西安電子科技大學(xué)出版社，2002.Ding Lufei，Geng Fulu.Radar principle［M］.Xi＇an:Xidian University Press，2002.

［2］ 俞小鼎，姚秀萍，熊廷南，等.多普勒天氣雷達(dá)原理與業(yè)務(wù)應(yīng)用［M］.北京:氣象出版社，2006.Yu Xiaoding，Yao Xiuping，Xiong Tingnan，et al.Principle and professional application of Doppler weather radar［M］.Beijing:China Meteorological Press，2006.

［3］ 張靄琛.現(xiàn)代氣象觀測［M］.北京:北京大學(xué)出版社，2008.Zhang Aichen.Modern meteorological observation［M］.Beijing:Peking University Press，2008.

［4］ 潘新民，熊 毅，柴秀梅，等.新一代天氣雷達(dá)退數(shù)據(jù)模糊方法探討［J］.氣象與環(huán)境科學(xué)，2010，33(1):17-23.Pan Xinming，Xiong Yi，Chai Xiumei，et al.Discussion on the CINRAD correcting data ambiguity method［J］.Journal of Meteorologyical and Environment Sciences，2010，33(1):17-23.

［5］ 周紅平，沙雪松，張光峰.批次處理解決氣象雷達(dá)距離模糊［J］. 現(xiàn)代雷達(dá)，2006，28(11):20-21，32.Zhou Hongping，Sha Xuesong，Zhang Guangfeng.Resolving range ambiguity of weather radar using batch processing method［J］.Modern Radar，2006，28(11):20-21，32.

［6］ Torres S M，Zrnic D S.Range and velocity ambiguity mitigation techniques for the WSR-88D weather radar［C］//2004 IEEE International Geoscience and Remote Sensing Symposium Anchorage，AK:IEEE Press，2004:1727-1729.

［7］ 熊 毅.多普勒天氣雷達(dá)中相位編碼退距離模糊的研究［D］.北京:北京郵電大學(xué)，2011.Xiong Yi.Research of phase coding algorithm to mitigate range ambiguity in doppler weather radar［D］.Beijing:Beijing University of Posts and Telecommunications，2001.

［8］ 張培昌，杜秉玉，戴鐵丕.雷達(dá)氣象學(xué)［M］.北京:氣象出版社，2001.Zhang Peichang，Du Bingyu，Dai Tiepi.Radar meteorology［M］.Beijing:China Meteorological Press，2001.