抽水蓄能電站潰壩過程模擬及應急響應措施研究
——以梅州抽水蓄能電站為例

周水兵，金 苗，江廷鋒，陳光澤，李東豪

(1.南方電網調峰調頻發電有限公司 運行分公司，廣州 514400；2.梅州蓄能發電有限公司，廣東 梅州 514600；3.南方電網調峰調頻發電有限公司 工程建設管理分公司，廣州 514400；4.珠江水利委員會 珠江水利科學研究院，廣州 510611)

1 概述

目前，大力發展清潔可再生能源成為全球應對環境污染和氣候變化的共同選擇。抽水蓄能作為經濟最優、最具大規模開發條件的一種綠色低碳可再生能源利用方式，在全球范圍內得到了廣泛應用[1-2]。根據《抽水蓄能中長期發展規劃(2021—2035 年)》，從“十四五”開始，未來我國抽水蓄能電站將持續增多。隨著近年來氣候變化的加劇，極端水文事件發生的概率顯著增加，抽水蓄能電站水庫潰壩、產生嚴重安全事故的風險增大。由于潰壩洪水的突發性和極大破壞性，大壩一旦失事潰決，將對下游地區人民群眾生命財產和生態環境造成毀滅性影響[3-5]。

正是由于潰壩洪水的巨大危害性，諸多學者對潰壩洪水進行了研究，其中，潰壩洪水過程模擬和應急響應措施研究是潰壩洪水風險分析和防災減災的重要手段之一。艾子欣[6]等以紅旗水庫為例，采用經驗公式方法計算潰口形式及流量過程，模擬分析不同潰壩工況下的洪水演進過程，編制洪水風險圖，為編制大壩安全管理應急預案提供技術支撐；沈洋等[7]以金牛山水庫為例，采用MIKE軟件耦合的方法，同時考慮不同水庫水位與入庫設計洪水組合情況，計算潰壩洪水及其演進過程；宋利祥[8]基于二維淺水方程和Godunnv型有限體積法，建立了適用于不規則計算域和地形上的潰壩洪水二維數學模型，以漳河水庫為例進行了潰壩洪水分析，繪制了最大淹沒水深、洪水到達時間等洪水風險圖。但上述研究均為單庫潰壩模擬和風險分析，而在具有上、下水庫的抽水蓄能電站中潰壩模擬研究較少[9]。主要是由于受到抽水蓄能電站雙庫水體輪換、調度利用計劃等因素影響，使得電站潰壩時機難控、傳遞和疊加機理復雜，進而較難進行數值模擬。

為此，文章以梅州抽水蓄能電站為研究對象(以下簡稱“梅蓄電站”)，結合電站調度計劃和電站大壩工程特征，考慮多種潛在潰壩風險，研究適合于雙庫抽水蓄能電站的潰壩數值模擬方法，并分析潰壩洪水影響，制定科學合理的應急響應措施，旨在提高電站應急響應能力，保障人民群眾生命財產安全。

2 研究區概況

梅蓄電站位于梅州市五華縣龍村鎮黃獅村內，距廣州市直線距離210 km，系韓江流域琴江支流先水河上游唯一抽水蓄能電站，控制集雨面積為36.38 km2，上下兩庫總庫容為9 180萬m3，裝機容量為2 400 MW(見圖1)。

梅蓄電站主要功能為發電，由上水庫、下水庫、輸水系統和發電廠房組成，按一等工程、一級建筑物設計，上下庫兩庫大壩設計標準和校核標準均分別為500 a和5 000 a。上下兩庫通過引水洞連接，為周調節水庫，無防洪庫容，總庫容控制在4 000萬～5 000萬m3。梅蓄電站上庫正常蓄水位為815.50 m，總庫容為4 254萬m3；主壩為鋼筋混凝土面板堆石壩，壩頂高程為820.00 m，壩頂寬度為10.00 m，最大壩高為60.00 m；副壩為均質土壩，設計壩頂高程為820.00 m，壩頂寬度為7.00 m，最大壩高為11.00 m，設有環形堰豎井式溢洪道，堰頂高程為815.50 m。下庫正常蓄水位為413.5 m，總庫容為4 926萬m3；主壩為碾壓混凝土重力壩，壩頂高程為419.00 m，最大壩高為85.0 m，壩頂長度為317.5 m。副壩右側為長約400 m的單薄分水嶺，防滲形式采用上墻下幕方式，防滲墻深入強風化花崗巖以下1.0 m。設有表孔溢洪道以及泄放管，溢洪道堰頂高程為413.5 m，正常蓄水位時泄放管泄流能力為67 m3/s。根據工程區域DEM數字高程模型分析結果，上庫潰壩后其潰壩水流都將流入下庫中。

3 建模理論及方法

潰壩洪水模擬研究主要包括大壩潰決模擬和下游洪水演進計算。其中，大壩潰決模擬包括潰口發展過程模擬和潰口流量過程計算，下游洪水演進計算包括一維河道洪水演算和二維地表洪水演算。通過大壩潰決模擬可獲得潰口流量過程，以此作為上邊界提供給下游洪水演進模型，最終實現特定水文條件和潰決形式下的潰壩洪水模擬。

3.1 工況設置

抽水蓄能電站潰壩受多類型因素影響，文章綜合考慮各潰壩因素影響，并結合梅蓄電站上下水庫總庫容控制計劃，設定多種潰壩場景。潰壩基本假設為：① 非汛期時，假定電站遭遇地震、戰爭等極端突發事件，壩高大、壩頂長、潰壩影響大的下庫混凝土主壩潰決；② 汛期時，假定水庫遭遇設計標準或超標準洪水，根據上下水庫運行條件不同，考慮上下庫連環潰壩和下庫單獨潰壩。其中，連環潰壩假定上庫初始為正常蓄水位815.5 m，豎井式溢洪道堵塞，無法正常下泄，以致庫內水位上漲，副壩達最高水位全潰，流入下庫后使得下庫副壩全潰，造成連環潰壩；下庫單獨潰壩假定下庫初始為正常蓄水位413.5 m，上水庫不下泄，下庫壩前水位達最高水位時主壩瞬時全潰/副壩全潰。共計設置6種潰壩工況場景(見表1所示)。

3.2 潰口模型

根據潰口形成歷時長短，可將潰壩方式分為瞬時潰決和逐漸潰決。混凝土壩主要考慮瞬時潰決由經驗公式方法計算潰口流量[10]，而土石壩的潰決過程是一個隨潰口寬度和深度增加而逐漸潰決的過程，主要通過基于參數的模型計算方法或基于物理過程的模型計算方法計算，如DAMBR模型、BREACH模型等，相關研究已較為深入[9]。但使用經驗公式方法計算混凝土壩瞬時潰決流量時難以考慮不規則潰口形狀和下游水位頂托的影響，因此文章采用零維水庫耦合二維河道地表網格的方法進行計算。該方法通過特定水庫的水位-庫容關系曲線對水庫進行零維概化，庫區下游河道及周圍地表則使用全二維模型進行概化，通過修正壩址處二維網格節點底高程概化任意潰口形狀，建立零維水庫和局部二維模型的耦合關系，使用基于圣維南方程特征線進行求解的HLL近似黎曼算子計算潰口流量大小，以適應上下游不同的水位、流速。

梅蓄電站上下庫主壩均為混凝土重力壩，副壩均為均質土壩。文章潰壩工況設置中，下庫主壩考慮瞬時全潰，下庫副壩壩長較長，考慮逐漸潰決影響，上庫副壩壩長較短，考慮瞬時全潰。

計算下庫主壩瞬時潰決的潰口流量時，使用零維水庫耦合二維河道地表網格的方法進行計算，零維水庫為下水庫，二維計算區域上游以下庫壩址為界，下游以先水河為界，先水河左、右岸兩側淹沒區以地勢較高的山區坡腳為界，下游采用自由出流邊界，零維水庫與二維模型上邊界進行耦合，模型示意見圖2。計算下庫副壩逐漸潰決的潰口流量時，則采用常用的BREACH模型進行計算。

圖2 零維—二維耦合潰口流量計算模型范圍示意

3.3 潰壩洪水演進模型

文中潰壩洪水演進模型為水庫下游一維河道和洪泛區二維地表側向耦合聯解模型，用于模擬潰壩洪水在河道中及漫溢至河道兩岸農田、村莊后的演進情況。一維河道水流模型采用圣維南方程組作為單一河道非恒定流控制方程，采用有限體積法進行方程離散，運用Riemann求解器計算控制體之間的數值通量以實現模型求解，可解決潰壩水流中干河床、逆坡、復雜流態模擬等問題，適應性強。二維地表洪水演進模型計算采用無結構網格有限體積法，用高性能格式建立適用于不規則計算域和地形上潰壩洪水淹沒區淺水流動的數學模型。數值求解方面，運用基于Godunov格式的有限體積法進行控制方程離散求解，詳細過程可參考文獻[8]。一維、二維模型通過耦合邊界的水力連接條件來實現模型聯解(見圖3)。

圖3 一維—二維模型側向耦合方式示意

文中一維河道模型范圍為梅蓄電站下水庫壩下、琴江(先水河匯口上)至琴江水口水文站以及五華河城區段，建模總河長為136.4 km，共布設109個斷面，斷面間距200～1 000 m，采用近期實測斷面地形資料。一維模型上邊界采用下庫潰口流量過程以及相應的琴江干支流洪水過程，下邊界采用琴江水口站實測水位流量關系曲線。二維地表模型計算區域上游以梅蓄電站下水庫壩址、琴江先水河匯口上游為界，下游以琴江水口水文站為界，干支流兩側淹沒區以地勢較高的上坡腳為界。采用三角網格進行剖分，以準確擬合任意不規則邊界。模型范圍面積約為327.2 km2，共計剖分161 037個單元、85 048個節點。基于1∶10 000DEM，對河道兩側可能淹沒區的網格進行高程插值。同時，結合土地利用類型數據，按照下墊面情況進行計算區域糙率賦值，糙率范圍為0.025～0.1。

將上述一維河道模型與二維地表模型進行側向耦合，耦合邊界共計105個，耦合后的水動力模型可用于模擬梅蓄電站不同潰壩場景下潰壩洪水在河道和河道兩岸地表的演進過程，一二維耦合水動力模型示意見圖4。

圖4 一二維耦合水動力模型示意

3.4 模型計算與結果分析

3.4.1潰口流量模擬結果

非汛期下庫主壩瞬時全潰工況，假定梅蓄電站非汛期遭遇地震、戰爭等極端突發事件，下庫主壩在正常蓄水位瞬時全潰，不考慮上庫下泄或發電。以F1為例，t=0 h時，庫水位最高為413.5 m，對應庫區水量為4 264萬m3，大壩潰決，由零維-二維耦合計算模型計算得到下庫最大潰口流量為42 650 m3/s，潰壩流量過程見圖5。

上下庫依次潰壩過程中(F2和F3)，假定上庫溢洪道堵塞，壩前水位達最高水位時副壩瞬時全潰，受上庫潰壩洪水沖擊影響，下庫達最高水位時副壩逐漸潰決。以F3為例：t=78 h，庫水位最高達到416.08 m，對應庫區水量為4 723萬m3，大壩潰決，使用BREACH模型計算得到下庫最大潰口流量為3 082 m3/s，潰壩流量過程見圖6。不同潰壩場景下的潰口流量計算結果見表2。

圖6 汛期下庫主壩逐漸潰決出入庫洪水過程(F3)示意

表2 各潰壩工況潰口流量計算結果

3.4.2潰壩洪水演進模擬結果

由上述潰口流量計算模型及水庫下游一二維耦合水動力模型，對不同潰壩場景下梅蓄電站潰壩洪水演進過程進行模擬。計算結果表明，潰壩后各特征點的水位過程線總趨勢是靠近壩址處波峰陡峭，沿流程逐漸坦化；由于潰壩洪水量級遠超河道正常行洪能力，下游兩岸地勢較陡的河段，地形對河道洪水擴散的約束作用明顯，導致河道洪水位漲落幅度大；地勢平坦區，潰壩洪水漫過河道向兩側擴散顯著，減小了洪水沿河道向下游傳播的強度和速度，其水位漲落幅度相對偏小。以5 000 a一遇上下庫副壩依次潰壩方案(F3)為例，當發生5 000 a一遇洪水過程時，上庫豎井式溢洪道堵塞，無法正常下泄，庫內水位上漲，水位達到最高818.40 m時副壩瞬間全潰，瞬時最大潰口流量910 m3/s，受潰口流量的作用，壩前水位迅速下降至最低。而受上水庫急速下落的潰壩洪水影響，下庫水位最高達到416.08 m，對應庫區水量為4 723萬m3，下水庫最大潰壩流量為3 082 m3/s。下庫潰壩后，下游河道各特征點的水位過程見圖7，由結果可知，各特征點的水位過程趨勢基本一致，在潰壩發生后水位突然上升，且距離水庫越近上升幅度越大，模型合理反映了潰壩洪水在河道中的傳播過程。先水河、琴江干流水位包絡線示意見圖8、圖9，由圖可知，該工況下，先水河中下游地區存在大片淹沒，尤其是車前寨、榕溪村、硿下等河道兩岸地勢較低的村社區域，淹沒較大，部分淹沒水深達到8 m以上。對于琴江，在河道兩岸的河灘地和低洼地區存在一定的淹沒情況，主要體現在龍村鎮、汕湛高速南北、安流鎮、橫陂鎮以及五華縣城下游地區。

圖7 特征點的水位過程計算結果(F3)示意

圖8 先水河最高水面線成果(F3)示意

圖9 琴江最高水面線成果(F3)示意

3.4.3洪水風險分析成果

本模型合理模擬了潰壩洪水傳播過程，相關計算結果符合洪水演進規律。依據模型計算結果，各潰壩方案場景下計算結論和洪水影響見表3所示。在此基礎上，文章建立了不同潰壩場景下的潰壩風險圖，對梅蓄電站潰壩洪水到達時間、淹沒水深、淹沒范圍及洪水影響等過程特征進行預測，以直觀展示電站潰壩后工程下游區域受淹沒影響情況和洪水災害程度，評估洪水災害可能造成的損失，更好地制定防災減災應急響應對策和措施，由于篇幅受限，文章不一一列出。

表3 各潰壩方案工況計算結論及洪水影響

當先水河流域發生超標準洪水導致梅蓄電站水庫潰壩時，潰壩洪水所影響的區域主要集中在先水河沿河的村莊，遭遇琴江干流洪水后，淹沒總面積呈現隨潰壩洪水量級、琴江干流洪水增大而增大的趨勢，電站下游琴江干支流兩岸龍村、河東、水寨、橫陂四鎮為主要的潰壩洪水受災區。以F6為例，最大淹沒水深風險示意見圖10。潰壩洪水總淹沒面積為67.36 km2，淹沒深度在1 m以上的面積占比達68%，洪水所到之處樹木房屋倒塌較多，受淹房屋35 000余棟，受影響人口175 299人，龍村、河東、水寨、橫陂四鎮中心遭受毀滅性打擊，受影響人口分別為31 164人、23 065人、61 291人、44 222人。此外，水庫潰壩后，下游G238國道、G325國道、S19興汕高速、G35濟廣高速等基礎設施嚴重受損，交通中斷，給先水河及琴江干流兩岸村莊、基礎設施等造成嚴重后果，對沿岸企業、居民生產生活造成較大影響。

4 應急響應措施

由上述梅蓄電站潰壩洪水計算分析成果可知，梅蓄電站在超標準洪水、地震等極端工況下潰壩后，將嚴重威脅下游人民群眾的生命和財產安全。為減少災害損失和保障人民群眾的生命財產安全，需要因地制宜地采取切實有效的應急響應措施。

1)預先分級，確定應急響應啟動條件

根據梅蓄電站上下水庫實際運行情況，將可能出現的超標準洪水、重大工程隱患、地震、水庫運行管理不當以及戰爭或恐怖襲擊等重大險情，按照險情的嚴重程度和影響范圍，以及可能發生的概率，從高到低分為四級：特別重大(Ⅰ級)、重大(Ⅱ級)、較大(Ⅲ級)和一般(Ⅳ級)，同時確定應急響應的啟動條件。

2)事先防范，做好險情監測及預警

梅蓄公司工作人員時刻監視上下水庫水位變化和降雨量，同時及時關注梅州市五華縣氣象臺發布的有關強降雨信息，及時準確地進行洪水預報，做好水雨情監測；加強對梅蓄電站大壩及各水工建筑物日常監測和定期巡視檢查，及時掌握大壩運行動態。根據預報的洪水量級和地震等級不同，啟動不同等級的預警，由電站防汛辦公室以各種通訊手段向公眾發布相應等級的預警信息，及時上報主管部門及地方政府等。

3)統一指揮，做好組織保障

梅蓄電站成立由應急指揮中心、應急辦公室、現場應急處置工作組組成的應急組織機構。發生突發事件時，由應急指揮中心統一指揮應急處置工作，統籌調配人力、物力資源，協調各相關部門開展緊急救援工作。應急辦公室負責收集、分析、傳遞和發布突發事件信息，并確保相關組織部門間聯絡通暢。現場應急處置工作組以現場處置組為核心，技術支持組、綜合后勤保障組及安全監督組各保障小組提供支持，負責落實梅蓄公司應急指揮中心的相關決策和要求，做好現場應急處理和救援工作。

4)應急處置，實施搶險救援

按照梅蓄電站工程險情分級啟動相應應急搶險措施。當遭遇超標準洪水或其他重大險情，可能發生水庫潰壩時，通過地方防汛搶險指揮部并協調當地政府，根據梅蓄電站潰壩洪水避險轉移圖，將工程下游沿河兩岸處在潰壩洪水危險區內的居民進行轉移，按照就近、安全等原則轉移到達相應的安置點，重點轉移居民區為龍村、河東、水寨、橫陂四鎮。安置點兼顧行政隸屬，以學校、政府單位或地勢較高的廣場、公園等作為下游群眾和企事業單位的主要避護場所，并按照“撤得出、住得下、不挨餓、能救援”的原則進行布置，確保人員安全和基本生活。潰壩事件發生后，及時進行災后恢復和重建工作，必要時進行損失賠償[11-14]。

5 結語

文章以梅州抽水蓄能電站為例，對潰壩洪水模擬方法進行研究和分析，采用零維水庫耦合二維河道地表網格法進行混凝土壩瞬時潰決流量潰口計算，采用BREACH模型法進行土石壩逐漸潰決潰口流量計算，建立基于Godunov格式有限體積法的一二維耦合的潰壩洪水演進模型，合理模擬了梅蓄電站單庫潰壩、上下水庫連環潰壩等工況下的洪水演進過程，分析了潰壩洪水風險并繪制洪水風險圖，在此基礎上制定了科學合理的應急響應措施。文章對抽水蓄能電站多種可能潰壩因素和場景進行了綜合考慮，具有一定的推廣應用價值，一些潰壩洪水模擬的關鍵問題的處理對其他類似潰壩過程洪水模擬也具有一定借鑒意義。