郭紅民，張 楊，柳 滔，楊百銀，馬 良，楊 瑩

(1.三峽大學(xué)水利與環(huán)境學(xué)院，湖北宜昌443002；2.廣西右江水利開發(fā)有限責(zé)任公司，廣西南寧530000；3.水電水利規(guī)劃設(shè)計(jì)總院，北京100120)

1 研究背景

混凝土壩由于材料強(qiáng)度較高且一致性較好，一般不易潰決，但在某些特殊條件下，如戰(zhàn)爭(zhēng)、地震、風(fēng)暴潮、質(zhì)量缺陷等，也可能發(fā)生潰決。大壩一旦潰決，巨大的水位落差，使水庫中的蓄水經(jīng)過潰口快速下泄到下游河道中，往往造成災(zāi)害性的突發(fā)洪水。

大壩不同的潰決方式及潰口形式對(duì)潰口的最大流量及相關(guān)水力特性影響很大，因此學(xué)者都較為重視，提出了許多潰壩潰決研究成果。如姚霄雯等[1]收集了國內(nèi)外74座混凝土壩潰壩資料，總結(jié)了混凝土壩的潰壩特點(diǎn)，統(tǒng)計(jì)了主要潰壩原因，進(jìn)而提出了壩基、壩肩破壞為混凝土壩最重要的潰壩模式；謝任之[2]系統(tǒng)的提出了可用于計(jì)算瞬時(shí)潰、逐漸潰、部分潰、全潰等情況的潰壩流量計(jì)算統(tǒng)一公式，并且還提出了適用于不同河道斷面的潰壩流量過程線計(jì)算方法；曲霞[3]利用MIKE11模型對(duì)大凌河的潰堤洪水過程進(jìn)行模擬分析，并選取實(shí)測(cè)洪水資料，對(duì)比模擬結(jié)果，表明能夠較好地反映潰口洪水實(shí)際狀況；姜治兵等[4]建立了壩體潰決過程與潰壩洪水演進(jìn)耦合的平面二維數(shù)值模型，通過實(shí)際算例表明模型合理地模擬了潰口的發(fā)展與洪水演進(jìn)過程；郭紅民等[5]根據(jù)混凝土壩結(jié)構(gòu)特征和受力特點(diǎn)，假定了3種可能的潰決方式，通過物理模型試驗(yàn)研究了3種潰決方式的潰口水力學(xué)特征，結(jié)果表明，3種潰決方式中滑動(dòng)式潰口峰值流量最大，危害最大。

以往大部分潰壩研究主要討論水庫大壩的潰決流量過程或下游洪水的演進(jìn)特征，而對(duì)于同一水庫大壩考慮不同潰決形式下的潰口水力特性研究較少。本文基于混凝土壩瞬時(shí)潰決的基本假設(shè)，采用數(shù)值模擬計(jì)算的方法，以某一實(shí)際水庫邊界及混凝土壩為研究背景，比較研究了大壩橫向局部潰決和垂向局部潰決下的潰口水力特性，為類似工程的大壩潰決及洪水傳播研究提供參考。

2 模型構(gòu)建

2.1 研究區(qū)域概況

某水電站為混凝土壩，壩高305 m，壩頂高程1 885 m，水庫正常蓄水位1 880 m，總庫容77.6億m3。壩址處多年平均流量1 220 m3/s，多年平均年徑流量385億m3/s，本文的研究區(qū)域包括水庫回水的干流河段及大壩下游的部分河段，模型中河道水深數(shù)據(jù)是根據(jù)實(shí)測(cè)斷面資料經(jīng)過插值加密處理后得到，斷面間距為200 m，陸上地形高程來自ASTER GDEM V2[6]，并由當(dāng)?shù)睾拥浪娓叱踢M(jìn)行了矯正處理，計(jì)算中在壩址上、下游各200 m處設(shè)置了水位觀測(cè)點(diǎn)，用于觀測(cè)庫水位及下游河道水位的變化過程。

2.2 控制方程

本文構(gòu)建的潰壩數(shù)值模型采用控制方程為

(1)

(2)

式中，Q為流量，m3/s；q為側(cè)向入流，m3/s；A為過水面積，m2；h為水位，m；R為水力半徑，m；C為謝才系數(shù)；α為動(dòng)量修正系數(shù)。

2.3 計(jì)算初始條件

上游入庫流量采用萬年一遇洪水作為入流條件，由于該超標(biāo)洪水過程的洪水歷時(shí)較長(zhǎng)，而在漫頂潰決條件下長(zhǎng)時(shí)段的過程線對(duì)潰口峰值流量影響不大，故本文取最大3天洪峰的平均值作為來流洪水進(jìn)行計(jì)算。

潰壩計(jì)算基于大壩漫頂時(shí)瞬時(shí)潰決，即大壩潰決時(shí)壩前水位為壩頂高程1 885 m；壩下起始水位為樞紐下泄上游來流量時(shí)的下游河道天然水位。

2.4 模型驗(yàn)證

由于所研究的大壩實(shí)際未發(fā)生過潰決，故模型的驗(yàn)證采用水庫正常運(yùn)行時(shí)段的水位特征進(jìn)行。驗(yàn)證時(shí)段為某年7月～8月，驗(yàn)證結(jié)果如圖1所示，計(jì)算值與實(shí)測(cè)值基本吻合，滿足進(jìn)一步研究的需要。

圖1 水文站水位驗(yàn)證

3 計(jì)算工況

為研究混凝土壩瞬潰條件下不同潰決形式對(duì)潰口水力特性的影響，本文根據(jù)實(shí)際工程可能產(chǎn)生的潰決形式，將混凝土壩的潰決形式概括為橫向局部潰決[7-9](在外力的作用下壩體發(fā)生結(jié)構(gòu)破壞，導(dǎo)致部分壩段整體失穩(wěn)或脫落)和垂向局部潰決[10-11](大壩潰決在某個(gè)高度發(fā)生，其下留有一定殘高)。

為方便后期計(jì)算成果的比較，在設(shè)定橫向局部潰決和垂向局部潰決的計(jì)算工況時(shí)，采用了潰口面積相等的潰口特征，即將大壩擋水橫剖面面積A分別以橫向、垂向局部潰決的潰口面積劃分為1/4A、2/4A和3/4A兩類共6種計(jì)算工況。圖2、圖3分別為橫向局部潰決(工況1、2、3)、垂向局部潰決(工況4、5、6)的示意。

4 計(jì)算結(jié)果及分析

4.1 潰口流量

各潰決工況下潰口流量[12-14]過程線如圖4所示，6種工況的潰口流量在潰壩發(fā)生后均在極短時(shí)間(約0.05～0.1 h)內(nèi)達(dá)到最大流量。表1為等面積潰決潰口最大流量。

當(dāng)橫向潰決時(shí)，潰口最大流量隨著潰口面積的增加而增加，但增速并不呈線性，當(dāng)潰口面積較小時(shí)，增速較大，當(dāng)潰口面積較大時(shí)，增速變緩；在庫容相同情況下，潰口流量降速的快慢與潰口的大小、庫容泄空的時(shí)間直接相關(guān)，如工況1潰口面積最小，庫容泄空時(shí)間長(zhǎng)，其潰口流量的下降速度就較緩，而工況2、3的潰口面積均較大，庫容泄空時(shí)間較短，其潰口流量的下降速度均較快，下降過程差別不大。

圖2 橫向局部潰決的3種工況

圖3 垂向局部潰決的3種工況

圖4 各工況下潰口流量過程

表1 等面積潰決潰口最大流量

垂向潰決由于存在殘高底坎，使各工況下泄的庫容不同，潰口流量總體仍隨著潰口面積的增大而增大，且最大潰口流量的大小與殘坎上的起始潰決水頭高度有關(guān)。工況6時(shí)，由于潰口面積最大，起始水頭也最大，潰口最大流量也最大。

由表1可知，當(dāng)潰決面積相同時(shí)，橫向潰決的最大潰口流量分別是垂向潰決的2.72、1.88和1.40倍，說明橫向潰決的下降趨勢(shì)比垂向潰決要快，水庫泄空時(shí)間要短，這種差別主要是由于橫向潰決的潰口部分是一潰到底，同面積條件下，潰口水頭較高，潰口流量較大；而垂向潰決時(shí)，潰口留有一定殘高，同面積潰決情況下，潰口水頭相對(duì)較小，潰口流量也相對(duì)較小，流量下降也明顯變緩。

4.2 流速

各工況潰口斷面流速變化過程如圖5所示，各潰決工況下潰口平均流速都是在潰決開始階段快速增大，在極短的時(shí)間(約0.06 h)內(nèi)達(dá)到最大值，然后再逐漸減小。表2為等面積潰決潰口最大平均流速。

圖5 各工況下潰口處平均流速

表2 等面積潰決潰口最大平均流速

當(dāng)橫向潰決時(shí)，潰口最大流速隨著潰口面積的增加而減小，但3種工況下潰口最大流速的差值并不大；在庫容相同情況下，潰口流速下降的快慢與潰口面積、庫容泄空時(shí)間有關(guān)，潰口面積較小時(shí)，泄空時(shí)間長(zhǎng)，流速下降慢，潰口面積較大時(shí)，泄空時(shí)間短，流速下降快；如工況1潰口面積最小，庫容泄空時(shí)間最長(zhǎng)，其潰口流速的下降速度就較緩，而工況3的潰口面積最大，庫容泄空時(shí)間最短，其潰口流速的下降速度也最快。

當(dāng)垂向潰決時(shí)，潰口最大流速總體仍隨著潰口面積的增大而增大，但3種工況下潰口最大流速的差別較大；垂向潰決由于存在殘高底坎，使各工況的下泄庫容不同，當(dāng)殘高底坎較高，下泄庫容較小時(shí)，潰口最大流速也較小，下降趨勢(shì)也更緩，殘高底坎較低，下泄庫容較大時(shí)，潰口最大流速也較大，下降趨勢(shì)也更快，如工況6時(shí)由于潰口面積最大，殘高底坎最低，下泄庫容也最大，潰口最大流速也最大，下降速率也更快。

由表2可知，等面積潰決時(shí)，橫向潰決的潰口最大流速分別是垂向潰決的2.69、1.78和1.30倍，橫向潰決比垂向潰決具有更大的流速，且流速變化更劇烈，這是因?yàn)橥娣e條件下，橫向潰決庫容泄空時(shí)間短，潰口水頭較高，潰口流速大。

4.3 水位

混凝土壩在各潰決工況時(shí)，上、下游水位測(cè)站處的水位變化情況如圖6、7所示，水位變化過程直接與水庫庫容及潰口流量變化過程相關(guān)。

圖6 上游水位變化

圖7 下游水位變化

壩體一旦潰決，上游庫區(qū)水位即快速下降，下降速率以工況3最大，工況4最小。當(dāng)橫向潰決時(shí)，隨著潰口面積的增加，水庫放空時(shí)間較快，水位下降則越快，趨于穩(wěn)定的耗時(shí)越短，由于潰決后的底部高程相同，潰決后3種工況的穩(wěn)定水位相同；當(dāng)垂向潰決時(shí)，3種計(jì)算工況(4、5、6)中，隨著潰口面積的增加，水庫潰決水頭也相應(yīng)增加，水位的降速也相應(yīng)增加，但最后水位趨于穩(wěn)定的耗時(shí)大致相同。潰決面積相同時(shí)，橫向潰決比垂向潰決水位的下降速率一般更快。

潰壩下游水位變化，工況1～6的水位上漲時(shí)間大致都在約0.1 h時(shí)達(dá)到最大，下游水位上升最高的峰值是工況3，最小的是工況4，水位達(dá)到峰值后均逐漸下降，之后趨于平穩(wěn)。橫向潰決時(shí)，隨著潰口面積的增加，潰口最大流量增加，水庫放空時(shí)間較快，水位變化呈“尖”“陡”的特點(diǎn)，趨于平衡的耗時(shí)較短；垂向潰決時(shí)，由于水庫放空時(shí)間相對(duì)較慢，水位的變化總體呈現(xiàn)快速上升緩慢下降的形態(tài)，趨于平衡的耗時(shí)越長(zhǎng)。

5 結(jié) 論

(1)混凝土壩橫向潰決和垂向潰決的潰口流量都是隨著潰口面積的增加而增大，而橫向潰決的總體潰口流量較大，下降較快，泄空時(shí)間較短，垂向潰決的總體下降趨勢(shì)比橫向潰決要緩；等面積潰決時(shí)，隨著潰口面積的增大，橫向潰決的最大流量分別是垂向潰決的2.72、1.88和1.40倍，即等面積潰決時(shí)橫向潰決比垂向潰決具有更大的危害，潰口面積越小，差別越明顯。

(2)橫向潰決由于潰口水頭大，潰口處平均流速整體較大，泄空時(shí)間隨著潰口面積的增大而縮短，流速下降也越快；垂向潰決的潰口水頭相對(duì)較小，流速相對(duì)較小，下降相對(duì)較慢；等面積潰決時(shí)，橫向潰決比垂向潰決具有更大的流速，且流速下降變化更劇烈。

(3)橫向潰決隨著潰口面積的增加，水庫放空時(shí)間加快，上、下游的水位變化也越快，趨于穩(wěn)定的耗時(shí)越短；垂向潰決隨著潰口面積的增加，水庫水位的下降速率也增加，但水庫的放空時(shí)間差別不大，下游水位在快速上升后，其下降過程均相對(duì)較緩慢，而趨于穩(wěn)定的耗時(shí)則大致相同。