曾明 陳瑜彬 鄒冰玉

摘要：“11·3”金沙江上游白格堰塞湖的潰決洪水破壞力強(qiáng)，屬于非常規(guī)、超歷史洪水，其中白格堰塞湖奔子欄江段洪水重現(xiàn)期超萬年。在缺乏高洪水文、河道地形資料的背景下，長江水利委員會水文局較準(zhǔn)確地預(yù)報了潰決洪水向下游演進(jìn)的沿程洪水特征，為應(yīng)急處置提供了決策支持。通過介紹此次堰塞湖潰決洪水演進(jìn)預(yù)報所采用的技術(shù)手段，歸納總結(jié)了預(yù)報經(jīng)驗和預(yù)報模型技術(shù)的優(yōu)缺點，以期為今后堰塞湖突發(fā)事件的應(yīng)急處置提供參考。

關(guān)鍵詞：潰決洪水;洪水演進(jìn);洪水預(yù)報;白格堰塞湖;金沙江

中圖法分類號：P338

文獻(xiàn)標(biāo)志碼：A

堰塞湖主要是由滑坡、泥石流等堵塞河道所形成的湖泊。堰塞湖形成后，水量不斷積蓄造成較高水頭差，潰決時由巨大能量形成的洪水常具有突發(fā)性強(qiáng)、流速快等特點。洪水特征主要表現(xiàn)為洪峰高、水量集中、洪水漲落迅速、傳播快等，對下游沿岸破壞力極強(qiáng)，遠(yuǎn)大于一般暴雨洪水，且傳播時間較一般洪水明顯縮短。準(zhǔn)確預(yù)報堰塞湖潰決洪水在下游沿程演進(jìn)的洪峰水位、流量以及峰現(xiàn)時間，可為提前轉(zhuǎn)移下游沿岸居民、保障沿岸村落安全提供重要指導(dǎo)依據(jù)。

“11·3”金沙江白格堰塞湖潰決洪水期間，長江水利委員會水文局根據(jù)實時報汛信息，基于模型適用性分析，不斷優(yōu)化參數(shù)，采用多模型計算成果比對，通過實時校正及滾動預(yù)報，準(zhǔn)確預(yù)報了堰塞湖下游沿程各站洪峰水位流量、起漲時間以及峰現(xiàn)時間，并向各級部門及時發(fā)布預(yù)報成果，為洪水預(yù)警、水文應(yīng)急監(jiān)測、人員轉(zhuǎn)移、水庫調(diào)度提供了重要技術(shù)支撐。

本文通過對比分析“11·3”堰塞湖應(yīng)急處置期間采用的洪水預(yù)報技術(shù)手段，探討總結(jié)了各預(yù)報模型的適用性和優(yōu)缺點，以期為今后應(yīng)急水文預(yù)報技術(shù)的選擇提供參考。

1堰塞湖及下游洪水演進(jìn)概況

2018年11月3日17：00許，西藏自治區(qū)昌都市江達(dá)縣波羅鄉(xiāng)白格村境內(nèi)金沙江右岸再次發(fā)生大規(guī)模山體滑坡，滑坡堵塞金沙江并形成堰塞湖。堰塞體上距崗?fù)兴恼炯s90km，距波羅站約15km，下距葉巴灘站約54km，距巴塘水文站約190km。由于堰塞體高、湖內(nèi)滯蓄水量大，潛在危險大，故采用開挖泄流槽方式對堰塞湖進(jìn)行引流處置，以降低堰塞湖的威脅。泄流槽于11月12日上午開始過流，經(jīng)過不斷的掏岸、溯源沖刷，過流斷面不斷擴(kuò)大，過流流量逐步增加，11月13日18：00過流洪峰流量達(dá)到31000m3/s。

此次潰壩洪水為非常規(guī)、超歷史洪水，對下游破壞力強(qiáng)。下游沿程巴塘、奔子欄站洪水重現(xiàn)期均超萬年，塔城站重現(xiàn)期超千年，然而演進(jìn)至石鼓站時洪水坦化明顯，轉(zhuǎn)為常遇洪水，其后經(jīng)金沙江中游梯級水庫調(diào)蓄，潰決洪水威脅得以解除。下游沿程各主要斷面洪峰特征值見表1。

2潰壩洪水演進(jìn)預(yù)報技術(shù)

白格堰塞湖形成于金沙江上游，地處高原山區(qū)，河道比降較大，河谷深切，斷面形態(tài)較穩(wěn)定。堰塞湖形成和潰決期間，流域內(nèi)無明顯降雨過程，區(qū)間來水穩(wěn)定，下游洪水演進(jìn)預(yù)報主要采用水力學(xué)模型、匯流曲線演算法、馬斯京根分段連續(xù)演算法進(jìn)行計算及綜合分析。還包括對下游沿程巴塘、奔子欄、石鼓、梨園水庫作了入庫預(yù)報。本文分述了3種河道演算方法在“11·3”金沙江白格堰塞湖下游洪水預(yù)報中的應(yīng)用。

2.1預(yù)報模型原理

2.1.1水力學(xué)模型

水動力學(xué)演算法主要基于圣維南（Saint-Venant）方程組的聯(lián)立求解，MIKE11模型可實現(xiàn)一維河道的水動力學(xué)模型演算，并可對連續(xù)急變流進(jìn)行較好的模擬。

MIKE11模型是基于一維非恒定流圣維南方程組來模擬河流或河口的水流狀態(tài)。方程組具體形式為式中，x、t分別表示空間坐標(biāo)和時間坐標(biāo);Q為斷面流量，m3/s;h為水位，m;A為斷面過流面積，m2;R為水力半徑，m;B.為河寬，m;g為單位河長的旁側(cè)人流量，M3/S;C為謝才系數(shù);a為動量修正系數(shù)。

MIKE11采用的有限差分格式為六點中心的Abbott-Ionescu格式。將圣維南方程組連續(xù)方程、動量方程分別寫成以水位點h和流量點Q為中心的形式進(jìn)行離散，采用追趕法對方程進(jìn)行求解。

2.1.2馬斯京根分段連續(xù)演算法

馬斯京根分段連續(xù)演算法是將演算河段劃分為Ⅳ個單元，經(jīng)過馬斯京根分段演算求得出流過程。計算參數(shù)有演算段數(shù)（單元河段數(shù)）Ⅳ，每個單元河段的馬斯京根法參數(shù)Xe與Ke。該演算法的具體公式為式中，Qi，j是第i河段末第j時段末的流量，m3/s;Qi-l，j是第i-1河段末第/時段末的流量，m3/s;Qi-1，j-1是第i-1河段末第j-l時段末的流量，m3/s;Qi，j.是第i河段末第j-l時段末的流量，m3/S;K為蓄量常數(shù);x為流量比重因素。

2.1.3 匯流曲線演算法

匯流曲線演算法是在加里寧河槽匯流曲線概念的基礎(chǔ)上，將流域（河網(wǎng)或河槽）沿程滯蓄作用概化為一系列等效的線性水庫和線性渠道，由連續(xù)方程和概化的蓄泄方程聯(lián)立推導(dǎo)演算出口流量的過程。流量演算采用卷積法，匯流曲線公式為式中，k為消退系數(shù);n為階數(shù)，即概化矩形人流的個數(shù);t為時段數(shù)。

2.2模型構(gòu)建

模型構(gòu)建范圍均為白格堰塞體至梨園水庫，考慮沿程區(qū)間人流，采用2018年“10·10”白格堰塞湖潰壩洪水演進(jìn)實況過程對預(yù)報模型進(jìn)行參數(shù)率定。各模型上邊界均采用堰塞體過流流量過程。

2.2.1水力學(xué)模型

由于研究區(qū)域無河道地形資料，水力學(xué)模型斷面采用葉巴灘、巴塘、奔子欄、石鼓水文站實測大斷面資料，對地形變化劇烈區(qū)域采用GIS等空間分析手段切割河道斷面，并結(jié)合專家經(jīng)驗進(jìn)行高程校核，同時對河段斷面形狀進(jìn)行插補(bǔ)處理。模型下邊界采用梨園水庫建庫前控制斷面水位流量關(guān)系，并對河道糙率進(jìn)行分段率定。

2.2.2 馬斯京根分段演算模型

馬斯京根分段連續(xù)演算法需率定的參數(shù)為河段單元數(shù)Ⅳ、單元演算系數(shù)Xe與Ke。分段演算模型將白格堰塞體至梨園水庫河段分為堰塞體一葉巴灘、葉巴灘一巴塘、巴塘一奔子欄、奔子欄一石鼓、石鼓一梨園5個計算區(qū)間，每個區(qū)間劃分子計算單元，采用金沙江“10 ·10”白格堰塞湖潰決過程試算得到子計算單元個數(shù)和單元演算系數(shù)。

2.2.3 匯流曲線演算模型

匯流曲線法根據(jù)匯流概化的條件不同，將匯流區(qū)域分段劃分為匯流單元，根據(jù)人出流洪水樣本分段率定，在確定各單元匯流曲線系數(shù)后，根據(jù)匯流曲線系數(shù)分段連續(xù)演算出流過程。根據(jù)報汛站網(wǎng)，將演算河段劃分為堰塞體一巴塘、巴塘一奔子欄、奔子欄一石鼓、石鼓一梨園4個計算單元。由于“10·10”白格堰塞湖潰壩洪水期間，巴塘一奔子欄河段存在蘇洼龍電站導(dǎo)流洞滯洪，為保障工程安全，“11·3”堰塞湖處置時將蘇洼龍電站拆除圍堰過流，使其基本接近天然河道。經(jīng)過分析巴塘一奔子欄匯流曲線參數(shù)，參考上河段并考慮洪水衰減坦化，綜合得出演算系數(shù)。

2.3預(yù)報檢驗及誤差分析

3種預(yù)報模型均采用堰塞體過流流量過程作為模型上邊界輸入進(jìn)行河道洪水演算。2018年11月13日18：00，推算堰塞體過流流量達(dá)到最大值后，依據(jù)水量平衡原理和“10·10”白格堰塞湖退水形態(tài)推求退水過程。采用3種模型并行演算下游洪水演進(jìn)，滾動實時校正，并綜合了3種模型計算結(jié)果。對下游巴塘、奔子欄、石鼓、梨園入庫起漲時間和洪峰水位、流量、峰現(xiàn)時間作了預(yù)報，預(yù)報誤差分析見表2。由表2可知：①巴塘、奔子欄站各項預(yù)報要素與實況基本一致。巴塘站峰現(xiàn)時間預(yù)報準(zhǔn)確，預(yù)報洪峰水位僅低0.01m，流量相對誤差0.5%。②奔子欄站預(yù)報峰現(xiàn)時間與實況一致，洪峰水位預(yù)報偏高0.02 m，流量相對誤差4.5%o。③石鼓、梨園入庫流量起漲時間較預(yù)報時間提前1～2h，峰現(xiàn)時間較預(yù)報提前8h左右，洪峰流量明顯偏大。主要是由于奔子欄以上沿程河段基本為穩(wěn)定深切斷面，洪水逐漸坦化但形態(tài)變化不大。奔子欄一石鼓河段存在多個灘地，“10·10”潰壩洪水量級不大，未達(dá)到漫灘水位;此次洪水峰高量大，該河段中出現(xiàn)多個灘地行洪，洪峰受灘地調(diào)蓄作用明顯坦化，洪水緩慢歸槽造成峰現(xiàn)時間較預(yù)報時間偏后。

采用3種模型對各站點流量過程模擬情況與實況對比見圖1。由圖可知，3種模型模擬結(jié)果與實況相比均有偏差，越向下游演進(jìn)，計算誤差越大，此次預(yù)報成果由3種模型結(jié)果綜合確定得出。通過比較，水力學(xué)模型模擬的峰現(xiàn)時間較實際略有提前，馬斯京根分段演算、匯流曲線演算模擬的峰現(xiàn)時間均偏后。原因為：①水力學(xué)模型沿程實測斷面較少，除地形變化劇烈處增加斷面外，其余斷面均為內(nèi)插，造成河道部分渠化，減小了河道調(diào)蓄能力，傳播加快;②馬斯京根分段演算法、匯流曲線演算法的模型參數(shù)綜合反映了河段傳播時間，率定采用的洪水樣本較此次潰壩洪水過程傳播時間長，造成這兩種模型計算的洪水時間較實際偏后。

2.4預(yù)報模型優(yōu)缺點比較

經(jīng)過此次堰塞湖非常規(guī)、超歷史洪水的預(yù)報實踐，3種洪水演進(jìn)模型在一定程度上能較好地反映出洪水演進(jìn)過程，3種模型方法的優(yōu)缺點對比分析如下。

（1）水力學(xué)模型屬于水力學(xué)方法，具有很強(qiáng)的物理意義，全動力波控制方程可滿足多種類型河道水力條件的計算需求，但模型資料要求高，沿程河道斷面地形資料是構(gòu)建模型、保障模型精度的基礎(chǔ)。在洪水過程模擬計算中，由于事發(fā)山區(qū)缺少實測資料，且高原峽谷區(qū)測量困難，對預(yù)報結(jié)果的精確度有一定影響。GIS與專家經(jīng)驗校正的切割斷面存在主觀誤差，造成奔子欄至石鼓地區(qū)未能較好反映灘槽調(diào)蓄作用。

（2）馬斯京根分段演算模型屬于水文學(xué)方法，對區(qū)間降雨徑流比重小的河道洪水演算具有很強(qiáng)的適應(yīng)性，但模型參數(shù)的率定需要多場次洪水資料，且河段傳播時間基本為綜合常數(shù)，對非常規(guī)洪水適用性一般。

（3）匯流曲線法同樣屬于水文學(xué)方法，與馬斯京根分段演算法相似，匯流曲線是多場洪水樣本的綜合參數(shù)，是洪水傳播時間的綜合反映，對傳播時間變化大的洪水適用性一般。

3結(jié)語

本文以“11·3”白格堰塞湖為例，分析了應(yīng)急處置時采用的預(yù)報技術(shù)手段和綜合預(yù)報過程。實踐檢驗表明，此次采用的3種洪水演進(jìn)預(yù)報模型，對分析計算上游山區(qū)河道的非常規(guī)超歷史洪水具有較高的準(zhǔn)確性，參數(shù)率定個數(shù)較少，能基本滿足突發(fā)事件下的洪水預(yù)報需求。其中，水力學(xué)模型能較為準(zhǔn)確地綜合考慮河道地形和洪水演進(jìn)坦化變形影響。由于各模型自身物理意義不同，各有優(yōu)缺點和適用范圍，應(yīng)急預(yù)報時仍需綜合考慮各模型結(jié)果的差異性和可靠性，多模型比較綜合預(yù)報成果，避免因采用單一模型計算而造成的較大誤差。

此次堰塞湖發(fā)生在山區(qū)，山區(qū)性河流漲落快速、水位流量關(guān)系單一，預(yù)報采用的均為適用山區(qū)性河道的洪水演進(jìn)模型。突發(fā)事件的應(yīng)急處置需快速響應(yīng)，未來應(yīng)整理歸類不同特性河道（如平原河網(wǎng)、湖泊、水庫等）和流域的應(yīng)急預(yù)報模型以及開展模型適用性等方面的研究工作，以便為堰塞湖的應(yīng)急處置、人員轉(zhuǎn)移避災(zāi)提供準(zhǔn)確及時的技術(shù)支撐。

參考文獻(xiàn)：

[1]康琳，朱軍，李維煉.基于潰壩洪水模型的山區(qū)堰塞湖避難場所選址[J].自然災(zāi)害學(xué)報，2018，27（5）：39-45.

[2]王云南，任光明，夏敏，等.滑坡堰塞壩穩(wěn)定性研究綜述[J].中國地質(zhì)災(zāi)害與防治學(xué)報，2016，27（1）：6-14.

[3]鄧宏艷，孔紀(jì)名，王成華，等，不同成因類型堰塞湖的應(yīng)急處置措施比較[J]山地學(xué)報，2011，29（4）：505-510.

[4]謝俊舉，李亞琦，基于洪水演進(jìn)數(shù)值模擬的潰壩危害性快速評估[J].災(zāi)害學(xué)，2011，26（2）：31-34.

[5]張莉.潰壩洪水特征及計算方法探析[J].地下水，2018，40（4）： 203-205.