陳祖煜，胡貴良，張 強(qiáng)，3，吳帥峰

(1.中國水利水電科學(xué)研究院，北京 100048；2.華電金沙江上游水電開發(fā)有限公司，四川 成都 610015；3.西安理工大學(xué)，陜西 西安 710048)

0 引 言

2018年10月和11月，我國金沙江上游白格村同一位置連續(xù)發(fā)生兩次大規(guī)模的山體滑坡，堵塞金沙江干流形成堰塞湖[1-2]。金沙江“11.03”白格堰塞湖潛在庫容達(dá)7.9億m3，現(xiàn)場應(yīng)急救災(zāi)指揮部啟動了人工干預(yù)開挖泄流槽、蘇洼龍水電站圍堰破拆、中游梯級水電站應(yīng)急騰庫等處置措施。11月13日堰塞湖潰決，形成超萬年一遇洪水，在各部門統(tǒng)一協(xié)調(diào)應(yīng)急搶險救災(zāi)過程中，及時轉(zhuǎn)移受威脅地區(qū)群眾，實現(xiàn)無一人傷亡，但由于此次潰決洪水超萬年一遇，對下游在建水電站和金沙江沿岸多個重要鄉(xiāng)鎮(zhèn)仍產(chǎn)生了巨大影響，直接經(jīng)濟(jì)損失約達(dá)135億元。

梯級開發(fā)是河流開發(fā)的主要形式[3]。梯級水庫群不僅具有增強(qiáng)河川徑流調(diào)節(jié)能力、提高河流水資源利用率，而且可在防洪抗旱、水力發(fā)電和生態(tài)保護(hù)等方面發(fā)揮更大的作用[4-6]。目前，我國已形成金沙江、雅礱江、大渡河、烏江、怒江、長江上游、南盤江、紅水河、瀾滄江等大江大河流域梯級水庫群。2014年云南省魯?shù)榭h“8.03”地震形成牛欄江紅石巖堰塞湖，堰塞湖發(fā)生后，上下游梯級水庫采取了上游德慶水庫蓄水?dāng)r截、下游天花板和黃角樹電站預(yù)泄騰庫等應(yīng)急處置措施，為堰塞湖搶險贏得了寶貴時間[7]。

金沙江中游梯級水電站是金沙江梯級水電基地的重要組成部分，“11.03”白格堰塞湖發(fā)生后，梨園、阿海、金安橋等水電站統(tǒng)一調(diào)度防洪庫容13.53億m3，阻斷了破壞性潰決洪水向大理、楚雄方向的延伸。相關(guān)監(jiān)測結(jié)果顯示，白格堰塞體高邊坡仍處于緩慢變形狀態(tài)，存在再次發(fā)生滑坡堵江的可能[8]。本文圍繞“11.03”白格堰塞湖實際潰決過程，運用一維洪水非恒定流演進(jìn)模型，對金沙江中游河段洪水演進(jìn)過程進(jìn)行反演分析，與梯級水庫應(yīng)急處置后實測洪水過程進(jìn)行了對比分析，研究了梯級水庫在堰塞湖處置過程中的作用，為以后此類堰塞湖應(yīng)急搶險中梯級水庫應(yīng)急處置提供借鑒。

1 金沙江梯級水電站概況

金沙江上游是指青海玉樹(巴塘河口)至云南石鼓區(qū)間段，河段長974 km，流域面積7.65萬km2，落差約1 715 m，河道平均坡降1.76%。金沙江上游是金沙江水電基地的重要組成部分，規(guī)劃實施“一庫八級”梯級水電站建設(shè)，分別為崗?fù)?、波羅、葉巴灘、拉哇、巴塘、蘇洼龍、昌波、旭龍、奔子欄。白格堰塞湖位于葉巴灘壩址上游54 km處，葉巴灘、拉哇、巴塘、蘇洼龍水電站已核準(zhǔn)在建。

金沙江中游西起云南麗江石鼓鎮(zhèn)，東至攀枝花市的雅礱江口，長564 km，落差838 m，河道平均坡降1.48%。河段規(guī)劃為“一庫八級”開發(fā)，目前投入運行并能夠發(fā)揮調(diào)洪作用的水電站有6座，分別為梨園水電站、阿海水電站、金安橋水電站、龍開口水電站、魯?shù)乩娬竞陀^音巖水電站。金沙江梯級水電站如圖1所示。

圖1 金沙江梯級水電站

兩次白格堰塞湖發(fā)生后，由于上游梯級多數(shù)在建，不能發(fā)揮調(diào)洪削峰作用，在長江委的統(tǒng)一調(diào)度下，中游已建梯級水電站分別采取降低庫水位，騰出庫容以容納堰塞湖潰決超大洪水。本文主要研究中游梯級調(diào)洪削峰作用，中游梯級水庫特性如表1所示。

表1 中游梯級水庫特性

2 金沙江中游水電站應(yīng)急處置過程

2.1 “11.03”洪水過程

2018年11月3日，西藏自治區(qū)昌都市江達(dá)縣波羅鄉(xiāng)和四川省甘孜藏族自治州白玉縣絨蓋鄉(xiāng)交界處，“10.10”白格滑坡體殘留體再次下滑，堵塞由第一次堰塞體潰決形成的泄流槽后，再次堵塞金沙江干流，形成“11.03”白格堰塞湖。堰頂埡口寬約195 m，長約273 m，堰頂高程約2 966 m，較“10.10”堰塞體頂高程高30余m，堰塞體高出水面58.24 m，潛在庫容將達(dá)7.9億m3，一旦潰決，將會對下游在建水電站和沿岸人民群眾生命財產(chǎn)安全產(chǎn)生巨大影響[9-10]。

2018年11月5日經(jīng)現(xiàn)場查勘和分析，現(xiàn)場應(yīng)急救災(zāi)指揮部啟動了人工干預(yù)開挖泄流槽，現(xiàn)場先后調(diào)用工程機(jī)械18臺，其中，挖掘機(jī)13臺，裝載機(jī)5臺。至11月11日下午泄流槽開挖完成，泄流槽長220 m，泄流槽底坎高程2 952.52 m，頂寬42 m，底寬3 m，最大深度15 m[1，10-11]。與此同時，為防止堰塞湖潰決洪水造成蘇洼龍圍堰潰決，產(chǎn)生潰決洪水疊加效應(yīng)，采取破拆蘇洼龍上游圍堰措施，極大限度地保護(hù)下游沿岸人民群眾的安全。

11月12日凌晨04∶45，泄流槽開始進(jìn)水，泄流槽河道在10∶50被淹沒；13日8∶00，泄流槽過流流量明顯增大，13日12∶00沖刷明顯加快；13日13∶40過流流量超過堰塞湖入湖流量，堰塞湖壩前最高水位2 956.40 m，相應(yīng)庫容5.78億m3；13日18∶00潰決流量達(dá)到最大流量31 000 m3/s。15日8∶00，泄流槽泄流量下降至天然來流量，水位穩(wěn)定，險情解除。下游各電站洪水過程如圖2所示。

圖2 “11.03”白格堰塞湖潰決下游各流量測點洪水過程

2.2 水電站應(yīng)急調(diào)度過程

金沙江“11.03”白格堰塞湖發(fā)生后，長江委根據(jù)金沙江白格堰塞湖人工干預(yù)情況，對金沙江中游梯級水電站實施聯(lián)合洪水調(diào)度，控制水位以應(yīng)對洪水。截至2018年11月13日16∶00，調(diào)度令梨園水電站以不高于1 592 m水位運行，最終實現(xiàn)騰庫最低水位1 590.08 m，梨園水庫防洪高水位1 618.00 m，正常蓄水位1 618.00 m，死水位1 605.00 m，實現(xiàn)3.8億m3的防洪庫容；阿海水電站以不高于1 493.3 m水位運行；金安橋水電站以不高于1 406 m水位運行；龍開口水電站、魯?shù)乩娬竞陀^音巖水電站分別以不高于1 294.00、1 220.00、1 130.00 m水位運行，實際調(diào)度防洪庫容13.53億m3。金沙江中游梯級水庫應(yīng)急調(diào)度情況如表2所示[12]。

表2 金沙江中游梯級水庫應(yīng)急調(diào)度情況[12]

2018年11月15日12∶30，白格堰塞湖潰決洪水以7 200 m3/s的洪峰流量入庫梨園水電站。梨園水電站最大泄流量為4 500 m3/s，經(jīng)過梨園、阿海、金安橋、龍開口、魯?shù)乩?、觀音巖等大型水電站進(jìn)行攔蓄調(diào)節(jié)，阿海電站水庫的出庫流量僅為854 m3/s，遠(yuǎn)低于11月3日堰塞湖發(fā)生斷流前的下游天然河道中的流量2 040 m3/s。

3 金沙江中游梯級水電站減災(zāi)效果評價

3.1 計算原理

洪水演進(jìn)計算的本質(zhì)是求解淺水運動方程，國內(nèi)外許多專家學(xué)者進(jìn)行了一維、二維、三維的求解方法研究，由于梨園至攀枝花市金沙江河道長350 km，本文洪水演進(jìn)采用MIKE11模型計算，MIKE11模型由丹麥水力研究所(DHI)開發(fā)[13]，在國內(nèi)外得到了廣泛應(yīng)用[14-16]。

模型的基本方程組為圣維南方程

(1)

式中，x、t分別為計算點空間和時間的坐標(biāo)；A為過水?dāng)嗝婷娣e；Q為過流流量；h為水位；q為旁側(cè)入流流量；C為謝才系數(shù)；R為水力半徑；α為動量校正系數(shù)；g為重力加速度。

MIKE11模型采用Abbott六點中心隱式差分方法圣維南方程進(jìn)行數(shù)值離散后應(yīng)用“追趕法”求解差分方程[13]。計算過程中在每個網(wǎng)格點不同時計算水位、流量值，而是按照順序交替進(jìn)行計算。

3.2 模型率定

“11.03”白格堰塞湖洪水發(fā)生后，在長江委的調(diào)度下，梨園、阿海、金安橋、龍開口、魯?shù)乩?、觀音巖等大型水電站進(jìn)行攔蓄調(diào)節(jié)，梨園下游無實測洪水資料，無法進(jìn)行洪水演進(jìn)模型的率定。巴塘至石鼓金沙江河段長974 km，落差約1 715 m，河道平均坡降1.76%，石鼓至攀枝花河段長564 km，落差838 m，河道平均坡降1.48%，采用蘇洼龍至梨園的實測流量過程資料進(jìn)行模型的率定，地形資料采用美國對地觀測衛(wèi)星Terra提供的30 m分辨率DEM地形資料。

蘇洼龍水電站至梨園水電站全長447 km，入流流量采用蘇洼龍?zhí)帉崪y“11.03”白格潰決洪水流量過程，如圖3所示。曼寧系數(shù)選取為0.06。圖4所示為蘇洼龍水電站至梨園水電站洪水演進(jìn)模擬結(jié)果。

圖3 蘇洼龍計算入流流量過程

模擬結(jié)果表明，曼寧參數(shù)取值為0.06時，模擬結(jié)果與實測結(jié)果擬合效果較好。奔子欄水電站實測洪峰流量15 700 m3/s，模擬洪峰流量14 090 m3/s，誤差10.25%；梨園水電站實測洪峰流量7 410 m3/s，模擬流量7 870 m3/s，誤差為6.21%。模擬洪峰到達(dá)時間與實測過程基本一致，滿足計算要求。

圖4 蘇洼龍水電站至梨園水電站模擬結(jié)果

3.3 洪水演進(jìn)模擬

“11.03”白格堰塞湖由于梨園水電站提前泄水調(diào)洪的作用，潰決洪水未對下游造成較大的影響，現(xiàn)假定梨園水電站未采取任何措施，計算白格堰塞湖潰決洪水自梨園水電站演進(jìn)至攀枝花的洪水過程，梨園水電站至攀枝花市全長350 km，計算地形資料采用30 m DEM資料。曼寧系數(shù)選用蘇洼龍水電站至梨園水電站率定的0.06。入庫流量采用“11.03”實測洪水梨園水電站入庫流量，入流流量過程如圖5所示。

圖5 計算入流流量過程

以梨園水電站至攀枝花市沿江9個典型城鎮(zhèn)為計算斷面，分別選取這些斷面計算洪水流量過程和水位變化過程，洪水演進(jìn)流量過程如圖6所示。

圖6 梨園-攀枝花市模擬結(jié)果

將斷面位置的流量與實際地形相結(jié)合，依據(jù)河道特征斷面水位—流量關(guān)系曲線計算水位抬升高度，水位抬升高度及各位置的淹沒面積如表3所示。

表3 淹沒信息統(tǒng)計

由圖6和表3可見，若無中游梯級水庫調(diào)蓄洪水作用，麗江市、楚雄市和攀枝花市將增加淹沒面積總計約160 km2，必將對沿江造成更大的損失。通過梯級水電站洪水調(diào)度及科學(xué)應(yīng)對，雖然梨園上游受洪水災(zāi)情較重，但梨園、阿海等大型電站水庫發(fā)揮了攔洪削峰調(diào)蓄的良好作用，將洪水災(zāi)害截斷于中游梯級電站水庫當(dāng)中。

4 結(jié) 論

本文對金沙江“11.03”白格堰塞湖的洪水過程進(jìn)行了回顧，采用一維洪水演進(jìn)分析水力學(xué)模型，對金沙江中游梯級水電站在防洪減災(zāi)中的作用進(jìn)行了分析，得出以下結(jié)論：

(1) “11.03”白格堰塞湖潰決在金沙江形成超萬年一遇洪水，通過開挖泄流槽、在建電站圍堰拆除和中游梯級水庫的調(diào)節(jié)作用，成功處置了此次超大洪水，中游梯級水庫在正確的調(diào)度作用下，減小了下游河道沿岸防洪壓力。

(2) 梨園、阿海、金安橋、龍開口、魯?shù)乩?、觀音巖等水電站經(jīng)聯(lián)合洪水調(diào)度，共調(diào)度防洪庫容13.53億m3，截斷了破壞性的潰壩洪水向大理、楚雄方向延伸，減輕了云南省的災(zāi)情。若無梯級水庫的作用，將造成額外的淹沒面積約160 km2。

(3) 金沙江上游梯級水庫正在建設(shè)期，在此次應(yīng)急處置過程中未發(fā)揮調(diào)洪作用，雖然中游梯級水電站成功攔截洪水，但梨園上游仍然受此次洪水災(zāi)情較重。應(yīng)加快金沙江上游梯級電站的建設(shè)，形成梯級規(guī)模，如全河段可進(jìn)行梯級調(diào)度，必將更多地減少金沙江沿岸損失。