陜西省涇、渭河流域短歷時暴雨頻率分布研究

王晶 王雙銀

摘 要：研究短歷時暴雨頻率分布可為小流域暴雨分析計算提供科學依據，對城市排澇和山洪災害防治具有重要意義。以陜西省涇河和渭河流域39個氣象站1981—2017年8個統計時段年最大值暴雨序列為基礎，以概率點據相關系數法（PPCC）和擬優絕對值準則法（MAE）為評價標準，結合計算暴雨量與實測暴雨量相關圖，以皮爾遜-Ⅲ型分布、指數分布、正態分布、廣義極值分布、廣義帕累托分布和耿貝爾分布為備選頻率分布對陜西省涇河和渭河流域短歷時暴雨頻率分布進行研究，同時對矩法、線性矩法和極大似然法3種參數估算方法進行對比分析。結果表明：各頻率分布采用線性矩法估算參數后的計算值與實測值之間誤差最小，表現出了很好的適用性，是頻率分布參數估算的最佳方法;PPCC法相較于MAE法的最優頻率分布集中，在區域頻率分布選擇時適用性更強;廣義極值分布為最優分布的暴雨序列數和氣象站點數均遠多于其他分布，表現出明顯優勢，可作為陜西省涇河和渭河流域短歷時暴雨頻率計算的理論頻率分布。

關鍵詞：頻率分布;短歷時暴雨;參數估算;陜西省涇、渭河流域

中圖分類號：P333.2 ? 文獻標志碼：A

doi：10.3969/j.issn.1000-1379.2020.07.003

Abstract：Studying on frequency distribution of short-duration rainstorm can provide scientific basis to small watershed rainstorm analysis and is also important to urban drainage and mountain flood disaster control. Based on eight series of the annual maximum rainstorm from 39 meteorological stations in Jinghe and Weihe basin， Shaanxi Province over 1981 to 2017， this paper analyzed the applicability of P-Ⅲ distribution， exponential distribution， normal distribution， generalized extreme value distribution， generalized Pareto distribution and Gumbel distribution， and their parameters assessed by combining the probability point correlation coefficient， mean absolute error and Q-Q methods， and the three parameter estimation methods including MOM， LM and MLE methods were also compared. The results show that the error between calculated and observed values of LM method are the smallest， suggesting that it may be the best method to parameter estimation; the optimal frequency distribution of PPCC method is more concentrated than that of MAE method， suggesting that it is more suitable for selection of regional frequency distribution; and the GEV distribution has obvious advantage in number of both rainstorm sequences and stations， indicating that it is optimal theoretical frequency distribution for short-duration rainstorm in Jinghe and Weihe basins of Shaanxi Province.

Key words： frequency distribution; short-duration rainstorm; parameter estimation; Jinghe and Weihe Basins of Shaanxi Province

短歷時暴雨極易導致城市內澇和山洪災害的發生，威脅著區域生態環境和人民的生命財產安全，對區域經濟社會的持續發展造成嚴重影響。研究短歷時暴雨的頻率分布，可為城鎮排水工程設計、小流域山洪災害防治和溝道綜合治理等提供設計依據，是合理確定工程規模和預警預報措施的基礎，對防治暴雨洪水災害具有重要意義。

對頻率分布的研究有兩派，一派以結合水文現象的物理成因和概率論的有關定理為準則（如耿貝爾分布、對數正態分布等）;另一派則以分布曲線與實測資料擬合最佳為準[1]。由于水文現象受氣象和地理等多種因素的綜合影響，有著錯綜復雜的成因過程，因此很難滿足統計學中的有關定理要求。目前應用最為廣泛的是結合頻率分布曲線統計特點，以實測資料為基礎，通過估計的頻率分布參數計算的“理論”值與實測值的誤差，在一定擬合優度準則下對不同頻率分布進行比較，以符合大多數水文系列的頻率分布作為“理論”頻率分布。我國很多學者研究了單個測站的暴雨頻率分布，發現常用的皮爾遜-Ⅲ型分布（P-Ⅲ）[2-4]、指數分布（EXP）[5]、廣義極值分布（GEV）[6-7]、廣義帕累托分布（GP）[8-9]和耿貝爾分布（GUM）[2，10]都有較好的適用性;也有學者針對區域進行過研究[4，11]，但依據的站點偏少或區域氣候分帶比較明顯，得出的結論可能存在一定的偶然性。

頻率分布參數的估計方法有兩大類，一類是以實測資料的矩為基礎的估計方法，如矩法[12]、線性矩法[11-13]、概率權重矩法[14]、權函數法[15]等，這類方法簡單實用;而另一類是以極大似然原理為基礎的極大似然法，在理論上被認為是最好的參數估算方法[16]，但計算過程復雜，有時候會出現似然方程無解，得不到參數的情況。近年來也有學者將群居蜘蛛優化算法[17]、模擬退火法[18]等智能優化算法用于頻率分布參數估計中，這類方法在計算頻率分布參數時優化函數多，需要編寫相應的計算程序才能實現，在水文分析中還未能廣泛應用。

筆者以正態分布（NORM）和上述常用頻率分布為備選分布對陜西省涇河和渭河流域短歷時暴雨頻率分布進行研究，以期為該區域的設計暴雨計算提供依據。人 民 黃 河 2020年第7期

1 研究區域概況

陜西省涇河和渭河流域介于東經106°19′—110°15′和北緯33°41′—35°31′之間，屬內陸性季風氣候區，暖溫帶半濕潤半干旱氣候帶，四季分明，雨熱同期。渭河干流以北多為樹狀水系，以南多為羽狀水系。區域內地勢南北高中部低，西部高東部低，是一個三面環山、向東敞開的河谷盆地，包括黃土臺塬、黃土高原溝壑和關中盆地3種地貌類型。受復雜地形的影響，暴雨、洪水、漬澇等水災一直是區域內主要的自然災害，尤其近年來，在人類活動和氣候變化的綜合影響下，暴雨導致的城市內澇、山洪災害等時有發生。

研究區為涇河和渭河陜西省部分，研究區域范圍及氣象站點分布見圖1。

基本資料為陜西省氣象中心提供的陜西省涇河和渭河流域39個氣象站1981—2017年歷年10、20、30、45、60、90、120、180 min統計時段的年最大暴雨量數據。

2 研究方法

（1）參數估算方法。結合數理統計理論，從簡單實用的角度考慮，選用矩法（MOM）、極大似然法（MLE）和線性矩法（L-M）作為各頻率分布的參數估算方法，具體計算方法見文獻[19]。

（2）擬合優度檢驗方法。通常采用擬優平方和準則法（RMSE）和擬優絕對值準則法（MAE）進行擬合優度檢驗，從實測資料的擬合情況看，MAE較為合理，且可避免個別點據對成果夸大的影響[20]。RMSE一般適用于長重現期的長歷時設計暴雨和設計洪水的頻率分析，而短歷時設計暴雨的重現期一般不超過20 a[21]，實測暴雨資料的整體擬合情況對其影響較大。筆者選用MAE，并結合概率點據相關系數法（PPCC）和圖形分析法（Q—Q圖）進行擬合優度檢驗。

3 結果與分析

3.1 參數估算方法對比分析

3種參數估算方法估算的各測站不同統計時段暴雨序列同一頻率分布的參數有一定差異，表1列出了乾縣站最大10 min暴雨序列3種參數估算方法所得各頻率分布的參數，可以看出，整體上MOM和L-M估算的參數比較接近，MLE估算的參數與MOM和L-M估算的參數差異較大。在用MLE估算P-Ⅲ分布參數時，會出現似然方程無解而得不到參數的情況（長武站180 min序列和陳倉站20～120 min序列等）。

為了對比各參數估算方法的優劣，分別計算3種參數估算方法下各假定頻率分布計算暴雨量與實測暴雨量的相關系數r與差值平均絕對值MAE*，采用概率點據相關系數法和擬優絕對值準則法優選最佳參數估算方法。從39個測站8個統計時段暴雨序列在3種參數估算方法下各頻率分布的最優方法站點數看（表2），采用PPCC檢驗法時，P-Ⅲ分布、EXP分布、GEV分布和GP分布在各序列中以L-M最優的站點數最多，均在17～26之間;GUM分布除180 min統計時段暴雨序列MOM的最優站點數最多外，其他7個統計時段暴雨序列也是L-M最優的站點數最多;對于NORM分布，以L-M、MLE和MOM為最優參數估算方法在各統計時段暴雨序列中均有一定數量，且最優站點數差異不大，但以L-M最優的總序列數最多，占計算序列數的38%，略微優于MLE和MOM;采用擬優絕對值準則法時，P-Ⅲ分布、EXP分布、NORM分布、GEV分布和GUM分布在各序列中以L-M最優的站點數最多，在19～38之間;對GP分布，L-M、MLE為最優參數估算方法在各歷時暴雨序列中均有一定數量，且最優站點數差異不大，但以L-M最優的總序列數最多，占計算序列數的52.9%，略微優于MLE。

綜合兩種優選方法的結果，L-M對6種頻率分布均有很好的適用性，可作為短歷時暴雨頻率分布對比分析的統一參數估算方法。當以GP分布為理論分布進行短歷時暴雨頻率計算時，L-M和MLE均可用于參數估算，而以其他分布為理論分布進行短歷時暴雨頻率計算時，建議采用L-M估算參數。

3.2 頻率分布適用性

根據前述結果，用L-M法估算各測站不同統計時段暴雨序列6種分布的參數，根據估算的參數分別計算6種頻率分布在各序列假定頻率分布計算值與實測值的r值和MAE*值，以此篩選各序列和各測站暴雨的最適用（最優）頻率分布。

3.2.1 各統計時段暴雨序列頻率分布

通過統計研究區域8個統計時段暴雨序列的最優頻率分布站點數（表3）可以看出，以PPCC為評價標準時，各序列最優頻率分布均以GEV分布最優站點數量最多，占比在44%～56%之間，且隨著統計時段長度的增加占比有增大趨勢，其次為GP分布，占比在13%～31%之間。以MAE為評價標準時，小于等于30 min序列中GP分布優勢站點最多，占比在33%～41%之間，其次為GEV分布和P-Ⅲ分布，且兩者占比接近;120 min序列中P-Ⅲ分布優勢站點最多，其次為GEV分布，二者占比差距僅為8%;其他時段以GEV分布最優站點數量最多，且從各頻率分布最優序列數量看，GEV分布、GP分布和P-Ⅲ分布為最優的序列數分別為103、93和86，GEV分布占比最大。綜合考慮PPCC和MAE的評價結果，認為GEV分布可作為各時段短歷時暴雨序列的最優頻率分布。

從評價結果的一致性來看，以PPCC為評價標準時，各時段間6種頻率分布最優站點數量穩定且各頻率分布間8個時段最優站點數量差距顯著，最優分布更集中，評價結果一致性強;以MAE為評價標準時，GEV分布、GP分布和P-Ⅲ分布在8個時段最優站點數量差距較小，沒有顯著差異。可見，PPCC檢驗法在區域頻率分布選擇時適用性更強。

3.2.2 各站點暴雨序列頻率分布

研究同一站點不同統計時段暴雨序列的頻率分布時，盡可能選用同一頻率分布。當以PPCC和MAE評價的某一暴雨序列結果不一致時，通過繪制兩種評價標準適用分布的計算暴雨量與實測暴雨量相關圖（Q—Q圖）選擇最優分布。表4列出了乾縣站各統計時段暴雨序列在兩種評價方法下的適用分布，對90、120、180 min暴雨序列，兩種評價方法的適用分布均為GEV分布，則這3個歷時的最優分布確定為GEV分布;而其他統計時段暴雨序列兩種方法的評價結果不一致，繪制10 min暴雨序列的Q—Q圖發現兩種適用分布在中值區很接近，如圖2所示，但GEV分布在低值區的離差大，高值區離差小，GP剛好相反，考慮到暴雨頻率計算的目的是推求設計暴雨，對中大值的精度要求較高，選取GEV分布作為乾縣站10 min的最優頻率分布;同樣選定的20、30、45、60 min暴雨序列的最優分布見表4。從8個統計時段暴雨序列選定的最優分布看，GEV分布有5個，相比其他兩個分布具有明顯優勢，據此以GEV作為乾縣站短歷時暴雨頻率計算的理論分布。

采用同樣的分析過程，所得研究區域內各測站最優頻率分布見表5，在39個測站中，最優分布為GEV分布的站點有22個，占總數的56%，其次是GP分布、P-Ⅲ分布、EXP分布和GUM分布，其站點數量分別為12、5、3、1個，其中寶雞站和千陽站對GP分布和P-Ⅲ分布同樣適用，武功站對EXP分布和GP分布同樣適用，興平站對GEV分布和EXP分布同樣適用，且GP分布、P-Ⅲ分布、EXP分布和GUM分布零散地分布在研究區域內。

同一測站不同統計時段暴雨序列和同一統計時段不同測站暴雨序列兩個角度的分析結果均表明，不同評價標準下適用分布不同，GEV分布和GP分布適用性均比較明顯，且GEV分布在8個短歷時暴雨序列和39個站點中最優的數量比例最高，推求的暴雨精度相對較高，其適用性更強。在對陜西省涇河和渭河流域進行短歷時暴雨頻率計算時，建議選用GEV分布作為理論頻率分布。

4 結 論

（1）各頻率分布采用L-M估算參數后的計算值與實測值之間的誤差最小，表現出了很好的適用性，是頻率分布的最優參數估算方法。

（2）PPCC法相較于MAE法的最優頻率分布集中，在區域頻率分布選擇時適用性更強。

（3）GEV分布為最優分布時的暴雨序列數和站點數均遠多于其他分布為最優時的數量，表現出明顯優勢，且據此推求的短歷時暴雨精度相對較高，建議在研究區域短歷時暴雨頻率分析計算時把GEV分布作為理論頻率分布。

參考文獻：

[1] 金光炎.水文統計理論與實踐[M].南京：東南大學出版社，2012：114-150.

[2] 朱齊艷，張端虎.韶關站年最大值法暴雨強度公式頻率曲線對比[J].人民珠江，2014，35（6）：70-72.

[3] 趙琳娜，王彬雁，白雪梅，等.北京城市暴雨分型及短歷時降雨重現期研究[C]//第33屆中國氣象學會年會S9水文氣象災害預報預警.北京：中國氣象學會， 2016：1-3.

[4] 毛慧琴，杜堯東，宋麗莉.廣州短歷時降水極值概率分布模型研究[J].氣象，2004，30（10）：3-6.

[5] 盧金鎖，程云，鄭琴，等.西安市暴雨強度公式的推求研究[J].中國給水排水，2010，26（17）：82-84.

[6] 樊利利，王艷永.廣義極值分布的參數估計及實例分析[J].首都師范大學學報（自然科學版），2017，38（3）：13-18.

[7] 金光炎.廣義極值分布及其在水文中的應用[J].水文，1998，18（2）： 10-16.

[8] 趙瑞星，翟宇梅.極端降水廣義帕累托分布參數的Pickands自助矩估計研究[J].水力發電學報，2015，34（10）：42-50.

[9] 馬京津，李書嚴，王冀.北京市強降雨分區及重現期研究[J].氣象，2012，38（5）：569-576.

[10] 楊譯舒，葉興成，王飛，等.基于不同頻率分布的樂山市暴雨參數分析[J].災害學，2017，32（4）：214-218.

[11] 牛晨，王雙銀，吳素娟.基于線性矩法的短歷時暴雨頻率分布研究[J].長江科學院院報，2019，36（8）：55-60.

[12] 梁玉音，劉曙光，鐘桂輝，等.線性矩法與常規矩法對太湖流域降雨頻率分析的比較研究[J].水文，2013，33（4）：16-21.

[13] HOSKING J R M， WALLIS J R. Regional Frequency Analysis： An Approach Based on L-Moments[M]. Cambridge：Cambridge University Press，1997：720-721.

[14] 雷冠軍，王文川，殷峻暹，等. P-Ⅲ型曲線參數估計方法研究綜述[J].人民黃河，2017，39（10）：1-7.

[15] 馬秀峰，阮本清.對權函數法求概率分布參數的討論[J].水利學報，2001，32（11）：34-40.

[16] 李航，宋松柏，石繼海.指數Gamma分布參數估計方法對比研究[J].水力發電學報，2019，38（4）：96-107.

[17] 吳光瓊.群居蜘蛛優化算法在水文頻率曲線參數優化中的應用[J].水資源與水工程學報，2015，26（6）：123-126，131.

[18] 李宏偉.水文頻率參數計算方法與應用研究[D].楊凌：西北農林科技大學，2009：37.

[19] 黃振平，陳元芳.水文統計學[M].北京：中國水利水電出版社，2011：168-186.

[20] 李松仕.幾種頻率分布線型對我國洪水資料適應性的研究[J].水文，1984，4（1）：1-7.

[21] 中華人民共和國住房和城鄉建設部.城市防洪工程設計規范：GB/T 50805—2012[S].北京：中國計劃出版社，2012：2.

[22] HEO J H， KHO Y W， SHIN H， et al. Regression Equations of Probability Plot Correlation Coefficient Test Statistics from Several Probability Distributions[J]. Journal of Hydrology， 2008， 355（1-4）：1-15.

[23] 李興凱，陳元芳.暴雨頻率分布線型優選方法的研究[J].水文，2010，30（2）：50-53.

【責任編輯 張 帥】