淺談北墅水庫洪水預報方法

【摘要】結合水庫洪水預報的實際情況，對前期影響雨量Pa、凈雨R、經驗單位線法進行探討分析

【關鍵詞】北墅水庫；洪水預報；方法

1、前言

北墅水庫是一座集防洪、供水、養殖于一體的中型水庫，汛期能否進行正確的洪水調度，關系到下游萊西市、平度市、青島市、即墨市12萬人， 20萬畝耕地的人民群眾生命及財產安全，及時準確的洪水預報，是進行洪水調度的前提和基礎。

2、前期影響雨量Pa

前期影響雨量Pa是反映雨前土壤蓄水量變化的參數，其主要影響因素是表層土壤的組成成分和厚度、降雨量的大小、蒸發量的多少、土壤入滲量的快慢等，較適用的計算方法大致分兩種：

2.1 扣損法

無雨日：Pa（t+1）=Pa（t）-E

有雨不產流日：Pa（t+1）=Pa（t）+P-E

有雨產流日：Pa（t+1）=Pa（t）+P-R-E≤Im

式是：Pa（t），Pa（t+1）—第t天，第（t+1）天的Pa值。

P—流域平均降雨量；R—凈雨量；E—第t日的陸面蒸發量，Im—土壤最大持水量。

北墅水庫水文站有E601蒸發皿觀測資料，就是陸面蒸發量。據多年資料統計，6月份平均日蒸發量為4.7毫米，7月份為4.5毫米，8月份為3.9毫米，9月份為3.6毫米，逐月逞遞減趨勢。

有雨產流日的公式中P-R-E指土壤入滲的部分降雨量。嚴格的講P-R-E的部分降雨量不能全部入滲土壤，總有一部分降雨量被坑、洼、庫、塘攔蓄，但計入攔蓄部分，使Pa的計算更接近實際，當土壤達到飽和，即土壤含水量達到Im，可以認為后期降雨量不再入滲土壤，全為產流凈雨，這時的Pa就達到最大值Im。

2.2 系數法

Pa（t+1）=K[Pa（t）+P]

但須控制Pa（t+1）≤Im

式中，K—流域土壤蓄水的日消退系數。北墅水庫預報時采用此法，精確與否取決于選用的系數K，K值取用范圍大致在0.85—1.0，當降雨小于5毫米時取K=0.9，當降雨在5—15毫米時取K=0.95，當降雨大于15毫米時取K=1.0。

3、凈雨量R

凈雨量是指毛降雨扣除損失后的產流雨量，損失量包括蒸發、植被截留、填洼、土壤入滲等，其中土壤入滲對凈雨量R影響較大，R的計算方法大致分為三種：

3.1 徑流系數法：R=ap，a—徑流系數，選擇水庫流域若干次降雨徑流資料，查出毛降雨量，由徑流量反推出產流凈雨，點繪P—R曲線，預報時查曲線即可。

3.2 降雨徑流相關法

水庫汛期控制運用方案中繪有降雨徑流關系P+Pa—R曲線，計算出P、Pa，然后查關系曲線得R值，并適當調整R值。

P+Pa—R關系曲線是按照一種典型雨型設計出來的，而實際發生的降雨徑流過程一次一個樣，對查出的凈雨R值調整多少憑經驗調整。例如2001年8月2日，北墅水庫發生的降雨徑流過程，流域P=51.6毫米，Pa=42毫米，查P+Pa—R關系曲線R值為30毫米，北墅水庫預報時將R值調整為4.7毫米，調整系數為0.156，實際發生徑流深為4.0毫米，預報與實際基本相符。

3.3 初損入滲法

P=P+Pa-Im=P-（Im-Pa）

蓄滿產流理論認為，一次降雨的損失量由流域的最大持水量Im減去降雨開始時的土壤含水量Pa求得，從毛降雨中扣除損失量即為凈雨R值，也就是形成的總徑流深。

2003年8月5日北墅水庫發生的降雨徑流過程，流域平均雨量84.7毫米，推算Pa為44.3毫米，Im為98毫米，計算R=84.7+44.3-98=31毫米，考慮降雨強度不同，雨時蒸發影響，取R值為29毫米，與實際發生的徑流深26.5毫米相符。

Im實際取用中考慮了流域植被、蓄水情況的影響，使R值的取用更接近實際，因此初損入滲法，理論與實際能很好的結合，有規律可循，容易掌握，是值得推廣的方法。

4、匯流的計算

預報匯流過程，我們一般采用以下法：

4.1 公式法：W=0.1RF

Qm=0.948F0.84J0.32R1.15tc-0.46

T=5W/QmW—洪水總量（萬米3） Qm—洪峰流量（米3/秒）

T—洪水歷時（小時） F—產流流域面積（平方公里）

Tc—凈雨歷時（小時） J—干流坡度

用公式計算，W、Qm、T式中各參數的選用比較固定，很難代表不同次降雨徑流過程的特殊規律。尤其是北墅水庫流域屬山丘區，洪水來得急并且快，還受萊洲橡皮壩的影響，加之近幾年連續干旱，各流域內城鄉競相挖水攔蓄，破壞了原流域的規律。因此，使用公式法計算的W、Qm、T都要大幅度調整。所以用公式法做為大致估算是可以的，若要提高精度，一般不易做到。

4.2 經驗單位線法

所謂單位線是指單位流域（100平方公里）在單位時段（1小時或2小時）內由單位凈雨深（10毫米）產生的入庫流量過程線。它是由實際的降雨徑流資料反推得出的，沒有抽象理論公式，所以有經驗之稱。

經驗單位線的分析查找，要經過一番系統的資料收集、數據分析。根據本流域發生的降雨徑流過程，找出實測流量過程，實際凈雨過程，分析出經驗單位線。

經驗單位線也可通過試算法得出：（1）用瞬時單位線法計算參數M，M=0.196×F0.33J-0.27tc0.17查瞬時單位線表，找出初步經驗單位線。（2）查出兩個或三個時段凈雨產生的徑流過程。

由實際凈雨過程，徑流過程初步經驗單位線進行試算修正，即可找出較為理想的經驗單位線。

單位線可以有多條，以凈雨的大小、歷時進行分列，當發生降雨時，選用情形較近的一條使用。例如北墅水庫分析得出時段為2小時的經驗單位線3條，分列如下：當R≤50mm，1-7時段入庫流量分別是33、61、28、10、2、1、0.1米3/秒 ；當50100mm，1-7時段入庫流量分別是36、64、29、10、2、1、0米3/秒。

預報時還可以將流域劃分成幾個區間，將區間洪水錯時段疊加。北墅水庫有水庫站、廟埠河站，郭家店站，將流域分成兩半流域，當整個流域發生不同步、不等量的降雨時，分區面疊加預報徑流過程要精確的多。例如在預報2003年7月31日洪水過程時，用經驗單位線法分兩個區間疊加預報，預報出洪水總量246萬方，洪峰流量120立方米/秒，洪峰時間為8月2日2時，而實際洪水總量為230萬方，洪峰流量109立方米每秒，洪峰出現時間為8月2日1時30分，精確度在90%以上。

參考文獻：

[1]《 工程水文學 》；

[2]《水庫管理指南》；

[3]《防汛抗旱專業干部培訓教材》