新一代天氣雷達(dá)標(biāo)準(zhǔn)格式數(shù)據(jù)可視化監(jiān)控系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)方法

謝曉林，胡 迪，羅宇昂，李 翔，景號(hào)然

(1.四川省氣象探測(cè)數(shù)據(jù)中心，成都 610072；2.高原與盆地暴雨旱澇災(zāi)害四川省重點(diǎn)實(shí)驗(yàn)室，成都 610072；3.四川省氣象臺(tái)，成都 610072)

0 引言

根據(jù)中國氣象局要求，2018—2020年全國新一代天氣雷達(dá)部署完善后已經(jīng)基本完成技術(shù)和標(biāo)準(zhǔn)統(tǒng)一，所有業(yè)務(wù)雷達(dá)實(shí)現(xiàn)了數(shù)據(jù)格式標(biāo)準(zhǔn)化[1]。雷達(dá)數(shù)據(jù)的傳輸方式由FTP文件傳輸升級(jí)為數(shù)據(jù)流傳輸，體掃的每1個(gè)仰角層都會(huì)形成1個(gè)單獨(dú)的文件，傳輸時(shí)效大大提高，由6 min/個(gè)文件升級(jí)為1 min/個(gè)單仰角文件和1個(gè)完全體掃文件[2]。

由于技術(shù)體制的升級(jí)，研究一種能夠?qū)崟r(shí)監(jiān)控雷達(dá)標(biāo)準(zhǔn)格式數(shù)據(jù)的可視化系統(tǒng)并應(yīng)用于業(yè)務(wù)，對(duì)提高天氣預(yù)報(bào)時(shí)效和監(jiān)測(cè)雷達(dá)數(shù)據(jù)質(zhì)量等方面具有十分重要的意義。文章闡述了一種針對(duì)新一代天氣雷達(dá)標(biāo)準(zhǔn)格式數(shù)據(jù)可視化監(jiān)控系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)方法。

1 總體設(shè)計(jì)

新一代天氣雷達(dá)標(biāo)準(zhǔn)格式數(shù)據(jù)可視化監(jiān)控系統(tǒng)能對(duì)天氣雷達(dá)標(biāo)準(zhǔn)格式數(shù)據(jù)的傳輸?shù)綀?bào)進(jìn)行監(jiān)控，能讀取文件內(nèi)容，進(jìn)行實(shí)時(shí)可視化展示，可以監(jiān)控雷達(dá)數(shù)據(jù)質(zhì)量，并進(jìn)行統(tǒng)計(jì)分析。系統(tǒng)可對(duì)四川12部新一代天氣雷達(dá)標(biāo)準(zhǔn)格式數(shù)據(jù)進(jìn)行可視化監(jiān)控，同時(shí)也支持對(duì)其他非業(yè)務(wù)天氣雷達(dá)數(shù)據(jù)的可視化，系統(tǒng)分為在線和本地兩種查看模式。

在線模式可自動(dòng)連接數(shù)據(jù)服務(wù)器，查看對(duì)應(yīng)站點(diǎn)的實(shí)時(shí)數(shù)據(jù)，可任意打開并讀取數(shù)據(jù)進(jìn)行可視化操作，如選擇不同數(shù)據(jù)產(chǎn)品、剖面、截圖等，可自動(dòng)可視化回放。在線模式也可實(shí)時(shí)可視化監(jiān)控?cái)?shù)據(jù)質(zhì)量，系統(tǒng)30 s自動(dòng)刷新，讀取最新數(shù)據(jù)并顯示，讀取數(shù)據(jù)后，能自動(dòng)分析數(shù)據(jù)質(zhì)量，如數(shù)據(jù)正常則生成LOG日志文件，異常則生成ALARM告警文件。系統(tǒng)能自動(dòng)統(tǒng)計(jì)一定時(shí)間區(qū)間內(nèi)的雷達(dá)數(shù)據(jù)質(zhì)量情況，并形成報(bào)表。

由于連接數(shù)據(jù)服務(wù)器需要內(nèi)網(wǎng)支持，為了方便用戶查看本地歷史數(shù)據(jù)，系統(tǒng)在線模式運(yùn)行時(shí)，監(jiān)控的所有數(shù)據(jù)會(huì)自動(dòng)保存到本地，并以站號(hào)和時(shí)間分類。打開本地?cái)?shù)據(jù)，可以對(duì)數(shù)據(jù)進(jìn)行常規(guī)分析，如切換產(chǎn)品或剖面等，支持對(duì)本地?cái)?shù)據(jù)的回放，生成數(shù)據(jù)質(zhì)量報(bào)表。

具體功能和流程圖如圖1所示。

2 關(guān)鍵技術(shù)

2.1 標(biāo)準(zhǔn)格式數(shù)據(jù)讀取

根據(jù)數(shù)據(jù)文件名中“_FMT”標(biāo)志可以判斷該文件為標(biāo)準(zhǔn)格式數(shù)據(jù)，如：

Z_RADR_I_Z9280_20210820150349_O_DOR_SC_CAP_FMT.bin.bz2表示成都雷達(dá)站2021-08-20T15：03：49生成的雷達(dá)標(biāo)準(zhǔn)格式完整體掃數(shù)據(jù)，雷達(dá)型號(hào)為CINRAD/SC，通過文件名中“_CUT”標(biāo)識(shí)可以判斷出該文件為標(biāo)準(zhǔn)格式數(shù)據(jù)的單層仰角數(shù)據(jù)，如：

Z_RADR_I_Z9280_20211018024441_O_DOR-CUT_SAD_CAP_43_7_FMT.bin.bz2表示該完整體掃數(shù)據(jù)中的第7個(gè)層仰角數(shù)據(jù)，仰角為4.3°。

圖1 系統(tǒng)功能流程

由于所有標(biāo)準(zhǔn)格式數(shù)據(jù)都被壓縮成bz2的壓縮文件，首先需要對(duì)文件進(jìn)行解壓。系統(tǒng)集成了Windows的封裝解壓程序bunzip2.exe文件，該文件可通過WinExec()函數(shù)進(jìn)行調(diào)用，達(dá)到解壓bz2文件的目的。調(diào)用方法為：WinExec(./bunzip2.exe-kf *.bz2，SW_HIDE)，其中‘-kf’表示解壓后保留原有壓縮文件并覆蓋同名文件，SW_HIDE表示隱藏調(diào)用方式，函數(shù)返回范圍值大于31表示調(diào)用成功。

解壓后的標(biāo)準(zhǔn)格式雷達(dá)數(shù)據(jù)為bin二進(jìn)制文件，分為多個(gè)區(qū)塊，每個(gè)區(qū)塊描述一組信息。數(shù)據(jù)可分為公共數(shù)據(jù)塊和徑向數(shù)據(jù)塊兩部分，其中公共數(shù)據(jù)塊用于提供站點(diǎn)信息和任務(wù)配置等公共信息。徑向數(shù)據(jù)塊用于存儲(chǔ)天氣雷達(dá)的探測(cè)資料，包括3個(gè)子塊:徑向頭、徑向數(shù)據(jù)頭和徑向數(shù)據(jù)。編寫程序時(shí)，每個(gè)數(shù)據(jù)塊對(duì)應(yīng)一個(gè)結(jié)構(gòu)體，再按照字節(jié)數(shù)整塊讀取。將所有的數(shù)據(jù)塊結(jié)構(gòu)體和讀取函數(shù)封裝成1個(gè)C++類即可。讀取完成后，所有產(chǎn)品數(shù)據(jù)都保存在對(duì)應(yīng)的1個(gè)三維數(shù)組中，以便處理數(shù)據(jù)和畫圖時(shí)使用。

2.2 數(shù)據(jù)質(zhì)量分析

天氣雷達(dá)數(shù)據(jù)質(zhì)量直接關(guān)系到雷達(dá)能否充分發(fā)揮其總體效益，直接影響短臨預(yù)報(bào)的及時(shí)性和準(zhǔn)確性。系統(tǒng)能夠?qū)崟r(shí)監(jiān)測(cè)雷達(dá)數(shù)據(jù)質(zhì)量，如發(fā)現(xiàn)異常則報(bào)警并生成告警文件。

2.2.1 同頻干擾識(shí)別

同頻干擾現(xiàn)象是由于本站周邊存在頻點(diǎn)相同或相近的電磁波發(fā)射設(shè)備造成電磁波干擾引起的[3]。四川全省共有9部S波段天氣雷達(dá)，其中7部雷達(dá)頻點(diǎn)都為2.88 GHz，相鄰最近距離為120 km左右。如成都雷達(dá)站北方會(huì)受到綿陽雷達(dá)的同頻干擾，回波圖形上會(huì)出現(xiàn)放射狀的弧形回波，強(qiáng)度能達(dá)到40 dBZ以上，嚴(yán)重影響了成都雷達(dá)的數(shù)據(jù)質(zhì)量和雷達(dá)業(yè)務(wù)運(yùn)行。

識(shí)別同頻干擾圖形需要對(duì)該異常回波進(jìn)行特征提取，經(jīng)統(tǒng)計(jì)發(fā)現(xiàn)，同頻干擾主要出現(xiàn)在低層仰角，一般在3.4°以下最為明顯；干擾回波圖形呈放射狀的弧形狀，類似于扇葉的形狀；干擾回波較強(qiáng)，最大值超過40 dBZ，徑向廓線雜亂，而氣象回波廓線具有連續(xù)性，同頻干擾回波和氣象回波的徑向廓線特征如圖2所示；同頻干擾回波的相關(guān)系數(shù)CC在強(qiáng)干擾帶為1.04，在弱干擾帶低于0.9，均值低于0.5。

圖2 氣象回波與同頻干擾徑向廓線特征圖

圖3給出了2021-09-01系統(tǒng)識(shí)別到的1次成都雷達(dá)同頻干擾回波，當(dāng)日成都無降水，成都雷達(dá)西北方向出現(xiàn)較強(qiáng)回波，反射率因子最大值為64 dBZ，強(qiáng)干擾帶呈現(xiàn)弧形分布，且強(qiáng)度高于40 dBZ，相關(guān)系數(shù)多為1.04。當(dāng)日系統(tǒng)識(shí)別到同頻干擾84次，雷達(dá)保障人員通過對(duì)硬件設(shè)備和濾波器的改良后，該異常回波問題得到解決。

圖3 成都雷達(dá)同頻干擾回波

2.2.2 地物雜波干擾

地物雜波干擾主要是由于雷達(dá)周邊山體、高樓、高塔，因反射雷達(dá)部分主波或者旁瓣信號(hào)而造成的。這是一個(gè)不可避免的現(xiàn)象，針對(duì)地物雜波而研究的抑制和消除方法已經(jīng)比較成熟，經(jīng)過處理后的雷達(dá)數(shù)據(jù)受到的地物雜波干擾主要是一些新型的地物，如風(fēng)力發(fā)電機(jī)等或由于系統(tǒng)故障未進(jìn)行地物濾波導(dǎo)致的[4]。

因此，系統(tǒng)識(shí)別地物雜波干擾的方法主要是針對(duì)上述原因：一是針對(duì)未被濾除的地物回波，雷達(dá)回波圖長(zhǎng)時(shí)間存在某一固定位置且面積較大的強(qiáng)回波區(qū)域，則判斷為地物雜波干擾；二是針對(duì)系統(tǒng)故障未進(jìn)行地物濾波的情況，該情況雷達(dá)回波圖上會(huì)出現(xiàn)大面積強(qiáng)度大，而速度為0 m/s左右的回波區(qū)域。圖4給出了2020-09-08系統(tǒng)識(shí)別到的1次廣元雷達(dá)地物雜波干擾，可見雷達(dá)北方10 km處存在1個(gè)強(qiáng)回波區(qū)域，強(qiáng)度超過60 dBZ，速度在-11 m/s左右，無法滿足固定地物速度在0 m/s左右的判定，所以雷達(dá)系統(tǒng)無法濾除，經(jīng)過與當(dāng)?shù)乩走_(dá)站聯(lián)系后，發(fā)現(xiàn)這是一片風(fēng)力發(fā)電機(jī)群，風(fēng)葉旋轉(zhuǎn)時(shí)雷達(dá)能觀測(cè)到有回波和速度，雷達(dá)系統(tǒng)無法濾波抑制，該問題暫時(shí)未得到解決。

圖4 廣元雷達(dá)地物干擾回波

2.2.3 雷達(dá)系統(tǒng)異常

系統(tǒng)能識(shí)別由于雷達(dá)系統(tǒng)異常而導(dǎo)致的數(shù)據(jù)質(zhì)量異常問題，包括因標(biāo)定問題導(dǎo)致的雷達(dá)回波偏強(qiáng)或者偏弱、因雷達(dá)伺服系統(tǒng)標(biāo)定異常導(dǎo)致的雷達(dá)仰角和方位角誤差過大、因雷達(dá)高壓未開導(dǎo)致的雷達(dá)數(shù)據(jù)為空回波等問題。系統(tǒng)能通過識(shí)別和監(jiān)控天氣雷達(dá)地物回波位置和面積的方法，對(duì)上述問題進(jìn)行識(shí)別，識(shí)別流程如圖5所示。

根據(jù)中國氣象局相關(guān)規(guī)定，在對(duì)天氣雷達(dá)進(jìn)行周、月、年維護(hù)時(shí)，需要對(duì)雷達(dá)參數(shù)進(jìn)行標(biāo)定，如雷達(dá)未進(jìn)行正確標(biāo)定或者參數(shù)存在異常跳變后，回波強(qiáng)度會(huì)隨之出現(xiàn)較大跳變。在雷達(dá)系統(tǒng)運(yùn)行正常的情況下，雷達(dá)的地物回波強(qiáng)度和面積不應(yīng)存在較大變化，系統(tǒng)通過算法識(shí)別出地物回波后，通過監(jiān)控地物回波的強(qiáng)度變化間接反映雷達(dá)回波強(qiáng)度的準(zhǔn)確性[5]。通過統(tǒng)計(jì)發(fā)現(xiàn)，雷達(dá)地物回波的強(qiáng)度變化范圍在3 dBZ內(nèi)，平均值低于1.4 dBZ。系統(tǒng)以3 dBZ為閾值，當(dāng)雷達(dá)地物回波突變量超過該閾值時(shí)進(jìn)行報(bào)警，提示監(jiān)控人員雷達(dá)回波強(qiáng)度可能出現(xiàn)偏差過大的問題。

圖5 雷達(dá)系統(tǒng)異常回波識(shí)別流程圖

當(dāng)雷達(dá)伺服系統(tǒng)標(biāo)定異常時(shí)，會(huì)導(dǎo)致雷達(dá)仰角和方位角存在較大偏差，雷達(dá)回波位置會(huì)出現(xiàn)偏移，影響預(yù)報(bào)員預(yù)報(bào)精確性。通過監(jiān)控雷達(dá)地物回波的位置和面積大小，可識(shí)別該異常回波。當(dāng)?shù)匚锘夭ǖ奈恢冒l(fā)生整體變化時(shí)，識(shí)別為雷達(dá)伺服系統(tǒng)方位角誤差過大；當(dāng)?shù)脱鼋堑匚锘夭娣e變大或者變小時(shí)，識(shí)別為雷達(dá)伺服系統(tǒng)仰角誤差過大。

當(dāng)雷達(dá)因高壓故障時(shí)，雷達(dá)還在繼續(xù)體掃的情況，這樣的雷達(dá)數(shù)據(jù)里并沒有掃描內(nèi)容，稱為空?qǐng)D異常，圖像和晴空雷達(dá)圖像相似，觀測(cè)員很難發(fā)現(xiàn)雷達(dá)故障，同時(shí)也影響預(yù)報(bào)員對(duì)天氣的判斷。系統(tǒng)監(jiān)測(cè)到雷達(dá)地物回波和其他所有產(chǎn)品都為空值時(shí)，則判斷為空?qǐng)D異常，發(fā)出告警信息。

2.3 剖面分析

由于雷達(dá)數(shù)據(jù)僅有9個(gè)仰角層，天線發(fā)射的電磁波脈沖寬度在1°左右，直接剖面出來的圖像結(jié)構(gòu)太過離散，所以應(yīng)對(duì)剖面數(shù)據(jù)進(jìn)行差值處理，讓剖面圖像形成一個(gè)連續(xù)的整體，這樣的剖面圖像更為直觀，更有利于分析回波的垂直結(jié)構(gòu)。

常用的差值方法有3種，分別為線性差值、雙線性差值和Barnes差值。線性差值較為簡(jiǎn)單，取上下兩個(gè)仰角的反射率因子插值。而雙線性差值則多了一組水平方向上的數(shù)據(jù)插值[6]。

Barnes插值法是一種指數(shù)權(quán)重插值法，十分適合在球坐標(biāo)系中使用，基本公式為：

(1)

式中，w表示相關(guān)插值點(diǎn)的權(quán)重系數(shù)；R0，θ0，Φ0表示相關(guān)插值點(diǎn)的徑向距離、方位角和仰角；Rg，θg，Φg表示待插值點(diǎn)的坐標(biāo)；Kr，Kθ，KΦ表示平滑系數(shù)，文章取1。

(2)

式中，fg表示待插值點(diǎn)的插值結(jié)果；f0(k)表示待插值點(diǎn)相關(guān)區(qū)域內(nèi)，第k個(gè)距離庫的反射率值；N表示待插值點(diǎn)的相關(guān)區(qū)域內(nèi)的距離庫數(shù)。

公式(1)表示各個(gè)插值相關(guān)點(diǎn)的權(quán)重系數(shù)，公式(2)決定待插值點(diǎn)的值。圖6給出了分別采用上述3種插值法得到的垂直剖面效果，可見線性插值和雙線性插值法在近處易出現(xiàn)空缺，在仰角層交界的區(qū)域容易出現(xiàn)空點(diǎn)，連續(xù)性和整體性都不如Barnes插值的結(jié)果，插值得到的圖像效果最好的是Barnes插值法。

圖6 RHI對(duì)比圖

2.4 格點(diǎn)化

由于在識(shí)別異常回波和配置地圖文件時(shí)，都需要對(duì)回波位置進(jìn)行準(zhǔn)確定位，需要按照經(jīng)緯度對(duì)雷達(dá)回波進(jìn)行格點(diǎn)化操作。

地圖文件格式為SHP二進(jìn)制文件，讀取方式與雷達(dá)數(shù)據(jù)類似，通過格式說明文件編寫對(duì)應(yīng)代碼即可，得到經(jīng)緯度信息后需要轉(zhuǎn)到屏幕坐標(biāo)才能進(jìn)行畫圖和定位[7]。另外雷達(dá)數(shù)據(jù)僅包含中心點(diǎn)的經(jīng)緯度信息，同樣需要將所有雷達(dá)距離庫的坐標(biāo)轉(zhuǎn)換到屏幕坐標(biāo)才能進(jìn)行畫圖。經(jīng)過上述轉(zhuǎn)換后，屏幕坐標(biāo)和經(jīng)緯度則有了對(duì)應(yīng)關(guān)系，即可對(duì)雷達(dá)回波進(jìn)行準(zhǔn)確定位。

由于經(jīng)緯度坐標(biāo)的起始原點(diǎn)和屏幕坐標(biāo)的原點(diǎn)不一致(圖7)，通過縱向坐標(biāo)和橫向坐標(biāo)的相似比，即可得到對(duì)應(yīng)關(guān)系，最后引入經(jīng)緯度和屏幕的大小，得到如下公式：

(3)

式中，x和y表示屏幕坐標(biāo)；lon和lat表示經(jīng)度和緯度；W和H表示屏幕的寬度和高度；文章中W=H=400；lonmax和lonmin表示經(jīng)度范圍的上限和下限；latmax和latmin表示緯度范圍的上限和下限。

圖7 地圖坐標(biāo)轉(zhuǎn)換屏幕坐標(biāo)示意圖

3 系統(tǒng)的實(shí)現(xiàn)

系統(tǒng)分為兩個(gè)部分，分別為數(shù)據(jù)獲取部分和數(shù)據(jù)顯示部分。

數(shù)據(jù)獲取部分可以選擇不同的站點(diǎn)和不同的數(shù)據(jù)獲取途徑，分為在線和本地。在線模式可以實(shí)時(shí)獲取近期服務(wù)器中存儲(chǔ)的標(biāo)準(zhǔn)格式雷達(dá)數(shù)據(jù)，并進(jìn)行實(shí)時(shí)更新顯示[8]；本地模式可展示計(jì)算機(jī)中存儲(chǔ)的所有雷達(dá)數(shù)據(jù)，并進(jìn)行歷史回放。

數(shù)據(jù)顯示部分，顯示用戶打開的標(biāo)準(zhǔn)格式雷達(dá)數(shù)據(jù)，可以選擇不同的數(shù)據(jù)產(chǎn)品內(nèi)容，不同仰角和距離圈。在圖上可以進(jìn)行垂直剖面、任意方向剖面、截圖和地圖加載等功能的操作。

4 系統(tǒng)運(yùn)行評(píng)估

通過對(duì)新一代天氣雷達(dá)標(biāo)準(zhǔn)格式數(shù)據(jù)可視化監(jiān)控系統(tǒng)的開發(fā)，實(shí)現(xiàn)了四川省12部新一代天氣雷達(dá)標(biāo)準(zhǔn)格式數(shù)據(jù)的傳輸時(shí)效和數(shù)據(jù)質(zhì)量的可視化監(jiān)控。天氣雷達(dá)站監(jiān)控人員通過該系統(tǒng)能及時(shí)有效地獲取到雷達(dá)數(shù)據(jù)的傳輸狀態(tài)、數(shù)據(jù)質(zhì)量和雷達(dá)回波剖面等信息，雷達(dá)標(biāo)準(zhǔn)格式數(shù)據(jù)通過流傳輸?shù)姆绞剑鄬?duì)于原有方式，傳輸及時(shí)性提到10 min以上。為了檢測(cè)系統(tǒng)的穩(wěn)定性，將系統(tǒng)部署到全省12個(gè)雷達(dá)站，并在本地增加50組虛擬用戶進(jìn)行1個(gè)月的系統(tǒng)穩(wěn)定性測(cè)試。網(wǎng)絡(luò)測(cè)試環(huán)境基于四川省氣象局內(nèi)網(wǎng)20～50 M的寬帶，操作系統(tǒng)基于WINDOWS-10，CPU性能不低于INTEL I3，內(nèi)存不低于4 G，測(cè)試結(jié)果見表1，系統(tǒng)運(yùn)行穩(wěn)定，用戶體驗(yàn)良好。

表1 系統(tǒng)頁面響應(yīng)測(cè)試記錄表

5 結(jié)束語

文章設(shè)計(jì)的新一代天氣雷達(dá)標(biāo)準(zhǔn)格式數(shù)據(jù)可視化監(jiān)控系統(tǒng)實(shí)現(xiàn)了對(duì)四川12部新一代天氣雷達(dá)標(biāo)準(zhǔn)格式數(shù)據(jù)傳輸?shù)谋O(jiān)控，并能實(shí)時(shí)監(jiān)控?cái)?shù)據(jù)質(zhì)量，系統(tǒng)自2021年初試運(yùn)行以來，僅出現(xiàn)極少BUG，并得到及時(shí)修復(fù)，如今運(yùn)行穩(wěn)定。

雷達(dá)數(shù)據(jù)質(zhì)量不高影響了天氣雷達(dá)數(shù)據(jù)傳輸?shù)臏?zhǔn)確性，文章設(shè)計(jì)開發(fā)的標(biāo)準(zhǔn)格式數(shù)據(jù)可視化監(jiān)控系統(tǒng)能夠發(fā)現(xiàn)一些常見的數(shù)據(jù)質(zhì)量問題，但雷達(dá)故障點(diǎn)較多，從數(shù)據(jù)上能反映出的問題有限，繼續(xù)開發(fā)天氣雷達(dá)數(shù)據(jù)質(zhì)量異常回波的識(shí)別方法，還需進(jìn)一步研究。