何旭乾，盧才武，李 萌，陳彥鈺，閆雪頌

(西安建筑科技大學(xué)資源工程學(xué)院，陜西 西安 710055)

露天礦山在開采或自然條件等因素影響下，一些邊坡會發(fā)生失穩(wěn)，從而引發(fā)滑坡，對礦山開采和人員生命安全造成嚴(yán)重威脅。因此，對滑坡風(fēng)險(xiǎn)定量分析進(jìn)行研究，對于滑坡災(zāi)害的預(yù)防治理具有重要意義。隨著計(jì)算機(jī)性能的不斷進(jìn)步，各種數(shù)值模擬軟件由于自身低成本、兼容強(qiáng)等優(yōu)點(diǎn)，在滑坡災(zāi)害模擬與致因分析過程中應(yīng)用廣泛。顆粒流(particle flow code，PFC)程序是基于CUNDALL等[1]提出的離散元理論，通過Itasca開發(fā)的二維模擬程序或三維模擬程序。TANG等[2-3]應(yīng)用PCD2D分別對1941年和1999年的滑坡的運(yùn)動(dòng)特性、裂紋擴(kuò)展、斷裂特性進(jìn)行分析；LO等[4]使用PCD3D對導(dǎo)致滑坡和破壞的運(yùn)動(dòng)學(xué)過程進(jìn)行模擬，評估了滑坡機(jī)理、泥石流運(yùn)動(dòng)、沉積程度。由此可知，PFC軟件能夠模擬剛性粒子的運(yùn)動(dòng)和相互作用，允許其旋轉(zhuǎn)和移動(dòng)，適用于滑坡形變破壞、滑移運(yùn)動(dòng)、堆積特征的分析。為了反映真實(shí)的滑坡運(yùn)動(dòng)，曹文等[5]指出了wall-ball法具備顆粒少、節(jié)省運(yùn)算時(shí)間、結(jié)果準(zhǔn)確度高的優(yōu)點(diǎn)，但是由于PFC軟件在滑坡建模這一前期處理方面不具備優(yōu)勢，因而受到一定的局限。

當(dāng)今無人機(jī)技術(shù)日趨成熟，其具備的高精度、高效率、低成本、低風(fēng)險(xiǎn)等優(yōu)點(diǎn)能夠很好地解決PFC軟件的不足。MA等[6]使用無人機(jī)劃分滑坡源區(qū)、重疊區(qū)和堆積區(qū)，結(jié)合DEM估算滑坡體積并分析DEM運(yùn)動(dòng)過程，定量描述了滑坡特征。因此，通過無人機(jī)建立高精度數(shù)字高程模型(digital elevation model，DEM)，能為PFC軟件的滑坡運(yùn)動(dòng)模擬提供準(zhǔn)確的三維離散元模型。

綜上所述，本文對于滑坡運(yùn)動(dòng)過程離散元數(shù)值模擬主要流程為：首先通過無人機(jī)航攝影像建立滑坡區(qū)DEM并復(fù)現(xiàn)原始地形，構(gòu)建出滑坡地形(wall)；然后以三軸壓縮試驗(yàn)與BP神經(jīng)網(wǎng)絡(luò)實(shí)現(xiàn)細(xì)觀參數(shù)的選取，對顆粒(ball)進(jìn)行參數(shù)賦值；最后構(gòu)建完整的ball-wall離散元三維滑坡模型，并從整體和局部的角度對滑坡區(qū)域的運(yùn)動(dòng)過程、堆積特征進(jìn)行反演與分析，為滑坡的定量評估提供一套快捷、有效、完整的解決思路。

1 研究區(qū)概述

1.1 露天礦區(qū)基本工程概況

本文選取位于陜西省境內(nèi)某露天煤礦作為研究對象(圖1)，該露天煤礦地處黃土高原，大部分地區(qū)為風(fēng)沙灘地地貌，總體地勢平坦，海拔為1 060～1 230 m。當(dāng)前采場非工作幫整體邊坡角為24°～32°，南端幫為15°，北端幫為17°。坡體巖性由上至下分別為松散沙層、粉質(zhì)黏土、砂巖，粉質(zhì)黏土顏色呈淺紅色或棕紅色，各種中砂、細(xì)砂、粉砂為淺黃色或褐黃色。雖然整體邊坡處于穩(wěn)定狀態(tài)，但是根據(jù)現(xiàn)場勘查和對地質(zhì)資料分析，部分位置在風(fēng)化剝蝕作用下仍出現(xiàn)了局部滑坡現(xiàn)象。

圖1 露天礦現(xiàn)場情況與巖性分布Fig.1 Open-pit mine site situation and lithology distribution

1.2 滑坡現(xiàn)場地質(zhì)情況

根據(jù)對采場現(xiàn)場的調(diào)查與分析，造成滑坡的主要影響因素包括節(jié)理裂隙、邊坡角設(shè)計(jì)、邊坡暴露時(shí)長等。該區(qū)域內(nèi)黃土層和紅土層垂直節(jié)理發(fā)育，自然風(fēng)化會引起局部淺層滑坡，邊坡巖體強(qiáng)度隨著時(shí)間增長呈下降趨勢，易產(chǎn)生滑坡。

該滑坡位于采場西北方向非工作幫的地表臺階上，臺階高度為20 m，整體地形主要由粉質(zhì)黏土構(gòu)成，采場區(qū)內(nèi)土層垂直節(jié)理發(fā)育，經(jīng)過長期暴露，邊坡巖體隨時(shí)間增長出現(xiàn)應(yīng)力松弛、強(qiáng)度下降的現(xiàn)象，部分坡體向下滑落產(chǎn)生滑坡并形成堆積，根據(jù)滑坡后堆積DEM可知(圖2)，主要滑動(dòng)幾乎是朝向東南方，滑移區(qū)域高約19 m，平均寬度約17 m，滑坡由頂端開始滑落，最終堆積層位于底端。

圖2 滑坡區(qū)域Fig.2 Landslide area

2 數(shù)值模型建立及參數(shù)確定

2.1 露天礦邊坡三維數(shù)值模型構(gòu)建

本文使用大疆精靈4 Pro無人機(jī)對露天礦區(qū)進(jìn)行航測，基于無人機(jī)影像數(shù)據(jù)，使用Agisoft Photoscan生成分辨率礦區(qū)的DEM數(shù)據(jù)。通過Global Mapper對進(jìn)行影像裁剪處理，選出不穩(wěn)定邊坡區(qū)域，并導(dǎo)入Rhinoceros軟件中，構(gòu)建STL網(wǎng)格滑坡模型。在PFC3D顆粒流程序中將模型作為Geometry進(jìn)行導(dǎo)入，Wall邊界設(shè)置為模型邊界。

在對滑坡區(qū)域進(jìn)行反演的過程中需要對原始地形進(jìn)行復(fù)現(xiàn)。根據(jù)DUMAN[7]的研究表明，通過對滑坡體積的估算來建立原始地形的方法具有可行性。在Rhinoceros軟件中導(dǎo)入滑坡模型后，生成滑動(dòng)后地形等高線。通過參考ALOS PALSAR 12.5 m的DEM數(shù)據(jù)來生成滑坡前地形等高線，在Rhinoceros軟件根據(jù)等高線劃分滑移區(qū)域。由圖3可知，滑坡區(qū)域地形發(fā)生了明顯改變，滑坡區(qū)域呈現(xiàn)下窄上寬，上部寬約23.3 m，下部寬約7.0 m，相對高差19.3 m，滑坡面積約360 m2，體積約為207 m3，1 150 m高程以下的等高線發(fā)生明顯的彎曲，根據(jù)實(shí)際調(diào)查情況可以看出其下方形成了堆積。

圖3 滑坡區(qū)域滑移前后的高程線Fig.3 Elevation lines before and after slippage in the landslide area

由于滑移區(qū)和堆積區(qū)存在重疊，因此在完成地形重建后，可以看出滑坡的頂端較陡，滑移區(qū)呈現(xiàn)倒梯形，邊坡滑動(dòng)面較淺，滑移區(qū)與堆積區(qū)有部分重疊，滑坡頂端與堆積體高差為14 m(圖4)。

圖4 滑坡區(qū)域三維地形Fig.4 3D topography of landslide area

2.2 巖土體的細(xì)觀力學(xué)參數(shù)確定

在PFC3D顆粒流軟件建立滑坡模型過程中，需要對模型賦予相應(yīng)的參數(shù)，由于材料的細(xì)觀參數(shù)與實(shí)際巖土體宏觀力學(xué)參數(shù)沒有直接的量化關(guān)系式，為了使得模擬結(jié)果接近實(shí)際情況，需要通過三軸壓縮試驗(yàn)，確定合適的細(xì)觀參數(shù)來準(zhǔn)確反映實(shí)際的巖石力學(xué)特征。一般來說，這些細(xì)觀參數(shù)包括顆粒彈性模量Ec、顆粒法向剛度kn、顆粒切向剛度ks、摩擦系數(shù)μ、法向黏結(jié)強(qiáng)度σ、切向黏結(jié)強(qiáng)度c，這些參數(shù)對彈性模量E、泊松比?、黏聚力等宏觀力學(xué)參數(shù)具有較大影響。

三軸壓縮試驗(yàn)的目的是對巖土體的宏細(xì)觀參數(shù)進(jìn)行標(biāo)定，以此獲得實(shí)驗(yàn)數(shù)據(jù)。在數(shù)值模擬實(shí)驗(yàn)中，首先根據(jù)以往對巖土體宏細(xì)觀參數(shù)的研究經(jīng)驗(yàn)，設(shè)置初始的細(xì)觀參數(shù)，然后在PFC3D軟件內(nèi)建立三軸試驗(yàn)?zāi)Ｐ?，設(shè)置符合滑坡模型的顆粒大小和接觸模型，并對初始細(xì)觀參數(shù)進(jìn)行模擬。通過對應(yīng)變曲線的分析，可以獲取到對應(yīng)數(shù)值的宏觀力學(xué)參數(shù)，將與實(shí)際試驗(yàn)結(jié)果相近的參數(shù)值作為輸出樣本。

該滑坡體主要為粉質(zhì)黏土，破壞誘發(fā)因素為自然風(fēng)化與長期暴露。首先，按照巖樣尺寸生成邊界墻體，依照墻體生成緊密顆粒，指定最大顆粒與最小顆粒的半徑比，完成設(shè)定后PFC3D系統(tǒng)會按照預(yù)設(shè)值，將最大值和最小值隨機(jī)均勻附加在模型的每個(gè)顆粒上，最終壓縮巖樣如圖5所示。

圖5 三軸試驗(yàn)數(shù)值模擬模型Fig.5 Numerical simulation model of triaxial test

為保證生成模型的顆粒尺寸和分布滿足實(shí)驗(yàn)要求，使得顆粒之間具備較高的接觸精度，需要對模型進(jìn)行均勻化處理。均勻化處理的目的是要保證生成的各相鄰顆粒之間的重疊度足夠小，讓顆粒與邊界接觸緊密，達(dá)到完全耦合的效果。根據(jù)圖6的均勻化處理，利用自編程序不斷調(diào)整模型邊界，使得顆粒之間接觸處于足夠理想的狀態(tài)。

圖6 三軸試驗(yàn)均勻化處理Fig.6 Homogenization of triaxial test

接觸模型選擇為平行黏結(jié)模型，該模型可以同時(shí)傳遞力和力矩，受到破壞后剛度下降，與滑坡破壞特性相近。完成接觸模型設(shè)定后，將一組細(xì)觀參數(shù)值進(jìn)行賦值，進(jìn)行軸向加載，完成數(shù)據(jù)測定后繪制應(yīng)力-應(yīng)變曲線，根據(jù)測定數(shù)據(jù)計(jì)算宏觀力學(xué)參數(shù)。

為了確定細(xì)觀力學(xué)參數(shù)對巖土體宏觀力學(xué)參數(shù)與力學(xué)特性的影響，需要對細(xì)觀力學(xué)參數(shù)做敏感性分析，以此定性分析宏細(xì)觀參數(shù)間的關(guān)系。根據(jù)以往學(xué)者[9-11]的研究，主要包括黏結(jié)強(qiáng)度、摩擦系數(shù)、剛度比、顆粒形狀的影響。黏結(jié)強(qiáng)度決定了細(xì)觀層面上顆粒間膠結(jié)強(qiáng)弱度，從而影響宏觀的峰值強(qiáng)度和應(yīng)力大小，若增大法向平行黏結(jié)強(qiáng)度與切向平行黏結(jié)強(qiáng)度，會引起彈性模量和峰值強(qiáng)度的增大。當(dāng)摩擦系數(shù)增大時(shí)，彈性模量和峰值強(qiáng)度也在增長但漲幅較小。而顆粒的法向剛度與切向剛度之比將對泊松比產(chǎn)生影響，隨著剛度比的增大，泊松比也會增大。模擬實(shí)驗(yàn)中顆粒形狀采用球形，不規(guī)則顆粒形狀雖然能夠更好地模擬現(xiàn)實(shí)中巖石的幾何形態(tài)，但考慮到該類塊體形態(tài)難以考慮到土石間凹凸不平引起的咬合作用，使得材料模擬結(jié)構(gòu)與實(shí)際差異較大引起實(shí)驗(yàn)結(jié)論不準(zhǔn)確，因此后續(xù)模擬將不考慮這一因素。

一般獲取細(xì)觀力學(xué)參數(shù)的方法需要假定細(xì)觀力學(xué)參數(shù)，通過PFC3D軟件內(nèi)三軸實(shí)驗(yàn)獲取宏觀力學(xué)參數(shù)并與實(shí)際室內(nèi)試驗(yàn)數(shù)據(jù)進(jìn)行對比來獲取符合參數(shù)。該方法較為耗時(shí)，因此本文采用機(jī)器學(xué)習(xí)方法中的BP神經(jīng)網(wǎng)絡(luò)，對上述三軸試驗(yàn)的數(shù)值模擬過程獲取數(shù)據(jù)進(jìn)行訓(xùn)練。具體流程是讓隨機(jī)組合的細(xì)觀參數(shù)，以及模擬相對應(yīng)宏觀力學(xué)參數(shù)作為輸出樣本和輸入樣本，通過機(jī)器學(xué)習(xí)模型進(jìn)行訓(xùn)練，根據(jù)訓(xùn)練樣本構(gòu)建出宏觀參數(shù)與細(xì)觀參數(shù)之間的非線性映射關(guān)系，對細(xì)觀力學(xué)參數(shù)進(jìn)行反演。將10組宏觀力學(xué)參數(shù)選做測試樣本來驗(yàn)證模型，最終將實(shí)際的宏觀力學(xué)參數(shù)帶入模型中對細(xì)觀力學(xué)參數(shù)進(jìn)行反演(表1)，分別獲取到最合適細(xì)觀標(biāo)定參數(shù)。

應(yīng)用反演精度公式[11]、平均殘差公式[12]、殘差均方差公式[13]來驗(yàn)證BP神經(jīng)網(wǎng)絡(luò)模型的性能，經(jīng)過驗(yàn)證其精度滿足建模要求。將表1獲得的細(xì)觀力學(xué)參數(shù)導(dǎo)入PFC3D軟件的三軸壓縮試驗(yàn)中，獲得應(yīng)力-應(yīng)變曲線(圖7)。

表1 模型細(xì)觀力學(xué)參數(shù)Table 1 Mesomechanical parameters of model

由圖7可知，依照摩爾-庫倫強(qiáng)度理論計(jì)算滑坡黏聚力和內(nèi)摩擦角[14-17]，計(jì)算同室內(nèi)三軸壓縮實(shí)驗(yàn)數(shù)據(jù)對比，由表2所得的對比值可知，其宏觀力學(xué)參數(shù)值非常相近，證明該方法可行，可以用于構(gòu)建露天礦滑坡的顆粒流模型。

圖7 應(yīng)力-應(yīng)變曲線示意圖Fig.7 Schematic diagram of stress-strain curve

表2 模擬值的宏觀力學(xué)參數(shù)對比Table 2 Comparison of macro-mechanical parameters of simulated values

3 邊坡滑移過程模擬與分析

3.1 滑移運(yùn)動(dòng)過程模擬與分析

建立的滑坡三維模型(圖8)，滑移區(qū)域共生成11 002顆顆粒，顆粒高程差約為19 m。在對顆粒附加細(xì)觀參數(shù)后，施加重力即可進(jìn)行滑坡運(yùn)動(dòng)特征的模擬?；履M過程中分別從整體區(qū)域滑移和區(qū)域特征點(diǎn)兩個(gè)方面分析，在滑坡區(qū)域分別選取9個(gè)特征點(diǎn)作為監(jiān)測點(diǎn)(圖9)，在模擬中對其速度和位移進(jìn)行監(jiān)測，監(jiān)測點(diǎn)的坐標(biāo)分別為A(-8，13，10)、B(-3，14，10)、C(3，14，10)、D(-7.5，10，7)、E(-2.5，11，7)、F(2.5，11，7)、G(-5，5，3.5)、H(-1.5，5.5，3.5)、I(2，5.5，3.5)。

圖8 原始區(qū)域PFC滑坡模型Fig.8 Original regional PFC landslide model

圖9 監(jiān)測點(diǎn)位置與編號Fig.9 Location and number of monitoring points

在將滑坡體的速度和位移清零后，完成顆粒與顆粒、顆粒與墻體間的細(xì)觀參數(shù)賦值，并刪除滑移坡面墻體，最后對滑坡模型施加重力，讓顆粒在重力作用下進(jìn)行滑移，其整體運(yùn)動(dòng)過程如圖10所示。由圖10可知，滑坡從失穩(wěn)到運(yùn)動(dòng)停止一共經(jīng)歷20 s。在邊坡滑移前期中，1 s時(shí)坡面發(fā)生明顯形變，坡面頂端首先出現(xiàn)形變，坡表面顆粒速度達(dá)到了約10 m/s；到2 s時(shí)坡體迅速下滑并開始形成堆積，整體坡面繼續(xù)進(jìn)行滑移，坡底端顆粒速度達(dá)到15 m/s，滑坡整體運(yùn)動(dòng)速度達(dá)到峰值；在3 s左右隨著堆積形成，滑移速度逐漸減慢，底端顆粒在堆積后繼續(xù)移動(dòng)，底端局部速度為17.5 m/s；到達(dá)5 s時(shí)，頂端沉降速度開始減慢，隨著滑移、摩擦對能量損耗，滑坡運(yùn)動(dòng)速度顯著減慢，底端局部顆粒仍舊向前移動(dòng)。隨著大部分顆粒的堆積，在10 s時(shí)滑坡體后端逐漸趨于穩(wěn)定，前中端局部顆粒保持2.0 m/s速度繼續(xù)移動(dòng)；在15 s時(shí)滑坡體前端整體開始穩(wěn)定，直至20 s時(shí)大部分顆?；眷o止，最終形成完整的堆積體，只有小部分土體以較小速度繼續(xù)運(yùn)動(dòng)。

圖10 邊坡滑移過程Fig.10 Slope slip process

應(yīng)用PFC程序自編FISH函數(shù)進(jìn)行滑坡整體的位移與速度變化監(jiān)測(圖11)。由圖11可知，滑坡平均速度在快速提升后又急速減緩，逐漸趨于平緩，在1.8 s時(shí)達(dá)到峰值約為5.6 m/s；平均位移不斷增加，但是增長速度逐漸減緩，最終位移距離約為13 m。

圖11 邊坡滑移平均速度和位移Fig.11 Average speed and displacement of slope slip

完成整體滑移運(yùn)動(dòng)分析后，分別對滑坡各區(qū)域速度和位移變化進(jìn)行分析，監(jiān)測點(diǎn)的顆粒運(yùn)動(dòng)速度變化見圖12。其中，圖12(a)中點(diǎn)C峰值速度最高，約為6.9 m/s，點(diǎn)A、點(diǎn)B的速度波動(dòng)較為接近，峰值速度略低于整體平均速度，三個(gè)顆粒經(jīng)過10 s都停止運(yùn)動(dòng)；圖12(b)中三個(gè)點(diǎn)峰值速度都與整體滑破平均速度接近，期間點(diǎn)F在3.5 s時(shí)有強(qiáng)烈波動(dòng)，同頂端點(diǎn)運(yùn)動(dòng)停止時(shí)間相近；圖12(c)的底端點(diǎn)運(yùn)動(dòng)時(shí)間相比于中頂端的點(diǎn)，運(yùn)動(dòng)時(shí)間和整體滑破運(yùn)動(dòng)時(shí)間相同，經(jīng)過20 s后點(diǎn)H、點(diǎn)I靜止，而點(diǎn)G還在緩慢移動(dòng)。

圖12 滑坡各部分監(jiān)測點(diǎn)速度變化曲線Fig.12 Velocity change curve of monitoring points in each part of the landslide

對比各點(diǎn)速度變化圖，對圖13中各部分監(jiān)測點(diǎn)位移曲線變化進(jìn)行分析。其中，圖13(a)中點(diǎn)B位移量最多為12.0 m；圖13(b)中點(diǎn)E位移量略高于點(diǎn)B，在10 s時(shí)達(dá)到12.2 m后靜止；圖13(c)中點(diǎn)H、點(diǎn)I經(jīng)過20 s后停止移動(dòng)，位移量分別為13.4 m和14.4 m，而點(diǎn)G經(jīng)過22 s后位移量還在增加。

圖13 滑坡各部分監(jiān)測點(diǎn)位移變化曲線Fig.13 Displacement change curve of monitoring points in each part of the landslide

通過對各部分顆粒速度位移變化分析可以看出，每個(gè)位置上監(jiān)測點(diǎn)的速度變化中，靠近滑坡右側(cè)的點(diǎn)C、點(diǎn)F、點(diǎn)I的速度波動(dòng)頻率最明顯，而點(diǎn)C、點(diǎn)F兩點(diǎn)對比相同位置上其他點(diǎn)位移量都最少，說明滑坡右側(cè)更為陡峭，而滑坡左側(cè)地勢較為平坦，所以底端點(diǎn)G位移量最多，并在滑坡整體基本靜止后仍有微量位移?；马敹撕椭卸私?jīng)過10 s靜止，底端顆粒在形成堆積后繼續(xù)移動(dòng)，監(jiān)測點(diǎn)速度也在快速減緩后再次升高。

完成模擬后，對滑坡整體滑移峰值平均速度進(jìn)行計(jì)算[15]，通過計(jì)算得到峰值為6.2 m/s，較高于模擬結(jié)果值?？紤]到該理論忽略了巖體顆粒間的摩擦、碰撞的動(dòng)能消耗，因此誤差可以近似忽略，可以看出滑坡的反演結(jié)果和實(shí)際滑坡堆積結(jié)果基本相似，證明了實(shí)驗(yàn)的準(zhǔn)確性。

3.2 滑移堆積過程分析

根據(jù)滑坡的滑移運(yùn)動(dòng)過程模擬分析，滑坡前邊坡最大寬度約為19 m，滑坡的滑移區(qū)域和堆積區(qū)域有部分重疊。為了更直觀展示滑坡的堆積特征，按照高程將這片區(qū)域的顆粒劃分為不同的組，對分組顆粒進(jìn)行染色。圖14對比了滑坡前后的分組顆粒，能夠觀察到頂端顆粒、中端顆粒排序并未因滑坡發(fā)生而產(chǎn)生明顯變化，只有低端顆粒在運(yùn)動(dòng)時(shí)發(fā)生明顯混合。其中，頂端的兩個(gè)堆積區(qū)域和滑移區(qū)域發(fā)生重疊，左側(cè)部分顆粒被掩埋并隆起，中后端形成堆積平面。通過剖面圖也可看出堆積主要產(chǎn)生在后端，這與滑坡后實(shí)地調(diào)查堆積特征基本一致。

圖14 分組顆粒堆積分布Fig.14 Grouped particle packing distribution

通過滑坡運(yùn)動(dòng)軌跡和堆積特征分析，證明了離散元能夠有效地模擬礦山滑坡特征與運(yùn)動(dòng)趨勢，其數(shù)值模擬預(yù)測結(jié)能為滑坡的防治提供依據(jù)。

4 結(jié) 論

1) 本文將無人機(jī)調(diào)查同PFC3D三維數(shù)值建模相結(jié)合，通過無人機(jī)航攝構(gòu)建的滑坡DEM與復(fù)現(xiàn)原始地形，以Geometry Import指令在PFC3D內(nèi)構(gòu)建出離散元三維數(shù)值模型，為后續(xù)滑坡的數(shù)值模擬提供了數(shù)據(jù)支撐。

2) 通過在三軸壓縮模擬實(shí)驗(yàn)中加入BP神經(jīng)網(wǎng)絡(luò)的方法，完成對宏細(xì)觀參數(shù)的反演與標(biāo)定。經(jīng)過驗(yàn)證，將獲取的細(xì)觀參數(shù)帶入PFC3D內(nèi)獲取數(shù)據(jù)與室內(nèi)試驗(yàn)對比，結(jié)果能夠滿足要求。

3) 滑坡整體歷時(shí)20 s，在1.8 s時(shí)平均速度達(dá)到峰值5.6 m/s?；轮许敹祟w粒下落過程中呈現(xiàn)聚集狀，底端趨向于向左側(cè)運(yùn)動(dòng)，對比監(jiān)測點(diǎn)速度位移曲線運(yùn)動(dòng)，滑坡體右端較為陡峭。堆積體主體主要為中頂端顆粒，滑坡與堆積體有部分重疊，與實(shí)際堆積特征基本相符，證明本研究的模擬方法可以應(yīng)用于滑坡的運(yùn)動(dòng)特征研究。