雷衛(wèi)延，路永平，王明輝，周嘉健，甘泉

（1.惠州市氣象局，廣東惠州 516003；2.廣東省氣象探測數(shù)據(jù)中心，廣東廣州 510080）

目前粵港澳大灣區(qū)已建成由25部X波段雙極化相控陣天氣雷達(dá)組成的災(zāi)害性天氣精細(xì)化監(jiān)測網(wǎng)，時空分辨率分別達(dá)到1 min和30 m，可有效彌補(bǔ)S波段天氣雷達(dá)對局地中小尺度強(qiáng)對流天氣監(jiān)測的不足，在防御暴雨、雷雨大風(fēng)、冰雹等發(fā)揮了積極作用，為防災(zāi)減災(zāi)爭取了更多的時間。但目前對X波段雙極化相控陣?yán)走_(dá)尚未健全定標(biāo)技術(shù)規(guī)范，偏振量數(shù)據(jù)質(zhì)量難以保障，因此迫切需要開展相控陣?yán)走_(dá)定標(biāo)試驗(yàn)，研究能滿足業(yè)務(wù)要求的定標(biāo)技術(shù)，盡早投入業(yè)務(wù)使用，發(fā)揮雷達(dá)效益最大化。國外研究人員Sachidananda等［1］于1985年提出Zdr測量精度不應(yīng)低于0.2 dB，理想精度最好是0.1 dB；Ryzhkov等［2］于2005年的研究結(jié)果表明，0.2 dB的Zdr偏差將導(dǎo)致雷達(dá)利用Zdr進(jìn)行降水估測產(chǎn)生18%的誤差，高精度探測是成功實(shí)現(xiàn)雙極化雷達(dá)氣象應(yīng)用的前提；李斐斐等［3］、胡東明等［4］研究分析了雙線偏振天氣雷達(dá)雙通道一致性測試及標(biāo)校方法；劉黎平等［5］、汪旭東等［6］、程張凡等［7］對X波段相控陣天氣雷達(dá)測試定標(biāo)方法開展了一定研究；李兆明等［8］對使用金屬球定標(biāo)X波段固態(tài)天氣雷達(dá)開展研究，取得初步結(jié)論。

目前國內(nèi)外雙極化雷達(dá)定標(biāo)方法主要有太陽法、金屬球法、小雨法、信號源法等，業(yè)務(wù)中使用的S波段雙極化雷達(dá)具有完善的自動定標(biāo)系統(tǒng)和技術(shù)規(guī)范［9-10］，由于X波段雙極化相控陣天氣雷達(dá)近幾年才投入使用，相對而言國內(nèi)外對其定標(biāo)研究文獻(xiàn)較少，本研究利用無人機(jī)懸掛金屬球開展雙通道強(qiáng)度一致性定標(biāo)研究，試驗(yàn)結(jié)果表明：使用無人機(jī)懸掛金屬球定標(biāo)方法切實(shí)可行，能夠滿足技術(shù)規(guī)范要求。

1 金屬球定標(biāo)原理

雷達(dá)定標(biāo)用的金屬球（圖1）對電磁波的散射是向空間各個方向傳播，后向散射截面與波長的關(guān)系是已知的，其半徑大小的選取與雷達(dá)波長有關(guān)，使后向散射位于幾何光學(xué)散射區(qū)域，即后向散射截面僅取決于金屬球半徑［6］。

圖1 雷達(dá)定標(biāo)用金屬球

根據(jù)天氣雷達(dá)方程，對于已知雷達(dá)目標(biāo)，回波強(qiáng)度可采用公式（1）計算。其中，Pr為回波信號功率（W）；Pt為雷達(dá)發(fā)射脈沖功率（W）；R為回波距離（m）；λ為雷達(dá)工作波長（m）；G為天線增益；Gr為接收通道增益；σ為目標(biāo)的雷達(dá)反射截面積RCS（m2）。金屬球作為探測目標(biāo)引起的差分反射率Zdr理論上為零值，因此使用金屬球定標(biāo)雷達(dá)實(shí)測Zdr值，即為雷達(dá)差分反射率系統(tǒng)誤差，檢驗(yàn)系統(tǒng)誤差是否滿足＜0.2 dB的指標(biāo)要求。如果超出指標(biāo)范圍，利用參數(shù)訂正使雷達(dá)Zdr小于0.2 dB，實(shí)現(xiàn)雷達(dá)從發(fā)射、接收、數(shù)據(jù)處理等全鏈路的系統(tǒng)定標(biāo)，確保雷達(dá)觀測數(shù)據(jù)質(zhì)量。

2 無人機(jī)飛行軌跡設(shè)計

相控陣?yán)走_(dá)的主瓣方向即3 dB波束寬度在垂直剖面上對應(yīng)一個橢圓區(qū)域，借助無人機(jī)上的定位設(shè)備，通過控制無人機(jī)使所攜帶的金屬球大致處于上述橢圓區(qū)域，由于雷達(dá)的目標(biāo)定位和無人機(jī)的位置定位有一定的誤差［11］，需要移動金屬球掃描一段區(qū)域，雷達(dá)觀測到金屬球反射信號的最大值時金屬球所在的位置［12］，即剛好處于雷達(dá)主瓣方向，并以該位置作為基準(zhǔn)點(diǎn)設(shè)計無人機(jī)飛行軌跡，如圖2所示。

圖2 無人機(jī)金屬球定標(biāo)示意圖

雷達(dá)定標(biāo)使用寬度為1μs的雷達(dá)短脈沖，考慮到雷達(dá)觀測盲區(qū)約為200 m，定標(biāo)時應(yīng)避開雷達(dá)盲區(qū)，標(biāo)定時金屬球距雷達(dá)天線選擇300 m左右為宜。當(dāng)雷達(dá)發(fā)射波束衰減25 dB時，主瓣寬度為±2°，第一旁瓣寬度為±2.6°，利用三角關(guān)系計算得到第一旁瓣距離雷達(dá)法線約13.6 m。為減少無人機(jī)對定標(biāo)的影響，利用20 m的系留細(xì)繩懸掛金屬球，使無人機(jī)遠(yuǎn)離雷達(dá)第一旁瓣，此時無人機(jī)的影響可以忽略不計。無人機(jī)帶動金屬球以基準(zhǔn)點(diǎn)為中心在12 m×12 m正方形范圍內(nèi)呈方波運(yùn)動，以5個等間距運(yùn)動為一個循環(huán)，最后一次上升到頂部后不做水平移動，原路下降重新回到最左面第一個起點(diǎn)位置，以此類推每個運(yùn)動周期約90 s，金屬球運(yùn)動軌跡示意圖如圖3所示。

圖3 金屬球運(yùn)動軌跡

3 定標(biāo)實(shí)施方案

選擇無風(fēng)或小風(fēng)的天氣開展金屬球定標(biāo)，以減少金屬球劇烈運(yùn)動對定標(biāo)的影響。按照雷達(dá)正常業(yè)務(wù)觀測模式，按照第2章無人機(jī)飛行軌跡策略，每個體掃周期90 s雷達(dá)共采樣7 200個點(diǎn)，每個采樣點(diǎn)包含32個脈沖累計。

雷達(dá)水平和垂直收發(fā)通道的阻抗均為固定值，因此回波功率可以通過電壓信號來反算，二者一一對應(yīng)，成正比關(guān)系。如圖4所示，縱坐標(biāo)代表電壓信號，橫坐標(biāo)代表采用點(diǎn)數(shù)，當(dāng)金屬球做垂直運(yùn)動時，上升和下降會穿過雷達(dá)主波束中心，電壓信號先增加后減小，產(chǎn)生①、③、⑤、⑦峰值功率極大值；圖形中兩個方波的寬度代表金屬球在兩個垂直運(yùn)動之間的平移，由于遠(yuǎn)離主波束中心此時產(chǎn)生②、④、⑥功率極小值；最后一個上升和下降是原路徑，無水平移動，所以兩個峰值⑨、○1間距較近。

圖4 實(shí)測電壓信號最大值曲線

當(dāng)金屬球經(jīng)過雷達(dá)主波束中心時，電壓大于特定閾值0.05 V為有效數(shù)據(jù)，分別計算水平和垂直通道所有有效數(shù)據(jù)的平均值，以此計算差分反射率Zdr。針對水平和垂直回波信號功率計算公式（1），回波距離可以通過無人機(jī)的定位系統(tǒng)計算出來，綜合考慮雷達(dá)盲區(qū)和無人機(jī)承載因素，金屬球到雷達(dá)天線陣面中心的距離應(yīng)大于200 m，但又不易太遠(yuǎn)；雷達(dá)的工作波長、天線增益在出廠前已經(jīng)過測試，為已知量；雷達(dá)觀測目標(biāo)金屬球直徑選定為13.8 cm，金屬球的RCS與極化方式無關(guān)，可采用公式（2）計算。其中，r為金屬球的半徑（m）。由此計算得知定標(biāo)用的金屬球的RCS為0.015 m2。

4 定標(biāo)過程

2020年1—10月在廣州多次開展金屬球定標(biāo)試驗(yàn)，試驗(yàn)數(shù)據(jù)如表1所示。

表1 金屬球定標(biāo)試驗(yàn)結(jié)果 dB

定標(biāo)方法：首先不升空金屬球，雷達(dá)依次選定0°、5°、10°、15°仰角對空氣進(jìn)行360°掃描，測試周圍環(huán)境噪聲和地物雜波情況；其次選擇不受外界干擾的方向，按照無人機(jī)預(yù)先設(shè)定飛行策略，金屬球距雷達(dá)天線約300 m的位置，依次對不同仰角開展定標(biāo)。如表1所示，在1—10月開展的多次定標(biāo)過程中花都、帽峰山、南海3部雷達(dá)在仰角為0°方向上存在遮擋，數(shù)據(jù)跳躍比較大，定標(biāo)數(shù)據(jù)不可信；而花都、帽峰山、南海、番禺4部雷達(dá)在仰角為5°、10°、15°方向上無外界干擾，Zdr平均值均表現(xiàn)穩(wěn)定，波動范圍＜0.2 dB，標(biāo)準(zhǔn)差均滿足＜0.1指標(biāo)要求。

根據(jù)表1可知番禺雷達(dá)10°仰角Zdr平均值為-0.308 8 dB，超出＜0.2 dB的指標(biāo)要求。一次雷達(dá)金屬定標(biāo)過程如下，修正前Zdr值偏差為-0.54 dB，分布如圖5a所示；經(jīng)過參數(shù)修正再次進(jìn)行定標(biāo)Zdr值偏差為-0.072 dB，分布如圖5b所示，表明通過金屬球定標(biāo)法能夠使雷達(dá)Zdr指標(biāo)滿足業(yè)務(wù)規(guī)范要求。

圖5 修正前（a）和修正后（b）Zdr概率分布

5 結(jié)論

通過開展X波段雙極化相控陣天氣雷達(dá)金屬球定標(biāo)試驗(yàn)，對天線、發(fā)射通道、接收通道、數(shù)據(jù)處理等全鏈路一致性進(jìn)行測試分析，對不同的仰角進(jìn)行定標(biāo)，對超出指標(biāo)要求的仰角進(jìn)行修正，得到如下結(jié)論：

1）開展金屬球定標(biāo)時，應(yīng)選擇無風(fēng)或小風(fēng)的天氣進(jìn)行，每個仰角均應(yīng)選擇無環(huán)境噪聲和地物雜波干擾的方位作為測試方向。

2）對廣州4部雷達(dá)現(xiàn)場定標(biāo)試驗(yàn)進(jìn)行分析，0°仰角容易受到地物干擾，數(shù)據(jù)波動大，定標(biāo)結(jié)果不符合要求，金屬球定標(biāo)存在一定的弊端。對5°、10°、15°等仰角開展金屬球定標(biāo)，Zdr標(biāo)準(zhǔn)差均小于0.1 dB，波動范圍均小于0.2 dB，表明定標(biāo)結(jié)論可信，超出指標(biāo)的仰角經(jīng)過參數(shù)修正，Zdr平均值均能滿足＜0.2 dB的指標(biāo)要求。

從定標(biāo)結(jié)果來看，絕大多數(shù)雷達(dá)定標(biāo)結(jié)果都能滿足指標(biāo)要求，僅個別雷達(dá)的某個仰角超出誤差范圍，建議每年開展一次金屬球定標(biāo)業(yè)務(wù)。