李桃英,王紅艷,曹逸希
(1.陜西省水文水資源勘測中心,陜西 西安 710068;2.陜西省西安水文水資源勘測中心,陜西 西安 710100;3.陜西省水務集團水生態(tài)綜合開發(fā)有限公司,陜西 西安 710000)
20世紀70年代陜西省編制了《陜西省可能最大暴雨圖集》《陜西省暴雨洪水圖集》以及20 世紀80 年代各地市編制了《地區(qū)水文手冊》,這些成果在水工程規(guī)劃、設計以及在各級防汛部門發(fā)揮了重要作用。隨著時間的推移,近30 年來發(fā)生了許多大暴雨,有的歷時實測暴雨量超過全省極值,也有暴雨量超過全國極值,原先的成果已不適應,因此陜西省于2020年重新修編了《陜西省水文手冊》,其中對暴雨頻率進行了重新計算,得出不同歷時的暴雨統(tǒng)計參數(shù),在此基礎上本文就暴雨統(tǒng)計參數(shù)特征進行分析。
受西太平洋副熱帶高壓與西風帶低槽相互作用的影響,暴雨是陜西省夏季多見的一種災害性天氣,具有明顯的季節(jié)變化,每年因暴雨引起山洪暴發(fā)和滑坡、泥石流,導致良田被毀,水庫跨壩,鐵路、公路、橋梁破壞,交通中斷等。陜西省暴雨時間跨度由南向北逐漸縮短,多集中于7、8、9三個月,陜南最早發(fā)生在3月27 日、最遲在11月19日,關中最早在4月5日、最遲在11月7日,陜北最早在4月23日、最遲在10月7日,連續(xù)暴雨天數(shù)由南向北越來越短。同時受山地地形條件影響,省內形成幾個暴雨高值區(qū)和暴雨低值區(qū),米倉山、大巴山是暴雨最多的地區(qū),秦嶺、鶻嶺、蟒嶺是次多暴雨區(qū),子午嶺、黃龍山及窟野河中下游是相對較多暴雨區(qū),白于山以北的風沙區(qū)是暴雨最少的地區(qū)。暴雨特性各地差異較大,局部地區(qū)災害性暴雨出現(xiàn)頻繁,暴雨量自南向北呈遞減趨勢。
結合工程使用方便和《中國暴雨參數(shù)圖集》中的標準歷時,確定本次采用的標準歷時為10 min、20 min、60 min、3 h、6 h、12 h、24 h及3 d等8 個時段。
(1)選取各歷時年最大值為樣本值。
(2)年最大時段雨量盡可能根據(jù)分段觀測間隔時間較短的雨量資料選取,以便可以滑動選取最大值。
(3)系列缺測年份的插補:中等以下暴雨不進行插補,按連續(xù)系列對待;對于大暴雨用相鄰站相關關系進行插補,取得缺測年份的暴雨數(shù)值。
(4)當系列中有個別或少數(shù)年份無法滑動選樣時,可采用下列方法改正。
①參考相鄰自記雨量站的分段雨量比例分配本站的雨量,使之可作滑動選樣。
②缺少24 h雨量值時,可利用日雨量資料近似估算。
③根據(jù)固定觀測段制、并記錄有降水起訖時間的雨量記錄,也可以判斷出個別標準歷時的年最大雨量。
(5)暴雨移置
若調查暴雨發(fā)生站點與臨近雨量站相距較近,暴雨天氣形勢相似,地形、地理條件基本相同,可直接將調查暴雨搬移到臨近雨量站使用;若調查暴雨發(fā)生站點與附近雨量站相距較遠,地形、地理條件差異顯著,則必須進行修正。
(6)特大暴雨及特大暴雨重現(xiàn)期
特大暴雨對頻率分析成果具有極為重要的影響。雖然大多數(shù)站系列中缺乏大暴雨資料,但有些站在不太長的系列中出現(xiàn)了重現(xiàn)期遠遠大于按系列年數(shù)求得的重現(xiàn)期,必須在進行頻率計算前對特大值先行處理和重現(xiàn)期的估算。
(7)系列代表性分析
由于暴雨資料為極值,因此,各雨量站資料系列不需要統(tǒng)一延長,資料系列長度采用從建站至2012 年。分析短系列暴雨(n=30 年)對長系列(n≥50 年)暴雨的代表性,共選209個長系列站的資料,對24 h雨量均值和頻率P=1%設計值,進行短系列(1983 年~2012 年)相對于長系列(建站至2012 年)計算結果誤差統(tǒng)計分析。24 h雨量均值相對誤差在±5%以內的142 站,占67.9%;相對誤差在±10%以內的站201站,占96.2%。頻率P=1%設計值相對誤差在±10%以內的187站,占89.5%;相對誤差在±20%以內的204站,占97.6%。因此,所選資料具有較好的代表性。
但是,短系列均值及設計值正負誤差不對稱,正誤差所占比例大于負誤差所占比例,但負誤差的最大值大于正誤差的最大值,并且正負誤差空間分布具有明顯的地域分布特征;黃河流域系統(tǒng)偏小,而且最大負誤差值出現(xiàn)在幾個區(qū)域,依次為北洛河中游劉家河、涇河北岸旬邑、北洛河下游狀頭一帶、榆林地區(qū)北部,長江流域主要在安康漢濱區(qū)一帶;正誤差最大值區(qū)出現(xiàn)在渭河流域寶雞一帶、北洛河中游延安銅川交界處一帶、無定河中游一帶;其他正負誤差交替出現(xiàn),總體上長江流域正誤差多于黃河流域正誤差。
《水利水電工程設計洪水計算規(guī)范》(SL 44-2006)(以下簡稱《規(guī)范》)和《中國暴雨統(tǒng)計參數(shù)圖集》規(guī)定,單站統(tǒng)計參數(shù)的計算按年最大值法選樣,以皮爾遜Ⅲ型曲線作為總體分布模型,CS/CV值采用3.5。采用8 個標準時段在同一張頻率紙上調試頻率曲線,其中10 min、20 min、60 min按分鐘滑動,3 h、6 h、12 h、24 h按小時滑動,3 d按日滑動統(tǒng)計計算。
單站各歷時暴雨的統(tǒng)計參數(shù),采用計算機約束準則適線與專家經(jīng)驗相結合的綜合適線方法初定;再利用設計暴雨公式參數(shù)約束8種歷時頻率曲線之間的間距,使之布局合理、不產(chǎn)生相交。這樣既保證了單站暴雨統(tǒng)計參數(shù)計算的高效率和客觀性,又融入了專家的寶貴經(jīng)驗。
單站多種歷時暴雨的適線,則要協(xié)調各頻率曲線的間距;使不同歷時的同一種統(tǒng)計參數(shù)服從“參數(shù)-歷時”關系的一般規(guī)律。適線時盡可能全部照顧各量級暴雨點據(jù),使頻率曲線通過點群中心,如無法全部照顧各量級點據(jù),則側重點通過大暴雨和中等暴雨的點群中心,并充分考慮特大暴雨的重現(xiàn)期、頻率曲線的中上部擬合、同一站點不同時段頻率曲線間的間距。
設計暴雨采用的頻率分布線型為皮爾遜Ⅲ型,包括3個統(tǒng)計參數(shù):均值、變差系數(shù)CV和偏態(tài)系數(shù)CS;均值直接計算所得、變差系數(shù)CV適線優(yōu)化而定,工程設計中,習慣上采用CS/CV取代CS。
(1)在歷時上,隨著歷時的增長,均值在增大。最大10 min點暴雨均值在10 mm~16 mm之間、20 min在14 mm~22 mm之間、60 min在22 mm~32 mm之間、3 h在30 mm~~45 mm之間、6 h在35 mm~50 mm之間、12 h在45 mm~65 mm之間、24 h在55 mm~80 mm之間、3 d在60 mm~75 mm之間。
(2)在空間分布上,各歷時均有從北往南逐漸增大的趨勢,歷時越長,這種趨勢越明顯。
(3)長歷時均值分布受山體影響特別明顯,6 h以上各歷時均值高值區(qū)出現(xiàn)在黃龍山、秦嶺脊線、米倉山和大巴山、鶻嶺和蟒嶺等山地。
(4)各歷時暴雨均值均有幾個低值區(qū),出現(xiàn)在毛烏素沙漠區(qū)、涇渭平原區(qū)、漢中盆地等區(qū)域。
(5)短歷時暴雨均值地區(qū)分布減弱,高值區(qū)減少,等值線分布稀疏。
單站變差系數(shù)CV隨歷時變化的規(guī)律表現(xiàn)為:左偏鈴形連續(xù)光滑曲線,極大值多出現(xiàn)在60 min、3 h或6 h處;或單調下降曲線;也有少數(shù)單調上升曲線,這個主要在漢中的米倉山暴雨高值區(qū)。10 minCV值在0.35~0.75之間、20 min在0.40~0.70 之間、60 min在0.40~0.60之間、3 h在0.45~0.60之間、6 h在0.45~0.60之間、12 h在0.40~0.60之間、24 h在0.40~0.55之間、3 d在0.40~0.55之間。
變差系數(shù)CV在區(qū)域分布上:與均值在地域分布上有明顯的不同,CV值渭河分界非常明顯,即渭河區(qū)域為全省高值區(qū);渭河以北CV值自北向南逐漸增大,渭河以南CV值自北向南逐漸減?。籆V值出現(xiàn)幾個低值區(qū),北洛河中游、延河中游一帶,米倉山、巴山和鶻嶺一帶。但在24 h以上歷時,渭河分界已不明顯,延河以北為高值區(qū),漢中盆地為低值區(qū)。
(1)統(tǒng)計最優(yōu)CS/CV值
采用不固定CS/CV值,對全省126~633 個雨量站實測資料系列8個歷時按統(tǒng)一準則CV和CS以經(jīng)驗頻率與P-Ⅲ型頻率曲線離差平方和最小準則,按CS/CV≥5.0、CS/CV≥4.0、CS/CV≥3.5、CS/CV≥3.0、CS/CV≥2.5、CS/CV≥2.0和CS/CV<2.0 7種情況分別統(tǒng)計,共得3233個CS/CV值,計算其各歷時的CS/CV均值為3.5。計算結果見統(tǒng)計見表1。
表1 全省最優(yōu)適線CS/CV值統(tǒng)計
(2)CS/CV取值3.5 對適線結果的影響
從633 處雨量站中選出92 處資料系列長、觀測完整(均有8個時段)、有代表性的典型雨量站,通過統(tǒng)一適線準則,使點據(jù)與P-Ⅲ型曲線擬合程度最高(即CS/CV不固定),以此計算結果與CS/CV采用3.5的頻率結算結果來比較分析,來量化其誤差大?。ò俜致剩?見表2。
從表2 可以看出,在歷時上誤差呈橫“S”型減小,即隨著歷時增長,設計值誤差逐漸減小,在3 h出現(xiàn)最小,隨后又增大,在3 d時又出現(xiàn)小值;在頻率設計值上呈單調減小,即頻率越小誤差越大。
(3)CS/CV取3.5 與3.0 設計成果分析
在一些地區(qū)或者一些學者對CS/CV取值不盡一致,比如《榆林市實用水文手冊》(2019年12月)CS/CV取值是3.0,因此,對篩選的92 站CS/CV分別取值3.5和3.0 計算其頻率,分析2種取值計算成果的差別(百分率),見表3。
從表3 可以看出,在歷時上和頻率設計值上均呈單調減小,即隨著歷時的增長,誤差越小,頻率越大(重現(xiàn)期越小)誤差越小,但各誤差均小于5%。充分說明CS/CV取值3.5與3.0設計成果較為接近,但CS/CV取值3.5時設計值較大,利于工程安全。因此,設計暴雨CS/CV取值3.5是合理的,符合《規(guī)范》要求。
表3 固定CS/CV取值3.5與3.0頻率計算結果誤差統(tǒng)計
陜西省暴雨統(tǒng)計參數(shù)特征在地域上有一定的分布規(guī)律,同一歷時均值自北向南逐漸增加,歷時越長增加越明顯。CV值和均值在地域上有明顯的不同,除24 h和3 d歷時外,CV值渭河分界非常明顯;渭河以北自北向南逐漸增大、渭河以南自北向南逐漸減??;渭河以北最小值在北洛河中游、延河中游一帶,渭河以南最小值在米倉山、大巴山和鶻嶺一帶。CS/CV均值為3.5;CS/CV取3.5和最優(yōu)準則適線取值頻率計算成果各歷時誤差在-7.2%~-2.6%之間;CS/CV取值3.5和3.0,頻率計算成果各歷時誤差在-4.2%~0.2%之間,說明CS/CV取3.5是最優(yōu)最安全的。