青藏高原因其复杂的地质构造背景、陡峻的地形地貌条件，高速岩质滑坡广泛发育于区内的高山峡谷中，且往往具有高位、高剪出口等特性^[1]。在地形陡峭的峡谷地区，高速岩质滑坡因其势能大、运动速度快的特点，还易造成堵江事件以及对上下游的淹没、洪水灾害^[1-5]。近年来，由于降雨、地震等原因，此类高速岩质滑坡频发，给社会带来了巨大的损失。由于高速岩质滑坡特殊复杂的发育条件和运动特征，无论是滑坡本身还是其引发的灾害链均严重威胁山区居民的生命财产安全，而目前对于该类滑坡运动特性并没有统一明确的分析方法，因此，急需开展高速岩质滑坡运动特性的深入研究。

2018年10月10日(“10·10”滑坡)和11月3日(“11·3”滑坡)西藏自治区江达县白格村在一个月内连续两次发生高位山体滑坡，滑坡碎屑体堵塞金沙江，形成堰塞湖。目前已有不少学者针对白格滑坡开展了研究工作^[6-18]，然而这些研究多集中于白格滑坡地质条件和形成机制的分析，或集中于白格滑坡的运动、堵江过程模拟，白格滑坡在各运动阶段的不同运动速度及状态尚未查清。

“10·10”白格滑坡先后经历了主体失稳破坏、牵引区启动、高速临空滑跃、冲击对岸、折返相撞及水砂射流、堆积坝次级滑移6个运动阶段，并于高速运动后在金沙江对岸造成严重冲击和水砂射流破坏，形成最高达176 m的明显冲击痕迹。因此，依据对白格滑坡的现场调查，提取和量测滑坡运动参数、冲击区几何参数、滑坡体积等关键数据，基于滑坡碎屑体冲击高度，反演和计算“10·10”白格滑坡在各阶段的运动速度。研究结果有助于加深高速岩质滑坡运动特性的认识，可为该类滑坡的预测和防治提供理论依据。

1 研究区概况 1.1 地质条件

白格滑坡(31.078887°N，98.698023°E)位于中国西藏自治区江达县波罗乡白格村，距波罗乡20 km，发生于金沙江右岸。滑坡所在区位于金沙江缝合带，区内主要断裂总体走向为NNW-SSE，波罗- 木协断裂是距离白格滑坡后缘最近的断裂(如图 1)。滑坡周边区域出露的地层岩性按昌都-江达地层分区主要有三叠系金古组(T₃jn)灰岩带，元古界雄松群(P_txn^a)斜长片麻岩组，燕山期戈坡超单元(ηγ₅^2b)和则巴超单元(γδ₅^2a)的花岗岩组，三叠系下逆松多组(T₃x)的碳酸岩盐和碎屑岩段(按义敦- 巴塘地层分区)，并在白格滑坡后缘的波罗- 木协断裂发育有华里西期金沙江超镁铁质和蛇纹岩带(φ_ω4)^[19]，燕山期的花岗闪长岩(γδ₅^2a)也出露于滑坡南侧，该酸性岩石也存在局部侵入。在滑动区左侧边界的局部地层中也发现可能为燕山期的花岗岩斑岩。研究区上下部岩性及风化程度不一，滑坡主要发生于上部全风化蚀变片麻岩带，中部以强风化片麻岩夹炭质板岩为主，下部以中风化的片麻岩为主，如图 1所示。

1.2 地形地貌条件

研究区为强烈侵蚀区，海拔差异超过1 720 m(2 880~4 600 m)。受河流的强烈侵蚀影响，金沙江以V形河谷为主。滑坡区地处金沙江河谷右岸，河道高程为2 880 m，滑坡后缘顶部高程3 720 m，相对高差840 m。滑坡区岸坡较为陡峻，平均坡度30°。滑坡表面由上到下发育3级平台：一级平台范围最广，包含白格村等区域，海拔高程3 450~3 550 m；二级平台次之，海拔高程3 100~3 200 m；三级平台靠近金沙江，位于滑坡剪出口右侧，形态特征不明显，海拔高程2 940~2 970 m。

1.3 气候条件

研究区气温在-15~28 ℃变化，年平均温度为4.5 ℃。据江达县气象站资料，该县多年平均降水量650 mm，最大年降水量1 067.7 mm，最大月降水量229.5 mm。年内降水分布不均，干湿季节分明。10月—次年5月为旱季，降水稀少，降雨量约占全年降雨量的12.8%。6—9月为雨季，降雨量约占全年降水量的87.2%。该县总体上雨量较为贫乏，但由于降水量时空分布不均匀，局地暴雨现象时有发生，成为泥石流、滑坡、崩塌等地质灾害的重要诱发因素。2018年滑坡区邻近3个乡镇降雨量显示，6—9月当地降雨量明显增加(如图 2)。因此，推断降雨对“10·10”白格滑坡的发育和前置斜坡变形应有一定贡献作用。

2 滑坡运动特征 2.1 滑坡运动分区

“10·10”白格滑坡主要分为启动区、堆积区、冲击区3区。“10·10”白格滑坡启动区沿坡面纵向距离较长，海拔高度分布于2 980~3 720 m。如图 3、5所示，滑坡启动区可进一步细分为3个亚区：1)牵引区(1-1)。位于一级平台之上，海拔为3 500~3 720 m，体积458万m³。该区靠近波罗- 木协断层，且风化非常严重，岩性主要以强蚀变片麻岩为主，性状近似为土体，滑坡体非常松散；2)主滑区(1-2)。该区主要分布于一、二级平台之间，海拔为3 100~3 500 m, 体积1 328万m³。主滑区主要夹于两侧基岩区之间，形成楔形体构造，岩体主要以片麻岩夹炭质板岩为主，风化程度相对弱于牵引区；3)锁固段(1-3)。主要位于白格滑坡三级平台区，海拔为2 980~3 100 m，体积80万m³。该区岩体主要由相对完整的片麻岩组成，结构面不发育。因此，在主滑区重力压力的作用下，锁固段首先发生失稳破坏，并诱导主滑区松散滑坡体一起向下滑动。

“10·10”白格滑坡堆积区具有左岸高、右岸低的堆积形态。左岸堆积体部分主要由两部分滑体组成，岩性分别为弱风化片麻岩块和碳质板岩细颗粒，对应着锁固段和主滑区的滑体；右岸堆积体部分由强蚀变片麻岩蛇纹岩颗粒组成，对应着牵引区滑体，在堆积区形成之初，其内部持续发生着次级滑移。根据现场精细调查，堆积区堰塞坝坝体长1 200 m，平均宽380 m，总面积38万m²。坝体溃口点海拔2 931 m，厚度47 m，最高点海拔为3 005 m，最高厚度为80 m，平均厚度32 m，堆积区滑坡方量达1 866万m³。

“10·10”白格滑坡冲击区主要影响金沙江左岸(四川侧)，如图 4、5所示，分为碎屑流冲击区(3-1)、水砂射流区(3-2)、淤泥污染区(3-3)3个亚区。碎屑流冲击区位于金沙江左岸上游段，金沙江左岸(四川侧)为碎屑体主要冲击区，其面积达到8.4万m²，分布形态为中间冲击高度向两侧递减，在冲击坡面上，残存有非常明显的滑坡碎屑物；金沙江右岸(西藏侧)为折返次级冲击区，面积仅为2.1万m²。水砂射流区位于金沙江左岸下游段，为碎屑体撞击激起的水砂造成的冲击破坏。

2.2 滑坡运动过程

基于现场调查，白格滑坡运动过程可分为主体失稳破坏、牵引区启动、高速临空滑跃、冲击对岸、折返相撞及水砂射流、堆积坝次级滑移6个阶段：1)主体失稳破坏, 滑坡启动区锁固段在重力作用下重压失稳，诱导上部主滑区整体启动发生滑移；2)牵引区启动, 主滑区滑移，牵引区下部临空面形成，由于失去主滑区滑体的支撑作用，最终失稳破坏；3)高速临空滑跃, 邻近江边，滑体于拔河高度约100 m处剪出飞跃河道；4) 冲击对岸阶段, 滑坡体冲击金沙江对岸四川侧，沿坡面爬升到海拔3 056 m；5)折返相撞及水砂射流, 滑体折返，与牵引区碎屑体对撞于金沙江上，激起水浪发生水砂射流；6)堆积坝次级滑移, 堆积坝形成后，由于高度差，堆积坝再次发生次级滑移。

3 滑坡运动速度计算 3.1 滑坡运动基本特征

滑坡体在高速运动的过程中，由于方量大、运动速度快，滑坡体在对岸山体坡面上留下明显的冲击痕迹。滑坡冲高在对面山体上留下长约943 m的弧形痕迹(碎屑流冲击区)，弧形最高处位于滑坡中部，冲起高度176 m(海拔3 056 m)，往冲击区两端走，冲起高度逐步降低，到两端时，冲起高度为0。另外，部分滑体滑入金沙江，造成高速水砂射流，在对岸(四川岸)也留下了明显痕迹。

3.2 滑坡启动速度计算

由于锁固段岩体主要由相对完整的片麻岩组成，结构面不发育。经现场调查发现该区滑床基本未残留碎屑，但是新鲜擦痕和基岩剪断破坏特征显著，呈现出较强的应力集中和剪切特征，并且发现剪出口下部存在未扰动区，见图 6。因此，推断锁固段以脆性剪断方式启动，和主滑区岩土体以一定初速度抛射而出，并且在重力作用下加速，高速与左岸(四川岸)相撞后逆坡爬行至顶部，随后折返与牵引区岩土体相撞于金沙江上。因此，白格滑坡速度计算必须考虑初速度。

为便于研究滑坡体起动机制，对白格滑坡的几何条件及边界条件进行简化(如图 7)，可见在重力及其他外力的长期作用下，滑坡体产生拉裂缝CD, AB段为锁固段，处于滑体剪出口。由主滑区DA段和阻滑区AB段构成的滑面，即将来产生的滑床面位置。白格滑坡在临界起动状态时，AB锁固段部位应力小于或接近该处岩石抗剪强度，此时AB锁固段存在变形但未剪滑，由AB段起到锁固作用。因此，在白格滑坡启动之前，边坡内的岩体应力逐渐增加，能量逐渐聚集。当锁固段AB段应力大于该处岩石抗剪强度的瞬间，抗剪应力骤然按该处岩体“峰残强降率”降低了50%~60%(有时甚至更大)，锁固段脆性剪断，释放出的巨大应力使滑坡体由静态转为动态，在滑体启动开始的瞬间出现相当高的“启程剧发速度”，向斜坡重力方向迅猛下滑，而使岩体获得初始速度，即为白格滑坡锁固段和主滑区的“启动速度”，这就是启程剧动的临床峰残强降加速机制^[20]。

根据临床峰残强降加速机制，应用断裂力学原理，分析斜坡累进性破坏过程中锁固段岩体的剪断释能效应，导出平面应力状态下峰残强降剧动速度公式^[20]为

$ v_{0}=\sqrt{\frac{{\rm{ \mathit{ π} }} g G_{V} \gamma}{2 E}} \cos \alpha\left(\tan \varphi_{\mathrm{p}}-\tan \varphi_{\mathrm{r}}\right) $

(1)

式中：根据白格滑坡工程地质特征及室内岩石力学试验结果，G_V为单位宽度滑体体积，此处取1 408万m³；α为滑动面倾角，此处取28°；E为岩石弹性模量，取24 GPa；φ_p、φ_r分别为启动区片麻岩峰值抗剪强度和残余抗剪强度，取为36.7°和30°; γ为岩石容重，取24.1 kg/m³。可得白格滑坡启动初速度为2.2 m/s。

3.3 滑坡撞击对岸高程计算

当主滑区和阻滑区滑体启动后，该部分滑体以一定速度从高程点2 980 m处飞跃河岸，撞击四川岸后逆坡爬行至高程点3 056 m后滑下。因此，要利用四川岸滑坡冲击高度计算滑坡行程速度，必须知道滑坡撞击四川岸高程。

目前，对于体积为百万至千万立方米巨型滑坡的速度估算，国际上普遍采用Scheidegger^[21]提出的一种滑坡计算公式：

$ v=\sqrt{2 g(\Delta h-f \Delta x)} $

(2)

式中：v为滑坡速率，Δh为滑坡的垂直高度，Δx为滑坡的水平距离，f为等效摩擦因数，g为重力加速度，一般取9.8 m/s²。

殷跃平^[22]、许强^[23]、郭长宝等^[1]分别采用式(2)对贵州关岭滑坡、四川茂县新磨村滑坡、理塘乱石包滑坡的滑动速率进行了计算。

式(2)为最后应用公式，未考虑初速度。其推导原理依据能量守恒定律，将全程的摩擦损失均匀化，定义为等效摩擦因数，即摩擦力和重力做功转化为动能的变化：

$ \Delta\left(\frac{1}{2} m v^{2}\right)=m g \Delta s \sin \beta-\Delta s m g f \cos \beta $

(3)

式中：Δs为滑坡斜距，β为西藏岸坡角。

如将Δssin β=Δh, Δscos β=Δx代入式(3)，可得

$ \frac{1}{g} \Delta\left(\frac{1}{2} v^{2}\right)=\Delta h-f \Delta x $

(4)

再经过变形，考虑初速度，可得滑体任意时刻所具有的瞬时速度：

$ v=\sqrt{v_{0}^{2}+2 g \Delta h(1-f \cot \beta)} $

(5)

如果滑坡发生在V型河谷，那么滑坡在对岸会有一定的冲起高度。如果考虑滑体在四川岸爬升过程，分析原理同上，那么滑坡冲起高度最顶点速度为0，则可利用滑坡冲起高度计算滑体爬升过程中任一高程点速度。

将滑坡末速度设为0，根据功能转化原理可得

$ 0-\frac{1}{2} m v^{2}=-m g \Delta h-m g f \Delta x $

(6)

$ \frac{1}{2} m v^{2}=m g \Delta h+m g f \Delta x $

(7)

$ v=\sqrt{2 g \Delta h\left(1+\frac{f}{\tan \alpha}\right)} $

(8)

如果已知滑体撞击四川岸时的速度，可以根据爬升顶点高程得出撞击点的高程：

$ \Delta h=\frac{\frac{v^{2}}{2 g}}{1+\frac{f}{\tan \alpha}} $

(9)

式中：Δh为滑体爬升高度；α为四川岸坡角，取26°；f为等效摩擦因数。由于谢德格尔假设滑体初速度为0时，f=h/x在白格滑坡不满足条件，选用经验公式^[24]lg f=alg V+b，a=-0.156 66，b=0.622 19，V为滑体体积，取1 408万m³，可得f=0.32。

已知白格滑坡主滑区和阻滑区滑体启动初速度为2.2 m/s，利用式(5)，Δh取为520 m(主滑区阻滑区滑体运动到海拔2 980 m)，β取为28°，f=0.32，可得

$ \begin{aligned} v=& \sqrt{2.2^{2}+2 \times 9.8 \times 520 \times\left(1-0.32 \times \cot 28^{\circ}\right)}=\\ & 63.8 \mathrm{~m} / \mathrm{s} \end{aligned} $

因此，得到滑体运动到剪出口位置速度为63.8 m/s，由于滑体从剪出口飞出时速度极大，飞跃时间极短，与四川岸坡碰撞时能量损失只是很少一部分，此处不进行考虑，将63.8 m/s记为滑体与四川岸坡撞击后的速度。代入式(9)可得

$ \Delta h=\frac{\frac{63.8^{2}}{2 \times 9.8}}{1+\frac{0.32}{\tan 26^{\circ}}}=125 \mathrm{~m} $

将滑体冲击海拔高度取为3 050 m，则大致得出主滑区和阻滑区滑体飞跃河岸撞击四川岸坡体海拔高度为2 925 m。

3.4 滑坡运动速度计算 3.4.1 主滑区及阻滑区滑坡运动速度计算

滑坡启动后，以一定初速度沿滑面下滑到海拔2 980 m时，从剪出口滑出，飞跃河岸撞击于四川岸海拔2 925 m处，滑坡碎屑体沿四川岸坡面向上爬高，而后回落堆积于河道之中，形成堰塞坝。通过现场调查，结合地形数据和遥感技术，获取滑坡的基本参数，可得滑坡冲击对岸(四川岸)痕迹弧上的5个特征点H₁、H₂、H₃、H₄、H₅(如图 8)，其冲击高度分别为130、170、176、172、139 m(相对金沙江河床最低海拔2 880 m)。因为已知滑体撞击四川岸海拔高度为2 925 m，则5个特征点H₁、H₂、H₃、H₄、H₅对应的滑体爬升高度分别为85、125、131、127、94 m。

Francis等^[25-26]针对岩质滑坡滑下并冲起的情况，提出了以下公式：

下降过程

$ \frac{1}{2} m v^{2}-m g H=-\frac{F_{1}}{100} m g H $

(10)

上升过程

$ m g h-\frac{1}{2} m v^{2}=-\frac{F_{2}}{100} m g H $

(11)

$ F_{1}+F_{2}=F_{\mathrm{T}} $

(12)

式中：m为滑体质量, kg; v为滑体速度, m/s; H为滑体下滑时下降的高度, m; h为滑体爬升时增加的高度, m; F₁、F₂分别为滑体下降和上升过程中因摩擦而损失的能量百分比; g为重力加速度，一般取9.8 m/s²。

为了估算滑动过程中的速度，将滑动路径分为N段，建立以下联立方程：

$ \begin{gathered} \left(\frac{1}{2} m v_{i}^{2}-\frac{1}{2} m v_{i-1}^{2}\right)+\left(m g h_{i}-m g h_{i-1}\right)= \\ -\frac{F_{i}}{100} m g H, i=1, \cdots, N \end{gathered} $

(13)

$ \frac{F_{i}}{W_{i}}=\frac{F_{\mathrm{T}}}{W_{\mathrm{T}}} $

(14)

$ W_{i}=-\mu m g \Delta l_{i} $

(15)

$ W_{\mathrm{T}}=-\mu m g L $

(16)

式中：H为滑坡原始高度，h_i为第i点高度，v_i为第i点速度，F_i为第i段的摩擦损失百分比，W_i为第i段损失的能量，F_T总摩擦损失百分比，W_T为总损失能量，μ为滑坡平均摩擦因数，Δl_i为滑坡体运动路径的水平距离，其余符号同上。

综合式(3)~(7)，可得滑坡路径中任一点速度计算公式：

$ v_{j}=\sqrt{2 g\left[H-h_{j}-(H-h) \frac{l_{j}}{L}\right]} $

(17)

其中 ${l_j} = \sum\limits_{i = 1}^j {\Delta {l_i}} , 1 \le j \le N$。

Oscar^[26]利用能量转化原理推导出的该速度计算公式，利用滑体产生的冲击高度计算滑体的速度，但是没有考虑初速度，必须在其基础上进行修改：

下降过程

$ \frac{1}{2} m v^{2}-\frac{1}{2} m v_{0}^{2}-m g H=-\frac{F_{1}}{100} m g H $

(18)

上升过程

$ m g h-\frac{1}{2} m v^{2}=-\frac{F_{2}}{100} m g H $

(19)

将滑动路径分为N段，仍有式(12)，取v_i-1=v₀，h_i-1=H, 按以上方法联立(18)、(19)、(14)~(16)，可得

$ v_{j}=\sqrt{v_{0}^{2}+2 g\left(H-h_{j}\right)-\frac{\left(2 g H+v_{0}^{2}-2 g h\right) l_{j}}{L}} $

(20)

其中 ${l_j} = \sum\limits_{i = 1}^j {\Delta {l_i}} , 1 \le j \le N$。

将5个特征点H₁、H₂、H₃、H₄、H₅对应的滑体爬升高度85、125、131、127、94 m，分别代入式(20)，在主滑区至对岸冲击区滑坡滑动方向上，选取7个滑坡高程点计算滑坡的滑动速度。由于对岸5个特征点高度不同，利用上述公式将滑坡体近似分为5部分在6个滑坡高程点进行计算。根据计算结果，白格滑坡从启动阶段开始，其滑动速度不断上升，在到达剪出口时达到顶值，为73.7 m/s，在四川岸撞击后其滑动速度为65.8 m/s，逆冲爬坡后滑动速度减小，直到停止。计算结果如表 1所示，滑坡主滑方向滑动速度特征如图 9所示。

Francis等^[25-26]结合滑坡运动的下降高度和上升高度，将摩擦损失的能量分散到运动全过程，在本研究计算滑体飞跃河岸过程中，同样计算了摩擦能量损失，由于滑体下降高差小于河床水平距离，速度有一定减少。实际上，滑体飞跃河岸过程中，会受到空气摩擦阻力而速度减小，更重要的是，滑体与四川岸撞击后，速度会进行衰减。在高速岩质滑坡动力学中，认为岩体材料性质在很大程度上决定了滑坡土石体与对岸陡壁撞击后能量的衰减系数，在数值上表现为恢复系数R。即用恢复系数来简化计算碰撞后岩体介质的速度^[27]：

$ R=\frac{v_{1}}{v} $

(21)

一般由坚硬完整基岩构成的岩坡表面具有较高的恢复系数，而软弱堆积物覆盖层覆盖的斜坡表面恢复系数较低。根据国内外主要斜坡岩体碰撞恢复系数的经验取值R=0.8^[27]，可计算得到白格滑坡主滑区阻滑区5部分滑体碰撞后速度分别为53.6、58.4、59.0、58.6、54.6 m/s，对应上述方法速度为54.3、54.6、65.8、65.0、56.5 m/s，计算结果相差不大，可认为合理。

3.4.2 牵引区滑坡运动速度计算

经现场调查，牵引区主要以强蚀变片麻岩为主，性状近似为土体，而且来自牵引区的滑坡碎屑堆积于堆积区靠右岸一侧，因此, 推断当锁固段脆性剪断后，主滑区阻滑区滑体以一定初速度先启动，牵引区后无初速度启动，主滑区阻滑区滑体高速飞跃河岸之后，逆冲到达冲击顶点，随后折返与牵引区岩土体相撞于金沙江上，堆积在金沙江右岸一侧。

已知牵引区体积为458万m³，按经验公式lg f= alg V+b可得f=0.38，β取30°，h=740 m，据式(5)可得

$ \begin{gathered} v=\sqrt{2 g h(1-f \cot \beta)}= \\ \sqrt{2 \times 9.8 \times 740 \times\left(1-0.38 \cot 30^{\circ}\right)}=70.2 \mathrm{~m} / \mathrm{s} \end{gathered} $

可知，牵引区土体在到达剪出口与返回的主滑区阻滑区滑体相撞时，速度极大，碰撞后的碎屑垂直向下运动，撞击河水，激起高速水砂射流与水雾，随后堆积于金沙江的右岸。

3.5 滑坡剧动时的能量计算

滑体在整体运动中，由于剪出口以上约790 m段内的滑面，相对上段较平缓，宽约500 m，构成了具有抗滑效应的锁固段，而滑体剧动正是由于克服了此段阻力，使之剪切破坏，相当于一次断裂破坏。在破坏的瞬间，滑坡体使长期积蓄的能量以应变弹性波的形式释放，从而产生震动效应。根据文献[28]，震级与破裂面积A(km²)之间关系为M=0.93log A+4.38；震级M与破裂长度L(km)经验公式为

$ M=(\log L+2.9) / 0.6 $

(22)

该滑坡破裂面长度为0.79 km，面积为0.395 km²，代入上述公式计算其震级分别为4.0和4.7。根据计算震级结果，利用Gutenberg等^[29-30]经验公式lg E=1.5M+4.8，计算滑体整体剧动时所释放的能量E至少为10^10.8J。因此，白格滑坡时引起的震动相当于4.0~4.7级表层基岩地震。

3.6 四川岸滑坡运动速度验证分析

根据以上分析可知，主滑区阻滑区滑体启动之后，是在高程2 980 m处剪出口滑出，在重力作用下滑体到达四川岸时，撞击点高程必然低于2 980 m，并且高于河床最低处高程2 880 m，若取滑体冲击高程为3 050 m，可得到滑体爬升高差最小为70 m，最大为170 m。利用上述公式得

$ v=\sqrt{2 g \Delta h\left(1+\frac{f}{\tan \alpha}\right)} $

经过计算可知，代入高差Δh₁=70 m时，得出v₁=47.7 m/s；代入高差Δh₁=170 m时，得出v₁=74.3 m/s。因此，主滑区阻滑区滑体在到达四川岸时，速度必然在47.7~74.3 m/s，这与本研究上述计算结果相符。

3.7 滑坡运动特征分析

“10·10”白格滑坡滑动运动过程具有如下特征：1)主滑区阻滑区滑体的主滑方向运动速度从2.2 m/s增加到65.8 m/s，变化极大，并且峰值速度较大，原因在于一方面，“10·10”白格滑坡相对高差巨大，从滑坡启动区到四川岸撞击点，高差达575 m，平均坡度为28°，地形陡峭，具有巨大的势能利于滑坡动能的释放；另一方面，白格滑坡的启动物质以岩体为主，并且以崩滑的形式启动，滑动速度快。2)根据滑坡冲击对岸的冲起高度，滑坡的滑动速度由中部向两端递减，原因在于滑坡的主滑方向滑动量大，运动距离最短，摩擦因数相对较小，其动能最大，冲起高度最大。3)牵引区岩土体由于依附于主滑区土体之上，主滑区和阻滑区启动之后，牵引区岩土体在重力作用下，无初速度下滑至金沙江之上。

4 结论

1) 白格滑坡锁固段脆性剪断，导致主滑区和锁固段滑体以2.2 m/s的初速度启动，高速飞跃河岸，逆冲到冲击顶点，随即折返与牵引区土体相撞于金沙江上。

2) 滑坡属高位高速岩质滑坡，由于本次滑坡动能巨大，在对岸留下明显的冲击痕迹，根据其对岸冲起高度，计算得到滑速在水平方向上的差异变化特性，即以主滑方向为中心，向两端延伸其滑速不断降低，其中主滑方向的速度最大，为73.7 m/s。牵引区无初速度启动，在到达剪出口时，速度为70.2 m/s。

3) 白格滑坡具有锁固段，剪切破坏瞬间，滑体整体剧动时所释放的能量E至少为10^10.8J，引起的震动相当于4.0~4.7级的表层基岩地震。