水汽風(fēng)速放大計算可能最大暴雨方法

梁忠民+季俊杰+卜慧+劉甜+楊靖

摘要：我國可能最大暴雨（PMP）計算中，通常采用兩種放大方法，即水汽效率放大和水汽風(fēng)速放大，但一般認為水汽風(fēng)速放大成果不穩(wěn)定，結(jié)果往往偏大，所以通常只作參考，最后基本上都選取水汽效率的放大結(jié)果作為最終的設(shè)計成果。借鑒港口工程設(shè)計中“風(fēng)玫瑰圖”概念，對傳統(tǒng)水汽風(fēng)速放大方法進行改進，即先根據(jù)研究區(qū)域常風(fēng)向?qū)崪y風(fēng)速資料進行風(fēng)向分組，對每組風(fēng)向的典型暴雨進行放大，再從中選取最大值作為水汽風(fēng)速的放大成果，為PMP計算提供了一種新的思路。以湖北省咸寧核電廠為例進行了應(yīng)用研究，并與水汽效率放大方法和傳統(tǒng)水汽風(fēng)速放大方法結(jié)果進行了比較，提供了更為合理的分析成果。

關(guān)鍵詞：可能最大暴雨P(guān)MP；水汽風(fēng)速放大；水汽效率放大；風(fēng)玫瑰圖

中圖分類號：P339文獻標(biāo)志碼：A文章編號：1672-1683（2017）01-0014-04

Abstract：Generally，moisture and efficiency factor maximization or moisture and wind maximization are two approaches used to estimate the probable maximum precipitation （PMP） in China.But only the value obtained from the moisture and efficiency factor maximization is used as the PMP estimation in most cases，because the moisture and wind maximization approach regularly overestimates and therefore is usually taken as reference only.This paper presents a modified moisture and wind maximization method based on the "Wind Rose Diagram" applied in port engineering.First，measured wind speed and direction data are grouped according to the direction of prevailing wind in the studied area.Then，the high-efficiency storm in each group of wind direction is maximized.Finally，the maximum value among all groups is chosen to be the final result of moisture and wind maximization，that is PMP value.The new approach was tested on the Xianning nuclear power plant in Hubei Province，and generated a more reasonable PMP value compared with the traditional approaches.

Key words：probable maximum precipitation （PMP）；moisture and wind maximization；moisture and efficiency factor maximization；Wind Rose Diagram

可能最大暴雨/洪水（PMP/PMF）作為重要水庫大壩的防洪設(shè)計標(biāo)準(zhǔn)，在世界范圍內(nèi)開展了大量的研究與應(yīng)用[1]。1973年，世界氣象組織（WMO）正式出版了可能最大降水估算手冊[2]，隨后分別于1986年和2009年對手冊進行了修訂完善[3-4]。我國開展PMP/PMF的研究始于20世紀(jì)50年代，“75.8”河南特大暴雨后，我國加快了PMP/PMF的研究和應(yīng)用工作，其中，理論研究成果以詹道江和鄒進上編著的《可能最大暴雨與洪水》[5]、王國安編著的《可能最大暴雨和洪水計算原理與方法》[6]最具代表性。目前，PMP/PMF已作為我國水利水電工程防洪設(shè)計的最高校核標(biāo)準(zhǔn)[7]，也是核電工程洪水的設(shè)防標(biāo)準(zhǔn)[8]。

近年來，國外學(xué)者在PMP的計算方法上取得了一些新進展。Ohara等利用區(qū)域氣候模式MM5模擬歷史極端暴雨，并通過調(diào)整模型的初始及邊界條件實現(xiàn)暴雨放大以推求PMP[9]；Ishida等采用MM5模式與相對濕度最大化相結(jié)合的途徑，計算了北加利福尼亞三個小流域的PMP[10]；Beauchamp等采用加拿大區(qū)域氣候模式CGCM計算了流域的可降水量[11]。

目前我國工程實際中的PMP計算方法包括水汽效率放大和水汽風(fēng)速放大兩類。相較而言，水汽效率放大法推求的PMP結(jié)果相對穩(wěn)定，而水汽風(fēng)速放大推求PMP，由于風(fēng)速的日變化較大，最大風(fēng)速與最小風(fēng)速的比值通常也較大，往往導(dǎo)致根據(jù)最大比值進行典型暴雨放大估算PMP的計算結(jié)果也偏大[12-13]。所以，很多實際工程中都以水汽效率放大方法的結(jié)果為主，水汽風(fēng)速放大方法結(jié)果通常只作參考，限制了該方法的應(yīng)用。水汽風(fēng)速放大法結(jié)果偏大的原因，除了風(fēng)速的日變化較大外，還由于該方法取所有方向上最大與最小風(fēng)速的比值進行典型暴雨放大之故，而這種做法與一個地區(qū)實際的風(fēng)向特點也不甚相符。本文借鑒海洋工程水文計算中“風(fēng)玫瑰圖”的概念[14]，先將一個地區(qū)的水汽風(fēng)速資料按照常風(fēng)向分為若干個主要的水汽入流方向，在每個方向上分別進行水汽風(fēng)速的極大化及典型放大，再取所有方向數(shù)值中的最大值（即外包值）作為PMP的估算結(jié)果。以湖北咸寧核電廠PMP的估算為例進行了應(yīng)用研究，通過與其它方法的對比分析，表明本文方法是合理可行的。

1 水汽風(fēng)速放大推求PMP方法

1.1 傳統(tǒng)的水汽風(fēng)速放大方法

水汽風(fēng)速放大方法分為按水汽W和風(fēng)速V分別放大后再相乘，以及水汽入流指標(biāo)（水汽與風(fēng)速乘積VW）放大兩種，但其基本假設(shè)是相同的，即認為降水P與V及W具有正相關(guān)關(guān)系[5]，該方法的主要步驟如下。

（1） 選取水汽入流方向，計算代表層風(fēng)速。先根據(jù)流域氣象條件規(guī)定水汽入流的大致方向，選取各氣象站暴雨當(dāng)天及前兩天（具體選取時段視情況而定）的風(fēng)向風(fēng)速，合成每一天0 h和12 h兩個時刻的風(fēng)速，最后挑選出三天中合成風(fēng)速最大且風(fēng)向合理的風(fēng)速作為典型暴雨的代表層風(fēng)速Vd。

（2） 計算代表性可降水量。選取水汽入流方向上的露點代表站暴雨當(dāng)天及前兩天（具體選取時段視情況而定）的露點資料，計算地面持續(xù)12 h最大露點，并按飽和假絕熱線化算到1 000 hPa，各站取平均即為代表露點，按此露點和流域平均高程推求典型暴雨可降水量Wd。

（3） 水汽風(fēng)速或水汽輸送率放大推求PMP。水汽風(fēng)速放大的極大化指標(biāo)為VmWm，水汽輸送率放大的極大化指標(biāo)為 （VW）m，放大公式如下：

1.2 改進的水汽風(fēng)速放大方法

本文根據(jù)海岸工程中“風(fēng)玫瑰圖”的概念，提出一種改進的水汽風(fēng)速放大計算PMP的思路?！帮L(fēng)玫瑰圖”也叫“風(fēng)向頻率玫瑰圖”，它是根據(jù)某一地區(qū)多年平均統(tǒng)計的各個風(fēng)向出現(xiàn)的頻率，并按一定比例繪制，通常分為8個方向或16個方向，見圖1。

現(xiàn)行的水汽風(fēng)速放大推求PMP方法中，在利用公式（1）或（2）時，是由可能最大風(fēng)速與所有方向上代表層風(fēng)速中的最小風(fēng)速相比，作為放大倍比，所以這個比值往往會很大，導(dǎo)致計算的PMP偏大。本文借鑒風(fēng)玫瑰圖概念，提出先將實測場次暴雨按常風(fēng)向進行分組，劃分成若干個水汽入流方向，再在每個方向上分別選取極大化指標(biāo)VmWm或（VW）m進行放大，最后再取各方向的外包值作為PMP的結(jié)果。由于實測場次暴雨資料中最大風(fēng)速與最小風(fēng)速出現(xiàn)在一個方向的概率較小，使得計算的放大倍比不至于過大，可以避免推算PMP結(jié)果過大現(xiàn)象。

2 應(yīng)用實例

2.1 流域基本情況

（1）基本情況。

咸寧核電廠位于湖北省咸寧市通山縣大畈鎮(zhèn)境內(nèi)大坑村、官塘村附近的“獅子巖”，坐落在富水水庫中部北岸之濱。廠址位于富水流域，富水為長江中游下段南岸一級支流，流域地勢西南高東北低。

（2）暴雨及水汽來源。

研究區(qū)的暴雨一般多為渦切變暴雨，也有冷鋒和臺風(fēng)暴雨，暴雨影響系統(tǒng)基本分為三類：一是盛夏臺風(fēng)在副熱帶高壓南緣西行影響湖北省，二是低渦沿副高北緣的切變線東移造成暴雨，三是蒙古槽由貝加爾湖帶來北路冷空氣結(jié)合低層低渦沿副高外圍形成暴雨帶。

據(jù)暴雨天氣分析，咸寧地區(qū)水汽來源及輸送的形式主要包括：①西南季風(fēng)。初夏時節(jié)，槽前西南氣流送來孟加拉灣的暖濕空氣，大暴雨常發(fā)生在鄂東和江漢平原東部；②東南水汽。盛夏，由于東風(fēng)帶系統(tǒng)活動頻繁，水汽多來自東海，暴雨區(qū)集中在鄂西；③形成于東?；蚰虾５呐_風(fēng)裹挾大量水汽從我國東南沿海登陸，一路向西北深入內(nèi)陸到達核電廠所在區(qū)域，并在地形作用下形成大暴雨。

因此，可將研究區(qū)的水汽入流方向大致分為東南、南和西南三個方向，見圖2。

2.2 資料選用

（1）歷年最大24 h雨量資料。

本次選取了選擇暴雨資料完整、代表性較好、離廠址較近的官塘站作為控制站，官塘站建站時間為1973年，其附近的陽辛站（距官塘約13 km）的建站年份比官塘站略早（但資料間斷）。根據(jù)《水利水電工程設(shè)計洪水規(guī)范》（SL 44-2006），當(dāng)雨量站相距較近時，可直接將鄰站資料移用于設(shè)計站。因此，將1961年-1972年陽辛站資料移用于官塘站，最終得到1961年-2013年共53年的歷年最大24 h雨量資料系列。設(shè)計流域雨量站位置及水系見圖3。

（2）氣象資料。

高空氣象站選取衢州、南昌、贛州、長沙和貴陽5個站，選定項目包括1961年-2013年的風(fēng)向風(fēng)速資料及1961年-1980年的高空露點資料，因核電廠廠址周圍露點平均高程約為85 m，因此選取850 hPa高空的風(fēng)向風(fēng)速作為代表層資料。其中，東南方向上的高空氣象代表站是衢州站和南昌站，正南方向上是南昌站和贛州站，西南方向上是長沙站和貴陽站，高空氣象站分布見圖4。地面露點站選南昌、修水、宜春、岳陽和長沙5個站，選定項目為1980年-2013年的地面露點，其中，東南方向上的代表站是南昌站，正南方向上是修水站和宜春站，西南方向上是岳陽站和長沙站，地面露點站分布見圖5。

2.3 水汽風(fēng)速放大推求PMP

（1）代表露點和最大可降水計算。

計算持續(xù)12 h最大露點，按飽和假絕熱化算到1 000 hPa，在水汽入流方向上各站取平均，即為各場降雨的代表露點。由于露點比較穩(wěn)定，一般在30年以上記錄中的持續(xù)最大露點所相應(yīng)的水汽含量就接近PMP時的水汽含量。所選擇的地面氣象站南昌、修水、宜春、岳陽、長沙均具有30年以上的地面露點資料，從各站歷年1000 hPa露點資料中選取持續(xù)12 h最大Tdm，南昌站為28.1 ℃、修水站為26.6 ℃、宜春站為27.2 ℃、岳陽站為27.9 ℃、長沙站為27.8 ℃，均小于本流域的水汽源地西太平洋的海洋表面水溫?？紤]工程安全，故可能最大露點取六站的最大值Tdm=28.1 ℃，按此露點計算出咸寧核電廠廠址周圍露點平均高程85 m至12 000 m高空的可降水101.7 mm，即為最大可降水Wm。

（2）代表層風(fēng)速及水汽風(fēng)速放大推求PMP。

根據(jù)前述分析，咸寧核電廠所在地區(qū)的水汽入流共分為東南、正南和西南三個方向，在每個方向上分別逐年確定年最大24 h暴雨當(dāng)天、前2 d共72 h，0時和12時的合成風(fēng)速，再從多年風(fēng)速系列中選取最大值作為該方向的歷史最大風(fēng)速。由于本次的資料系列較長（53年），可以將歷史最大風(fēng)速近似作為該方向上的可能最大風(fēng)速。據(jù)此，可以求得各水汽入流方向上的可能最大風(fēng)速。

分別選取極大化指標(biāo)VmWm和（VW）m，按傳統(tǒng)水汽風(fēng)速放大方法和本文的改進方法推求咸寧核電廠可能最大24 h暴雨，計算結(jié)果見表1。

結(jié)果表明，改進的水汽風(fēng)速或水汽輸送率放大方法計算的PMP比傳統(tǒng)方法的結(jié)果降低較多。根據(jù)頻率分析結(jié)果，咸寧核電廠址24 h萬年一遇設(shè)計暴雨為505.5 mm，本文改進方法推求的可能最大24 h暴雨為576.7 mm，其與萬年一遇的比值為1.14，該比值與我國大部分地區(qū)PMP與萬年一遇比值介于1.1～1.2的認識較為一致；而根據(jù)水汽效率放大方法推求的24 h PMP成果為391 mm，小于萬年一遇設(shè)計值。因此，建議采用改進的水汽輸送率放大計算的PMP結(jié)果作為當(dāng)?shù)乇┯攴糯蠓ǖ腜MP成果，即24 h PMP為576.7 mm。當(dāng)然，咸寧核電廠址最終PMP值還需結(jié)合移置暴雨放大后綜合確定。

3 結(jié)論

本文借鑒“風(fēng)玫瑰圖”的概念，對可能最大降水PMP計算中的水汽風(fēng)速放大方法進行改進。先根據(jù)暴雨天氣系統(tǒng)形成機制，確定研究區(qū)的主要水汽入流方向，并對實測場次暴雨資料進行水汽入流方向分組；在每個水汽入流方向上，分別確定水汽風(fēng)速放大指標(biāo)VmWm或（VW）m，再取各方向上的最大值以推求PMP。由于是先在各水汽入流上進行極大化然后再取外包，所以改進方法可避免目前水汽風(fēng)速方法結(jié)果往往偏大的問題。在湖北咸寧核電廠的應(yīng)用表明，本文方法推求的結(jié)果與目前相應(yīng)方法相比，更具合理性。

采用水汽風(fēng)速放大推求PMP，其前提條件是降水量與風(fēng)速或水汽輸送率應(yīng)具有正相關(guān)關(guān)系，但在實際工作中這種正相關(guān)達到多大時才能使用水汽風(fēng)速方法，尚無定論，需要對相關(guān)問題進行深入研究，以提高PMP計算的可靠性。

參考文獻（References）：

[1] 王國安.國內(nèi)外PMP /PMF的發(fā)展和實踐[J].水文，2004，24（5）：5-9，47.（WANG Guo-an.Worldwide development and practice of PMP /PMF[J].Journal of China Hydrology，2004，24（5）：5-9，47.（in Chinese））

[2] World Meteorological Organization （WMO）.Manual for estimation of probable maximum precipitation[M].Geneva.No.332，1973.

[3] World Meteorological Organization （WMO）.Manual for Estimation of Probable Maximum Precipitation[M].2nd Edition，Geneva.No.332，1986.

[4] World Meteorological Organization （WMO）.Manual for estimation of probable maximum precipitation[M].Geneva.No.1045，2009.

[5] 詹道江，鄒進上.可能最大暴雨與洪水[M].北京：水利電力出版社，1983.（ZHAN Dao-jiang，ZOU Jin-shang.Probable maximum precipitation and flood[M].Beijing：Water Resources and Electric Power Press，1983.（in Chinese））

[6] 王國安.可能最大暴雨和洪水計算原理和方法[M].北京：中國水利水電出版社，鄭州：黃河水利出版社，1999.（WANG Guo-an.Principles and methods of PMP/PMF calculations[M].Beijing：China Water & Power Press，Zhengzhou：The Yellow River Water Conservancy Press，1999.（in Chinese））

[7] 中華人民共和國水利部.《防洪標(biāo)準(zhǔn)》-GB50201-94[M].北京：中國計劃出版社，1994.（The Ministry of Water Resources of the People′s Republic of China.Flood control standard-GB50201-94[M].Beijing：China Planning Press.1994.（in Chinese））

[8] 國家核安全局.濱河核電廠廠址設(shè)計基準(zhǔn)洪水的確定 HAD101/08，核安全導(dǎo)則匯編 （上冊）[M].北京：中國法制出版社，2000.（National Nuclear Safety Administration.The design basis flood Determination of Riverside Nuclear Power Plant site HAD101/08，Nuclear Safety Guide compilation （Vol 1）[M].Beijing：China Legal Publishing House，2000.（in Chinese））

[9] N Ohara，M L Kavvas，S Kure，et al.Physically based estimation of maximum precipitation over American river watershed[J].Journal of Hydraulic Engineering.2011，16 （4）：351-361.

[10] K Isdida，M L Kavvas，S Jiang，et al.Physically Based Estimation of Maximum Precipitation over Three Watersheds in Northern California：Relative Humidity Maximization Method[J].Journal of Hydraulic Engineering.2015，10.1061/（ASCE）HE.1943-5584.0001175，04014052.

[11] J Beauchamp，R Leconte，M Trudel，et al.Estimation of the summer-fall PMP and PMF of a northern watershed under a changed climate[J].Water Resources Research.2013，49 （6）：3852-3862.

[12] 華家鵬，楊輝，王達雨，等.江坪河水電站可能最大暴雨研究[J].河海大學(xué)學(xué)報：自然科學(xué)版，2004（5）：523-525.（HUA Jia-peng，YANG Hui，WANG Da-yu，et al.Study on probable maximum storm at Jiangping River hydropower station[J].Journal of Hohai University：Natural Sciences，2004（5）：523-525.（in Chinese））

[13] 梁忠民，鐘平安，華家鵬.水文水利計算（第2版）[M].北京：中國水利水電出版社，2008.（LIANG Zhong-min，ZHONG Pin-an，Hua Jia-peng.Hydrological design and water conservancy calculation （Second Edition）[M].Beijing：China Water & Power Press，2008.（in Chinese））

[14] 邱大洪.工程水文學(xué)[M].北京：人民交通出版社，1999.（QIU Da-hong.Engineering Hydrology[M].Beijing：China Communications Press，1999.（in Chinese））