基于Python的西雙版納州人工增雨效果評估平臺的設(shè)計與實現(xiàn)

摘 要：為實現(xiàn)人工增雨地面作業(yè)效果評估結(jié)果的客觀快速上報，提高西雙版納州人工增雨作業(yè)效果評估能力，根據(jù)中國氣象局人工影響天氣中心編制的《人工增雨作業(yè)效果檢驗技術(shù)指南》（2016版），基于人工增雨地面作業(yè)、智能網(wǎng)格實況、天氣雷達、T-lnP等數(shù)據(jù)，利用Python語言將效果評估算法方案轉(zhuǎn)化成自動化生成的西雙版納州人工增雨作業(yè)效果評估平臺。

關(guān)鍵詞：西雙版納；人工增雨；地面作業(yè)；效果評估；評估平臺

中圖分類號：P481 文獻標志碼：B 文章編號：2095–3305（2024）05–0-03

人工增雨作業(yè)效果評估既是人工影響天氣科研和業(yè)務(wù)工作的重要組成部分，也是人工影響天氣最關(guān)鍵的核心問題。人工增雨效果的定量評估在為農(nóng)服務(wù)、改善生態(tài)環(huán)境和防災(zāi)減災(zāi)中具有重要意義，科學(xué)、準確地評估作業(yè)效果是政府和公眾對人工影響天氣活動支持和投入的依據(jù)[1]。當前，人工增雨效果評估方法基本上分為物理分析方法、統(tǒng)計檢驗方法和數(shù)值模擬方法。采用區(qū)域歷史回歸統(tǒng)計檢驗技術(shù)確定催化影響區(qū)，根據(jù)對比區(qū)確定原則選擇影響區(qū)的匹配對比區(qū)，計算相對增雨和絕對增雨，從而獲得人影作業(yè)的增雨效果[2-7]。

針對西雙版納州地面作業(yè)（火箭、煙爐作業(yè)），將效果評估算法方案轉(zhuǎn)化成自動化生成的西雙版納州人工增雨作業(yè)效果評估平臺（以下簡稱平臺）以投入業(yè)務(wù)使用，提高工作效率。在自動生成人工增雨作業(yè)效果評估報前，西雙版納人工增雨效果評估工作一般采用手工計算的方式，具有一定人為主觀性且還會影響作業(yè)效果上報速度。該自動平臺的應(yīng)用有利于改變?nèi)斯ぴ鲇晷Чu估手工計算的現(xiàn)狀，可在每次作業(yè)完成后較便捷、快速地計算出人工增雨效果并存檔，且可統(tǒng)計計算多次增雨效果，從而實現(xiàn)效果評估結(jié)果的客觀快速上報，提高西雙版納州人工增雨作業(yè)效果評估能力。

1 數(shù)據(jù)準備和算法方案

1.1 數(shù)據(jù)準備

（1）實際申請空域增雨作業(yè)信息，包括作業(yè)開始時間、經(jīng)緯度、彈藥類型和用彈量信息。將收集到的火箭增雨站點申請空域作業(yè)信息作為系統(tǒng)的基礎(chǔ)數(shù)據(jù)庫，整合至電腦文件。

（2）使用近10年區(qū)域站實況數(shù)據(jù)、普洱多普勒天氣雷達數(shù)據(jù)（Z9879）、思茅站（56964）高空實況資料及人工增雨作業(yè)最近氣象站點獲取的云南省天擎系統(tǒng)氣象站點的雨量信息等多源觀測數(shù)據(jù)。

1.2 算法方案

效果評估算法方案采用區(qū)域歷史回歸統(tǒng)計檢驗方法，基于日雨量資料分析該作業(yè)過程增雨效果，計算相對增雨和絕對增雨，從而獲得人影作業(yè)的增雨效果。算法方案根據(jù)中國氣象局人工影響天氣中心編制的《人工增雨作業(yè)效果檢驗技術(shù)指南》（2016版）。

（1）利用西雙版納區(qū)域自動站的小時降水量資料對統(tǒng)計檢驗分析該地面增雨過程作業(yè)效果，所用小時降水量數(shù)據(jù)集的時間序列是2014—2023年。將近10年人工增雨作業(yè)時段的降水數(shù)據(jù)剔除，計算出一個多年來小時降水值資料，使用Cressman插值法將其進行網(wǎng)格化插值處理，作為歷史對比降水量。

（2）借助距離較近的思茅探空站（56964）各時次探空數(shù)據(jù)及普洱多普勒天氣雷達（ZZ9879）雷達回波組合反射率對該增雨作業(yè)過程進行合理性分析。

（3）基于地面降水量資料對該作業(yè)過程增雨效果進行統(tǒng)計檢驗。

（4）根據(jù)實際作業(yè)情況，結(jié)合作業(yè)高度高空風(fēng)的風(fēng)向風(fēng)速和雷達回波移向和移速，大致確定作業(yè)影響區(qū)即目標區(qū)站點。采用半定量拋物線方法，來測算人工作業(yè)影響區(qū)范圍：

y2=p[x+（a2/p）]（1）

x=u×t（2）

式（1），式（2）中，x軸為水平風(fēng)向，y軸垂直于水平風(fēng)向，u為水平風(fēng)速，a為火箭水平射距，p為拋物線形狀的系數(shù)，采用的p系數(shù)值為1.5，坐標原點取作業(yè)點，則拋物線所包圍的區(qū)域就作為人工作業(yè)影響區(qū)。

（5）作業(yè)影響區(qū)站點確定后，根據(jù)對比區(qū)的選擇原則，選擇對比區(qū)站點，求得歷史期和作業(yè)期目標區(qū)和對比區(qū)的區(qū)域降水量、區(qū)域平均日雨量的相關(guān)系數(shù)。

（6）利用區(qū)域歷史回歸統(tǒng)計檢驗方法進行統(tǒng)計分析時要求統(tǒng)計變量滿足正態(tài)分布。

（7）將目標區(qū)作業(yè)期實測降水量與目標區(qū)作業(yè)期自然降水量的估計值進行比較確定絕對增雨量和相對增雨率。

（8）采用t-檢驗法對上述求得的絕對增雨量和相對增雨率進行統(tǒng)計顯著性檢驗。

2 評估報告設(shè)計與實現(xiàn)

根據(jù)西雙版納人工增雨作業(yè)的特點，充分利用氣象資源，建立一個操作簡單、運行穩(wěn)定、功能實用，基于Python實現(xiàn)西雙版納州人工增雨地面作業(yè)效果評估報告的自動生成，評估報告的設(shè)計與實現(xiàn)流程圖見圖1。

3 評估報告自動生成的實現(xiàn)

服務(wù)產(chǎn)品所用數(shù)據(jù)均來自氣象大數(shù)據(jù)云平臺·天擎，使用Pathon進行編譯，下載作業(yè)期間的實況數(shù)據(jù)，繪制雷達組合反射率圖、作業(yè)前T-lnP圖及作業(yè)后三小時降水實況圖，根據(jù)實況及作業(yè)信息計算作業(yè)面積、增雨量、增雨率。其計算部分代碼如下：

def XYrange（u，d，a，p，lat0，lon0）：

u為水平風(fēng)速；d為風(fēng)；a為火箭水平射距；火箭彈5 km；煙爐2 km；p為拋物線形狀系數(shù)，取1.5；lat0、lon0為作業(yè)點經(jīng)緯度；XXX、YYY、S分別為作業(yè)面積經(jīng)度、緯度、作業(yè)面積。

x=symbols（‘x’）

y=（p*（x+a**2/p））**0.5

A=integrate（y，（x，0，x2））

S=A*2

X1=np.linspace（0，u*3*60*60/1 000，100）

Y1=（p*（X1+a**2/p））**0.5

Y2=-Y1

Y2=[u*3*60*60/1 000]*100

Y3=np.linspace（Y1[-1]， Y2[-1]，100）

X3=[0]*20

Y4=np.linspace（Y2[0]，Y1[0]，20）

X=list（X1）+list（X2）+list（X1[：：-1]）+list（X3）

Y=list（Y1）+list（Y3）+list（Y2[：：-1]）+list（Y4）

angle=d-90

if anglegt;=0：

XX=[a+b for a，b in zip（[item*math.cos（angle） for item in X]，[item*math.sin（angle） for item in Y]）]

YY=[a-b for a，b in zip（[item*math.cos（angle） for item in Y]，[item*math.sin（angle） for item in X]）]

else：

angle=abs（angle）

XX=[a-b for a，b in zip（[item*math.cos（angle） for item in X]，[item*math.sin（angle） for item in Y]）]

YY=[a+b for a，b in zip（[item*math.sin（angle） for item in Y]，[item*math.cos（angle） for item in X]）]

XXX=[item/111.1+lon0 for item in XX]

YYY=[item/111.1+lat0 for item in YY]

return XXX，YYY，S

4 評估報告的主要內(nèi)容

針對地面作業(yè)，利用Python語言將效果評估算法方案轉(zhuǎn)化為自動計算并輸出人工增雨作業(yè)效果評估報告，報告使用Word文檔的方式，自動給出作業(yè)效果檢驗報告的文件格式，以2024年3月21—22日西雙版納人工增雨作業(yè)為例，西雙版納州人工增雨作業(yè)效果評估報告內(nèi)容包括：

（1）增雨作業(yè)過程效果評估詳細信息表包括單次作業(yè)的作業(yè)時間、地點、經(jīng)緯度、作業(yè)用彈量、作業(yè)增雨率、增雨面積和增雨量報表（表1）；分縣區(qū)的作業(yè)次數(shù)、用彈量、增雨率、增雨面積和增雨量報表。

（2）增雨作業(yè)過程進行合理性分析：在地圖上標注作業(yè)點位置，應(yīng)用T-lnP探空數(shù)據(jù)和多普勒天氣雷達組合反射率，以圖文形式給出作業(yè)云和對比云變化特征的對比結(jié)果，分析作業(yè)過程是否充分（圖2）。

（3）增雨作業(yè)前后降雨量對比圖：生成作業(yè)前和作業(yè)后3 h降水實況圖，可以更清晰、更簡便、更直觀地得出增雨作業(yè)效果的相關(guān)判斷依據(jù)，且直觀、簡便，使增雨作業(yè)效果評估工作方便、快捷（圖3）。

5 結(jié)論

（1）基于Python語言的西雙版納州人工增雨地面作業(yè)效果評估平臺，集數(shù)據(jù)讀取、處理、展示為一體效果評價自動生成平臺，實現(xiàn)作業(yè)效果評價方法和計算過程的統(tǒng)一。

（2）平臺一鍵式生成，操作過程簡單，為西雙版納各縣市人工增雨作業(yè)提供了一套操作簡便、實用性強的人工增雨效果評價平臺，提高工作人員的工作效率。

（3）人工增雨效果報告不僅給出了作業(yè)前后的雷達回波參數(shù)及作業(yè)前后降水實況變化圖，對作業(yè)效果進行數(shù)值統(tǒng)計檢驗，為提升作業(yè)的準確性、科學(xué)性提供了數(shù)據(jù)。

（4）平臺投入使用后，不僅充實了西雙版納人工影響天氣業(yè)務(wù)平臺和業(yè)務(wù)工作流程，提供了快速計算增雨效果的方法，提高了工作效率，還規(guī)范了西雙版納人工增雨效果評估工作，為開展人工增雨作業(yè)效果評估工作提供了更加科學(xué)的依據(jù)。

（5）該自動生成平臺的應(yīng)用改變?nèi)斯ぴ鲇晷Чu估手工計算的現(xiàn)狀，可以在每次作業(yè)完成后較便捷、快速地計算出人工增雨效果并存檔，且可對多次增雨效果作統(tǒng)計分析，從而實現(xiàn)效果評估結(jié)果的客觀快速上報，提高西雙版納人工增雨作業(yè)效果評估能力。

參考文獻

[1] 王飛，李集明，姚展予，等.我國人工增雨作業(yè)效果定量評估研究綜述[J].氣象，2022，48（8）：945-962.

[2] 汪玲，韋增岸，程鵬，等.湖南人工增雨作業(yè)效果統(tǒng)計檢驗與分析[J].氣象研究與應(yīng)用，2019，40（3）：85-89.

[3] 翟晴飛，敖雪，袁健，等.基于區(qū)域歷史回歸法的遼寧地區(qū)一次人工增雨作業(yè)效果檢驗[J].氣象與環(huán)境學(xué)報，2017，33 （6）：96-104.

[4] 尹先龍，王鑫凱.臨海市人影作業(yè)的實施和效果分析[J].科技通報，2020，36（7）：35-38.

[5] 李斌，鄭博華，史蓮梅，等.利用區(qū)域回歸分析法對阿克蘇地區(qū)人工防雹作業(yè)效果再分析[J].新疆農(nóng)業(yè)科學(xué)，2017，54 （9）：1756-1764.

[6] 高建飛，吳進忠，熊凱，等.務(wù)川一次人工增雨作業(yè)效果分析評估[J].科技創(chuàng)新與應(yīng)用，2023，13（19）：101-104，109.

[7] 王霄，孫建印，趙宇，等.氣象雷達在人工增雨效果物理檢驗中的應(yīng)用[J].河南科學(xué)，2022，40（11）：1842-1849.

作者簡介：洪群艷（1987—），女，云南大理人，工程師，研究方向為氣象科技服務(wù)。#通信作者：高婷婷（1994—），女，云南楚雄人，研究方向為氣象科技服務(wù)，E-mail：599394620@qq.com。