河谷寬度對高速遠程滑坡碎屑化運動影響機制研究

【摘 要】基于物理模型實驗與數(shù)值模擬對不同河谷寬度條件下滑坡的運動破碎特征、速度特征與堆積結(jié)構(gòu)特征進行研究。該成果為致災(zāi)范圍研究提供了一種新思路，有利于防災(zāi)減災(zāi)的開展。

【關(guān)鍵詞】高速遠程滑坡；碎屑化運動；河谷寬度

【中圖分類號】P642.22【文獻標志碼】A

0 引言

高速遠程滑坡［1］是大規(guī)?；鶐r失穩(wěn)后快速或極快速的以碎屑流形式進行似流體、遠距離運動現(xiàn)象，主要集中在高山峽谷地帶，滑體方量多達數(shù)百萬方以上［2-6］，巨大的能量使得滑坡多堆積于河谷甚至爬升至對面山坡形成滑坡壩和堰塞湖，經(jīng)地震或大規(guī)模降雨極易失穩(wěn)造成壩體潰決、洪水泛濫，引發(fā)次生災(zāi)害。高速遠程滑坡的方量、溝道條件和溝道流量的變化形成不同規(guī)?；聣危?-8］，從而形成不同滑坡壩破壞模式與潰決模式［9-11］。Stefanelli等［12］在對300個滑坡分析基礎(chǔ)上提出兩個指數(shù)為滑坡壩形成提供可靠預(yù)測。年廷凱等［13］開展室內(nèi)物理實驗?zāi)M溝谷不同截面形狀來研究滑坡壩穩(wěn)定性。鄧明楓等［14］以茂縣宗渠堰塞壩為原型，通過漫頂潰決模型實驗研究壩體組成結(jié)構(gòu)對潰決速率和潰口形式的影響機制，分析壩體密度和級配特征對潰決過程的動力學(xué)作用。

雖然過往研究已經(jīng)取得豐碩成果，但在高山峽谷地區(qū)河谷寬度作為滑坡碎屑化運動過程的重要影響因素未得到系統(tǒng)研究與理論總結(jié)。本文以高速遠程滑坡為研究對象，基于滑坡物理模型實驗和離散元數(shù)值模擬開展不同河谷寬度對滑坡碎屑化運動的影響機制研究，為防災(zāi)減災(zāi)提供參考。

1 研究方法

1.1 物理模型實驗

1.1.1 實驗裝置

本文實驗裝置由斜槽、水平槽以及前緣阻擋結(jié)構(gòu)組成，如圖1所示。斜槽和水平槽固定連接置于鋼架上，斜槽坡度45°，長寬2.4 m×0.1 m，水平槽尺寸2.4 m×0.1 m。斜槽的最外端至0.4 m范圍為物源區(qū)，用于相似塊體的堆積釋放。前緣阻擋結(jié)構(gòu)由G型夾穩(wěn)固安裝在側(cè)板上，呈15°角狀結(jié)構(gòu)。監(jiān)測裝置由一臺高速相機、一臺高清攝像機和一組傳感器組件組成進行監(jiān)測與數(shù)據(jù)采集。

1.1.2 實驗材料

本文按一定比例（表1）配置相似材料，尺寸0.05 m×0.05 m×0.02 m。相似材料的力學(xué)參數(shù)由單軸抗壓強度試驗測得，得到單軸抗壓強度0.281 MPa，彈性模量0.035 7 GPa，泊松比0.11。

1.1.3 實驗方案

本文物理模型實驗部分主要考慮不同河谷寬度（即前緣阻擋結(jié)構(gòu)距轉(zhuǎn)折端的距離）對滑體運動堆積特征的影響，因此設(shè)置5個工況，具體參數(shù)設(shè)置見表2。

1.2 離散元數(shù)值模擬

本文離散元數(shù)值模型還原物理模型尺寸，設(shè)置斜槽和水平槽兩部分，長度均為2.4 m，槽寬0.1 m，槽深0.23 m，斜槽坡度45°，采用落雨法生成模型物料顆粒。設(shè)置節(jié)理層間顆粒膠結(jié)強度值為零來實現(xiàn)物理模型實驗中節(jié)理結(jié)構(gòu)巖體的堆疊形式，設(shè)置節(jié)理面之間顆粒數(shù)量是構(gòu)成節(jié)理顆粒數(shù)的5倍及以上來降低尺寸效應(yīng)。前緣阻擋結(jié)構(gòu)采用分段式建墻的方法以剛性體在指定位置生成位移邊界條件來實現(xiàn)。

2 結(jié)果

2.1 滑坡運動破碎特征分析

為充分了解塊體破碎演變以及速度變化，通過高速攝影圖像對滑坡在水平段運動停積過程逐幀分析以獲取運動學(xué)特征。

圖2為工況1的系列高速攝影圖像，以塊體到達轉(zhuǎn)折端前一幀為初始時刻t=0.0 s?；瑝K在轉(zhuǎn)折端處（t=0.0 s）呈現(xiàn)不同堆疊形式且存在分離，撞擊水平槽時發(fā)生破碎：完整塊體在t=0.142 35 s時從塊體中部破碎并向四周碎裂開來，上部塊體被迫向上彈起，塊體之間空隙增大，整體密實度減小。塊體撞擊破碎在前端產(chǎn)生的明顯飛霧是一種明顯的氣浪現(xiàn)象。當大部分塊體沖擊水平槽時（t=0.302 35 s）產(chǎn)生第二次氣浪現(xiàn)象。

在塊體撞擊水平槽以及前后部塊體由于速度差異產(chǎn)生的撞擊、塊體與槽面摩擦后，滑體已經(jīng)有了較大程度的破碎。在t=0.34 s時，可以觀察到碎屑顆粒沖擊至水平槽側(cè)壁，且在此后的0.30 s內(nèi)未掉落，說明塊體破碎賦予了碎屑顆粒極大的能量，使得破碎后的碎屑顆粒能向四周遠距離飛躍。對比模型實驗工況1，數(shù)值模擬也呈現(xiàn)出類似的運動破碎現(xiàn)象，如圖3所示。

圖4為工況1物理模型實驗和數(shù)值模擬的運動破碎過程對比圖。運動過程中的解體、氣浪、架空現(xiàn)象在物理模型實驗高速影像和數(shù)值模擬過程圖中都可以觀察到，反映了數(shù)值模型在描述塊體運動及破碎過程中具有很好的準確性。

2.2 滑坡速度特征分析

除了對高速相機影像直接進行分析外，還對系列影像進行處理，得到不同工況下運動速度的變化趨勢，如圖5所示。從圖中可以看出速度整體分布趨勢是上層塊體速度略大于底層塊體，斜槽段塊體速度略大于水平槽段塊體速度。從云圖和矢量圖中也可以明顯看出，塊體間存在速度差異時產(chǎn)生的霧狀飛濺顆粒也有不小的速度，約3.5 m/s左右，且存在一定的覆蓋范圍。隨著河谷寬度不斷增大，可以看出塊體爬上阻擋結(jié)構(gòu)的速度也在降低，并且產(chǎn)生的霧狀飛濺顆粒的速度也在降低。

圖6為工況1在不同時刻下數(shù)值模擬的滑塊運動速度圖。t=0.985 s時塊體前后部產(chǎn)生較為顯著的速度差異；塊體前部撞擊底板后迅速產(chǎn)生較大程度的破碎，前部塊體整體速度驟降，但角塊顆粒速度有所增加。t=0.106 s時前部塊體底部已經(jīng)因碰撞損失大量動能，速度驟降至1 m/s以下。t=1.1 s時后部塊體高速跟進與前部塊體發(fā)生劇烈碰撞，產(chǎn)生顯著的能量傳遞使前部塊體速度驟增至1.25 m/s，而后部塊體在短時間內(nèi)迅速停積，在t=1.19 s時幾乎就停止了運動，中后部分依次停積。從放大圖中可以明顯看出，運動過程中不同位置的塊體有較為明顯的速度差異，整體呈現(xiàn)前部、上部速度較快，后部、下部速度較慢的規(guī)律。

2.3 滑坡堆積結(jié)構(gòu)特征分析

滑坡堆積結(jié)構(gòu)特征是滑坡碎屑化運動的結(jié)果，也從側(cè)面反映了滑坡的運動特征。圖7為模型實驗和數(shù)值模擬的堆積體對比圖，三角形標注的為質(zhì)心運動距離，五角星標注為標準化最遠運動距離。從圖中可以看出，不同工況下物理模型實驗與數(shù)值模擬結(jié)果在質(zhì)心運動距離的變化趨勢方面保持高度一致性，標準化運動距離雖保持增加趨勢但存在一定誤差，這可能是因為塊體在物理模型實驗中運動至最遠運動距離后由于破碎程度增加而與槽底的摩擦作用減小后在重力作用下發(fā)生持續(xù)滑動從而產(chǎn)生差異。

圖8為滑坡壩堆積體側(cè)視對比圖。從圖中可以看出，各工況均排列規(guī)律，隨著河谷寬度增加滑坡堆積厚度逐漸降低，架空現(xiàn)象逐漸明顯，堆積體密實度也逐漸降低。從爬高的角度來看，模型實驗中河谷越寬，塊體的爬高能力越差，且在前緣阻擋結(jié)構(gòu)上的堆積范圍越小，這是因為塊體在爬上前緣阻擋結(jié)構(gòu)之前經(jīng)過的水平段距離越長，動能消耗越大，爬坡能力就越差，因此爬坡高度就越低。而在數(shù)值模擬中，滑坡運動距離較遠且呈離散形態(tài)堆積，這是因為塊體破碎程度較高時不會呈現(xiàn)出模型實驗中的粉末或粉塵，而是解體破碎成較小顆粒以滾動形式為主要運動方式，區(qū)別于模型實驗中以滑動為主的運動方式。

3 討論

本節(jié)分別對實驗數(shù)據(jù)和數(shù)值模擬數(shù)據(jù)進行歸一化處理，得到各工況運動學(xué)參數(shù)和相對破碎率的對比圖，如圖9所示。整體而言，各參數(shù)之間保持高度一致的規(guī)律性，并且標準化參數(shù)值之間差值較小。

對比分析實驗結(jié)果與數(shù)值模擬結(jié)果發(fā)現(xiàn)：河谷寬度越大，標準化最遠爬高距離越短、標準化質(zhì)心高度越低、標準化最遠運動距離越遠、標準化最遠爬高距離越短、標準化爬高展布面積越小、破碎率越小。原因在于：滑塊在河谷寬度較大的工況中，有更充分的時間在水平槽段運動，在這期間不斷發(fā)生摩擦耗能，因此會攜帶更少能量爬上前緣阻擋結(jié)構(gòu)，經(jīng)過短暫運動就完成停積，導(dǎo)致其標準化最遠爬高距離較小，塊體也就有更小的標準化爬高展布面積，即河谷寬度與標準化質(zhì)心運動距離、標準化最遠運動距離呈正相關(guān)，與其他運動學(xué)參數(shù)和破碎率呈負相關(guān)關(guān)系。

就相對破碎率而言，物理模型實驗和數(shù)值模擬表現(xiàn)出相同的規(guī)律，但是從數(shù)值來看，數(shù)值模擬中的相對破碎率值明顯高于實驗值，主要有兩方面原因：其一是由于顆粒數(shù)量的限制導(dǎo)致發(fā)生破碎后數(shù)值模擬對于破碎的表現(xiàn)相較于實驗略顯粗糙，破碎后的顆粒無法進一步破碎，大大提升了堆積體的相對破碎率；另一方面，離散元程序采用的是剛性顆粒集合體來進行建模和運算，以彈性行為為主，對于實驗中因劇烈碰撞而發(fā)生塑性變形的現(xiàn)象無法較好還原，也就導(dǎo)致塊體在碰撞過程中受到相對較大的沖擊力，黏結(jié)鍵處產(chǎn)生較大局部應(yīng)力，黏結(jié)鍵斷裂較多，堆積體相對破碎率增大。

綜上所述，離散元數(shù)值模擬在模擬物理模型實驗中具有較好的表現(xiàn)，但是受限于顆粒數(shù)目和顆粒單元的剛性特征，模型對于塊體受強沖擊力下的塑性應(yīng)變和相對破碎率這方面的再現(xiàn)存在一定局限性。

4 結(jié)論

（1）通過對物理模型實驗和數(shù)值模擬結(jié)果的分析討論，揭示了河谷寬度越大，滑坡的標準化質(zhì)心運動距離、標準化最遠運動距離呈上升趨勢，堆積體的標準化質(zhì)心高度、標準化爬高距離和標準化爬高展布面積呈下降趨勢的規(guī)律。

（2）對比數(shù)值模擬結(jié)果與物理模型實驗中的運動破碎特征、速度特征和堆積特征結(jié)果，發(fā)現(xiàn)數(shù)值模擬可以很好地再現(xiàn)物理模型實驗，并且對于氣浪現(xiàn)象、抑阻隆起現(xiàn)象、后部塊體反方向坐落停積和堆積體中的架空現(xiàn)象等都可以很好地還原。在運動學(xué)參數(shù)方面，數(shù)值模擬能夠很好地反映出實驗結(jié)果的變化趨勢，且差值較小，精度較高。

（3）分析了離散元數(shù)值模擬在研究高速遠程滑坡碎屑化運動中具備的適用性和存在的局限性。

參考文獻

［1］ Heim. Landslides and human lives［M］. Vancouver， B C： Bitech Publishers， 1932： 80-88.

［2］ 程謙恭， 張倬元， 黃潤秋. 高速遠程崩滑動力學(xué)的研究現(xiàn)狀及發(fā)展趨勢［J］. 山地學(xué)報， 2007， 25（1）： 72-84.

［3］ 陳達， 許強， 鄭光，等. 基于顆粒識別分析系統(tǒng)的碎屑流堆積物顆粒識別和統(tǒng)計方法研究［J］. 水文地質(zhì)工程地質(zhì)， 2021， 48（1）： 60-69.

［4］ 陳祿俊， 邢愛國， 陳龍珠，等. 高速遠程滑坡飛行數(shù)值分析［J］. 水文地質(zhì)工程地質(zhì)， 2008， 35（5）： 7-12.

［5］ 溫銘生， 陳紅旗， 張鳴之，等. 四川茂縣\"6·24\"特大滑坡特征與成因機制分析［J］. 中國地質(zhì)災(zāi)害與防治學(xué)報， 2017， 028（3）： 1-7.

［6］ 戴興建， 殷躍平， 邢愛國. 易貢滑坡-碎屑流-堰塞壩潰壩鏈生災(zāi)害全過程模擬與動態(tài)特征分析［J］. 中國地質(zhì)災(zāi)害與防治學(xué)報， 2019， 30（5）： 1-8.

［7］ 鐘啟明，單熠博. 堰塞壩穩(wěn)定性快速評價方法對比［J］. 人民長江， 2019， 50（4）： 20-24.

［8］ 柴賀軍，劉漢超，張倬元. 中國堵江滑坡發(fā)育分布特征［J］.山地學(xué)報， 2000， 18（B02）： 51-54.

［9］ 劉惠軍， 沈軍輝， 聶德新，等. 5·12汶川地震青川縣木魚鎮(zhèn)滑坡壩穩(wěn)定性分析［J］. 山地學(xué)報， 2008，26（3）： 263-266.

［10］ 黃來源， 季為. 四川什邡青牛沱崩塌群特征及其形成的堰塞壩體穩(wěn)定性評價［J］. 中國地質(zhì)災(zāi)害與防治學(xué)報， 2012， 23（1）： 39-47.

［11］ Sattar A， Konagai K， Ikeda T， et al. A quantitative approach to assess landform changes of Hattian Ballah landslide dam formed by the 2005 Kashmir Earthquake［J］. Bulletin of Earthquake Resistant Structure Research Center， 2009 （42）： 19-20.

［12］ Stefanelli C T， Segoni S， Casagli N， et al. Geomorphic indexing of landslide dams evolution［J］. Engineering Geology， 2016， 208： 1-10.

［13］ 年廷凱， 吳昊， 陳光齊. 堰塞壩穩(wěn)定性評價方法及災(zāi)害鏈效應(yīng)研究進展［J］. 水利水電快報， 2018， 39（8）： 1.

［14］ 鄧明楓， 陳寧生， 廖世春， 等. 堰塞壩漫頂潰決過程及其受組成結(jié)構(gòu)的影響［J］. 人民長江， 2012， 43（2）：30-35.

［基金項目］川藏鐵路重大基礎(chǔ)科學(xué)問題專項（項目編號：41941017）

［作者簡介］蘇婕（1998—），女，碩士，主要從事高速遠程滑坡研究工作。