地震动力作用触发的斜坡崩滑响应研究现状及展望

杨旭朝，崔芳鹏

（1.天津市市政工程设计研究院，天津 300201；

2.中国矿业大学（北京）地球科学与测绘工程学院，北京 100083；

3.成都理工大学 地质灾害防治与地质环境保护国家重点实验室，四川 成都 610059）

0 引言

近年来，强震触发的斜坡崩、滑、流等地质灾害因其巨大的致灾力引起了国内外广泛关注与研究。据多次典型地震的致灾统计数据看［1-3］，强烈地震触发的斜坡崩、滑、流地质灾害造成的人员伤亡数量占到地震伤亡总人数的近三分之一，其危害程度“不亚于一场大地震”［4］。据统计，举世瞩目的汶川8.0级强震更是触发了15000多处滑坡、崩塌、泥石流地质灾害，估计直接造成2万人死亡。地质灾害隐患点达12700多处，以崩塌体灾害增加最为显著，威胁近百万人的生命安全［5］。

总体而言，国内外对地震作用触发的斜坡崩、滑、流灾害研究主要集中在以下几个方面：一是地震诱发的斜坡崩滑灾害特征及影响因素统计学研究，即着重从分析地震诱发的斜坡崩滑灾害实例入手，从统计学角度论证地震崩滑灾害发育特征同地震参数（如地震震级、地震烈度、震源深度、震中距等）及斜坡体特征（如坡高、坡角、岩性、构造、水文地质条件等）之间的关系；二是地震诱发的斜坡崩滑灾害形成机制及动力响应特征研究；三是地震作用下斜坡稳定性评价及致灾预测研究。

1 地震诱发的斜坡崩滑灾害特征及影响因素统计学分析

国外方面，Keefer依据地震滑坡运动的特性、内部破裂、地质环境等，将地震滑坡分为三大类：Ⅰ类滑坡为高度破裂或严重高度破裂，滑坡物质含大量的大石块、小一些的岩石碎块、小土块和土粒。该类滑坡典型的特征是起源于陡的斜坡；Ⅱ类滑坡主要指存在分离的剪切面，该类滑坡起源于坡度中等的斜坡，厚度大于3m；Ⅲ类滑坡指滑动机理主要以流相为主，表面块体由内部裂缝造成破裂而成，并且由于相互分离而形成地堑，块体之间可能充填液相物质。Keefer还通过对历史地震资料的统计分析，得出了震级与滑坡体积的统计公式：

log V=1.44（±0.21）Mw-2.34（±1.5）5.3≤Mw≤8.6

式中：V——滑坡总体积（m3）；

Mw——矩震级［6］。

Richter研究了地震滑坡与地震烈度之间的关系［7］，Prestininzi也做了类似的工作，他对意大利从公元前461年到1992年间的地震目录进行了分析，给出了地震烈度与地震滑坡之间的图形关系［8］。Papadopoulos等根据1000年到1995年希腊境内的47个地震事件：（Ms=5.3～7.9），统计出地震滑坡与震中距的关系为：

log Rc=0.75Ms-2.98

式中：R c——地震触发滑坡的最大震中距离（km）［9］。

Keefer根据1989年 LomaPrieta，California地震（M=6.9）震后的调查，得出了滑坡分布（密度）与震中距、岩性、坡角等的关系［10］。Rodriguez统计了1980～1997年全世界范围内地震造成的岩土边坡滑坡灾害情况，分析了地震造成滑坡的最大密度、影响滑坡范围、滑坡数量以及滑坡类型。分析了滑坡面积、滑坡震中距与震级之间的关系。分析了不同裂度地区的地震数量，得到Ⅶ和Ⅷ级烈度下滑坡数量较多的结论［11］。Parise和Jibson根据1994年美国Northridge地震情况的调查，分析了滑坡面积、形式，以及滑坡分布、影响因素（地震强度、坡角）等与滑坡的关系，指出地震滑坡的坡角多处于 20°～35°之间［12］。Julian和Carlos于2002年统计了中美洲地区1898～2001年地震造成的滑坡情况，分析了地震作用的影响面积、震中距以及影响情况［13］。

国内方面，周本刚等根据我国西南地区的地震滑坡特点，将地震滑坡按其运动方式划分为推移式滑坡、牵引式滑坡、溜滑性滑坡和崩塌性滑坡四大类，并对每一类进行了细致的描述［14］。薛守义等在考虑了动荷载作用特点的基础上，将边坡的变形破坏形式分为滑动型、崩塌型、塑性变形破坏、层体弯折四大类，各大类之下又分若干亚类［15］。李天池根据我国区域地质和地貌条件的特性，通过回归计算，得出了单个地震Ⅶ度以上烈度区的面积和震级的近似关系［16］。辛鸿博等通过对我国过去近800年间发生的125次地震及其触发的285例典型边坡崩滑实例进行分析研究，提出了以场地烈度、震级、震中距为标准的三参数边坡地震崩滑初判准则［17］。孙崇绍等通过对中国历史地震资料（1500～1949年）统计分析，得出以下结论：地震引起的崩塌、滑坡多在5级以上的震区，6级以上的震区内崩、滑的数量显著增大，根据历史地震记载中所反映的烈度，崩、滑大多发生在Ⅶ度及Ⅶ度以上的地区，Ⅷ度以上发生的可能性急剧增大。几次规模特大、破坏特重的滑坡都发生在Ⅸ度及Ⅸ度以上的地区［18］。李忠生依据《中国地震目录》及其他主要地震目录，对5级地震以下的斜坡崩滑灾害做了统计分析，得出了诱发地震崩滑灾害的最小震级统计结果［19］。尹光华等利用野外考察资料，对1812年新疆尼勒克地区的滑坡进行了统计分析和研究，认为尼勒克地震产生的滑坡规模大，数量多，分布集中，其规模和分布密度与烈度成正比，有由震中向外变小、变稀的特点［20］。王余庆等提出了地震滑坡最大震中距与震级之间的关系，并初步建立了分析地震滑坡的综合指数法（综合指数考虑坡高系数、岩性条件、坡角影响、地震强度以及降水系数），得到边坡危害指数，并结合云南丽江地震后的调查进行了验证［21-22］。杨涛等通过对四川地区3.5级以上地震诱发的崩塌滑坡事件进行分析，指出了四川地区地震滑坡的分布特点。还分析了崩塌面积与震级、震中距、震源深度之间的关系，提出了诱发边坡崩塌的最小震级［23］。崔芳鹏等通过对重庆小南海地震崩滑体的研究，着重分析了地震作用诱发的斜坡崩滑体抛射方向与发震构造破裂机制、斜坡体构造、地震作用力及崩滑堆积体粒径分布之间的关系［24］。

以上多为对多次地震诱发的众多斜坡崩滑灾害统计分析得出的统计结论，对于认识地震作用诱发的斜坡崩滑灾害总体特征及相关影响因素具有指导和参考作用，但其对认识特定地质条件下特定地震诱发的斜坡崩滑灾害参考价值则有所降低。

2 地震诱发的斜坡崩滑灾害形成机制及动力响应特征研究

近年来，对地震诱发的斜坡崩滑灾害形成机制有较为深入的研究。

地震动载对岩质边坡的破坏作用主要有：①地震荷载通过岩体本身不同层面与结构不连续的界面起反射和折射作用，导致超压增大；或者致使欠稳定的岩坡发生掉块、局部崩塌和呈台阶滑坡；②地震和其它外界因素共同作用造成边坡失稳，特别是降雨、洪水和地下水位的变化，加剧了破坏的形成［25］。张倬元等认为地震对边坡稳定性的影响表现为累积和触发效应两个方面。前者主要表现为地震作用引起边坡岩体结构松动、破裂面、弱面错位和孔隙水压力累积上升等。后者则主要表现为地震的作用造成边坡中软弱层的触变液化、砂层液化以及使处于临界状态的边坡瞬间失稳［26］。胡广韬等提出了边坡动力失稳机制的坡体波动振荡加速效应假说。当天然边坡处于或接近于极限状态时，地震产生的惯性力会激发滑坡体骤然下滑，出现启程剧动，产生“启程剧发速度”，并把这种效应称之为坡体振荡加速效应［27］。他根据地震纵横波的不同特性（即地震纵波为拉压波，质点振动方向与波前进方向一致，只在斜坡体中产生体积变形而无畸变，形成斜坡体的近竖直抛起和拉离效应；而横波为剪切波，质点振动方向与波前进方向垂直，在斜坡体中只产生形状改变而无体积改变，形成斜坡体的水平震荡移动效应），分别推导出坡体在垂直波动震荡的抛起作用、拉离作用和水平波动震荡产生的水平震荡移动作用下斜坡体的位移量公式［28］，但他忽略了地震作用过程中导致的孔隙水压力增大及累积作用对斜坡体变形破坏所产生的巨大破坏。毛彦龙等认为地震对滑坡形成的影响，主要是通过坡体波动振荡来产生；坡体波动振荡在边坡岩土体变形破坏过程中产生三种效应：累进破坏效应、启动效应和启程加速效应［29］。祁生文等认为地震边坡的失稳是由于地震惯性力的作用以及地震产生的超静孔隙水压力迅速增大和累积作用这两个方面造成的，这两种作用使得边坡沿着某一滑动面的下滑力增大，抗滑力减少，致使边坡失稳破坏。不同的边坡破坏类型，导致边坡动力失稳的主导因素也不同，一般来讲，塑性流动失稳破坏是孔隙水压力的累积作用起主导引起的，崩塌型、层体弯折型则是地震惯性力起决定作用，而对于滑动型破坏则视具体条件而定。对于后发型边坡失稳，则是由于地震的累积效应导致边坡岩体力学参数的降低，为后来的外地质营力创造了条件，最终导致了边坡的失稳破坏［109］。岩质边坡中不可避免地存在断层、节理、软弱夹层等结构面，当地震波通过这些结构面时，将发生复杂的反射、折射现象，由于岩体的抗拉强度远远低于它的抗压强度，反射拉伸波将引起岩体受拉破坏；而地震荷载的往复作用，将导致岩体的抗剪强度降低、下滑力增大，岩体可能发生剪切破坏。但这两种破坏形式决不会单独发生，而是拉剪共同作用的结果，其破坏程度与岩体力学参数、结构面空间展布等有着密切的关系［110］。

“5.12”汶川地震诱发的斜坡崩滑灾害引起了国内外学者的广泛关注，殷跃平等在对震区大量斜坡崩滑灾害调查研究基础上，认为汶川地震极震区斜坡崩滑灾害一般不具连续完整的破坏带（面），而呈现明显的“尖点突起”或“边缘突出”特征，反映出上部滑体被地震力振动解体，甚至抛掷后与下部滑床边缘撞击的动力学特性。依据强震地面运动记录和大量实例调查表明，触发斜坡崩滑的地震竖向力作用非常明显，大量崩滑灾害经历了初始斜坡（风化碎裂岩体）→地震抛掷→撞击崩裂→高速滑流的动力学特征。在该过程中，斜坡崩滑的动力学过程具有鲜明的抛掷效应、碰撞效应、刮铲效应及气垫效应［30，31］。黄润秋等通过对震区大量斜坡崩滑灾害进行现场调查、遥感解译和研究后认为汶川地震地质灾害具有以下鲜明特征：①在区域上具有沿地震断裂带呈带状分布和沿河流水系成线状分布的特点；②具有明显的上盘效应，即发震断裂上盘地质灾害发育密度明显大于下盘，且上盘强发育带宽度约为10km；③地形坡度是地震地质灾害发育的控制性因素之一，绝大部分的灾害集中在坡度20°～50°的范围内；④地震地质灾害与高程和微地貌具有很好的对应关系，大部分灾害发生在高程1500～2000m以下的河谷峡谷段，尤其是峡谷段的上部（即宽谷向峡谷的转折部位），单薄的山脊、孤立或多面临空的山体对地震波最为敏感，具有显著的放大效应，此部位崩塌滑坡最为发育；⑤不同岩性与地质灾害的发育虽没有显著对应关系，但决定了地质灾害的类型，通常情况下，滑坡多发生在软岩中，而硬岩中多发生的是崩塌［32-34］。黄润秋等还在对汶川地震触发的最大规模滑坡：大光包—红洞子沟巨型滑坡研究基础上，指出其形成机理和动力过程大致分为坡体震裂、松弛和解体→高速溃滑→震动堆积→二次抛射和碎屑流堆积四个阶段［35］。许强等通过对汶川地震斜坡崩滑灾害大量的现场调查和分析认为，本次地震诱发的大型崩滑灾害具有震裂溃屈、临空抛射和碎屑流化等独特的动力学特征［36］，为强震诱发斜坡崩滑灾害成因机理提供了新认识。

此外，伍法权［37］、付碧宏［38］、崔鹏［39］、唐川［40］、王运生［41］、赵纪生［42］、张永双［43］、祁生文［44］等学者都对汶川地震触发的斜坡崩滑地质灾害的动力学特征和形成机制做了较深入研究，取得了可喜的研究成果。

事实上在地震发生时，对斜坡体起破坏作用的纵、横波由于传输速度的差异导致其到达特定斜坡体存在时间差，且由于纵、横波对穿越介质力学作用（纵波从断裂处以同等速度向所有方向外传，交替地挤压和拉张其穿过的岩土体，其颗粒在波传播的方向上向前和向后运动，即这些颗粒的运动是垂直于波前的；而横波使岩土体颗粒的运动横过运移方向，交替地剪切其穿过的介质）的差异，导致其对斜坡岩土体的破坏效应存在差别，即纵波周期拉压其穿越的斜坡岩土体，而横波周期剪切其穿越的斜坡岩土体。但令人遗憾的是，在分析地震作用下斜坡损伤和崩滑破坏成灾机理时，已有的研究鲜有从纵、横波时差（由于传输速度差异导致的到达特性斜坡体存在时间差）耦合作用（当其都到达斜坡体后，纵波的周期拉压破坏和横波的周期剪切耦合作用于斜坡体）角度入手的，而该角度正是比较或是最符合地震作用触发斜坡损伤和崩滑效应实际的，所以这无疑是该研究较具前景的发展方向之一。

3 地震作用下斜坡稳定性评价及致灾预测研究

3.1 地震作用下斜坡稳定性评价方法研究

安全系数法是进行地震作用下斜坡稳定性评价的一个较常规方法。包括两种，一种是极限平衡法，即把计算得到的滑动面的抗滑力与滑动力之比作为边坡的安全系数，其中根据合力计算方法的差异又分为几种不同的方法，如 Bishop 法［45］、Janbu法［46］、Morgenstern-Price法［47］、Sarma法［48］等。另一种是强度折减法，其是在外荷载不变的条件下逐渐降低岩土体的抗剪强度，直到坡体达到破坏或临界稳定状态，定义安全系数为岩土体的实际剪切强度与临界状态时对应的折减后强度比值。随后，国内外很多学者利用强度折减法对边坡的安全系数进行了分析，如：Naylor［49］、Donald ＆ Giam［50］、Ugai［51-52］、Matsui et al［53-54］、郑颖人［55-56］、葛修润［57］、郑宏［58］。

目前，评价地震作用下斜坡稳定性的方法还有物理模型试验、数值模拟分析等。

物理模型试验主要包括拟静力试验、拟动力试验和地震模拟振动台试验，它们都是实验室常用的测试结构抗震性能的有效手段。拟静力试验以较低的加载速率实现单调或周期的加载，可以最大限度地获得构件刚度、承载力、变形和耗能等信息；拟动力试验是将计算机的计算和结构试验相结合的一种试验方法，最大的优点是结构的恢复力特性不是来自数学模型，而是直接从被测结构上实时取得，但是试验的加载过程还是拟静力的；地震模拟振动台试验是真正意义上的地震模拟试验，台面上可以真实地再现各种形式的地震波，结构在地震作用下的破坏机理也可以直观的被了解，是目前研究结构抗震性能最直接也是较准确的试验方法［59］。早在20世纪40年代，地震模拟振动台首次在土木结构上被用来模拟地震作用［60］，60年代以后地震模拟振动台开始广泛被建设［61］。

近年来，振动台试验在岩质边坡动力响应特征及稳定性评价研究上得到了广泛应用。王思敬等较早地通过振动模拟试验及相关理论推导对上述问题进行了研究，得出了以下主要成果：①应用振动台和加速度测量系统，测得滑动面在块体运动时的动摩擦系数及其与运动速度的关系曲线；②建立了边坡块体运动的动力微分方程，通过数值积分求得块体滑动的动力学特征，即块体运动加速度和块体相对基岩的运动加速度、运动速度和位移曲线；③进行了层状边坡岩体滑动稳定的三维动力学分析，建立了块体运动的三维动力学方程，并探讨了层状边坡层间振动加速度在传递过程中的加速度变化情况，从而可以计算出在不同震源方向作用下边坡各层间的稳定性［62-65］。王存玉在二滩拱坝动力模型试验中发现，岩石边坡对地震加速度不仅存在铅直向的放大作用，而且还存在水平向的放大作用［66］。徐光兴等用大型振动台模型试验，探讨了地震作用下土体模型边坡的动力特性与动力响应规律，以及地震动参数对动力特性和动力响应的影响，结果表明：边坡对输入地震波存在垂直放大和临空面放大作用，边坡土体对输入地震波低频部分存在放大作用，对高频部分存在滤波作用。坡面加速度峰值放大系数随输入地震波振幅、频率的增加而减小，持时对加速度峰值响应的影响不大；坡体位移随振幅、持时的增加而显著增大，随频率的增大而减小［67-68］。土耳其学者 Kumsar H等人对楔形体的动力稳定性，利用模型试验和理论相结合的方法进行过一些研究［69］。Lin M L在岩土边坡的动力特性和动力响应研究方面，也做了类似的工作［70］。此外，很多学者［71-73］都进行了节理岩体的动力特性试验研究，深化了对于节理岩体动力特性的认识。

此外，随着计算机技术发展的突飞猛进，多种数值模拟技术应运而生并成功应用于斜坡动力分析。相比物理模拟试验，数值模拟技术具有简便、易行、快速、可行性强等特点，在一定程度上弥补了物理模拟试验耗时、耗财、耗力等缺陷，因而具有良好的发展前景。目前主要的数值模拟技术有：有限单元法、有限差分法、离散单元法、拉格朗日元法、非连续变形分析方法［74］、流形元法［75］、边界元法、无界元法以及几种半解析元法［76］。

数值方法引入边坡的动力响应分析始于20世纪60年代，Clough R W［77］和 Chopra A K［78］第一次把有限单元法引入分析坝坡的地震反应问题。随后，两例有限差分研究［79-80］较早地注意了无限长坝坡在垂直入射S波作用下的动力反应问题。20世纪70年代以来，有限单元法已经被广泛地应用于土体的动力反应分析中。目前动力有限单元法的发展趋势是应用“真非线性”的应力—应变关系，进而引入反映土体弹塑性、流变性质的动力本构关系。由于有限单元法不但可以采用总应力法，而且可以采用有效应力法，可以考虑复杂地形、土的非线性、非均质性、弹塑性及土中孔隙水等因素的影响，能够深入分析土的自振特性及土体各部分的动力反应，因此有限单元法已经成为土体动力分析中最重要的分析方法［81］。差分法的基本思想是用差分网格离散求解域，用差分公式将科学问题的控制方程（常微分方程或偏微分方程）转化为差分方程，然后结合初始条件和边界条件，求解线形方程组。由于这种方法比较直观，容易编制程序，所以从20世纪40年代以来，至今仍被广泛采用［76］。有限差分法是最早被引入研究边坡动力响应的数值方法。后来由于有限元法的灵活性以及边界元法的异军突起，使得有限差分法在土工动力分析中（包括静力分析）趋于停滞。然而在最近几年，差分法又有了新进展。任意网格的差分，使这个老方法又可以和有限元法相匹敌。

自20世纪70年代开始，国内外研究人员提出了多种考虑结构面特征以不连续变形分析方法为基础的数值分析方法，如离散单元法（DEM）、非连续变形分析法（DDA）和刚块—弹簧法（RBSM）等。离散单元法是Cundall P A［82］提出的，于20世纪80年代引入我国［83］。该法特别适合运用于节理岩体的应力应变分析，在采矿、隧道、边坡及基础工程方面均有重要应用。Cundall发表了三维动态松弛离散元基本算法［84-85］，由ITASCA咨询集团开发了三维商业化程序3DEC（Three-Dimensional Element Code）。3DEC成为目前岩土力学领域中功能最全也最强大的三维非连续体程序之一，在采矿、地下核废料处理、冰川力学等诸多方面取得了巨大的成功。国内在离散元研究上起步较晚，但发展迅速。

在将以上数值计算方法应用于边坡动力稳定性分析方面，国内外做了大量的工作。何蕴龙在大量有限元动力计算的基础上，分析了地震荷载作用下边坡动力系数的分布规律以及与坡高、坡角以及坡体弹模之间的关系［86］。毛彦龙论证了地震导致斜坡岩土体变形破坏的启移动力学过程。研究认为，地震时的坡体波动振荡在斜坡岩土体变形破坏过程中产生三种效应：累进破坏效应、启动效应和启程加速效应。并采用离散元对斜坡在地震作用下的变形破坏过程进行了数值模拟分析［87-88］。彭定采用有限元 ADINA分析了地震荷载作用下三峡船闸高边坡的位移与应力分布情况，并通过SARMA法分析了稳定性系数［89］。徐卫亚采用有限元 FINAL分析了地震荷载作用下水布垭马崖高陡边坡的应力分布情况［90］。蔡德所应用非线性显式有限元程序DYNA3D，对三峡工程左岸厂房坝段在地震动力作用下的应力、位移进行了研究［91］。陶连金对节理岩体边坡的动力稳定性进行了离散元模拟，并推导了离散元的计算公式［92］；赵坚用 UDEC模拟了爆炸波在节理岩体中的传播［93］；张丙印采用离散元方法对某边坡的稳定性进行了数值模拟［94］，崔芳鹏［95-99］等对具地域性和空间非均质性的地震纵横波时差耦合作用导致的斜坡体初期的崩滑破坏及后期的抛射、堆积过程进行了离散元模拟。

国外方面，Lorig在博士论文中将离散元法与边界元法进行了耦合［100］。Wang B L与Ghaboussi在静态松弛的离散单元法的基础上提出了他们的刚块—弹簧法模型［101-102］，编制了 Fortran程序 BOLSMER（Block Spring Modelfor Excavation in Rocks）。Devin P等做了岩坡倾倒破坏实验［103］，并与离散元模拟结果进行了对比，发现二者并不完全吻合，认为其重要原因是实际接触点附近损伤区的存在和裂纹的扩展导致几何条件与计算模型产生差异。此外，Bouckovalas［104］、刘汉龙［105］等也采用数值模拟对地震作用下斜坡的动力响应特性做了较深入分析，使该项研究进一步走向深入。

3.2 地震动力作用下斜坡崩滑致灾预测研究

地震作用下斜坡崩滑破坏的形成机制及动力学特征等的研究，其落脚点都对准了致灾预测这一难题。在致灾预测研究上，有通过统计学分析既有斜坡崩坏灾害得出其一般规律的，也有通过理论研究通过理论公式计算出斜坡崩滑运动的永久位移的，通过数值模拟技术对特定崩滑灾害研究得出其一般特征的研究也不在少数，当然，物理模型试验作为一种可靠手段也被广泛地应用在此研究中。

Jibson等以1994年美国加利福尼亚的Northridge地震（Ms=6.7）为实例，对地震附近一个区域制作了数字化地震滑坡灾害图，该图所显示地震滑坡危险区域与震后的滑坡实际情况吻合较好［106］。国内方面，张楚汉等［107］采用动力离散元研究了地震动输入对三峡船闸及高岩石边坡所产生永久位移的影响，并与拟静力法进行了对比。陈云敏等［108］结合 Newmark的有限滑动位移思想，在动力有限元计算的基础上，提出了基于圆弧滑动的永久位移的计算方法。

在地震崩滑灾害的致灾预测上，鲜有见到成功的实例，究其原因一方面是大多数斜坡体主客观资料获取存在难度，另一方面是现有的技术手段较难实现斜坡崩滑距离、影响范围等因素在内的预测内容，故该研究还有待于进一步继续深入。

4 地震动力作用触发斜坡崩滑响应研究展望

目前，国内外学者对地震诱发的滑坡、崩塌地质灾害已有了较深入和细致的研究并取得了可喜成果，但依然存在不少问题。就研究内容而言，对于地震诱发斜坡崩滑灾害的形成机制及主控触发因素仍需深入研究；就研究对象而言，由于受到试验条件等因素的限制，岩质边坡的研究仍显薄弱，而对于土质边坡的动力响应问题研究则相对较多；此外，多数研究集中在对斜坡本身动力稳定性研究上，但对边坡输入的动力条件研究却较少，在进行动力条件输入时没有考虑远场的地质环境对近场动荷载施加的影响（主要表现为地震波的衰减及由此引起的地震峰值加速度比的变化趋势等），这与地震动力条件地面监测资料的有限性有较大关系；最后，针对地震动力条件下的斜坡崩滑早期识辨技术鲜有见到，因此包括斜坡坡度、岩土体风化程度及动力条件的地域性和空间非均质性特征等因素在内的斜坡体崩滑早期识辨技术体系的建立便具有重要的理论和现实意义。

［1］祁生文，伍法权，严福章，等.岩质边坡动力反应分析［M］.北京：科学出版社，2007.QI Shengwen，WU Faquan，YAN Fuzhang，et al.Rock slope dynamic response analysis［M］.Beijing：Science Press，2007.

［2］殷跃平.汶川八级地震地质灾害研究［J］.工程地质学报，2008，16（4）：433-444.YIN Yueping.Research on the geo-hazards triggered by Wenchuan earthquake，Sichuan ［J］.Journal of Engineering Geology，2008，16（4）：433-444.

［3］胡聿贤.地震工程学［M］.北京：地震出版社，1988.HU Yuxian.Earthquake engineering［M］.Beijing：Seismological Publishing House，1988.

［4］郑颖人，叶海林，黄润秋.地震边坡破坏机制及其破裂面的分析探讨［J］.岩石力学与工程学报，2009，28（8）：1714-1723.ZHENG Yingren，YE Hailin，HUANG Runqiu.Analysis and discussion of failure mechanism and fracture surface of slope under earthquake［J］.Chinese Journal of Rock Mechanics and Engineering，2009，28（8）：1714-1723.

［5］殷跃平.汶川地震地质与滑坡灾害概论［M］.北京：地质出版社，2009.YIN Yueping.Introduction of seismic geology and landslide hazard triggered by Wenchuan earthquake［M］.Beijing：Geological Publishing House，2009.

［6］Keefer D K and Wilson R C.Predicting earthquakeinduced landslides with emphasis on arid and semi-arid environments［A］.in Sadler P M，and Morton D M.eds.Landslides in a semi-arid environment with emphasis on the inland Valleys of Southern California［C］.California：Publications of the Inland Geological Society，1989，2：118-149.

［7］Richter C F.Elementary seismology［M］.San Francisco：W H.Freeman and Co.，1958.

［8］Prestininzi A and Romeo R.Earthquake—induced ground failures in Italy［J］.Engineering Geology，2000，58（3-4）：387-397.

［9］Papadopoulos G A，and Plessa A.Magnitude-distance relations for earthquake-induced landslide in Greece［J］.Engineering Geology，2000，58（3-4）：377-386.

［10］Keefer D K.Statistical analysis of an earthquake-induce landslide distribution-the 1989 Loma Prieta，California event［J］.Engineering Geology，2000，58（3-4）：231-249.

［11］Rodriguez C E，Bommer JJ，Chandler R J.Earthquakeinduced landslides：1980-1997［J］.Soil Dynamics and Earthquake Engineering，1999，18：325-346.

［12］Mario Parise，Randall W，Jibson.A seismic landslide susceptibility rating of geologic units based on analysis of characteristics of landslides triggered by the 17 January，1994 Northridge，California earthquake［J］.EngineeringGeology，2000，58：251-270.

［13］Julian J B，Carlos E，Rodriguez.Earthquake-induced landslides in Central America［J］.Engineering Geology，2002，63：189-220.

［14］周本刚，王裕明.中国西南地区地震滑坡的基本特征［J］.西北地震学报，1994，16（1）：95-103.ZHOU Bengang，WANG Yuming.Some characteristics of earthquake-induced landslide in Southwestern China.Northwestern Seismological Journal，1994，16（1）：95-103.

［15］薛守义.岩体边坡动力性研究［D］.北京：中国科技大学，1989.XUE Shouyi.Study on dynamic response of rock slope［D］.Beijing：University of Science and Technology of China，1989.

［16］李天池.地震与滑坡的关系及地震滑坡预测探讨［A］.滑坡文集（第二集）［C］.北京：中国铁道出版社，1979，127-132.LI Tianchi.Discussion on relationship between earthquake and landslide and earthquake-induced landslide predication［A］.Landslide Collected Works（Volume No.2）［C］.Beijing： China Railway Publishing House，1979，127-132.

［17］辛鸿博，王余庆.岩土边坡地震崩滑及其初判准则［J］.岩土工程学报，1999，21（5）：591-594.XIN Hongbo，WANG Yuqing.Earthquake-induced landslide and avalanche［J］.Chinese Journal of Geotechnical Engineering，1999，21（5）：591-594.

［18］孙崇绍，蔡红卫.我国历史地震时滑坡崩塌的发育及分布特征［J］.自然灾害学报，1997，6（1）：25-30.SUN Chongshao，CAI Hongwei.Developing and distributing characteristics of collapses and landslides during strong historic earthquake in China［J］.Journal of Natural Disasters，1997，6（1）：25-30.

［19］李忠生.国内外地震滑坡灾害研究综述［J］.灾害学，2003，18（4）：64-70.LIZhongsheng.The state of the art of the research on seismic landslide hazard at home and abroad［J］.Journal of Catastrophology，2003，18（4）：64-70.

［20］尹光华，李军，张勇，等.尼勒克地震滑坡的统计分析及初步研究［J］.内陆地震，2001，15（1）：56-63.YIN Guanghua，LI Jun，ZHANG Yong，et al.Statistical analysis and research of landslides caused by Nileke Earthquake［J］.Inland Earthquake，2001，15（1）：56-63.

［21］丁彦慧，王余庆，孙进忠，等.地震崩滑预测方法及其工程应用研究［J］.工程地质学报，2000，8（4）：475-480.DING Yanhui，WANG Yuqing，SUN Jinzhong，et al.Research on the method for prediction of earthquake induced landslides and its application to engineering projects［J］.Journal of Engineering Geology，2000，8（4）：475-480.

［22］王余庆，辛鸿博，高艳平，等.预测岩土边坡地震崩滑的综合指标法研究［J］.岩土工程学报，2001，23（3）：311-314.WANG Yuqing，XIN Hongbo，GAO Yanping，et al.Study on comprehensive index method for predicting earthquake-induced landslides［J］.Chinese Journal of Geotechnical Engineering，2001，23（3）：311-314.

［23］杨涛，邓荣贵，刘小丽，等.四川地区地震崩塌滑坡的基本特征及危险性分区［J］.山地学报，2002，20（4）：456-460.YANG Tao，DENG Ronggui，LIU Xiaoli，et al.The distributing and subarea character of the seismic landslides in Sichuan［J］.Journal of Mountain Science，2002，20（4）：456-460.

［24］崔芳鹏，胡瑞林，张明，等.重庆小南海地震崩滑体的基本特征及形成机制研究［J］.工程地质学报，2009，17（2）：167-174.CUI Fangpeng，HU Ruilin，ZHANG Ming，et al.Basic characteristics and formation mechanism of collapsed and slided debris mass induced by xiaonanhai earth-quake in Chongqing city［J］.Journal of Engineering Geology，2009，17（2）：167-174.

［25］周维垣.高等岩石力学［M］.北京：水利电力出版社，1990.ZHOU Weiyuan.Higher rock mechanics［M］.Beijing：Water Resources and Electric Power Publishing，1990.

［26］张倬元，王士天，王兰生.工程地质分析原理［M］.北京：地质出版社，1993.ZHANG Zhuoyuan，WANG Shitian，WANG Lansheng.Engineering geology analysis principle［M］.Beijing：Geological Publishing House，1993.

［27］胡光韬.滑坡动力学［M］.北京：地质出版社，1995.HU Guangtao.Dynamic landslide［M］.Beijing：Geological Publishing House，1995.

［28］胡广韬，毛彦龙，程谦恭，等.地震时滑坡体波动震荡的启程加速动力学问题［R］.西安：西安工程学院，1997.HU Guangtao，MAO Yanlong，CHENG Qiangong，et al.Set-out violent-slide dynamic problems of slopeshaking during an earthquake［R］.Xi'an： Xi'an Engineering University，1997.

［29］毛彦龙，胡广韬，毛新虎，等.地震滑坡启程剧动的机理研究及离散元模拟［J］.工程地质报，2001，9（1）：74-80.MAO Yanlong，HU Guangtao，MAO Xinhu.et al.Mechanism of set-out violent-slide of slope mass during earthquake and its simulation by using discrete element method［J］.Journal of Engineering Geology，2001，9（1）：74-80.

［30］殷跃平.汶川八级地震滑坡特征分析［J］.工程地质学报，2009，17（1）：29-38.YIN Yueping.Fearures of landslides triggered by the Wenchuan earthquake［J］.Journal of Engineering Geology，2009，17（1）：29-38.

［31］殷跃平.汶川八级地震滑坡高速远程特征分析［J］.工程地质学报，2009，17（2）：153-166.YIN Yueping.Rapid and long run-out features of landslides triggered by the Wenchuan Earthquake［J］.Journal of Engineering Geology，2009，17（2）：153-166.

［32］黄润秋，李为乐.“5.12”汶川大地震触发地质灾害的发育分布规律研究［J］.岩石力学与工程学报，2008，17（12）：2585-2592.HUANG Runqiu，LIWeile.Research on development and distribution rules of geohazards induced by wenchuan earthquake on 12th may，2008［J］.Chinese Journal of Rock Mechanics and Engineering，2008，17（12）：2585-2592.

［33］黄润秋，李为乐.汶川大地震触发地质灾害的断层效应分析［J］.工程地质学报，2009，17（1）：19-28.HUANG Runqiu，LIWeile.Fault effect analysis of geohazard triggered by the Wenchuan Earthquake［J］.Journal of Engineering Geology，2009，17（1）：19-28.

［34］黄润秋，许强，等.中国典型灾难性滑坡［M］.北京：科学出版社，2008.HUANG Runqiu，XU Qiang.Catastrophic landslides in China［M］.Beijing：Science Press，2008.

［35］黄润秋，裴向军，李天斌.汶川地震触发大光包巨型滑坡基本特征及形成机理分析［J］.工程地质学报，2008，16（6）：730-741.HUANG Runqiu，PEI Xiangjun，LI Tianbin，et al.Basic characteristics and formation mechanism of the largest scale landslide at Daguangbao occurred during the Wenchuan Earthquake［J］.Journal of Engineering Geology，2008，16（6）：730-741.

［36］许强，黄润秋.5.12汶川大地震诱发大型崩滑灾害动力特征初探［J］.工程地质学报，2008，16（6）：721-729.XU Qiang，HUANG Runqiu.Kinetics charateristics of large landlides triggered by may 12th Wenchuan Earthquake［J］.Journal of Engineering Geology，2008，16（6）：721-729.

［37］伍法权，胡瑞林，岳中琦，等.5.12汶川地震地质灾害［M］.北京：地质出版社，2009.WU Faquan，HU Ruilin，YUE Zhongqi，et al.Geohazard triggered by the Wenchuan Earthquake［M］.Beijing：Geological Publishing House，2009.

［38］付碧宏，王 萍，孔 屏，等.中国汶川“5.12”8.0级大地震地震地质灾害图集［M］.北京：地震出版社，2009.FU Bihong，WANG Ping，KONG Ping，et al.Atlas of geo-hazard triggered by the Wenchuan Earthquake［M］.Beijing：Seismological Publishing House，2009.

［39］崔 鹏，韦方强，何思明，等.5.12汶川地震诱发的山地灾害及减灾措施［J］.山地学报，2008，26（3）：280-282.CUI Peng，WEI Fangqiang，HE Siming，et al.Mountain Disasters Induced by the Earthquake of May 12 in Wenchuan and the Disasters Mitigation［J］.Journal of Mountain Science，2008，26（3）：280-282.

［40］唐 川，梁京涛.汶川震区北川9.24暴雨泥石流特征研究［J］.工程地质学报，2008，16（6）：751-758.TANG Chuan，LIANG Jingtao.Characteristics of debris flows in Beichuan epicenter of the Wenchuan Earthquake triggered by rainstorm on Septermber 24，2008［J］.Journal of Engineering Geology，2008，16（6）：751-758.

［41］王运生，罗永红，吉峰，等.汶川大地震山地灾害发育的控制因素分析［J］.工程地质学报，2008，16（6）：759-763.WANG Yunsheng，LUO Yonghong，JI Feng，et al.Analysis of the controlling factors on geo-hazard in mountainous epicenter zones of the Wenchuan Earthquake［J］.Journal of Engineering Geology，2008，16（6）：759-763.

［42］赵纪生，魏景芝，吴景发，等.汶川8.0级地震滑坡、崩塌机制［J］.震灾防御技术，2008，3（4）：379-383.ZHAO Jisheng，WEI Jingzhi，WU Jingfa，et al.Mechanism of landslide and rockfall caused by Wenchuan Ms8.0 Earthquake［J］.Technology for Earthquake Disaster Prevention，2008，3（4）：379-383.

［43］张永双，雷伟志，石菊松，等.四川5.12地震次生地质灾害的基本特征初析［J］.地质力学学报，2008，14（2）：109-116.ZHANG Yongshuang，LEIWeizhi，SHI Jusong，et al.General characteristics of 5.12 earthquake-induced geohazards in Sichuan［J］.Journal of Geomechanics，2008，14（2）：109-116.

［44］祁生文，许 强，刘春玲，等.汶川地震极重灾区地质背景及次生斜坡灾害空间发育规律［J］.工程地质学报，2009，17（1）：39-49.QIShengwen，XU Qiang，LIU Chunling，et al.Slope instabilities in the severest disaster areas of 5·12 Wenchuan Earthquake［J］.Journal of Engineering Geology，2009，17（1）：39-49.

［45］Bishop A W.The Use of the Slip Circle in the Stability Analysis of Slopes［J］.Geotechnique，1955，5（1）：7-17.

［46］Janbu N.Slope Stability Computations［J］.Embankment Dam Engineering，New York：John Wiley and Sons，1973：47-86.

［47］Morgenstern N R，Price V E.The Analysis of the Stability of General Slip Surfaces［J］.Geotechnique，1965，15（1）：79-93.

［48］Sarma S K.Stability Analysis of Embankments and Slopes［J］.Journal of the Geotechnical Engineering Division，ASCE，1979，105：1511-1524.

［49］Naylor D J.Finite elements and slope stability［A］.Numerical method in geo-mechanics，Proc.NATO advanced study institute［C］，Lisbon，Portugal，1981：229-244.

［50］Donald I B，Giam S K.Application of the nodal displacement method to slope stability analysis［A］.Proc.Fifth Australia-New Zealand Conf.on Geomech［C］.Sydney，Australia，1988：456-460.

［51］Ugai K.A method of calculation of total factor of safety of slope by elasto-plastic FEM ［J］.Soils and Foundations，1989，29（2）：190-195.

［52］Ugai K，Leshchinsky D.Three-dimensional limit equilibrium and finiteelement analysis：a comparison of results［J］.Soils and Foundations，1995，35（4）：1-7.

［53］Matsui T，San K C.Slope stability analysis by shear strength reduction technique［A］.Proceeding of the International Conference on Computer Methods and Advances in Geomechanics［C］，Balkema，Rotterdam，1991：499-504.

［54］Matsui T，San K C.Finite element slope stability analysis by shear strength reduction technique［J］.Soils and Foundation，1992，32（1）：59-70.

［55］郑颖人，赵尚毅.有限元强度折减法在土坡与岩坡中的应用［J］.岩石力学与工程学报，2004，23（19）：3381-3388.ZHENG Yingren，ZHAO Shangyi.Application of strength reduction FEM in soil and rock slope［J］.Chinese Journal of Rock Mechanics and Engineering，2004，23（19）：3381-3388.

［56］郑颖人，赵尚毅，张鲁渝.用有限元强度折减法进行边坡稳定分析［J］.中国工程科学，2002，4（10）：57-61，78.ZHENG Yingren，ZHAO Shangyi，ZHANG Luyu.Slope stability analysis by strength reduction FEM ［J］.Engineering Science，2002，4（10）：57-61，78.

［57］葛修润，任建喜，李春光，等.三峡左厂3#坝段深层抗滑稳定三维有限元分析［J］.岩土工程学报，2003，25（4）：389-394.GE Xiurun，REN Jianxi，LI Chunguang，et al.3DFEM analysis of deep sliding stability of 3# dam foundation of left power house of the Three Gorges［J］.Chinese Journal of Geotechnical Engineering，2003，25（4）：389-394.

［58］郑宏，李春光，李焯芬，等.求解安全系数的有限元法［J］.岩土工程学报，2002，24（5）：626-628.ZHENG Hong，LIChunguang，LIZhuofen，et al.Finite element method for solving the fator of safety［J］.Chinese of Geotechnical Engineering，2002，24（5）：626-628.

［59］王燕华，程文瀼，陆 飞，等.地震模拟振动台的发展［J］.工程抗震与加固改造，2007，29（5）：53-56.WANG Yanhua，CHENG Wenrang，LU Fei，et al.Development of the shaking table［J］.Earthquake Resistant Engineering and Retrofitting，2007，29（5）：53-56.

［60］Williams M S，Blakeborough A.Laboratory testing of structures under dynamic loads：an introductory review［J］.Phil.Trans R.Soc.Lond.A，2001，359（1786）：1651-1669.

［61］黄浩华.地震模拟振动台的设计与应用技术［M］.北京：地震出版社，2008.HUANG Haohua.Design and application of the shaking table［M］.Beijing：Seismological Publishing House，2008.

［62］王思敬，张菊明.岩体结构稳定性的块体力学分析［J］.地质科学，1980，1：19-33.WANG Sijing，ZHANG Juming.Block mechanical analysis of structure stability for rock mass［J］.Geology Science，1980，1：19-33.

［63］王思敬，张菊明.边坡岩体滑动稳定的动力学分析［J］.地质科学，1982，2：162-170.WANG Sijing，ZHANG Juming.Kinetics analysis of slope rock mass stability［J］.Geology Science，1982，2：162-170.

［64］Wang S J，Xue S Y，Maugeri M，et al.Dynamic stability of the left abutment in the Xiaolangdi project on the Yellow River［A］.In：Proc.Int.Symp.on Assessment and Prevention of Failure Phenomena in Rock Engineering［C］.Ankara：1993，5-7.

［65］Wang S J，Zhang JM.On the dynamic stability of block sliding on rock slopes［A］.In：Proc.Int.Cont.on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics［C］.St.Louis：1984，431-434.

［66］王存玉.地震条件下二滩水库库坡稳定性研究［M］.北京：科学出版社，1987.WANG Cunyu.Study on seismic stability of Ertan reservoir slope［M］.Beijing：Science Press，1987.

［67］徐光兴，姚令侃，高召宁，等.边坡动力特性与动力响应的大型振动台模型试验研究［J］.岩石力学与工程学报，2008，27（3）：624-632.XU Guangxing，YAO Lingkan，GAO Zhaoning，et al.Large-scale shaking table model test study of dynamic characteristics and dynamic response of slope［J］.Chinese Journal of Rock Mechanics and Engineering，2008，27（3）：624-632.

［68］徐光兴，姚令侃，李朝红，等.边坡地震动力响应规律及地震动参数影响研究［J］.岩土工程学报，2008，30（6）：918-923.XU Guangxing，YAO Lingkan，LI Chaohong，et al.Dynamic response of slopes under earthquakes and influence of ground motion parameters［J］.Chinese Journal of Geotechnical Engineering，2008，30（6）：918-923.

［69］Kumsar H，Aydan O，Ulusay R.Dynamic and static stability assessment of rock slopes against wedge failures［J］.Rock Mechanics and Rock Engrg，2000，33（1）：31-51.

［70］Lin M L，Wang K L.Seismic slope behavior in a largescale shaking table model test［J］.Engineering Geology，2006，86（2/3）：118-133.

［71］Crawford A W，Curran J H.The influence of shear velocity on the frictional resistance of rock discontinuities［J］.Int JRock Mech Min Sci＆ Geomech Abstr，1981（18）：505-515.

［72］Fox D J，Kana D D，Hsiung SM.Influence of interface roughness on dynamic shear behavior in jointed rock［J］.Int JRock Mech and M in Sci，1998，35（7）：923-940.

［73］Yoshinaka R，Tran T V，Osada M.Non-linear，stress and strain dependent behavior of soft rocks under cyclic triaxial conditions［J］.Int J Rock Mech and Min Sci，1998，35（7）：941-955.

［74］石根华.数值流形方法与非连续变形分析［M］.北京：清华大学出版社，1993.SHI Genhua.Numerical manifold method and discontinuous deformation analysis［M］.Beijing：Tsinghua University Press，1993.

［75］裴觉民，石根华.岩石滑坡体的动态稳定和非连续变形分析［J］.水利学报，1993，24（3）：28-34.PEI Juemin，SHIGenhua.Analysis of dynamic stability and discontinuous deformation of rock slope［J］.Journal of Hydraulic Engineering，1993，24（3）：28-34.

［76］龚晓南.土工计算机分析［M］.北京：中国建筑工业出版社，2000.GONG Xiaonan.Geotechnical computer analysis［M］.Beijing：China Architecture and Building Press，2000.

［77］Clough R W and Chopra A K.Earthquake stress analysis in earth dams［J］.J.Engrg.Mech.，ASCE 1966，92（EM 2）：51-60.

［78］Chopra A K.Earthquake response of earth dams［J］.J.Soil Mech and Found Div，ASCE 1967，93（SM 2）：66-82.

［79］Ishizaki Hatekeyama.Consideration on the dynamical behavior of earth dams［R］.Bul.No.52，Disaster Prevention Research Inst.，Kyoto Unvi.，1963.

［80］Medevedev S and Sinitsym A.Seismic effects on earth fill dams［A］.In：Van Roekel J.H.eds.Proc.3rdworld Conf.Earthq.Engrg.［C］，R.E.Owen government printer，Wellington，New Zealand，1965，Paper III：373-390.

［81］肖克强.地震荷载作用下顺层岩体边坡变形特征及稳定性研究［D］.武汉：中国科学院武汉岩土力学研究所，2006.XIAO Keqiang.Study on deformation characteristics and stability of rock bedding slopes under earthquake［D］.Wuhan：Institute of Rock and Soil Mechanics，Chinese Academy of Sciences，2006.

［82］Cundall P A.A computer model for simulating progressively large scale moments in blocky rock system［A］.In：Proceedings of the International Symposium on Rock Mech［C］.ISRM，Nancy，1971，II-8：129-136.

［83］王泳嘉.离散单元法—一种适用于节理岩石力学分析的数值方法［A］.第一届全国岩石力学数值计算及模型试验讨论会论文集［C］.成都：西南交通大学出版社，1986.WANG Yongjia.DEM-a numerical method for analysis of jointed rock mass mechanics［A］.Chengdu：Southwest Jiaotong University Press，1986.

［84］Cundall P A.Formulation of a three-dimensional distinct element model-Part I.A scheme to detect and represent contacts in system composed of many polyhedral blocks［J］.Int.J.Rock Mech.Sci.＆ Geotech.，Abstr.，1988，25（3）：107-116.

［85］Hart R，Cundall P A＆Lemos J.Formulation of a threedimensional distinct element model-Part II.Mechanical calculation for motion and interaction of a system composed of many polyhedral blocks［J］.Int.J.Rock Mech.M in.Sci.＆ Geotech.Abstr.1988，25（3）：107-116.

［86］何蕴龙，陆述远，段亚辉.岩石边坡地震作用计算方法研究［J］.长江科学院院报，1998，15（4）：35-38.HE Yunlong，LU Shuyuan，DUAN Yahui.Study on calculation method of seismic effect on rock slope［J］.Journal of Yangtze River Scientific Research Institute，1998，15（4）：35-38.

［87］毛彦龙，胡广韬，赵法锁，等.地震导致斜坡变形破坏的启程动力学过程初探［J］.地下空间，1999，19（5）：416-419.MAO Yanlong，HU Guangtao，ZHAO Fasuo，et al.Study on set-out kinetics process of landslide under earthquake［J］.Underground Space，1999，19（5）：416-419.

［88］毛彦龙，胡广韬，赵法锁，等.地震触发滑坡体滑动的机理［J］.西安工程学院学报，1998，20（4）：45-48.MAO Yanlong，HU Guangtao，ZHAO Fasuo，et al.Formation mechanism of landslide under earthquake［J］.Journal of xi’an engineering institute，1998，20（4）：45-48.

［89］彭定，陈玲玲.三峡船闸高边坡在地震作用下的稳定性分析［J］.长江科学院院报，2001，18（6）：47-49.PENG Ding，Chen Lingling.Stability analysis on high slopes of TGP lock under seismic action［J］.Journal of Yangtze River Scientific Research Institute，2001，18（6）：47-49.

［90］徐卫亚，谢守益，徐瑞春，等.水布垭马崖高陡边坡三维数值模拟研究［J］.武汉水利电力大学（宜昌）学报，1998，20（3）：11-15.XU Weiya，XIE Shouyi，XU Ruichun，et al.Study on three dimensional numerical simulation of Maya high slope on Shuibuya［J］.Journal of Wuhan University of Hydraulic and Electrical Engineering，1998，20（3）：11-15.

［91］蔡德所.地震输入波形对三峡工程左岸坝段的动力影响研究［J］.岩石力学与工程学报，1998，17（2）：183-187.CAI Desuo.Effects of earthquake input waves on dynamic response of dam segment of the left bank of Three Gorges project［J］.Chinese Journal of Rock Mechanics and Engineering，1998，17（2）：183-187.

［92］陶连金.节理岩体边坡的动力稳定性分析［J］.工程地质学报，2001，9（1）：23-38.TAO Lianjin.Dynamic stability analysis of jointed rock slope［J］.Journal of Engineering Geology，2001，9（1）：23-38.

［93］赵坚.用 UDEC模拟爆炸波在节理岩体中的传播［J］.中国矿业大学学报，2002，31（2）：111-115.ZHAO Jian.Simulation of blast wave propagation in jointed rock mass using UDEC［J］.Journal of China University of Mining and Technology，2002，31（2）：111-115.

［94］张丙印.昌马水库枢纽工程右岸岩石边坡稳定性的离散元分析［J］.水力发电学报，2002，1：73-80.ZHANG Bingyin.Stability analysis of the rock slope of Changma Reservoir’s right bank with discrete element method［J］.Journal of Hydrauelectric Engineering，2002，1：73-80.

［95］崔芳鹏.地震P-S波时差耦合作用的斜坡崩滑效应研究［D］.北京：中国科学院地质与地球物理研究所，2009.CUI Fangpeng.Study of collapsing and sliding response of slope triggered by single and combined action with time difference of P and S seismic waves［Ph.D.Thesis］［D］.Beijing：Institute of Geology and Geophysics，Chinese Academy of Science，2009.

［96］崔芳鹏，胡瑞林，殷跃平，等.地震纵横波时差耦合作用的斜坡崩滑效应研究［J］.工程地质学报，2009，17（4）：455-462.CUI Fangpeng，HU Ruilin，YIN Yuepin，et al.Numerical simulation of slope collapsing and sliding due to single or combined actions of seismic p-and s-waves with time difference［J］.Journal of Engineering Geology，2009，17（4）：455-462.

［97］崔芳鹏，胡瑞林，殷跃平，等.纵横波时差耦合作用的斜坡崩滑效应离散元分析——以北川唐家山滑坡为例［J］.岩石力学与工程学报，2010，29（2）：319-327.CUI Fangpeng，HU Ruilin，YIN Yueping，et al.Discrete element analysis of collapsing and sliding response of slope triggered by time difference coupling effects of p and s seismic waves-taking Tangjiashan landslide in Beichuan county for example［J］.Chinese Journal of Rock Mechanics and Engineering，2010，29（2）：319-327.

［98］崔芳鹏，殷跃平，许 强，等.地震纵横波时差耦合作用的斜坡崩滑响应研究［J］.岩土工程学报，2010，32（8）：1266-1273.CUI Fangpeng，YIN Yueping，XU Qiang，et al.Collapsing and sliding response of slopes triggered by single and combined action with time difference of P and S seismic waves［J］.Chinese Journal of Geotechnical Engineering，2010，32（8）：1266-1273.

［99］崔芳鹏，胡瑞林，殷跃平，等.地震动力作用触发的斜坡崩滑高差效应研究［J］.水文地质工程地质，2010，37（4）：38-43.CUI Fangpeng，HU Ruilin，YIN Yueping，et al.Study on difference in elevation response of slope collapsing and sliding triggered by seismic action ［J］.Hydrogeology and Engineering Geology，2010，37（4）：38-43.

［100］Lorig L，Brady B H G，Cundall P A.Hybrid distinct element-boundary element analysis of jointed rock［J］.Int.J.Rock Mech.Min.Sci.＆ Geomech，Abstr.，1986，23（4）：303-312.

［101］Wang B L，Vinod K Garga.A numerical method for modeling large disp lacement of jointed rocks.I.Fundamentals［J］.Can Geotech J，1993，30：96-108.

［102］Brady B H G，Cramer M L，Hart R D.Technical note preliminary analysis of a loading test on a large basalt block［J］.Int.J.Rock Meck.Sc.＆ Geomech.，Abstr.，1985，22（5）：345-348.

［103］Devin P，Zaninetti A.Validation of a distinct element model for topping rock slopes［A］.7thInt.Cont.，On Rock Mech［C］，1995，1：417-421.

［104］Bouckovalas G D，Papadimitrioul A G.Numerical evaluation of slope topography effects on seismic ground motion［J］.Soil Dynamics and Earthquake Engineering，2005，25（7-10）：547-558.

［105］刘汉龙，费 康，高玉峰.边坡地震稳定时程分析方法［J］.岩土力学，2003，24（4）：553-556.LIU Hanlong，FEI Kang，GAO Yufeng.Time history analysismethod of slope seismic stability［J］.Rock and Soil Mechanics.2003，24（4）：553-556.

［106］Jibson R W，Harp E l and Michael J A.A method for producing digital probabilistic seismic landslide hazard maps［J］.Engineering Geology，2000，58（3-4）：271-289.

［107］Zhang Chunhan，Pekau O A，Jin Feng，et al.Application of distinct element method in dynamic analysis of high rock slopes and blocky structures［J］.Soil dynamics and earthquake engineering，1997，16（6）：385-394.

［108］陈云敏，柯瀚，凌道盛.城市垃圾填埋体的动力特性及地震响应［J］.土木工程学报，2002，35（3）：66-72.KE Yunmin，KE Han，LING Daosheng.Dynamic properties and seismic response of municipal solid waste［J］.China Civil Engineering Journal，2002，35（3）：66-72.

［109］祁生文，伍法权，刘春玲，等.地震边坡稳定性的工程地质分析［J］.岩石力学与工程学报，2004，23（16）：2792-2796.QI Shengwen，WU Faquan，LIU Chunling，et al.Engineering geology analysis on stability of slope under earthquake［J］.Chinese Journal of Rock Mechanics and Engineering，2004，23（16）：2792-2796.

［110］杜晓丽.地震荷载作用下岩质边坡稳定性研究［D］.西安：西安科技大学，2008.DU Xiaoli.Study on rock slope stability in earthquake load［D］.Xi’an：Xi’an University of Science and Technology，2008.