易侵蝕層厚度對(duì)滑坡運(yùn)動(dòng)過(guò)程與堆積特征影響的模型試驗(yàn)研究

段釗 張慶 李明 朱錕 吳延斌

摘 要：為研究滑坡與易侵蝕層相互作用下產(chǎn)生的地質(zhì)構(gòu)造現(xiàn)象，利用砂盒試驗(yàn)裝置對(duì)滑坡的運(yùn)動(dòng)過(guò)程進(jìn)行研究，揭示滑體運(yùn)動(dòng)過(guò)程中速度場(chǎng)與位移場(chǎng)規(guī)律，剖析滑坡堆積體幾何形態(tài)及內(nèi)部構(gòu)造特征。結(jié)果表明：滑體運(yùn)動(dòng)過(guò)程分為3個(gè)階段，初始階段，滑體沿斜板下滑，速度和加速度快速增加；加速、穩(wěn)定階段，易侵蝕層會(huì)極大地抑制滑體的運(yùn)動(dòng)，滑體與易侵蝕層碰撞后速度與加速度瞬間降低，使得滑體位移量顯著減??；隨著易侵蝕層厚度的增加，滑體在撞擊易侵蝕層后速度減小的幅度增大，堆積體面積、長(zhǎng)度隨厚度的增加而逐漸減??；在易侵蝕層厚度增加的過(guò)程中，滑體對(duì)易侵蝕層的作用形式逐漸由擠向覆轉(zhuǎn)變，作用過(guò)程變得復(fù)雜；隨著厚度的增加，相互作用變得復(fù)雜，堆積體內(nèi)部發(fā)育直立褶皺、平臥褶皺、逆沖剪切、水平剪切、包卷等結(jié)構(gòu)，表面形成褶皺、隆起等構(gòu)造，反映滑體從壓縮到剪切運(yùn)動(dòng)的應(yīng)力狀態(tài)。研究結(jié)果為滑坡運(yùn)動(dòng)路徑中存在易侵蝕時(shí)的災(zāi)害范圍計(jì)算提供了一定的理論參考。

關(guān)鍵詞：滑坡；易侵蝕層；模型試驗(yàn)；運(yùn)動(dòng)過(guò)程；堆積形態(tài)；內(nèi)部構(gòu)造

中圖分類號(hào)：P 642.22

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

文章編號(hào)：1672-9315（2023）04-0733-10

DOI：10.13800/j.cnki.xakjdxxb.2023.0410

Model test study on the influence of the depth of erodible layer

on landslide movement process and deposits morphology

DUAN Zhao1，2，ZHANG Qing3，LI Ming4，ZHU Kun5，WU Yanbin6

（1.College of Geology and Environment，Xian University of Science and Technology，Xian 710054，China；

2.Institute of Ecological Enviromental Restoration in Mine Areas of West China，

Xian University of Science and Technology，Xian 710054，China；

3.Xian Building Materials Geological Engineering Survey Institute Co.，Ltd.，Xian 710003，China；

4.Xian Survey，Design and Research Institute Co.，Ltd.，of CREC，Xian 710054，China；

5.Xian Geological Environmental Monitoring Station，Xian 710018，China；

6.Ocean College，Zhejiang University，Hangzhou 310058，China）Abstract：To study the geological structure phenomenon under the interaction between landslide and erodible layer，the movement process of landslide is studied by using sand box test device，and the law of velocity field and displacement field in the movement process of landslide is revealed，the geometric shape and internal structure characteristics of landslide deposits are analyzed.The results show that the motion process of the sliding mass can be divided into three stages：initial stage，when the sliding mass slides along the inclined plate，the velocity and acceleration increase rapidly.In the acceleration and stabilization stages，the existence of the erodible layer will greatly inhibit the movement of the sliding mass，and the velocity and acceleration of the sliding mass immediately decrease after the collision with the erodible layer，resulting in the significant decrease of the sliding body position shift.With the increase of the depth of the erodible layer，the velocity decrease of the sliding mass after impacting the erodible layer increases.Meanwhile，the area and length of the deposits will gradually decrease with the increase of the depth.As the depth of the erodible layer increases，the action form of the sliding mass on the erodible layer gradually changes from squeezing to overlying，and the action process becomes complicated.With the increase of the depth，the interaction becomes complicated.The internal structure of the deposits develops upright fold，recesses fold，thrust shear，horizontal shear，enveloping and other structures，and the surface forms fold and uplift，which reflects the stress state of the sliding mass from compression to shear movement.The research results provide a theoretical reference for the calculation of disaster scope when

the erasion tends to emerge

in the landslide movement path.

Key words：landslide；erodible layer；model test；movement process；deposits morphology；internal structure

0 引 言

滑坡作為一種常見的地質(zhì)現(xiàn)象，運(yùn)動(dòng)時(shí)在復(fù)雜的地形和各種基底層（易侵蝕層）上快速擴(kuò)散，不可避免地對(duì)運(yùn)動(dòng)路徑上的易侵蝕層造成動(dòng)力侵蝕，經(jīng)歷鏟刮、推擠后其體積顯著增加，運(yùn)動(dòng)過(guò)程變得更為復(fù)雜［1-3］。易侵蝕層指存在于滑坡運(yùn)動(dòng)路徑上的堆積層，滑坡的運(yùn)動(dòng)過(guò)程和堆積特征會(huì)受到易侵蝕層的強(qiáng)烈影響［4-5］。

野外調(diào)查是研究滑體運(yùn)動(dòng)過(guò)程與堆積特征的主要方法之一。ZENG和WANG發(fā)現(xiàn)青藏高原尼續(xù)滑坡堆積體下伏的易侵蝕層中形成了底辟侵入體、卷曲層理、斷層和平臥不對(duì)稱褶皺等構(gòu)造，同時(shí)受擾動(dòng)易侵蝕層中顆粒具有較好的定向性，在地層內(nèi)部存在大量的破裂浮石，認(rèn)為該滑坡在運(yùn)動(dòng)過(guò)程中與下伏易侵蝕層存在強(qiáng)烈的相互作用［6-7］。DUFRESNE在新西蘭圓頂滑坡中發(fā)現(xiàn)滑坡下伏易侵蝕層存在底辟侵入體、旋卷層理和局部小斷層等構(gòu)造現(xiàn)象，并在維多利亞蘭普盧滑坡研究中發(fā)現(xiàn)了類似的現(xiàn)象，認(rèn)為這些滑坡在運(yùn)動(dòng)過(guò)程中對(duì)易侵蝕層造成了劇烈的擾動(dòng)［8-9］。DUAN對(duì)涇陽(yáng)南塬黃土滑坡研究發(fā)現(xiàn)，在滑坡堆積體內(nèi)部有從階地層中向上躍層流動(dòng)的流態(tài)砂質(zhì)粉土，在地表形成“砂沸”現(xiàn)象，在滑體前緣地表形成了橫向脊、徑向裂縫等特征，認(rèn)為滑體在運(yùn)動(dòng)過(guò)程中對(duì)階地飽和易侵蝕層造成了強(qiáng)烈沖擊和推擠［10］。

滑坡與沿途路徑上易侵蝕層的相互作用是普遍存在的，其作用下的形態(tài)特征與構(gòu)造跡象可作為推斷滑坡運(yùn)動(dòng)過(guò)程與受力特征的重要依據(jù)［11-12］?；峦ǔ１灰暈橐环N單次性事件，導(dǎo)致在野外調(diào)查中所獲取的數(shù)據(jù)并不能系統(tǒng)地反映易侵蝕層對(duì)滑坡運(yùn)動(dòng)特征的影響［13-14］?；碌男纬杀尘安町愋院艽螅瑑H通過(guò)地質(zhì)調(diào)查難以有效地分析該類滑坡的共性特征。而物理模型試驗(yàn)的引入對(duì)易侵蝕層與滑坡相互作用過(guò)程的研究提供了定量研究的方法和思路。MANGENEY進(jìn)行了斜板上鋪設(shè)有厚度可變易侵蝕層的滑坡模擬試驗(yàn)，表明滑體在與易侵蝕層碰撞后表現(xiàn)出高速擴(kuò)展和緩慢細(xì)流2個(gè)階段，在后一個(gè)階段侵蝕作用強(qiáng)烈，易侵蝕層的存在對(duì)滑坡的運(yùn)動(dòng)起到促進(jìn)作用，隨著厚度的增加，滑體運(yùn)移距離遠(yuǎn)且耗時(shí)更長(zhǎng)［15］。IVERSON進(jìn)行了在斜面上鋪設(shè)易侵蝕層的大型滑坡物理模型試驗(yàn)，發(fā)現(xiàn)當(dāng)易侵蝕層處于高飽和狀態(tài)下時(shí)，會(huì)產(chǎn)生較大超孔隙壓力，同時(shí)發(fā)生夾帶和動(dòng)量增加現(xiàn)象［16］。DUFRESNE在斜槽底部鋪設(shè)了厚度分別為15 mm和30 mm的PVC球和面粉制成的易侵蝕層［17］，易侵蝕層距坡腳有一定距離且表面與底板面齊平，試驗(yàn)結(jié)果表明在這種試驗(yàn)條件下易侵蝕層對(duì)滑坡運(yùn)動(dòng)能力有顯著的抑制作用，同時(shí)滑坡造成易侵蝕層發(fā)生水平剪切和擠壓，并將其夾帶至較遠(yuǎn)距離。SHEA將3 mm厚的易侵蝕層鋪設(shè)于光滑的鋁質(zhì)斜曲面上以研究滑體與易侵蝕層之間的相互作用［18］，試驗(yàn)結(jié)果顯示經(jīng)歷碰撞后的堆積體后緣廣泛發(fā)育正斷層，而前緣發(fā)育大量的逆沖斷層，認(rèn)為滑體與易侵蝕層相互作用時(shí)后緣表現(xiàn)出強(qiáng)烈的拉伸作用，而前緣主要以推擠作用為主。

在滑坡路徑上普遍存在易侵蝕層，這對(duì)于滑坡的運(yùn)移以及堆積形態(tài)有著顯著影響，滑體與易侵蝕層碰撞時(shí)存在強(qiáng)烈的相互作用，造成一系列復(fù)雜的堆積體形態(tài)與構(gòu)造變化。目前對(duì)于易侵蝕層所在的位置和材料類型往往是研究的重點(diǎn)，而對(duì)于易侵蝕層的厚度變化沒有被系統(tǒng)研究。通過(guò)砂盒物理模型試驗(yàn)在滑坡運(yùn)動(dòng)路徑上設(shè)置不同厚度的易侵蝕層，研究滑坡運(yùn)動(dòng)過(guò)程和堆積特征，以期為計(jì)算滑坡運(yùn)動(dòng)路徑中存在易侵蝕時(shí)的影響范圍提供一定的理論參考。

1 材料與方法

1.1 試驗(yàn)儀器和材料

采用砂盒試驗(yàn)裝置研究滑坡沖擊易侵蝕層的過(guò)程，裝置包含物理模型系統(tǒng)和監(jiān)測(cè)系統(tǒng)。物理模型系統(tǒng)包括斜板、底板和砂盒。斜板1.5 m×1.2 m（長(zhǎng)×寬），通過(guò)支架調(diào)節(jié)角度；底板1.5 m×1.2 m（長(zhǎng)×寬），水平放置，兩側(cè)有滑軌；砂盒用于裝滑體，為對(duì)開門結(jié)構(gòu)。監(jiān)測(cè)系統(tǒng)包括1臺(tái)3D掃描儀（8 fps）用于獲取滑體的運(yùn)動(dòng)過(guò)程以及堆積體的數(shù)字高程模型數(shù)據(jù)（DEM）；2臺(tái)高速相機(jī)（120 fps），分別用于記錄滑體運(yùn)動(dòng)過(guò)程中的正視和俯視影像（圖1）。

滑體與易侵蝕層均采用干燥的中細(xì)石英砂作為材料（圖2），該材料能很好的模擬出由于孔隙水壓力積累超過(guò)有效應(yīng)力而發(fā)生滑坡的地質(zhì)現(xiàn)象，在試驗(yàn)中能表現(xiàn)出良好的流態(tài)化特征［14，19］。砂的粒徑在0.075～0.5 mm范圍內(nèi)累積百分比為87.71%，粒徑分布曲線如圖3所示。通過(guò)直剪試驗(yàn)測(cè)得砂的內(nèi)摩擦角φ為34°，黏聚力為0，摩擦系數(shù)為0.675。在直剪試驗(yàn)中用有機(jī)玻璃替換下剪切盒中的砂樣，測(cè)得有機(jī)玻璃與砂的界面摩擦系數(shù)為0.474。試驗(yàn)所用石英砂為染色砂，目的是為更好地區(qū)分滑體與易侵蝕層，不同顏色石英砂的物理、力學(xué)性質(zhì)一致。

1.2 試驗(yàn)方法

尺度效應(yīng)的存在，使得野外實(shí)際滑坡的一些重要物理力學(xué)過(guò)程無(wú)法在模型試驗(yàn)中反映出來(lái)，包括動(dòng)力學(xué)、靜電現(xiàn)象［20］、巖石破碎和大型滑坡中的地震效應(yīng)［21］。野外滑坡復(fù)雜的邊界條件和多樣的材料性質(zhì)也難以在模型試驗(yàn)中再現(xiàn)［22］。這使得很難闡明單一因素如何影響滑坡的演變。在本研究中，簡(jiǎn)化了邊界條件和材料選擇，便于測(cè)試的重復(fù)性［23］。

在60°時(shí)，滑動(dòng)質(zhì)量可以從以往的研究中識(shí)別出運(yùn)動(dòng)學(xué)參數(shù)［10，24］。在進(jìn)行預(yù)試驗(yàn)時(shí)發(fā)現(xiàn)，當(dāng)易侵蝕層為24 mm時(shí)將不能完全侵蝕，易侵蝕層底部存在1～2 mm原始地層（圖14（c））。砂盒重心高度設(shè)置為1 m，使得滑體在坡面轉(zhuǎn)折點(diǎn)處可以擁有足夠的動(dòng)能，與易蝕層相互作用，產(chǎn)生更多可見的內(nèi)外沉積物結(jié)構(gòu)。本次選用的砂體試驗(yàn)材料具有與自然事件相似的低流變強(qiáng)度和沉積行為［25］，因此被作為試驗(yàn)材料。

試驗(yàn)設(shè)置斜板與水平板之間的角度α為60°，砂盒重心高度H為1 m，滑體體積為3.6×103cm3。在砂盒中裝砂時(shí)控制砂的質(zhì)量與體積，使每次試驗(yàn)時(shí)滑體的初始密度均達(dá)到1.5 g/cm3。設(shè)置易侵蝕層厚度分別為0，4，8，12，16，20 mm和24 mm，密度為1.5 g/cm3，其中0 mm為水平底板未鋪設(shè)易侵蝕層試驗(yàn)組。易侵蝕層由鋪砂器進(jìn)行鋪設(shè)。將彩砂裝入鋪砂器內(nèi)，隨著鋪砂器沿著底板兩側(cè)導(dǎo)軌向前移動(dòng)，抽出擋板使鋪砂器內(nèi)的彩砂向下流出，并通過(guò)鋪砂器后側(cè)的刮刀將彩砂表面整平。通過(guò)高度調(diào)節(jié)卡口調(diào)整鋪砂厚度，每4 mm更換彩砂顏色（圖4）。

打開砂盒底部的對(duì)開門，將滑體從砂盒釋放，利用監(jiān)測(cè)系統(tǒng)記錄滑坡運(yùn)動(dòng)過(guò)程和堆積特征。為了明確試驗(yàn)條件和滑坡堆積體特征，定義了滑體重心高度、斜板角度、滑體滑動(dòng)距離、易侵蝕層厚度、滑坡堆積體長(zhǎng)度、滑坡堆積體寬度、堆積厚度、堆積體面積以及堆積體周長(zhǎng)（圖5）。設(shè)滑坡運(yùn)移方向?yàn)閅軸，垂直方向?yàn)閄軸，滑坡位移為水平投影位移。在滑坡運(yùn)動(dòng)停止后，利用有機(jī)玻璃切板對(duì)滑坡堆積體進(jìn)行切割，觀察并分析內(nèi)部構(gòu)造特征。

2 結(jié)果與分析

2.1 滑體運(yùn)動(dòng)過(guò)程

根據(jù)滑體與易侵蝕的相互作用過(guò)程，將滑體的運(yùn)動(dòng)過(guò)程劃分為3個(gè)階段。初始階段滑體啟動(dòng)下滑，在與易侵蝕層碰撞之前，滑體在斜板上快速滑動(dòng)，在該階段，各試驗(yàn)組下滑體運(yùn)動(dòng)特征一致；加速階段滑體與易侵蝕層猛烈碰撞且互相混入，滑體滑動(dòng)受阻，長(zhǎng)度被顯著壓縮，滑速降低；穩(wěn)定階段滑坡夾帶易侵蝕層以低速滑移最終停止，過(guò)程中前緣滑體與易侵蝕層不會(huì)產(chǎn)生更多的相對(duì)位移。相比之下，未鋪設(shè)易侵蝕層的試驗(yàn)組，滑體在與水平底板碰撞后，滑體長(zhǎng)度壓縮幅度小，滑體在加速階段至穩(wěn)定階段滑動(dòng)速度的階段性變化不顯著（圖6）。

2.1.1 速度特征

初始階段（未與易侵蝕層發(fā)生碰撞前）滑體速度迅速增加，最快速度達(dá)到2.7×103～3×103 mm/s。在該階段，所有試驗(yàn)組下的速度場(chǎng)特征是一致的。當(dāng)滑體與易侵蝕層（或易侵蝕層厚度為0 mm）碰撞瞬間，滑體運(yùn)動(dòng)進(jìn)入加速階段，加速度迅速降低到負(fù)值，速度因慣性小幅增加后快速降低。易侵蝕層厚度越大，易侵蝕層對(duì)滑體運(yùn)動(dòng)的抑制作用越強(qiáng)，滑體加速度的下降幅度越大，碰撞后滑體減速越顯著。易侵蝕層厚度為24 mm時(shí)，易侵蝕層對(duì)滑體的阻力達(dá)到最大，在撞擊后的0.125 s內(nèi)速度下降80%。而易侵蝕層厚度為0 mm時(shí)，在撞擊后的0.125 s內(nèi)速度僅下降48%。在穩(wěn)定階段，滑體的加速度由負(fù)值趨于0，減速運(yùn)動(dòng)逐漸減緩?；w與易侵蝕層相互作用減弱，滑體速度逐漸降低，直至運(yùn)動(dòng)停止。當(dāng)易侵蝕層厚度為4，12，16 mm和20 mm時(shí)，滑體加速度出現(xiàn)短暫的正值，滑體做加速運(yùn)動(dòng)，表現(xiàn)為二次加速。主要是由于受前側(cè)滑體和被擾動(dòng)的易侵蝕層的阻礙，后側(cè)滑體在慣性力作用下在深度上出現(xiàn)了速度差，一部分滑體沿滑坡表層滑動(dòng)并快速翻越阻力較大的隆起區(qū)，在下坡過(guò)程中出現(xiàn)二次加速（圖7、圖8）。

2.1.2 位移特征

試驗(yàn)累計(jì)位移從砂盒底部開始計(jì)算（位移為0），以任意時(shí)刻滑體前緣位置與起點(diǎn)的差值作為該時(shí)刻滑體的累計(jì)位移（圖9）。初始階段，滑體的位移大致相同，但在后2個(gè)階段，隨著易侵蝕層厚度的增加，滑體位移顯著減小。當(dāng)易侵蝕層厚度為0時(shí)，滑體位移為1 533.9 mm，而易侵蝕層厚度為4 mm和24 mm時(shí)，位移分別下降至1 414.3 mm和1 312.7 mm，比易侵蝕層厚度為0時(shí)比例分別下降了7.7%和14.4%。

水平底板上易侵蝕層的存在對(duì)滑坡的運(yùn)動(dòng)能力有顯著的抑制作用。易侵蝕層在與滑坡相互作用下密度增加，厚度變大，滑體前緣的地形隆高，在滑體的推動(dòng)作用下易侵蝕層抗剪強(qiáng)度充分發(fā)揮，滑體運(yùn)動(dòng)阻力增大。易侵蝕層通過(guò)與滑體間顆粒的碰撞和摩擦來(lái)實(shí)現(xiàn)對(duì)速度和位移的抑制［26-29］。

2.2 堆積體平面形態(tài)

隨著易侵蝕層厚度的增加，滑坡堆積體范圍明顯縮小，幾何形態(tài)由橢圓形至舌形過(guò)渡，且X軸長(zhǎng)Y軸短（圖10）。相比于未鋪設(shè)易侵蝕層的試驗(yàn)組，存在易侵蝕層時(shí)堆積體前緣有明顯的拱形隆脊［30］。以隆脊為界可將堆積體劃分為前坡和后坡，所有試驗(yàn)組下前坡均比后坡的面積更大且坡角更陡，在未鋪設(shè)侵蝕層和侵蝕層較薄的試驗(yàn)組下，堆積體的后緣形成不連續(xù)洼地（圖11）。

提取堆積體數(shù)字高程模型（DEM），分析其形態(tài)參數(shù)（圖12）。在易侵蝕層厚度為0 mm時(shí)，堆積體長(zhǎng)、寬、厚以及面積分別為503.35，775.90，32.51 mm和2.9×105 mm2；在底板存在易侵蝕層時(shí)，堆積體長(zhǎng)度和面積都顯著下降，但厚度和寬度增加。隨著易侵蝕層厚度的增加，堆積體長(zhǎng)度和厚度的變化量相對(duì)于寬度更為顯著。當(dāng)易侵蝕層厚度為4～24 mm，堆積體長(zhǎng)度從368.94 mm下降至291.34 mm；寬度從626.0 mm增加至650.6 mm；厚度由56.42 mm增加至74.60 mm；面積由1.9×105 mm2下降至1.7×105 mm2；堆積體前坡角度由27.87°增加至36.32°；后坡角度由23°增加至29°（圖12（a）（b）（c）（d）和（e））。

隨著易侵蝕層厚度增加，堆積體周長(zhǎng)面積比呈現(xiàn)出先增大后減小的變化趨勢(shì)，在20 mm達(dá)到最大（圖12（d））。相對(duì)于鋪設(shè)易侵蝕層的試驗(yàn)組，未鋪設(shè)易侵蝕層的試驗(yàn)組下堆積體周長(zhǎng)面積比數(shù)值要小很多。堆積體周長(zhǎng)面積比數(shù)值越大，表明堆積體邊界越復(fù)雜，也說(shuō)明了滑體與易侵蝕層相互作用過(guò)程越強(qiáng)烈，作用形式越復(fù)雜。

易侵蝕層厚度越大，對(duì)滑坡運(yùn)動(dòng)的抑制作用越顯著，尤其表現(xiàn)在平面形態(tài)上。易侵蝕層厚度在20 mm及以下的條件下，均符合該規(guī)律。當(dāng)易侵蝕層厚度增加至24 mm時(shí)，堆積體長(zhǎng)和寬卻表現(xiàn)出增加趨勢(shì)。通過(guò)斷面可知在易侵蝕層厚度達(dá)到24 mm時(shí)，滑體未將易侵蝕層完全侵蝕，即侵蝕深度未至底板（圖14）?；w體積和易侵蝕材料一定，滑體的侵蝕深度存在最大限度，當(dāng)易侵蝕層厚度大于這個(gè)深度，滑體的沖擊力則不能將其擊穿。滑體在撞擊易侵蝕層后并不能沿光滑底板向前滑動(dòng)，侵蝕面上移，前緣易侵蝕層隆起程度低，在后部主滑體的推動(dòng)下前部表層滑體輕易翻過(guò)隆起區(qū)，形成運(yùn)動(dòng)距離增大的現(xiàn)象。

2.3 堆積體內(nèi)部構(gòu)造

通過(guò)分析滑坡堆積體的內(nèi)部構(gòu)造特征，推演滑體與易侵蝕層的相互作用形式是砂箱試驗(yàn)中的常見方法［31-32］。為記錄堆積體內(nèi)部構(gòu)造特征，試驗(yàn)利用透明切割板沿滑坡堆積體Y軸（中軸）進(jìn)行了切割（圖13）。除易侵蝕層厚度為24 mm外，易侵蝕層在深度上均被完全侵蝕。在滑體和易侵蝕層的相互作用下，堆積體內(nèi)部出現(xiàn)了直立褶皺、平臥褶皺、逆沖剪切、水平剪切、包卷等構(gòu)造現(xiàn)象（圖14）。

直立褶皺分布于堆積體前緣的易侵蝕層中，表明該區(qū)域受到滑體推力的影響，形成垂直于推力方向的擠壓變形。直立褶皺在所有試驗(yàn)組下普遍存在，隨著易侵蝕層厚度的增加，直立褶皺的范圍偏移至滑體下方，厚度有所增加。

平臥褶皺發(fā)育于堆積體前坡下伏的堆積體中。其成因?yàn)榛w在運(yùn)動(dòng)過(guò)程中將易侵蝕層完全鏟起并通過(guò)推動(dòng)和覆蓋使其發(fā)生折疊形成地層重復(fù)對(duì)稱。平臥褶皺在所有試驗(yàn)組下均有發(fā)育。當(dāng)易侵蝕層厚度較薄時(shí)（4 mm和8 mm），平臥褶皺的核部在水平方向上遠(yuǎn)于隆脊位置，隨著易侵蝕層厚度的增加，平臥褶皺核部和隆脊位置均向斜坡坡腳運(yùn)移。

逆沖剪切主要發(fā)育于平臥褶皺的核部和上側(cè)翼部，該區(qū)域所受到的剪應(yīng)力最集中。逆沖剪切在易侵蝕層厚度為4 mm和8 mm時(shí)不發(fā)育，隨著易侵蝕層厚度增加至20 mm，剪切面角度逐漸變陡。平臥褶皺的形成，使得滑體下伏的擾動(dòng)易侵蝕層厚度增加，導(dǎo)致滑坡水平方向運(yùn)動(dòng)阻力增大。受滑坡推力作用，擾動(dòng)的易侵蝕層會(huì)沿著褶皺的軸線發(fā)生相應(yīng)的變形和剪切位移。因剪切面傾向與滑體運(yùn)動(dòng)方向相反，表現(xiàn)為逆沖剪切構(gòu)造特征。

水平剪切主要分布于平臥褶皺的上側(cè)翼部，僅在易侵蝕層厚度為20 mm和24 mm的試驗(yàn)組中發(fā)育。滑體在運(yùn)動(dòng)的過(guò)程中，下部滑體（及被擾動(dòng)的易侵蝕層）受到的阻力大于上部滑體二者之間產(chǎn)生速度差，在二者接觸面附近形成水平剪切。

包卷主要發(fā)育于前坡下部，僅在易侵蝕層厚度為16 mm和20 mm試驗(yàn)組下發(fā)育。在滑體將易侵蝕層鏟起的同時(shí)，下部滑坡阻力增大，而上部滑體在慣性力作用下繼續(xù)向前運(yùn)動(dòng)并卷入平臥褶皺核部，最終被后側(cè)滑體覆蓋形成包卷構(gòu)造。

3 結(jié) 論

1）滑體運(yùn)動(dòng)過(guò)程可以分為3個(gè)階段：初始階段滑體沿斜板下滑，速度和加速度快速增加，該階段所有試驗(yàn)組下滑體速度場(chǎng)，位移場(chǎng)的特征都是類似的；在加速、穩(wěn)定階段，滑體與易侵蝕層碰撞后速度迅速降低，加速度瞬間降低至負(fù)值。

2）在滑體體積一定的條件下，易侵蝕層的存在，在一定厚度的范圍內(nèi)會(huì)對(duì)滑體的運(yùn)動(dòng)產(chǎn)生抑制作用。隨著易侵蝕層厚度的增加，滑坡堆積體范圍明顯縮小。隆脊將堆積體分為前坡和后坡，所有試驗(yàn)組下前坡均比后坡的面積更大且坡角更陡。

3）滑體體積一定，所具有的侵蝕能力是一定的。當(dāng)超過(guò)滑體的影響范圍時(shí)，深層次的易侵蝕層將難以受到擾動(dòng)和侵蝕。在受到擾動(dòng)和侵蝕的地層內(nèi)部會(huì)出現(xiàn)直立褶皺、平臥褶皺、逆沖剪切、水平剪切、包卷等現(xiàn)象。這表明在滑體內(nèi)部，隨著易侵蝕層厚度的增加，滑體對(duì)易侵蝕層的作用形式逐漸由擠向覆轉(zhuǎn)變。

參考文獻(xiàn)（References）：

［1］ LUCAS A，MANGENEY A，AMPUERO J P.Frictional velocity-weakening in landslides on Earth and on other planetary bodies［J］.Nature Communications，2014，5（1）：3417-3426.

［2］王玉峰，林棋文，李坤，等．高速遠(yuǎn)程滑坡動(dòng)力學(xué)研究進(jìn)展［J］.地球科學(xué)與環(huán)境學(xué)報(bào)，2021，43（1）：164-181.

WANG Yufeng，LIN Qiwen，LI Kun，et al.Review on rock avlanche dynamics［J］.Journal of Earth Sciences and Environment，2021，43（1）：164-181.

［3］殷躍平，王文沛.高位遠(yuǎn)程滑坡動(dòng)力侵蝕犁切計(jì)算模型研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2020，39（8）：1513-1521.

YIN Yueping，WANG Wenpei.A dynamic erosion plowing model of long run-out landslides initialized at high locations［J］.Chinese Journal of Rock Mechanics and Engineering，2020，39（8）：1513-1521.

［4］高楊，李濱，高浩源，等.高位遠(yuǎn)程滑坡沖擊鏟刮效應(yīng)研究進(jìn)展及問(wèn)題［J］.地質(zhì)力學(xué)學(xué)報(bào)，2020，26（4）：510-519.

GAO Yang，LI Bin，GAO Haoyuan，et al.Progress and issues in the research of impact and scraping effect of high-elevation and long-runout landslide［J］.Journal of Geomechanics，2020，26（4）：510-519.

［5］YIN Y，WANG W，ZHANG N，et al.Long runout geological disaster initiated by the ridge-top rockslide in a strong earthquake area：A case study of the Xinmo landslide in Maoxian County，Sichuan Province［J］.Geology in China，2017，44（5）：827-841.

［6］ZENG Q L，YUAN G X，MCSAVENEY M，et al.Timing and seismic origin of Nixu rock avalanche in southern Tibet and its implications on Nimu active fault［J］.Engineering Geology，2020，268（2）：105522-105536.

［7］WANG Y F，CHENG Q G，SHI A W，et al.Sedimentary deformation structures in the Nyixoi Chongco rock avalanche：implications on rock avalanche transport mechanisms［J］.Landslides，2018，16（3）：523-532.

［8］DUFRESNE A，WOLKEN G，HIBERT C，et al.The 2016 Lamplugh rock avalanche，Alaska：Deposit structures and emplacement dynamics［J］.Landslides，2019，16（12）：2301-2319.

［9］DUFRESNE A，DAVIES T R.Longitudinal ridges in mass movement deposits［J］.Geomorphology，2009，105（3-4）：171-181.

［10］DUAN Z，CHENG W C，PENG J B，et al.Interactions of landslide deposit with terrace sediments：Perspectives from velocity of deposit movement and apparent friction angle［J］.Engineering Geology，2021，280（5）：105913-105931.

［11］HUNGR O，LEROUEIL S，PICARELLI L.The Varnes classification of landslide types，an update［J］.Landslides，2013，11（2）：167-194.

［12］ZHANG M，TING Y P，WU S R，et al.Development status and prospects of studies on kinematics of long runout rock avalanches［J］.Engineering Geology，2010，18（6）：805-817.

［13］XU Y，ALLEN M B，ZHANG W，et al.Landslide cha-racteristics in the Loess Plateau，northern China［J］.Geomorphology，2020，359：107150-107188.

［14］DUAN Z，WU Y B，TANG H，et al.An analysis of factors affecting flowslide deposit morphology using Taguchi method［J］.Advances in Civil Engineering，2020（9）：1-14.

［15］MANGENEY A，ROCHE O，HUNGR O，et al.Erosion and mobility in granular collapse over sloping beds［J］.Journal of Geophysical Research，2010，115（F3）：3040-3061.

［16］IVERSON R M，REID M E，LOGAN M，et al.Positive feedback and momentum growth during debris-flow entrainment of wet bed sediment［J］.Nature Geoscience，2010，4（2）：116-121.

［17］DUFRESNE A.Granular flow experiments on the interaction with stationary runout path materials and comparison to rock avalanche events［J］.Earth Surface Processes and Landforms，2012，37（14）：1527-1541.

［18］SHEA T，VAN WYK DE VRIES B.Structural analysis and analogue modeling of the kinematics and dynamics of rockslide avalanches［J］.Geosphere，2008，4（4）：657-686.

［19］DUAN Z，CHENG W C，PENG J B，et al.Investigation into the triggering mechanism of loess landslides in the south Jingyang platform，Shaanxi province［J］.Bulletin of Engineering Geology and the Environment，2018，78（7）：4919-4930.

［20］IVERSON R M，LOGAN M，DENLINGER R P.Granular avalanches across irregular three-dimensional terrain：2.Experimental tests［J］.Journal of Geophysical Research：Earth Surface，2004，109（F1）：1015-1031.

［21］DAVIES T R H，MCSAVENEY M J.Runout of dry granular avalanches［J］.Canadian Geotechnical Journal，1999，36（2）：313-320.

［22］DELANNAY R，VALANCE A，MANGENEY A，et al.Granular and particle-laden flows：from laboratory experiments to field observations［J］.Journal of Physics D：Applied Physics，2015，50（5）：53001-53057.

［23］CROSTA G，IMPOSIMATO S，RODDEMAN D.Granular flows on erodible and non erodible inclines［J］.Granular Matter，2015，17（5）：667-685.

［24］LI K，WANG Y F，LIN Q W，et al.Experiments on granular flow behavior and deposit characteristics：implications for rock avalanche kinematics［J］.Landslides，2021，18（5）：1779-1799.

［25］MANZELLA I，LABIOUSE V.Flow experiments with gravel and blocks at small scale to investigate parameters andmechanisms involved in rock avalanches［J］.Engineering Geology，2009，109（1-2）：146-158.

［26］LIU W，WANG D，ZHOU J，et al.Simulating the Xinmo landslide runout considering entrainment effect［J］.Environmental Earth Sciences，2019，78（19）：585-601.

［27］ZHOU Y，SHI Z，ZHANG Q，et al.3D DEM investigation on the morphology and structure of landslide dams formed by dry granular flows［J］.Engineering Geology，2019，258（3）：105151-105162.

［28］CROSTA G B，DE BLASIO F V，DE CARO M，et al.Modes of propagation and deposition of granular flows onto an erodible substrate：experimental，analytical，and numerical study［J］.Landslides，2016，14（1）：47-68.

［29］YUAN R M，TANG C L，HU J C，et al.Mechanism of the Donghekou landslide triggered by the 2008 Wenchuan earthquake revealed by discrete element modeling［J］.Natural Hazards and Earth System Sciences，2014，14（5）：1195-1205.

［30］范奉鑫，林美華，江榮華，等.東海陸架前緣斜坡（沖繩海槽西坡）北部的斷塊隆脊地貌［J］.青島海洋大學(xué)學(xué)報(bào)（自然科學(xué)版），2000，30（1）：173-176.

FAN Fengxin，LIN Meihua，JIANG Ronghua，et al.The rift block ridge in north continental shelf front slope of the East China Sea［J］.Journal of Ocean University of Qingdao，2000，30（1）：173-176.

［31］PAOLA C，STRAUB K，MOHRIG D，et al.The “unreasonable effectiveness” of stratigraphic and geomorphic experiments［J］.Earth Science Reviews，2009，97（1-4）：1-43.

［32］DENG B，LIU S，LI Z.Structural geometry and kinematic processes at the intracontinental Daloushan mountain chain［J］.Comptes Rendus Geosciences，2015，348（2）：159-168.

（責(zé)任編輯：李克永）