摘要 為了提高簡(jiǎn)便、快捷、自動(dòng)識(shí)別冰雹云的能力，利用雷達(dá)拼圖組合反射率因子（Combined Reflectivity，CR）數(shù)據(jù)，在冰雹回波超級(jí)單體四要素特征分析的基礎(chǔ)上，提出了簡(jiǎn)約快捷的自動(dòng)識(shí)別冰雹云的方法。結(jié)果表明：聚類算法（Clustering Algorithm）、散點(diǎn)輪廓算法（Scatter Contour Algorithm）能夠較好地識(shí)別出冰雹云回波中心強(qiáng)度和強(qiáng)回波面積;強(qiáng)回波梯度算法（Strong Echo Gradient Algorithm）、云砧回波算法（Cloud Anvil Echo Algorithm）計(jì)算快捷。雷達(dá)拼圖CR、強(qiáng)回波面積（Strong Echo Area，SEA）、強(qiáng)回波梯度（Strong Echo Gradient，SEG）和云砧回波（Cloud Anvil Echo，CAE）四要素被用來(lái)確定回波與冰雹云的關(guān)系。江西冰雹大多數(shù)發(fā)生在超級(jí)單體（Supercell）中，當(dāng)CR≥60 dBZ、SEA≥100 km2、SEG≤8 km、CAE比值在1∶2～1∶3時(shí)，就可能發(fā)生冰雹;有些微型超級(jí)單體（Micro Supercell）在合適的天氣背景和環(huán)境條件下，即使SEA=18 km2也會(huì)發(fā)生冰雹。自動(dòng)識(shí)別冰雹云的方法在2022、2023年各3次冰雹過程中得到實(shí)踐驗(yàn)證，其識(shí)別出的冰雹云區(qū)域與冰雹實(shí)況區(qū)域吻合，但也存在10％～20％的空?qǐng)?bào)率。本研究結(jié)果為簡(jiǎn)便、快捷、自動(dòng)識(shí)別冰雹天氣提供了有效依據(jù)。

關(guān)鍵詞冰雹云;聚類算法;散點(diǎn)輪廓算法;強(qiáng)回波梯度;云砧回波

2023-08-01收稿，2024-03-29接受

北極閣基金（BJG202208）;國(guó)家重點(diǎn)研發(fā)計(jì)劃課題（2022YFC3003904）;中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)（A類-XDA19040202）;國(guó)家自然科學(xué)基金項(xiàng)目（41975001）;江西省氣象局重點(diǎn)科研項(xiàng)目（JX2022Z04）;2022年江西省氣象局面上項(xiàng)目（JX2022M03）;景德鎮(zhèn)市科技計(jì)劃項(xiàng)目（2022SF003）

引用格式：馬中元，王金鑫，張林才，等，2024.基于雷達(dá)拼圖CR產(chǎn)品四要素識(shí)別冰雹云的方法［J］.大氣科學(xué)學(xué)報(bào)，47（5）：826-840.

Ma Z Y，Wang J X，Zhang L C，et al.，2024.A method for identifying hail clouds based on four elements of radar mosaic CR products［J］.Trans Atmos Sci，47（5）：826-840.doi：10.13878／j.cnki.dqkxxb.20230801001.（in Chinese）.

冰雹天氣是一種短時(shí)強(qiáng)對(duì)流天氣，冰雹發(fā)生時(shí)常伴隨狂風(fēng)、雷電、強(qiáng)降水、急劇降溫等天氣現(xiàn)象，不僅對(duì)農(nóng)業(yè)、電力、交通等造成巨大的損失，還對(duì)人民生命財(cái)產(chǎn)造成極大的威脅。冰雹發(fā)生過程歷時(shí)較短且具有不確定性和局地性等特點(diǎn)。因此，簡(jiǎn)便、快捷、自動(dòng)識(shí)別冰雹云的算法研究及便于普及的監(jiān)測(cè)預(yù)警產(chǎn)品平臺(tái)的搭建，一直是雷達(dá)氣象重點(diǎn)研究課題。

1981年Forsyth等提出了WSR-88D冰雹云探測(cè)方法，由于WSD-88D算法是半經(jīng)驗(yàn)統(tǒng)計(jì)算法，所以在不同的氣候和地形影響下算法會(huì)出現(xiàn)誤差（王瑾和劉黎平，2011）。隨著雷達(dá)探測(cè)技術(shù)的發(fā)展，新一代天氣雷達(dá)CINRAD／SA可以獲得大量數(shù)據(jù)，對(duì)冰雹天氣的監(jiān)測(cè)和分析起到了重要作用。目前國(guó)內(nèi)冰雹識(shí)別方法研究較多，如：利用雷達(dá)回波強(qiáng)度和平均零度層高度等特點(diǎn)進(jìn)行識(shí)別（李金輝和樊鵬，2007）;根據(jù)強(qiáng)冰雹發(fā)生的特征指標(biāo)，將垂直積分液態(tài)水含量（Vertical Integral Liquid water content，VIL）、回波頂輻散、三體散射長(zhǎng)釘（吳劍坤和俞小鼎，2009）、負(fù)溫區(qū)回波厚度、強(qiáng)回波梯度（李德俊等，2011）、回波組合反射率（段鶴等，2014）等作為冰雹預(yù)警指標(biāo)。上述方法對(duì)冰雹云的識(shí)別都有重要影響，但距離實(shí)現(xiàn)冰雹云自動(dòng)識(shí)別和識(shí)別產(chǎn)品自動(dòng)處理的要求還有一定差距。

隨著科技發(fā)展，人工智能、機(jī)器學(xué)習(xí)和圖像識(shí)別等方法被不斷運(yùn)用于冰雹云的識(shí)別。路志英等（2015）基于圖像處理方法，對(duì)雷達(dá)反射率因子剖面圖進(jìn)行了研究。張文海和李磊（2019）應(yīng)用貝葉斯分類法、機(jī)器學(xué)習(xí)和人工智能方法來(lái)識(shí)別冰雹云。李冰村等（2022）根據(jù)支持向量機(jī)、決策樹和樸素貝葉斯方法對(duì)冰雹進(jìn)行分類識(shí)別及預(yù)報(bào)，而劉新偉等（2021）則根據(jù)隨機(jī)森林模型對(duì)冰雹進(jìn)行分類識(shí)別及預(yù)報(bào)。張秉祥等（2014）將模糊邏輯算法應(yīng)用于冰雹天氣的識(shí)別，其識(shí)別命中率、虛警率和臨界成功指數(shù)分別為73.9％、36.4％和51.9％。李博勇等（2021）利用貝葉斯方法改進(jìn)了華南地區(qū)冰雹的識(shí)別效果，并比較了WSR-88D冰雹識(shí)別算法和貝葉斯方法對(duì)冰雹識(shí)別的效果;結(jié)果表明，這兩種方法都能較準(zhǔn)確地識(shí)別出冰雹云，其中貝葉斯方法識(shí)別范圍較大。盡管這些方法已在冰雹云識(shí)別中得到應(yīng)用，但這些研究基于不同的天氣雷達(dá)數(shù)據(jù)，其選取的特征量存在較大差異。這些方法識(shí)別范圍較大、識(shí)別準(zhǔn)確率較低、算法效率較低，遇到小范圍的冰雹區(qū)域時(shí)，其識(shí)別效果并不理想。冰雹回波的傳統(tǒng)識(shí)別方法是從回波強(qiáng)度、強(qiáng)回波面積、強(qiáng)回波梯度等方面進(jìn)行分析。在雷達(dá)產(chǎn)品終端（PUP）上，還可以進(jìn)行回波頂高度、強(qiáng)回波頂高、特殊特征（懸掛、回波墻、勾指狀、三體等）、垂直積分液態(tài)水含量VIL、中氣旋、中渦旋、速度對(duì)、中層輻合區(qū)等分析。雙偏振雷達(dá)還可以根據(jù)反射率因子差分反射率（Differential Reflectance，ZDR）、差分傳播相移（Propagation Differential Phase shift，PDP）、相關(guān)系數(shù)（co-polar cross-Correlation Coefficient，CC）等產(chǎn)品進(jìn)行分析。但這些分析存在一定的門檻：一是需要雷達(dá)產(chǎn)品終端條件支持（硬件）;二是對(duì)分析人員的雷達(dá)資料分析能力的要求較高（門檻），同樣資料、不同人員的分析結(jié)果往往不同。這兩點(diǎn)限制了廣大基層臺(tái)站（市、縣兩級(jí)）預(yù)報(bào)人員的使用以及方法的普及。如何僅使用普及的雷達(dá)拼圖組合反射率因子（Combined Reflectivity factor，CR）產(chǎn)品，就能實(shí)現(xiàn)對(duì)冰雹云的自動(dòng)識(shí)別？這是本文研究的目的和意義。

本文擬在對(duì)以往江西冰雹雷達(dá)產(chǎn)品終端（PUP）回波特征和雷達(dá)拼圖組合反射率因子CR回波特征的認(rèn)知基礎(chǔ)上，根據(jù)聚類算法和散點(diǎn)輪廓算法生成回波候選區(qū)域，并根據(jù)不同候選區(qū)域的回波特征與冰雹云之間的內(nèi)在聯(lián)系，建立自動(dòng)識(shí)別方法，形成自動(dòng)化的冰雹云識(shí)別產(chǎn)品平臺(tái)（10 min滾動(dòng)更新），從而實(shí)現(xiàn)為廣大基層臺(tái)站（市、縣兩級(jí)）預(yù)報(bào)人員提供監(jiān)測(cè)預(yù)警冰雹天氣的分析依據(jù)。

1 資料和概況

1.1 資料

雷達(dá)資料來(lái)源于江西雷達(dá)拼圖產(chǎn)品（江西雷達(dá)拼圖有24部S波段多普勒雷達(dá)，其中江西8部、湖南5部、湖北1部、安徽2部、浙江2部、福建3部和廣東3部;2012年完成開發(fā)并上網(wǎng)共享）。雷達(dá)拼圖組合反射率因子（CR，單位：dBZ）是雷達(dá)反射率因子數(shù)據(jù)（9層）在垂直方向上的累計(jì)最大值。雷達(dá)拼圖數(shù)據(jù)文件經(jīng)緯度網(wǎng)格點(diǎn)范圍為110.2°～121.2°E、22°～33°N，空間分辨率為1 km×1 km，時(shí)間分辨率為10 min（簡(jiǎn)稱為1100格式數(shù)據(jù)）;24部雷達(dá)極坐標(biāo)數(shù)據(jù)轉(zhuǎn)換為經(jīng)緯度網(wǎng)格點(diǎn)數(shù)據(jù);全文使用北京時(shí)間（BST）。

冰雹實(shí)況信息來(lái)源于江西基本國(guó)家氣象站記錄，相關(guān)業(yè)務(wù)微信群照片、視頻（需經(jīng)過雷達(dá)回波強(qiáng)度訂正確認(rèn)），冰雹發(fā)生地氣象信息員、愛好者的信息反饋，以及冰雹災(zāi)后氣象部門組織的調(diào)查報(bào)告等。其中最快最有效的途徑是微信群信息。

基于雷達(dá)拼圖CR產(chǎn)品四要素識(shí)別冰雹云的方法可簡(jiǎn)寫成“CR-識(shí)別”法，它根據(jù)雷達(dá)拼圖經(jīng)緯度網(wǎng)格點(diǎn)CR產(chǎn)品數(shù)據(jù)，對(duì)冰雹云的結(jié)構(gòu)和特征通過四要素進(jìn)行量化處理和分析，進(jìn)而識(shí)別出有可能出現(xiàn)冰雹的指示圖和表格數(shù)據(jù)文件?！癈R-識(shí)別”法的優(yōu)點(diǎn)是算法能有效快速地識(shí)別出冰雹云，10 min雷達(dá)拼圖間隔，每次計(jì)算（視回波多少）僅需1～2 min，滿足雷達(dá)拼圖出圖的時(shí)間要求。其不足之處是：存在一定比例的空?qǐng)?bào)率（10％～20％）;對(duì)于少數(shù)CR回波強(qiáng)度小于60 dBZ出現(xiàn)小冰粒的回波單體，可能存在漏報(bào)現(xiàn)象。

1.2 概況

江西冰雹95％以上發(fā)生在超級(jí)單體回波結(jié)構(gòu)中，只有不到5％的冰雹發(fā)生在江西春季（3—4月）有利的環(huán)境條件下，即發(fā)生在55 dBZ強(qiáng)回波單體（非超級(jí)單體）中，且都是1 cm以下小冰粒（應(yīng)冬梅等，2007;陳鮑發(fā)和馬中元，2019;黃中根等，2022;段和平等，2023）。江西超級(jí)單體分為經(jīng)典超級(jí)單體（簡(jiǎn)稱超級(jí)單體）和微型超級(jí)單體兩種，超級(jí)單體出現(xiàn)冰雹的機(jī)會(huì)最大，微型超級(jí)單體只是在3—4月產(chǎn)生1 cm以下較小冰雹。微型超級(jí)單體與強(qiáng)回波單體（55 dBZ）的回波強(qiáng)度差不多，具有60 dBZ時(shí)稱為微型超級(jí)單體，小于60 dBZ時(shí)稱為強(qiáng)回波單體。在合適季節(jié)（江西3—4月）、合適天氣條件下，小部分強(qiáng)回波單體也會(huì)產(chǎn)生小冰雹。超級(jí)單體與微型超級(jí)單體的主要區(qū)別是：回波強(qiáng)度不同，前者60 dBZ回波強(qiáng)度中存有65 dBZ或70 dBZ“回波核”，而后者沒有“回波核”;60 dBZ強(qiáng)回波面積的大小不同，前者≥100 km2，后者≤100 km2。

江西冰雹還會(huì)發(fā)生在多單體風(fēng)暴（在同一個(gè)卷云罩下的多個(gè)超級(jí)單體群）或颮線回波帶上的超級(jí)單體中，只要單體回波達(dá)到超級(jí)單體四要素閾值，都可以被快速識(shí)別出來(lái)。識(shí)別程序每10 min連續(xù)滾動(dòng)運(yùn)行，通過比較前后回波的演變趨勢(shì)，如兩個(gè)時(shí)次強(qiáng)回波面積的增加率，從而提前預(yù)測(cè)冰雹回波的形成。

江西3—4月春季冰雹和7—8月副高冰雹的天氣形勢(shì)不同，識(shí)別冰雹的閾值不同，如春季面積18 km2的微型超級(jí)單體強(qiáng)回波也能產(chǎn)生小冰雹（2022年3月14日安義縣老觀測(cè)站冰雹），因此按季節(jié)確定不同閾值是提高識(shí)別冰雹云識(shí)別率的一個(gè)有效途徑。

利用CR雷達(dá)拼圖對(duì)江西冰雹進(jìn)行識(shí)別其實(shí)有5個(gè)要素：回波強(qiáng)度和回波核、強(qiáng)回波面積（水平面積）、強(qiáng)回波垂直厚度（垂直面積）、強(qiáng)回波梯度、云砧前伸弱回波。但由于江西雷達(dá)拼圖上暫時(shí)沒有實(shí)現(xiàn)等高面顯示，所以強(qiáng)回波垂直面積無(wú)法被計(jì)算。2024年將利用“天衍”雷達(dá)拼圖提供的垂直剖面信息，開展強(qiáng)回波垂直面積識(shí)別實(shí)踐。根據(jù)雷達(dá)拼圖CR四要素原則，識(shí)別冰雹的準(zhǔn)確率可達(dá)80％以上，因存在近20％的空?qǐng)?bào)，所以若增加強(qiáng)回波垂直厚度（垂直面積）綜合識(shí)別，則可有效減少空?qǐng)?bào)。垂直積分液態(tài)水含量VIL大值區(qū)也是識(shí)別冰雹的一個(gè)指標(biāo)，江西冰雹大多數(shù)處在45～60 kg／m2量級(jí)，加入VIL判別，能有效識(shí)別冰雹云和減少空?qǐng)?bào)。

2 研究方法

按照雷達(dá)拼圖上超級(jí)單體識(shí)別四要素來(lái)區(qū)別，超級(jí)單體和微型超級(jí)單體四要素基本相同，都要達(dá)到超級(jí)單體的基本閾值，即：CR圖上具有60 dBZ的回波強(qiáng)度、60 dBZ強(qiáng)回波面積、強(qiáng)回波梯度（30～60 dBZ之間的距離）、具有云砧形成的“前伸”弱回波。超級(jí)單體和微型超級(jí)單體主要區(qū)別在于前2項(xiàng)：1）超級(jí)單體60 dBZ中還有65、70 dBZ的“回波核”，微型超級(jí)單體沒有回波核;2）超級(jí)單體強(qiáng)回波面積≥100 km2，甚至達(dá)到700 km2，降雹時(shí)間較長(zhǎng)，冰雹直徑在2 cm以上，微型超級(jí)單體強(qiáng)回波面積≤100 km2，觀測(cè)到最小18 km2出現(xiàn)小冰雹（冰雹直徑<1 cm）的個(gè)例，降雹時(shí)間不足10 min。

本方法僅使用雷達(dá)拼圖組合反射率因子CR產(chǎn)品，從雷達(dá)拼圖1 100 km×1 100 km經(jīng)緯度網(wǎng)格數(shù)據(jù)著手，研究產(chǎn)生冰雹時(shí)的回波特征。根據(jù)回波四個(gè)顯著特征，選取≥60 dBZ的強(qiáng)回波區(qū)域，通過聚類算法得到強(qiáng)回波中心強(qiáng)度（dBZ）和強(qiáng)回波面積（SEA），通過散點(diǎn)輪廓算法得到強(qiáng)回波梯度（SEG），通過高空風(fēng)尋找特征點(diǎn)來(lái)確定云砧“前伸”弱回波（CAE）比值，最后根據(jù)不同天氣背景、環(huán)境條件和季節(jié)等制定閾值，確定是否識(shí)別為冰雹。

盡管數(shù)據(jù)單一，但雷達(dá)拼圖CR產(chǎn)品四要素揭示了冰雹回波的基本特性，識(shí)別冰雹回波的準(zhǔn)確率較高。我們對(duì)2022—2023年江西6次冰雹過程中20個(gè)冰雹回波個(gè)例的冰雹識(shí)別率達(dá)到100％，即出現(xiàn)冰雹的回波全部被識(shí)別出來(lái)，但也看到空?qǐng)?bào)率為10％～20％，即部分識(shí)別為冰雹的卻沒有出現(xiàn)。空?qǐng)?bào)率高是四要素識(shí)別方法的局限性。事后分析總結(jié)并尋找空?qǐng)?bào)原因，我們發(fā)現(xiàn)有2個(gè)方面需要考慮：考慮不同天氣背景、環(huán)境條件和不同季節(jié)等因素，可以減少部分空?qǐng)?bào);要考慮雷達(dá)拼圖上冰雹的強(qiáng)回波垂直厚度（強(qiáng)回波垂直面積）。冰雹識(shí)別空?qǐng)?bào)都表現(xiàn)為強(qiáng)回波厚度不夠，雖然回波強(qiáng)度達(dá)到標(biāo)準(zhǔn)，但伸展高度不夠，強(qiáng)回波垂直面積達(dá)不到要求。因此，2024年準(zhǔn)備加入“天衍”雷達(dá)拼圖等高顯示產(chǎn)品數(shù)據(jù)，以增加強(qiáng)回波垂直面積的識(shí)別判斷，達(dá)到降低冰雹回波空?qǐng)?bào)率的目的。

冰雹回波的識(shí)別方法很多，如單部雷達(dá)PUP產(chǎn)品、雙偏振雷達(dá)產(chǎn)品等。但這些產(chǎn)品PUP終端的門檻高，不是廣大基層臺(tái)站（縣級(jí)）都能擁有的（江西93個(gè)臺(tái)站只有11個(gè)市級(jí)臺(tái)站擁有）。加上雷達(dá)PUP產(chǎn)品分析方法、基礎(chǔ)知識(shí)、技術(shù)經(jīng)驗(yàn)等方面要求高，也不是廣大基層臺(tái)站在短時(shí)間內(nèi)能掌握的。短臨預(yù)報(bào)講究的是快，能否有一種簡(jiǎn)便、快捷、明確的冰雹識(shí)別產(chǎn)品成為當(dāng)務(wù)之急，廣大基層臺(tái)站預(yù)報(bào)員非常迫切需要這類產(chǎn)品。因此，我們才有了在雷達(dá)拼圖CR產(chǎn)品上識(shí)別冰雹的想法。冰雹識(shí)別平臺(tái)每10 min更新一次，手機(jī)、電腦隨時(shí)隨地都可以瀏覽，沒有門檻，一目了然，滿足了廣大基層臺(tái)站預(yù)報(bào)員的需求。

2.1 強(qiáng)回波值

雷達(dá)反射率因子Z為單位體積內(nèi)所有小球粒子直徑的6次方之和，如式（1）所示，單位為mm6/m3。

Z=∑Ni=1D6i。 （1）

由于粒子直徑的變化對(duì)雷達(dá)反射率因子Z的影響非常敏感，可利用式（2）將雷達(dá)反射率因子Z轉(zhuǎn)化為ZdB，ZdB的變化與粒子直徑的變化幅度趨于一致。

ZdB=10lgZ1 mm6/m3。（2）

雷達(dá)反射率因子Z為單位體積內(nèi)所有小球粒子直徑的6次方之和：1）假設(shè)單位體積（m3）內(nèi)3 mm的小球有一個(gè)，1 mm的小球一共729個(gè)，根據(jù)公式Z=∑Ni=1D6i，Z=（3 mm）6／m3+729×（1 mm）6／m3，得：Z=1 458 mm6／m3;2）將Z=1 458 mm6／m3代入公式2），ZdB=31.6 dBZ。

通過式（1）、（2）可得到單位體積內(nèi)單個(gè)粒子直徑與ZdB的對(duì)應(yīng)關(guān)系（表1）。由表可見，雷達(dá)回波強(qiáng)度的大小與降水粒子的直徑關(guān)系密切。當(dāng)雷達(dá)反射率因子達(dá)到55 dBZ時(shí)，降水粒子直徑為8.25 mm;當(dāng)雷達(dá)反射率因子達(dá)到60 dBZ時(shí)，降水粒子直徑為10 mm。一般冰雹粒子的直徑≥6 mm，因此，當(dāng)回波強(qiáng)度≥50 dBZ時(shí)就有可能出現(xiàn)冰雹。

雷暴大風(fēng)與短時(shí)強(qiáng)降水的多普勒雷達(dá)產(chǎn)品對(duì)比分析結(jié)果表明：雷雨大風(fēng)與短時(shí)強(qiáng)降水過程雷達(dá)反射率因子達(dá)到45～55 dBZ（應(yīng)冬梅等，2007;夏文梅等，2021）。江西冰雹雷達(dá)拼圖回波特征是：回波中心強(qiáng)度達(dá)60 dBZ，超級(jí)單體回波往往發(fā)展為65～70 dBZ（陳鮑發(fā)和馬中元，2019;陳鮑發(fā)等，2022;段和平等，2023;劉懿樞等，2023），一些大冰雹超級(jí)單體回波強(qiáng)度超過70 dBZ，≥60 dBZ的強(qiáng)回波面積≥300 km2（黃中根等，2022）。因此，上述研究結(jié)果表明，江西冰雹CR需要達(dá)到60 dBZ。根據(jù)雷達(dá)反射率與降水粒子的關(guān)系以及實(shí)際觀察結(jié)果，將雷達(dá)反射率因子≥60 dBZ的回波認(rèn)定為強(qiáng)回波。

降水粒子的直徑和粒子數(shù)濃度之間存在一定的關(guān)系，對(duì)回波強(qiáng)度有著顯著的影響。一般來(lái)說(shuō)，降水粒子的直徑越小，粒子數(shù)濃度就越高;粒子直徑越大且數(shù)濃度越高，則回波越強(qiáng)。降水粒子直徑和粒子數(shù)濃度之間的關(guān)系可以通過測(cè)量和監(jiān)測(cè)降水中的顆粒物來(lái)研究。粒子譜分布是指在特定環(huán)境下粒子的直徑或體積在不同尺寸范圍內(nèi)的分布情況。在體積粒子譜中，粒子的直徑和濃度分布關(guān)系到回波的強(qiáng)度分布。雨滴和冰雹數(shù)濃度相差甚大，在降雹過程中，冰雹分鐘數(shù)濃度最大為5 m-3，雨滴分鐘數(shù)濃度最大為1 423 m-3（岳治國(guó)和梁谷，2018）。

2.2 強(qiáng)回波面積

強(qiáng)回波區(qū)域的面積大小對(duì)冰雹生成非常關(guān)鍵，因?yàn)樵谧R(shí)別開始前強(qiáng)回波個(gè)數(shù)是未知的，但強(qiáng)回波都是密集聚在一起，所以強(qiáng)回波面積的識(shí)別可利用密度聚類DBSCAN（Density-Based Spatial Clustering of Applications with Noise）算法。密度聚類將簇定義為密度相連的樣本的最大集合，能夠?qū)⒚芏茸銐蚋叩膮^(qū)域劃分為簇，不需要給定簇?cái)?shù)量，并可在有噪聲的空間數(shù)據(jù)集中發(fā)現(xiàn)任意形狀的簇。樣本分布的緊密程度是一組鄰域參數(shù)ε和MinPts來(lái)描述的，要求聚類空間一定區(qū)域內(nèi)所包含對(duì)象（點(diǎn)或其他空間）的數(shù)目不小于某一給定閾值（馮少榮和肖文俊，2008）。

給定數(shù)據(jù)集D=｛x1，x2，…，xm｝，ε-鄰域（Eps）對(duì)于xj∈D，其ε-鄰域包含樣本集D中與xj的距離不大于ε的樣本。式（3）中，Nε（xj）為子樣本集的個(gè)數(shù)，dist（xi，xj）為xi、xj之間的歐氏距離。

Nε（xj）=｛xi∈D|dist（xi，xj）≤ε｝。 （3）

DBSCAN算法速度快，能夠發(fā)現(xiàn)任意形狀的空間聚類，對(duì)處理的數(shù)據(jù)不需要提前劃分聚類個(gè)數(shù)。但是該算法調(diào)試參數(shù)比較復(fù)雜，主要需要對(duì)距離閾值Eps，鄰域樣本數(shù)閾值MinPts進(jìn)行聯(lián)合調(diào)參，不同的參數(shù)組合對(duì)最后的聚類效果有較大影響。若Eps設(shè)置得非常小，則意味著沒有點(diǎn)是核心樣本，可能會(huì)導(dǎo)致所有點(diǎn)被標(biāo)記為噪聲;若Eps設(shè)置得非常大，則可能會(huì)導(dǎo)致所有的點(diǎn)形成單個(gè)簇。雖然不需要顯示設(shè)置簇的個(gè)數(shù)，但是設(shè)置Eps可以隱式地控制簇的個(gè)數(shù)。

通過對(duì)雷達(dá)拼圖經(jīng)緯度1100格點(diǎn)數(shù)據(jù)進(jìn)行預(yù)處理，選取雷達(dá)回波強(qiáng)度≥60 dBZ的范圍，并記錄≥60 dBZ回波出現(xiàn)的位置和回波強(qiáng)度值，通過DBSCAN算法，可以在≥60 dBZ的回波中提取回波連續(xù)且回波面積≥10 km2的回波區(qū)域。在實(shí)際情況中回波區(qū)域通常是不連續(xù)的，但若是距離非常近的情況，則可將距離相近的回波看作一個(gè)整體，將距離小于3 km的回波歸為同一個(gè)簇，因此對(duì)Eps參數(shù)的設(shè)定以3 km為標(biāo)準(zhǔn)。因?yàn)楫?dāng)≥60 dBZ的回波面積超過10 km2時(shí)極有可能出現(xiàn)冰雹，所以對(duì)閾值Minpts參數(shù)的設(shè)定以10 km2為標(biāo)準(zhǔn)。由圖1可見，圖1a中紫色區(qū)域?yàn)榛夭ā?0 dBZ的范圍，圖1b中綠色點(diǎn)則是自動(dòng)選取出來(lái)的回波區(qū)域。

強(qiáng)回波面積由選定回波區(qū)域中點(diǎn)的總數(shù)確定，回波強(qiáng)度則是選取回波區(qū)域中CR最大值，因?yàn)榛夭◤?qiáng)度的大小和強(qiáng)回波面積對(duì)冰雹過程的影響很大，有些強(qiáng)冰雹發(fā)生時(shí)，回波強(qiáng)度會(huì)超過65 dBZ甚至70 dBZ，所以還需分別記錄選定回波區(qū)域中回波強(qiáng)度≥65 dBZ和回波強(qiáng)度≥70 dBZ的強(qiáng)回波面積。通過記錄不同的回波強(qiáng)度面積大小，可以對(duì)冰雹過程的強(qiáng)弱進(jìn)行區(qū)分。

強(qiáng)回波中心的確定需要考慮回波每個(gè)位置強(qiáng)度的不同，回波最強(qiáng)的位置不一定是幾何中心，所以回波中心不能簡(jiǎn)單地設(shè)定為幾何中心。因此在選定好的回波區(qū)域中，計(jì)算此區(qū)域回波強(qiáng)度的均值，選取此區(qū)域中回波強(qiáng)度大于均值的回波區(qū)域，在這些大于均值的區(qū)域中選取一個(gè)點(diǎn)，此點(diǎn)到所有均值的點(diǎn)距離最短，可以確定此點(diǎn)為回波中心。在回波區(qū)域中還會(huì)遇到多個(gè)強(qiáng)度很大的區(qū)域，而這些區(qū)域的距離又較遠(yuǎn)，這里對(duì)每個(gè)不同回波強(qiáng)度區(qū)域進(jìn)行單獨(dú)計(jì)算，并確定其多個(gè)強(qiáng)回波中心。

2.3 強(qiáng)回波梯度

強(qiáng)回波梯度是指強(qiáng)回波區(qū)域內(nèi)強(qiáng)度變化的大小，以60 dBZ邊緣線至30 dBZ邊緣線的距離為值，這個(gè)距離越短，強(qiáng)回波梯度就越大。強(qiáng)回波梯度的大小和方向可反映回波系統(tǒng)的強(qiáng)度和演變趨勢(shì)。

強(qiáng)回波梯度的計(jì)算方法通常采用數(shù)值微分，通過計(jì)算相鄰兩個(gè)回波強(qiáng)度的差值，再除以它們之間距離的算術(shù)平均值來(lái)估計(jì)強(qiáng)回波梯度。本方法定義強(qiáng)回波梯度為60 dBZ回波區(qū)域邊界點(diǎn)到回波強(qiáng)度30 dBZ邊界之間最短的距離，因此需要知道強(qiáng)回波區(qū)域的外圍輪廓。用離散點(diǎn)提取外圍輪廓時(shí)，利用凸包算法對(duì)平面中的N個(gè)點(diǎn)，找出一個(gè)由其中某些點(diǎn)作為頂點(diǎn)組成的凸多邊形，而這個(gè)凸多邊形恰好能圍住所有的N個(gè)點(diǎn)，凸包雖然和輪廓近似相似，但是凸包算法中凸形曲線總是凸出來(lái)的，在凹進(jìn)去的地方檢測(cè)時(shí)并不精確，會(huì)遇到凸形缺陷。而散點(diǎn)輪廓算法利用一個(gè)半徑為r的圓在點(diǎn)集外圍滾動(dòng)，當(dāng)r的選取適當(dāng)時(shí)，此圓就不能滾動(dòng)到散點(diǎn)的內(nèi)部，此圓經(jīng)過的點(diǎn)就是選取的外圍輪廓點(diǎn)，可以有效避免凸包中凸性缺陷，對(duì)散點(diǎn)邊界進(jìn)行精確提?。▌⒔?，2013）。

利用散點(diǎn)輪廓算法確定回波邊界，散點(diǎn)輪廓算法中的參數(shù)主要是選取圓的半徑，半徑大小對(duì)輪廓的精細(xì)度非常重要，半徑太大時(shí)選取的輪廓外圍太大，而半徑太小時(shí)，有些區(qū)域無(wú)法選擇出來(lái)，因此半徑的選取尤為重要。因?yàn)榇饲癉BSACN算法中選取的MinPts參數(shù)為3 km，當(dāng)圓半徑設(shè)置為1 km時(shí)，直徑是2 km，比較接近MinPts參數(shù)。經(jīng)過驗(yàn)證知道，參數(shù)選取為1 km時(shí)，選取的外圍輪廓效果較好。圖2中白色點(diǎn)是通過散點(diǎn)輪廓算法得到的回波邊界。

由于江西地區(qū)3—5月的強(qiáng)回波梯度通常會(huì)出現(xiàn)在強(qiáng)回波區(qū)域的東南和西南方向，所以選取回波邊界上東南方向和西南方向的點(diǎn)，圖2中綠色點(diǎn)即為選取點(diǎn)。提取選取點(diǎn)周圍回波強(qiáng)度小于30 dBZ的回波區(qū)域，再計(jì)算選取點(diǎn)到小于30 dBZ點(diǎn)之間的距離，選取其中最短的距離，并將此距離確定為回波梯度。圖2中A點(diǎn)是60 dBZ輪廓上的點(diǎn)，B點(diǎn)是60 dBZ到30 dBZ輪廓上的點(diǎn)，A點(diǎn)到B點(diǎn)的距離為強(qiáng)回波梯度。

2.4 云砧回波

云砧回波是以強(qiáng)回波中心點(diǎn)為起點(diǎn)沿著100 hPa和200 hPa高空風(fēng)的平均方向，尋找此方向上強(qiáng)回波中心（≥60 dBZ）到回波強(qiáng)度小于30 dBZ點(diǎn)的距離，與回波強(qiáng)度30 dBZ到小于10 dBZ的距離之間的比值，這個(gè)比值設(shè)定為云砧回波。

通常情況下沿著高空風(fēng)的方向?qū)ふ視r(shí)，回波強(qiáng)度是依次減小的，實(shí)際情況中會(huì)遇到一些回波強(qiáng)度突然變大或突然變小的異常情況，遇到這種情況時(shí)，可以對(duì)風(fēng)向上的回波強(qiáng)度進(jìn)行數(shù)據(jù)平滑處理，在算法的設(shè)計(jì)上使用滑動(dòng)平均方法。根據(jù)式（4），滑動(dòng)平均法把回波強(qiáng)度前后一共2n+1個(gè)值做平均，得到當(dāng)前位置的回波強(qiáng)度，xi為當(dāng)前位置的回強(qiáng)度，分別選取前后兩個(gè)位置點(diǎn)的回波強(qiáng)度做平均，這里選取n為1，求出當(dāng)前的回波強(qiáng)度Z。

Z=∑ni=1（xi-1+xi+1）+xi2n+1。 （4）

高空風(fēng)方向并不完全是云砧回波的方向，可設(shè)置一個(gè)誤差范圍，即在高空風(fēng)上下30°的范圍內(nèi)去選取最長(zhǎng)的云砧。沿著高空風(fēng)的方向如果出現(xiàn)另外一個(gè)強(qiáng)回波區(qū)域，此強(qiáng)回波區(qū)域已完全遮擋了云砧的方向且回波區(qū)域非常大，則在此情況下設(shè)定此時(shí)遇到回波遮擋，可以通過圖3中A點(diǎn)到B點(diǎn)距離與B點(diǎn)到C點(diǎn)距離的比值為云砧回波的大小。

3 “CR-識(shí)別”法的應(yīng)用

我們將“CR-識(shí)別”法應(yīng)用于江西2022年3次冰雹過程（20個(gè)例）和2023年3次冰雹過程（9個(gè)例）。結(jié)果表明，該方法對(duì)這6次冰雹過程都識(shí)別出了冰雹，效果較好，且具有一定的提前量，冰雹識(shí)別率達(dá)到100％，但也有10％～20％的空?qǐng)?bào)率。雷達(dá)拼圖CR產(chǎn)品由于使用簡(jiǎn)單、快捷、方便（手機(jī)、電腦都可以瀏覽），所以對(duì)于廣大基層臺(tái)站（縣級(jí)）預(yù)報(bào)員來(lái)說(shuō)更具有參考和推廣意義。對(duì)于10％～20％的空?qǐng)?bào)率，后期在考慮不同天氣背景、環(huán)境條件、季節(jié)等因素，并融入垂直積分液態(tài)水含量VIL產(chǎn)品和強(qiáng)回波垂直面積要素，可望提高識(shí)別準(zhǔn)確率。造成10％～20％空?qǐng)?bào)率的原因是，雖然四要素都達(dá)到標(biāo)準(zhǔn)，但垂直積分液態(tài)水含量VIL和強(qiáng)回波垂直面積沒有達(dá)到閾值，后期增加這兩個(gè)要素，可以大大減少空?qǐng)?bào)率。這是“CR-識(shí)別”法改進(jìn)的方向。

下面給出6次冰雹過程的識(shí)別效果圖及識(shí)別數(shù)據(jù)表（需要說(shuō)明的是，一次冰雹過程可有多次降雹個(gè)例）。

3.1 2022年3月14日江西及周邊冰雹

2022年3月14日，江西及其周邊出現(xiàn)較大范圍的冰雹天氣，冰雹基本上產(chǎn)生于超級(jí)單體回波。有些是較孤立的經(jīng)典超級(jí)單體，而有些則是微型超級(jí)單體，這兩種超級(jí)單體的區(qū)別簡(jiǎn)單說(shuō)就是強(qiáng)弱差別。一些超級(jí)單體分布在雷暴回波帶上或雷暴回波群中，孤立存在的經(jīng)典超級(jí)單體表現(xiàn)更強(qiáng)（圖4）。圖4是16個(gè)超級(jí)單體產(chǎn)生冰雹的瞬間，可見超級(jí)單體回波強(qiáng)度CR都≥60 dBZ，最強(qiáng)達(dá)到70 dBZ（浙江上甘鄉(xiāng)冰雹），60 dBZ強(qiáng)回波面積都≥25 km2，受天氣系統(tǒng)、環(huán)境條件和季節(jié)影響，3月超級(jí)單體極易產(chǎn)生冰雹。

由2022年3月14日16次冰雹的識(shí)別信息匯總表（表2）可以看到：冰雹發(fā)生時(shí)回波≥60 dBZ的強(qiáng)回波面積基本超過25 km2，只有13：40安義老觀測(cè)站觀測(cè)到冰雹發(fā)生時(shí)的強(qiáng)回波面積為18 km2，回波強(qiáng)度值為61 dBZ，說(shuō)明18 km2也能產(chǎn)生冰雹，但冰雹較小、影響范圍小、時(shí)間短，屬于回波帶上微型超級(jí)單體結(jié)構(gòu)。

由表2可見，在16次冰雹過程中，13次冰雹出現(xiàn)了≥65 dBZ的強(qiáng)回波核，其中上甘鄉(xiāng)≥65 dBZ的強(qiáng)回波面積超過51 km2，并出現(xiàn)5 km2的70 dBZ回波核。其余3次冰雹沒有出現(xiàn)≥65 dBZ的強(qiáng)回波區(qū)域。

是否出現(xiàn)≥65 dBZ的強(qiáng)回波區(qū)域，是區(qū)分經(jīng)典超級(jí)單體與微型超級(jí)單體的判斷條件之一。當(dāng)≥60 dBZ的強(qiáng)回波面積超過40 km2時(shí)，隨之≥65 dBZ的強(qiáng)回波面積增多。因此，60 dBZ強(qiáng)回波面積越大越有利于更強(qiáng)回波的發(fā)展，從而更有利于冰雹的形成與發(fā)展。根據(jù)回波識(shí)別圖的時(shí)間與冰雹出現(xiàn)時(shí)間的比對(duì)，經(jīng)典超級(jí)單體回波基本上能提前10～50 min識(shí)別出冰雹，但微型超級(jí)單體生命史很短，提前量為0 min，也就是說(shuō)無(wú)法提前發(fā)現(xiàn)。當(dāng)然這也與雷達(dá)提掃時(shí)間延遲有關(guān)。

3.2 2022年7月29日婺源冰雹

2022年7月29日17：00（圖5a），婺源出現(xiàn)多個(gè)對(duì)流單體回波發(fā)展，強(qiáng)度只有55 dBZ，并出現(xiàn)回波合并趨勢(shì)。同時(shí)，德興附近識(shí)別出超級(jí)單體冰雹（但沒有冰雹記錄），回波強(qiáng)度為60 dBZ，強(qiáng)回波面積為51 km2。17：10（圖5b），婺源識(shí)別出現(xiàn)冰雹，回波強(qiáng)度為64 dBZ，強(qiáng)回波面積為21 km2，強(qiáng)回波梯度為3.61 km，移動(dòng)十分緩慢，城區(qū)出現(xiàn)大風(fēng)災(zāi)害（婺源縣城景觀橋被吹倒），屬于微型超級(jí)單體回波結(jié)構(gòu)，婺源城內(nèi)出現(xiàn)短時(shí)間≤2 cm的小冰雹。17：20（圖5c），婺源冰雹減弱為55 dBZ，德興冰雹也明顯減弱。這次過程說(shuō)明，微型超級(jí)單體的生命史很短，只維持10 min左右。婺源冰雹識(shí)別信息匯總表如表3所示。

3.3 2022年8月7日石城冰雹

2022年8月7日16：00（圖6a），石城北部識(shí)別出60 dBZ強(qiáng)回波，強(qiáng)回波面積為49 km2，冰雹正在生成中。16：10（圖6b），回波強(qiáng)度增至66 dBZ，強(qiáng)回波面積達(dá)91 km2，強(qiáng)回波梯度為5.1 km，屬于典型超級(jí)單體回波結(jié)構(gòu)，石城北部木蘭鄉(xiāng)和小松鎮(zhèn)出現(xiàn)冰雹;冰雹一直持續(xù)到16：20（圖6c）至16：30（圖6d），回波強(qiáng)度分別達(dá)到67和66 dBZ。16：40（圖6e）和16：50（圖6f），石城冰雹停止，回波逐漸減弱。石城冰雹識(shí)別信息匯總表如表4所示。

3.4 2023年3月22日寧都、石城冰雹

2023年3月22日，江西南部出現(xiàn)一次較強(qiáng)經(jīng)典超級(jí)單體風(fēng)暴，冰雹維持時(shí)間較長(zhǎng)，其中石城出現(xiàn)≥5 cm的大冰雹。14：00（圖7a），中尺度輻合線在贛南萬(wàn)安、興國(guó)、寧都、廣昌一帶，沿輻合線形成一條由多個(gè)對(duì)流單體組成的不連續(xù)回波帶，帶上強(qiáng)回波單體中心強(qiáng)度為55 dBZ。14：30（圖7b），回波帶上強(qiáng)單體不斷發(fā)展，中心強(qiáng)度達(dá)到60 dBZ。15：00（圖7c），回波帶上3個(gè)超級(jí)單體發(fā)展形成，中心強(qiáng)度都達(dá)到60 dBZ，其中最西端興國(guó)超級(jí)單體強(qiáng)回波面積達(dá)到55 km2，是最強(qiáng)的一個(gè)，這時(shí)，由3個(gè)超級(jí)單體組成的中尺度對(duì)流回波群已經(jīng)形成，共一個(gè)卷云罩;此時(shí)，興國(guó)部分鄉(xiāng)鎮(zhèn)出現(xiàn)冰雹。15：30（圖7d），中尺度對(duì)流回波群發(fā)展壯大，興國(guó)超級(jí)單體強(qiáng)回波面積達(dá)到150 km2，回波強(qiáng)度最大達(dá)66 dBZ，出現(xiàn)65 dBZ回波核，強(qiáng)回波梯度為5.39 km，寧都開始出現(xiàn)冰雹。16：00（圖7e），寧都超級(jí)單體發(fā)展最旺盛，回波強(qiáng)度達(dá)到68 dBZ，強(qiáng)回波面積為250 km2，65 dBZ強(qiáng)回波面積為99 km2，強(qiáng)回波梯度為4.12 km，造成寧都、石城出現(xiàn)大冰雹。16：30（圖7f），超級(jí)單體影響石城縣城，回波強(qiáng)度為67 dBZ，強(qiáng)回波面積為85 km2，有65 dBZ回波核，強(qiáng)回波梯度為5.1 km;與前10 min相比，中尺度對(duì)流回波群有所減弱，慢慢移出江西進(jìn)入福建。寧都、石城冰雹識(shí)別信息匯總表如表5所示。

3.5 2023年4月13日峽江、永豐冰雹

2023年4月13日，江西中部出現(xiàn)一次超級(jí)單體冰雹過程，冰雹維持時(shí)間較長(zhǎng)（約2 h）。15：00（圖8a），新余至峽江間有小塊對(duì)流回波發(fā)展，強(qiáng)度為35 dBZ。15：30（圖8b），峽江西部快速生成強(qiáng)回波單體，中心強(qiáng)度達(dá)到60 dBZ，但面積較小（≤10 km2），因此沒有被識(shí)別出來(lái)。16：00（圖8c），回波單體迅速發(fā)展成超級(jí)單體，強(qiáng)回波面積為67 km2，中心強(qiáng)度達(dá)到66 dBZ，出現(xiàn)了65 dBZ回波核，強(qiáng)回波梯度為2.24 km，峽江西部鄉(xiāng)鎮(zhèn)出現(xiàn)冰雹。16：30（圖8d），超級(jí)單體移到峽江縣城，強(qiáng)度為66 dBZ，回波面積為99 km2，存在65 dBZ強(qiáng)回波面積27 km2的回波核，強(qiáng)回波梯度為1.41 km，峽江縣城出現(xiàn)冰雹。17：00（圖8e），超級(jí)單體東移，回波強(qiáng)度為65 dBZ，強(qiáng)回波面積為64 km2，強(qiáng)回波梯度為1 km，有65 dBZ回波核，峽江、永豐部分鄉(xiāng)鎮(zhèn)出現(xiàn)冰雹。17：30（圖8f），永豐超級(jí)單體回波強(qiáng)度為67 dBZ，回波面積為113 km2，強(qiáng)回波梯度為5 km，存在65 dBZ、52 km2回波核，永豐出現(xiàn)冰雹。18：00（圖8g），永豐超級(jí)單體繼續(xù)向東移動(dòng)，回波強(qiáng)度為64 dBZ，回波面積為91 km2，強(qiáng)回波梯度為2.83 km，永豐東部鄉(xiāng)鎮(zhèn)仍有冰雹。18：30（圖8h）和19：00（圖8i），超級(jí)單體開始減弱，地面降雹停止。這次冰雹過程出現(xiàn)了中尺度對(duì)流回波群特征（圖8f）：在同一個(gè)卷云罩下由3個(gè)超級(jí)單體組成的中尺度對(duì)流回波群。峽江永豐冰雹識(shí)別信息匯總表如表6所示。

3.6 2023年11月9日上高冰雹

2023年11月9日，江西上高縣城出現(xiàn)一次小冰雹事件，產(chǎn)生冰雹的超級(jí)單體維持時(shí)間不長(zhǎng)，冰雹直徑≤2 cm。從天氣形勢(shì)來(lái)看，中高緯兩槽一脊，東北冷渦位于黑龍江北部，低渦后部有冷空氣從西北向東南輸送。副高位于華南呈東北-西南走向，受青藏高原繞流影響形成南支槽，有多個(gè)短波系統(tǒng)東移影響長(zhǎng)江中下游地區(qū)。07：30（圖9a），識(shí)別出湖南東部2處冰雹云的強(qiáng)回波朝東-東北方向移動(dòng)。強(qiáng)回波單體中心強(qiáng)度為60 dBZ，但強(qiáng)回波面積≤29 km2。08：00（圖9b），接近萍鄉(xiāng)的強(qiáng)回波單體不斷發(fā)展，中心強(qiáng)度達(dá)到60 dBZ，移入江西境內(nèi)，湖南東北另一處強(qiáng)回波減弱。08：30（圖9c），強(qiáng)回波有所減弱，60 dBZ強(qiáng)回波基本消失。09：00（圖9d），回波重新發(fā)展，其中上高西側(cè)由兩塊對(duì)流小回波合并發(fā)展成超級(jí)單體，60 dBZ強(qiáng)回波中有65 dBZ“回波核”，60 dBZ強(qiáng)回波面積為188 km2，上高縣城開始出現(xiàn)小冰雹，冰雹直徑為1～2 cm。同時(shí)，萬(wàn)載西部回波也在發(fā)展壯大，強(qiáng)回波中心強(qiáng)度為55 dBZ。09：10（圖9e），上高超級(jí)單體的強(qiáng)回波面積減小為111 km2，上高城區(qū)冰雹維持。萬(wàn)載強(qiáng)回波發(fā)展成微型超級(jí)單體，中心強(qiáng)度為60 dBZ，強(qiáng)回波面積為28 km2，地面沒有出現(xiàn)冰雹信息。09：20（圖9f），上高超級(jí)單體明顯減弱，中心強(qiáng)度雖然有60 dBZ，但強(qiáng)回波面積只有11 km2，冰雹停止。同時(shí)，萬(wàn)載微型超級(jí)單體在發(fā)展，中心強(qiáng)度為60 dBZ，強(qiáng)回波面積為66 km2。09：30（圖9g），上高回波減弱為50 dBZ，萬(wàn)載微型超級(jí)單體發(fā)展成超級(jí)單體，中心強(qiáng)度為60 dBZ，強(qiáng)回波面積為114 km2。雖然有65 dBZ“回波核”，但65 dBZ強(qiáng)回波面積很?。ㄖ挥? km2）;雖然60 dBZ強(qiáng)回波面積達(dá)到114 km2，但強(qiáng)回波垂直面積較?。?0 dBZ回波到達(dá)高度較低），萬(wàn)載縣城無(wú)冰雹信息反饋。10：00（圖9h），上高、萬(wàn)載強(qiáng)回波東移減弱。由上高冰雹識(shí)別信息匯總表（表7）可見，09：00，≥60 dBZ的強(qiáng)回波面積為188 km2，≥65 dBZ的強(qiáng)回波面積為30 km2，強(qiáng)回波中心強(qiáng)度為68 dBZ，強(qiáng)回波梯度為1.14 km。

綜上所述，在2022—2023年江西出現(xiàn)的6次冰雹過程中，冰雹識(shí)別率達(dá)到100％，也就是說(shuō)所有冰雹回波都能被識(shí)別出來(lái)，但有10％～20％的冰雹回波識(shí)別出現(xiàn)空?qǐng)?bào)，其原因在于雷達(dá)拼圖四要素識(shí)別冰雹方法存在局限。識(shí)別出來(lái)的冰雹回波，都有可能伴隨短時(shí)強(qiáng)降水、雷暴大風(fēng)天氣，但以冰雹為主。

從回波強(qiáng)度等條件來(lái)看，冰雹要求最高（四要素都具備），雷暴大風(fēng)次之，短時(shí)強(qiáng)降水最低。例如：江西能產(chǎn)生≥30mm／h短時(shí)強(qiáng)降水的回波強(qiáng)度不要求達(dá)到60 dBZ，≥45 dBZ即可，45 dBZ回波面積越大，越有利于短時(shí)強(qiáng)降水的形成。江西雷暴大風(fēng)可以發(fā)生在冰雹回波中，也可以發(fā)生在比冰雹回波弱的較強(qiáng)回波中（55 dBZ）。因此，冰雹的回波強(qiáng)度都比短時(shí)強(qiáng)降水、雷暴大風(fēng)的回波更強(qiáng)。

雷達(dá)拼圖四要素識(shí)別冰雹方法的四要素包含了許多信息，如：強(qiáng)回波梯度的距離，包含了回波墻的緊密程度。當(dāng)然，誤判肯定有，主要是“空?qǐng)?bào)率”較高，初步統(tǒng)計(jì)空?qǐng)?bào)率在10％～20％之間。如果增加雷達(dá)拼圖垂直積分液態(tài)水含量VIL和強(qiáng)回波垂直面積參加計(jì)算，那么空?qǐng)?bào)率將大大降低。

4 結(jié)論與討論

1）在冰雹識(shí)別方法研究中，采用聚類算法和散點(diǎn)輪廓算法可以確定回波強(qiáng)度（回波核）、強(qiáng)回波面積、強(qiáng)回波梯度和云砧回波比值4個(gè)關(guān)鍵識(shí)別特征;通過這些回波特征值，設(shè)定不同天氣背景和環(huán)境條件下的冰雹閾值，來(lái)判斷是否出現(xiàn)冰雹和預(yù)估冰雹直徑的大小。

2）聚類算法和散點(diǎn)輪廓算法計(jì)算效率非常高，對(duì)一些小范圍（≤10 km2）的冰雹也可以準(zhǔn)確地識(shí)別出來(lái)。在江西早春季節(jié)（3—4月）冰雹過程中，該方法可以提前預(yù)判出冰雹可能發(fā)生的前兆，提高預(yù)警時(shí)效;可以對(duì)冰雹進(jìn)行識(shí)別追蹤，提高監(jiān)測(cè)預(yù)警能力。

3）“CR-識(shí)別”法的應(yīng)用實(shí)踐表明：產(chǎn)生冰雹的回波都是超級(jí)單體或微型超級(jí)單體，回波強(qiáng)度CR≥60 dBZ，60 dBZ強(qiáng)回波面積≥100 km2（但早春≥18 km2就有可能出現(xiàn)冰雹），強(qiáng)回波梯度≤8 km，云砧“前伸”回波比值為1∶2～1∶3。

雷達(dá)拼圖CR上識(shí)別冰雹云四參數(shù)的確定，一是根據(jù)多年分析冰雹特征的經(jīng)驗(yàn)值，二是根據(jù)2013—2023年冰雹歷史個(gè)例。強(qiáng)回波垂直面積受到雷達(dá)拼圖沒有CAPPI數(shù)據(jù)的限制，“CR-識(shí)別”法還沒有應(yīng)用它，2024年將使用“天衍”雷達(dá)拼圖中回波垂直結(jié)構(gòu)分析進(jìn)行嘗試。

本文四要素冰雹回波識(shí)別提供兩種產(chǎn)品，一是識(shí)別圖像產(chǎn)品，二是識(shí)別數(shù)據(jù)表格，且每10 min滾動(dòng)更新。后續(xù)研究工作將融合垂直積分液態(tài)水含量VIL和CAPPI數(shù)據(jù)，在垂直強(qiáng)回波面積上進(jìn)行深入探討和實(shí)踐，以進(jìn)一步完善冰雹云識(shí)別方法。

參考文獻(xiàn)（References）

陳鮑發(fā)，馬中元，2019.江西局地冰雹WebGIS雷達(dá)拼圖回波特征分析［J］.氣象與環(huán)境科學(xué)，42（2）：104-114. Chen B F，Ma Z Y，2019.Echo characteristic analysis of WebGIS radar mosaic on hailstone in Jiangxi［J］.Meteor Environ Sci，42（2）：104-114.doi：10.16765／j.cnki.1673-7148.2019.02.015.（in Chinese）.

陳鮑發(fā)，馬中元，王立志，等，2022.江西昌江流域特大暴雨短時(shí)強(qiáng)降水回波特征分析［J］.氣象，48（11）：1418-1427. Chen B F，Ma Z Y，Wang L Z，et al.，2022.Echo characteristics of short-term heavy rainfall due to severe rainstorm in Changjiang River Basin，Jiangxi Province［J］.Meteor Mon，48（11）：1418-1427.doi：10.7519／j.issn.1000-0526.2022.041002.（in Chinese）.

段鶴，嚴(yán)華生，馬學(xué)文，等，2014.滇南冰雹的預(yù)報(bào)預(yù)警方法研究［J］.氣象，40（2）：174-185. Duan H，Yan H S，Ma X W，et al.，2014.Study of the method for hail forecasting and warning in southern Yunnan［J］.Meteor Mon，40（2）：174-185.doi：10.7519／j.issn.1000-0526.2014.02.005.（in Chinese）.

段和平，馬中元，陳鮑發(fā)，等，2023.江西兩次強(qiáng)降雹過程的雷達(dá)回波特征對(duì)比分析［J］.沙漠與綠洲氣象，17（2）：137-145. Duan H P，Ma Z Y，Chen B F，et al.，2023.Comparative analysis of radar echo characteristics of two severe hails in Jiangxi Province［J］.Desert Oasis Meteor，17（2）：137-145.doi：10.12057／j.issn.1002-0799.2023.02.018.（in Chinese）.

馮少榮，肖文俊，2008.DBSCAN聚類算法的研究與改進(jìn)［J］.中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，37（1）：105-111. Feng S R，Xiao W J，2008.An improved DBSCAN clustering algorithm［J］.J China Univ Min Technol，37（1）：105-111.doi：10.3321／j.issn：1000-1964.2008.01.021.（in Chinese）.

黃中根，馬中元，諶蕓，等，2022.江西兩次超級(jí)單體風(fēng)暴雷達(dá)回波特征分析［J］.氣象與環(huán)境科學(xué)，45（5）：39-48. Huang Z G，Ma Z Y，Chen Y，et al.，2022.Analysis of radar echo characteristics of two supercell storms in Jiangxi Province［J］.Meteor Environ Sci，45（5）：39-48.doi：10.16765／j.cnki.1673-7148.2022.05.006.（in Chinese）.

李冰村，唐曉文，何建新，等，2022.基于機(jī)器學(xué)習(xí)的冰雹天氣識(shí)別研究［J］.氣象科學(xué)，42（5）：581-590. Li B C，Tang X W，He J X，et al.，2022.Hail identification based on machine learning methods［J］.J Meteor Sci，42（5）：581-590.doi：10.12306／2021jms.0106.（in Chinese）.

李博勇，胡志群，鄭佳鋒，等，2021.利用貝葉斯方法改進(jìn)華南地區(qū)冰雹識(shí)別效果［J］.熱帶氣象學(xué)報(bào)，37（1）：112-125. Li B Y，Hu Z Q，Zheng J F，et al.，2021.Using Bayesian metho8b333720381c92dcd447ba9b4ba9363ef307b11ed383e98ebb3ceb66ae6d6b2fd to improve hail identification in South China［J］.J Trop Meteor，37（1）：112-125.doi：10.16032／j.issn.1004-4965.2021.011.（in Chinese）.

李德俊，唐仁茂，熊守權(quán)，等，2011.強(qiáng)冰雹和短時(shí)強(qiáng)降水天氣雷達(dá)特征及臨近預(yù)警［J］.氣象，37（4）：474-480. Li D J，Tang R M，Xiong S Q，et al.，2011.Radar features and nowcasting of severe hail and short-time heavy rainfall［J］.Meteor Mon，37（4）：474-480.（in Chinese）.

李金輝，樊鵬，2007.冰雹云提前識(shí)別及預(yù)警的研究［J］.南京氣象學(xué)院學(xué)報(bào)，30（1）：114-119. Li J H，F(xiàn)an P，2007.Investigation on early identification and warning of hail clouds［J］.J Nanjing Inst Meteor，30（1）：114-119.doi：10.13878／j.cnki.dqkxxb.2007.01.016.（in Chinese）.

劉杰，2013.點(diǎn)云邊界點(diǎn)提取算法及其在散點(diǎn)式幾何圖案矢量化方法中的應(yīng)用［D］.無(wú)錫：江南大學(xué). Liu J，2013.Algorithm for extracting boundary points from point clouds and its application in vectorization of scattered geometric patterns［D］.Wuxi：Jiangnan University.（in Chinese）.

劉新偉，蔣盈沙，黃武斌，等，2021.基于雷達(dá)產(chǎn)品和隨機(jī)森林算法的冰雹天氣分類識(shí)別及預(yù)報(bào)［J］.高原氣象，40（4）：898-908. Liu X W，Jiang Y S，Huang W B，et al.，2021.Classified identification and nowcast of hail weather based on radar products and random forest algorithm［J］.Plateau Meteor，40（4）：898-908.doi：10.7522／j.issn.1000-0534.2020.00063.（in Chinese）.

劉懿樞，馬中元，魏鳴，等，2023.一次副熱帶高壓控制下的龍卷雷達(dá)回波特征分析［J］.沙漠與綠洲氣象，17（3）：35-43. Liu Y S，Ma Z Y，Wei M，et al.，2023.Echo characteristics of tornado radar under the control of a subtropical high［J］.Desert Oasis Meteor，17（3）：35-43.（in Chinese）.

路志英，朱俊秀，田碩，等，2015.雷達(dá)回波反射率垂直剖面圖的冰雹識(shí)別方法［J］.天津大學(xué)學(xué)報(bào)（自然科學(xué)與工程技術(shù)版），48（8）：742-749. Lu Z Y，Zhu J X，Tian S，et al.，2015.Hail recognition using radar echo reflectivity cross section［J］.J Tianjin Univ Sci Technol，48（8）：742-749.（in Chinese）.

王瑾，劉黎平，2011.WSR-88D冰雹探測(cè)算法在貴州地區(qū)的評(píng)估檢驗(yàn)［J］.應(yīng)用氣象學(xué)報(bào)，22（1）：96-106. Wang J，Liu L P，2011.The evaluation of WSR-88D hail detection algorithm over Guizhou region［J］.J Appl Meteor Sci，22（1）：96-106.doi：10.3969／j.issn.1001-7313.2011.01.010.（in Chinese）.

吳劍坤，俞小鼎，2009.強(qiáng)冰雹天氣的多普勒天氣雷達(dá)探測(cè)與預(yù)警技術(shù)綜述［J］.干旱氣象，27（3）：197-206. Wu J K，Yu X D，2009.Review of detection and warning methods for sever hail events by Doppler weather radars［J］.J Arid Meteor，27（3）：197-206.doi：10.3969／j.issn.1006-7639.2009.03.001.（in Chinese）.

夏文梅，馬中元，慕瑞琪，等，2021.江西副高邊緣雷暴大風(fēng)雷達(dá)拼圖回波特征分析［J］.大氣科學(xué)學(xué)報(bào)，44（5）：717-726. Xia W M，Ma Z Y，Mu R Q，et al.，2021.Analysis of characteristics of mesoscale thunderstorm gale echo on the edge of subtropical high［J］.Trans Atmos Sci，44（5）：717-726.doi：10.13878／j.cnki.dqkxxb.20201024001.（in Chinese）.

應(yīng)冬梅，許愛華，黃祖輝，2007.江西冰雹、大風(fēng)與短時(shí)強(qiáng)降水的多普勒雷達(dá)產(chǎn)品的對(duì)比分析［J］.氣象，33（3）：48-53，130-131. Ying D M，Xu A H，Huang Z H，2007.A comparative analysis of Doppler radar products in hail，gale and short-time heavy rainfall in Jiangxi［J］.Meteor Mon，33（3）：48-53，130-131.doi：10.3969／j.issn.1000-0526.2007.03.007.（in Chinese）.

岳治國(guó)，梁谷，2018.陜西渭北一次降雹過程的粒子譜特征分析［J］.高原氣象，37（6）：1716-1724. Yue Z G，Liang G，2018.Characteristics of precipitation particles in a hailstorm process in Weibei area of Shaanxi Province［J］.Plateau Meteor，37（6）：1716-1724.doi：10.7522／j.issn.1000-0534.2018.00023.（in Chinese）.

張秉祥，李國(guó)翠，劉黎平，等，2014.基于模糊邏輯的冰雹天氣雷達(dá)識(shí)別算法［J］.應(yīng)用氣象學(xué)報(bào)，25（4）：415-426. Zhang B X，Li G C，Liu L P，et al.，2014.Identification method of hail weather based on fuzzy-logical principle［J］.J Appl Meteor Sci，25（4）：415-426.doi：10.11898／1001-7313.20140404.（in Chinese）.

張文海，李磊，2019.人工智能在冰雹識(shí)別及臨近預(yù)報(bào)中的初步應(yīng)用［J］.氣象學(xué)報(bào)，77（2）：282-291. Zhang W H，Li L，2019.A preliminary application of artificial intelligence on the detection and nowcasting of hail weather［J］.Acta Meteor Sinica，77（2）：282-291.（in Chinese）.

·ARTICLE·

A method for identifying hail clouds based on four elements of radar mosaic CR products

MA Zhongyuan1，2，WANG Jinxin3，4，ZHANG Lincai3，MU Ruiqi2，5，CHEN Baofa6，ZHENG Yuanyuan7，WANG Lizhi8，DUAN Heping9，HUANG Zhikai4，DONG Ling10，ZHANG Qijie3

1Jiangxi Institute of Meteorological Sciences，Nanchang 330096，China;

2Key Laboratory of Transportation Meteorology of China Meteorological Administration，Nanjing 210041，China;

3Nanchang Yunyiran Technology Co.，Ltd.，Nanchang 330000，China;

4School of Information Engineering，Nanchang University of Engineering，Nanchang 330099，China;

5Jiangsu Meteorological Observatory，Nanjing 210041，China;

6Jingdezhen Meteorological Bureau，Jingdezhen 333000，China;

7Nanjing Meteorological Science and Technology Innovation Research Institute，Nanjing 210041，China;

8Institute of Atmospheric Physics，Chinese Academy of Sciences，Beijing 100029，China;

9Jiangxi Climate Center，Nanchang 330096，China;

10Fuzhou Meteorological Bureau，F(xiàn)uzhou 334400，China

Abstract Radar is the most effective tool for detecting hail.In the 1960s and 1970s，the widespread use of rain measurement radar （e.g.，Danka 41 in the UK and domestic radars such as 711，713，etc.） enabled the identification of hail clouds through radar echo features like hook-shaped，finger-shaped，and V-shaped notches.The introduction of Doppler weather radar in the late 1980s and early 1990s provided more accurate data，including radial velocity，for effective hail cloud identification.However，four factors affect the accuracy of hail cloud recognition： 1）single radar limitations such as detection range，distance attenuation，blind spots，and Earth’s curvature;2）the requirement for forecasters to possess high echo analysis skills;3）the variability of classic hail cloud characteristics with orientation，elevation，and distance;4）the limited number of PUP terminals for single radar use，insufficient for county-level forecasters.Radar mosaics can effectively compensate for some of these limitations，especially through web-based radar mosaic CR （Combined Reflectivity） products，which are accessible via computers，tablets，and mobile phones for simple and convenient operation.

The radar mosaic CR product gathers data from multiple radars simultaneously （within ± 3 min）.The blind spot of one radar is covered by another，and the inter-radar distance of 100—150 km is optimal for detection，minimizing issues related to angle blind spots，Earth’s curvature，and distance attenuation.The four key elements for identifying hail clouds on the CR product chart are： 1）echo intensity of 60 dBZ and strong echo core≥65 dBZ，2）strong echo in horizontal and vertical areas，3）strong echo gradient of 30—60 dBZ，and 4）weak echo length formed by cloud anvils.For example，in Jiangxi，hail echo intensity is typically≥60 dBZ，with larger hail having strong echo nuclei above 65 dBZ.A strong echo area of 60 dBZ should be ≥100 km2，although smaller hail may represent a smaller area.The vertical thickness of the strong echo （≥6 km） is also significant，though Jiangxi radar mosaics lack CAPPI products for this measurement.The strong echo gradient indicates hail echo walls，with a steep gradient suggesting a shorter distance.The weak echo formed by cloud anvils reflect the high-altitude winds “pumping” effect.Using the radar mosaic CR product，identifying hail clouds based on these four elements is nearly 100％ successful，with a false report rate below 20％，primarily due to seasonal variability in element thresholds.Adding the vertical area of the strong echo can reduce the false alarm rate.

Automated identification of hail-inducing echoes based on these four elements involves specific algorithms： 1）echo intensity identification through comparison of adjacent points;2）strong echo area identification using clustering and scatter contour algorithms;3）strong echo gradient determination by comparing the distance between 30 dBZ and 60 dBZ;4）cloud anvil echo calculation by measuring the 10 dBZ distance from the 30 dBZ edge along the high-altitude wind direction.Results indicate that hail may occur when the radar mosaic CR is ≥60 dBZ，and the Strong Echo Area （SEA） is ≥100 km2，the Strong Echo Gradient （SEG） is ≤8 km，and the Cloud Anvil Echo （CAE） ratio is between 1∶2 and 1∶3.Most hail in Jiangxi occurs in supercells，though some micro supercells with SEA=18 km2 may also produce hail under suitable conditions.

A method for identifying hail clouds based on these four elements was verified through six hail processes in Jiangxi from 2022 to 2023.The identified hail cloud areas matched actual hail areas，with a false alarm rate of 10％—20％.Future efforts should focus on reducing false alarm rates by incorporating strong echo vertical area and vertical integrated liquid water content elements.This research provides practical experience for simple，fast，and automatic identification of hail weather.

Keywords hail clouds;clustering algorithm; scatter contour algorithm;strong echo gradient;cloud anvil echo

doi：10.13878／j.cnki.dqkxxb.20230801001

（責(zé)任編輯：倪東鴻）