王 云，張 鑫，王文亞，蘇夏羿

(西北農(nóng)林科技大學水利與建筑工程學院，陜西 楊凌 712100)

0 引 言

近年來，我國山洪災(zāi)害頻發(fā)，造成了大量人員傷亡和財產(chǎn)損失，已成為我國防災(zāi)減災(zāi)工作中亟待解決的突出問題。許多國家在山洪預(yù)警指標分析方面開展了大量的研究工作，美國研發(fā)了在歐美國家和地區(qū)廣泛應(yīng)用的山洪預(yù)警指標確定方法(FFG)及平臺[1,2]，該方法在東南亞和南美等地區(qū)也逐漸得到推廣；澳大利亞經(jīng)充分調(diào)研世界各國山洪預(yù)警成果后建立了自己的預(yù)警平臺[3]；日本提出了土壤雨量指數(shù)法、實效雨量法、多重判別分析統(tǒng)計法等方法來確定臨界雨量[4]。目前國內(nèi)針對山洪災(zāi)害預(yù)警指標研究采用的主要方法包括統(tǒng)計歸納法和流域水文模型法，如實測雨量統(tǒng)計法、單站臨界雨量法、P-Ⅲ頻率分析法、新安江模型法[5]和降雨徑流API模型法[6]等。

但是統(tǒng)計歸納法和流域水文模型法對下墊面和水文氣象資料要求較高，使得其很難在基層預(yù)警指標分析中廣泛應(yīng)用[8]。自新中國成立以來，雨量站基本建在大江大河以及水庫上，對于山區(qū)村落雨量站分布稀少且年份資料短，部分地區(qū)甚至沒有[7]。水位流量反推法是山洪災(zāi)害預(yù)警指標確定的主要方法之一[9]。水位流量反推法的前提是假定暴雨與洪水同頻率，在這樣的假設(shè)條件下，所需資料僅包括預(yù)警地點斷面地形資料、設(shè)計暴雨和設(shè)計洪水資料等[10,11]，可應(yīng)用在資料匱乏的山區(qū)小流域。本文通過對比分析研究，獲得適宜于該地區(qū)的設(shè)計洪水計算方法，基于計算所得設(shè)計洪水，采用考慮不同土壤含水量情況下的水位反推法計算典型流域的臨界雨量，為陜南地區(qū)的山洪災(zāi)害防治提供技術(shù)支撐。

1 研究方法

1.1 水位反推法

水位反推法是基于山洪災(zāi)害發(fā)生的物理機制，假設(shè)洪水與降雨同頻率，應(yīng)用水文、水力學方法推求臨界水位對應(yīng)的臨界流量，再計算臨界流量對應(yīng)的臨界頻率，暴雨頻率曲線中臨界頻率對應(yīng)的雨量即為臨界雨量。因假定洪水與降雨同頻率，故該方法只適用于無資料小流域地區(qū)。具體步驟如下：

(1)繪制各時段降雨量與頻率關(guān)系曲線圖。

(2)進行洪水計算，確定各設(shè)計頻率P下各時段設(shè)計洪水洪峰流量，并繪制頻率與設(shè)計洪水洪峰流量的關(guān)系曲線圖P～Qmax,p。

(3)確定控制斷面水位流量關(guān)系，由成災(zāi)水位查水位流量關(guān)系曲線得到相應(yīng)預(yù)警流量Ql，從P～Qmax,p關(guān)系曲線圖上查得對應(yīng)的頻率值Pl。根據(jù)Pl查設(shè)計暴雨頻率曲線圖查得各時段臨界雨量，即立即轉(zhuǎn)移預(yù)警雨量。

1.2 設(shè)計洪水推求方法

1.2.1 經(jīng)驗公式法

陜南地區(qū)計算洪峰流量的經(jīng)驗公式基本形式為[12]：

Qp=CpFn

(1)

式中：Qp為頻率為p的設(shè)計洪峰流量，m3/s；F為流域面積，km2；Cp、n為頻率為p的經(jīng)驗參數(shù)(由表查得)。

1.2.2 推理公式法

(2)

式中：hτ、hR分別為τ所對應(yīng)的最大凈雨量、總產(chǎn)流量；

繪制Qt～t和Qτ～t兩組曲線，兩曲線交點所對應(yīng)的坐標即為所求的設(shè)計洪峰流量和匯流時間。

1.2.3 瞬時單位線法

瞬時單位線就是歷時趨于無限小時流域內(nèi)單位凈雨在出流斷面所形成的徑流過程線。實際應(yīng)用時，用S曲線法將瞬時單位線轉(zhuǎn)換為時段單位線，其基本形式如下：

(4)

式中：U(0，t)為t時刻瞬時單位線縱高；n為相當于線性水庫的個數(shù)或調(diào)節(jié)次數(shù)；K為水庫型線性蓄泄方程的匯流歷時；T為時間；Δt為時段長，h；S(t-Δt)為移后Δt時段的S(t)。

如果i≤iK，則i=IK代入求m1;若i>iK，則直接用i代入求m1。

2 應(yīng)用實例

2.1 研究區(qū)概況

選取勉縣將臺村、廟坪村和小砭河村為典型防災(zāi)對象，如圖1所示。勉縣為典型的內(nèi)陸性季風氣候，屬亞濕潤區(qū)。氣候溫暖濕潤，四季分明，雨熱同季，時空分布不均，垂直差異較大，土壤類型主要為棕壤與黃棕壤。勉縣地區(qū)常降雷暴雨，降雨歷時短，籠罩面積小。本文基于ArcGIS水文分析模塊生成流域河網(wǎng)，并量算出所選典型小流域的集水面積、主溝道長度、比降等特征參數(shù)，如表1所示。

圖1 勉縣水系圖Fig.1 Drainage map of Mianxian

典型對象L/mF/km2J/×10-3n將臺村16．0842．2347．40．035廟坪村24．0692．8557．00．035小砭河村31．99129．7226．50．030

2.2 設(shè)計暴雨頻率曲線

根據(jù)小流域面積大小確定設(shè)計點暴雨歷時分別為1、3、6、12 h，結(jié)合《漢中地區(qū)水文手冊》和《陜西省暴雨圖集》由圖表查算法求得設(shè)計暴雨。并繪制各時段設(shè)計暴雨頻率曲線，如圖2所示。

2.3 不同方法設(shè)計洪水推求

2.3.1 經(jīng)驗公式法

由《漢中地區(qū)水文手冊》可知陜南地區(qū)經(jīng)驗公式各參數(shù)Cp、n取值如表2所示，根據(jù)經(jīng)驗公式得設(shè)計洪峰流量，結(jié)果見表3。

圖2 典型防災(zāi)對象設(shè)計暴雨頻率曲線Fig.2 Design storm frequency curve of the typical areas of disaster prevention

表2 經(jīng)驗公式參數(shù)Tab.2 Parameter of empirical formula

表3 經(jīng)驗公式法設(shè)計洪峰流量 m3/s

2.3.2 推理公式法

陜南地區(qū)屬于濕潤地區(qū)，以“蓄滿產(chǎn)流”為主，用初損一次扣除法計算凈雨。采用流域最大蓄水量Im的2/3為設(shè)計情況下的前期影響雨量，即Pa=2/3Im。評價對象均位于漢江北，故取Im=82 mm，潛流量按產(chǎn)流總量的20%計算[12]。

根據(jù)公式(2)計算得設(shè)計洪峰流量如表4所示。

表4 推理公式法設(shè)計洪峰流量 m3/s

2.3.3 瞬時單位線法

由典型小流域面積與比降得n、k，查S(t)曲線表[12]，計算得各評價對象的設(shè)計洪峰流量表5所示。

表5 瞬時單位線法設(shè)計洪峰流量 m3/s

由表3～5可知三種方法計算結(jié)果存在差異，經(jīng)驗公式法計算結(jié)果偏小，瞬時單位線法計算結(jié)果偏大。經(jīng)驗公式法與推理公式法計算結(jié)果相對誤差多數(shù)小于20%，推理公式法與瞬時單位線法計算結(jié)果相對誤差多數(shù)小于15%。經(jīng)綜合分析可知經(jīng)驗公式法僅將流域面積作為變量，考慮影響因素過少，計算結(jié)果可靠性較差；瞬時單位線法適用于面積較大的流域[12]，而山區(qū)洪水流域面積偏小，且無實測系列資料，故選用推理公式法計算典型防災(zāi)對象的設(shè)計洪水。

2.4 不同土壤含水量的設(shè)計洪水

因土壤含水量對預(yù)警結(jié)果影響較大，故采用雨前土壤較干、一般和較濕三種典型情況進行預(yù)警指標的確定，三種情況下土壤含水量分別為0.2Im、2/3Im和0.8Im。根據(jù)推理公式法及各項參數(shù)可計算得典型防災(zāi)對象不同土壤含水量下的設(shè)計洪峰流量，如表6所示。并繪制設(shè)計洪水洪峰流量頻率曲線P～Omax,p，如圖3所示。

2.5 臨界雨量

根據(jù)實地調(diào)查情況，將住戶宅基地高程按河道比降轉(zhuǎn)換至控制斷面，并選取最低宅基地高程作為成災(zāi)水位。根據(jù)各防災(zāi)對象斷面具體形狀由曼寧公式確定水位流量關(guān)系并繪制水位流量關(guān)系曲線，如圖4所示。由成災(zāi)水位查水位流量關(guān)系曲線得到相應(yīng)預(yù)警流量，結(jié)果見表7。

圖3 設(shè)計洪峰流量頻率曲線P～Qmax,pFig.3 Design peak flow frequency curve，P～Qmax,p

圖4 水位流量關(guān)系曲線Fag.4 Relation curve between water level and discharge

表6 不同土壤含水量下的設(shè)計洪峰流量 m3/s

由水位反推法計算三種土壤含水量情況下臨界雨量結(jié)果如表7所示。因?qū)⒆畹驼馗叱套鳛槌蔀?zāi)水位，故所得臨界雨量即為山洪預(yù)警時立即轉(zhuǎn)移預(yù)警雨量。

廟坪村與小砭河村均位于秦嶺山區(qū)，下墊面條件相似，且小流域集水面積相近，臨界雨量計算結(jié)果符合實際情況，故由水位反推法得到的臨界雨量預(yù)警指標較為可靠。

經(jīng)分析可知，相同時段下，臨界雨量隨著土壤含水量的增大而減??；土壤含水量對臨界雨量的影響程度隨著時段的增長而增大。

表7 典型對象臨界雨量值 mm

3 結(jié) 論

本研究采用水位反推法計算陜南地區(qū)無資料典型防災(zāi)對象在不同土壤含水量情況下的臨界雨量，得到如下結(jié)論：

(1)經(jīng)驗公式法、推理公式法和瞬時單位線法三種方法對于計算資料缺乏地區(qū)的設(shè)計洪水均有一定的適用性，且各有其優(yōu)缺性。推理公式法在陜南地區(qū)山區(qū)小流域適用性較好，計算結(jié)果更符合實際。

(2)水位反推法計算原理簡單，對水文氣象資料要求較低，在陜南地區(qū)資料匱乏小流域有一定的適用性。

(3)三種土壤含水量對應(yīng)的臨界雨量值 相差甚大，故在預(yù)警時必須考慮前期土壤含水量的影響。臨界雨量值隨著前期土壤含水量的增大而減小。前期土壤含水量對臨界雨量的影響隨著時段的增長而增大。

(4)水位反推法是基于暴雨、洪水同頻率的假定，而該假定在某些情況下并不成立，因此在應(yīng)用時，還需要結(jié)合流域水文模型法等其他方法，經(jīng)過對比分析來確定最終的臨界雨量成果。

□

[1] Carpenter T M, Georgakakos K P. GIS-based procedures in support of flash flood guidance (IIHR Report No. 366,lowaInstitute of Hydraulic Research, The University of Iowa，Iowa City)[R].1993.

[2] U.S.Department of Commerce，National Oceanic and Atmospheric Administration， National Weather Service. Modernized Areal FlashFlood Guidance[R].1992.

[3] Prasantha Hapuarachchi H A，Wang Q J. A review of methods andsystems available for flash flood forecasting[R]. 2008.

[4] Japan National Institute for Land and Infrastructure Management.Operating methods of Critical Rainfall for Warning and Evacuationfrom Sediment-Related Disasters[R]. 2001.

[5] 葉金印，李致家，常 露.基于動態(tài)臨界雨量的山洪預(yù)警方法研究與應(yīng)用[J].氣象,2014,(1):101-107.

[6] 陳瑜彬，楊文發(fā)，許銀山.不同土壤含水量的動態(tài)臨界雨量擬定方法研究[J].人民長江,2015,(12):21-26.

[7] 王新宏,唐永鵬,張美洋,等.基于水位/流量反推法的資料匱乏地區(qū)臨界雨量研究[J].水資源與水工程學報,2016,27(4):125-128.

[8] 盧燕宇，謝五三，田 紅.基于水文模型與統(tǒng)計方法的中小河流致洪臨界雨量分析.自然災(zāi)害學報,2016,25(3);38-47.

[9] 李昌志,孫東亞.山洪災(zāi)害預(yù)警指標確定方法[J].中國水利,2012,(9):54-56.

[10] 程衛(wèi)帥.山洪災(zāi)害臨界雨量研究綜述[J].水科學進展,2013,24(6):901-908.

[11] 陳桂亞等. 山洪災(zāi)害臨界雨量分析計算方法研究 [J]. 水資源研究, 2005,36(12):36-38.

[12] 陜西省漢中地區(qū)水電局. 漢中地區(qū)實用水文手冊[Z].1988.

[13] 李克先.基于推理公式的山洪預(yù)警雨量計算方法研究[J].水文,2016,36(1):84-87.