有效降水指數在暴雨洪澇監測和評估中的應用*

秦鵬程，劉 敏，李 蘭

（武漢區域氣候中心，武漢 430074）

?

有效降水指數在暴雨洪澇監測和評估中的應用*

秦鵬程，劉 敏，李 蘭

（武漢區域氣候中心，武漢 430074）

摘要：科學有效的監測和評估是防范和減輕暴雨洪澇災害的重要基礎。基于有效降水指數（EP）構建單站和區域暴雨洪澇監測、評估指標，利用1961-2014年湖北省76站逐日氣象觀測資料及相關災情資料，確定降水衰減參數及致澇閾值，在此基礎上分析EP指數在歷史暴雨洪澇評估及實時暴雨洪澇過程監測中的應用效果。結果表明：經參數率定后的EP指數對農作物洪澇受災面積的解釋方差達78.1%，對年際間暴雨洪澇強度差異反應敏感，能識別歷史典型大澇年和嚴重洪澇年，在2014年實時暴雨洪澇過程監測中能直觀診斷出一般性暴雨洪澇的起止時間和過程動態變化，但對局地性和間歇性發生的暴雨洪澇過程刻畫不足。創建EP指數所需數據資料少、計算簡便，可用于洪澇災害歷史排位、年景評價、災情預評估、風險區劃以及作物產量建模等。

關鍵詞：有效降水指數；參數率定；暴雨洪澇；應用檢驗

秦鵬程，劉敏，李蘭.有效降水指數在暴雨洪澇監測和評估中的應用[J].中國農業氣象，2016，37（1）:84-90

洪澇災害是全球最為頻繁的自然災害之一，其造成的生命財產和經濟損失居各類自然災害的前列*[1]。中國季風氣候顯著，降水集中，且地形復雜，是全球洪澇災害頻發和重發區域之一[2-4]。因此，對暴雨洪澇的監測預警和評估具有重要的現實意義。

國內外針對洪澇災害的監測評估方法和案例研究較多，其復雜程度各異，對資料的要求也不同，如依據氣象指標、地形地貌、水文模型、衛星遙感及災情統計資料等開展的洪澇災害監測、預警、災情評估、風險區劃等[1]。目前，對洪澇災害評估的研究進展迅速，從微觀、簡單系統到中觀復雜及宏觀巨系統，洪澇災情評估范圍不斷拓展，評估手段不斷完善，評估結果的客觀性和科學性不斷提高。然而，隨著專業化和多元化評估的不斷深入，其對信息獲取、軟硬件設施及計算時間的要求也越來越高[5]，相反，其普適性和應用范圍卻越來越有限。

盡管形成暴雨洪澇的災害系統異常復雜，但其致災因子主要是過強或過于集中的降水[6]，因此，以降水為主導因子建立暴雨洪澇危險性的評估指標，在洪澇災害監測預警及災前和災中快速評估中仍有良好的應用前景。溫泉沛等[7]基于10個降水因子開展了中國中東部地區暴雨氣候及其農業災情風險評估，郭廣芬等[6]基于過程雨量建立了湖北省暴雨洪澇等級劃分模型，其不足之處在于對地表水分收支考慮不足，同時也不便于確定暴雨洪澇過程的起止時間。Byun等[8]提出了有效降水的概念，即降水經蒸發、滲漏、徑流等物理過程后的剩余量，Lu[9]通過推導證明有效降水隨著時間的推移呈指數衰減，并示范了基于有效降水指數的旱澇監測應用，Deo等[10]基于有效降水指數建立了洪澇監測評估指，并在澳大利亞進行了應用檢驗，張國平[11]基于有效雨量建立了滑坡泥石流災害預測模型，趙一磊等[12]基于有效降水指數建立了干旱監測指標，但在暴雨洪澇監測評估方面的應用還未見報道。本文以湖北省為例，利用歷史災情資料對有效降水指數進行參數率定，并對其在暴雨洪澇監測和評估中應用的可行性進行分析，以期為暴雨洪澇災害監測預警服務和評估業務提供依據。

1 資料與方法

1.1 資料來源

氣象資料為1961年1月1日-2014年12月31日湖北省76個氣象站的逐日觀測資料，由湖北省氣象信息與技術保障中心提供，數據均經過質量檢驗，站點分布及高程信息見圖1。1961-2014年湖北省洪澇災情資料來自歷年《湖北省農村統計年鑒》、《中國氣象災害大典（湖北卷）》及湖北省民政廳災情快報。

圖1 研究區域高程及站點(·）分布Fig. 1 Distribution of digital elevation and 76 meteorological stations(·) in study area

1.2 利用有效降水指數監測與評估暴雨洪澇過程的方法

1.2.1 有效降水指數

本文采用文獻[9]定義的有效降水指數，即

式中，EP為有效降水指數，a為降水衰減參數，取值0～1，P（t）為t時刻降水量，t為距離檢測日的日數，t＝0表示檢測當日，t＝1表示前一日，以此類推，N為前期降水對當前影響的有效時長，理論上可以取無窮大（即考慮距當前無窮日前降水的影響），但由于隨著距離當前日數的增加，降水的權重衰減迅速，距當前14d時降水權重已不足5%（a≤0.8時），因此，在暴雨洪澇的監測評估中取14d既可滿足監測需求，同時也降低了資料收集和計算要求。從式（1）可以看出，EP實際上相當于加權累積降水，因此與降水量具有相同的量綱，EP數值越小，表示有效降水越少，偏旱；EP數值越大，表示有效降水越多，偏澇。

1.2.2 單站暴雨洪澇過程識別與評估

通過對降水資料序列進行滾動計算，建立EP指數的時間序列，根據暴雨洪澇致災閾值，當EP指數超過暴雨洪澇致災閾值時確定發生一次暴雨洪澇過程，過程起始日期為第1天EP指數大于致災閾值的日期，結束日期為最后1次EP指數大于致災閾值的日期。過程強度以過程內的EP指數累積值表示，即

式中，S（t）為t時刻暴雨洪澇強度，EP（i）為自過程起始日開始第i日的有效降水，EPthr為暴雨洪澇致災閾值。

1.2.3 區域暴雨洪澇過程識別與評估

區域性暴雨洪澇的確定通常要求發生暴雨洪澇的站點數達到一定數量或百分比，如福建省區域性暴雨的界定要求發生暴雨的站點數不少于3個（或至少5%）[13]，湖北省區域性暴雨天氣過程判定的最低站點數量為7～10個，約占總站點數的10%[14]。為了綜合考慮暴雨洪澇過程的發生范圍和強度，首先對單站建立EP指數序列，依據單站暴雨洪澇致災閾值，確定單站的暴雨洪澇過程，然后對每日發生暴雨洪澇的單站EP指數進行累加，作為區域暴雨洪澇的監測指標，為了避免監測站數不同及過程內數據缺失造成的誤差，將累加值除以總站數，參照單站暴雨洪澇過程的識別和評估方法確定起止時間和強度，其中區域暴雨洪澇過程的閾值為單站暴雨洪澇致災閾值與區域性過程判定最低站數（百分比）的乘積。

區域暴雨洪澇監測指數用公式表示為

其中

式中，RFI（t）為t時刻區域暴雨洪澇監測指數，n為區域內觀測站的總個數，EPthr為暴雨洪澇致災閾值。

區域暴雨洪澇的綜合強度RSI表示為式中，Si（t）為第i個站點t時刻的暴雨洪澇強度。區域暴雨洪澇站次比RPI定義為當日發生暴雨洪澇的站數（nflood）與區域內總站數（n）的比值，即

1.3 參數率定及洪澇等級劃分

降水衰減參數和洪澇致災閾值通過遺傳算法優化確定，遺傳算法是借鑒生物界自然選擇思想和自然遺傳機制的一種全局優化算法，在水文模型參數率定中具有廣泛應用[15]。以湖北省歷年農作物洪澇受災面積與該年暴雨洪澇過程累積強度的相關系數作為遺傳算法的目標函數，在預先設定的參數取值空間，當降水衰減參數和洪澇致災閾值的組合使洪澇受災面積與過程累積強度的相關系數達到最大時即為最優參數。優化過程通過R語言genalg包實現。暴雨洪澇等級劃分為一般性洪澇、嚴重洪澇和特大洪澇3級。一般性洪澇等級劃分閾值直接采用上述優化后的致災閾值，嚴重洪澇和特大洪澇等級通過對76站1961-2014年歷次暴雨洪澇過程的最大EP值進行概率分布擬合計算重現期，分別以1a一遇和5a一遇對應的EP指數值作為等級劃分依據。概率分布擬合基于超定量法的廣義帕累托分布[16]，通過R語言extRemes包實現。

2 結果與分析

2.1 降水衰減權重及洪澇等級劃分閾值的確定

從式（1）可以看出，如果降水衰減參數a取值0.5，則降水量減弱50%需要1d，取值0.7則需3d，取值0.9需7d，通常一般性暴雨洪澇的衰退需要3～10d，據此，參數a可近似取值0.7～0.9。由于不同地理區域及不同季節降水在地表的滯留時間具有明顯差異，故采用統一的降水衰減參數難免存在偏差，為此，利用湖北省歷年農作物洪澇受災面積資料，基于遺傳算法對降水衰減參數和致澇閾值同時進行率定。其中參數a取值空間設置為0～1，EPthr取值空間設置為10～150，種群規模設置為200，迭代次數設置為100次。優化結果顯示（圖2），降水衰減參數a取值0.825，暴雨洪澇致災閾值EPthr取值70較為合理。基于率定后的EP指數統計年內暴雨洪澇累積強度與歷年農作物洪澇災害受災面積百分比具有較好的線性關系，解釋方差達78.1%。基于超定量法的廣義帕累托分布，以70mm為門限值，對76站1961-2014年歷次暴雨洪澇過程的最大EP值進行概率分布擬合，計算得到1a一遇和5a一遇的重現期對應的EP值分別為140和220mm，由此確定嚴重洪澇和特大洪澇的劃分閾值分別為140和220mm。

圖2 基于遺傳算法的降水衰減參數a（a）和洪澇致災閾值EPthr(b)優化結果頻率統計Fig. 2 Frequency distribution of parameter values selected by the genetic algorithm, for parameter a(a) and for parameter EPthr(b)

2.2 利用有效降水指數評估暴雨洪澇過程方法的檢驗

2.2.1 歷史暴雨洪澇災害的評估與檢驗

利用湖北省1961-2014年76站逐日氣象觀測資料，計算各站逐日EP指數，依據式（2）統計各站暴雨洪澇過程。由于缺乏詳細、可靠的單站驗證資料，故針對全省統計歷年暴雨洪澇過程的累積強度，并與《中國氣象災害大典（湖北卷）》的災情記載及湖北省歷年農作物洪澇受災面積進行對比驗證。

圖3為基于EP指數的湖北省歷年暴雨洪澇累積強度指數，為了區分不同等級類型的洪澇過程，依據一般性洪澇、嚴重洪澇和特大洪澇的劃分閾值分別進行洪澇強度統計，其中一般性洪澇累積強度在統計時包含了嚴重洪澇和特大洪澇過程，同理，嚴重洪澇強度包含特大洪澇。從圖中可以看出，1961年以來的1964、1969、1980、1983、1991、1996、1998、1999、2010年暴雨洪澇累積強度指數相對較高，反映的洪澇程度較為嚴重，這與災情記載中的大澇年一致，對應年份農作物洪澇災害受災面積均在30%以上，嚴重洪澇年主要集中在20世紀90年代，均與實際情況吻合，表明EP指數能客觀反映洪澇災害程度及年際間的差異。從圖還可以看出，1980 和1996年特大洪澇等級強度較弱，主要以一般性洪澇和嚴重洪澇為主，而其它大澇年份均有特大洪澇發生，表明基于不同致澇閾值能夠有效診斷出洪澇過程強度特征。

圖4為基于一般性洪澇臨界閾值和EP指數的湖北省歷年暴雨洪澇過程累積強度空間分布格局，從圖可以看出，EP指數反映的湖北省域范圍暴雨洪澇強度具有明顯的空間差異，頻發重發區域主要位于鄂西南、鄂東南及鄂東北地區，幾近每年發生，鄂西北地區發生頻率較低，程度也相對較輕，這與湖北省的地形和降水分布特征一致。通過對典型洪澇年災害發生范圍驗證比較，EP指數反映的暴雨洪澇空間格局與實際災情相符，洪澇發生的范圍和強度與農作物受災面積也有較好的對應關系。

圖3 基于EP指數的1961-2014年湖北省歷年暴雨洪澇累積強度Fig.3 Annual accumulated severity of rainstorm flood identified by EP index

2.2.2 實時暴雨洪澇過程監測和診斷

以2014年為例，利用EP指數計算湖北省4-10月逐日的區域暴雨洪澇監測指數、站次比及綜合強度指數，繪制時間序列曲線如圖5所示。從圖可以看出，2014年湖北省共發生5次明顯的區域性暴雨洪澇過程，分別在4月中下旬、5月中旬、7月上中旬及9月上旬，其中最強過程為7月上旬，其次為9月上旬，最大過程的站次比達到1/3左右，強度達到全省平均165mm，此外，在5月下旬、6月下旬、8月上旬、9月中下旬還出現了數次局地性暴雨洪澇，以上過程均在民政廳災情統計資料（表1）中得到印證，其中5月9-10日過程監測結果范圍偏大，程度偏重，這是由于降雨過程發生在鄂東南地區，該地區常年多暴雨，當地抗災能力較強，9月10-19日監測結果顯示出現兩次輕度暴雨洪澇過程，但實際災情較重，這與9月上旬以來連續數次暴雨洪澇過程的累積影響有關，表明EP指數對間歇性發生的洪澇過程監測存在不足。總體而言，基于EP指數建立的區域暴雨洪澇監測和評估指標對于暴雨洪澇起止時間和過程強度的動態變化診斷基本合理，但由于對日內降水集中度考慮不足及監測站點密度有限，對局地性過程的診斷存在一定偏差，在監測中有必要結合單站指標進行綜合分析。

圖4 1961-2014年湖北省歷年暴雨洪澇累積強度空間分布格局Fig.4 Spatial pattern of the accumulated severity of rainstorm flood in Hubei province from 1961 to 2014注：圖中標注百分數為農作物洪澇受災面積占當年耕地面積的百分比Note:Percentage annotated on the map represent the ratio of flooded area to total agricultural acreage in that year

圖5 基于EP指數的湖北省2014年暴雨洪澇過程監測Fig. 5 Monitoring of the rainstorm flood process in Hubei province in 2014注：RPI、RFI、RSI分別為區域暴雨洪澇站次比、監測指數和強度指數Note: RPI， RFI and RSI represent flood stations proportion， monitoring index and severity index， respectively

表1 2014年湖北省暴雨洪澇過程及災情Table 1 Rainstorm flood processes and relevant disaster losses in Hubei province in 2014

3 結論與討論

利用洪澇災情資料對經驗參數進行率定后的EP指數，可以解釋農作物洪澇受災面積78.1%的變異，針對1961-2014年歷史暴雨洪澇的評估檢驗，能有效識別出典型大澇年和嚴重洪澇年，合理刻畫出空間分布的差異性，2014年實時暴雨洪澇過程監測結果表明，其對大部分災情信息具有較為一致的響應關系，能動態反映出過程起始時間、結束時間及強度變化，但對局地性較強和間歇性發生的暴雨洪澇過程監測存在一定偏差。暴雨洪澇引發的次生災害種類繁多，限于災情資料收集的困難，本文僅利用湖北省域年尺度的農業災情資料對EP指數進行了初步應用檢驗，表明其在農業洪澇損失評估、年景評價、風險區劃以及作物產量建模中具有一定的應用前景。對于暴雨洪澇引發泥石流、山體滑坡以及城市內澇等的應用檢驗需進一步收集災情資料系統分析。

暴雨洪澇災害的發生遠不止區域內的自然降水一個因素，地理、地形、土壤質地、江湖水位、外來洪水以及防災抗災能力對災害的形成都具有重要的影響[17]。本文在對EP指數進行參數率定中采取通過災情資料反演的方法，未考慮地形因子對地表水分平衡的影響，湖北省地形復雜，山地、丘陵分布廣泛，對洪澇災害的形成具有不同的影響，這是導致EP指數在洪澇實時監測中驗證效果不佳的重要原因之一。因此，在降水衰減參數及致災臨界雨量的確定上，有必要結合地形或分區域考慮更加準確的參數。此外，本文基于日雨量建立EP指數，在資料允許的條件下，可以考慮采用小時雨量，以提高降水集中度對暴雨洪澇形成的反映能力，以及使用更多的加密氣象站資料，從而提高EP指數在暴雨洪澇監測評估中的應用效果。在農業災情的評估中還應盡可能結合農作物土地利用圖以提高評估精度。

參考文獻References

[1]Hirabayashi Y，Mahendran R，Koirala S，et al.Global flood risk under climate change[J].Nature Climate Change， 2013，3（9）: 816-821.

[2]李翠金.中國暴雨洪澇災害的統計分析[J].災害學，1996，11（1）:59-63. Li C J.A statistical analysis of the storm flood disasters in China[J].Journal of Catastrophology，1996，11（1）:59-63.（in Chinese）

[3]馮強，王昂生，李吉順.我國降水的時空變化與暴雨洪澇災害[J].自然災害學報，1998，7（1）:87-93. Feng Q，Wang A S，Li J S.The variation of precipitation in time and space and heavy rain flood disaster in China[J].Journal of Natural Disasters，1998，7（1）:87-93.（in Chinese）

[4]李茂松，李森，李育慧.中國近50年洪澇災害災情分析[J].中國農業氣象，2004，25（1）:40-43. Li M S，Li S，Li Y H.Analysis of flood disaster in the past 50 years in China[J].Agricultural Meteorology，2004，25（1）:40-43. （in Chinese）

[5]周月華，郭廣芬，邵末蘭，等.基于水位和雨量的洪澇受災面積評估模型研究[J].暴雨災害，2007，26（4）:323-327. Zhou Y H，Guo G F，Shao M L，et al.Study on index of heavy rain flood disaster grade in Hubei Province[J].Torrential Rain and Disasters，2007，26（4）:323-327.（in Chinese）

[6]郭廣芬，周月華，史瑞琴，等.湖北省暴雨洪澇致災指標研究[J].暴雨災害，2009，28（4）:357-361.Guo G F，Zhou Y H，Shi R Q，et al.Study on index of heavy rain flood disaster grade in Hubei Province[J].Torrential Rain and Disasters，2009，28（4）:357-361.（in Chinese）

[7]溫泉沛，霍治國，馬振峰，等.中國中東部地區暴雨氣候及其農業災情的風險評估[J].生態學雜志，2011，30（10）:2370-2380. Wen Q P，Huo Z G，Ma Z F，et al.Risk assessment of rainstorm climate and its induced agricultural disasters in east-central China[J].Chinese Journal of Ecology，2011，30（10）:2370-2380. （in Chinese）

[8]Byun H R，Wilhite，D A.Objective quantification of drought severity and duration[J].Journal of Climate，1999，12（9）:2747-2756.

[9]Lu E.Determining the start，duration，and strength of flood and drought with daily precipitation:rationale[J].Geophysical Research Letters，2009，36（12）:L12707.

[10]Deo R，Byun H R，Adamowski J，et al.A real-time flood monitoring index based on daily effective precipitation and its application to Brisbane and Lockyer valley flood events[J]. Water Resources Management，2015，29（6）:1-19.

[11]張國平.有效雨量和滑坡泥石流災害概率模型[J].氣象，2014，40（7）:886-890. Zhang G P.Study on the relation between effective precipitation and landslide/debris-flow with probabilistic model[J]. Meteorological Monthly，2014，40（7）:886-890.（in Chinese）

[12]趙一磊，任福民，李棟梁，等.基于有效降水干旱指數的改進研究[J].氣象，2013，39（5）:600-607. Zhao Y L，Ren F M，Li D L，et al.Study on improvement of drought index based on effective precipitation[J]. Meteorological Monthly，2013，39（5）:600-607.（in Chinese）

[13]鄒燕，葉殿秀，林毅，等.福建區域性暴雨過程綜合強度定量化評估方法[J].應用氣象學報，2014，25（3）:360-364. Zou Y，Ye D X，Lin Y，et al.A quantitative method for assessment of regional heavy rainfall intensity[J].Journal of Applied Meterological Science，2014，25（3）:360-364.（in Chinese）

[14]邵末蘭，張寧，岳陽，等.基于距離函數的區域性暴雨災害風險預估方法研究[J].暴雨災害，2010，29（3）:268-273，278. Shao M L，Zhang N，Yue Y，et al.Study on a disaster risk prediction method for regional heavy rain events based on Euclidean distance function[J].Torrential Rain and Disasters，2010，29（3）: 268-273，278.（in Chinese）

[15]Dong S H.Genetic algorithm based parameter estimation of Nash model[J].Water Resources Management，2008，22（4）: 525-533.

[16]丁裕國，李佳耘，江志紅，等.極值統計理論的進展及其在氣候變化研究中的應用[J].氣候變化研究進展，2011，7（4）: 248-252. Ding Y G，Li J Y，Jiang Z H，et al.Advances in extremes statistics and their application to climate change study[J]. Advances in Climate Change Research，2011，7（4）:248-252. （in Chinese）

[17]謝立華，趙寒冰.洪澇災害與地形的相關性研究:以肇慶市為例[J].自然災害學報，2013，22（6）:240-245. Xie L H，Zhao H B.Correlation between flood disaster and topography:a case study of Zhaoqing City[J].Journal of Natural Disasters，2013，22（6）:240-245.（in Chinese）

Application of Effective Precipitation Index in Rainstorm Flood Disaster Monitoring and Assessment

QIN Peng-cheng, LIU Min, LI Lan
（Wuhan Regional Climate Center， Wuhan 430074， China）

Abstract:Scientific and effective monitoring and assessment are significant to prevent and mitigate the rainstorm flood disasters. Based on the Effective Precipitation Index （EP）， the indices for flood monitoring and assessment at both site and regional scales were established. Using daily meteorological data from 76 stations in Hubei province and disaster loss data related during 1961 to 2014， the decay parameter and critical rainfall were determined， and the performance of the index in historical rainstorm flood assessment and real-time monitoring was examined. The results showed that the calibrated EP index explained 78.1% variances of damaged crop area， and it was sensitive to identify typical extreme and severe flood year from 1961 to 2014. During the real-time rainstorm flood in 2014， the EP index showed a certain capability to determine the start， duration， and strength of each flood process， while there was disagreement for local and discontinuous flood processes. Generally， the EP index was simple in calculation，reasonable in analysis and easy in spread. So it had more advantages in rainstorm flood evaluation， such as historical ranking， annual climatic assessment， disaster loss pre-evaluation， risk zoning， as well as crop yield prediction.

Key words:Effective Precipitation Index; Parameter calibration; Rainstorm flood; Application validation

作者簡介：秦鵬程（1986-），工程師，主要從事氣象災害監測與評估研究。E-mail：qinpengcheng027@163.com

基金項目：中國工程院重大咨詢項目“我國旱澇事件集合應對戰略研究（2012-ZD-13）”；公益性行業科研專項（GYHY201406028）

* 收稿日期：2015-06-14

doi:10.3969/j.issn.1000-6362.2016.01.011