朗鄉河流域2013.8.12暴雨洪水計算分析

李金楠

摘 要：2013年8月12日，湯旺河右岸二級支流朗鄉河遭受短歷時強降雨災害，造成山災害，導致居民受災，耕地、牲畜、人民財產受損。因該流域屬小面積山溪性河流且無水文控制斷面，災后洪水調查分析采用實地過水斷面測量和水文比擬以及水文圖集瞬時單位線分析計算等方法對該流域洪水洪峰流量進行了還原計算，得出了流域暴雨和洪水的出現頻率量級，并對計算方法和結果以及成災原因進行了分析。

關鍵詞：暴雨洪水；洪峰流量；洪水頻率；暴雨頻率；流域匯流單位線

1 流域及洪水概況

朗鄉河屬湯旺河右岸二級支流，河源經度128°54′55.5″，緯度46°48′50.9″。河口經度128°52′24.7″，緯度46°57′14.6″。河流長度21.5km，流域面積160km2，河道比降11.6‰，2013年8月12日18時至24時該流域發生強降雨，6小時內降雨總量達130mm以上，河流水位迅速起漲，致使該地區形成山洪，12時后洪水漸漸消退。流域位置見圖1，因流域內無流量控制站，本次洪水沒有實測洪水過程，洪峰流量及出現頻率采用理論和調查方法推算。

2 成災原因分析

（1）本流域山洪災害3小時臨界雨量為42毫米，英山站3小時降雨為142.8毫米，大大超過山洪災害臨界雨量，歷時短、強度大是形成山洪的直接原因。（2）本河流為山區性河流，來水陡漲陡落，來勢猛、突發性強。（3）本次洪水屬于歷史上特大洪水，洪水量級已經超過現有工程防洪標準。

3 流域洪峰流量計算方法介紹

3.1 相鄰流域水文比擬法

水文比擬法是指相鄰或流域自然地理特性相似的流域之間，無控制斷面流域借用有控制斷面流域的洪水預報方案，按流域面積或其他要素的倍比關系將計算結果進行移用的一種水文分析計算方法，一般常用在流域徑流或洪峰流量的估算方面。該方法的基本依據就是認為相似流域的降雨徑流特性類似，產流規律也類似，主要差別僅在流域面積不同上。本次分析計算就是將朗鄉河相鄰的位于湯旺河右岸二級支流永翠河流域新曙光站洪水預報方案用于朗鄉河洪水預報。兩流域自然地理特征類似、流域面積、流域產匯流規律也基本相同。

3.2 流域瞬時單位線法

采用黑龍江省水文圖集中概化瞬時單位過程線法，該方法主要用來推求工程設計中針對無資料流域的設計洪水。其各種參數和備選流域匯流單位線都是通過歷史洪水整體規律概化分析而來，使用比較方便，各種參數可由圖集直接查得，基本要求事先確定降雨扣除損失后的產流凈雨量。實際工作中該凈雨量一般按百分比預估或由相鄰流域產流方案中的流域降雨～徑流相關圖查得。

3.3 曼寧公式法

該方法需在對河道進行實地洪水調查的基礎上進行，需要測量河道過水斷面，確定最高洪痕水位，測量河道縱斷面水位線洪痕點來確定斷面處洪水比降，然后利用曼寧公式計算出斷面平均流速，進而求出洪峰流量。該方法外業工作量較大，成果可信度一般較高。

4 流域降雨、洪水計算分析

4.1 降雨分析

流域主要降雨始于8月12日18時，至12日21時各站降雨均超100mm，降雨強度較大，之后降雨強度相對減弱，21日24時降雨基本停止，各站6h累計雨量均在130mm以上。流域平均雨量達100mm。

可見，英山最大6小時降雨152.6毫米，其中最大3小時降雨142.8毫米，占6小時雨量的94%；朗鄉站同時段6小時降雨132.5毫米。暴雨中心位于流域中上游，暴雨走向自上游向下游移動。

4.2 洪峰流量及出現頻率計算分析

相鄰流域水文比擬法借用新曙光站洪水預報匯流單位線成果，利用朗鄉河流域降雨量在新曙光站進行洪峰流量預報，然后將計算成果按面積倍比法來移用至梅花河和馬場河，成果移用時計算公式為：Q朗鄉河=（F朗鄉河/F新曙光）0.67*Q新曙光。借用流域預報時所需雨量直接使用朗鄉河流域實際雨量，新曙光站假定徑流深由朗鄉河流域實際降雨利用新曙光站降雨徑流相關圖查得，借用單位線選擇新曙光站歷史洪水單位線中與本次暴雨雨型相似單位線，利用查得的徑流深逐段乘所選單位線流量，即得出洪水過程和洪峰流量。

水文圖集瞬時單位線法采用降雨歷時12h，徑流深由伊春站產流方案求得，圖集所需匯流參數計算公式為：m1=B*（a）-0.31，B=C3*F0.16*J-0.35，上式中a為徑流深和降雨歷時的比值，當有效降雨歷時小雨24小時時按24小時計算。C3為匯流參數，由圖集中圖36查得。參數N=0.2*m10.74，參數k=m1/n，F為流域面積，J為河道平均坡降。利用參數n、k即可在圖集中查得S（t）曲線，轉換該曲線后即可求出相應時段洪水過程和洪峰流量。

曼寧公式法根據流域調查數據，采用公式V=（1/n）R（2/3）I（1/2），式中V為斷面平均流速（m/s），n為河床糙率，R為平均水深（m），I為水面比降。斷面洪峰流量由公式Q=VF求出，Q為流量（m3/s），V為斷面平均流速（m/s），F為斷面面積（m2）。因馬場河面積較小，本次洪水所占比重較小，僅以梅花河為例來計算分析。

以上方法計算出的洪峰流量頻率根據水文圖集來確定，計算公式為：Qmp=（kp/k5%）*Cp*F0.67。式中Qmp為流域不同頻率設計洪峰流量，kp為不同頻率模比系數，k5%為5%頻率的模比系數，Cp為最大洪水參數，F為流域面積。查圖，Cv=1.1，Cs=2.25Cv，Cp=7，計算不同頻率的最大流量見下表2、3。由表可見402立方米/秒的洪峰流量接近200年一遇。

5 結束語

經分析可知，水文比擬法與曼寧公式計算法計算結果比較接近，瞬時單位線法計算結果相對偏小，因水文比擬法和曼寧公式法所選要素都與流域實況比較接近，計算結果可信程度高些。水文圖集瞬時單位線法因各項計算參數概化程度較大，對面積較小流域實用性可能偏低，因此綜合分析確定朗鄉河流域洪峰流量為360m3/s，頻率為百年一遇洪水標準。目前對無流量測站控制的流域進行洪水預報還缺少比較權威有效的方法，文章僅以最基本的幾種方法來進行了實例分析，方法觀點難免有偏差之處，僅供參考交流。

參考文獻

[1]黑龍江省水利廳.黑龍江省水文圖集[M].哈爾濱：黑龍江省水利廳，1996.

[2]葛守西.現代洪水預報技術[M].北京：中國水利水電出版社，1999：125-145.

[3]水利部.水文情報預報規范[M].北京：中國水利水電出版社，2000：9-11.

[4]揚州水利學校.水文預報[M].中國水利水電出版社，1983：23-40.