段 煉，郭鵬宇，張 楊

(中國(guó)民用航空飛行學(xué)院 空中交通管理學(xué)院，四川 廣漢 618307)

中國(guó)是一個(gè)強(qiáng)對(duì)流災(zāi)害天氣多發(fā)國(guó)家，強(qiáng)對(duì)流天氣預(yù)報(bào)難度較大，上沖云頂?shù)戎行〕叨葟?qiáng)對(duì)流天氣預(yù)報(bào)難度很大。目前強(qiáng)對(duì)流的監(jiān)測(cè)主要是利用天氣雷達(dá)和氣象衛(wèi)星這兩類遙感資料。天氣雷達(dá)靈敏度高，對(duì)流回波有明顯的特征，但探測(cè)范圍有限，一般有效探測(cè)距離是200～300 km，且容易受測(cè)站周圍障礙物、大氣層結(jié)和天氣類型的影響，不能宏觀探測(cè)對(duì)流系統(tǒng)；氣象衛(wèi)星探測(cè)范圍大，時(shí)空分辨率高，探測(cè)對(duì)流天氣具有直觀、形象和不受地形影響等優(yōu)點(diǎn)，再加上氣象衛(wèi)星技術(shù)的發(fā)展，衛(wèi)星在時(shí)間分辨率和空間分辨率上都得到了顯著的提高[1]，由于上沖云頂生命周期短，相較于雷達(dá)，氣象衛(wèi)星探測(cè)面積大且可以全天候進(jìn)行檢測(cè)的特點(diǎn)可以很好地檢測(cè)和識(shí)別上沖云頂。

美國(guó)氣象學(xué)會(huì)對(duì)上沖云頂?shù)亩x為:“卷云砧上的穹頂狀突起，上升氣流通過它侵入平流層”。上沖云頂比周圍砧云高出2 km左右，在可見光通道觀測(cè)形狀如同花椰菜，在紅外光云圖上表現(xiàn)為低亮溫的成簇圖像。通常直徑小于15 km，當(dāng)上沖云頂上升至較低平流層時(shí)，以7～9 K/km的速率持續(xù)降溫，使它們明顯比周圍云砧溫度更低[2-3];且有明顯突起的紋理。上沖云頂生命期短，有強(qiáng)烈的對(duì)流現(xiàn)象出現(xiàn)，通常發(fā)生在熱帶陸地和海洋地區(qū)上空的強(qiáng)對(duì)流風(fēng)暴中，伴隨災(zāi)害性天氣發(fā)生，如強(qiáng)降水、災(zāi)害性大風(fēng)、大冰雹和龍卷等，是強(qiáng)雷暴天氣的重要指示者，影響飛機(jī)的飛行安全[4]。因此對(duì)于上沖云頂?shù)难芯烤哂惺种匾囊饬x。

1 數(shù)據(jù)

1.1 葵花8氣象衛(wèi)星介紹

數(shù)據(jù)來(lái)自日本氣象廳的Himawari-8(葵花8)對(duì)地靜止氣象衛(wèi)星?？?氣象衛(wèi)星搭載的葵花成像儀(Advanced Image，AHI)可以提供高時(shí)間頻次的區(qū)域和全盤掃描功能，并且單次掃描時(shí)間很短，有利于對(duì)上沖云頂?shù)戎行〕叨葟?qiáng)對(duì)流天氣的觀測(cè)與監(jiān)控。FY-2系列衛(wèi)星紅外通道的空間分辨率為5 km，完成一次全圓盤數(shù)據(jù)掃描需要30 min，相較而言，葵花8號(hào)氣象衛(wèi)星在時(shí)間分辨率和空間分辨率都有顯著的優(yōu)勢(shì)，完成一次全圓盤數(shù)據(jù)掃面僅需要10 min，其中紅外通道空間分辨率為2 km，而可見光通道的空間分辨率為0.5 km，因此，利用葵花8號(hào)氣象衛(wèi)星可以對(duì)上沖云頂進(jìn)行更好的識(shí)別。

1.2 數(shù)據(jù)處理

Himawari-8靜止衛(wèi)星使用的數(shù)據(jù)存儲(chǔ)格式是HDF5，Python具有功能強(qiáng)大、簡(jiǎn)單易學(xué)、運(yùn)用靈活等特點(diǎn)，只需要幾個(gè)函數(shù)命令就可以讀取數(shù)據(jù)，同時(shí)Python有豐富而強(qiáng)大的庫(kù)，利用h5py，即可實(shí)現(xiàn)衛(wèi)星數(shù)據(jù)的讀取。Matplotlib是Python的繪圖庫(kù)，利用此庫(kù)可以將衛(wèi)星數(shù)據(jù)可視化，利用cmaps可以對(duì)不同數(shù)據(jù)進(jìn)行不同顏色的一個(gè)映射，可以利用顏色上的不同區(qū)分識(shí)別物體。本文采用Himawari-8靜止衛(wèi)星的VIS通道、IR3通道、IR8通道對(duì)上沖云頂進(jìn)行識(shí)別，同時(shí)利用ERA-5中的逐小時(shí)降雨數(shù)據(jù)驗(yàn)證上沖云頂?shù)挠绊憽?/p>

2 識(shí)別方法

最初對(duì)于上沖云頂?shù)淖R(shí)別，由Berendes等[5]提出，利用可見光通道和紅外通道對(duì)云紋理進(jìn)行分析，從而探測(cè)深對(duì)流和上沖云頂。雖然該方法在白天能夠很好地執(zhí)行，但是在夜晚或能見度低的條件下則并不能有效地執(zhí)行。此后又有學(xué)者提出下列方法。

2.1 水汽通道和紅外通道的雙通道差值法

在平流層，溫度隨高度的增高而增高，水汽隨著上沖云頂?shù)纳仙龑饬鲙У礁哂谠祈數(shù)奈恢?，所以水汽通道接收到的水汽溫度高于紅外通道的云頂溫度，如果發(fā)生上沖現(xiàn)象，則差值應(yīng)為正數(shù)，Martin等[6]用此來(lái)探測(cè)上沖云頂?shù)拇嬖?。但是這個(gè)方法可能會(huì)受到水汽水平對(duì)流的影響使水汽發(fā)生水平移動(dòng)使識(shí)別出來(lái)的上沖云頂不一定在實(shí)際位置的正上方，并且十分依賴衛(wèi)星的性能，當(dāng)衛(wèi)星性能過低時(shí)不能有效識(shí)別出來(lái)同時(shí)非常冷的卷云砧邊緣可能被錯(cuò)誤地分類為上沖云頂。

2.2 IRW-texture法

Kristopher等[7]利用過沖云頂與砧狀云平均亮度溫度差值判斷上沖云頂?shù)臈l件，提出的IRW-texture(通過小群的紅外窗非常冷亮溫區(qū)域的存在推測(cè)上沖云頂)法。由統(tǒng)計(jì)1979—2003年美國(guó)的平均對(duì)流層頂溫度[7]，使用-58 ℃或當(dāng)?shù)貙?duì)流層頂溫度兩者較小值作為基準(zhǔn)點(diǎn)。亮溫低于-48 ℃或者當(dāng)?shù)貙?duì)流層頂溫度加10 ℃兩者較小值作為上沖云頂?shù)暮蜓a(bǔ)區(qū)域，計(jì)算候補(bǔ)像元的平均亮溫值，根據(jù)上沖云頂?shù)闹睆叫∮?5 km，若最冷點(diǎn)低于-58 ℃，直徑小于15 km，且低于候補(bǔ)區(qū)域的平均亮溫值6.5 ℃則為識(shí)別出來(lái)的真上沖云頂。但是該閾值的選取并不一定適用于所有情況，此外有的上沖云頂會(huì)出現(xiàn)漏判的現(xiàn)象。

2.3 局部最小值法

本文選取局部最小值法[8]，結(jié)合上文的IRW-texture法和水汽通道和紅外通道的雙通道差值法，以11.2 μm通道識(shí)別的亮溫值劃分成小塊區(qū)域，以測(cè)算到最小值為上沖云頂?shù)膮⒖级皇且詫?duì)流云頂溫度，這樣可以根據(jù)具體情況來(lái)選取閾值。方法如下：將11.2 μm通道利用Python畫出紅外通道衛(wèi)星云圖，將圖像劃分為若干個(gè)小區(qū)域，計(jì)算出每一個(gè)小區(qū)域的最小亮溫值，如果一個(gè)像元它的11.2 μm亮溫小于亮溫最小值加20 ℃，及11.2 μm亮溫小于-48 ℃，則被劃分為候補(bǔ)區(qū)域，最后利用水汽通道和紅外通道的雙通道差值法進(jìn)行篩選，當(dāng)差值大于1 ℃且直徑小于15 km的區(qū)域則被認(rèn)為是真上沖云頂。

3 案例分析

受切變?cè)葡涤绊懀?019年6月21日,湖北、安徽、重慶等地出現(xiàn)強(qiáng)對(duì)流天氣，08:00 時(shí)，低渦切變線云系呈東西帶狀分布，其中湖北東部、安徽南部、重慶東南部有強(qiáng)對(duì)流云團(tuán)發(fā)展，云頂亮溫較低，最低亮溫達(dá)到-65 ℃以下,上沖云頂明顯，對(duì)流發(fā)展較強(qiáng)。圖1的高層散度圖顯示，從西南地區(qū)東部至江漢、江淮、江蘇北部等地維持一個(gè)高層輻散區(qū)，有利于上述地區(qū)出現(xiàn)上升運(yùn)動(dòng)。該對(duì)流過程給湖北、安徽、重慶、浙江等地帶來(lái)降雨。

圖1 高層散度圖

本文利用Python將Himawari-8靜止氣象衛(wèi)星數(shù)據(jù)可視化后，利用上文所介紹的水汽通道和紅外通道雙通道差值法、IRW-texture法、局部最小值法對(duì)上沖云頂進(jìn)行識(shí)別。圖2(a)～圖2(f)分別為7：10—8：00可見光云圖，右邊圖例為圖像反射率。在圖中可以看見，在重慶中部以及湖北的中東部存在若干個(gè)濃白色成簇的云團(tuán)同時(shí)伴有陰影，相較于周圍云團(tuán)亮度更高，反射率更大，是上沖云頂因?yàn)楦哂谠茍F(tuán)在可見光照射下所存在的特征。

圖2 7:10—8:00每10 min可見光云圖

首先運(yùn)用雙通道插值法將IR3(6.9 μm)通道減去IR8(11.2 μm)通道，當(dāng)差值大于1 ℃的時(shí)候就是需要識(shí)別的上沖云頂位置。圖3中，紅點(diǎn)為識(shí)別出的上沖云頂，但是黑色箭頭所指位置，在可見光云圖上并沒有上沖云頂突起的特征，因此該方法將很多溫度較低云砧識(shí)別為上沖云頂。

圖3 IR3-IR8雙通道差值法

利用IRW-texture法，根據(jù)郭佳賓等[9]統(tǒng)計(jì)的中國(guó)地區(qū)對(duì)流層頂參數(shù)變化，本文選取的目標(biāo)區(qū)域30°N附近對(duì)流層頂溫度為-68 ℃，如圖4所示，其中黃色區(qū)域?yàn)楹蜻x區(qū)域，紅點(diǎn)為識(shí)別的上沖云頂，我們能看到在可見光云圖的許多有上沖云頂特征的突起并成功識(shí)別，如圖4中黑色箭頭所指，存在很多漏判的現(xiàn)象。

圖4 IRW-texture法

最后利用局部最小值法，將圖像劃分為若干份小區(qū)域，首先圖5中藍(lán)點(diǎn)為識(shí)別出的局部亮溫最低點(diǎn)，以此為基準(zhǔn)點(diǎn)，黃色區(qū)域?yàn)楹Y選出來(lái)的候選區(qū)域，最后利用雙通道差值，用IR3(6.9 μm)通道減去IR8(11.2 μm)通道，以1 ℃為閾值進(jìn)行篩選，落在篩選區(qū)域內(nèi)的則被判定為真上沖云頂區(qū)域。

圖5 局部最小值法

對(duì)比3種方法發(fā)現(xiàn)，水汽通道和紅外通道差值法識(shí)別的上沖云頂最多，IRW-texture 法較少，局地最小法介于中間，三者識(shí)別的位置，基本上而言是一致的，識(shí)別的位置落于對(duì)流較強(qiáng)的云系中，根據(jù)ERA5數(shù)據(jù)得到7：00—8：00一小時(shí)降雨量(圖6)可以看到，在上沖云頂出現(xiàn)的位置大都出現(xiàn)了較強(qiáng)的降雨?？梢姽庠茍D上較白且紅外線云圖上亮溫較低的地方，使用水汽通道減去紅外通道差值法和IRW-texture法會(huì)出現(xiàn)誤判，IRW-texture法會(huì)將一些上沖云頂漏判，這可能是由于閾值以及對(duì)流層頂溫度的選取有關(guān)。而水汽通道減去紅外通道差值法識(shí)別位置比較零亂，且會(huì)將一些中尺度的對(duì)流系統(tǒng)以及比較冷的云砧誤判成上沖云頂，造成的原因可能是對(duì)流持續(xù)的時(shí)間比較長(zhǎng)，從而造成大多數(shù)水汽被帶到對(duì)流層之上。局部最小值法可以將兩者進(jìn)行結(jié)合，既利用IRW-texture法上沖云頂高于對(duì)流層頂?shù)奶攸c(diǎn)，同時(shí)由結(jié)合上沖云頂將水汽帶到更高處的情況，以局部溫度最小值作為判定依據(jù)，通過最小值來(lái)確定閾值避免了認(rèn)為選定閾值所帶來(lái)的誤差，有效地識(shí)別上沖云頂。

圖6 7:00—8:00一小時(shí)降雨量

4 結(jié)語(yǔ)

利用水汽通道和紅外通道差值法、IRW-texture法、局地最小法來(lái)識(shí)別上沖云頂，這3種方法所使用的程序應(yīng)用靈活方便，運(yùn)算時(shí)間較短，能夠顯示出上沖云頂與對(duì)流之間的關(guān)系。其中局部最小值法相較于其他兩種方法能夠更準(zhǔn)確地識(shí)別上沖云頂。在未來(lái)的工作中將收集足夠的上沖云頂案例，對(duì)3種方法進(jìn)行統(tǒng)計(jì)分析，得到三者偵測(cè)的準(zhǔn)確率、虛報(bào)率及漏報(bào)率，找到缺陷，提高方法的可靠性，能更好地對(duì)上沖云頂?shù)葟?qiáng)對(duì)流天氣進(jìn)行監(jiān)測(cè)。