基于RAMMS的冰湖潰決型泥石流演進模擬及危害性

劉府生, 周瑞宸, 孫紅林, 胡卸文, 黎尤, 徐玉龍

(1. 中鐵第四勘察設計院集團有限公司, 武漢 430063; 2.西南交通大學地球科學與環(huán)境工程學院, 成都 610031)

由于溫室效應不斷積累導致全球氣候變暖,冰川消融速率呈明顯上升趨勢,高海拔或高緯度山區(qū)冰湖規(guī)模及數(shù)量急速增加,冰湖潰決事件高發(fā)頻發(fā)[1-2]。近50年喜馬拉雅地區(qū)至少記載有20次大型冰磧湖潰決災害事件,其中3/4發(fā)生在西藏[3]。在高原鐵路規(guī)劃建設的背景下,冰湖潰決型泥石流作為目前西藏境內(nèi)最主要的地質(zhì)災害之一,研究其形成演化特征具有重要意義[4-5]。西藏冰湖潰決泥石流主要分布于藏東南地區(qū),其造成的后果往往比常見的暴雨泥石流更為嚴重,不僅直接對下游的房屋建筑、公路、橋梁等工程及人員安全造成巨大危害,還極易誘發(fā)沿岸崩塌、滑坡等次生災害[6-7]。例如,1981年7月11日聶拉木縣樟藏布溝發(fā)生特大規(guī)模冰湖潰決型泥石流,致使中尼公路聶友段多數(shù)公路及橋梁被毀,200余人罹難[8]。2013年嘉黎縣然則日阿錯發(fā)生潰決,導致周邊居民和基礎設施不同程度受災,直接經(jīng)濟損失達2.7億元[9]。2020年6月嘉黎縣吉翁錯冰湖發(fā)生潰決,泥石流淤積溝口、堵塞主河道,對下游村鎮(zhèn)構(gòu)成嚴重威脅[10]。

冰湖潰決型泥石流通常暴發(fā)突然、能量巨大,具有持續(xù)時間短、洪峰流量大、波及范圍廣等特點。國內(nèi)外許多學者對冰湖潰決的致災誘因及機理進行了研究,總的來說分為兩類:第一類是冰崩(冰滑坡)、強降水、地震等外部誘因造成水位上漲、激發(fā)涌浪,引起冰磧湖自身失衡而導致的冰湖潰決;第二類是冰磧壩內(nèi)死冰消融、堤壩管涌擴大等內(nèi)部誘因,導致壩體失效、引發(fā)洪水,沿途沖刷侵蝕形成冰湖潰決型泥石流[11-12]。冰湖潰決誘因眾多而復雜,且災害一旦發(fā)生后果嚴重,因此對冰湖潰決進行風險評價十分重要,這也是冰湖潰決型泥石流演進模擬的必要前提。冰湖潰決風險評價方法大致分為定性、半定量和定量三類,隨著研究的深入與技術的進步,逐步向評價指標定量化的方向發(fā)展。根據(jù)歷史上已潰決冰磧湖的不同特點,國內(nèi)外學者提出了相應的冰湖潰決風險評價指標體系。Mckillop等[13]以遙感數(shù)據(jù)為基礎,根據(jù)冰磧湖、冰磧壩、母冰川和冰湖盆的特點及其相互關系歸納出18項評價指標,并篩選出4個參數(shù)、基于邏輯回歸方法建立了冰湖潰決風險概率方程。Ermini等[14]基于對84座堰塞壩的統(tǒng)計分析(阿爾卑斯山脈和亞平寧山脈36座,日本17座,美國和加拿大20座,新西蘭和印度等國家11座),選取壩體體積、流域面積和壩高作為評價因素,提出采用無量綱堵塞指數(shù)(dimensionless blockage index,DBI)來評估堰塞壩穩(wěn)定性。位宏等[15]選取母冰川、冰湖和終磧壩3類指標共8個因子構(gòu)建了新疆冰湖潰決風險模型。余斌等[16]基于冰磧湖潰決機理,提出采用6個無量綱影響因子的冰湖潰決易發(fā)性評價模型。

冰湖潰決數(shù)值模擬旨在通過數(shù)學公式或物理模型估算潰決洪水災害動力學特征參數(shù)、重現(xiàn)或預測冰湖潰決演進過程及下游危險性范圍。Wang等[17]利用BREACH模型對西藏龍巴薩巴湖和皮達湖進行潰壩風險評價及兩湖全潰模擬分析,并驗證了該模型的可靠性。韓松林等[18]基于HEC-RAS一維非恒定流水動力模型,建立了上游冰湖至下游河道的洪水演進數(shù)學模型,模擬了桑旺錯潰壩、洪水下泄的演進過程。劉建康等[19]通過RAMMS數(shù)值模擬分析了西藏波密丹卡弄巴冰湖在全潰條件下的演進過程及致災風險。劉波等[20]利用RAMMS內(nèi)置泥石流模塊的Vollemy-Salm單相流模型,模擬了西藏洛隆巴曲冰湖潰決型泥石流在4種極端工況下的演進過程。王翔等[21]采用FLO-2D水文動力學模型對藏東南地區(qū)典型冰湖強宗克措和吉萊普措進行潰決泥石流模擬,為下游地區(qū)災害防治措施提供依據(jù)。

總的來說,國內(nèi)外學者有關冰湖潰決洪水及泥石流研究取得了較好的進展,但由于基礎資料的缺失及調(diào)查方法的局限,對冰湖潰決風險評價仍缺乏完善的體系[22-23]。現(xiàn)以西藏洛隆凍錯曲冰湖為例,在查明泥石流孕災條件的基礎上,對泥石流的動力學特征參數(shù)及進行計算;采用無量綱堵塞指數(shù)(DBI)對冰湖堰塞體進行潰決風險判別;采用經(jīng)驗公式分別計算2種潰決模式下的潰口峰值流量;基于三維動態(tài)模擬軟件RAMMS分析冰湖泥石流潰決演進過程,定量評價潛在威脅區(qū)危害性,以期為研究區(qū)公路、鐵路等重大工程的風險管理與防治設計提供參考。

1 流域概況

凍錯曲冰湖泥石流溝位于西藏昌都市洛隆縣臘久鄉(xiāng),擬建工程位于該泥石流溝溝口、凍錯曲主河下游3 km處(圖1)。凍錯曲冰湖泥石流溝流域面積226.27 km2,主溝長度31.44 km,流域范圍內(nèi)最高海拔約5 180 m,最低海拔約3 815 m,相對高差1 365 m,主溝平均縱坡降27.27‰。

圖1 凍錯曲流域全貌圖Fig.1 The view of the Dongcuoqu watershed

1.1 地質(zhì)環(huán)境條件

凍錯曲冰湖泥石流溝位于伯舒拉嶺山脈東北側(cè),為高山峽谷地貌,受區(qū)域環(huán)境因素和地勢影響,冰湖發(fā)育,該區(qū)域內(nèi)泥石流以冰川-暴雨混合型泥石流為主[19]。該流域地勢東高西低,溝道呈寬淺“V”形,支溝較為發(fā)育,谷坡25°～70°,谷寬20～100 m。溝口與巴曲、義俄4號溝交接,附近堆積扇由于人工改造較少,基本保持原地貌形態(tài)。研究區(qū)在區(qū)域構(gòu)造上屬隆格爾至納木錯斷裂帶,區(qū)內(nèi)延伸60 km,產(chǎn)狀203°～210°∠45°～50°,表現(xiàn)為逆推覆構(gòu)造,斜切凍錯曲冰湖上游匯水區(qū)。地表出露基巖主要為石炭-二疊系來姑組石英砂巖、白堊紀花崗巖、侏羅紀二長花崗巖(圖2)。流域內(nèi)支溝發(fā)育,主溝道縱坡降存在明顯的分段特性,可分為上游匯水區(qū)、形成流通區(qū)和溝口堆積區(qū)(圖3)。

圖2 凍錯曲冰湖泥石流溝物源分布圖Fig.2 Distribution of material sources in the Dongcuoqu glacier-lake debris flow watershed

圖3 凍錯曲冰湖泥石流溝主溝道縱斷面圖Fig.3 Longitudinal section of the main channel of the Dongcuoqu glacier-lake debris flow

1.2 物源條件

凍錯曲冰湖泥石流溝受常年積雪影響,流域內(nèi)松散物源以溝道物源、坡面物源、崩滑物源、凍融物源為主(圖2)。溝域內(nèi)可啟動物源共264處,總靜儲量約9.5 × 107m3,總動儲量約1.7 × 106m3。凍融物源及冰磧物穩(wěn)定性較好,只有在大型泥石流鏟刮下啟動,相比之下崩滑物源和溝道物源在暴雨條件下更易于啟動[24]。

1.3 水文條件

研究區(qū)所在區(qū)域?qū)贉貛О敫珊导撅L氣候區(qū),日照充足,干濕季節(jié)分明,多年平均降雨量490.23 mm,主要集中在6—9月,最大日降雨量39.2 mm。該流域年平均氣溫4.5～6.4 ℃,最高氣溫30.1 ℃,多出現(xiàn)在6—8月,多年平均無霜期128 d。凍錯曲流域?qū)倥?水源主要來自大氣降水、冰川融水和地下水,溝內(nèi)流水常年匯入凍錯曲,流量隨季節(jié)變化較大。

2 泥石流動力學特征參數(shù)

以擬建工程位置作為計算斷面,根據(jù)規(guī)范公式計算凍錯曲冰湖泥石流在不同降雨頻率下的動力學特征參數(shù)。

2.1 泥石流流速

根據(jù)《泥石流防治工程勘察規(guī)范(試行)》(T/CAGHP 006—2018),凍錯曲泥石流的易發(fā)性評分為83分,根據(jù)查表法得泥石流容重為1.572 t/m3,屬稀性泥石流,故采用式(1)計算泥石流流速。

(1)

式(1)中:VC為泥石流斷面平均流速,m/s;γH為泥石流固體物質(zhì)容重,t/m;φ為泥沙修正系數(shù),據(jù)查表法取0.54;1/n為清水河床糙率系數(shù),取1.2;R為水力半徑,一般可用平均泥深H代替,m;I為泥石流水力坡度,取全流域平均縱坡降26.12‰。

計算結(jié)果如表1所示。

表1 凍錯曲泥石流擬建工程處流速Table 1 Flow velocity of the Dongcuoqu debris flow at the location of the proposed construction site

2.2 泥石流峰值流量

假設泥石流與暴雨同步發(fā)生,且暴雨流量完全轉(zhuǎn)化為泥石流流量。研究表明,冰湖潰決型泥石流對比暴雨型泥石流具有流量放大效應,可將冰湖潰決考慮為溝道堵塞潰決的一部分,采用放大堵塞系數(shù)的雨洪法[25-26]。首先獲取不同降雨頻率下計算斷面的暴雨洪水設計流量,再通過采用合適的堵塞系數(shù)對泥石流流量進行快速評估。根據(jù)式(2)計算泥石流峰值流量。

QC=(1+φ)QDDC

(2)

式(2)中:QC為泥石流峰值流量,m3/s;QD為暴雨洪水設計流量,m3/s;DC為泥石流堵塞系數(shù),取5.5。

計算結(jié)果如表2所示。

表2 凍錯曲泥石流擬建工程處峰值流量Table 2 Peak flow of the Dongcuoqu debris flow at the location of the proposed construction site

2.3 一次泥石流過程總量

由于泥石流比一般洪水更具暴漲暴落的特點,一次泥石流過程一般均比較短,泥石流過程線可以概化成五角形,故一次泥石流過程總量按照式(3)進行計算。

Q=0.264TQC

(3)

式(3)中:Q為一次泥石流的總量,m3;T為泥石流歷時,s。

一次泥石流沖出的固體物質(zhì)總量按式(4)計算。

QH=Q(γC-γw)/(γH-γw)

(4)

式(4)中:QH為一次泥石流固體物質(zhì)沖出量,m3;γc為泥石流容重,t/m3;γw為清水容重,t/m3。

計算結(jié)果如表3所示。

表3 凍錯曲泥石流一次沖出總量Table 3 Total volume of run-out materials of the Dongcuoqu debris flow

2.4 泥石流沖擊力

泥石流沖擊力是泥石流防治工程設計的重要參數(shù),包括(流體)整體沖壓力和(石塊)單塊最大沖擊力。采用式(5)計算泥石流整體沖壓力。

(5)

式(5)中:δ為泥石流整體沖壓力,kPa;λ為受力物形狀系數(shù),取1.33;g為重力加速度,m/s2;α為受力面與泥石流沖壓力方向的夾角。

考慮石塊正面撞擊下游橋墩,單塊最大沖擊力按式(6)計算。

(6)

式(6)中:F為石塊沖擊力,kN;γ為動能折減系數(shù),取0.3;W為石塊質(zhì)量,t;α為石塊運動方向與受力面的夾角;C1、C2為石塊、橋墩的彈性變形系數(shù),取C1+C2=0.005。

計算結(jié)果如表4所示。

表4 凍錯曲泥石流整體沖擊力及單塊最大沖擊力Table 4 Impact pressure and maximum impact force of rock mass of the Dongcuoqu debris flow

3 凍錯曲冰湖發(fā)育特征

3.1 冰湖基本特征

凍錯曲冰湖距泥石流溝溝口約8.3 km,東側(cè)湖口地理坐標:30°18′18.26″N,96°4′0.41″E。湖面呈長條狀,軸向長度7 100 m,平均寬度 290 m。非汛期湖面高程為 4 055 m,汛期湖面高程為 4 060 m,常年處于溢流狀態(tài)。經(jīng)過對Landsat-8衛(wèi)星影像解譯,凍錯曲冰湖面積約2.05 km2。

冰湖庫容是反映冰湖蓄能大小與冰湖潛在危害性的重要指標,既能為冰湖危險程度初步評價提供參考,也是峰值流量及潰決洪水演進的估算基礎。目前尚不能通過遙感等手段直接測量冰湖庫容,一般采用經(jīng)驗公式間接估算。Huggel經(jīng)驗公式樣本取自與青藏高原地質(zhì)條件相近的阿爾卑斯山脈,且公式誤差波動較小[27-28],因此采用式(7)估算凍錯曲冰湖平均深度,進而計算冰湖庫容,計算結(jié)果見表5。

表5 凍錯曲冰湖基本參數(shù)Table 5 Basic parameters of the Dongcuoqu glacier lake

(7)

式(7)中:V為冰湖庫容,m3;A為冰湖面積,m2;D為冰湖平均深度,m。

3.2 冰湖堰塞體穩(wěn)定性分析

2018年11月及2019年8月兩次現(xiàn)場踏勘結(jié)果表明,凍錯曲冰湖堰塞體為兩岸大型崩滑堆積及右岸支溝泥石流堆積綜合形成,壩體主要以孤、塊石為主,其左岸大型崩塌堆積體結(jié)構(gòu)松散、植被稀少或無,而右側(cè)崩塌及局部泥石流堆積體上有低矮灌木分布。野外調(diào)查結(jié)果及遙感影像顯示,堰塞體高度約45 m,順溝長度約458 m,壩體頂寬約320 m,估算其總體積為6.59×106m3(圖4)。

圖4 凍錯曲冰湖堰塞體概況Fig.4 Overview of the Dongcuoqu glacier-lake barrier dam

凍錯曲冰湖常年處于穩(wěn)定滲流狀態(tài),在壩頂處無水漫流,但在壩體內(nèi)存在一個較為穩(wěn)定的滲流通道,溝內(nèi)可以觀察到在壩頂下游約400 m開始出現(xiàn)穩(wěn)定水流溝道,說明該壩體具有穩(wěn)定的滲透浸潤線,因此對凍錯曲冰湖的潰決風險評價尤其重要。冰湖潰決風險評價的主要內(nèi)容是對冰湖堰塞體的穩(wěn)定性評價[12],在缺乏壩體材料土石強度參數(shù)以及土力學參數(shù)的情況下,常采用堰塞壩穩(wěn)定性快速評價方法[29-30]。鐘啟明等[31]總結(jié)了國內(nèi)外學者基于形態(tài)學指標,提出的5種快速評價堰塞壩穩(wěn)定性的經(jīng)驗公式,并對世界范圍內(nèi)421 個堰塞壩實例進行了評價,結(jié)果表明DBI和邏輯回歸模型指標Ls(AHV)的準確性較好。DBI較Ls(AHV)錯判率更低,綜合考慮了包括漫頂和管涌等各種形式壩體潰決,且樣本具有廣泛性和代表性,因此采用DBI來判別凍錯曲冰湖堰塞體穩(wěn)定性。DBI表達式為

(8)

式(8)中:Ab為堰塞壩控制流域面積,km2;Hd為壩體高度,m;Vd為堰塞壩體積m3。

通常認為當DBI<2.75時,壩體穩(wěn)定;當2.753.08時,壩體不穩(wěn)定。

基于遙感影像劃定堰塞體控制流域面積,結(jié)合現(xiàn)場調(diào)查結(jié)果對其高度和體積進行估算。計算得到凍錯曲冰湖堰塞體的DBI為3.01,位于非穩(wěn)定區(qū)與穩(wěn)定區(qū)之間,從圖5中也可以發(fā)現(xiàn)有為數(shù)不少的不穩(wěn)定堰塞壩在緩沖區(qū)甚至穩(wěn)定區(qū),同樣對于穩(wěn)定堰塞壩也有此類情況出現(xiàn),因此不排除凍錯曲冰湖有潰決的可能性[14]。

圖5 凍錯曲冰湖堰塞體DBI判別圖[14]Fig.5 DBI discriminant diagram of the Dongcuoqu glacier-lake barrier dam[14]

3.3 潰口峰值流量計算

冰湖潰決潰口峰值流量的計算結(jié)果通常與堰塞壩的形態(tài)特征、壩前水深等因素有關,計算方法則主要根據(jù)潰決模式有所不同。對于凍錯曲冰湖,分別在瞬時部分潰決和瞬時全部潰決2種潰決模式下,對其潰口峰值流量進行計算。

劉建康等[32]統(tǒng)計早期西藏地區(qū)冰湖潰決記錄,認為洪水一般在下泄1/3庫容時出現(xiàn)峰值流量。因此亦采用1/3潰決模式對凍錯曲冰湖進行瞬時部分潰決計算,并且考慮極端工況下凍錯曲冰湖全部潰決模式作為參考。對于峰值流量的計算,國內(nèi)外不少學者對此進行研究,并提出了許多計算公式,但是這些公式中絕大部分都是基于冰湖庫容這一單一因子對峰值流量估算[33]。這里分別采用Froehlich針對潰口高度提出的式(9)計算凍錯曲冰湖瞬時部分潰決潰口峰值流量qp[34],以及謝任之[35]利用特征線法歸納出的式(10)計算瞬時全部潰決潰口峰值流量Qp,表達式為

qp=0.607V0.295h1.24

(9)

(10)

式中:V為冰湖庫容,m3;h為潰決口高度,m;λ為流量參數(shù),λ=m2m-0.5[2/(2m+1)]2m+1,m為溝道斷面指數(shù),取1.25;B為堰塞體寬度,m;g為重力加速度,m/s2;H1為潰壩前最大水深,取冰湖平均深度,m。

表6 不同潰決模式潰口峰值流量Table 6 Peak flow of different dam failure modes

(11)

(12)

式中:k為洪峰流量系數(shù),k=1+(γC-γw)/(γH-γc);γc為泥石流容重,t/m3;γw為清水容重,t/m3;γH為泥石流固體容重,t/m3。

4 冰湖潰決泥石流數(shù)值模擬

4.1 RAMMS

RAMMS采用改進的VS(Voellmy-Salm)單相流連續(xù)介質(zhì)模型,泥石流流體被假設為一種非穩(wěn)定以及非均質(zhì)流體,通過雙參數(shù)Voellmy摩擦模型來描述流動巖屑之間的摩擦行為[37]。研究表明,VS模型能夠準確模擬泥石流的運動,并已被廣泛運用于泥石流的動力學研究[38-39]。此外,RAMMS還能夠?qū)⒔Y(jié)果導出到GIS、修改地形數(shù)據(jù)(例如設置構(gòu)筑物)、增加附加參數(shù)(屈服應力等),并增強泥石流的預測效果[40]。

VS模型假設剪切變形集中在流動基面附近[37],將摩擦阻力分為兩部分:一是靜態(tài)摩擦阻力,用與法向應力成比例的干庫侖型摩擦因數(shù)μ表示,不依賴于速度;二是運動阻力,與速度平方阻力或黏性湍流摩擦因數(shù)ξ相關。這兩項都提供了流動阻力,可以綜合模擬高速區(qū)域和速度較小的尾部的流動行為。改進后的總摩擦阻力S表示為

(13)

式(13)中:N為法向應力,由ρhgcosφ表示;ρ為密度,kg/m3;g為重力加速度,kg/m2;φ為流體材料內(nèi)摩擦角,(°);h為流體高度,m;u為流體速度,m/s;N0為流體材料屈服應力。

當法向應力N較小(低流體高度h)時,剪應力S的數(shù)值將迅速從0變?yōu)镹0;當法向應力N較大時,屈服應力N0會使剪應力S增大,從而使泥石流提前停止。

4.2 模型設置

RAMMS為泥石流模擬提供物源釋放和水力釋放兩種方式。水力釋放基于流量和時間的關系,通過在給定位置處釋放一個流量對泥石流進行模擬,計算準確度較好,該方法可以顯著減少模擬時間(只需計算冰湖下游而不是全流域),且潰決洪水流量曲線恰好可以作為數(shù)值模型的輸入以進行下一步模擬,相比物源釋放更具優(yōu)勢。因此,基于凍錯曲冰湖堰塞體穩(wěn)定性評價結(jié)果,將不同潰決模式下冰湖潰口峰值流量輸入RAMMS,分別得到2種場景下的冰湖泥石流潰決演進過程:瞬時部分潰決和瞬時全部潰決。

考慮計算機的最佳計算能力,將ASTER GDEMV3獲取的30 m分辨率DEM處理后作為地形數(shù)據(jù),基于網(wǎng)格進行數(shù)值分析。根據(jù)RAMMS水文手冊,可以將泥石流流量-時間曲線概化為三點式水文過程圖,即通過定義泥石流峰值流量Qmax及其對應時刻t1、結(jié)束時間t2進行水力釋放,具體參數(shù)見表7。根據(jù)野外調(diào)查結(jié)果和物理原理計算進行參數(shù)率定,最終模型基本參數(shù)設置為:密度為1 572 kg/m3,重力加速度取9.8 m/s2,干庫侖型摩擦因數(shù)為0.18,黏性湍流摩擦因數(shù)為192 m/s2。

表7 RAMMS水力釋放參數(shù)Table 7 Parameters of hydrology in RAMMS

4.3 潰決演進過程分析

根據(jù)RAMMS模擬結(jié)果,凍錯曲冰湖泥石流在2種潰決模式下的演進過程相似,根據(jù)時序和流體運動特點,可歸納為4個階段:初始潰決階段、加速運動階段、減速運動階段、溝口停淤階段。

在初始潰決階段,冰湖堰塞體在極端地震或強降雨等外動力地質(zhì)作用下發(fā)生破壞,凍錯曲冰湖開始發(fā)生潰決,湖水傾瀉而下,此時瞬時部分潰決和瞬時全部潰決潰口處最大流深可分別達到34.62 m和50.96 m[圖6(a),圖7(a)]。在加速運動階段,泥石流鏟刮沿程溝道,夾帶溝道內(nèi)泥沙、碎石等松散固體物質(zhì)向下游運移,泥石流在下游溝道狹窄段流速相對更大;對比到達溝口的時刻,部分潰決相對于全部潰決具有一定的滯后性[圖6(b),圖7(b)]。凍錯曲冰湖下游主溝縱坡降僅有21.9‰,沖出溝口后,泥石流流速逐漸變小,進入減速運動階段,由于潰決洪水總量較大,存在一個持續(xù)沖出溝口的過程[圖6(c),圖7(c)]。 隨著潰決過程的結(jié)束,沖出的泥石流在溝口漫流后自然減速停淤,RAMMS會通過泥石流流體低通量平衡自動停止,凍錯曲冰湖瞬時部分潰決和瞬時全部潰決演進過程分別歷時5 h和3.5 h[圖6(d),圖7(d)]。

圖6 凍錯曲冰湖泥石流瞬時部分潰決演進過程Fig.6 Evolution process of instantaneous partial outburst of the Dongcuoqu glacier lake debris flow

圖7 凍錯曲冰湖瞬時全部潰決演進過程Fig.7 Evolution process of instantaneous complete outburst of the Dongcuoqu glacier lake debris flow

4.4 危害性評價

根據(jù)2種潰決模式下泥石流的影響范圍及特征參數(shù)可以評價其對下游工程的影響。如圖6(e)、圖7(e)所示,凍錯曲冰湖瞬時部分潰決和瞬時全部潰決堆積范圍均經(jīng)過工程位置,其中全部潰決的影響范圍更大;潛在威脅區(qū)泥石流平均深度分別達4.87 m和8.26 m,最大流速分別為5.29 m/s和5.74 m/s[圖6(f)、圖7(f)],將對擬建工程構(gòu)成一定影響。考慮到擬建工程僅從泥石流堆積區(qū)邊緣通過,可通過提高路基和橋梁標高以保留一定的安全高度,并采取排導槽或?qū)Я鞯痰确雷o措施對泥石流進行引流和攔截,避免其直接對擬建工程直接沖擊。由于冰湖潰決泥石流致災的復雜性與隨機性,在進行防治工程設計時,需特別注意對模擬結(jié)果的解讀和使用。

5 結(jié)論

(1)凍錯曲冰湖泥石流溝流域面積226.27 km2,主溝長度31.44 km,相對高差1 365 m,平均縱坡降27.27‰。流域內(nèi)松散物源豐富,溝內(nèi)流水常年匯入凍錯曲冰湖。考慮冰湖潰決流量放大效應,100年一遇凍錯曲泥石流在擬建工程處洪峰流量為2 509.10 m3/s,流速為4.99 m/s,整體沖擊力為53.18 kPa。

(2)凍錯曲冰湖面積2.05 km2,平均深度46.56 m,庫容9.54 × 107m3。無量綱堵塞指數(shù)(DBI)分析結(jié)果顯示凍錯曲冰湖堰塞體位于非穩(wěn)定區(qū)與穩(wěn)定區(qū)之間,存在發(fā)生潰決的可能性。運用RAMMS的Voellmy-Salm單相流模型對凍錯曲冰湖泥石流進行模擬,結(jié)果表明其潰決演進過程可歸納為初始潰決、加速運動、減速運動、溝口停淤4個階段。

(3)在2種潰決模式下,凍錯曲冰湖泥石流在潛在威脅區(qū)平均深度分別為4.87 m和8.26 m;在瞬時全部潰決場景下,擬建工程處最大流速為5.74 m/s,峰值流量為2 843.38 m3/s,與規(guī)范公式計算結(jié)果接近。通過提高路基和橋梁標高以保留一定的安全高度,并采取排導槽或?qū)Я鞯痰确雷o措施對泥石流進行引流和攔截,避免其直接對擬建工程直接沖擊。