東坑水庫除險加固設計洪水分析

鄭 毅

(福建省浦城縣東坑水庫管理處，福建 浦城 353400)

1 工程概況

東坑水庫位于福建浦城縣管厝鄉巖步溪上游約3 km處，所在河流為南浦溪支流巖步溪，距浦城縣城約35 km。水庫壩址以上集雨面積為38.1 km2，主河道長16.32 km，河道平均坡降25.1‰，是一座以灌溉為主兼有發電、防洪、養殖等綜合利用的中型水庫[1]。工程由水庫樞紐和灌區兩大部分組成，水庫現狀總庫容1552.2萬m3，其中興利庫容1206萬m3，為年調節中型水庫。水庫現狀設計洪水標準為50 a一遇設計，1000 a一遇校核，校核洪水位為363.24 m，水庫正常蓄水位為360.00 m。水庫設計灌溉面積1.02萬畝，灌溉保證率90%。保護下游人口0.96萬人，耕地1.7萬畝。水庫于1980年6月建成開始蓄水運行，至今已運行40余年，發揮了巨大的工程效益。2015年6月8日以來，東坑水庫內連降大雨和暴雨，6月11日因大量洪水沖刷導致水庫溢洪道底部靠廠房側山體產生大面積滑坡，嚴重影響了大壩和電站的安全。為確保水庫的安全運行，充分發揮水庫的灌溉、防洪、發電等綜合效益，對東坑水庫進行應急除險加固。本文主要對除險加固中設計洪水計算和調洪計算分析，將得到的設計成果用于本次除險加固設計。

2 設計洪水分析

2.1 水文測站及暴雨資料

2.1.1 測站分布基本情況

離東坑水庫壩址15 km處有黨溪雨量站，黨溪雨量站有1964年～2012年共49年實測的6 h、24 h的暴雨成果，考慮到黨溪雨量站缺乏1 h的暴雨資料，增加城關雨量站的資料進行分析比較，同時采用2002年修正后的福建省暴雨等值線圖查得各短歷時暴雨資料進行比較，東坑水庫區域水系見圖1。

圖1 浦城水系圖

2.1.2 基本資料

(1)水庫鄰近的黨溪雨量站暴雨資料

距離東坑水庫壩址15 km處有黨溪雨量站，黨溪雨量站有1964年～2012年共49 a實測的6 h、24 h的暴雨系列資料，詳見表1。

由于黨溪雨量站缺乏1 h的暴雨資料，增加城關雨量站資料進行分析比較。城關雨量站距水庫壩址33.6 km，有1950年～2013年共64 a實測的1 h、6 h和24 h 的暴雨系列資料，詳見表1。根據2002年修正后的福建省暴雨等值線圖查得各短歷時暴雨資料，見表1。

(2)本工程防洪標準采用50 a一遇設計，1000 a一遇校核洪水標準，應根據地區水文手冊，查得此處的點面系數，然后才能由點雨量推出面雨量，見表2。

表1 水庫鄰近的雨量站及福建省暴雨等值線圖查得各短歷時暴雨資料

(a)黨溪站年最大6 h暴雨頻率曲線圖

(b)浦城城關雨量站年最大6 h暴雨頻率曲線圖

表2 各頻率面暴雨量成果表

根據表2的比較可知，暴雨等值線圖的暴雨參數均大于各雨量站的成果，因此本次復核，從偏安全的角度考慮，采用暴雨等值線的成果進行洪水復核。

2.2 設計洪水分析

2.2.1 設計洪水計算方法選擇

由于東坑水庫所在流域內無水文站，考慮到水庫的集雨面積為38.1 km2，本次計算洪水采用推理公式法、瞬時單位線法，和地區綜合法，將三者得到的成果進行分析比較，選擇更加合理的成果作為東坑水庫設計洪水成果[2-3]。

(1)推理公式法

由小溪流域特征參數F、L、J；暴雨參數Sp、n；產流參數μ；匯流參數m等，用推理公式法求得設計洪峰流量，計算公式如下：

(1)

(2)

當tc≥τ時，即全面匯流情況下：

(3)

當tc<τ時，即部份匯流情況下：

(4)

式中：Qm為地表洪峰流量，m3/s；F為流域面積，km2；Τ為流域匯流歷時，hr；tc為地表產流歷時，hr；Rτ為匯流歷時內的最大地表凈雨量，mm；Rtc為產流歷時內的地表凈雨量，mm；m為匯流參數，內陸區取m=0.039θ0.712(θ≥1.5 h)，m=0.045θ0.335(θ<1.5時)。

(2)瞬時單位線法

采用納希瞬時單位線，假設流域是由n個相同的線性水庫串聯而成，降落在流域上瞬時的單位凈雨，經過n個相同的線性水庫的調節，調蓄后在流域出口斷面處所形成的出流過程線就是相應于該凈雨的瞬時單位線。計算公式如下：

(5)

式中：Gn為伽馬函數，其值為n-1的階乘；n為線型水庫的調節次數，n=3.87F0.066J-0.151；t為匯流時間；e為自然對數的底；k為蓄泄系數，相當于流域匯流時間的參數，k=m1/n，其中：

(3)地區綜合法

為滿足無實測流量資料推求設計洪水及成果合理性檢查的需要，本次設計擬采用臨近小安下、赤橋、柳溪、渡里水文站的洪水分析成果，見表3。

表3 各水文站的洪水特征值表

根據地區綜合法求得的地區綜合線(見圖3)，計算得到東坑水庫壩址處的Q=97.6 m3/s、Cv=0.81、Cs=3.5Cv，求得壩址處的各頻率洪水成果。

圖3 流域及鄰近各站年最大流量均值Cv-F關系線

2.2.2 設計洪水推求

各頻率設計洪水成果見表4。

表4 設計洪水成果表

從表4中可以看出推理公式法和瞬時單位線法的計算結果十分接近，地區綜合法的計算成果偏大。其中推理公式法和瞬時單位線法適用于集雨面積大于30 km2，小于500 km2的流域，而地區綜合法更適用于集雨面積大于500 km2的流域，本次計算的東坑水庫壩址以上集雨面積為38.1 km2，推理公式法和瞬時單位線法更加適用，而從工程偏安全的角度考慮，選擇計算成果更大的瞬時單位線法更加合理可靠。即本次計算的東坑水庫50 a一遇設計洪峰流量為267 m3/s，1000 a一遇校核洪峰流量為424 m3/s。

2.3 合理性分析

將水庫的歷次洪水成果進行列表比較分析本次計算成果的合理性，結果見表5。

表5 歷次洪水成果比較表 單位：m3/s

表5為東坑水庫歷次洪水計算的成果，都采用了推理公式法，但本次除險加固的推理公式法中的暴雨參數采用最新版的“福建省暴雨等值線圖”，故計算結果偏大于2004年和2006年的成果(校核情況)。其中2004年安全評價最后采用水文比擬法的計算成果，選取的參證站為小安下水文站，由于小安下水文站系列僅有29年(1971年～1999年)的流量資料，存在較大的抽樣誤差且兩者的就面積相差較大，故成果存在一定的誤差。2006年除險加固最后采用的是地區綜合法的計算成果，地區綜合法分析的是一個區域內的洪水特征情況，綜合線的成果受采用的水文站的影響很大，不能很好地體現本次計算的洪水特征，通常只作為洪水合理性分析的驗證手段。本次洪水成果計算采用的推理公式法和瞬時單位線法，能夠反映流域的水文氣象、地形地貌等多因素的差別，計算公式中的暴雨參數均使用最新版本的查算圖表和計算公式，成果合理可靠，可以很好地體現水庫流域洪水特性。

2.4 洪水總量計算和洪水過程線

水庫缺乏歷年歷次暴雨洪量的實測資料，本次復核依據《福建省產匯流計算方法》，以24 h設計暴雨不扣損失計算洪水總量：

(6)

經計算，各頻率設計洪水成果見表6。

表6 設計洪水總量計算成果表(24 h)

設計洪水過程線的計算根據福建省暴雨洪水概化過程線系數概化，計算公式如下：

(7)

式中：Wp為地表洪水總量，m3；Qm為設計洪峰流量，m3/s；T=50 h該洪水過程線為50 h過程線，選取其中最大24 h洪水過程線，根據W24/(F×H24)比值進行過程線放大，成果見表7。

表7 設計洪水過程線 單位：m3/s

3 結論及建議

設計洪水成果在東坑水庫除險加固設計中十分重要，本文根據現有水文資料和成果，選取推理公式法、瞬時單位線法和地區綜合法3種方法進行洪水成果計算，最終選擇出適合該工程水文計算的瞬時單位線法洪水成果，為本次除險加固設計提供了具體的數據資料。目前，工程已完工，且運行良好。