2014年鲁甸MS6．5地震震前与同震滑坡空间分布规律对比分析

田颖颖　许　冲　徐锡伟　吴赛儿　陈　剑

1）中国地震局地质研究所，活动构造与火山重点实验室，北京　100029

2）中国地质大学（北京），工程技术学院，北京　100083

2014年鲁甸MS6．5地震震前与同震滑坡空间分布规律对比分析

田颖颖1，2）许冲1）*徐锡伟1）吴赛儿2）陈剑2）

1）中国地震局地质研究所，活动构造与火山重点实验室，北京100029

2）中国地质大学（北京），工程技术学院，北京100083

2014年8月3日，云南省昭通市鲁甸县发生了MS6.5地震，地震诱发了大量滑坡。文中以牛栏江沿线鲁甸县、巧家县和会泽县交界处面积为44.13km2的区域为研究区，开展地震震前与同震滑坡的空间分布规律对比分析。根据震前Google Earth高分辨率影像与震后0.2m分辨率的超高分辨率航片数据，分别建立了震前滑坡与同震滑坡数据库。结果表明，研究区内震前有284处滑坡，本次地震触发1 053处滑坡。借助10m×10m分辨率的数字高程模型（DEM）数据，基于GIS平台提取研究区的高程、坡度、坡向、曲率、岩性、烈度、河流共7个主要因子，并利用滑坡的面积百分比（Landslide Areas Percentage，LAP）和点密度（Landslide Number Density，LND）对比分析震前与同震滑坡的空间分布规律。结果表明，震前与同震滑坡的易发高程区间分别为＜1 200m与1 200～1 300m。坡度越大越容易发生滑坡，其中坡度＜10°的区域由于距离河流很近，也为滑坡易发区。震前与同震滑坡发育的优势坡向都是近S向。当斜坡为凹坡时（曲率值为负值），滑坡易发性较高。地震烈度越大，越易发生同震滑坡。灰岩夹白云质灰岩分布区很容易发生滑坡，玄武岩和火山角砾岩分布区在地震力的作用下边坡的稳定性也大大降低。震前、同震滑坡的发生与到河流的距离大致呈现正相关性。震前滑坡LAP的峰值大多数都与震前已经存在的大型滑坡有密切的对应关系。

鲁甸地震震前滑坡同震滑坡空间分布滑坡密度

0　引言

2014年8月3日16时30分，云南省昭通市鲁甸县（27.1°N，103.3°E）发生6.5级地震。地震造成云南省昭通市鲁甸县、巧家县、昭阳区、永善县、曲靖市会泽县共55个乡镇遭受不同程度破坏，617人死亡、112人失踪、3 143人受伤，22.97万人紧急转移安置。

表1　震源机制解参数Table 1　Focal mechanism solution of the Ludian earthquake

震后利用地震波形记录，通过P波初动方向快速求得震源机制解（表1）。结果显示2个节面的走向分别为NE向和NW向。由于此次震中位于NE向昭通－鲁甸断裂附近，初期做出了地震为NE向破裂的误判。后期，通过余震空间优势分布方位、强地震动记录、现场地震灾性讨论（程佳等，2014；徐锡伟等，2014a，b；赵静等，2014）和地震地质灾害（李西等，2014；许冲等，2014a）等方面。地震触发滑坡往往数量多、分布密集、致灾严重，如许冲等（2013a）利用目视解译的方法针对2008年汶川MS8.0地震建立了地震滑坡编目数据库，共解译出197 481处滑坡，总面积达1 160km2。统计结果表明，汶川地震滑坡直接导致2万余人死亡，上千亿元经济损失（殷跃平，2009）。对地震诱发的滑坡进行空间分析，有利于了解未来再次发生地震时比较容易发生滑坡的区域，从而减少人员伤亡和经济损失（Keefer，2000；Dai et al．，2011）。Keefer（1984）对全球40多次历史地震滑坡进行了编目，并讨论了滑坡的空间分布和地震参数的关系。之后，Keefer（2000）用滑坡密度分析了1989年Loma Prieta地震滑坡的分布与到震中的距离、坡度和岩石类型的关系。Wang等（2007）利用滑坡面积百分比与滑坡密度分析2004年Chuetsu地震诱发滑坡与研究区内地质、地形地貌、地震等参数的关系。Qi等（2010）在汶川地震后也利用滑坡密度对2008年汶川地震触发的滑坡进行了研究，对活动构造、岩性以及地形等因素进行了空间分析。Dai等（2011）利用面积密度和点密度对汶川地震诱发的滑坡进行了研究，对震源、地形和岩性等参数进行了分析。类似研究的实例还有很多，如针对1994年美国北岭地震（Harp et al．，1995）、1999年台湾集集地震（Liao et al．，2000）、2005年克什米尔地震（Sato et al．，2007）、2010年海地太子港地震（Xu et al．，2014c）、2010年青海玉树地震（许冲等，2011；许冲等，2012；Xu et al．，2013a，2014a）、2013年四川芦山地震（许冲等，2013b，c；Cui et al．，2014；Liu et al．，2014；Wang，2014；Xu et al．，2014b；周庆等，2014）、2013年甘肃岷县漳县地震（田颖颖等，2013；许冲等，2013d，e；Xu et al．，2014d）等震例所进行的地震滑坡空间分布规律研究。然而，以上研究主要针对同震滑坡进行空间分布与地形、地质和地震等因素的分析，而关于震区同震滑坡与震前滑坡分布规律差异对比的研究却非常少。

2014年8月3日鲁甸地震发生后，许冲等（2014a）根据地震前后的高分辨率（2m）的卫星影像，开展了覆盖整个鲁甸震区的地震滑坡解译。结果表明，鲁甸地震触发了至少1 024处面积＞100m2的滑坡，这些滑坡分布在一个面积＞250km2的区域内，滑坡面积约为5.19km2，体积约为2.2×107m3。尽管如此，通过与覆盖有限震区的震后超高分辨率（0.2m）航片比较，一些在0.2m分辨率的航片上能够清楚识别的滑坡在分辨率为2m的SJ9A或TH01－02卫星影像上却难以识别。因此，利用航片解译的滑坡要比基于卫片解译得到的滑坡详细得多。但是，考虑到震后的航片仅仅覆盖了部分地震区，因此本文选择航片覆盖良好的牛栏江一段流域作为研究区，结合Google Earth上震前影像与震后超高分辨率航片，借助GIS平台，基于地震、地质和地貌因子，从对比地震前后滑坡的空间分布规律的角度出发，开展鲁甸地震震前滑坡与同震滑坡的分布规律对比，分析了滑坡在高程、坡度、坡向、曲率、岩性、烈度和河流7个因素影响下的易发参数区间，为震区滑坡防灾减灾提供科学参考。

1　研究区与区内滑坡

鲁甸地震发生在云南省鲁甸县，牛栏江北侧。此次地震震中位于龙头山镇，距离鲁甸县城约25km，距震中30km范围内平均海拔约为2 500m。该地区属于低纬度山地季风气候，干湿季分明，湿季主要集中在5—10月，雨量集中，年降雨量达1 000mm（唐立梅，2007）。本文所选的研究区位于鲁甸县、巧家县和会泽县交界处，沿牛栏江对称分布。研究区面积44.13km2，高程分布范围1 000～2 000m，坡度范围20°～50°，牛栏江是研究区的主要河流，此外其支流沙坝河也是研究区内的一条主要河流。牛栏江和沙坝河的强烈下蚀形成了该区V型高山峡谷地貌（胡金等，2008）。研究区主要位于Ⅶ、Ⅷ和Ⅸ烈度区（http：∥www．cea．gov．cn／publish／dizhenj／464／478／20140807085249557322083／index．html）。根据区域地震构造、余震分布、主震震源机制解、烈度分布长轴方位和滑坡空间分布样式等资料，确定此次地震的发震断裂为震区内NW向的包谷垴－小河断裂（常祖峰等，2014；房立华等，2014；李西等，2014；徐锡伟等，2014b；张盛峰等，2014）。鲁甸地震震源特征研究结果（赵旭等，2014）表明，此次地震为一次左旋走滑型兼有少量正断倾滑分量的事件。起始破裂点开始以震源为中心向NW和SE两侧同时扩展，2s后破裂主要向SE方向扩展（赵旭等，2014）。

图1　鲁甸MS6.5地震区域构造与研究区分布图Fig．1　Regional tectonic setting of the Ludian MS6.5 earthquake affected area and the study area．

基于GIS和遥感平台，分别对Google Earth上的震前高分辨率卫星影像（2014年1月30日）和震后0.2m高分辨率航片（2014年8月4日）开展滑坡解译工作。滑坡识别方法与准则主要参考许冲等（2014b）提出的地震滑坡编录原则。结果表明，解译得到的震前284个滑坡主要沿河流和道路分布，最小的滑坡面积约为48m2（图2a）。鲁甸地震触发滑坡1 053处，是震前滑坡数量的3.7倍。同震滑坡最小面积仅为2m2左右（图2b），表明了超高分辨率（0.2m）航片强大的地震滑坡辨识能力。

选取高程、坡度、坡向、曲率、岩性、烈度和河流7个控制滑坡发生的因素，在GIS平台下对这些因素进行分类，基于由面数据提取出的滑坡点数据，并以滑坡点密度（LND）和面积百分比（LAP）2个指标对地震前和同震滑坡分布的空间分布规律进行分析，各个因素的分类图见图4。其中，LND参考Keefer（2000）提出的滑坡密集程度LC来定义，LAP主要参考Xu等（2014e）提到的滑坡的面积百分比LAP定义，即：滑坡的点密度（LND）＝对应分类滑坡的个数／分类区间总面积；滑坡的面积百分比（LAP）＝对应分类滑坡的面积／分类区间的面积× 100%。

图2　研究区的滑坡分布图Fig．2　The spatial distribution of the landslides in the study area．

图3　震前、同震典型滑坡图示Fig．3　The typical pre-earthquake and coseismic landslides．

2　震前、同震滑坡空间分布规律研究

2.1高程

图4　研究区各因素的分类分布图Fig．4　The classification and distribution of the influence factors of landslide in the study area．

表1　研究区地层和岩性分类Table 1　The stratum and lithological classification of the study area

地震滑坡往往在某个高程范围内较易发生。研究区内高程范围是1 001．74～1 953．85m。基于ArcGIS与10m×10m分辨率的DEM，将研究区分为9类，高程范围分别为＜1 100m，1 100～1 200m，1 200～1 300m，1 300～1 400m，1 400～1 500m，1 500～1 600m，1 600～1 700m，1 700～1 800m，＞1 800m。震前滑坡的LAP在1 300～1 400m和1 100～1 200m范围内较大，分别占该分类区域总面积的17.93%和11.37%；＞1 400m范围内滑坡的LAP很小。＜1 100m和1 100～1 200m分类区域内滑坡的LND较大，分别为22.89个／km2和26.31个／km2。其中在1 300～1 400m范围内，滑坡的LAP最大而滑坡的LND却较小。同震滑坡的LAP大致先增后减，在1 200～1 300m分区内滑坡达到最大，约为11.72%；LND也先增后减，在1 100～1 200m高程范围内达到最大，约为56.50个／km2。在各高程范围内，震前1 100～1 200m的范围内LAP和LND同时达到峰值，这说明，虽然该区域分级面积不是很大，但较容易发生滑坡。震前LAP的最大值出现在1 300～1 400m范围内，但此区域内滑坡的LND值较小，这说明该区域内发育大型滑坡，通过对比震前影像，震前红石岩大滑坡即位于该分类区域内。同震滑坡LAP最大的区域为1 200～1 300m，该区域内LND值也较大，这说明在地震力的作用下，1 200～1 300m范围内岩土体很容易失稳发生大量小规模的滑坡。震前、同震滑坡的LND在1 100～1 200m范围内均达到最大，这说明了无论震前、同震，该区域内都比较容易发生滑坡，这是由于该分级不仅距离河流较近，还有道路分布，人类对其地形的改造等作用也很明显。＜1 100m的范围主要位于牛栏江西侧，此处河流的切割作用比较明显，河谷狭窄，分级面积较小，所以滑坡数量和面积也较小。＞1 400m的范围内滑坡的LND和LAP都较小，这是由于这些区域位于河谷较宽阔的部分（Qi et al．，2010），切割作用较弱，地形较平缓，不容易发生滑坡。

图5　滑坡分布与高程关系图Fig．5　The relationship between landslide distribution and elevation．

2.2坡度

坡度也是影响滑坡分布的因素之一。一般坡度越大，越易发生滑坡。通过对研究区的DEM进行坡度提取，得到研究区滑坡坡度范围0°～70.26°，本文以10°为分类间隔对坡度进行重分类：＜10°，10°～20°，20°～30°，30°～40°，40°～50°，＞50°，统计每个分类滑坡的面积和数量，制作滑坡LND和LAP图（图6）。震前滑坡在30°～40°以及50°以上坡度范围内LAP较大，分别为10.82%和11.33%，30°以下的范围内滑坡的LAP较小。滑坡的LND先减小后增加，在＞10°的范围内随着坡度的增大而增大，最大值为10.20个／km2，最小值为10°～20°区间的3.55个／km2。同震滑坡LAP随着坡度的增加呈增大趋势，当坡度＞50°时，LAP突变至最大，为38.35%；滑坡的LND先减小后增加，最后平稳发展，在＜10°的范围内最大，为28.36个／km2；20°～30°范围内最小，为19.96个／km2。无论震前、震后，＞50°范围内滑坡的LAP都达到了极值，即该分级区域内较易发生滑坡，这是由于＞50°坡度分类区间坡度较陡，而坡度越陡，地震波的地形放大效应越显著（刘洪兵等，1999），加之其分类面积很小，使得该分区内滑坡的LAP和LND都较大。震前30°～40°分级范围的LAP值也达到了峰值，这是由于该分级内包含了几处大型滑坡，如震前红石岩滑坡。无论震前、同震条件下＜10°的范围内，由于接近于岸边，岩土体受河流的冲刷侵蚀和孔隙水压力的作用，也很容易发育大量滑坡，所以LND线在此范围有所上升。同震滑坡的LND与坡度的相关性不是很大，这可能是由于研究区较小。

图6　滑坡分布与坡度关系图Fig．6　The relationship between landslide distribution and slope．

2.3坡向

基于GIS和DEM，提取研究区的坡向分布并把其重分类为平坦、N、NE、E、SE、S、SW、W和NW 9类，经过统计发现：当坡向平坦时，不发育滑坡。由此制得的滑坡随坡向的LND和LAP分布图见图7。由震前滑坡与坡向分布图可以很明确地看到，正N向滑坡的LAP最大，为21.41%；其次为NE向和ES向，其LAP值为5.60%和4.67%。而LND在SE向最大，为11.65个／km2。对于同震滑坡，LAP在正S向达到最大，约为19.65%；而LND在S、SE和SW向均较大，其值分别为37.78、32.41和36.02个／km2。不同坡向影响下滑坡的易发性和规模也是不同的（Xu et al．，2014c）。滑坡的易发性与坡向的关系，能在一定程度上反映地震的发震断裂的位移方位（Xu et al．，2014e）。通过对比震前、震后坡向统计图，发现不管震前还是在地震力作用下，接近S方向的斜坡都比较容易发生大型滑坡，这一方面可能是由于牛栏江在研究区内共有4处急转弯，其北岸正好为河流凹岸，受河流冲刷侵蚀严重，位于北岸的沙坝河也由N向S汇入牛栏江，使牛栏江北岸坡向接近S向的岩土体很容易发生滑坡；另一方面同震滑坡在S、SW和SE向很发育，也可能与此次地震发震断裂的活动特征有一定的联系。震前某些大型滑坡，如震前红石岩滑坡，倾向大致为N，这也使震前N方向LAP达到了极值。另外，震中位于研究区的NNW向，在地震产生的纵波传播过程中对近S向斜坡产生背离斜坡的拉应力，而对近N向的斜坡产生挤压应力（Xu et al．，2013a，b），这也是导致震后滑坡在正S向数量多且规模较大的原因之一。

2.4曲率

斜坡曲率也是影响滑坡分布的一个地形因素。曲率为负值表示凹坡，曲率为正值表示凸坡，曲率为0或接近0为平直坡。对DEM进行曲率提取并重分类为6类：＜－2，－2～－1，－1～0，0～1，1～2，＞2，分别统计各个分类的LND和LAP。震前滑坡曲率在－1～0的滑坡LAP最大，为7.97%。－2～－1滑坡的LND最大，为1.57个／km2。同震滑坡在曲率＜－2的LAP和LND都达到了最大，分别为78.14%和42.63个／km2。分析结果表明，震前滑坡最易发生在凹坡处（曲率值为负值），即被河流冲刷严重的一侧。－1～0的峰值也是由于震前红石岩滑坡的存在造成的。祁生文等（2004）认为在地震动力作用下，凹坡比凸坡和直线型斜坡更容易发生滑坡、崩塌等现象，同震滑坡统计图更明显地证明了这一点。即在地震力的作用下，凹坡的稳定性更差。

图7　滑坡分布与坡向关系图Fig．7　The relationship between landslide distribution and aspect．

图8　滑坡分布与曲率关系图Fig．8　The relationship between landslide distribution and curvature．

2.5岩性

岩性是影响滑坡发育和分布的地质因素。依据表1对研究区岩性进行重分类，进而统计每个类别内滑坡的LAP和LND。震前第5类（泥质白云岩与泥页岩间层）区内滑坡的LAP最大，为105.75%；而其LND却较小。在第2类（灰岩夹白云质灰岩）区的LND较大，10.98个／km2；但其LAP很小。另外，第1类（玄武岩、火山角砾岩）、第4类（白云岩、灰岩夹砂、页岩）区的LAP和LND都很小。对比同震滑坡分布的统计图，第4类区内（白云岩、灰岩夹砂、页岩）的LAP最大，为22.12%；但其LND值也较小。滑坡的LND在第1，2类区内（玄武岩、火山角砾岩和灰岩夹白云质灰岩）较大，分别为24.76个／km2和39.85个／km2。在第5类区内（泥质白云岩）滑坡的LAP和LND都较小。对比震前、同震的统计图，震前、同震滑坡LAP最大值分别位于第5类（泥质白云岩、页岩）和第4类（白云岩、灰岩夹砂页岩）区域，而其LND则正好相反，这是由于震前的红石岩滑坡位于泥质白云岩和页岩区；在地震力作用下，震前红石岩滑坡上仅发育了几处浅层小型滑坡，而在其对岸，白云岩、灰岩夹砂页岩区形成了同震的红石岩滑坡，该滑坡体阻塞了牛栏江，形成了堰塞湖。震前、同震的LND最大值都位于第2分类（灰岩夹白云质灰岩区），说明该类区内发育了大量的小型滑坡；经分析发现该类区主要沿牛栏江分布，不仅面积最大，其受河流的影响也是很显著的。其中，玄武岩、火山角砾岩分布区由于面积很大，且地势较高、地形平缓，震前很少发育滑坡；而震后在沙坝河汇入牛栏江附近的河流拐弯处发育了大量的滑坡，使得震后该类区内滑坡的数量和面积有所增加。

图9　滑坡分布与岩性关系图Fig．9　The relationship between landslide distribution and lithology．

2.6烈度

烈度是影响滑坡分布的地震因素。烈度越大，越容易发育滑坡。由图1可以看出，研究区处于Ⅶ、Ⅷ、Ⅸ3个烈度等级内，烈度为Ⅶ的区域所占的面积很小。通过ArcGIS的空间分析功能，统计这3个烈度区内滑坡的LAP和LND。在无地震力作用下，烈度为Ⅸ的区域内的滑坡LAP最大为6.19%，在Ⅷ烈度区内最小为3.40%；LND在Ⅷ烈度区间内最大，为7.93个／km2，烈度为Ⅸ的区域滑坡LND最小，约为2.70个／km2。同震滑坡的LAP和LND在Ⅶ、Ⅷ、Ⅸ3个烈度区间内呈增长的趋势，在烈度为Ⅸ的区域内LAP和LND均达到最大，分别为8.95%，37.27个／km2；在Ⅶ烈度区内LAP和LND均最小，为1.59%和1.57个／km2。烈度只对本次地震触发的滑坡的发育和分布产生影响，本文所统计的震前滑坡的发育不受本次地震烈度的影响，所以震前滑坡LAP和LND与烈度没有呈现出正相关的关系。震前Ⅸ烈度所覆盖的范围位于沙坝河与牛栏江的交会处，滑坡LAP较大，LND较小，该区域内滑坡的规模较大，如震前红石岩大滑坡；而Ⅷ烈度所分布的区域滑坡的规模较小。同震滑坡随烈度的变化呈现出较好的相关性，即滑坡的LAP和LND随着烈度值的增大，呈增长趋势，这也验证了烈度越大，越容易发生滑坡。

图10　滑坡分布与烈度关系图Fig．10　The relationship between landslide distribution and seismic intensity．

2.7河流

河岸两侧的岩土体具有侧向临空面，且长期受河流的冲刷侵蚀和水的浸泡作用，岩土体的成分和结构都会发生变化，所以河流对滑坡的发育也具有一定的控制作用。另外，河流的存在常常伴随着人类的活动，这对滑坡的发育产生的影响也是不容忽视的。在研究区内道路常沿河流分布，河流和道路附近的特殊地形为滑坡的发育提供了良好的条件。本文选取对研究区影响最大的河流——牛栏江和沙坝河为研究对象，以200m为缓冲距离建立缓冲区，然后对缓冲区进行重分类，按到河流的距离分为8类：0～200m，200～400m，400～600m，600～800m，800～1 000m，1 000～1 200m，1 200～1 400m，＞1 400m。通过统计得到的不同缓冲区内的滑坡的LAP和LND见图11。震前滑坡的LND总体上呈随到河流的距离的增加而减小的趋势。LAP在400～600m的缓冲区内最大为15.31%。LND在＜200m的分类达到最大为17.79个／km2。同震条件下，距离河流越远，滑坡的LAP越小；LND也大体呈递减趋势。LAP和LND在＜200m的分区内同时达到最大，分别为11.81%和53.83个／km2。河流与滑坡分布的关系统计图显示：震前、同震＜200m区域内滑坡的LAP、LND都较大，说明此区域比较容易发生大量小规模滑坡。震前400～600m范围内的滑坡的LAP最大，这与该分区内震前红石岩大滑坡的存在有着直接关系。通过对比发现，同震滑坡的LAP和LND与到河流的距离的相关性较震前更好，这说明地震力使得到河流的距离与滑坡易发性的相关性得到了加强。

3　结论

通过分析研究区内震前和同震滑坡与高程、坡度、坡向、曲率、岩性、烈度和河流7个因子的LAP和LND统计图，进而进行鲁甸地震震前与同震滑坡空间分布规律研究。结果表明，震前与同震滑坡受高程、坡度、坡向、曲率、岩性、烈度和河流的影响，表现出不同的空间分布特征。由于某些震前已经存在的滑坡的影响，震前滑坡LAP和LND随各因子的变化相关性不是很好。地震力的作用使滑坡的LAP和LND变化较震前的相关性有不同程度的加强。其中，河流的冲刷侵蚀作用在各个影响因子的统计结果中都有所体现。

震前滑坡的易发高程区间为＜1 200m，且在1 100～1 200m范围内滑坡的规模比较大，由于震前某些大型滑坡的存在，使震前1 300～1 400m范围内的LAP达到极大值；同震条件下在1200～1300m容易发育大量小规模的滑坡。震前滑坡的LAP和LND随着坡度的增大先减小后增大，表明无地震力作用时，坡度越大越易发生滑坡；震前滑坡在＜10°区由于河流的作用也较容易发生滑坡；同震滑坡的分布与坡度的相关性不是很大；由于地震波的地形放大效应，＞50°的坡度范围内也很容易发生滑坡。由于牛栏江的侧向剥蚀作用，震前近S向的边坡较容易失稳；同震滑坡的优势发育方向也大致为近S向。曲率在震前、同震都表现出了凹坡（曲率值为负）更易发生滑坡；但在地震力的作用下，凹坡表现出了极大的不稳定性，且滑坡易发性随着斜坡曲率值的减小而升高。灰岩夹白云质灰岩分布的区域由于位于河流附近，所以震前和同震都比较容易发育大量滑坡；玄武岩和火山角砾岩分布区的边坡在震前较稳定，在地震力的作用下其稳定性大大降低，进而发育大量滑坡；泥质白云岩分布区和白云岩、灰岩夹砂、页岩分别发育了震前的红石岩大滑坡和同震的红石岩滑坡。这次地震的烈度仅对地震的同震滑坡的分布产生影响；同震滑坡的分布与烈度呈现出了很好的相关性，且最易发区间为烈度最大的区间，即Ⅸ烈度区。震前、同震滑坡的分布随到河流距离的增加，滑坡的数量和规模都逐渐减小。其中，不管震前还是同震条件下，到河流的距离＜200m的区域都比较容易发生大量小规模滑坡。

图11　滑坡分布与河流关系图Fig．11　The relationship between landslide distribution and rivers．

致谢感谢国家测绘局地理信息局提供的超高分辨率航片与Google Earth平台，感谢审稿专家对本文详细的审阅以及所提出的宝贵意见。

常祖峰，周荣军，安晓文，等．2014．昭通－鲁甸断裂晚第四纪活动及其构造意义［J］．地震地质，36（4）：1260—1279．doi：10．3969／j．issn．0253－4967.2014.04.025．

CHANG Zu-feng，ZHOU Rong-jun，AN Xiao-wen，et al．2014．Late-Quaternary activity of the Zhaotong-Ludian Fault zone and its tectonic impliction［J］．Seismology and Geology，36（4）：1260—1279（in Chinese）．

程佳，刘杰，徐锡伟，等．2014．大凉山次级块体内强震发生的构造特征与2014年鲁甸6.5级地震对周边断层的影响［J］．地震地质，36（4）：1128—1242．doi：10．3969／j．issn．0253－4967.2014.04.023．

CHENG Jia，LIU Jie，XU Xi-wei，et al．2014．Tectonic characteristics of strong earthquakes in Daliangshan sub-blockand impact of the MS6.5 Ludian earthquake in 2014 on the surrounding faults［J］．Seismology and Geology，36（4）：1128—1242（in Chinese）．

房立华，吴建平，王未来，等．2014．云南鲁甸MS6.5地震余震重定位及其发震构造［J］．地震地质，36（4）：1173—1185．doi：10．3969／j．issn．0253－4967．2014．04．019．

FANG Li-hua，WU Jian-ping，WANG Wei-lai，et al．2014．Relocation of the aftershock sequence of the MS6.5 Ludian earthquake and its seismogenic structure［J］．Seismology and Geology，36（4）：1173—1185（in Chinese）．

冯文凯，许强，黄润秋．2009．斜坡震裂变形力学机制初探［J］．岩石力学与工程学报，28（增1）：3124—3130．

FENG Wen-kai，XU Qiang，HUANG Run-qiu．2009．Preliminary study on mechanical mechanism of slope earthquake-induced deformation［J］．Chinese Journal of Rock Mechanics and Engineering，28（Suppl）：3124—3130（in Chinese）．

胡金，李波，杨艳锋．2008．GIS在云南鲁甸县地质灾害易发性分区中的应用［J］．灾害学，23（1）：73—75．

HU Jin，LI Bo，YANG Yan-feng．2008．The application of GIS in geological hazard susceptibility zonation in Ludian County，Yunnan［J］．Journal of Catastrophology，23（1）：73—75（in Chinese）．

李西，张建国，谢英情，等．2014．鲁甸MS6.5地震地表破坏及其与构造的关系［J］．地震地质，36（4）：1280—1291．doi：10．3969／j．issn．0253－4967.2014.04.026．

LI Xi，ZHANG Jian-guo，XIE Ying-qing，et al．2014．Ludian MS6.5 earthquake surface damage and its relationship with structure［J］．Seismology and Geology，36（4）：1280—1291（in Chinese）．

李永华，徐小明，张恩会，等．2014．青藏高原东南缘地壳结构及云南鲁甸、景谷地震深部孕震环境［J］．地震地质，36（4）：1204—1216．doi：10．3969／j．issn．0253－4967.2014.04.021．

LI Yong-hua，XU Xiao-ming，ZHANG En-hui，et al．2014．Three-dimensional crust structure beneath SE Tibetan plateau and its seismotectonic implications for the Ludian and Jinggu eatthquakes［J］．Seismology and Geology，36（4）：1204—1216（in Chinese）．

刘洪兵，朱晞．1999．地震中地形放大效应的观测和研究进展［J］．世界地震工程，15（3）：20—25．

LIU Hong-bing，ZHU Xi．1999．Advance on topographic amplification effects of seismic response［J］．World Information on Earthquake Engineering，15（3）：20—25（in Chinese）．

祁生文，伍法权，刘春玲，等．2004．地震边坡稳定性的工程地质分析［J］．岩石力学与工程学报，23（16）：2792—2797．

QI Sheng-wen，WU Fa-quan，LIU Chun-ling，et al．2004．Engineering geology analysis on stability of slope under earthquake［J］．Chinese Journal of Rock Mechanics and Engineering，23（16）：2792—2797（in Chinese）．

石磊，贾晓东，陈石，等．2014．2014年云南鲁甸6.5级地震前重力变化特征与3维反演［J］．地震地质，36（4）：1217—1227．doi：10．3969／j．issn．0253－4967.2014.04.022．

SHI Lei，JIA Xiao-dong，CHEN Shi，et al．2014．Gravity changes before the Ludian，Yunnan，MS6.5 earthquake of 2014 and the three-dimensional inversion research［J］．Seismology and Geology，36（4）：1217—1227（in Chinese）．

唐立梅．2007．鲁甸县地质灾害特征、成因及防治区划研究［D］：［学位论文］．昆明：昆明理工大学．

TANG Li-mei．2007．Study on the characteristics，causes and prevention regionalization of geological hazards in Ludian County［D］．Dissertation．Kunming University of Science and Technology，Kunming（in Chinese）．

田颖颖，许冲，徐锡伟，等．2013．2013年岷县漳县MS6.6地震滑坡特征参数分析［J］．地震工程学报，35（4）：761—767．

TIAN Ying-ying，XU Chong，XU Xi-wei，et al．2013．Analysis of characteristic factors of landslides triggered by the July 22，2013 MS6.6 Minxian-Zhangxian earthquake［J］．China Earthquake Engineering Journal，35（4）：761—767（in Chinese）．

许冲，徐锡伟，戴福初，等．2011．2010年4月14日玉树地震滑坡空间分布与控制变量分析［J］．工程地质学报，19（4）：505—510．

XU Chong，XU Xi-wei，DAI Fu-chu，et al．2011．Analysis of spatial distribution and controlling parameters of landslides triggered by the April 14，2010 Yushu earthquake［J］．Journal of Engineering Geology，19（4）：505—510（in Chinese）．

许冲，徐锡伟，于贵华．2012．玉树地震滑坡分布调查及其特征与形成机制［J］．地震地质，34（1）：47—62．doi：10．3969／j．issn．0253－4967.2012.01.006．

XU Chong，XU Xi-wei，YU Gui-hua．2012．Study on the characteristics，mechanism，and spatial distribution of Yushu earthquake triggered landslides［J］．Seismology and Geology，34（1）：47—62（in Chinese）．

许冲，徐锡伟，吴熙彦，等．2013a．2008年汶川地震滑坡详细编目及其空间分布规律分析［J］．工程地质学报，21（1）：25—44．

XU Chong，XU Xi-wei，WU Xi-yan，et al．2013a．Detailed catalog of landslides triggered by the 2008 Wenchuan earthquake and statistical analysis of their spatial distribution［J］．Journal of Engineering Geology，21（1）：25—44（in Chinese）．

许冲，肖建章．2013b．2013年芦山地震滑坡空间分布分析：以太平镇东北方向的一个典型矩形区为例［J］．地震地质，35（2）：436—451．doi：10．3969／j．issn．0253－4967.2013.02.021．

XU Chong，XIAO Jian-zhang．2013b．Spatial analysis of landslides triggered by the 2013 MS7.0 Lushan earthquake：A case stusy of a typical rectangle area in the northeast of Taiping town［J］．Seismology and Geology，35（2）：436—451（in Chinese）．

许冲，徐锡伟，郑文俊，等．2013c．2013年四川省芦山“4.20”7.0级强烈地震触发滑坡［J］．地震地质，35（3）：641—660．doi：10．3969／j．issn．0253－4967.2013.03.018．

XU Chong，XU Xi-wei，ZHENG Wen-jun，et al．2013c．Landslides triggered by the April 20，2013 Lushan，Sichuan Province MS7.0 strong earthquake of China［J］．Seismology and Geology，35（3）：641—660（in Chinese）．

许冲，徐锡伟，郑文俊．2013d．2013年7月22日岷县漳县MS6.6级地震滑坡编录与空间分布规律分析［J］．工程地质学报，21（5）：736—749．

XU Chong，XU Xi-wei，ZHENG Wen-jun．2013d．Compiling inventory of landslides triggered by Minxian-Zhangxian earthquake of July 22，2013 and their spatial distribution analysis［J］．Journal of Engineering Geology，21（5）：736—749（in Chinese）．

许冲，徐锡伟，郑文俊，等．2013e．2013年甘肃岷县漳县6.6级地震触发滑坡及其构造分析［J］．地震地质，35（3）：616—626．doi：10．3969／j．issn．0253－4967.2013.03.015．

XU Chong，XU Xi-wei，ZHENG Wen-jun，et al．2013e．Landslides triggered by the 2013 Minxian-Zhangxian，Gansu Province MS6.6 eartqhuake and its tectonic analysis［J］．Seismology and Geology，35（3）：616—626（in Chinese）．

许冲，徐锡伟，沈玲玲，等．2014a．2014年鲁甸MS6.5地震触发滑坡编录及其对一些地震参数的指示［J］．地震地质，36（4）：1186—1199．doi：10．3969／j．issn．0253－4967.2014.04.020．

XU Chong，XU Xi-wei，SHEN Ling-ling，et al．2014a．Inventory of landslides triggered by the 2014 MS6.5 Ludian earthquake and its implications on several earthquake parameters［J］．Seismology and Geology，36（4）：1186—1199（in Chinese）．

许冲，徐锡伟．2014b．21世纪初几次大地震事件触发滑坡基础数据建设［J］．地震地质，36（1）：90—104．doi：10．3969／j．issn．0253－4967．2014．01．008．

XU Chong，XU Xi-wei．2014b．Construction of basic earthquake-triggered landslides dataset for several large earthquake events at the beginning of the twenty-first century［J］．Seismology and Geology，36（1）：90—104（inChinese）．

徐锡伟，程佳，许冲，等．2014a．青藏高原块体运动模型与地震活动主体地区讨论：鲁甸和景谷地震的启示［J］．地震地质，36（4）：1116—1134．doi：10．3969／j．issn．0253－4967．2014．04．015．

XU Xi-wei，CHENG Jia，XU Chong，et al．2014a．Discussion on block kinematic model and future themed areas for earthquake occurrence in the Tibetan plateau：Inspiration from the Ludian and Jinggu earthquakes［J］．Seismology and Geology，36（4）：1116—1134（in Chinese）．

徐锡伟，江国焰，于贵华，等．2014b．鲁甸6.5级地震发震断层判定及其构造属性讨论［J］．地球物理学报，57（9）：3060—3068．

XU Xi-wei，JIANG Guo-yan，YU Gui-hua，el al．2014b．Discussion on seismogenic fault of the Ludian MS6.5 earthquake and its tectonic attribution［J］．Chinese Journal of Geophysics，57（9）：3060—3068（in Chinese）．

杨光亮，申重阳，谈洪波，等．2014．云南鲁甸MS6.5地震震区地壳密度结构特征［J］．地震地质，36（4）：1145—1156．doi：10．3969／j．issn．0253－4967.2014.04.017．

YANG Guang-liang，SHEN Chong-yang，TAN Hong-bo，et al．2014．Crustal density structure of Yunnan Ludian MS6.5 earthquake area［J］．Seismology and Geology，36（4）：1145—1156（in Chinese）．

殷跃平．2009．汶川八级地震滑坡特征分析［J］．工程地质学报，17（1）：29—38．

YIN Yue-ping．2009．Features of landslides triggered by the Wenchuan earthquake［J］．Journal of Engineering Geology，17（1）：29—38（in Chinese）．

张盛峰，吴忠良，房立华．2014．双差（DD）定位地震目录能用于地震序列的统计地震学参数计算吗？［J］．地震地质，36（4）：1145—1259．doi：10．3969／j．issn．0253－4967．2014．04．024．

ZHANG Sheng-feng，WU Zhong-liang，FANG Li-hua．2014．Can the DD-relocated earthquake catalogue be used for the statistical parameters of an earthquake sequence？：A case study of the spatial distribution of B－values for the aftershocks of the 2014 Ludian MS6.5 earthquake［J］．Seismology and Geology，36（4）：1145—1259（in Chinese）．

赵静，刘杰，牛安福，等．2014．大凉山次级断块周边断层的闭锁特征［J］．地震地质，36（4）：1135—1144．doi：10．3969／j．issn．0253－4967.2014.04.016．

ZHAO Jing，LIU Jie，NIU An-fu，et al．2014．Study on fault locking characteristic about the Daliangshan sub-block［J］．Seismology and Geology，36（4）：1135—1144（in Chinese）．

赵旭，刘杰，冯蔚．2014．2014年云南鲁甸MS6.5地震震源运动学特征［J］．地震地质，36（4）：1157—1172．doi：10．3969／j．issn．0253－4967．2014．04．018．

ZHAO Xu，LIU Jie，FENG Wei．2014．The kinematic characteristics of the MS6.5 Ludian Yunnan earthquake in 2014［J］．Seismology and Geology，36（4）：1157—1172（in Chinese）．

周庆，江亚风，吴果，等．2014．芦山地震崩滑灾害空间分布及相关问题探讨［J］．地震地质，36（2）：344—

357．doi：10．3969／j．issn．0253－4967．2014.02.006．

ZHOU Qing，JIANG Ya-feng，WU Guo，et al．2014．Distribution of coseismic landslides in Lushan earthquake and discussion on related problems［J］．Seismology and Geology，36（2）：344—357（in Chinese）．

Cui P，Zhang J Q，Yang Z J，et al．2014．Activity and distribution of geohazards induced by the Lushan earthquake，April 20，2013［J］．Natural Hazards，73（2）：711—726．

Dai F C，Xu C，Yao X，et al．2011．Spatial distribution of landslides triggered by the 2008 MS8.0 Wenchuan earthquake，China［J］．Journal of Asian Earth Sciences，40（4）：883—895．

Harp E L，Jibson R W．1995．Inventory of landslides triggered by the 1994 Northridge，California earthquake［EB／OL］［Z］．http：∥geo-nsdi．er．usgs．gov／metadata／open-file／95－213／metadata．faq．html．

Keefer D K．1984．Landslides caused by earthquakes［J］．Geological Society of America Bulletin，95（4）：406—421．Keefer D K．2000．Statistical analysis of an earthquake-induced landslide distribution—The 1989 Loma Prieta，California event［J］．Engineering Geology，58（3）：231—249．

Liao H W，Lee C T．2000．Landslides triggered by the Chi-Chi earthquake［C］．In：Proceedings of the 21st Asian Conference on Remote Sensing，Taipei，1—2：383—388．

Liu L N，Xu C，Xu X W，et al．2014．Interactive statistical analysis of predisposing factors for earthquake-triggered landslides：A case study of the 2013 Lushan，China MS7.0 earthquake［J］．Environmental Earth Sciences，1—10．doi：10．1007／s12665－014－3758－8．

Qi S W，Xu Q，Lan H X，et al．2010．Spatial distribution analysis of landslides triggered by 2008.5．12 Wenchuan earthquake，China［J］．Engineering Geology，116（1）：95—108．

Sato H P，Hasegawa H，Fujiwara S，et al．2007．Interpretation of landslide distribution triggered by the 2005 Northern Pakistan earthquake using SPOT 5 imagery［J］．Landslides，4（2）：113—122．

Wang G L．2014．Comparison of the landslides triggered by the 2013 Lushan earthquake with those triggered by the strong 2008 Wenchuan earthquake in areas with high seismic intensities［J］．Bulletin of Engineering Geology and the Environment：1—13．

Wang H B，Sassa K，Xu W Y．2007．Analysis of a spatial distribution of landslides triggered by the 2004 Chuetsu earthquakes of Niigata Prefecture，Japan［J］．Natural Hazards，41（1）：43—60．

Xu C，Xu X W，Yao Q，et al．2013a．GIS-based bivariate statistical modelling for earthquake-triggered landslides susceptibility mapping related to the 2008 Wenchuan earthquake，China［J］．Quarterly Journal of Engineering Geology and Hydrogeology，46（2）：221—236．

Xu C，Xu X W，Yu G H．2013b．Landslides triggered by slipping-fault-generated earthquake on a plateau：An example of the 14 April 2010，MS7.1，Yushu，China earthquake［J］．Landslides，10（4）：421—431．

Xu C，Xu X W．2014a．Statistical analysis of landslides caused by the MW6.9 Yushu，China，earthquake of April 14，2010［J］．Natural Hazards，72（2）：871—893．

Xu C，Xu X W．2014b．The spatial distribution pattern of landslides triggered by the 20 April 2013 Lushan earthquake of China and its implication to identification of the seismogenic fault［J］．Chinese Science Bulletin，59（13）：1416—1424．

Xu C，Shyu J B，Xu X W．2014c．Landslides triggered by the 12 January 2010 Port-au-Prince，Haiti，MW＝7.0 earthquake：Visual interpretation，inventory compiling，and spatial distribution statistical analysis［J］．Natural Hazards and Earth System Sciences，14（7）：1789—1818．

Xu C，Xu X W，Shyu J B H，et al．2014d．Landslides triggered by the 22 July 2013 Minxian-Zhangxian，China，MW5.9 earthquake：Inventory compiling and spatial distribution analysis［J］．Journal of Asian Earth Sciences，92：125—142．

Xu C，Xu X W，Yao X，et al．2014e．Three（nearly）complete inventories of landslides triggered by the May 12，2008 Wenchuan MW7.9 earthquake of China and their spatial distribution statistical analysis［J］．Landslides，11（3）：441—461．

Abstract

On August 3，2014，an MS6.5 earthquake happened in the Ludian County，Zhaotong City of Yunnan Province．This earthquake caused a large number of landslides．In order to study the spatial distribution of the coseismic and pre-earthquake landslides，a 44.13km2area at the junction of Ludian County，Qiaojia County and Huize County along the Niulanjiang River is selected in this study．By visual interpretation of the Google Earth pre-earthquake high resolution images and the coseismic aerial data of 0.2m resolution of this area，the landslide databases of pre-earthquake and coseismic are established．The result shows that there are 284 pre-earthquake landslides，and 1053 earthquake-induced landslides．Then with the help of 10m×10m resolution DEM data and the GIS，the extracted factors of elevation，slope angle，slope aspect，curvature，lithology，earthquake intensity and drainages are used to analyze the spatial distribution of the coseismic and pre-earthquake landslides by adopting LAP（Landslide Areas Percentage）and LND（Landslide Number Density）．The results show that areas with elevation＜1 200m and 1 200～1 300m are prone to landslides whatever pre-earthquake or coseismic．With the slope gradient increasing，it is much more prone to landslide，and the area of＜10°，close to the rivers，is also much susceptible．The advantage slope aspect is almost near S direction．Concave slope（when the curvature is negative）is much susceptible to landslides，and with the curvature decreasing，the landslide susceptibility gets higher．The region of limestone with dolomitic limestone is sensitive to landslide；in the areas consisting of basalt and volcanic breccia，the slope stability is greatly reduced under the effect of seismic force．The larger the intensity is，the more landslides happened．For either pre-earthquake or coseismic landslides，there is a positive correlation between landslide spatial distribution and the distance to rivers．The large pre-earthquake landslides have effective influence on LAP．

SPATIAL DISTRIBUTION ANALYSIS OF COSEISMIC AND PRE-EARTHQUAKE LANDSLIDES TRIGGERED BY THE 2014 LUDIAN MS6.5 EARTHQUAKE

TIAN Ying-ying1，2）XU Chong1）XU Xi-wei1）WU Sai-er2）CHEN Jian2）
1）Key Laboratory of Active Tectonics and Volcano，Institute of Geology，China Earthquake Administration，Beijing100029，China
2）School of Engineering and Technology，China University of Geosciences，Beijing100083，China

Ludianearthquake，pre-earthquakelandslides，coseismiclandslides，spatial distribution，landslide density

P315．9

A文献标识码：0253－4967（2015）01－0291－16

10．3969／j．issn．0253－4967．2015．01．023

田颖颖，女，1989年生，中国地质大学（北京）地质工程专业在读硕士研究生，研究方向为滑坡危险性评价，E－mail：t190403098＠163.com。

2014－10－28收稿，2015－01－23改回。

国家自然科学基金（41472202，41202235）和中国地震局“云南鲁甸6.5级地震专题研究”共同资助。

许冲，副研究员，E－mail：xuchong＠ies．ac．cn；xc11111111＠126．com。