基于侵蚀循环理论的地震触发崩塌滑坡灾势评价

段书苏 姚令侃 郭沉稳

摘要：构建了一种与选线原则方案精度相匹配的、区域性的预测地震触发崩塌滑坡灾势的方法.论述了在侵蚀循环（幼年期-壮年期-老年期）的一个地貌发育周期内，流域谷坡将相应经历向临界坡发展、达到临界坡、偏离临界坡的过程，提出了根据流域演化发育阶段预测山地灾害危险性的原理；基于河谷与谷坡的反馈机制，提出通过点坡度变化确定纵剖面特征点，进而区分不同的发育阶段，解决了斯特拉勒分析方法只能适应于小流域的问题.通过汶川地震、玉树地震、芦山地震这3次21世纪以来我国震级在Ms7.0以上的大地震实震资料，检验了本方法的正确性.流域的斯特拉勒积分在0.5～0.6之间时，崩塌滑坡灾势最严重；距上值越远，灾势越轻.

关键词：侵蚀循环理论； 地震触发崩塌滑坡； 灾势预测； ArcGIS

中图分类号：U211.9 文献标识码：A

21世纪以来，我国发生了3次Ms.7.0以上的大地震（汶川地震、玉树地震、芦山地震）.在相同的地震烈度区内，震后次生山地灾害的严重性差异却很大.地震触发的崩塌滑坡灾害态势的评价引起了各国学者的重视.

区域性研究地震触发崩塌滑坡严重性的方法，从定性分析法发展至综合指标法，到现在流行的Newmark累计位移法.定性分析法和综合指标法，都是通过统计历史地震中滑坡与各种指标之间的关系，得出结论.例如：美国学者Keefer[1]，Rodríguezpeces等[2]对世界范围内的滑坡进行统计，归纳了地震滑坡分布与震中距、岩性等之间的关系； 2000年，我国学者乔建平、王余庆[3-5]等运用综合指标法对影响地震滑坡的各种因素进行了区分和定量化，进而对区域地震滑坡进行评价分级；汶川地震后，黄润秋[6]等建立了强震作用下斜坡危险性评价的指标体系；许冲[7-8]利用ArcGIS技术平台，采用综合指标法对玉树、汶川地震滑坡危险度进行了评价.研究结论反映了历史地震崩塌滑坡的统计特性，成果的普适性受统计样本代表性的限制.目前国际上常用的地震滑坡区划模型是Newmark累计位移法[9]，但Newmark模型所需要的参数较多（包括滑坡深度及滑动面形状、滑动面黏聚力、内摩擦角等），更适用于岩石工程数据详实的小区域.

本文提出，地貌侵蚀循环的过程就是流域内坡体在各发育阶段依次经历向临界坡发展、达到临界坡、偏离临界坡的演变过程.利用区域内坡体偏离临界坡的程度，可以定性地评定流域内斜坡重力灾害的严重性.从而建立起一种基于侵蚀循环理论评判地震触发崩塌滑坡灾势的方法.

1基于谷坡状态的地震触发崩塌滑坡灾势预测原理

基于发生学的理论，Davis在1899年首次创立了侵蚀循环学说（Theory of the Cycle of Erosion），认为地块从开始上升到被逐渐剥蚀夷平，直至降低到起伏不大的地面或者接近基准面的准平原之间，存在着连续的、同时又有阶段性的剥蚀过程和地表形态.在地表发育的过程中，Davis强调构造、作用和时间（侵蚀阶段）这3个要素之间的相互作用影响[10]，进而将循环过程中的地形发展分为3个阶段：地形起伏不大，河间地广阔平坦的幼年期，地面主要由谷坡和狭窄的分水岭组成的壮年期，和具有残丘的准平原的老年期.

按照地貌循环理论，侵蚀轮回中坡地发育过程如图1所示.由图可知，当地貌处于幼年期时，坡度最陡，在谷坡之上存在大量的上个地貌循环留存下来的夷平面.从幼年期发展到壮年期的阶段，谷底强烈下切，海拔持续下降；夷平面不断被侵蚀，面积不断减小，高程却没有明显变化；谷坡坡度持续减小.从壮年期发展到老年期阶段，谷底与谷坡的高程都在缓慢下降，但谷坡下降的速度比谷底快.这导致流域内谷坡坡度持续变缓，夷平作用持续进行，直至准平原.

地貌的形成和发展是内、外营力相互作用的结果.内营力趋向于使山体隆升，是使流域内坡度增加的过程；外营力作用趋向于使山体高度降低、削平，是使流域内坡度减小的过程.由于二者对山地地貌塑造的反向效应，山体坡度最大只能达到一个特定值，即所谓的临界坡度.那么，在一个侵蚀循环内，地貌经历幼年期、壮年期、老年期的同时，流域内的坡体相应经历了向临界坡度发展、达到临界坡、偏离临界坡的演变过程.处于临界坡的斜坡系统，在地震、降雨等外界扰动下，极易发生失稳破坏，造成大规模的崩塌滑坡；反之，坡度小于临界坡时，斜坡系统就具有一定的安全裕度，而且偏离临界坡度越远，安全裕度越大，发生崩塌滑坡可能性越小.因此可以以区域内坡体偏离临界坡的程度，作为评估流域山地灾害发生危险性的依据.这就是基于侵蚀循环理论预测地震触发崩塌滑坡灾势的原理.

3地震触发崩塌滑坡灾势评估案例分析

现以21世纪以来我国发生的7级以上3次大地震为例，通过对地震触发崩塌滑坡实震资料的对比分析，对该理论进行验证，同时也对具体操作方法进行说明.

3.1玉树地震震区斯特拉勒积分

2010年4月14日，玉树县发生Ms7.1地震.境内主要水系包括通天河、扎曲、巴曲等，地形以高海拔、低起伏为主.震区所在计算流域的斯特拉勒积分计算步骤如下：

第一步，划定研究区域.震区位于通天河和金沙江的交汇地段，那么我们圈定的研究区域的流域应该是包括沱沱河、通天河、金沙江上游在内的长江上游河段，研究区域的水系格局如图3所示.

沱沱河（长江正源）周围发育着色林错、赤布张错、乌兰乌拉湖、西金兰湖等高原内陆湖，扎加藏布江、曾松曲等内流河等；楚玛尔河（长江北源）周围发育着格尔木河等内流河，可可西里湖、库赛湖、盐湖、达布逊湖等内陆湖；当曲（长江南源）周围发育着怒江的上游那曲、澜沧江的上游扎曲；通天河的东侧与黄河的上游约古列宗渠毗邻；在金沙江河段，东西两侧基本平行发育着澜沧江和雅砻江.根据上述的流域格局圈定研究区域周界为：雅砻江-古宗列渠-格尔木河-库赛湖-可可西里湖-乌兰乌拉湖-赤布张错-扎加藏布江-澜沧江，提取相应区域的DEM.

第二步，划定流域边界.长江的源头及通天河属于广大的藏北腹地，区域内内流河与外流河之间往往没有明显的分水岭，河流弯转曲折，分水线在平面上呈犬牙交错状；金沙江、澜沧江、雅砻江之间存在宽广的夷平面，地貌特征不明显.本文采用了较高精度的DEM（30 m×30 m），生成研究区域的河网及集水区如图4所示.针对流域边界复杂情况，需注意采取以下措施以保证边界划定精度：①圈定分析区域时，区域范围务必圈定至相邻的独立水系，②在困难区域，选择高精度的DEM分析，在流域划分的过程中应该与现有的其他数据资料进行相互验证校核，从而精准划定流域边界.

第三步，提取河谷纵剖面特征点，确定计算区域.通过提取河流点坡度和纵剖面复合图（图5），可以看出，在玉树下游的不远处，海拔4 000 m的地方，坡度突然降低至0附近，疑似夷平面.在夷平面的上游，河谷纵剖面变化不大，均是比较缓的坡度；而在夷平面的下游，河谷坡面变化剧烈.这是由于青藏高原的分阶段、非均一隆升，使高原东缘的外流河产生溯源侵蚀，在侵蚀到达的地方，形成高山峡谷，在侵蚀未达到的地方，盆地面内平坦开阔，切割微弱.夷平面上、下游2个地貌单元差异明显[19].选择这一夷平面为流域计算的出口点，如图4所示.

3.2汶川地震、芦山地震震区的斯特拉勒积分

“5·12”汶川大地震和“4·20”芦山地震2次地震的主震区均属于龙门山地区，分别位于青衣江流域、岷江流域范围内.划定研究区域后，该两流域与相邻流域分水岭地貌明显，可通过流域分析直接划定流域边界.龙门山、岷山等一系列的山脉组成了青藏高原东缘的地形陡变带，其间的河流走向多与龙门山垂直，以深切河谷为主要特征，进入四川盆地后河流弯转曲折，流速缓慢，流域河道控制点即为山区河段和平原河段的分界点.通过提取青衣江和岷江的全河段的河流点坡度和纵剖面复合图（图7），发现岷江在都江堰处存在拐点，青衣江在芦阳镇处存在拐点，可将都江堰和芦阳镇分别作为岷江和青衣江的河道控制点（计算流域的出口点），计算流域的斯特拉勒曲线和斯特拉勒积分如图8所示.

根据以上计算可知，青衣江和岷江的Strahler积分值分别为0.46和0.52.

距河源的距离/km

3.3计算结果与实震坡体资料对比

芦山地震震中位于北纬30.3°， 东经103°，震源深度约13 km，主震区沿龙门山断裂带分布，呈东北—西南走向，南起凉山州甘洛县，往东北经汉源、荥经、芦山至大邑县，主要位于雅安市境内.震中烈度达Ⅸ度.

利用震后遥感影像资料进行人工目视解译是大面积获取震区崩塌滑坡信息的主要方法.芦山地震后的半年时间里，我们通过多种途径收集了震区航空、航天遥感影像资料，包括：1）中国科学院遥感与数字地球研究所提供的三批航空遥感数据，覆盖芦山、宝兴、邛崃等县市约5 000 km2.第一批航片获取时间为2013年4月20日10：30～12：40，分辨率0.6 m；第二批航片获取时间为4月20日15：00～17：00，包括0.4 m和2 m两种分辨率；第三批航片获取时间为4月21日上午，包括0.4 m和2 m两种分辨率；2）四川省测绘地理信息局芦山地震信息发布平台公布的芦山震中区的航片影像资料（影像获取截止时间为4月25日，分辨率0.5 m）；3）国家测绘局公布的龙门乡、太平镇、宝盛乡三地的航片影像资料（影像获取时间为4月20日18：28，分辨率为0.16 m） .对以上遥感影像资料进行几何纠正、融合、拼接、图像增强等数据处理，获得芦山震区的遥感覆盖面积208 km2，覆盖整个Ⅸ度区.区域内的灾害点分布如图9所示.

表1所列是3次大地震所在流域地震崩塌滑坡发生条件的评判结果，以及3次地震中Ⅸ度烈度区实震资料统计.

据表1可知，芦山地震和汶川地震震区均处于壮年中期，谷坡处于临界坡状态.玉树地震震区处于老年期，谷坡已偏离临界坡.考虑到地震烈度条件相同时，3场地震崩塌滑坡的严重程度才具有可比性，因此均选取Ⅸ度烈度区的崩塌滑坡密度作为比较指标（玉树、芦山地震最大烈度均为Ⅸ度）.实震资料表明，芦山地震和汶川地震均触发了大量的崩塌滑坡，而玉树地震次生灾害相对弱得多.这与芦山地震、汶川地震震区具备地震触发大规模崩塌滑坡的条件，而玉树地震震区不具备此条件的评判结论是吻合的.

至此，基于侵蚀循环理论的地震触发崩塌滑坡灾势预测的理论，得到了21世纪以来3次7级以上大地震实震资料的检验.

4结论

1）基于地貌循环理论，论证了在一个地貌发育循环周期内，流域谷坡将经历向临界坡发展、到达临界坡，偏离临界坡的演变阶段，在降雨、地震等触发条件具备时，山地灾害相应地呈现从发展到旺盛再到衰退的规律，这就是根据流域演化发育阶段预测地震触发崩塌滑坡危险性的原理；提出了根据河床特征点划分流域地貌单元的理念，从而解决了面积高程分析理论应用于复杂大流域的适用性问题；归纳了利用DEM和ArcGIS技术计算斯特拉勒积分的程式.从而在仅利用社会公共资料的条件下，建立起一种与原则选线阶段精度要求相匹配的区域性地震触发崩塌滑坡危险性预测方法，并通过了21世纪以来3次7级以上大地震实震资料的检验.本文提出的方法亦可为其他诱因的区域性山地灾害危险性预测提供借鉴.

2）高烈度地震山区的铁路选线，是在节省工程投资和减轻未来地震灾害风险的矛盾中起着统筹规划作用的多目标决策过程，线路原则方案的选择是风险调控的首要环节.对于谷坡不具备地震触发大规模崩塌滑坡条件的流域，可以采用直接通过方案；对于预测地震触发崩塌滑坡灾势严重的流域，若采用通过方案，宜以隧道为主穿越，而这一般会导致工程造价大幅度增加，因此还应与大范围绕避方案进行经济技术比较后确定合理方案.显然，本文工作将为上述方案论证提供重要依据.

参考文献

[1]KEEFER D K. Landslides caused by earthquakes [J]. Geological Society of America Bulletin， 1984， 95 （4）： 406-421.

[2]RODRGUEZPECES M J， PREZ J L， GARCAMAYORDOMO J，et al. Applicability of Newmark method at regional， subregional and site scales： seismically induced Bullas and La Paca rockslide cases （Murcia， SE Spain）[J]. Natural Hazards， 2011， 59（2）： 1109-1124.

[3]乔建平，蒲晓虹. 川西南滇北接壤带地震滑坡概述[J]. 山地研究， 1987（3）： 181-186.

QIAO Jianping，PU Xiaohong. An introduction to earthquake landslides in the contiguous area between Southwest Sichuan and North Yunnan [J].Mountain Research， 1987（3）： 181-186.（In Chinese）

[4]王余庆，辛鸿博，高艳平，等. 预测岩土边坡地震崩滑的综合指标法研究[J]. 岩土工程学报， 2001， 23（3）： 311-314.

WANG Yuqing，XIN Hongbo，GAO Yanping，et al. Study on comprehensive index method for predicting earthquakeinduced landslides [J]. Chinese Journal of Geotechnical Engineering， 2001， 23（3）： 311-314. （In Chinese）

[5]丁彦慧，王余庆，孙进忠，等. 地震崩滑预测方法及其工程应用研究[J]. 工程地质学报， 2000， 8（4）： 475-480.

DING Yanhui，WANG Yuqing，SUN Jinzhong，et al. Research on the method for prediction of earthquakeinduced landslides and its application to engineering projects[J].Journal of Engineering Geology， 2000， 8（4）： 475-480. （In Chinese）

[6]黄润秋. 汶川地震地质灾害研究[M]. 北京： 科学出版社， 2009： 203-257.

HUANG Runqiu. Geohazard assessment of the Wenchuan earthquake [M]. Beijing： Science Press， 2009： 203-257. （In Chinese）

[7]许冲，徐锡伟，于贵华. 玉树地震滑坡分布调查及其特征与形成机制[J]. 地震地质， 2012， 34（1）： 47-62.

XU Chong，XU Xiwei，YU Guihua. Study on the characteristics， mechanism， and spatial distribution of Yushu earthquake triggered landslides[J]. Seismology and Geology， 2012， 34（1）： 47-62. （In Chinese）

[8]许冲，戴福初，陈剑，等. 汶川Ms8.0地震重灾区次生地质灾害遥感精细解译[J]. 遥感学报， 2009， 13（4）： 745-762.

XU Chong，DAI Fuchu，CHEN Jian， et al. Identification and analysis of secondary geological hazards triggered by a magnitude 8.0 Wenchuan Earthquake [J]. Journal of Remote Sensing， 2009， 13（4）： 745-762. （In Chinese）

[9]NEWMARK N. Effects of earthquakes on dams and embankments [J]. Geotechnique， 1965， 15（2）： 139.

[10]DAVIS W M. The geographical cycle[J]. The Geographical Journal， 1899， 14（5）： 481-504.

[11]STRAHLER A N. Hypsometric （areaaltitude） analysis of erosional topography [J]. Geological Society of America Bulletin， 1952， 63 （11）： 1117-1142.

[12]WEISSEL J K， PRATSON L F， MALINVERNO A. The lengthscaling properties of topography[J]. Journal of Geophysical Research， 1994，99（B7）：799-808.

[13]艾南山. 侵蚀流域系统的信息熵[J]. 水土保持学报， 1987（2）： 1-8.

AI Nanshan. Comentropy in erosion drainagesystem[J]. Journal of Soil and Water Conservation， 1987（2）： 1-8. （In Chinese）

[14]刘怀湘，王兆印，陆永军，等.山区下切河流地貌演变机理及其与河床结构的关系[J].水科学进展，2011，22（3）：367-372.

LIU Huaixiang， WANG Zhaoyin， LU Yongjun， et al. Mechanism of morphological evolution in incised river and its relationship to streambed structures[J]. Advances in Water Science， 2011，22（3）：367-372. （In Chinese）

[15]左大康. 现代地理学辞典[M].北京：商务印书馆，1990：182.

ZUO Dakang. Dictionary of modern geography[M]. Beijing： Commercial Press，1990：182. （In Chinese）

[16]杨景春，李有利. 地貌学原理[M]. 北京： 北京大学出版社， 2005： 24-26.

YANG Jingchun，LI Youli. Theory of geomorphology [M]. Beijing： Beijing University Press， 2005： 24-26. （In Chinese）

[17]HAYAKAWA Y S， OGUCHI T. DEMbased identification of fluvial knickzones and its application to Japanese mountain rivers [J]. Geomorphology， 2006，78（1）：90-106.

[18]吴立新，史文中. 地理信息系统原理与算法[M]. 北京： 科学出版社， 2003： 441.

WU Lixin，SHI Wenzhong. Geographic information system principle and algorithm [M]. Beijing： Science Press， 2003： 441. （In Chinese）

[19]中国科学院青藏高原综合科学考察队. 西藏河流与湖泊[M]. 北京： 科学出版社， 1984： 5-215.

Scientific Expedition of the Chinese Academy of Sciences to the QinghaiXizang Plateau. Rivers and lakes of Xizang （Tibet）[M]. Beijing： Science Press， 1984： 5-215. （In Chinese）

[20]殷跃平，张永双，马寅生，等. 青海玉树M_S7.1级地震地质灾害主要特征[J]. 工程地质学报， 2010， 18（3）： 289-296.

YIN Yueping，ZHANG Yongshuang，MA Yansheng，et al. Reserch on major charcteristics of geohazards induced by the Yushu Ms7.1 earthquake[J]. Journal of Engineering Geology， 2010， 18（3）： 289-296. （In Chinese）