王 振

(新疆水利水電勘測(cè)設(shè)計(jì)研究院，烏魯木齊 830000)

潰壩洪水具有突發(fā)性強(qiáng)、歷時(shí)短、峰高量大、難以預(yù)測(cè)的特點(diǎn)，潰壩洪水突發(fā)事件是破壞力強(qiáng)、危害范圍廣、損失嚴(yán)重的突發(fā)事件。新疆洛浦縣布爾庫(kù)木水庫(kù)位于灌區(qū)內(nèi)部，為均質(zhì)土壩小型水庫(kù)，水庫(kù)下游居民點(diǎn)、耕地密布，水庫(kù)修建年代早、設(shè)計(jì)標(biāo)準(zhǔn)低、防洪能力差、管理不到位，一旦發(fā)生潰壩，將對(duì)下游造成極為嚴(yán)重的生命、財(cái)產(chǎn)損失。截至目前，國(guó)內(nèi)外專家學(xué)者對(duì)土石壩水庫(kù)潰壩洪水?dāng)?shù)值模擬方面取得了豐碩成果[1-6]。本文選用BREACH模型和簡(jiǎn)化計(jì)算法模擬布爾庫(kù)木水庫(kù)潰壩洪水，模擬計(jì)算成果可為水庫(kù)突發(fā)洪水事件應(yīng)急預(yù)案制定和風(fēng)險(xiǎn)評(píng)估提供技術(shù)支持。因此，開(kāi)展水庫(kù)潰壩洪水研究具有重要的現(xiàn)實(shí)意義。

1 工程概況

布爾庫(kù)木水庫(kù)位于新疆和田地區(qū)洛浦縣多魯鄉(xiāng)境內(nèi)，水庫(kù)利用天然洼地筑壩而成，是一座以灌溉為主的引水注入式平原水庫(kù)，水庫(kù)水源為玉龍喀什河河水。壩址距離洛浦縣多魯鄉(xiāng)政府約2.5 km，地理坐標(biāo)為E80°16′08″，N37°5′45″。

布爾庫(kù)木水庫(kù)主要由大壩、放水涵洞、引水渠和放水渠組成。水庫(kù)正常蓄水位1 324.98 m，正常蓄水位相應(yīng)庫(kù)容150×104m3；死水位1 323.00 m，死庫(kù)容10×104m3；調(diào)節(jié)庫(kù)容140×104m3。水庫(kù)大壩為均質(zhì)土壩，全長(zhǎng)2 300 m，壩頂高程1 326.5 m，最大壩高4.18 m，工程為小(Ⅰ)型Ⅳ等工程，水庫(kù)大壩為4級(jí)，地震設(shè)防烈度為Ⅶ度。

2 BREACH模型模擬成果

BREACH模型是由美國(guó)國(guó)家氣象局Fread研制的一種有關(guān)土壩失事的模型，用來(lái)預(yù)測(cè)由管道失事或洪水漫溢引起的土壩/堆石壩潰口出流過(guò)程和決口大小。

2.1 水量平衡

BREACH模型采用水量平衡關(guān)系模擬水庫(kù)水位的變化，公式如下：

(1)

其中：ΔH為計(jì)算時(shí)段內(nèi)水庫(kù)水位變化；Δt為計(jì)算時(shí)段，較小值；Sa為水面面積；Qi為水庫(kù)入庫(kù)流量；Qb為流經(jīng)決口的寬頂堰流量；Qsp為溢洪道出流流量；Qo為流經(jīng)壩頂寬頂堰流量。

2.2 潰口流量計(jì)算公式

漫頂潰口流量公式采用寬頂堰公式計(jì)算：

Qb=3Ab(H-hb)0.5(寬頂堰水流)

(2)

如果決口因管道失事所致，用短管、孔流方程代替寬頂堰水流方程，則：

Qb=Ab[2g(H-hb)/(1+fL/D)]0.5(孔流)

(3)

其中：Ab為水流所流過(guò)的堰或潰口的面積；hb為壩下游面決口底部高程；f為達(dá)西摩擦系數(shù)，取決于D50粒徑；L為管道長(zhǎng)度；D為管道直徑或?qū)挾取?/p>

2.3 決口寬度

決口最初可認(rèn)為是寬度為B0的矩形，B0=BrY，Y是決口入口處的臨界水深，Y=2/3(H-hb)。Br為無(wú)量綱因子，漫溢決口時(shí)Br為2，管道失事決口時(shí)Br為1。

有效決口深度超過(guò)臨界水深dc發(fā)生坍塌，決口由原來(lái)的矩形可以變?yōu)樘菪巍?/p>

dc=4Ccosφsinθ′/[γ-γcos(θ′-φ)]

(4)

其中：C為土壤內(nèi)聚性；Φ為壩料內(nèi)部摩擦角；γ為土壤容重；θ′為壩坡坡腳。

坍塌引起的角度θ′減小到新角度，公式為：

θ′=(θ′+φ)/2

(5)

2.4 決口侵蝕

除側(cè)面坍塌之外，假定決口侵蝕沿底部和側(cè)面同時(shí)進(jìn)行。一旦決口沖刷到規(guī)定的壩底，侵蝕只在決口側(cè)面繼續(xù)進(jìn)行。決口侵蝕速度取決于水流輸送侵蝕材料的能力。泥沙輸送關(guān)系用下式表示：

(6)

其中：Qs為壩料沖蝕率；D90、D50、D30分別為小于或等于某粒徑的沙重占總沙重90%、50%和30%的泥沙粒徑；D為用曼寧方程求得的任一瞬時(shí)決口水流的水力深度；S為決口底部坡降,假定與壩下游面平行；P為決口總周長(zhǎng)；n為曼寧系數(shù)；τc為產(chǎn)生侵蝕所必須超過(guò)的Shield臨界剪切力。

在計(jì)算時(shí)間內(nèi)，模型假設(shè)決口底部和兩側(cè)沖蝕率相等，底部和側(cè)面產(chǎn)生的增量相等，計(jì)算公式如下：

ΔHc=QsΔt/PL(1-P)

(7)

其中：L為穿過(guò)大壩的決口長(zhǎng)度；P為決口材料的孔隙度。

2.5 潰壩流量過(guò)程線

本次采用BREACH模型，對(duì)布爾庫(kù)木水庫(kù)潰壩洪水模擬，計(jì)算結(jié)束條件為大壩水位與潰口底部高程差小于0.1 m、潰口基本不再出流。經(jīng)計(jì)算，壩址處最大潰壩流923 m3/s。

圖1 布爾庫(kù)木水庫(kù)BREACH模型模擬潰壩洪水

3 簡(jiǎn)化計(jì)算法模擬結(jié)果

簡(jiǎn)化計(jì)算法為保守地利用地形與水文過(guò)程線資料和憑經(jīng)驗(yàn)確定的潰決流量表達(dá)式，適用于壩高小于等于12 m的均質(zhì)土壩[7]。

壩址處的最大潰壩流量Qm可以按下式計(jì)算：

Qm=2.5FV0.76H0.1

(8)

式中：Qm為最大潰壩流量，m3/s；F為簡(jiǎn)化評(píng)估特征參數(shù)，F(xiàn)=1.3；V為下泄水的總體積，103m3；H為最大水深，m。

布爾庫(kù)木水庫(kù)大壩為均質(zhì)土壩，最大壩高4.18 m。從大壩壩高和材料類型方面考慮，可采用簡(jiǎn)化計(jì)算公式計(jì)算最大潰壩流量。經(jīng)計(jì)算，布爾庫(kù)木水庫(kù)大壩壩址處潰壩洪水洪峰為955 m3/s。

4 兩種模擬的分析比較

1) BREACH模型模擬壩址最大潰壩流量923 m3/s，簡(jiǎn)化計(jì)算法計(jì)算最大潰壩流量955 m3/s，兩者相差較小，僅3.5%。

2) BREACH模型模擬結(jié)果為潰壩洪水過(guò)程線，為完整的潰決流量過(guò)程，可分析最大潰壩流量出現(xiàn)時(shí)間。但需綜合水利學(xué)、泥沙、土力學(xué)等學(xué)科知識(shí)，參數(shù)多，計(jì)算過(guò)程復(fù)雜，在資料具備條件下，應(yīng)盡量采用BREACH模型。

3) 簡(jiǎn)化計(jì)算公式模型簡(jiǎn)單，未考慮實(shí)際潰壩機(jī)理，僅能計(jì)算壩址最大潰壩流量。但所需計(jì)算參數(shù)少，計(jì)算簡(jiǎn)單，在資料缺乏條件下，可采用該方法進(jìn)行初步計(jì)算及突發(fā)情況應(yīng)急估算。

5 結(jié) 論

BREACH模型和簡(jiǎn)化計(jì)算法均能進(jìn)行布爾庫(kù)木水庫(kù)潰壩模擬，模擬成果可為制定水庫(kù)突發(fā)洪水事件應(yīng)急預(yù)案和風(fēng)險(xiǎn)評(píng)估提供技術(shù)支持，對(duì)其他均質(zhì)土壩小型水庫(kù)潰壩洪水模擬具有一定指導(dǎo)意義，但應(yīng)根據(jù)水庫(kù)資料情況合理選擇計(jì)算模型。