21世纪以来，中国进入地震高发期，在短短5 a间四川境内就相继发生了两场大地震 (Ms 8.0汶川大地震、Ms 7.0芦山地震)^[1-3].在这两场地震中，均触发了大量的崩塌、滑坡灾害，对灾区道路造成了严重损毁.据统计，汶川地震所引发的崩塌滑坡灾害造成公路受损近2.8×10⁴ km (21条高速公路、15条国省干线公路、2 795条农村公路)，同时也造成铁路受损 (成灌铁路、广岳铁路、德天铁路、宝成铁路绵阳至广元段等).山区公路、铁路作为生命线工程，在震后的抗震救灾中起着巨大的作用，而在震后的实地调研中发现，对线路损害最严重的是崩塌滑坡等次生山地灾害，少到数十方的崩塌落石、大到数百万方的滑坡体掩埋线路并阻断交通.在线路工程抢通过程中，极大地阻碍了救援队伍和工程机械进入灾区，增大了救援难度^[4].

5.12汶川地震触发的次生灾害结果表明，大量的高位崩塌滑坡是造成线路工程震害的主要威胁，如崔鹏等对Ms 8.0汶川地震后所触发的崩塌滑坡到岭脊与到河网的距离进行了归一化处理，并得到了地震触发的崩塌滑坡与河流、山脊位置分布关系.位于坡面距山脊归一化距离小于0.4的崩塌滑坡占总崩塌滑坡的65%，这是汶川地震作用下斜坡破坏的一个突出的特点，即崩塌滑坡发生部位一般都是位于山脊附近或山坡的坡肩，这与山脊处和坡肩处应力集中释放及对地震波的放大作用直接相关^[5].

汶川地震前，有关高位崩塌滑坡问题并未引起人们足够的重视，而在汶川地震后，就个体高位崩塌滑坡的形成机理、影响范围等研究已成为热点课题^[6]；对于铁路、公路而言，最关心的是高位崩滑坡体能否冲到路域内对线路工程造成危害.目前有关线路两侧边坡崩塌滑坡的判识已较成熟，但对来自于路域以外复杂坡体条件下的高位崩塌滑坡体能够威胁到线路工程的判识，还涉及到崩滑坡体的运动轨迹、运动距离的确定等问题，尚属甚少有人涉足的领域.

针对以上不足，本文提出了一种地震触发高位崩塌滑坡灾势预测方法，该方法在ArcGIS (地理信息系统软件) 技术支持下，按照地质条件及水文地质条件划分边坡单元，并建立区分威胁线路的高位边坡模型.该模型可以用于指导高烈度区线路选线以及既有线改建工程.

1 地震触发高位崩塌滑坡灾势预测的原理方法

地震触发高位崩塌滑坡灾势预测方法的原理是基于卫星遥感数据，应用GIS的空间分析技术与地震触发崩塌滑坡统计方法相结合，对预建线路或既有线改建地区划分边坡单元，然后建立高位边坡单元判识模型，该模型根据边坡能否同时满足能量条件及路径条件来判断威胁线路安全的高位边坡.最后基于综合指标法快速生成区域地震触发高位崩塌滑坡灾害的灾势区划图，其判别程式如图 1所示.

1.1 区域边坡单元化

目前划分边坡单元的方式主要分为规则网格单元和复杂边坡单元两种.规则网格单元的计算和算法相对简单，但与实际的边坡破坏单元及力学机制无任何联系，与地质、地形及水文也毫不相干.边坡单元定义为一块与邻近区域具有明显不同的地形特征的区域，而地形的形成本身反映了地质条件及水文地质条件的长期作用效果^[7-9].如图 2边坡划分示意图所示，通过ArcGIS的水文分析功能找出数字地形图中的山脊线与山谷线，并通过山谷线与山脊线将集水流域划分出边坡A与边坡B.边坡单元的大小选定取决于研究区域内历史滑坡的平均尺寸 (如果缺少历史滑坡的记录可参考临近地区历史滑坡资料来确定边坡单元的大小平均尺寸)，即对于同一区域内具有相同或类似的地形、地质及水文条件的研究对象，可以认为其将来可能产生的滑坡破坏与过去的滑坡破坏在规模上是一致的.在实际应用时，需要进行多个阈值的对比、试算确定合适的阈值水平^[9].

水文特征与区域地表高程的起伏状况有关，降落在山体表面上的雨水，除去下渗部分，其余在地表径流的水体则遵循一定的规律在重力作用下从坡顶高处向坡底低洼处流动并汇集.使用ArcGIS软件的水文分析工具能够根据水流流动规律，完成对区域内的水文分析.利用水文分析模块分别对DEM (数字高程模型) 进行无洼地DEM提取、水流方向提取、水流汇集、河网生成、集水流域分析，得到山脊线.而山谷线的提取则是通过将原来的山谷变为山脊，然后采用相同的方法对反向DEM进行无洼地DEM提取、水流方向提取、水流汇集、河网生成和积水流域分析等操作就可以得到山谷线.经栅格数据转换成面域，在此运算中会产生一些不正确的多边形边坡单元和面积很小或不协调的多边形边坡单元，通过使用GIS的融合 (dissolve) 和统合 (integrate) 功能消除不合理的多边形边坡单元.

1.2 地震触发高位崩塌滑坡模型判别方法

高位崩塌滑坡具有滑动距离远、破坏强度大等特点，一直都是众多学者研究的热点.目前，国际上通常采用滑坡体质心的最大高度H^*与质心的最大水平距离L^*(或滑坡顶部的最大高度H_max和最大水平距离L_max) 之比 (μ=H^*/L^*) 作为描述滑坡运动性的特征参数 (见图 3高位滑坡运动特征)，这一特征参数在物理意义上相当于整个滑道上的平均摩擦系数^[10]. μ值只与动能和动摩擦系数有关，能够较好地表征滑坡的运动性. μ值越大，平均摩擦系数越大，高位滑坡滑动时越困难；μ值越小，平均摩擦系数越小，高位滑坡滑动时越容易. Legros等对世界上203处在各种环境下产生的滑坡、运动距离、下落高度、滑坡面积以及滑坡运动特征值μ进行统计分析，结果表明发生在陆地上的滑坡μ值大约在0.1~0.6之间，发生在火山口的滑坡μ值约为0.1，发生在海底滑坡μ值约为0.05~0.1^[11].许强等为了研究汶川地震区内大型高位滑坡的运动特征，根据遥感影像从112处大型滑坡中筛选出由滑动转化为碎屑流并做远程运动的16处高位滑坡作为重点分析研究对象^[12].

根据这16处高位滑坡体可得到了L^*/H^*与滑坡面积S的对数关系，见图 4，对其进行线性拟合，拟合公式为

$\begin{align} &{{L}^{*}}/{{H}^{*}}=0.5684\ln s-4.6235\Rightarrow \\ &\ \ \ \ \ \ \ \ \ \ \ \ \ {{H}^{*}}=\frac{{{L}^{*}}}{0.5684\ln s-4.6235}, \\ &\ \ \ \ \ \ \ \ \ \ \ \ {{R}^{2}}=0.6427. \\ \end{align}$

(1)

由式 (1) 可知，当面积为S的高位边坡在地震作用下发生滑动的最远距离为L^*时，其质心点临界高度约为H^*.故可以通过此公式来判别高位边坡在地震作用下发生崩塌滑坡是否会对线路工程造成威胁.当高位边坡的质心点高程H < H^*时，则不会威胁到线路工程，如图 5所示，尽管在地震作用下发生崩塌但是由于其质心点高程不够，致使其无法达到坡脚并威胁线路工程；而当高位边坡的质心点高程H > H^*时，则有可能威胁线路工程.需要说明的是，式 (1) 既可以采用滑坡体的表面面积S也可以采用总体体积V，两者都是灾害规模的表征方式并且之间存在一定相关关系.但是较难准确预测地震作用下滑坡可能的滑动深度，而表面面积S却可以准确得到.因此，本文在对地震触发高位崩塌滑坡的灾势预测时只考虑表面面积S.该方法能够满足在选线设计中指导线路工程的方案比选精度.

结合以上分析及实地调查发现，在地震作用下，高位边坡须同时满足能量和路径两种必要条件，才能判定其是否会对线路造成危害.所谓满足能量条件是指使高位崩坡体滑动到路域以内时所需具有足够的重力势能.而所谓满足路径条件是指其滑动时的运动路径.段书苏等对芦山地震实地调查发现，在Ⅶ度及以上地震烈度区确定的1 754个地震触发崩塌滑坡点中，崩塌滑坡个数的优势方向与自然坡体的优势坡向相吻合^[13-14]，因此本文将高位边坡的主坡向做为滑动的优势方向.但是高位边坡并非都会沿着边坡主坡方向一滑到底，如图 6所示，青川马鹿乡老鹰岩危岩首先沿边坡坡向滑动到冲沟1中，然后由冲沟1滑向村庄.经过这一路径不仅增长了其滑动路径的长度，同时在冲向冲沟1转向时也消耗了部分动能，因此路径是必须要考虑的要素.

基于以上分析提出了高位崩塌滑坡是否会威胁到线路工程的判别模型，判别方法如下.

满足能量条件判别方法：当面积为S的一处高位边坡质心点到线路距离L^*已知时，就可以根据式 (1) 判断其所对应的临界质心高度H^*.若此高位边坡的质心点高程H > H^*时，则可以认为该滑坡体满足能量条件，故可以判定其为可能威胁到线路工程的高位边坡.

满足路径条件判别方法：针对如上满足能量条件的高位边坡，根据其主倾向优势方向及冲沟位置，确定滑坡路径的水平投影长度，并重新将其作为L₁^*，根据式 (1) 判断调整后的临界质心高度H₁^*，如果此时高位边坡的质心点高程H>H₁^*，则认为该滑坡体同时满足能量条件和路径条件，故可以判定其为威胁线路工程的高位边坡.

在此基础上结合综合指标法预判高位边坡在地震作用下发生滑坡的可能性，完成地震触发高位崩塌滑坡的灾势预测.

1.3 基于综合指标法预测研究区地震崩塌滑坡

基于威胁线路工程边坡的区划图，采用综合指标法，对地震触发的崩塌滑坡进行预测.该方法是在对中国已有地震资料研究分析基础上，将地层岩性构造、坡高、坡角、地震、降雨作为影响地震崩滑的主要因子，并对其进行分级赋值，适用于我国西部地区地震崩塌滑坡的预测^[16]，计算式为

$H=\left( {{S}_{i}}\times {{S}_{a}}\times {{S}_{h}} \right)\times \left( {{T}_{s}}+{{T}_{p}} \right).$

(2)

式中：S_i为岩性条件系数；S_a为坡角影响系数；S_k为坡高影响系数；T_s为地震强度系数；T_p为降水强度系数；H为边坡危害程度.当H < 264时，所有边坡都是稳定的；当H>396时，边坡均失稳；H值处于二者之间时，边坡既有可能稳定又可能失稳.

式 (2) 表明了地震崩滑的产生是外部因素 (地震、降雨) 相叠加作用于内部因素 (坡高、岩性、坡角) 的结果.最后，计算研究区内稳定边坡和失稳边坡的综合指标值，进而确定地震崩滑发生临界值.

2 案例分析

选取西藏地区通麦-105道班段作为研究区域，该区域是川藏公路的重要组成部分，是一既有改建工程.区内地震活动频繁，新构造运动强烈，地质构造错综复杂，岩体破碎，松散失稳，降水丰沛，海洋性冰川广泛分布，河流水系发育，河谷、沟谷侵蚀下切强烈，使该段道路成为川藏公路全线地质病害病种齐全，灾害集中、频繁，规模最大及损失最为严重的路段，是典型的地震触发高位崩塌滑坡研究区域.为此，购买了关于该区域的资源3号，高分1号卫星影像数据，影像精度为1:25 000.

2.1 研究区内的边坡划分

该研究区通过水文分析方法共判别出804处边坡.然后，通过C#与ArcGIS软件结合，二次开发编程统计各边坡的平均倾向.

2.2 研究区内威胁线路工程的高位边坡判识 2.2.1 判断满足能量条件的高位边坡

假定边坡在地震作用下发生崩塌滑坡时，其滑坡体的质心与形心重合.首先，求取各边坡单元的面积S及形心点；其次，利用ArcGIS的Near工具，判断各边坡平面形心点到线路的垂直距离，将其作为最远滑动距离L^*；然后，根据式 (1) 求取每一边坡形心点所对应的临界高度H^*；再次，利用ArcGIS的Extract values to points工具求取各个边坡形心点的高差H，并判别H > H^*的边坡.以通麦-105道班地区为例，判别结果如图 7所示，图中威胁线路的高位边坡是满足能量条件的高位边坡.

2.2.2 判断同时满足能量和路径条件的高位边坡

基于图 7满足能量条件的高位边坡判断满足路径条件的高位边坡.根据高位边坡的主坡方向及冲沟位置，判断其滑动时可能的滑动路径，将满足能量条件的高位边坡分为3类：第1类高位边坡是指滑坡滑动路径既不存在缓坡也不经过支沟；第2类高位边坡是指滑坡滑动路径不经过支沟但存在较长缓坡，如图 8所示；第3类高位边坡是指滑坡滑动路径需经过支沟，然后由冲沟滑向线路工程.对于第1类高位边坡可认为满足能量条件即可，而对于第2、3类高位边坡还需进行单独判断.

对于第2类高位边坡又可分为先陡后缓型和先缓后陡型高位边坡.先陡后缓地形对高位滑坡的滑动不构成制约因素，滑坡受地震作用启动后就可利用陡坡迅速将重力势能转化成动能，而先缓后陡地形则对高位边坡滑动构成制约因素，滑坡受地震作用启动后，需要克服较长的缓坡地段，无法迅速将重力势能转化成动能，此时需要考虑先缓后陡型高位边坡是否能够滑过缓坡区.由于已经经过能量条件判断，故只需要再判断第2类的先缓后陡型高位滑坡体是否能滑过缓坡区即可.首先确定滑坡体质心到缓坡终点水平投影距离L₁^*，然后根据式 (1) 判断两点间的临界质心高度H₁^*，如果此时高位边坡的质心点高程H>H₁^*，则认为该滑坡体能够通过缓坡地段，是威胁线路的高位边坡. 图 9为第2类高位边坡平面图，该图中的高位边坡属于第2类高位边坡，其边坡面积S=69 266 m²，L₁^*=900.0 m，H₁^*=525.8 m，H=500.0 m，H < H₁^*，故此高位边坡满足能量条件，不满足路径条件.

对于第3类高位边坡，首先需要根据滑动路径确定滑坡路径的水平投影长度，重新将其作为L₁^*，然后根据式 (1) 判断调整后的临界质心高度H₁^*，如果此时高位边坡的质心点高程H>H₁^*，则认为该滑坡体同时满足能量和路径条件，可以判定其为威胁线路工程的高位边坡. 图 10为第3类高位边坡滑动路径平面图，该图中的高位边坡属于第3类高位边坡，需要确定其滑动路径.滑动路径沿该边坡的主坡向优势方向滑动到冲沟后，沿冲沟滑向线路，如图 10中黑线所示.其边坡面积S=192 184.5 m²，L₁^*=3 032.0 m，H₁^*=1 322.9 m，H=1 275.0 m，H < H₁^*，故该高位边坡满足能量条件，不满足路径条件.此外，高位滑坡沿主坡优势方向滑向冲沟后会造成一定的动能损失，动能损失大小与主坡优势方向和冲沟间的夹角大小有关.本文未考虑动能损失，对预测结果更趋于安全.

综合第1、2、3类高位边坡，得到了威胁线路工程的高位边坡图，如图 11所示.该图中威胁线路工程的高位边坡比满足能量条件的高位边坡少了32处，进一步缩小了范围.

2.3 研究区地震触发高位崩塌滑坡判别

1) 区域气候特征.研究区域气候温和湿润，年平均气温为8~13 ℃，年降水为1 000~1 400 mm；一年干湿两季气候，4月份至10月份都属雨季，长达7个月，11月份至次年3月份为旱季，仅5个月.雨季降水占了全年降水总量的85.5%，旱季仅占14.5%；一年中雨季虽长，雨日也多，降水比率很高，但绝大多数都是小到中雨，大雨很少，没有日降水大于50.0 mm的暴雨发生.

2) 区域边坡平均坡度及边坡高度.基于C#与ArcGIS软件二次开发编程统计划分好的每一边坡单元的平均坡度及边坡单元的相对高度，其中边坡平均坡度如图 12所示.

3) 区域地层岩性及近场区地震构造.根据已有地质资料可知，滑坡上覆底层主要为第四系全新统松散堆积体 (第四系全新统崩坡积层Q₄^c+dl、第四系全新统滑坡堆积层Q₄^del、第四系上更新统冲洪积层Q₃^al+pl第四系全新统泥石流堆积层Q₄^sef、第四系全新统冲洪积层Q₄^2al+pl).下伏基岩主要为前震旦系冈底斯岩群 (AnZgd) 及前震旦系南迦巴瓦岩群阿尼桥片岩组 (AnZa).区域内地震活动水平较高，在空间上分布极不均匀，大震多位于断裂带上.估计未来百年在墨脱8.6级地震附近地区仍存在发生7级地震的可能性，但发生7.5级以上地震的可能性不大.近场区内历史地震较多，部分地震震级较小，但距离场址较近，对场区造成一定的影响.现代地震活动频度高，但震级较小，对场区的影响强度较弱.对工程场区影响烈度最大者应当为1950年墨脱8.6级地震^[17-18]，距场址的最小距离仅为125.4 km，其近场区地震构造.

4) 地震崩滑预测结果.基于综合指标法对研究区域内的边坡进行地震触发高位崩塌滑坡灾势区划，其判别情况如图 13所示.

3 结论

1) 提出一种线路工程地震触发高位崩塌滑坡灾势预测方法.该方法是根据Ms 8.0汶川地震高位滑坡的特征建立起来的，可以在路域以外的复杂坡体上，判识由地震触发的高位崩塌滑坡是否会滑动到路域内威胁线路工程.

2) 汶川实震表明，大量的高位崩塌滑坡是造成线路工程震害的主要威胁，因此当线路工程通过高烈度区时，需考虑防范来自线路用地界以外的高位崩塌滑坡.所提出的地震触发高位崩塌滑坡灾势预测方法可以为线路通过方案与绕避方案的比选提供决策依据.鉴于样本数仅有16处，研究工作仍属第一阶段，随着样本数据的不断增加，地震作用下威胁线路工程的高位崩塌滑坡灾势预测的精度将会进一步增加.