汪旭東,胡志群,劉 浩,郝立勇
(1.安徽四創(chuàng)電子股份有限公司,安徽合肥230088;2.中國氣象科學(xué)研究院災(zāi)害天氣國家重點實驗室,北京100081)
雙線偏振雷達(dá)通過同時或交替發(fā)射水平和垂直偏振波,并接收兩個偏振方向的回波信號的方法,可同時探測到降水系統(tǒng)的回波強度ZH、差反射率因子ZDR、差傳播相移率KDP和水平垂直信號相關(guān)系數(shù)ρHV,這些量直接反映了降水系統(tǒng)粒子相態(tài)、滴譜分布等微物理結(jié)構(gòu)的變化規(guī)律,從而可以明顯提高雷達(dá)資料質(zhì)量、降水估測和相態(tài)識別的能力[1-6]。
但是,雙線偏振雷達(dá)探測的偏振量(如差分反射率因子ZDR、差傳播相移率KDP)是水平通道回波強度和傳播相位的差別,一般來說,ZH的取值在-10 dBz到65 dBz之間,而ZDR在-5 dB到5 dB之間變化,KDP通常也在-2°/km到5°/km之間變化,相對于ZH而言,它們的取值范圍小很多;另一方面,這些量非常小的誤差就會引起降水強度和相態(tài)識別結(jié)果很大的誤差。這就要求雙線偏振雷達(dá)探測的精度非常高。為了得到準(zhǔn)確可靠的雙線偏振雷達(dá)資料,對系統(tǒng)的測試和訂正就顯得非常重要。Doviak[1]對研制的S波段雙線偏振雷達(dá)在工作模式、天線結(jié)構(gòu)、數(shù)據(jù)質(zhì)量等方面進行了深入的分析。美國在WSR-88D上發(fā)展的KOUN雙線偏振雷達(dá)系統(tǒng)在應(yīng)用前,也對其進行了嚴(yán)格的測試,并將其測試結(jié)果應(yīng)用到資料的分析中[4]。
我國早在上世紀(jì)80年代就利用數(shù)字化雷達(dá)改造的雙線偏振雷達(dá)開展了降水估測、降水粒子相態(tài)識別等方面的研究工作[7-8]。在新一代天氣雷達(dá)建設(shè)的同時,帶有多普勒功能的雙線偏振雷達(dá)研制及其應(yīng)用工作得到了各方面的重視,雷達(dá)廠家研制了不同體制的雙線偏振雷達(dá)系統(tǒng),我國已經(jīng)有多種波長、多種工作模式的雙線偏振雷達(dá)先后投入外場試驗和科研工作,并在應(yīng)用研究方面取得了不少研究成果[9-11]。但迄今為止,雙線偏振天氣雷達(dá)目前在我國尚無統(tǒng)一的測試內(nèi)容、測試方法和性能評估的標(biāo)準(zhǔn),傳統(tǒng)的方法是以我國新一代多普勒天氣雷達(dá)的測試項目和測試方法的標(biāo)準(zhǔn)為依據(jù),在此基礎(chǔ)上增加有關(guān)偏振方面的測試內(nèi)容。如對天線和接收機分別測兩個通道的特性,以兩個通道一致性的好壞作為評價指標(biāo),這種方法無疑是正確的,偏振雷達(dá)出廠前都作這方面測試,接收單位也是以此作為驗收的依據(jù)。但由于出廠測試結(jié)果多是各個分機性能測試結(jié)果數(shù)據(jù)的集成,缺乏對雷達(dá)總的性能指標(biāo)的測試(如包括天線在內(nèi)的性能),且由于運輸和現(xiàn)場安裝原因,雷達(dá)在現(xiàn)場安裝起來后,有些分機的性能可能有所變化,再加上天線罩的影響等,這些都會引起雷達(dá)參數(shù)的改變,如在雷達(dá)氣象方程里還采用出廠前的數(shù)據(jù),必然會帶來誤差,這可能就是目前我國出現(xiàn)不同雷達(dá)探測同一個目標(biāo),其強度信號差別較大的一個重要原因。
為此,安徽四創(chuàng)電子股份有限公司在新一代天氣雷達(dá)測試定標(biāo)的基礎(chǔ)上,綜合研究了雙線偏振雷達(dá)的測試、定標(biāo)和偏振信息的訂正和預(yù)處理等方面的工作,并在其為北京市氣象局研制的3836雙線偏振多普勒雷達(dá)中進行了實際的驗證,其研究結(jié)果可供今后雙線偏振雷達(dá)研制工作參考。
在多普勒天氣雷達(dá)中,為保證雷達(dá)自身的工作參數(shù)穩(wěn)定,給出的結(jié)果數(shù)據(jù)更為可信,必須隨時檢測雷達(dá)的主要參數(shù),以便對探測結(jié)果進行必要的修正,這就要求雷達(dá)系統(tǒng)具有對主要工作參數(shù)進行自動檢測和對探測結(jié)果進行修訂的功能,這個過程就是雷達(dá)的標(biāo)定。雷達(dá)可進行自動或人工干預(yù)檢測和標(biāo)定,其檢測和標(biāo)定內(nèi)容為:發(fā)射機脈沖功率在線檢測、天饋駐波的檢測、雙接收通道一致性標(biāo)校、接收特性曲線標(biāo)校、強度測試、速度驗證、噪聲系數(shù)測試、相位噪聲檢測以及天線波束指向標(biāo)定。
根據(jù)不同類型降水粒子的雙線偏振雷達(dá)探測量的取值范圍和雷達(dá)探測能力,以及各國的雙線偏振雷達(dá)實際的探測精度,一般認(rèn)為ZH,ZDR和KDP的誤差應(yīng)分別小于1.0 dB,0.2 dB和0.2°/km。為了保證這一精度要求,除多普勒雷達(dá)測試和標(biāo)定內(nèi)容外,還要根據(jù)雙線偏振雷達(dá)探測的偏振量的種類和特點采用其他不同的測試方法。在實際的雷達(dá)測試中,本文提出了結(jié)合機內(nèi)機外儀表測試、標(biāo)定和利用氣象觀測目標(biāo)相結(jié)合的綜合測試方法,其中以機內(nèi)機外儀表測試、標(biāo)定為主,氣象和非氣象目標(biāo)測量為輔。目的是確定雙線偏振雷達(dá)偏振量的測量標(biāo)定方法,提高探測精度,減小系統(tǒng)誤差對氣象產(chǎn)品反演的影響,提高資料的可靠性。
雙線偏振雷達(dá)的天饋線系統(tǒng)由天線和饋線兩部分組成,天線部分由共軸雙模喇叭口輻射器、正交模耦合器、圓拋物面反射體組成,饋線部分由定向耦合器、隔離器、諧波濾波器、三態(tài)極化開關(guān)、環(huán)流器、方位向雙通道高功率旋轉(zhuǎn)關(guān)節(jié)、俯仰向高功率旋轉(zhuǎn)關(guān)節(jié)、TR管、固態(tài)噪聲源、PIN開關(guān)組件、波導(dǎo)同軸轉(zhuǎn)換及若干波導(dǎo)組成。對天饋線系統(tǒng)的標(biāo)定測試主要是對天線系統(tǒng)水平、垂直通道一致性測試,即通過天線發(fā)射出去的水平、垂直偏振波的波瓣寬度、主要旁瓣的一致性。兩通道的增益一致性是保證水平垂直偏振波的后向散射能力來自同一降水散射體的重要因素,是偏振量可靠性的重要保證。在雷達(dá)出場前,都要在微波暗室里對天線進行測試。
目前利用外部信號源的測試方法逐漸應(yīng)用到雷達(dá)安裝現(xiàn)場的測試中,由于該方法利用雷達(dá)接收外部標(biāo)準(zhǔn)信號源的信號,所以得到的接收數(shù)據(jù)更能反映雷達(dá)的整體性能和探測能力,這對要求有更高探測數(shù)據(jù)精度的雙線偏振雷達(dá)尤為重要。其測試的原理就是將信號源和標(biāo)準(zhǔn)喇叭安放在合適的陣地上,開機后將信號源發(fā)射喇叭對準(zhǔn)雷達(dá)天線,調(diào)整信號源喇叭使其能準(zhǔn)確地發(fā)射水平(或垂直)線偏振波,斷開雷達(dá)接收機前端與混頻器的連線,將其接到與信號源發(fā)射頻率相同的頻譜儀的輸入端。測試過程中,通過調(diào)整雷達(dá)天線,以雷達(dá)對準(zhǔn)信號源時的位置為中心,在方位角和仰角偏離一定范圍內(nèi)測量雷達(dá)接收到的信號功率。然后計算水平和垂直通道的電軸方向、波瓣寬度、旁瓣電平、隔離度、方向圖的一致性和相關(guān)性等性能指標(biāo)。圖1~4給出了2012年5月在北京3836雙線偏振雷達(dá)的一次測試結(jié)果。測試結(jié)果如下(頻率5 350 MHz):雷達(dá)天線水平極化波束寬度為0.92°(H面)、0.95°(E面);垂直極化波束寬度為0.91°(H 面)、1.00°(E面);旁瓣電平24~28.5 d B;水平通道E面水平發(fā)垂直收交叉極化隔離度為-28.17 d B,其他優(yōu)于-30 d B。
圖2 水平通道交叉極化方向圖對比
圖3 垂直通道交叉極化方向圖對比
圖4 3 dB寬度波束
該雷達(dá)發(fā)射機發(fā)射的功率通過功分器分到水平和垂直通道,進入到兩個通道的功率的比值是計算ZH、ZDR的必要因子。對“雙發(fā)雙收”體制來說,功分開關(guān)的安裝位置越靠近饋源位置越容易實現(xiàn)發(fā)射通道的幅相平衡,理論上安裝在天線的背面為最佳方案,因為功分器至饋源之間的饋線連接器件,短時間內(nèi)幅相起伏很小,但選用該方式,就必須采用三通道旋轉(zhuǎn)鉸鏈,然而目前三通道鉸鏈的電性能,尤其是其在旋轉(zhuǎn)過程中插損和相位的起伏很難保證雙線偏振天氣雷達(dá)的探測精度要求;該雷達(dá)只采用雙通道鉸鏈,將功分開關(guān)組件安裝在方位鉸鏈與發(fā)射管之間,這對雙通道饋線系統(tǒng)的插損和相位平衡的控制造成很大的難度,一般采用對稱性設(shè)計來保證,則此時影響偏振參數(shù)測量結(jié)果的主要器件為三極化開關(guān)、雙通道高功率鉸鏈、正交模耦合器等。
對兩個通道功率差的測試采用在發(fā)射口外接信號輸入源,信號經(jīng)過極化開關(guān)、高功率鉸鏈(方位、俯仰)和正交模耦合器后在饋源的出口處測量其兩個通道的差值。測試結(jié)果:兩個通道的功率差為0.3 d B(5 350 MHz)。
該雷達(dá)采用同時收發(fā)體制,兩個通道一致性測試包括通道對功率的損耗差和兩個通道雷達(dá)波的相位差。通道損耗差與ZDR的計算有關(guān),而從發(fā)射機到天線饋源水平、垂直通道的相位差關(guān)系到發(fā)射的電磁波是否為與水平方向成45°夾角的線偏振波。
兩個通道一致性測試的方法是將接收機內(nèi)部的DDS測試信號放大后接入數(shù)控衰減器,數(shù)控衰減器可同時輸出兩路功率相等的經(jīng)過衰減后信號分別至兩個接收前端??刂茢?shù)控衰減器的衰減碼,使注入到接收前端的測試信號逐漸降低,兩路測試信號分別經(jīng)過接收前端、數(shù)字中頻和信號處理器后送入數(shù)據(jù)處理與顯示分系統(tǒng),雷達(dá)終端監(jiān)控軟件對信號處理輸出信號進行采集并計算,最終在雷達(dá)終端軟件上記錄并顯示出通道增益差曲線。圖5給出了雙通道一致性標(biāo)校原理框圖。
圖5 雙通道一致性標(biāo)校原理框圖
測量結(jié)果:水平垂直通道的整體損耗差0.39 dB,饋源處兩個通道發(fā)射雷達(dá)波的相位差小于2°。
偏振雷達(dá)的數(shù)據(jù)由于受硬件條件的限制,對于弱回波的探測往往誤差較大,為了考慮判別偏振數(shù)據(jù)精度情況,可以此作為判別偏振數(shù)據(jù)精度的門限。實際的測試方法是將雷達(dá)饋源處使用金屬擋板擋住,在發(fā)射端耦合器(耦合度53 dB)處用信號源送入信號,在線性區(qū)域以5 dB、在高端和低端以1 d B的步長改變信號源的輸出功率,由數(shù)據(jù)處理終端讀取ZH、ZDR、ΦDP、ρHV等數(shù)據(jù)并記錄。以各種偏振參量數(shù)據(jù)為縱軸,以輸入功率為橫軸,畫出曲線圖,從圖中即可看出偏振參量與輸入功率(對應(yīng)與信噪比)的關(guān)系。
圖6是一次隨機測試給出的結(jié)果圖。從圖中可以看出,對于ZDR和ΦDP在輸入信號為-20 dB(對應(yīng)的信噪比為25 dB)以下數(shù)據(jù)的抖動誤差較大,參數(shù)不可信。實際應(yīng)用中為保證偏振數(shù)據(jù)的精度,可以此作為判別偏振數(shù)據(jù)精度的門限。
圖6 偏振數(shù)據(jù)與回波功率大小的關(guān)系
理論上對于固定強度的信號輸入而言,雷達(dá)在不同方位和不同仰角上測量到的值應(yīng)該也是固定的,但是在實際的觀測過程中,由于各種誤差的存在,各參數(shù)在方位和仰角上會有變化,這種變化如果超出一定的范圍,就需要根據(jù)實際情況對雷達(dá)作出調(diào)整。
實際的測試方法是將饋源喇叭用金屬板堵住喇叭口,使其發(fā)射信號產(chǎn)生全反射,信號源信號從發(fā)射機功分器前面耦合入波導(dǎo)??刂评走_(dá)在最低仰角上作PPI掃描,然后以0.2°的步長在垂直方向控制雷達(dá)天線在不同的仰角上作PPI掃描,在終端上保存其輸出的PPI數(shù)據(jù)文件,這時得到的ZH、ZDR、ΦDP、ρHV即為該雷達(dá)的各偏振參數(shù)對這方位和仰角的系統(tǒng)誤差。測試結(jié)果如表1所示,圖7給出了ZDR、ΦDP隨方位起伏變化圖。
表1 CPDRW各偏振參數(shù)對方位、仰角的系統(tǒng)誤差
圖7 偏振參數(shù)隨方位起伏變化圖
雷達(dá)在實際工作時,水平、垂直兩個通道(包括饋線通道和接收通道)會存在一定的通道差,該差值可能會隨著工作環(huán)境的變化而變化,如果其中一個通道發(fā)生較大變化時,仍然用固定的通道差值代入雷達(dá)參數(shù)進行相關(guān)計算,可能會造成有些雙線偏振參數(shù)發(fā)生變化,影響探測的準(zhǔn)確性。因此需要對雙接收通道的幅相一致性進行實時標(biāo)校。標(biāo)校的方法就是在水平、垂直通道上通過信號源輸入兩個相位一致的信號,通過移相器改變兩個通道的信號的相位差,然后從信號處理器上輸出實際測量的相位差。
在雷達(dá)天線俯仰轉(zhuǎn)臺上設(shè)置一個穩(wěn)定的幅相標(biāo)定信號源,該信號源同時輸出兩路等幅同相的測試信號至水平和垂直通道上的定向耦合器測試口。雷達(dá)正常工作時,關(guān)閉該測試信號源電源,以防止泄漏;當(dāng)雷達(dá)需要標(biāo)定時,通過控制電路控制信號源輸出兩路功率、相位相同的測試信號,通過穩(wěn)幅穩(wěn)相電纜分別接入水平、垂直兩個通道的耦合器,經(jīng)過饋線通道后進入雙通道接收機,中頻處理完成后送入信號處理器,信號處理器將強度和相位的檢測結(jié)果送至雷達(dá)終端顯示,并由終端將該結(jié)果送至信號處理器,由其完成一致性實時標(biāo)校。
標(biāo)校信號源在數(shù)字處理分機送來的同步信號控制下,產(chǎn)生用于系統(tǒng)幅相一致性標(biāo)校的兩路標(biāo)定信號,該標(biāo)定信號經(jīng)穩(wěn)幅穩(wěn)相電纜送至波導(dǎo)耦合器,饋入雷達(dá)系統(tǒng)H、V饋線通道,經(jīng)饋線傳輸至對應(yīng)的接收通道,輸出的中頻信號送至數(shù)字處理分機處理,完成幅相一致性標(biāo)校。這種方法的局限性在于:1)不能標(biāo)定發(fā)射通道引起的系統(tǒng)差;2)受標(biāo)校信號源、耦合器和連接電纜性能影響較大。現(xiàn)在已經(jīng)開始考慮直接在發(fā)射端加標(biāo)校源的方式,這樣可以一次性地標(biāo)校出整個硬件通道(包括發(fā)射通道)的系統(tǒng)差,但同樣也存在著器件損耗的影響。
對ZDR的定標(biāo)測試是先分別測量水平和垂直偏振波的反射率因子,然后計算ZDR。具體的測試方法是:從天線后端水平、垂直兩個通道上分別饋入相同強度的輸入信號(逐漸線性遞增),考慮到兩個通道的衰減差,得到理想的ZDR。然后從終端上讀出ZDR的測量值,分析ZDR的測量誤差及其信噪比對ZDR測量的影響。測試結(jié)果表明,當(dāng)對應(yīng)的信噪比 >30 dB時,雙通道的ZDR系統(tǒng)誤差在0.8 dB左右,這部分差值需要在后續(xù)資料處理中給予消除。
為了檢驗利用機內(nèi)和機外儀表標(biāo)定的結(jié)果,采用天線垂直指向的探測方法來進行驗證。理論上,當(dāng)被測量的氣象目標(biāo)滿足方位對稱性時,則其固有ZDR為0 dB,因此在天線波束垂直入射時的測量值就能確定全系統(tǒng)的“修正”因子。另一種方法是假設(shè)風(fēng)暴頂層的冰粒子基本滿足偏振面各向同性條件,與天線仰角無關(guān),從而固有ZDR≈0 dB。由于這樣的雷達(dá)數(shù)據(jù)更容易收集,此方法比較容易實施。這種情況下多時刻、不同距離測量的平均ZDR≈0 dB,KDP≈0°/km。
因為太陽輻射是隨機偏振的,也可以利用太陽信號來模擬一個輸入到兩個通道輻射量相等的外部參考信號源的方式來直接標(biāo)定兩部接收機之間的增益差和相位。當(dāng)雷達(dá)天線對準(zhǔn)太陽進行掃描時,就相當(dāng)于對兩個通道提供相等“噪聲”的功率激勵,這時探測到的回波數(shù)據(jù)也可以作為系統(tǒng)誤差的一個參考。實際的處理方法是通過在不同時次上獲取的幾百個太陽信號數(shù)據(jù)進行統(tǒng)計處理,計算其平均值和標(biāo)準(zhǔn)差,然后將其結(jié)果作為系統(tǒng)誤差用來修正雷達(dá)實際獲取的氣象回波數(shù)據(jù)。由于太陽噪聲是隨機的,理論上其ρHV應(yīng)該為0附近,實際對其觀測獲得的ρHV在0.1附近擺動。
需要說明的是,在輸入的信噪比較小時,雷達(dá)探測到偏振參數(shù)的值擾動較大(如圖8所示),偏振數(shù)據(jù)不可信。因此,當(dāng)偏振數(shù)據(jù)精度的門限較高,大于太陽信號自身的信噪比(一般不超過20 dB)時,此時探測到的偏振數(shù)據(jù)由于抖動較大,對其進行統(tǒng)計處理得到的結(jié)果就沒有明顯的意義,用其結(jié)果來訂正實際氣象回波數(shù)據(jù)得到的結(jié)果就是不正確的。
圖8 Z DR太陽目標(biāo)定標(biāo)過程中獲取的雷達(dá)采樣數(shù)據(jù)(單位:dB)
本文以安徽四創(chuàng)電子股份有限公司和北京市氣象局合作研發(fā)的3836雙線偏振多普勒天氣雷達(dá)升級后的狀態(tài)為主要平臺,根據(jù)雙線偏振雷達(dá)探測能力的要求,在新一代天氣雷達(dá)測試方法的基礎(chǔ)上,提出了雙線偏振雷達(dá)的測試方法,進行了實際的雷達(dá)測試分析。同時利用氣象回波和太陽回波資料,驗證了該雷達(dá)的測試結(jié)果。測試結(jié)果和初步資料分析結(jié)果表明:雙線偏振雷達(dá)的測試方法合理可行,該雙線偏振雷達(dá)探測精度能基本滿足雙線偏振雷達(dá)探測降水和識別粒子系統(tǒng)的需求,資料可信度較高。
[1]DOVIAK R J,BRINGI V,RYZHKOV A,et al.Consideration for Polarimetric Upgrades to Operational WSR-88D Radars[J].Journal of Atmospheric and Oceanic Technology,2000,17(3):257-278.
[2]LIU L,HU Z,FANG W,et al.Calibration and Data Quality Analysis with Mobile C-Band Polarimetric Radar[J].Journal of Meteorological Research,2010,24(4):501-509.
[3]BRANDES E A,ZHANG G,VIVEKANANDAN J.Drop Size Distribution Retrieval with Polarimetric Radar:Model and Application[J].Journal of Applied Meteorology,2004,43(3):461-475.
[4]GORGUCCI E,CHANDRASEKAR V,BRINGI V N,et al.Estimation of Raindrop Size Distribution Parameters from Polarimetric Radar Measurements[J].Journal of the Atmospheric Sciences,2002,59(15):2373-2384.
[5]MELNIKOV V M,ZRINIC D S,DOVIAK R J,et al.Calibration and Performance Analysis of NSSL’s Polarimetric WSR-88D[R].Norman,OK:National Severe Storms Laboratory,2003.
[6]RYZHKOV A V,ZRNIC D S.Assessment of Rainfall Measurement that Uses Specific Differential Phase[J].Journal of Applied Meteorology,1996,35(11):2080-2090.
[7]殷加鵬,李永禎,曹學(xué)軍,等.同時極化雷達(dá)波形隔離度分析與選擇[J].雷達(dá)科學(xué)與技術(shù),2013,11(4):379-384.YIN Jia-peng,LI Yong-zhen,CAO Xue-jun,et al.Analysis and Design of Waveform Isolation for Simultaneous Polarimetric Radar[J].Radar Science and Technology,2013,11(4):379-384.(in Chinese)
[8]劉黎平,錢永甫,王致君.用雙線偏振雷達(dá)研究云內(nèi)粒子相態(tài)及尺度的空間分布[J].氣象學(xué)報,1996,54(5):590-599.
[9]王致君,楚榮忠.X波段雙通道同時收發(fā)式多普勒偏振天氣雷達(dá)[J].高原氣象,2007,26(1):135-140.
[10]曹俊武,劉黎平,葛潤生.模糊邏輯法在雙線偏振雷達(dá)識別降水粒子相態(tài)中的研究[J].大氣科學(xué),2005,29(5):827-836.
[11]曹俊武,胡志群,陳曉,等.影響雙線偏振雷達(dá)相位探測精度的分析[J].高原氣象,2011,30(3):817-823.