魯峻麟，黃惺惺，顧桃峰，張志堅(jiān)，賈寧遠(yuǎn)

（1.廣州市突發(fā)事件預(yù)警信息發(fā)布中心，廣州 511430；2.廣東省氣象公共安全技術(shù)支持中心，廣州 510641）

引言

氣象領(lǐng)域傳統(tǒng)的測風(fēng)手段有風(fēng)杯風(fēng)速計(jì)、常規(guī)探空氣球等［1］，隨著測風(fēng)手段的更新以及對(duì)氣象設(shè)備認(rèn)同度的提升，風(fēng)廓線雷達(dá)也逐漸被氣象部門所接納［2－3］，成為日常的風(fēng)場測量設(shè)備。傳統(tǒng)的風(fēng)速計(jì)和風(fēng)向標(biāo)只能實(shí)現(xiàn)單點(diǎn)測量，借助測風(fēng)塔后實(shí)現(xiàn)對(duì)應(yīng)高度層的風(fēng)場信息檢測，但這類傳統(tǒng)裝置易受冰凍天氣影響，測風(fēng)塔的搭建和維護(hù)也需要花費(fèi)大量的人力物力，還存在移動(dòng)困難和前期征地手續(xù)復(fù)雜等問題［4－7］。微波雷達(dá)常用波長主要為厘米波，與大氣中的大尺寸粒子（如云、雨、冰等）相互作用產(chǎn)生回波，無法與大氣中的分子或氣溶膠顆粒產(chǎn)生作用，而晴空時(shí)大氣中大尺寸粒子較少，因此微波雷達(dá)在晴空天氣條件下將出現(xiàn)探測盲區(qū)，此外微波雷達(dá)還具備龐大的收發(fā)系統(tǒng)也導(dǎo)致其移動(dòng)的困難［8－9］。因此，迫切需要補(bǔ)充新型的風(fēng)場測量手段替代傳統(tǒng)測風(fēng)裝置實(shí)現(xiàn)大氣風(fēng)場信息的測量。

掃描型多普勒測風(fēng)激光雷達(dá)系統(tǒng)是一種小型、全自動(dòng)、無環(huán)境電磁干擾的多普勒激光雷達(dá)，實(shí)現(xiàn)全天候、不間斷的三維低空風(fēng)場探測［10－11］，通過配置多臺(tái)激光測風(fēng)雷達(dá)組網(wǎng)觀測業(yè)務(wù)軟件平臺(tái)，實(shí)現(xiàn)風(fēng)廓線、垂直氣流及PPI／RHI／水平風(fēng)場的組網(wǎng)探測及業(yè)務(wù)分析功能。測風(fēng)激光雷達(dá)采用人眼安全1550nm激光波段、近紅外激光探測等技術(shù)，具有時(shí)空分辨率高、自動(dòng)化程度高、安裝簡單易維護(hù)、移動(dòng)便攜性好等優(yōu)勢，因此成為了最具前景的風(fēng)場信息測量手段。近年來激光測風(fēng)雷達(dá)已在氣象、風(fēng)電、民航和環(huán)保等領(lǐng)域有越來越多的成熟應(yīng)用案例［12－14］，在多種不同的氣象環(huán)境中驗(yàn)證了雷達(dá)風(fēng)場探測的可靠性和準(zhǔn)確性。

在廣州超大城市綜合觀測建設(shè)中，需連續(xù)空間上觀測城市上空邊界層內(nèi)風(fēng)分布狀況，獲取典型天氣過程特征風(fēng)數(shù)據(jù)。本文通過對(duì)比分析3 部不同測風(fēng)雷達(dá)的探測效能，為綜合觀測建設(shè)提供依據(jù)，提高廣州超大城市綜合氣象觀測資料的本地化應(yīng)用水平，為探索建立業(yè)務(wù)化應(yīng)用提供參考依據(jù)。

1 觀測設(shè)備及數(shù)據(jù)來源

1.1 觀測試驗(yàn)及設(shè)備性能

自2019 年廣州超大城市綜合氣象觀測試驗(yàn)啟動(dòng)以來，為驗(yàn)證超大城市綜合觀測試驗(yàn)設(shè)備觀測數(shù)據(jù)的可靠性，選擇3 部不同的測風(fēng)雷達(dá)與清遠(yuǎn)探空站資料進(jìn)行對(duì)比觀測，為保護(hù)廠家信息，將3部測風(fēng)激光雷達(dá)分別以L、W、Y 表示，各性能參數(shù)見表1。

表1 三部測風(fēng)激光雷達(dá)性能參數(shù)

1.2 數(shù)據(jù)來源及處理

選取2020 年6 月1 日—7 月6 日的觀測資料作為有效對(duì)比資料。在對(duì)比評(píng)估中主要基于探空資料，把三部測風(fēng)雷達(dá)的可探測高度、垂直分辨率作為參考，分析評(píng)估觀測高度上的風(fēng)向、風(fēng)速，包括數(shù)據(jù)可用率（非缺測率）、風(fēng)向／風(fēng)速準(zhǔn)確率、風(fēng)向／風(fēng)速平均絕對(duì)誤差以及風(fēng)向／風(fēng)速均方根誤差。具體計(jì)算公式如下：

式中，xi表示某一高度測風(fēng)雷達(dá)觀測值，yi表示對(duì)應(yīng)高度的探空觀測值，m、n 表示該高度上有效樣本數(shù)。

2 觀測試驗(yàn)數(shù)據(jù)分析

2.1 數(shù)據(jù)可用率

根據(jù)觀測資料，統(tǒng)計(jì)觀測日期內(nèi)探空施放球時(shí)間（北京時(shí)07 時(shí)、19 時(shí)）內(nèi)3 部測風(fēng)雷達(dá)的數(shù)據(jù)可用率，得出100m 以下L 雷達(dá)探測效能最優(yōu)，數(shù)據(jù)可用率達(dá)100%；W 雷達(dá)次之，數(shù)據(jù)可用率為80%；Y測風(fēng)雷達(dá)數(shù)據(jù)可用率僅75%。100～1000m，L 和Y 雷達(dá)相當(dāng)，數(shù)據(jù)可用率均為100%，W 雷達(dá)從80%逐漸下降到70%。1000～2000m，Y 雷達(dá)最優(yōu)，數(shù)據(jù)可用率在90%以上，L 雷達(dá)數(shù)據(jù)可用率從90%逐漸下降，W雷達(dá)從70%下降到55%左右。2000m 以上，L 雷達(dá)數(shù)據(jù)可用率從70%逐漸下降，W 雷達(dá)從50%逐漸下降，Y 雷達(dá)數(shù)據(jù)可用率為0。

2.2 風(fēng)向數(shù)據(jù)分析

基于3 部測風(fēng)雷達(dá)探測高度，提取探空數(shù)據(jù)接近高度上的風(fēng)向數(shù)據(jù)作為對(duì)比值，對(duì)測風(fēng)雷達(dá)風(fēng)向數(shù)據(jù)的準(zhǔn)確率、絕對(duì)平均誤差以及均方根誤差進(jìn)行分析。如果測風(fēng)雷達(dá)所觀測風(fēng)向與探空風(fēng)向偏差小于10°則認(rèn)為測風(fēng)雷達(dá)數(shù)據(jù)準(zhǔn)確。

2.2.1 風(fēng)向觀測數(shù)據(jù)準(zhǔn)確率分析

根據(jù)3 部測風(fēng)雷達(dá)風(fēng)向觀測數(shù)據(jù)和清遠(yuǎn)站探空數(shù)據(jù)，分析3 部測風(fēng)雷達(dá)風(fēng)向觀測數(shù)據(jù)準(zhǔn)確率，如圖1 所示。

從圖1 可看出，L 測風(fēng)雷達(dá)風(fēng)向觀測數(shù)據(jù)準(zhǔn)確率分布較不穩(wěn)定，在統(tǒng)計(jì)樣本量較小的高度范圍（1500～3000m）易出現(xiàn)準(zhǔn)確率為0 的探測結(jié)果；在統(tǒng)計(jì)樣本量較大的高度上，200m 以下高度準(zhǔn)確性較差，未超過60%，高層準(zhǔn)確率基本維持在80%以上，個(gè)別高度風(fēng)向探測準(zhǔn)確率能達(dá)到100%。W 測風(fēng)雷達(dá)風(fēng)向觀測數(shù)據(jù)準(zhǔn)確率較差，2000m 高度以上層次探測數(shù)據(jù)準(zhǔn)確率基本為0，2000m 高度以下層次準(zhǔn)確率分布不穩(wěn)定，大部分高度上準(zhǔn)確率低于50%。Y測風(fēng)雷達(dá)風(fēng)向觀測數(shù)據(jù)準(zhǔn)確率分布較穩(wěn)定，除在個(gè)別層次由于統(tǒng)計(jì)樣本量較小出現(xiàn)了準(zhǔn)確率為0 的探測結(jié)果（800m 和1000m 高度）之外，在統(tǒng)計(jì)樣本量較大的高度上，風(fēng)向準(zhǔn)確率基本能達(dá)到60%以上。其中，400～1300m 之間是數(shù)據(jù)準(zhǔn)確率較高，某些高度可達(dá)到90%以上的準(zhǔn)確率，個(gè)別高度風(fēng)向探測準(zhǔn)確率能達(dá)到100%。

圖1 三部測風(fēng)雷達(dá)風(fēng)向數(shù)據(jù)準(zhǔn)確率

2.2.2 風(fēng)向觀測數(shù)據(jù)誤差分析

為進(jìn)一步分析測風(fēng)雷達(dá)風(fēng)向觀測數(shù)據(jù)的誤差分布情況，分別統(tǒng)計(jì)3 部測風(fēng)雷達(dá)風(fēng)向觀測數(shù)據(jù)的平均絕對(duì)偏差和均方根誤差隨高度的分布（圖2）。

圖2 三部測風(fēng)雷達(dá)風(fēng)向數(shù)據(jù)平均絕對(duì)誤差和均方根誤差分布

從風(fēng)向探測數(shù)據(jù)看，3 部測風(fēng)雷達(dá)對(duì)中層的風(fēng)向探測誤差最小，L 測風(fēng)雷達(dá)平均絕對(duì)誤差和均方根誤差控制在10°左右；W 測風(fēng)雷達(dá)平均絕對(duì)誤差維持在15°左右。低層、高層探測誤差均較大。其中Y 測風(fēng)雷達(dá)風(fēng)向探測的平均絕對(duì)誤差基本上控制在20°以內(nèi)，除800m 和1000m 高度以外，其余層次平均絕對(duì)誤差基本在10°以內(nèi)，個(gè)別高度只有1～2°的誤差；風(fēng)向探測的均方根誤差除60m 高度均方根誤差接近40°之外，其余高度均方根誤差基本控制在20°以內(nèi)，大部分高度均方根誤差都在10°左右。

比較3 部測風(fēng)雷達(dá)的風(fēng)向準(zhǔn)確率以及誤差分布可以發(fā)現(xiàn)，Y 雷達(dá)除了近地面層準(zhǔn)確率誤差穩(wěn)定性較低以外，其余高度準(zhǔn)確率和誤差穩(wěn)定性均優(yōu)于L測風(fēng)雷達(dá)和W 測風(fēng)雷達(dá)，在高層2000m 以上，L 雷達(dá)較優(yōu)，W 雷達(dá)次之，Y 雷達(dá)最差。

2.3 風(fēng)速數(shù)據(jù)分析

基于3 部測風(fēng)雷達(dá)探測高度，提取探空數(shù)據(jù)接近高度上的風(fēng)速數(shù)據(jù)作為對(duì)比值，對(duì)測風(fēng)雷達(dá)風(fēng)速數(shù)據(jù)的準(zhǔn)確率、絕對(duì)平均誤差以及均方根誤差進(jìn)行了分析。如果測風(fēng)雷達(dá)所觀測風(fēng)速與高度接近的探空風(fēng)速偏差小于1m·s-1則認(rèn)為測風(fēng)雷達(dá)數(shù)據(jù)準(zhǔn)確。

2.3.1 風(fēng)速觀測數(shù)據(jù)準(zhǔn)確率分析

根據(jù)3 部測風(fēng)雷達(dá)風(fēng)速觀測數(shù)據(jù)和清遠(yuǎn)站探空數(shù)據(jù)，分別分析測風(fēng)雷達(dá)風(fēng)速觀測數(shù)據(jù)準(zhǔn)確率（圖3）。

由圖3 可知，L 測風(fēng)雷達(dá)風(fēng)速觀測數(shù)據(jù)準(zhǔn)確率分布與風(fēng)向準(zhǔn)確率基本相同，呈現(xiàn)較不穩(wěn)定的趨勢。在統(tǒng)計(jì)樣本量較小的高度范圍（1500～3000m）易出現(xiàn)準(zhǔn)確率為0 的探測結(jié)果。在統(tǒng)計(jì)樣本量較大的高度上，200m 以下高度準(zhǔn)確性較差，未超過60%，高層準(zhǔn)確率基本維持在80%以上，個(gè)別高度風(fēng)速探測準(zhǔn)確率能達(dá)到100%；W 測風(fēng)雷達(dá)風(fēng)速觀測數(shù)據(jù)準(zhǔn)確率較差，2000m 高度以上層次探測數(shù)據(jù)準(zhǔn)確率基本為0，2000m 高度以下層次準(zhǔn)確率分布不穩(wěn)定，大部分高度上準(zhǔn)確率低于50%；Y 測風(fēng)雷達(dá)風(fēng)速觀測數(shù)據(jù)準(zhǔn)確率分布較穩(wěn)定，除在個(gè)別層次由于統(tǒng)計(jì)樣本量較小出現(xiàn)了準(zhǔn)確率為0%的探測結(jié)果（1100m和1360m 高度）之外，在統(tǒng)計(jì)樣本量較大的高度上，風(fēng)速準(zhǔn)確率基本能維持在60%左右，個(gè)別高度風(fēng)向探測準(zhǔn)確率能達(dá)到100%。

圖3 三部測風(fēng)雷達(dá)風(fēng)速數(shù)據(jù)樣本數(shù)與準(zhǔn)確率

2.3.2 風(fēng)速觀測數(shù)據(jù)誤差分析

為進(jìn)一步分析3 部測風(fēng)雷達(dá)風(fēng)速觀測數(shù)據(jù)的誤差分布情況，分別統(tǒng)計(jì)了3 部測風(fēng)雷達(dá)風(fēng)速觀測數(shù)據(jù)的平均絕對(duì)偏差和均方根誤差所高度的分布，如圖4 所示。

圖4 三部測風(fēng)雷達(dá)風(fēng)速數(shù)據(jù)平均絕對(duì)誤差和均方根誤差分布

從風(fēng)速探測數(shù)據(jù)來看，3 部測風(fēng)雷達(dá)對(duì)中低層的風(fēng)速探測誤差最小，L 測風(fēng)雷達(dá)平均絕對(duì)誤差和均方根誤差分別為0.69m·s-1和0.80m·s-1；W 測風(fēng)雷達(dá)平均絕對(duì)誤差和均方根誤差分別為2.08m·s-1和2.37m·s-1；Y 測風(fēng)雷達(dá)平均絕對(duì)誤差和均方根誤差分別為0.83m·s-1和0.99m·s-1。高層探測誤差均較大，高層探測誤差較大的原因有可能由于探空氣球飄移距離較遠(yuǎn)所致。

比較三部測風(fēng)雷達(dá)的風(fēng)速準(zhǔn)確率以及誤差分布可以發(fā)現(xiàn)，在中低層（2000m 以下）Y 測風(fēng)雷達(dá)和L測風(fēng)雷達(dá)準(zhǔn)確率和探測誤差較為接近，W 雷達(dá)較差。

2.4 逐日誤差分析

通過分析3 部測風(fēng)雷達(dá)逐日所有觀測高度的平均絕對(duì)誤差（MAE）評(píng)估測風(fēng)雷達(dá)在觀測期間的穩(wěn)定性。從風(fēng)向的絕對(duì)平均誤差值來看，W 測風(fēng)雷達(dá)誤差明顯高于L 和Y 測風(fēng)雷達(dá)，L 和Y 測風(fēng)雷達(dá)風(fēng)向觀測誤差較為接近。

從風(fēng)速的絕對(duì)平均誤差值來看，W 測風(fēng)雷達(dá)誤差明顯高于L 和Y 測風(fēng)雷達(dá)，L 和Y 測風(fēng)雷達(dá)風(fēng)速觀測誤差較為接近，但在觀測前期Y 測風(fēng)雷達(dá)風(fēng)速觀測還是略低于L 測風(fēng)雷達(dá)風(fēng)速觀測。

3 結(jié)論與討論

（1）在數(shù)據(jù)可用率方面，Y 測風(fēng)雷達(dá)數(shù)據(jù)可用率除60m 高度較低外，其余高度數(shù)據(jù)可用率均達(dá)到90%以上，在同樣高度層次上優(yōu)于L 測風(fēng)雷達(dá)數(shù)據(jù)可用率和W 測風(fēng)雷達(dá)數(shù)據(jù)可用率；其中W 測風(fēng)雷達(dá)數(shù)據(jù)可用率表現(xiàn)最差。2000m 以下Y 雷達(dá)數(shù)據(jù)最優(yōu)，L 雷達(dá)次之，W 雷達(dá)最次，2000m 以上L 雷達(dá)數(shù)據(jù)最優(yōu)，W 雷達(dá)次之，Y 雷達(dá)基本為0。

（2）從測風(fēng)雷達(dá)的風(fēng)向探測數(shù)據(jù)來看，3 部測風(fēng)雷達(dá)對(duì)中層的風(fēng)向探測誤差最小，低層次之，高層探測誤差均較大。比較3 部測風(fēng)雷達(dá)的風(fēng)向準(zhǔn)確率可以發(fā)現(xiàn)，在中層（2000m 以下），L 雷達(dá)觀測數(shù)據(jù)的平均絕對(duì)誤差和均方根誤差控制在10°左右，L 雷達(dá)與Y 雷達(dá)較優(yōu)，W 雷達(dá)較差。

（3）從測風(fēng)雷達(dá)的風(fēng)速準(zhǔn)確率以及誤差分布來看，在中層（2000m 以下），L 雷達(dá)和Y 雷達(dá)較優(yōu)，W雷達(dá)次之。在低層（1000m 以下）Y 測風(fēng)雷達(dá)和L 測風(fēng)雷達(dá)準(zhǔn)確率和探測誤差較為接近，探測效果最好；W 測風(fēng)雷達(dá)對(duì)風(fēng)速的探測效果最差。

（4）從逐日三部測風(fēng)雷達(dá)的風(fēng)向、風(fēng)速平均絕對(duì)誤差分布可以看出，L 和Y 雷達(dá)的逐日風(fēng)速、風(fēng)向誤差較為接近，但Y 雷達(dá)觀測誤差還是略低于L 雷達(dá)，W 測風(fēng)雷達(dá)誤差明顯高于L 和Y 測風(fēng)雷達(dá)。

本研究所取的時(shí)段較短，樣本較少，研究結(jié)果可能有一定的局限性。且由于探空氣球在上升過程中會(huì)存在飄移情況，而測風(fēng)雷達(dá)進(jìn)行垂直探測，因此在對(duì)高層數(shù)據(jù)的評(píng)估中無法剔除由于探空氣球飄移動(dòng)所產(chǎn)生的誤差。