磨子潭水庫(kù)潰口洪水分析

1、    磨子潭水庫(kù)潰口洪水分析    蔡黎明+麻天為摘 要：簡(jiǎn)要介紹了磨子潭水庫(kù)的基本情況，遵循“最不利，最可能”的原則，提出了磨子潭水庫(kù)可能潰壩的4種典型方案，并根據(jù)這些方案分析了磨子潭水庫(kù)潰口洪水的流量過(guò)程線，從而為潰壩洪水演進(jìn)模擬入流邊界的設(shè)置提供必要的參考。關(guān)鍵詞：磨子潭水庫(kù)；潰壩；潰口；洪水分析：tv122+.4 ：a doi：10.15913/ki.kjycx.2016.05.018在實(shí)際工作中，假設(shè)混凝土壩的潰決方式和潰口形態(tài)，根據(jù)水庫(kù)運(yùn)行水位和不同上游來(lái)流量值進(jìn)行組合計(jì)算，比較和分析了潰口最大流量、潰口流量過(guò)程線。潰口發(fā)展和潰口流量的計(jì)算

2、、模擬是研究潰壩問(wèn)題的基礎(chǔ)。潰口流量過(guò)程線作為潰壩洪水波演進(jìn)的入流邊界，其計(jì)算精確度關(guān)系著潰壩洪水波的演進(jìn)效果，而且還會(huì)影響對(duì)潰壩風(fēng)險(xiǎn)的分析和地區(qū)應(yīng)急預(yù)案的制訂。1 工程概況磨子潭水庫(kù)位于淠河上游安徽省霍山縣境內(nèi)，水庫(kù)下游25 km是佛子嶺水庫(kù)壩址，它與佛子嶺水庫(kù)形成串聯(lián)式的梯級(jí)樞紐，是一座以防洪為主，集灌溉、發(fā)電等功能為一體的綜合水利樞紐。其控制流域面積570 km2，總庫(kù)容3.47×108 m3。樞紐工程主要包括攔河壩、溢洪道、泄洪隧洞和發(fā)電廠。攔河壩主體為混凝土支墩大頭壩，它是由中部的12個(gè)雙支墩（ 212號(hào)） 壩段、兩壩頭各3個(gè)單支墩壩段、重力壩段、左岸2號(hào)垛以左（壩軸線向上

3、游轉(zhuǎn)折42°）的轉(zhuǎn)折重力壩段組成。其壩頂全長(zhǎng)331.0 m，最大壩高83.1 m，壩頂高程 202.0 m，防浪墻頂高程204 m，正常蓄水位187.0 m，設(shè)計(jì)洪水位201.19 m，對(duì)應(yīng)壩前水深79.19 m，對(duì)應(yīng)庫(kù)容為3.47×108 m3，校核洪水位203.93 m，對(duì)應(yīng)壩前水深81.93 m，對(duì)應(yīng)庫(kù)容3.47×108 m3。2 潰口洪水計(jì)算2.1 確定潰口形態(tài)潰口是大壩潰決失事時(shí)形成的口門(mén)，其形態(tài)主要取決于壩型和筑壩材料。目前，對(duì)實(shí)際潰壩機(jī)理的研究尚不明確，因此，潰口的形態(tài)多用近似假定的方式確定。一般情況下，潰口形態(tài)是由4個(gè)參數(shù)確定的，即潰決歷時(shí)、潰口寬

4、度b、潰口底部高程hb、潰口邊坡的坡比m。第一個(gè)參數(shù)可以確定大壩的潰決形式是瞬時(shí)潰壩還是逐漸潰壩，后面的3個(gè)參數(shù)則可以確定潰口的斷面形狀矩形、三角形、梯形，進(jìn)而判斷大壩是局部潰決還是全部潰決。2.2 計(jì)算潰壩最大流量潰壩壩址的峰值流量與壩址潰壩前上下游水深、壩址斷面形狀和尺寸有關(guān)，而其計(jì)算方法則取決于潰壩形式瞬時(shí)全潰、瞬時(shí)部分潰、逐漸潰決。由于筑壩材料和結(jié)構(gòu)形式不同，所以，對(duì)于混凝土重力壩或拱壩等剛性體壩，經(jīng)常出現(xiàn)一次性瞬時(shí)全潰或瞬時(shí)部分潰的情況。此次計(jì)算結(jié)合工程的實(shí)際情況，遵循“最不利，最可能”的原則計(jì)算壩體瞬時(shí)全潰、瞬時(shí)部分潰2種情況下的相關(guān)值。對(duì)于瞬時(shí)潰壩，水流呈非恒定、不連續(xù)的運(yùn)動(dòng)狀態(tài)

5、。目前，國(guó)內(nèi)外有不少潰壩潰口流量的計(jì)算模型，比如謝任之提出的統(tǒng)一公式，sigh和scarlato等通過(guò)實(shí)驗(yàn)觀察修正的寬頂堰流量計(jì)算公式，還有美國(guó)國(guó)家氣象局fldwav模型中使用的潰口流量公式等。本文擬采用謝任之提出的統(tǒng)一公式計(jì)算瞬時(shí)潰壩的最大流量，即：. （1）式（1）中：為流量系數(shù)；b0為潰口平均寬度，m；g為重力加速度，m/s2；h0為潰壩時(shí)壩前平均水頭，m。當(dāng)大壩瞬時(shí)全潰時(shí)：. （2）當(dāng)大壩瞬時(shí)部分潰時(shí)：. （3）式（2）（3）中：m為斷面形狀系數(shù)；為沉溺系數(shù)；e為矩形斷面、自由出流、口門(mén)拉至河底、堰寬比為e的流量系數(shù)；f為堰高比；n2，n4，n6分別為指數(shù)。2.3 壩址斷面潰壩流量過(guò)程

6、線潰壩流量過(guò)程線可以概化為4次拋物線型或2.5次拋物線型，即當(dāng)潰壩初瞬流量陡增至qmax時(shí)，在短時(shí)間內(nèi)流量迅速降為潰壩前入庫(kù)流量q0所形成的下凹曲線。在工程中，多選擇4次拋物線型。其具體工作流程是先擬定過(guò)程線總歷時(shí)t，再驗(yàn)算過(guò)程線與q=q0直線間的水量是否等于潰壩庫(kù)容。如果不相等，則需要調(diào)整t值，直到二者相等為止。3 計(jì)算方案的選取和結(jié)果分析3.1 計(jì)算方案方案潰壩條件：瞬時(shí)全潰，庫(kù)水位取校核洪水位，即203.93 m；潰口底部高程為122 m；下泄庫(kù)容為3.47×108 m3；潰口寬度取壩址斷面平均河谷寬度，約為165.5 m；水庫(kù)上游來(lái)流量取與校核洪水位對(duì)應(yīng)的洪峰流量，即12 4

7、00 m?/s。方案潰壩條件：瞬時(shí)部分潰，庫(kù)水位取校核洪水位，即203.93 m；潰口底部高程為163 m；下泄庫(kù)容為2.94×108 m3；潰口寬度取壩址斷面潰口底部以上平均河谷的寬度，約為220 m；水庫(kù)上游來(lái)流量取與校核洪水位對(duì)應(yīng)的洪峰流量，即12 400 m?/s。方案潰壩條件：瞬時(shí)全潰，庫(kù)水位取設(shè)計(jì)洪水位，即201.19 m；潰口底部高程為122 m；下泄庫(kù)容為3.47×108 m3；潰口寬度取壩址斷面平均河谷寬度，約為165.5 m；水庫(kù)上游來(lái)流量取與設(shè)計(jì)洪水位對(duì)應(yīng)的洪峰流量，即6 640 m3/s。方案潰壩條件：瞬時(shí)部分潰，庫(kù)水位取設(shè)計(jì)洪水位，即201.19

8、m；潰口底部高程為163 m；下泄庫(kù)容為2.665×108 m3；潰口寬度取壩址斷面潰口底部以上平均河谷寬度，約為220 m；水庫(kù)上游來(lái)流量取與設(shè)計(jì)洪水位對(duì)應(yīng)的洪峰流量，即6 640 m3/s。3.2 結(jié)果分析潰壩流量計(jì)算結(jié)果詳見(jiàn)表1，4種方案潰壩流量過(guò)程線如圖1、圖2、圖3、圖4所示。由表1可知，通過(guò)計(jì)算4種方案的相關(guān)值，可得到壩址處潰壩峰值流量和下泄庫(kù)容放空用時(shí)。方案潰壩峰值流量最大，方案潰壩峰值流量最小，方案下泄庫(kù)容放空用時(shí)最長(zhǎng)，方案下泄庫(kù)容放空用時(shí)最短。由圖1、圖2、圖3、圖4可知，在4種潰壩條件下，在分析壩址斷面處潰口流量過(guò)程線的變化時(shí)，并沒(méi)有考慮到岸邊溢洪道、新老泄洪洞對(duì)

9、潰口流量過(guò)程線的影響，而潰口流量在潰壩初迅速達(dá)到了潰壩峰值流量qmax，之后又隨時(shí)間迅速展平，最后達(dá)到潰壩前入庫(kù)流量q0.比較圖1和圖3可知，壩前水頭越高，潰壩洪水波攜帶的重力勢(shì)能就越高，潰壩峰值流量就越大，破壞能力就越強(qiáng)。4 結(jié)論磨子潭水庫(kù)多年平均徑流量為17.6 m?/s，歷年最大洪峰流量為9 350 m?/s，大壩采用千年一遇的校核標(biāo)準(zhǔn)，與其對(duì)應(yīng)的流量為12 400 m?/s。當(dāng)水庫(kù)遇到入庫(kù)洪水量超過(guò)設(shè)計(jì)洪水量或是其他突發(fā)情況，比如地震、戰(zhàn)爭(zhēng)等時(shí)，存在一定的潰壩風(fēng)險(xiǎn)。即使方案的潰壩峰值流量為66 039 m?/s，其真正流量也遠(yuǎn)大于設(shè)計(jì)流量。文中潰口流量過(guò)程線的計(jì)算中并沒(méi)有考慮到潰口形成

10、過(guò)程對(duì)其的影響，但是，在“最不利，最可能”的原則下，這種計(jì)算方式也是比較合理的，它能滿足工程的實(shí)際需要。在此過(guò)程中，較為準(zhǔn)確的潰口流量過(guò)程線的計(jì)算方法對(duì)計(jì)算潰壩洪水波演進(jìn)、分析潰壩風(fēng)險(xiǎn)、制訂地區(qū)應(yīng)急預(yù)案是非常有利的。因此，后續(xù)探尋較為準(zhǔn)確的潰口流量過(guò)程線計(jì)算方法是很有必要的。參考文獻(xiàn)1楊志，馮民權(quán).潰口近區(qū)二維數(shù)值模擬與潰壩洪水演進(jìn)模擬j.水利水運(yùn)工程學(xué)報(bào)，2015（1）.2劉慶紅，付湘，何海源，等.梯級(jí)水庫(kù)群洪水潰壩模擬j.中國(guó)農(nóng)村水利水電，2011（3）.3陳春燕，李遠(yuǎn)強(qiáng)，吳春耀，等.梯級(jí)水庫(kù)連續(xù)潰壩計(jì)算j.廣東水利水電，2011（9）.4張建國(guó)，王乘，徐學(xué)軍.唐家山堰塞湖潰口流量數(shù)值模擬j.人民長(zhǎng)江，2009，40（22）.5謝任之.潰壩水力學(xué)m.濟(jì)南：山東科學(xué)技術(shù)出版社，1993.6李煒