范根昌 樊嶸 王梧熠

(河北省人工影響天氣辦公室， 石家莊 050021)

引 言

雙偏振多普勒天氣雷達(dá)能在水平和垂直方向上發(fā)出偏振波，不僅可以得到回波強(qiáng)度(ZH)、徑向速度(V)、譜寬(W)等常用雷達(dá)參量，還可以通過比較其雙通道回波功率大小以及相移等獲取差分反射率因子(ZDR)、差分傳播相移(ΦDP)、差分傳播相移率(KDP)和零滯后相關(guān)系數(shù)(ρHV)等多個(gè)雷達(dá)偏振參量[1-3]，常用來研究降水粒子結(jié)構(gòu)和微物理機(jī)制。與普通雷達(dá)相比，該雷達(dá)在識(shí)別粒子相態(tài)等方面擁有較大優(yōu)勢(shì)。但在雷達(dá)運(yùn)行過程中，數(shù)據(jù)質(zhì)量會(huì)因地物雜波、雷達(dá)內(nèi)部及外部噪聲、電磁波衰減等影響，在使用雙偏振多普勒雷達(dá)數(shù)據(jù)之前，評(píng)估資料質(zhì)量便顯得尤為重要，為數(shù)據(jù)有效應(yīng)用奠定基礎(chǔ)。

雙偏振多普勒雷達(dá)資料質(zhì)量評(píng)估和控制是保證數(shù)據(jù)可信和產(chǎn)品可靠性的前提，因此長(zhǎng)期受到國(guó)內(nèi)外學(xué)者的關(guān)注。ZDR作為關(guān)鍵的偏振參量，表示水平偏振和垂直偏振的反射率因子之比。Ryzhkov, et al[4]發(fā)現(xiàn)，當(dāng)ZDR的誤差為0.2 dB時(shí)，將會(huì)產(chǎn)生18%的降水估測(cè)誤差，所以在應(yīng)用ZDR前，需要對(duì)其訂正系統(tǒng)誤差。在對(duì)ZDR的系統(tǒng)誤差訂正方面，國(guó)外許多學(xué)者使用了垂直天頂掃描、太陽(yáng)標(biāo)定法和氣象回波法等訂正方法[5-6]。國(guó)內(nèi)吳林林等[7]將自然目標(biāo)物法應(yīng)用在訂正ZDR系統(tǒng)誤差之中，結(jié)果表明使用該方法訂正后的ZDR比訂正之前的數(shù)據(jù)質(zhì)量有了明顯提高。胡志群等[8]在將許多ZDR標(biāo)定法進(jìn)行比較后，發(fā)現(xiàn)微雨滴法效果最好。杜牧云等[1]將垂直掃描數(shù)據(jù)應(yīng)用于對(duì)C波段雙偏振雷達(dá)的ZDR系統(tǒng)誤差訂正中，對(duì)ZDR系統(tǒng)誤差進(jìn)行了有效訂正。在雷達(dá)探測(cè)過程中，ΦDP會(huì)出現(xiàn)抖動(dòng)起伏的情況。很多國(guó)外學(xué)者采用低通濾波器將高頻噪聲過濾掉，從而保留曲線平均趨勢(shì)，這方面的國(guó)內(nèi)研究較少，主要內(nèi)容是關(guān)于數(shù)據(jù)預(yù)處理部分[9-10]。胡明寶等[11]通過比較雷達(dá)衰減訂正的很多方法后，得出了KDP法效果最佳的結(jié)論。

河北于2016年在邢臺(tái)布設(shè)了1部724XSP雙偏振雷達(dá)。為了解和評(píng)估數(shù)據(jù)的質(zhì)量狀況，首先選取偏振雷達(dá)上空層狀云弱降水過程，對(duì)雙偏振雷達(dá)的ZDR、ρHV、ΦDP等偏振參量進(jìn)行質(zhì)量以及穩(wěn)定性分析，并對(duì)ZH和ZDR進(jìn)行簡(jiǎn)單衰減訂正并和系統(tǒng)誤差訂正，并將訂正結(jié)果與位于石家莊新樂站(38.351 9°N，114.711 9°E，海拔134.8 m)的CINRAD/SA多普勒天氣雷達(dá)(簡(jiǎn)稱SA天氣雷達(dá))相對(duì)比，希望能夠深入了解724XSP雙偏振雷達(dá)對(duì)云、降水的觀測(cè)機(jī)制以及探測(cè)能力。

1 724XSP雙偏振雷達(dá)參數(shù)

本文所使用的724XSP雙偏振雷達(dá)是一部車載X波段全固態(tài)雙線偏振多普勒天氣雷達(dá)。該雷達(dá)架設(shè)在河北邢臺(tái)市內(nèi)丘縣柳林氣象站內(nèi)(37.290 7°N，114.386 9°E，海拔67.9 m)，已經(jīng)進(jìn)行多次降水個(gè)例觀測(cè)。該雷達(dá)采用新一代固態(tài)雷達(dá)發(fā)射機(jī)，主要技術(shù)參數(shù)見表1。X波段全固態(tài)雙線偏振多普勒天氣雷達(dá)長(zhǎng)短脈沖一致性將影響到數(shù)據(jù)質(zhì)量和后續(xù)應(yīng)用，采用高仰角數(shù)據(jù)交界處以外的長(zhǎng)短脈沖距離庫(kù)進(jìn)行比較，求二者隨方位的均值來進(jìn)行長(zhǎng)短脈沖訂正，并對(duì)交界處的距離庫(kù)采用5點(diǎn)滑動(dòng)平均進(jìn)行平滑處理。該雷達(dá)終端軟件已經(jīng)對(duì)長(zhǎng)短脈沖不連續(xù)進(jìn)行了訂正，基數(shù)據(jù)輸出為訂正后產(chǎn)品。

表1 724XSP雙偏振雷達(dá)相關(guān)性能參數(shù)Table 1 Performance parameters of 724XSP dual-linear polarization radar

2 偏振數(shù)據(jù)質(zhì)量分析

2.1 系統(tǒng)初始相位分析

圖1 (a) 不同時(shí)刻初始相位隨方位角變化；(b) 不同仰角初始相位隨方位角變化；(c) 2017年7月28日08時(shí)09分(北京時(shí)，下同)PPI圖Fig.1 (a) Initial phase changes with azimuth angle at different time;(b) initial phase changes with azimuth angle at different elevation angles;(c) the PPI diagram at 08∶09 BST on 28 July 2017

圖2 2017年7月28日08時(shí)56分 90°仰角層狀云降水資料中小雨天氣SNR-ZDR(a,單位：dB)和SNR-CC(b,單位：dB)的散點(diǎn)分布Fig.2 Scatter distribution of SNR-ZDR (a, unit:dB) and SNR-CC (b, unit:dB) in light rain in stratified cloud precipitation data of90° elevation angle at 08∶56 BST on July 28， 2017

724XSP雙偏振雷達(dá)為雙發(fā)雙收模式，在實(shí)際探測(cè)過程中，兩個(gè)通道的不完全匹配將會(huì)引起ΦDP產(chǎn)生偏差。首先分析初始相位的穩(wěn)定性和偏差。本文參考杜牧云等[1]提出的沿徑向向外搜索連續(xù)10個(gè)庫(kù)CC大于0.9、且ΦDP的標(biāo)準(zhǔn)差小于3°，那么這些庫(kù)的平均值便被定義成初始相位。以2017年7月28日體掃數(shù)據(jù)為例，獲得每根徑向上的ΦDP初始值，圖1a、b分別為不同時(shí)刻、不同仰角求取的0～360°方位角初始相位分布，其中圖1a是仰角為2.4°得到的相位變化。圖1a的ΦDP方差為0.302，圖1b的ΦDP方差為0.285，圖1c 是2017年7月28日08時(shí)09分(北京時(shí)，下同)時(shí)刻對(duì)應(yīng)的PPI圖。從圖1中可以看出不同時(shí)刻、不同仰角的初始相位在徑向上波動(dòng)不大且變化趨勢(shì)基本一致，系統(tǒng)初始相位隨時(shí)間的變化率較小。將2017年6月21—22日、7月6—7日類似的層狀云降水過程分析后，也可得出相同結(jié)果。

2.2 信噪比對(duì)ZDR和CC影響分析

在信躁比(Signal Noise Ratio, SNR)小于15 dB時(shí)，ZDR離散值較多，且波動(dòng)較大，其對(duì)應(yīng)的相關(guān)系數(shù)CC也偏小很多，且隨SNR降低呈減小趨勢(shì)(圖2)。這表明當(dāng)SNR較低時(shí)所取得的數(shù)據(jù)質(zhì)量不高。當(dāng)信噪比變得越來越大時(shí)，ZDR的分布也越發(fā)穩(wěn)定，CC也大體上超過0.9，而且快速接近1。在SNR大于30 dB時(shí)，也可以看到CC有一段異常數(shù)據(jù)區(qū)域(圖2中紅框區(qū))，下面討論不同高度下的SNR-ZDR和SNR-CC分布。

圖3為3個(gè)不同高度區(qū)間下SNR-ZDR和SNR-CC的散點(diǎn)分布。從圖中可以看出，低層(H<3.6 km)時(shí)，數(shù)據(jù)信噪比較高，且ZDR數(shù)據(jù)質(zhì)量較好，但圖2紅框數(shù)據(jù)異常區(qū)域出現(xiàn)在3.6 km以下，這是由于層狀云降水主要的微物理過程是冰晶聚合成為干雪，到融化層后變?yōu)闈裱?，濕雪進(jìn)一步融化成為雨引起。不過圖3中大部分?jǐn)?shù)據(jù)CC均大于0.8；當(dāng)高度偏高(H>5 km)時(shí)，因?yàn)楹屠走_(dá)相距較遠(yuǎn)，SNR整體偏低(<15 dB)，CC也明顯偏下，大部分都小于0.9，這是由于噪聲對(duì)數(shù)據(jù)的干擾作用。此外，還可能是冰晶在高空風(fēng)取向的作用下，導(dǎo)致其離散度變大[12]；處于中間層(3.60.9。其他個(gè)例也得到相同結(jié)論。

圖3 圖2中不同高度區(qū)間的SNR-ZDR和SNR-CC散點(diǎn)分布：(a、b) 3.6 km以下；(c、d) 3.6～5 km；(e、f)5 km以上Fig.3 Scatter distribution of SNR-ZDR and SNR-CC in different height intervals in fig.2： (a，b) less than 3.6 km；(c，d)3.6-5 km；(e，f) more than 5 km

2.3 ZDR系統(tǒng)誤差穩(wěn)定度分析

利用雷達(dá)垂直掃描數(shù)據(jù)對(duì)探測(cè)到的小雨資料ZDR進(jìn)行系統(tǒng)誤差分析，為避免信噪比、融化層等的影響，選擇的數(shù)據(jù)高度在1 000 m附近，結(jié)合前文分析為保證數(shù)據(jù)的穩(wěn)定性，樣本庫(kù)選擇SNR>20 dB、ZH<28 dBZ和CC>0.97。

圖4為小雨天氣下SNR-ZDR散點(diǎn)分布。信噪比SNR以0.5 dB為間隔求取每個(gè)信噪比區(qū)間內(nèi)對(duì)應(yīng)的ZDR均值。該圖說明在信噪比變大的過程中，小雨的ZDR系統(tǒng)誤差訂正曲線變化幅度不大(其中>32 dB由于樣本數(shù)較小波動(dòng)較大)，表明在雷達(dá)運(yùn)行時(shí)系統(tǒng)較為穩(wěn)定，當(dāng)信噪比波動(dòng)時(shí)，ZDR系統(tǒng)誤差波動(dòng)不明顯。表2為2017年7月28日的降水過程資料得到的不同時(shí)刻ZDR系統(tǒng)誤差均值，可以看出ZDR系統(tǒng)誤差均值大都處于-0.2～-0.1 dB，且標(biāo)準(zhǔn)差穩(wěn)定在0.4～0.5 dB。根據(jù)一次小雨過程共18 970個(gè)樣本數(shù)統(tǒng)計(jì)，ZDR系統(tǒng)誤差均值為-0.18 dB，標(biāo)準(zhǔn)差為0.44 dB。

圖4 2017年7月28日08時(shí)56分90°仰角層狀云降水資料中小雨天氣零度層以下SNR-ZDR散點(diǎn)分布(其中紅色點(diǎn)線為均值，藍(lán)線為ZDR系統(tǒng)誤差訂正曲線)Fig.4 Scattered distribution of SNR-ZDR below zero layerin light rain in stratified cloud precipitation data of90° elevation angle at 08∶56 BST on July 28， 2017

表2 2017年7月28日不同時(shí)刻觀測(cè)個(gè)例中ZDR系統(tǒng)誤差、標(biāo)準(zhǔn)差統(tǒng)計(jì)Table 2 Statistics of ZDR system errors and standard deviationsin individual cases observed at different times on July 28， 2017

3 衰減訂正及結(jié)果評(píng)估

雨滴會(huì)吸收穿過降水粒子的電磁波，進(jìn)而引起電磁波衰減情況的出現(xiàn)。對(duì)于波長(zhǎng)越短的電磁波，被衰減的程度越大。由于ΦDP是相位變化，因此衰減對(duì)ΦDP的影響較小。所以采用ΦDP來對(duì)ZH及ZDR進(jìn)行衰減訂正效果較好。

圖5 2017年7月16日08時(shí)58分：(a) 訂正前724XSP雙偏振雷達(dá)反射率因子、(b) 訂正后724XSP雙偏振雷達(dá)反射率因子和(c) SA天氣雷達(dá)3 000 m高度反射率CAPPI Fig.5 (a) 724XSP dual polarization radar reflectivity factor before revision;(b) 724XSP dual polarization radar reflectivity factor after revision;(c) 3 000 m altitude reflectance CAPPI chart of SA weather radar at 08∶58 BST on July 16， 2017

本文用ΦDP訂正ZH及ZDR的公式為：

(1)

(2)

為了客觀評(píng)估訂正前后724XSP雙偏振雷達(dá)回波反射率因子資料，需要把724XSP雙偏振雷達(dá)和位于新樂SA天氣雷達(dá)二者極坐標(biāo)下空間分辨率不等的回波數(shù)據(jù)統(tǒng)一插值到空間分布均勻的笛卡爾坐標(biāo)系下。本文參考張志強(qiáng)等[14]提出的格點(diǎn)對(duì)比法進(jìn)行對(duì)比，格點(diǎn)對(duì)比法是把不同雷達(dá)的數(shù)據(jù)插值到同一高度層，然后比較它們的相交區(qū)域，該方法可以較好的找到其對(duì)應(yīng)點(diǎn)，但缺點(diǎn)在于不同雷達(dá)的庫(kù)長(zhǎng)和波束寬度存在差異，或許會(huì)出現(xiàn)插值誤差。在2017年7月16日對(duì)流天氣過程中，選取(36.9°～37.7°N，114.1°～114.8°E)區(qū)域，作出該區(qū)域3 000 m高度處兩部雷達(dá)CAPPI圖像，并用雙線性差值法[14]插值成120×100網(wǎng)格，格點(diǎn)分辨率約為3 km。3 000 m高度CAPPI數(shù)據(jù)避開了地物及遮擋的影響，且回波較為連續(xù)，圖5a、b空白圓圈為雷達(dá)站所在位置。

由圖5可見，724XSP雙偏振雷達(dá)與SA天氣雷達(dá)回波強(qiáng)度分布及結(jié)構(gòu)特征大致相似。將系統(tǒng)誤差導(dǎo)致的明顯不合理點(diǎn)剔除后得到的訂正后的反射率因子相關(guān)系數(shù)達(dá)到0.78，與訂正前相比，訂正后的擬合結(jié)果得到明顯提高，如圖6所示。從訂正后的擬合情況來看，在回波較弱時(shí)724XSP雙偏振雷達(dá)存在低估而回波較強(qiáng)時(shí)存在高估現(xiàn)象，但整體偏差不大，后續(xù)工作將結(jié)合雨滴譜對(duì)選取更優(yōu)化的α值進(jìn)行研究。訂正前724XSP雙偏振雷達(dá)探測(cè)到的強(qiáng)回波(>40 dBZ)面積要小于SA天氣雷達(dá)，雷達(dá)所缺失的回波主要集中在遠(yuǎn)離724XSP雙偏振雷達(dá)一側(cè)，主要原因是724XSP雷達(dá)波束穿過回波強(qiáng)中心后被衰減，影響其探測(cè)強(qiáng)回波后的區(qū)域。訂正后，ZH得到改善，與SA天氣雷達(dá)對(duì)比強(qiáng)回波中心位置及強(qiáng)度基本相同，而且在強(qiáng)回波區(qū)域后部的回波區(qū)衰減也得到了了明顯改善。

通過ZH和ZDR兩者間的理論表達(dá)式來更好的說明訂正前和訂正后的區(qū)別以及效果，表達(dá)式為Bringi[18]使用散射模式計(jì)算出的理論結(jié)果。圖7為訂正前、后ZDR-ZH的散點(diǎn)分布。從圖中可以看出經(jīng)過訂正后，其離散度減小了很多。當(dāng)ZH<25 dBZ時(shí)，其ZDR值均較小，集中在0值波動(dòng)。當(dāng)ZH值處于25～45 dBZ時(shí)，訂正后的ZDR值大都處于0.2～2 dB，并且當(dāng)ZH值增大時(shí)，訂正后的ZDR具有上升趨勢(shì)，這種變化同Ryzhkov， et al[3]的結(jié)論一致。通過對(duì)該雷達(dá)探測(cè)到的其它過程數(shù)據(jù)的分析，也可以得出同樣的結(jié)論。

4 結(jié)論

(1)724XSP雙偏振雷達(dá)初始相位ΦDP(0)隨仰角、方位的變化都是較為穩(wěn)定的，波動(dòng)較小。分析垂直掃描數(shù)據(jù)，發(fā)現(xiàn)ZDR系統(tǒng)偏差較為穩(wěn)定，變化不大。說明724XSP雙偏振雷達(dá)系統(tǒng)較為穩(wěn)定，分析結(jié)論可以作為一段時(shí)間內(nèi)的系統(tǒng)誤差進(jìn)行訂正[14-16]。

(2)小雨天氣，ZDR、CC隨信噪比SNR變化穩(wěn)定，數(shù)據(jù)可信度高，隨信噪比SNR增大數(shù)據(jù)越穩(wěn)定，尤其當(dāng)SNR>10 dB時(shí)。

(3)724XSP雙偏振雷達(dá)與SA天氣雷達(dá)對(duì)同一回波區(qū)域，觀測(cè)結(jié)果大體相似，經(jīng)過衰減訂正后724XSP雙偏振雷達(dá)反射率因子ZH衰減訂正結(jié)果與S波段雷達(dá)更為接近，尤其強(qiáng)回波區(qū)改善明顯。經(jīng)過衰減訂正和系統(tǒng)誤差訂正后ZDR-ZH散點(diǎn)分布更接近理論值，說明訂正效果好。

圖6 SA天氣雷達(dá)和訂正前雙偏振雷達(dá)反射率因子散點(diǎn)對(duì)比：(a)訂正前；(b)訂正后Fig.6 The scatter plot of reflectivity factor between SA weather radar and unrevised dual polarization radar：(a) before revision;(b) after revision

圖7 2017年8月12日18時(shí)37分4.3°仰角下ZDR-ZH散點(diǎn)分布：(a) 訂正前；(b) 訂正后Fig.7 Scatter distribution of ZDR-ZH at 4.3° elevation angle at 18∶37 BST on 12 August， 2017 ： (a) before revision;(b) after revision