冉林，馬鵬輝，彭建兵，孔嘉旭，張作蓬，楊炬，李澤坤，馬哲

（1. 長(zhǎng)安大學(xué)地質(zhì)工程與測(cè)繪學(xué)院，陜西 西安 710054；2. 西部礦產(chǎn)資源與地質(zhì)工程教育部重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710054；3. 中國(guó)有色金屬工業(yè)西安勘察設(shè)計(jì)研究院有限公司，陜西 西安 710054）

0 引言

我國(guó)是世界上最大的黃土分布區(qū)，總面積約6.3×105km2，占國(guó)土總面積的6.6 %[1]。由于黃土的類卡斯特結(jié)構(gòu)，導(dǎo)致黃土具有強(qiáng)烈的水敏性，天然狀態(tài)下的黃土具有極高的強(qiáng)度，遇水后，強(qiáng)度急劇降低，易發(fā)生軟化變形，進(jìn)而形成黃土地質(zhì)災(zāi)害。據(jù)統(tǒng)計(jì)，90 %以上的黃土地質(zhì)災(zāi)害發(fā)生與水有關(guān)[2?3]。甘肅黑方臺(tái)是我國(guó)最著名的水致黃土滑坡研究場(chǎng)地之一，被稱為黃土滑坡研究的天然試驗(yàn)場(chǎng)。從1968 年至今，幾十年農(nóng)業(yè)漫溉使黑方臺(tái)地下水位大幅抬升，塬邊黃土滑坡頻繁發(fā)生[4]，引起了國(guó)內(nèi)外學(xué)者廣泛關(guān)注，由此開展了大量的研究。例如：在滑坡類型及分布規(guī)律方面，任曉虎[5]基于多種滑坡識(shí)別指標(biāo)對(duì)黑方臺(tái)黃土滑坡發(fā)育早期特征識(shí)別進(jìn)行研究。許領(lǐng)等[6]等利用IKONOS 影像與DEM 疊加重現(xiàn)了滑坡場(chǎng)景，對(duì)黑方臺(tái)滑坡進(jìn)行了系統(tǒng)的分類；在成因機(jī)理方面，張茂省等[7?8]、王家鼎等[9]、潘攀[10]、藺曉燕[11]對(duì)灌溉水誘發(fā)黃土滑坡的形成機(jī)理進(jìn)行了較系統(tǒng)研究；在風(fēng)險(xiǎn)評(píng)估方面，彭大雷[12]在研究歸納了各類型災(zāi)害的形成條件、致災(zāi)因子、成災(zāi)模式等的基礎(chǔ)上，提出黃土滑坡潛在隱患地質(zhì)識(shí)別方法。薛強(qiáng)等[13]基于多期DEM 數(shù)據(jù)，利用GIS 技術(shù)，通過滑坡體積與速度乘積計(jì)算滑坡強(qiáng)度，對(duì)黑方臺(tái)南塬典型滑坡進(jìn)行了危險(xiǎn)性定量評(píng)估。然而隨著滑坡的頻繁發(fā)生，當(dāng)?shù)鼐用癜踩珕栴}異常突出，而黑方臺(tái)黃土滑坡的運(yùn)動(dòng)特征方面的研究直接關(guān)乎著滑體最終的致災(zāi)范圍和災(zāi)害效應(yīng)，影響當(dāng)?shù)鼐用裆?cái)產(chǎn)。

針對(duì)黃土滑坡的運(yùn)動(dòng)特征研究，數(shù)值模擬是最常用也是最可視化的手段。例如：王思源[14]運(yùn)用FLAC3D有限差分模擬軟件，分析了黑方臺(tái)焦家滑坡的實(shí)際變形破壞過程；胡煒等[15]在實(shí)際測(cè)試及工程地質(zhì)測(cè)試的基礎(chǔ)上，基于有限差分法的Sassa 運(yùn)動(dòng)學(xué)模型模擬焦家崖頭滑坡運(yùn)動(dòng)全過程，并分析滑坡的速度變化規(guī)律；賈俊等[16]以離散元模型對(duì)焦家崖頭13 號(hào)滑坡為例，研究地下水對(duì)滑坡作用機(jī)制及灌溉型黃土滑坡破壞過程和高速遠(yuǎn)程運(yùn)動(dòng)學(xué)特征;趙紀(jì)飛等[17]以黑方臺(tái)焦家滑坡為例，基于野外調(diào)查和室內(nèi)試驗(yàn)，通過Sassa 運(yùn)動(dòng)學(xué)模型模擬該滑坡泥流運(yùn)動(dòng)過程?？紤]到黑方臺(tái)地區(qū)的黃土滑坡具有很強(qiáng)的連續(xù)性和流動(dòng)性，而Massflow 基于深度積分的連續(xù)介質(zhì)力學(xué)方法，將運(yùn)動(dòng)物體視作連續(xù)流體，可用來模擬碎屑流、滑坡、泥石流、山洪等地質(zhì)災(zāi)害動(dòng)力演化全過程，且建模效率高、計(jì)算速率快，故采用該軟件反演黑方臺(tái)黃土滑坡較為適用，已在多個(gè)案例[12，18?21]得到應(yīng)用。

文中以黑方臺(tái)2019 年10 月5 日發(fā)生在黨川村附近的一起黃土滑坡為研究對(duì)象，通過野外科學(xué)考察，分析該滑坡的基本特征，并采用Massflow 軟件反演該滑坡運(yùn)動(dòng)全過程，為黑方臺(tái)地區(qū)滑坡風(fēng)險(xiǎn)評(píng)估及防災(zāi)減災(zāi)提供一定的科學(xué)參考。

1 研究區(qū)區(qū)域概況

黑方臺(tái)位于甘肅省永靖縣鹽鍋峽鎮(zhèn)境內(nèi)，地處黃河與湟水河交匯處上游，距蘭州市約40 km（圖1），屬于黃河Ⅳ級(jí)階地，臺(tái)塬總面積約14 km2，塬面東西向長(zhǎng)約7.7 km，南北向?qū)捈s2.5 km。臺(tái)塬被虎狼溝切割為黑臺(tái)和方臺(tái)，東側(cè)為黑臺(tái)，面積約9 km2，西側(cè)為方臺(tái)，面積約2 km2。臺(tái)塬出露地層自上而下為：全新統(tǒng)崩滑堆積物（）；上更新統(tǒng)馬蘭黃土（）；中更新統(tǒng)粉質(zhì)黏土（）；中更新統(tǒng)砂卵石層（）；下白堊統(tǒng)河口群粉砂質(zhì)泥巖（K1hk）（表1）。

表1 研究區(qū)地層巖性Table 1 Stratigraphic lithology of the study area

研究區(qū)內(nèi)黃土滑坡具有群體性分布特征[6]（圖2），根據(jù)黑方臺(tái)地區(qū)內(nèi)物質(zhì)組成和運(yùn)動(dòng)特征，將黑方臺(tái)滑坡分為四大類：（1）黃土—錯(cuò)落型；（2）黃土—碎屑流型；（3）黃土—泥流型；（4）黃土基巖—滑移型。

圖2 黑方臺(tái)滑坡分布圖Fig.2 Distribution of landslides in the Heifangtai platform

2 滑坡基本特征

“10·5”滑坡發(fā)生于2019 年10 月5 日，位于甘肅黑方臺(tái)黨川村，據(jù)當(dāng)?shù)鼐用窠榻B，該滑坡運(yùn)動(dòng)時(shí)間極短，不超過1 min。該滑坡為黃土—碎屑流型滑坡，基本沒有水的影響，多為崩塌型的黃土內(nèi)部滑動(dòng)，發(fā)生于塬邊裂縫發(fā)育和內(nèi)部存在軟弱結(jié)構(gòu)面的Q3黃土中，通常情況下前期次滑坡增加坡體的臨空條件，并在坡體后緣形成一定數(shù)量的拉張裂縫，為黃土—碎屑流型滑坡的發(fā)生創(chuàng)造條件，坡體失穩(wěn)破壞之后，在運(yùn)動(dòng)過程中逐漸破碎解體?！?0·5”滑坡存在多期性、隱蔽性、突發(fā)性及影響范圍廣等特點(diǎn)（圖3）。

圖3 滑坡基本特征Fig.3 The characteristics of the landslide

2.1 地形地貌

滑坡區(qū)在地貌單元上屬于黃土高原丘陵溝壑區(qū)，整體位于黃河Ⅳ級(jí)基座階地，滑坡所處黨川段為目前以來滑坡發(fā)育最活躍的區(qū)域?；聟^(qū)斜坡坡體上陡下緩，滑坡后壁坡度范圍為60°～85°，滑坡前緣為較平緩開闊農(nóng)田，坡度范圍20°～25°，平均坡度約50°，這種后緣高陡、前緣開闊的地形特點(diǎn)，為滑坡的發(fā)生提供了良好的臨空條件。

2.2 滑體物質(zhì)及坡體結(jié)構(gòu)

“10·5”滑坡的滑體物質(zhì)成分主要由風(fēng)成黃土組成，位于滑坡上部，厚度約35 m，巖性為馬蘭黃土，淺黃色，土體松散，遇水后土體體積含水量迅速增大，導(dǎo)致黃土黏聚力大幅降低，發(fā)生軟化，坡體穩(wěn)定性降低，滑動(dòng)面在黃土層非飽和帶，剪出口在斜坡中上部，多在坡腳形成堆積。坡體及周邊出露地層主要有馬蘭黃土、粉質(zhì)黏土、砂卵石層、粉砂質(zhì)泥巖[22?24]（圖4）。

圖4 “10·5”滑坡剖面圖Fig.4 The profile of the “10·5” landslide

2.3 滑坡形態(tài)與規(guī)模

根據(jù)野外實(shí)地調(diào)查，該滑坡整體平面形態(tài)呈長(zhǎng)椅型，后壁近似圓弧狀，后壁寬約110 m，順坡長(zhǎng)約355 m，后緣高程1 710 m，前緣高程1 605 m，相對(duì)高差105 m，坡體后緣后退約11 m，滑源區(qū)滑體平均厚度30 m，滑體在坡腳形成堆積區(qū)，呈扇形展布，堆積區(qū)地形較為平緩，滑體極為松散，滑坡前緣堆積區(qū)寬度不一，范圍在70～140 m（圖3），滑坡區(qū)域面積為700 m2，體積約2.7×104m3，屬于小型黃土滑坡。

2.4 滑坡變形特征

坡體后緣發(fā)育數(shù)條拉張裂縫，平均寬度5～20 cm，走向平行于坡邊，豎向上呈上寬下窄的特點(diǎn)，深度受坡高、坡度等因素控制，分布形態(tài)為直線型，呈現(xiàn)出漸近式后退的趨勢(shì)，有明顯錯(cuò)臺(tái)，其錯(cuò)臺(tái)最大下錯(cuò)量可達(dá)90 cm，平均下錯(cuò)大于20 cm。選取一條具有明顯變化的裂縫1，通過幾期不同時(shí)間段的對(duì)比，裂縫橫向長(zhǎng)度與寬度逐漸增大，在201912 期前，斜坡土體沿裂縫1 貫通發(fā)生崩塌滑移，并在坡頂后部發(fā)育多條拉張裂縫（圖5）。

圖5 滑坡后緣錯(cuò)臺(tái)及局部裂縫變化過程Fig.5 Dislocation of the trailing edge of the landslide and the change of local fractures

2.5 滑坡運(yùn)動(dòng)過程特征

據(jù)野外考察發(fā)現(xiàn)，滑源區(qū)土體高度臨空，后緣受拉張應(yīng)力影響，裂縫不斷擴(kuò)張，最終形成滑動(dòng)面發(fā)生剪切破壞（圖6），滑源區(qū)土體在重力作用下失穩(wěn)破壞以235°方向開始滑動(dòng)，坡下地形無約束，滑坡在運(yùn)動(dòng)過程中未發(fā)生較大的改道行為，滑坡東臨山脊，部分滑體受到右側(cè)山脊的碰撞，滑至底部平緩地區(qū)時(shí)，向西南方向發(fā)展?；录舫隹谖恢幂^高，滑體儲(chǔ)有較大勢(shì)能，滑坡前緣到達(dá)斜坡底部之后繼續(xù)向前運(yùn)動(dòng)，穿過底部農(nóng)田和灌溉渠道直至停止，滑坡前緣最遠(yuǎn)處距公路60 m（圖7）。

圖6 滑坡發(fā)育過程Fig.6 The developping process of landslide

圖7 滑坡運(yùn)動(dòng)過程路徑Fig.7 Motion path of the landslide

3 “10·5”滑坡運(yùn)動(dòng)過程反演

3.1 運(yùn)動(dòng)模型建立

采用Massflow 軟件進(jìn)行滑坡特征反演，該軟件是基于深度積分的連續(xù)介質(zhì)力學(xué)方程而建立的有限差分法數(shù)值軟件，能夠很好地考慮運(yùn)動(dòng)過程中地形地貌對(duì)其運(yùn)動(dòng)的影響。該軟件把復(fù)雜地質(zhì)模型簡(jiǎn)化為上下界面，計(jì)算公式如下[18,21]：

式中：ρ1——流動(dòng)層密度/（kg·m?3）；

τ1——上部流動(dòng)層對(duì)底部的剪力/kPa；

z10——頂部流動(dòng)層的基底邊界/m；

E1——底部穩(wěn)定層的夾帶率；

u1、v1——材料邊界速率/（m·s?1）；

h1——流體厚度/m；

g——重力加速度，取9.8 m/s2。

該軟件內(nèi)有Coulomb、Manning、Voellmy 三種計(jì)算模型，其中Coulomb 模型適用于滑坡、巖崩等碎屑流（固相）、Manning 模型適用于潰壩、洪水等水力學(xué)災(zāi)害（液相）、Voellmy 模型適用于泥漿，泥石流等災(zāi)害（固液相）。文中采取Coulomb 計(jì)算模型，進(jìn)行初步模擬計(jì)算，Coulomb 模型其主要計(jì)算方程如式（2），關(guān)于軟件的基本力學(xué)方程，可參考?xì)W陽(yáng)朝軍等的相關(guān)研究[18]：

式中：τb——基底摩檫力/kPa；

c、φ——土體黏聚力和內(nèi)摩擦角/（°）；

g——重力加速度，取9.8 m/s2；

h——滑坡豎向厚度/m；

λ——基底液化系數(shù)，取值為0～1。

關(guān)于滑坡前后地形數(shù)據(jù)及影像的獲取，文中通過無人機(jī)航測(cè)獲得，將所獲DEM 數(shù)據(jù)導(dǎo)入ArcGIS10.7 中進(jìn)行處理，得到滑坡物源區(qū)范圍，將滑坡地形和物源區(qū)柵格數(shù)據(jù)轉(zhuǎn)換為Massflow 能夠識(shí)別的ASCII 格式。

3.1.1 參數(shù)選取

模擬所需參數(shù)有土體密度（ρ）、黏聚力（c）、基底液化系數(shù)（λ）、巖土體內(nèi)摩擦角（φ），參考黑方臺(tái)地區(qū)已有室內(nèi)實(shí)驗(yàn)數(shù)據(jù)及參數(shù)[12,25]，將模擬所需非關(guān)鍵參數(shù)土體密度與內(nèi)摩擦角視為定值。對(duì)于黏聚力與基底液化參數(shù)，這兩者主要影響著滑坡的堆積形態(tài)和堆積厚度，滑坡的基底液化系數(shù)和黏聚力通過控制變量法確定，主要參數(shù)取值見表2。

表2 模擬計(jì)算參數(shù)Table 2 Simulation parameters

通過反演計(jì)算得到黏聚力和基底液化系數(shù)最佳取值，滑體物質(zhì)主要是黃土，滑床基本為老的滑坡堆積體，因此滑坡黏聚力可簡(jiǎn)化為一個(gè)定值，黃土浸水后其黏聚力大幅下降[26?27]，故黏聚力取值范圍設(shè)為0～7 kPa。第一輪反演將基底液化系數(shù)視為定值（取0.4），黏聚力在0.5 kPa 的梯度取值條件進(jìn)行反演試算，模擬結(jié)果與滑坡實(shí)際范圍進(jìn)行對(duì)比，確定黏聚力值；第二輪反演在黏聚力為5.5 kPa 時(shí)，基底液化系數(shù) λ以0.1 的梯度進(jìn)行反演試算，發(fā)現(xiàn)在0.6～0.8 區(qū)間結(jié)果較好，進(jìn)而將參數(shù)以0.01 的梯度細(xì)化，得到更精確的基底液化系數(shù)λ。

3.2 反演結(jié)果分析

根據(jù)表2 中模擬計(jì)算參數(shù)進(jìn)行計(jì)算，模擬滑坡不同時(shí)間段的運(yùn)動(dòng)特征及滑坡體的堆積范圍和厚度，直至運(yùn)動(dòng)停止。該滑坡為多期次滑坡，所獲滑坡前后地形數(shù)據(jù)時(shí)間分別為2019 年5 月及2019 年12 月，通過前后DEM數(shù)據(jù)所獲滑源區(qū)土體大于“10·5”滑坡實(shí)際滑源區(qū)，并且期間滑坡上部土體存在多次崩落的可能性，崩落土體因重力作用落于坡面上，改變了滑坡下墊面微地形，使得堆積區(qū)右側(cè)地形高于左側(cè)，滑體向下運(yùn)動(dòng)過程中向兩側(cè)擴(kuò)展，因滑坡右側(cè)為一山脊，故在滑坡模擬過程中，滑源區(qū)土體在重力作用下向下運(yùn)動(dòng)時(shí)，東側(cè)滑體受到右側(cè)山脊的碰撞，發(fā)生一定程度的破碎、折返運(yùn)動(dòng)，隨即向西南方向發(fā)展運(yùn)動(dòng)堆積，模擬所得結(jié)果與實(shí)際堆積范圍出現(xiàn)些許偏差。

3.2.1 堆積厚度分析

通過對(duì)滑坡各時(shí)間段運(yùn)動(dòng)過程的模擬，得到了不同時(shí)間段的堆積范圍（圖8），滑坡體初始厚度最大約16 m，滑坡啟動(dòng)后，快速?zèng)_向坡腳和斜坡兩側(cè)，并最終堆積在滑坡前緣耕地處。圖8 中（b）可看出t=5 s 時(shí)滑坡已啟動(dòng)，滑源逐漸向坡腳滑動(dòng)，此時(shí)滑坡最大堆積厚度為7.08 m，t=10 s 時(shí)滑坡主體部分已基本抵達(dá)坡腳，最大堆積厚度為3.36 m，t=15 s 時(shí)，滑坡在滑坡前緣耕地處堆積，t>15 s 時(shí)滑坡體在前緣耕地處緩慢流動(dòng)至停止，45 s 左右時(shí)滑坡已基本停止運(yùn)動(dòng)。

圖8 滑坡不同時(shí)刻堆積范圍Fig.8 The simulation results of landslide accumulation boundary at different times

滑坡運(yùn)動(dòng)是一個(gè)動(dòng)態(tài)運(yùn)動(dòng)過程，為定量分析滑坡堆積厚度分布特征，選取t=45 s 的模擬結(jié)果，在滑坡的主滑方向、滑源區(qū)與滑坡前緣三個(gè)位置提取堆積厚度，與真實(shí)堆積情況對(duì)比。根據(jù)滑坡前后高程數(shù)據(jù)對(duì)比，獲得滑坡真實(shí)堆積厚度分布（圖9）。對(duì)于主滑方向的1-1′剖面，滑移距離在約355 m 處基本達(dá)到收斂，滑坡運(yùn)動(dòng)方向的反演堆積厚度與實(shí)際堆積特征相符，堆積區(qū)堆積厚度整體趨勢(shì)為先增大減小，平均堆積厚度逐漸變小，最大堆積厚度約為0.82 m。對(duì)于3-3′剖面，因滑坡前緣地形無約束且滑體流動(dòng)性較強(qiáng)，在滑坡物源體積一定的情況下，滑體向前運(yùn)動(dòng)的同時(shí)向兩側(cè)擴(kuò)散堆積，主要以散狀堆積在坡腳，呈扇形展布，實(shí)際平均堆積厚度約0.63 m，模擬平均堆積厚約0.49 m。由2-2′剖面看出，在滑坡滑源區(qū)附近，坡體崩塌下滑，故滑源區(qū)堆積厚度為負(fù)值，2-2′剖面真實(shí)下滑最大深度約14.86 m，反演結(jié)果與真實(shí)下滑深度一致。

圖9 滑坡剖面位置及堆積厚度圖（t=45 s）Fig.9 Location of typical profile and accumulate thickness (t=45 s)

通過三個(gè)剖面堆積厚度的比較，可發(fā)現(xiàn)在主滑方向與滑源方向上的堆積厚度變化和實(shí)際情況一致，滑坡前緣模擬堆積厚度與實(shí)際情況存在誤差，這是因?yàn)榛露逊e區(qū)的整體厚度均小于1 m，在小范圍厚度上誤差體現(xiàn)比較明顯，同時(shí)也因滑坡運(yùn)動(dòng)路徑受地形和前期堆積體的限制，運(yùn)動(dòng)過程中發(fā)生了碰撞折返，故在模擬滑坡運(yùn)動(dòng)時(shí)，滑體到達(dá)實(shí)際堆積區(qū)范圍后，繼續(xù)向西南方向運(yùn)動(dòng)堆積，模擬所得堆積區(qū)范圍略大于實(shí)際堆積區(qū)范圍，使得滑坡前緣處模擬堆積厚度小于實(shí)際堆積厚度。

3.2.2 運(yùn)動(dòng)速度分析

圖10 為該滑坡不同時(shí)刻運(yùn)動(dòng)速度模擬結(jié)果。為獲取滑坡不同部位在不同時(shí)刻的運(yùn)動(dòng)速度，在滑坡主滑方向上布置速度監(jiān)測(cè)點(diǎn)1?5 號(hào)，得出1?5 號(hào)監(jiān)測(cè)點(diǎn)滑坡運(yùn)動(dòng)速度變化圖。

圖10 不同時(shí)刻滑坡運(yùn)動(dòng)速度Fig.10 The movement speed of the landslide at different times

滑坡自啟動(dòng)到停止運(yùn)動(dòng)總用時(shí)45 s，可將其運(yùn)動(dòng)堆積過程分為啟動(dòng)加速-穩(wěn)定加速-減速堆積三個(gè)階段（圖11）：0～7.2 s 內(nèi)為啟動(dòng)加速階段；7.2～11.3 s 為穩(wěn)定加速階段；11.3～45 s 為減速堆積階段，滑坡動(dòng)能逐漸耗盡，最終停止。其中啟動(dòng)加速階段與穩(wěn)定加速階段時(shí)間較短，僅占運(yùn)動(dòng)總時(shí)間的25.1 %，說明“10·5”滑坡具有突發(fā)性，坡體上部迅速崩塌滑動(dòng)，在極短的時(shí)間內(nèi)速度達(dá)到最大值，減速堆積階段占據(jù)滑坡運(yùn)動(dòng)主要部分，約74.9 %，說明滑坡具有較強(qiáng)的流動(dòng)性。

圖11 運(yùn)動(dòng)速度-時(shí)間變化曲線Fig.11 Relationship between the movement speed and time

從5 個(gè)速度監(jiān)測(cè)點(diǎn)獲得的數(shù)據(jù)來看，到達(dá)每個(gè)監(jiān)測(cè)點(diǎn)的最大速度數(shù)值相差較小，滑坡各點(diǎn)運(yùn)動(dòng)速度均表現(xiàn)為先快速增大后迅速減小。坡下橫向地形無約束，能量消散快，故滑坡運(yùn)動(dòng)速度變化迅速。對(duì)于1 號(hào)監(jiān)測(cè)點(diǎn)，坡體上部土體開始崩塌失穩(wěn)，滑動(dòng)開始約1 s 運(yùn)動(dòng)速度達(dá)到最大值16 m/s，因坡體上部存在些許土體掉落，在減速階段速度未能減至0；2 號(hào)監(jiān)測(cè)點(diǎn)上，滑坡滑源物質(zhì)經(jīng)4 s 左右達(dá)到最大速度23 m/s，17 s 時(shí)速度趨于0；3號(hào)監(jiān)測(cè)點(diǎn)上，滑坡歷經(jīng)7 s 達(dá)到最大速度26 m/s，25 s 時(shí)速度趨于0；4 號(hào)監(jiān)測(cè)點(diǎn)上，滑坡歷經(jīng)11 s 達(dá)到最大速度24.6 m/s，30 s 時(shí)速度趨于0；5 號(hào)監(jiān)測(cè)點(diǎn)位于堆積區(qū)前緣，滑坡動(dòng)能已耗散部分，故此處最大速度較低為18.5 m/s，在45 s 時(shí)趨于0。結(jié)合運(yùn)動(dòng)速度-時(shí)間變化曲線圖可以看出，滑坡從后緣至前緣處逐步達(dá)到靜止?fàn)顟B(tài)，且在滑坡不同部位，滑坡運(yùn)動(dòng)速度從加速直至降低為0 的時(shí)間間隔都比較小，進(jìn)一步說明了“10·5”滑坡運(yùn)動(dòng)速度快。

4 結(jié)論

文中通過野外考察研究，對(duì)黑方臺(tái)黨川村“10·5”滑坡的基本特征及運(yùn)動(dòng)特征進(jìn)行了一定的討論，并利用Massflow 數(shù)值模擬軟件反演滑坡運(yùn)動(dòng)全過程。得到以下結(jié)論：

（1）“10·5”滑坡體積約為2.7×104m3，主滑方向?yàn)?35°，最大滑距355 m，具有突發(fā)性，影響范圍廣的特點(diǎn)。

（2）根據(jù)滑坡運(yùn)動(dòng)速度可將“10·5”滑坡運(yùn)動(dòng)堆積過程分為啟動(dòng)加速-穩(wěn)定加速-減速堆積三個(gè)階段，自啟動(dòng)到停止整個(gè)過程僅持續(xù)45 s，其中滑坡最大運(yùn)動(dòng)速度可達(dá)26 m/s，堆積區(qū)最大堆積厚度為0.82 m。

（3）通過1—5 號(hào)速度監(jiān)測(cè)點(diǎn)的數(shù)據(jù)可看出，“10·5”滑坡運(yùn)動(dòng)速度呈先快速增大后迅速減小的趨勢(shì)，所用時(shí)間極短。

（4）經(jīng)過Massflow 數(shù)值模擬獲得的滑坡危險(xiǎn)范圍較大于真實(shí)范圍，存在些許偏差，結(jié)果偏于保守。