新疆射月溝水庫潰壩過程數(shù)值模擬

馬 黎，鐘啟明，楊 蒙，吳 迪，李大成，梅勝堯

(1.中國電建集團(tuán)貴陽勘測設(shè)計(jì)研究院有限公司，貴州 貴陽 550081； 2.南京水利科學(xué)研究院，江蘇 南京210029； 3.水利部水庫大壩安全重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210029)

近年來，受全球氣候變化影響，我國極端暴雨事件呈增加趨勢[1]。2000 年1 月至2021 年12 月，我國共發(fā)生99 起潰壩事故，多起潰壩事故發(fā)生在干旱和半干旱地區(qū)，其中新疆和內(nèi)蒙古自治區(qū)分別發(fā)生9 起和6 起由極端暴雨導(dǎo)致的潰壩事故[2-3]。如2018 年7 月31 日，位于新疆自治區(qū)哈密市的射月溝水庫大壩因壩頂溢流發(fā)生潰決[4]。合理預(yù)測潰壩過程對潰壩洪水致災(zāi)后果的評價(jià)具有重要意義，采用數(shù)值分析方法模擬潰口流量過程和潰口尺寸演化過程是預(yù)測潰壩過程的有效手段。目前常用的潰壩過程數(shù)學(xué)模型可分為經(jīng)驗(yàn)?zāi)Ｐ秃突跐Q機(jī)理的模型[5-7]。經(jīng)驗(yàn)?zāi)Ｐ鸵话阃ㄟ^考慮壩體形態(tài)、水庫特征、壩料特性等信息，基于回歸分析或機(jī)器學(xué)習(xí)方法建立經(jīng)驗(yàn)計(jì)算式，獲取潰壩時(shí)的潰口峰值流量、最終潰口尺寸和潰壩歷時(shí)等潰壩參數(shù)[5]。但經(jīng)驗(yàn)?zāi)Ｐ蜔o法獲取潰壩時(shí)的潰口流量過程和潰口發(fā)展過程?；谕潦瘔螡Q機(jī)理的潰壩過程數(shù)學(xué)模型[8-14]，一般采用寬頂堰流或孔流公式模擬漫頂或滲透破壞潰壩時(shí)的潰口流量，并利用各種沖蝕公式來描述壩料的輸移，一些模型還可以考慮潰口邊坡的失穩(wěn)[13-14]，通過按時(shí)間步長迭代的數(shù)值計(jì)算方法獲取各時(shí)間步長的潰口流量和潰口尺寸。

本文通過對射月溝水庫潰壩案例的調(diào)研，總結(jié)潰壩水流的沖蝕特性和潰口發(fā)展過程。在此基礎(chǔ)上建立基于壩體潰決機(jī)理的簡化漫頂潰壩過程數(shù)學(xué)模型，反演射月溝水庫的潰壩過程，并與國外典型的NWS BREACH 模型的計(jì)算結(jié)果進(jìn)行對比分析。

1 射月溝水庫潰壩事件調(diào)研

射月溝流域地處哈爾里克山南坡，由6 條支流匯流形成射月溝水庫（圖1）[15]。射月溝水庫集水面積約408 km2，徑流補(bǔ)給主要由季節(jié)性冰雪融水、夏季降雨和泉水組成，年平均降雨強(qiáng)度約為150 mm，年平均徑流量為1.58×107m3。

圖1 射月溝流域水文概況Fig.1 Hydrological survey of Sheyuegou Basin

射月溝水庫大壩距離二宮村頭宮隊(duì)2 km，距離哈密市115 km。大壩為瀝青混凝土心墻砂礫石壩，3 層瀝青混凝土心墻厚度均為0.4 m，壩頂高程為1 496.65 m，防浪墻頂高程為1 497.65 m；大壩最大壩高41.15 m，壩頂長403 m、寬6 m，上、下游坡比（垂直/水平）分別為1∶2.25 和1∶2.00[16]，射月溝水庫大壩典型斷面見圖2。

圖2 射月溝水庫大壩典型斷面示意Fig.2 Schematic diagram of the typical section of Sheyuegou reservoir dam

1.1 2018 年“7·31”洪水產(chǎn)匯流過程

2018 年7 月31 日01:00 至13:00，射月溝流域發(fā)生特大暴雨，小堡、射月溝、秦城3 個(gè)氣象站12 h的降水強(qiáng)度分別為116、105 和79 mm（3 個(gè)氣象站位置見圖1），均打破了當(dāng)?shù)氐臉O值記錄[17]。2018 年7 月31 日22:00，新疆自治區(qū)水文局基于坡面面積法，通過現(xiàn)場洪水痕跡調(diào)查，得出射月溝水庫入庫洪峰流量達(dá)到1 848 m3/s，遠(yuǎn)超設(shè)計(jì)標(biāo)準(zhǔn)，導(dǎo)致水庫發(fā)生漫頂潰壩。趙懸濤等[15]通過模擬獲取了射月溝水庫2018 年7 月31 日的降雨-入庫流量過程（圖3），入庫峰值流量相對誤差僅為3.6%，在本文的潰壩過程模擬中采用此入庫洪水流量過程線。

圖3 射月溝水庫2018 年7 月31 日降雨-入庫流量過程線Fig.3 The rainfall-inflow hydrograph of Sheyuegou reservoir on July 31, 2018

1.2 射月溝水庫漫頂潰壩過程

2018 年7 月31 日8:00 至11:40，一名目擊者用手機(jī)錄制了射月溝水庫潰壩過程的多段視頻。雖然視頻是臨時(shí)拍攝的，但獲取的信息對研究人員直觀了解壩體沖蝕特性和潰口形態(tài)演變過程具有重要意義。

根據(jù)現(xiàn)場視頻，可將射月溝水庫潰壩過程分為4 個(gè)階段：（1）初始沖坑出現(xiàn)。水流越過壩頂后表現(xiàn)為全斷面溢流，由于壩頂處水流流速較小，頂部未出現(xiàn)明顯沖蝕；受重力作用，水流在下游坡行進(jìn)時(shí)流速逐漸增大，導(dǎo)致下游坡某一位置的土體發(fā)生沖蝕，形成初期沖坑（圖4(a)）。（2）溯源沖蝕。在漫頂水流的作用下，沖坑以上下游坡逐漸變陡，直至接近垂直；隨后發(fā)生溯源沖蝕，在此過程中下游坡面幾乎保持垂直，沖蝕不斷向上游發(fā)展（圖4(b)～(c)）。（3）初始潰口出現(xiàn)。當(dāng)溯源沖蝕發(fā)展至上游水庫后，由于壩頂高程突然下降，初始潰口出現(xiàn)，漫頂水流匯聚到潰口處，其他部位的漫頂沖蝕結(jié)束（圖4(d)）。（4）潰口擴(kuò)展。隨著漫溢水頭的逐漸增大，潰口在橫斷面和縱斷面的沖蝕逐漸加快；由于潰口形態(tài)的快速變化，潰口邊坡發(fā)生失穩(wěn)；當(dāng)潰口處水流的剪應(yīng)力小于壩料的臨界啟動剪應(yīng)力時(shí)，沖蝕過程停止（圖4(e)～(f)）。

圖4 射月溝水庫潰壩過程照片F(xiàn)ig.4 Photos of the Sheyuegou dam breach process

根據(jù)潰壩視頻和現(xiàn)場調(diào)查，大壩內(nèi)部的瀝青混凝土心墻無法抑制潰口的發(fā)展，其原因主要有兩個(gè)方面：一是瀝青混凝土心墻均位于壩體上游部位（圖2），主要起防滲作用，對初期潰口的發(fā)展影響微弱；二是溯源沖蝕發(fā)展至瀝青混凝土心墻位置時(shí)，潰口處的漫溢水頭已達(dá)到較大值，僅有40 cm 厚度的瀝青混凝土心墻無法阻止?jié)⒖诘陌l(fā)展。因此在數(shù)值模擬中不考慮心墻。

2 基于潰決機(jī)理的簡化漫頂潰壩過程數(shù)學(xué)模型

基于射月溝水庫漫頂潰壩過程的視頻和現(xiàn)場調(diào)查資料，建立基于物理機(jī)制的簡化潰壩數(shù)學(xué)模型。由于壩體潰決時(shí)間相對較短，滲流的影響可以忽略不計(jì)。模型主要包括4 個(gè)部分：潰口流量計(jì)算、初始沖坑位置確定、縱斷面潰口發(fā)展和橫斷面潰口發(fā)展。

2.1 潰口流量計(jì)算

水量平衡是潰壩模擬中必須遵循的基本原理，利用水量平衡關(guān)系描述水庫水位變化關(guān)系：

式中：As為與庫水位相對應(yīng)的庫面面積；zs為庫水位；t為時(shí)間；Qin為入庫流量；Qb為潰口流量；Qs為溢洪道泄流量；QT為閘門泄流量。

潰口流量可采用寬頂堰流量方程[18]表示：

式中：c1和c2分別為流量系數(shù)，其中c1= 1.7 m0.5/s，c2= 1.3 m0.5/s[18]；b為潰口底寬；H=zs－zb，其中zb為潰口底高程；m為潰口邊坡坡比（垂直/水平）的倒數(shù)；ksm為尾水修正系數(shù)[18-19]。

2.2 初始沖坑位置確定

根據(jù)相關(guān)文獻(xiàn)[20-21]，土石壩潰壩時(shí)初始沖坑的位置計(jì)算如下：

式中：ln為從下游坡頂順坡至初始沖坑的距離；β0為大壩下游坡角；dn為下游坡法向水深；Frn為初始沖坑處的弗勞德數(shù)；Un為水流在初始沖坑處及越過沖坑后的法向平均流速，Un=C(Rnsinβ0)1/2=C(Hsinβ0)1/2，其中C為謝才系數(shù)，Rn為水力半徑；g為重力加速度；Bw為壩頂處漫頂水流的頂部寬度；Bn為下游壩坡處漫頂水流的頂部寬度；當(dāng)漫頂水流越過初始沖坑后，dn=Qb/(UnBn)。由于潰壩第一階段為全斷面溢流，因此假設(shè)Bw=Bn。

2.3 縱斷面潰口發(fā)展

根據(jù)射月溝潰壩案例的視頻記錄，在數(shù)學(xué)模型中可將溯源沖蝕過程分為兩個(gè)階段：（1）初始沖坑出現(xiàn)后，水流繼續(xù)掏刷沖坑，下游壩坡坡角逐漸變大至基本垂直，形成陡坎；（2）陡坎保持直立狀向上游發(fā)生運(yùn)移，直至進(jìn)入水庫。

潰口的底床沖蝕率采用基于剪應(yīng)力的沖蝕公式[22]表示：

式中：E為底床沖蝕率；kd為壩料沖蝕系數(shù)，可通過試驗(yàn)測量[23]或經(jīng)驗(yàn)公式計(jì)算[24]，即kd= 10(ρw/ρd)exp，其中ρw為水的密度取1 000 kg/m3，ρd為壩料的干密度，cr為壩料的黏粒含量；τb為潰壩水流的剪應(yīng)力，可通過曼寧公式計(jì)算[10]；τc為壩料顆粒的臨界起動剪應(yīng)力，可由希爾茲曲線[25]獲取。

對于下游壩坡變陡階段，坡角增量計(jì)算如下：

式中：dβ/dt為下游坡角變化率；E1為沖坑處的底床沖蝕率；E0為壩頂處的底床沖蝕率。

對于陡坎運(yùn)移階段，利用時(shí)間平均的沖蝕率公式來計(jì)算溯源沖蝕率[26]：

式中：dx/dt為陡坎運(yùn)移速率；CT為陡坎運(yùn)移系數(shù)；q為潰口單寬流量；He為陡坎高度。

隨著陡坎的運(yùn)移，陡坎上游側(cè)的壩體在縱斷面方向逐漸變薄，在水動力作用下可能沿水平面破壞，可通過比較壩體在水平破壞面上的驅(qū)動力和抵抗力來評價(jià)陡坎上游側(cè)壩體的穩(wěn)定性。驅(qū)動力由上游水壓力和陡坎底部的水流剪應(yīng)力組成，而抵抗力則包括水平破壞面和陡坎兩側(cè)的黏聚力和摩擦力及尾水引起的水壓力。當(dāng)驅(qū)動力大于抵抗力時(shí)，陡坎上游側(cè)壩體會發(fā)生失穩(wěn)破壞，并假定失穩(wěn)的壩體物質(zhì)在瞬間被漫頂水流沖走，至此溯源沖蝕發(fā)展至水庫。

2.4 橫斷面潰口發(fā)展

當(dāng)溯源沖蝕發(fā)展至水庫后，漫頂水頭突然增大，壩體橫斷面出現(xiàn)初始潰口，根據(jù)射月溝潰壩視頻假設(shè)初始潰口的形狀為矩形。隨著水流的不斷侵蝕，初始潰口迅速加深加寬，當(dāng)潰口深度達(dá)到臨界高度時(shí)，潰口邊坡沿平面發(fā)生失穩(wěn)破壞，初始潰口由矩形變?yōu)榈固菪?，此后漫頂水流繼續(xù)沖蝕潰口直至潰壩結(jié)束。

假設(shè)漫頂水流對潰口底部和側(cè)壁具有同等的沖蝕能力，則潰口頂寬增量為：

式中：ΔBt為潰口頂寬增量；nloc為潰口位置參數(shù)（nloc= 1 代表單側(cè)沖蝕，nloc= 2 代表雙側(cè)沖蝕）；α為潰口邊坡坡角。

潰口底寬增量表示為：

式中：ΔBb為潰口底寬增量。

利用極限平衡法評價(jià)潰口邊坡的穩(wěn)定性，當(dāng)失穩(wěn)塊體自重引起的驅(qū)動力大于失穩(wěn)塊體自重和土體黏聚力引起的抵抗力時(shí)，潰口邊坡發(fā)生失穩(wěn)破壞。視頻記錄和其他模型試驗(yàn)表明，邊坡失穩(wěn)時(shí)的破壞面近乎平面，并假定失穩(wěn)土體被潰壩水流瞬間沖走。潰口邊坡失穩(wěn)的驅(qū)動力Fd和抵抗力Fr分別為：

式中：W為失穩(wěn)塊體重力；θ為失穩(wěn)后的潰口邊坡坡角；γs為邊坡土體重度；Hs為潰口邊坡的臨界高度；φ為土體內(nèi)摩擦角；C為土體黏聚力。

采用按時(shí)間步長迭代的數(shù)值計(jì)算方法模擬潰壩過程中的水土耦合效應(yīng)，每個(gè)時(shí)間步長可輸出潰口流量、潰口尺寸、庫水位等潰壩參數(shù)。

3 射月溝水庫潰壩案例數(shù)值模擬

采用上述建立的數(shù)學(xué)模型對射月溝水庫潰壩過程進(jìn)行反演分析，通過視頻記錄和現(xiàn)場實(shí)測數(shù)據(jù)，驗(yàn)證模型的合理性，并對關(guān)鍵輸入?yún)?shù)進(jìn)行敏感性分析。

3.1 輸入?yún)?shù)

從拍攝的視頻中可以看到，射月溝水庫大壩于7 月31 日9:10 開始漫頂溢流，并在當(dāng)日11:00 左右潰壩結(jié)束。模擬持續(xù)時(shí)間為3 h，時(shí)間步長為1 s。根據(jù)視頻記錄，將潰壩計(jì)算時(shí)的初始水位設(shè)置為1 497.75 m，超過防浪墻頂高程0.1 m，入庫流量過程如圖3 所示。壩料沖蝕系數(shù)（kd）和陡坎運(yùn)移系數(shù)（CT）是潰壩模擬的重要指標(biāo)。根據(jù)設(shè)計(jì)報(bào)告，ρd=2 240 kg/m3，cr= 0.04，由前文經(jīng)驗(yàn)公式可得，kd=3.0 cm3/N/s。另外，基于潰壩過程視頻確定的陡坎移動時(shí)間，推導(dǎo)出CT= 0.001 /m1/6/s2/3。表1 給出了射月溝水庫和壩體的物理力學(xué)指標(biāo)[4,16]。

表1 模型輸入?yún)?shù)Tab.1 Input parameters of the numerical model

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

由于初始漫頂溢流的水位為0.1 m，由式（4）可以得出，初始沖坑位于下游頂部順坡向下8.93 m 處。計(jì)算結(jié)果顯示，陡坎在潰壩后0.90 h（對應(yīng)時(shí)間為10:04）進(jìn)入水庫，現(xiàn)場觀測時(shí)間為10:00，計(jì)算結(jié)果與實(shí)測數(shù)據(jù)吻合較好。另外，根據(jù)計(jì)算結(jié)果，在潰壩后1.38 h（對應(yīng)時(shí)間為10:33）出現(xiàn)潰口峰值流量，這也與實(shí)測數(shù)據(jù)10:25 出現(xiàn)峰值流量相吻合。

表2 給出了射月溝水庫潰壩案例的計(jì)算和實(shí)測潰壩參數(shù)，包括潰口峰值流量（Qp）、潰口最終頂寬（Btf）、潰口最終底寬（Bbf）、潰口最終深度（Df）和潰口峰值流量出現(xiàn)時(shí)刻（Tp）。從表2 可以看出，潰口峰值流量、潰口最終頂寬和底寬、潰口峰值流量出現(xiàn)時(shí)刻的計(jì)算誤差均較小。由于沒有潰口流量過程和潰口形態(tài)發(fā)展過程的實(shí)測數(shù)據(jù)，僅有推求的峰值流量和實(shí)測最終潰口尺寸，圖5 和圖6 主要展示了潰口流量過程和潰口形態(tài)發(fā)展過程的計(jì)算值。另外，從圖5 和圖6 的變化趨勢看出，壩體橫斷面和縱斷面潰口形態(tài)的演變趨勢與現(xiàn)場視頻記錄相吻合，這表明本文所提出的數(shù)值模擬方法能夠較為準(zhǔn)確地描述射月溝水庫的潰壩過程。

表2 射月溝水庫計(jì)算與實(shí)測潰壩參數(shù)值對比Tab.2 Comparison of numerical and measured breaching parameters of Sheyuegou reservior

圖5 潰口流量過程計(jì)算值與潰口峰值流量實(shí)測值Fig.5 Numerical and measured results of breach flow discharge

圖6 潰口形態(tài)發(fā)展過程計(jì)算值與潰口最終尺寸實(shí)測值Fig.6 Numerical and calculated results of breach development

3.3 參數(shù)敏感性分析

在射月溝水庫潰壩過程的數(shù)值模擬中，壩料沖蝕系數(shù)（kd）和陡坎運(yùn)移系數(shù)（CT）是具有高度不確定性的關(guān)鍵輸入?yún)?shù)。kd對潰壩過程中的沖蝕率有較大影響，是決定潰壩過程和潰口尺寸的關(guān)鍵參數(shù)；CT直接決定陡坎的運(yùn)移速率，也對潰壩過程具有重要影響。通過將數(shù)學(xué)模型中的相關(guān)參數(shù)（kd和CT）分別乘以0.5 和2.0，重新計(jì)算射月溝水庫的潰壩過程，分析模型分別對這2 個(gè)參數(shù)的敏感性。表3 和表4 分別顯示了在不同kd和CT條件下，模型計(jì)算出的潰口峰值流量、潰口最終尺寸（長度、寬度和深度）和潰口峰值流量出現(xiàn)時(shí)刻的變化情況。圖7 和圖8 還分別比較了不同kd和CT條件下的潰口流量過程。

表3 不同壩料沖蝕系數(shù)（kd）對潰壩過程的影響Tab.3 Influence of different dam material erodibility coefficients (kd) on dam breach process

表4 不同陡坎運(yùn)移系數(shù)（CT）對潰壩過程的影響Tab.4 Influence of different headcut migration coefficients (CT) on dam breach process

圖7 不同壩料沖蝕系數(shù)kd 條件下的潰口流量過程Fig.7 Breach hydrographs under different kd

圖8 不同陡坎運(yùn)移系數(shù)CT 條件下的潰口流量過程Fig.8 Breach hydrographs under different CT

對于沖蝕系數(shù)（kd），各潰壩參數(shù)的總體敏感性由強(qiáng)到弱依次為：潰口峰值流量，潰口峰值流量出現(xiàn)時(shí)刻，潰口最終尺寸。當(dāng)沖蝕系數(shù)增大1 倍或縮小至50%時(shí)，潰口峰值流量將分別產(chǎn)生57.5%和?41.7%的變化，而潰口峰值流量對潰壩致災(zāi)后果評價(jià)具有重要意義，因此在選取參數(shù)時(shí)應(yīng)注重其合理性。對于陡坎運(yùn)移系數(shù)（CT），由于原始模型中CT的取值較小，當(dāng)CT縮小至50%時(shí)，對計(jì)算結(jié)果的影響較小；當(dāng)CT增大1 倍時(shí)，各潰壩參數(shù)的敏感性有所增強(qiáng)，尤其是潰口峰值流量出現(xiàn)時(shí)間，較原始模型提前了23 min。各潰壩參數(shù)的總體敏感性由強(qiáng)到弱依次為：潰口峰值流量出現(xiàn)時(shí)刻、潰口最終尺寸、潰口峰值流量?？傮w而言，kd主要影響潰口峰值流量和潰口峰值流量出現(xiàn)時(shí)刻，而CT主要影響潰口峰值流量出現(xiàn)時(shí)刻。

4 與典型潰壩模型比較

選擇美國NWS BREACH 模型[8]模擬射月溝水庫的潰壩過程，然后通過計(jì)算結(jié)果比較其與本文模型的優(yōu)缺點(diǎn)。選擇NWS BREACH 模型的原因在于，該模型也是基于潰決機(jī)理的簡化漫頂潰壩過程數(shù)學(xué)模型，且目前在國際上應(yīng)用較廣泛。表5 列舉了NWS BREACH 模型和本文模型的特點(diǎn)。

表5 NWS BREACH 模型和本文模型的特點(diǎn)總結(jié)Tab.5 Summary of characteristics of NWS BREACH model and the present model

NWS BREACH 模型采用修正Meyer-Peter-Mueller 公式[27]模擬壩料的沖蝕過程，其中起控制作用的參數(shù)是d90/d30（d90和d30分別為小于該粒徑的顆粒質(zhì)量占總質(zhì)量30%和90%對應(yīng)的顆粒粒徑），反映了壩料顆粒的分布特征。根據(jù)射月溝壩料的顆分曲線[4]，d90/d30約為45。模型其他輸入?yún)?shù)可參考表1。

表6 對比了NWS BREACH 模型和本文模型預(yù)測的潰壩參數(shù)結(jié)果及相對誤差。圖9 為2 個(gè)模型預(yù)測的射月溝水庫潰壩過程中的潰口流量過程線。

表6 不同基于潰決機(jī)理的簡化潰壩模型計(jì)算結(jié)果對比Tab.6 Comparison of numerical results of different simplified physically-based dam breach models

圖9 不同基于潰決機(jī)理的簡化潰壩模型預(yù)測的潰口流量過程Fig.9 Predicted breach hydrographs of different simplified physically-based dam breach models

由表6 和圖9 可以看出，NWS BREACH 模型預(yù)測的潰壩參數(shù)和潰口流量過程與本文模型相比表現(xiàn)出更大的相對誤差。由于壩料的沖蝕特性對潰壩過程至關(guān)重要，而NWS BREACH 模型中決定沖蝕系數(shù)的關(guān)鍵參數(shù)是d90/d30，射月溝水庫的堆石材料表現(xiàn)出寬級配的特征，且d90/d30= 45，導(dǎo)致壩料的沖蝕率較大，因而預(yù)測得到的潰口峰值流量偏大，峰值流量出現(xiàn)時(shí)刻提前，潰壩歷時(shí)縮短。雖然NWS BREACH模型預(yù)測的潰口流量過程線呈現(xiàn)出較大的誤差，但由于潰壩歷時(shí)較短，潰口形態(tài)后期無明顯擴(kuò)展，因此潰口最終形態(tài)誤差較小。通過射月溝水庫潰壩實(shí)例的比較發(fā)現(xiàn)，本文模型的預(yù)測精度總體高于NWS BREACH 模型。

5 結(jié) 語

根據(jù)射月溝水庫潰壩過程的視頻記錄和現(xiàn)場調(diào)查，本文建立了基于潰決機(jī)理的描述水庫潰壩過程的簡化數(shù)學(xué)模型。利用射月溝水庫潰壩案例的實(shí)測資料對模型的準(zhǔn)確性進(jìn)行驗(yàn)證，并通過與其他典型潰壩模型預(yù)測結(jié)果的對比，證明了本文模型預(yù)測結(jié)果的合理性。主要結(jié)論如下：

（1）射月溝水庫的潰壩機(jī)制主要表現(xiàn)為縱斷面方向初始沖坑的形成與溯源沖蝕，以及橫斷面方向潰口的連續(xù)下切、側(cè)向沖蝕和潰口邊坡的間歇性失穩(wěn)坍塌。此外，隨著沖蝕的不斷加劇，縱斷面方向大壩下游坡角逐漸減小，而橫斷面潰口的邊坡角基本保持不變，直至發(fā)生潰口邊坡失穩(wěn)。

（2）參數(shù)敏感性分析表明，壩料沖蝕系數(shù)對潰壩過程有顯著影響，潰口峰值流量和潰口峰值流量出現(xiàn)時(shí)刻對壩料沖蝕系數(shù)的敏感性大于潰口最終尺寸，且潰口峰值流量對壩料沖蝕系數(shù)的敏感性最高。另外，潰口峰值流量出現(xiàn)時(shí)刻比其他潰壩參數(shù)對陡坎移動系數(shù)更為敏感。因此，在進(jìn)行潰壩模擬時(shí)，壩料沖蝕系數(shù)和陡坎移動系數(shù)的選取至關(guān)重要。

（3）因壩料沖蝕率公式的選擇和對縱、橫斷面潰口發(fā)展過程計(jì)算方法的過度簡化，目前廣泛應(yīng)用的NWS BREACH 未能合理描述射月溝潰壩過程的特點(diǎn)。但由于對土石壩潰決過程的認(rèn)識不足，且測量數(shù)據(jù)存在較大的不確定性，因此僅通過單個(gè)潰壩案例難以全面展示各模型的功能。本文對模型預(yù)測能力的對比僅作為參考，未來仍需對各模型的精度開展進(jìn)一步驗(yàn)證。