Discontinuities and Engineering Geological Significances of Strong Earthquake-Induced Daguangbao Landslide
-
摘要: 大光包滑坡是2008年汶川地震触发的最大规模滑坡,也是世界上百年罕见的大型滑坡,引起国内外学者持续争论和广泛研究. 本文为研究地质构造对大光包滑坡影响,基于现场照片、遥感图及岩体结构面产状数据,首先开展大光包滑坡断壁、滑带及滑坡堆积体不连续地质特征调查,进而分析该不连续地质特征成因,进一步讨论强地震动、不连续地质界面组合及地下水对大光包滑坡启动影响. 研究结果表明:震前随大水闸背斜构造演化发育的断层、层间错动带及优势岩体结构面等不连续地质结构控制了大光包滑坡边界的形成和滑体堆积过程;汶川强震为这些不连续地质结构破坏提供了动力条件;地下水伴随着强震过程不连续地质界面贯通促使了大光包滑坡快速启动.Abstract: The Daguangbao (DGB) landslide is the largest landslide triggered by the 2008 Wenchuan earthquake and considered as one of the largest of these events in the world. In this paper, the geological cause of the DGB landslide was analysed on the basis of rich photographs, remote sensing images, and orientation data. The discontinuity features of the landslide scarp, slip zone, and deposition were investigated firstly, and then the formations of these discontinuities were studied. The effects of strong seismic shaking, geological discontinuities, and ground water on the DGB landslide were discussed. The results showed that the initiation and deposition of the DGB landslide were dominated by geological discontinuities such as local faults, bedding faults, and joint sets, which were developed during the formation of an anticline. Strong ground shaking of the Wenchuan earthquake induced the failure of these geological discontinuities, and ground water might have finally initiated the DGB landslide during the earthquake-induced discontinuity connectivity.
-
Key words:
- rock structure /
- geological structure /
- Wenchuan earthquake /
- sliding zone /
- ground water
-
有砟轨道作为一种常用的轨道结构形式,散体道床是其重要组成部分,在长期服役过程中,不可避免会出现病害,其中道床板结是有砟轨道最主要的设备病害之一[1]. 道床板结是道砟颗粒破碎粉化产物、行车坠落及风吹降落物、路基层向上迁移的细颗粒等混合在道砟中,并在长期作用下形成的混合体. 道床板结后会降低轨道弹性,加剧列车的振动冲击作用,加速轨道结构的状态恶化,缩短有砟轨道线路的服役寿命,增大线路的养护维修工作量,严重影响铁路运营经济效益,相关问题值得深入探究.
目前,国内外针对道床板结问题已开展了一定的研究工作:在试验方面,郭善辛[2]对板结道床的成因进行了分析,并对道床的脏污率、孔隙率、变形模量等参数进行了测试,提出了道床板结的参考技术指标;赵成江等[3]开展了板结道床的室内冻胀试验,研究了不同含泥量与含水量下板结道床的冻胀特性;Fortunato等[4]开展了三轴试验指出,脏污道床的弹性模量受脏污干湿状态的影响;Tang等[5]通过三轴试验,分析了路基土侵入有砟道床形成的脏污板结道床集料的力学特性,当脏污板结道床的含水率减小时道床的抗剪强度会提高;Paiva等[6]进行了有砟道床的室内渗透率试验,指出随着道床中土壤含量的增加,道床的透水性迅速降低;Ishikawa等[7]开展了不同含水率条件下洁净道床和脏污道床的室内三轴试验,结果表明,土和水的共同作用会降低有砟道床的抗剪强度. 既有的关于板结道床的试验研究多采用现场原位试验和室内三轴试验,测试分析了板结道床在水等因素作用下的冻胀、抗剪和透水性能等。但由于试验测试大多只是对宏观现象的分析,对道床内部道砟颗粒细观角度下的力学特性的揭示方面存在不足.
在理论方面,金寿延等[8]采用二自由度钢轨振动理论定性地分析道床板结对钢轨振动的影响,结果表明,道床板结会增大钢轨和轨枕的振动,降低钢轨和轨枕的位移,减弱轨道系统的弹性功能;赵宁[9]基于车轨耦合理论和有限单元法建立动力学模型,通过改变材料刚度参数模拟道床板结,其指出道床板结会增强轨道结构的振动,会使轨道结构的位移显著降低. 既有的研究手段多采用简化理论和有限单元法,多通过直接改变材料参数模拟道床板结,对于道床板结的细观力学行为的表征方面存在不足,并且也无法有效地考虑有砟道床的散粒体特性.
综上,既有针对道床板结的研究大多没有从细观角度出发分析板结对有砟道床力学特性的影响,并且在考虑有砟道床的散体特性方面存在不足. 有砟道床由不同粒径的道砟颗粒堆积而成,具有散体特性. 板结道床是由于道砟颗粒之间的脏污材质相互胶结作用而形成的硬块,具有一定的整体性,但是在列车荷载作用下道砟颗粒之间的胶结作用可能被破坏,或者板结程度较低时道砟颗粒之间并没有完全粘接牢固,导致板结道床处于介于散粒体和连续体之间的状态,用有限单元法模拟散体道砟存在一定的不足,未有效考虑道砟颗粒及脏污板结材质的散体特性. 基于此,本文将采用离散单元法建立脏污道床的仿真分析模型分析列车荷载作用下道床板结对有砟道床动力行为的影响,并从细观角度揭示板结使道床弹性降低的内在原因,相关结论可为有砟道床的养护维修提供参考.
1. 模型建立
大规模离散元数值模拟对计算设备具有极高的要求,直接建立三维模型开展分析会消耗大量的计算资源. 考虑到实际中有砟道床主要承受列车垂向荷载的作用,可主要关注道砟颗粒沿竖向的运动情况,为此,可将铁路线路简化为二维问题进行数值模拟[10]. Zhang等[11]验证了有砟道床的二维离散元模型适用于有砟道床的动力特性分析. 另外,板结道床中的细小脏污颗粒数量巨大,会极大地增加仿真计算量,所以经比选后采用二维离散元模型进行仿真分析.
道砟颗粒形状对颗粒力学行为有重要影响[12],实现对道砟颗粒真实外形的模拟是仿真计算的关键所在. 本文首先利用数码相机对道砟颗粒进行图像采集,然后运用图像处理技术获得道砟颗粒的外轮廓信息. 并根据道砟颗粒的轮廓信息,基于离散元软件PFC用clump单元模拟道砟颗粒,如图1所示.
按照《铁路碎石道砟》(TB/T 2140—2018)[13]中特级碎石道砟要求生成散体道床模型. 首先通过Fish语言编程,设置每个粒径区间道砟质量的上下限,结合clump distribute命令和bin关键词控制道砟的级配,通过落雨法并配合分层压实最终得到密实的道床. 另外,参考已有文献[11, 14],钢轨和扣件采用ball单元模拟,在单元之间建立平行黏结接触模型,实现对钢轨的连续性模拟. 轨枕则采用clump单元模拟,轨枕及道砟颗粒之间采用线性接触模型.
线性接触模型的计算式[15]为
{Fn=KnUn,Fs=−∑ KsΔUs, (1) 式中:Fn为法向接触力,以压力为正;Kn为接触点的法向刚度;Un为单元之间的法向重叠量;Fs为切向接触力;Ks为接触点的切向接触刚度;ΔUs为单元之间切向位移增量.
当Fs达到最大摩擦力Fsmax=μ|Fn|(其中μ为摩擦系数)时,道砟颗粒之间发生滑动.
平行黏结的断裂准则[15]为
{Fn<−Fn[C],|Fs|>Fs[C], (2) 式中:Fn[C]为法向黏结强度;Fs[C]为切向黏结强度. 当单元之间的接触力满足式(2)中任一不等式时,黏结键断裂,平行黏结失效,单元之间的接触转化为线性接触.
建立有砟轨道离散元模型如图2所示. 将轨枕从左至右依次进行编号. 所建模型中道砟颗粒的clump数量为1 227个,内部共含151791个pebble单元. 模型全长为50.4 m,对中间的9跨轨枕建立完整的道床模型,为减小计算量,对其余区域进行简化仅建立钢轨和扣件模型,约束模型两端钢轨沿水平方向的位移.
板结道床是由于道砟的粉化、脏污与道砟间胶合物质的相互作用使道砟粘接成块而形成的. 从板结道床的微观结构出发,可将板结道床视为混有粘胶、具有一定黏结力的道砟体[16]. 本文将道砟颗粒视为块体单元,在洁净道床模型中道砟颗粒的空隙之间投放小颗粒单元模拟脏污介质,通过定义脏污颗粒之间的黏结作用模拟道床板结现象. 道砟颗粒之间采用线性接触模型进行模拟.
采用直径2 mm的圆盘单元模拟脏污颗粒. 板结程度近似采用文献[17]中的脏污指标进行评估,即认为对于脏污较严重的道床,在定义脏污颗粒之间的黏结作用后用于模拟板结较程度较严重的道床. 所建立的板结道床的局部模型如图3所示.
通过控制板结道床中脏污颗粒数量的多少,共设置4种计算工况,以分析道床板结以及道床板结程度对有砟道床动力特性的影响,如表1所示.
表 1 计算工况Table 1. Calculation conditions工况编号 道床 板结程度/% 脏污颗粒数量/粒 1 洁净道床 0 0 2 板结道床 25 17605 3 板结道床 50 35213 4 板结道床 100 70427 参考文献[16],离散元模型的计算参数如表2所示,忽略脏污颗粒的材质影响,取脏污颗粒的刚度参数与道砟颗粒一致.
表 2 离散元模型计算参数Table 2. Parameters of DE model参数 取值 道砟颗粒法向刚度/(N•m−1) 5.0 × 108 道砟颗粒切向刚度/(N•m−1) 5.0 × 108 道砟颗粒摩擦系数 0.55 道砟颗粒密度/(kg•m−3) 2 500 脏污颗粒密度/(kg•m−3) 2500 钢轨密度/(kg•m−3) 1029.5 轨枕密度/(kg•m−3) 3025.8 脏污颗粒黏结强度/MPa 0.5 黏结半径比例系数 1.0 2. 模型验证
开展现场动测试验以校验模型,选定某有砟道床试验段开展动态测试. 根据《轨道检测道床应变、振动测试方法》(TB/T 2488—2017)[18],现场测试中道床加速度传感器的埋设位置位于枕底以下15 cm位置处,测试过程中列车的运行速度为40 km/h. 在离散元模型中通过在钢轨上施加经过平面应变变换的移动荷载模拟列车作用,并对枕底以下15 cm处的道砟颗粒的加速度进行记录. 道床加速度实测结果、仿真结果以及频谱变换如图4所示.
由图4可以看出:现场实测的道床加速度与离散元模型仿真得到的时程曲线的幅值及波形相近. 因现场线路存在不平顺,导致实测结果存在明显的毛刺现象,但现场实测的道床振动加速度的频谱曲线与离散元计算的道床振动加速度的频谱曲线波形及幅值相近,振动的主频成分也相吻合. 其他不同速度的列车荷载下道床加速度测试和计算结果相对比后均能得到同样的验证结果. 从而验证了本文所建立模型的正确性.
3. 计算结果分析
在钢轨上施加移动的轴重荷载进行动力仿真. 考虑到道床的振动和受力是导致道砟粉化破碎的最主要因素[19],本文将对列车荷载作用下道砟颗粒的振动加速度、接触力进行监测. 另外,还将分析列车荷载作用下道砟颗粒之间的相对滑动情况,以期揭示板结道床弹性降低、刚度增大的细观机理.
3.1 道砟颗粒振动
在仿真过程中对道砟颗粒的振动加速度进行监测. 限于篇幅,仅给出工况1中当列车的第二个轮载位于 ⑤ 号轨枕上方时道砟颗粒振动加速度的整体矢量图以及中间 ④ 号、⑤ 号和 ⑥ 号轨枕的局部矢量图,如图5所示. 图中:P为列车荷载;v为列车速度.
由图5可看出:在列车荷载作用下,道砟颗粒的振动加速度矢量分布较为混乱,同一时刻不同位置处道砟颗粒的瞬时加速度存在显著差异,这是道床散体特征的表现. 道砟颗粒的振动加速度受道砟颗粒自身重量、形状以及与周围道砟颗粒的接触状态等诸多因素影响. 分析发现其他计算工况也存在类似的规律,这也反映出即便是板结道床,也具有明显的散体特性. 板结会使得道床的力学特性更趋向于连续介质的表现,但在列车荷载作用下,道砟颗粒之间的板结介质并未将道砟颗粒黏结成一个牢固的整体,板结道床的问题仍然可以划定为散体力学范畴.
进一步量化加载过程中道砟颗粒的振动情况,选定各轨枕正下方15 cm附近的道砟颗粒进行分析. 考虑到粒径较大的道砟颗粒的配位数更多,其振动状态更为稳定,故所选取道砟颗粒的粒径大于4.5 cm. 对选定的道砟颗粒振动加速度进行监测. 限于图幅,给出工况1、2中 ⑤ 号轨枕下方道砟颗粒的振动加速度时程结果,如图6所示.
由图6可以看出:工况1和工况2中的 ⑤ 号轨枕下方道砟颗粒的振动波形较为相似,均是在转向架荷载经过时道砟颗粒的振动加速度增大并出现最大值;工况1中 ⑤ 号轨枕下方道砟颗粒的振动加速度幅值为1.89g,工况2中 ⑤ 号轨枕下方道砟颗粒的振动加速度幅值为2.08g.
各工况中不同轨枕下方所选定道砟颗粒的振动加速度幅值,结果如图7所示.
由图7可知:同一工况中不同道砟颗粒的振动加速度幅值有所不同,这是由于有砟道床的散体特性导致,道砟颗粒受其自身重量、外形及与周围道砟颗粒的接触状态等因素的影响;工况2~4中道砟颗粒振动加速度的幅值均大于工况1,且有依次增大的趋势;对4种计算工况中各自道砟颗粒振动加速度幅值取平均,工况1~4依次为1.94g、2.30g、2.54g和2.69g,工况2~4比工况1分别增大了18.56%、30.93%和38.66%. 因此可知,道床板结会增大道床中道砟颗粒的振动水平,且影响程度随着道床板结程度的增大而增大.
为使结果更具代表性,选取不同工况中 ⑤ 号轨枕下方15 cm~20 cm区域的多个道砟颗粒,统计其各自的加速度,结果表明不同道砟颗粒之间的振动存在差异,但对不同道砟颗粒的振动加速度峰值取平均后,经对比发现道砟颗粒振动仍然存在随着板结程度增大而增大的趋势.
3.2 道砟接触力
为量化道床板结对道砟颗粒接触力的影响,对列车荷载作用下各轨枕下方道砟颗粒之间的接触力进行监测,在仿真过程中对每根轨枕中心位置两侧各0.3 m范围内的道砟颗粒接触力进行遍历并取平均,得到每根轨枕所在区域处的道砟颗粒平均接触力. 限于篇幅,仅给出4种计算工况的 ⑤ 号轨枕处道砟颗粒之间的平均接触力时程曲线,如图8所示.
由图8可看出:各工况中道砟颗粒平均接触力的时程曲线较为相似,均是在列车轮载经过时逐渐增大;4种计算工况中道砟颗粒平均接触力的峰值有所不同,工况1~4的道砟颗粒平均接触力的峰值依次为1.21、1.23、1.28和1.38 kN,工况2~4比工况1分别增大了1.65%、5.79%和14.05%. 由此表明,板结因弹性减低而导致列车荷载对道床的冲击作用加强,从而增大道砟颗粒的接触受力.
4种计算工况的各道砟箱中道砟颗粒平均接触力的幅值如图9所示.
由图9可知:同一工况的不同轨枕下方道砟颗粒平均接触力的幅值有所不同,且道床板结会放大此差异;工况1中的 ③ 号轨枕下方道砟颗粒的受力较 ④ 号轨枕下方的道砟颗粒的受力大7.34%,而工况4中差异增大到了12.47%. 由于有砟道床是散粒体材料,不同轨枕位置处道床的密实程度不尽相同,导致道砟颗粒的受力略有差异. 由此表明,板结会增大道床中不同区域处道床服役状态的差异性,增加了道床刚度的不均匀性.
3.3 道砟颗粒相对滑动
板结道床在宏观方面的表现是使道床弹性丧失,引起轨道刚度急剧增加. 为从细观角度揭示板结道床引起道床刚度增大的机理,本文记录仿真计算过程中道砟颗粒的相对滑动情况. 定义道砟颗粒之间发生滑动的接触数量与总接触数量的比值为滑动分数. 图10 为仿真过程中工况1和工况4的 ④ 号轨枕下方道砟颗粒的滑动分数时程曲线.
由图10可知:工况1和工况4中道砟颗粒滑动分数的时程波形较为相似,均是在列车经过时逐渐增长;工况1中道砟颗粒滑动分数的幅值为0.263,工况4中道砟颗粒滑动分数的幅值为0.122,工况4的比工况1的减少了53.61%. 因而可知,板结会降低列车荷载作用下道砟颗粒之间接触发生相对滑动的几率,使更多的接触由发生滑动的活跃状态转化为不发生滑动的稳定状态. 另外还可看出,虽然工况4板结道床中道砟颗粒的滑动分数相比工况1的洁净道床有所减低,但并不为0,这也表明在列车荷载作用下板结道床内部道砟颗粒之间仍然存在滑动现象,板结道床仍具有散粒体的特征.
图11为列车荷载作用于 ⑤ 号轨枕上方时道砟颗粒之间所有接触在每一时步中的切向滑动位移量. 由于本文所建模型中含有大量粒径较小的脏污颗粒,故设定了较小的计算时步,从而导致道砟颗粒在单一时步中的滑动位移量值较小.
由图11可知:工况1中道砟颗粒在单一时步中瞬时滑动位移的最大值为21.12 nm,工况4的瞬时滑动位置的最大值为7.23 nm,比工况一的降低了65.77%. 可以看出,板结降低了道砟颗粒之间滑动位移的量值.
板结减小道砟之间滑动程度的原因是由于板结介质改变了道砟颗粒之间的接触状态,如图12所示. 为便于表述,将图中所示的道砟颗粒分别编号为1号道砟和2号道砟. 由图12可知:脏污板结填充了道砟颗粒之间的空隙,由于有砟道床为散粒体,外部荷载作用下相邻道砟颗粒之间速度方向不尽相同,可能存在相向运动的速度分量,例如图12中的1号和2号道砟颗粒. 在洁净的道床中,1号和2号道砟颗粒不存在直接相互作用,而是通过彼此相邻的道砟传力,容易引起与相邻道砟之间产生滑动位移;而脏污道床中1号和2号道砟颗粒之间的脏污板结形成了新的传力路径,起到了分担荷载的作用,从而减低了这两颗道砟与周围道砟之间的接触力,导致道砟之间不再滑动或降低滑动位移.
有砟道床在列车荷载作用下道床中道砟颗粒之间会出现相对滑动和错位,消耗了列车荷载传递至道床的能量,并在宏观上表现出道床整体的弹塑性变形,保证了有砟轨道的弹性充足. 而板结道床中的脏污板结材质降低了道砟之间发生相对滑动的几率,同时还会降低滑动程度(相对滑动位移). 这反映出板结抑制了道床中细观角度下相互咬合的道砟颗粒之间的相对运动,导致有砟道床的力学特性更趋向于连续介质材料,使得道床在宏观上表现得更不易被压缩,从而削弱了有砟道床的弹性机能,导致道床刚度增大. 因此对于道床板结严重的地段应及时进行清筛养护作业,以避免因道床刚度过大而引发其他次生病害.
4. 结 论
本文建立了板结道床的离散元分析模型,分析了列车荷载作用下道砟颗粒的振动、受力及颗粒滑动等宏细观力学特性,结果显示:
1) 列车荷载作用下,板结道床中不同道砟颗粒的振动状态存在显著差异,且道砟颗粒之间存在相对滑动现象,说明板结道床仍具有散体特性. 道床板结会提高道砟颗粒的振动水平,并且道床板结越严重对道砟颗粒振动水平的影响越大.
2) 道床板结会增加增大列车荷载对道床的冲击作用,从而增大道砟颗粒的接触受力,增加了道砟颗粒发生破碎劣化的风险;由于脏污颗粒在道床中的位置分布较为随机,导致不同区域处道砟颗粒受力状态的差异性变大,增加了道床刚度的不均匀性.
3) 脏污板结材质会减小道砟颗粒之间发生相对滑动的几率,同时会减小道砟颗粒之间的滑动量值,脏污板结在细观角度下对道砟颗粒相对移动的抑制是引起宏观上的道床弹性丧失、刚度变大的重要原因.
-
黄润秋, 唐川, 李勇, 等. 汶川地震地质灾害研究[M]. 北京: 科学出版社, 2009: 1-30 刘静,张智慧,文力,等. 汶川8级大地震同震破裂的特殊性及构造意义——多条平行断裂同时活动的反序型逆冲地震事件[J]. 地质学报,2008,82(12): 1707-1722 doi: 10.3321/j.issn:0001-5717.2008.12.007LIU Jing, ZHANG Zhihui, WEN Li, et al. The Ms8.0 Wenchuan earthquake co-seismic rupture and tis tectonic implications—an out-of-sequence thrusting event with slip partitioned on multiple faults[J]. Acta Geologica Sinica, 2008, 82(12): 1707-1722 doi: 10.3321/j.issn:0001-5717.2008.12.007 XU C, XU X W, YAO X, et al. Three (nearly) complete inventories of landslides triggered by the May 12,2008 Wenchuan Mw 7.9 earthquake of China and their spatial distribution statistical analysis[J]. Landslides, 2014, 11(3): 441-461 doi: 10.1007/s10346-013-0404-6 XU Q, ZHANG S, LI W L. Spatial distribution of large-scale landslides induced by the 5.12 Wenchuan earthquake[J]. Journal of Mountain Science, 2011, 8(2): 246-260 doi: 10.1007/s11629-011-2105-8 许强, 裴向军, 黄润秋. 汶川地震大型滑坡研究[M]. 北京: 科学出版社, 2009: 53-91 黄润秋,裴向军,崔圣华. 大光包滑坡滑带岩体碎裂及其形成机制研究[J]. 岩石力学与工程学报,2016,35(1): 1-15HUANG Runqiu, PEI Xiangjun, CUI Shenghua. Cataclastic characteristic and formation mechanism of rock mass sliding zone of Daguangbao landslide[J]. Chinese Journal of Rock Mechanism and Engineering, 2016, 35(1): 1-15 黄润秋,裴向军,张伟锋,等. 再论大光包滑坡特征与形成机制[J]. 工程地质学报,2009,17(6): 725-736 doi: 10.3969/j.issn.1004-9665.2009.06.001HUANG Runqiu, PEI Xiangjun, ZHANG Weifeng, et al. Further examination on characteristics and formation mechanism of Daguangbao landslide[J]. Journal of Engineering Geology, 2009, 17(6): 725-736 doi: 10.3969/j.issn.1004-9665.2009.06.001 黄河清,赵其华. 汶川地震诱发文家沟巨型滑坡-碎屑流基本特征及成因机制初步分析[J]. 工程地质学报,2010,18(2): 168-177 doi: 10.3969/j.issn.1004-9665.2010.02.003HUANG Heqin, ZHAO Qihua. Basic characteristics and preliminary mechanism analysis of large scale rock slide-sturzstrom at Wenjiagou triggered by Wenchuan earthquake[J]. Journal of Engineering Geology, 2010, 18(2): 168-177 doi: 10.3969/j.issn.1004-9665.2010.02.003 DAI F C, XU C, YAO X, et al. Spatial distribution of landslides triggered by the 2008 Ms 8.0 Wenchuan earthquake,China[J]. Journal of Asian Earth Sciences, 2011, 40(4): 883-895 doi: 10.1016/j.jseaes.2010.04.010 黄润秋,李为乐. " 5.12”汶川大地震触发地质灾害的发育分布规律研究[J]. 岩石力学与工程学报,2008,27(12): 2585-2592 doi: 10.3321/j.issn:1000-6915.2008.12.028HUANG Runqiu, LI Weile. Research on development and distribution rules of geohazards induced by Wenchuan earthquake on 12th[J]. Chinese Journal of Rock Mechanism and Engineering, 2008, 27(12): 2585-2592 doi: 10.3321/j.issn:1000-6915.2008.12.028 WANG G, HUANG R, LOURENÇO S D N, et al. A large landslide triggered by the 2008 Wenchuan (M8.0) earthquake in Donghekou area:Phenomena and mechanisms[J]. Engineering Geology, 2014, 18(2): 148-157 殷跃平. 汶川八级地震滑坡高速远程特征分析[J]. 工程地质学报,2009,17(2): 153-166 doi: 10.3969/j.issn.1004-9665.2009.02.002YI Yueping. Features of landslides triggered by the Wenchuan earthquake[J]. Journal of Engineering Geology, 2009, 17(2): 153-166 doi: 10.3969/j.issn.1004-9665.2009.02.002 许强,李为乐. 汶川地震诱发大型滑坡分布规律研究[J]. 工程地质学报,2010,18(6): 818-826 doi: 10.3969/j.issn.1004-9665.2010.06.002XU Qiang, LI Weile. Distribution of large-scale landslide induced by the Wenchuan earthquake[J]. Journal of Engineering Geology, 2010, 18(6): 818-826 doi: 10.3969/j.issn.1004-9665.2010.06.002 殷跃平. 汶川八级地震地质灾害研究[J]. 工程地质学报,2008,16(4): 433-444 doi: 10.3969/j.issn.1004-9665.2008.04.001YIN Yueping. Rapid and long run-out feature of landslide ides triggered by the Wenchuan earthquake[J]. Journal of Engineering Geology, 2008, 16(4): 433-444 doi: 10.3969/j.issn.1004-9665.2008.04.001 黄润秋,李为乐. 汶川大地震触发地质灾害的断层效应分析[J]. 工程地质学报,2015,17(8): 19-28HUANG Runqiu, LI Weile. Fault effect analysis of geohazard triggered by Wenchuan earthquake[J]. Journal of Engineering Geology, 2015, 17(8): 19-28 HUANG R Q, LI W L. Analysis of the geo-hazards triggered by the 12 May 2008 Wenchuan earthquake,China[J]. Bulletin of Engineering Geology and the Environment, 2009, 68(3): 363-371 doi: 10.1007/s10064-009-0207-0 黄润秋,张伟锋,裴向军. 大光包滑坡工程地质研究[J]. 工程地质学报,2014,22(4): 557-585HUANG Runqiu, ZHANG Weifeng, PEI Xiangjun. Engineering geological study on Daguangbao landslide[J]. Journal of Engineering Geology, 2014, 22(4): 557-585 HUANG R Q, FAN X. The landslide story[J]. Nature Geoscience, 2013, 5(6): 325-326 HUANG R, PEI X, FAN X, et al. The characteristics and failure mechanism of the largest landslide triggered by the Wenchuan earthquake,May 12,2008,China[J]. Landslides, 2012, 1(9): 131-142 黄润秋,裴向军,李天斌. 汶川地震触发大光包巨型滑坡基本特征及形成机理分析[J]. 工程地质学报,2008,16(6): 730-741 doi: 10.3969/j.issn.1004-9665.2008.06.002HUANG Runqiu, PEI Xiangjun, LI Tianbin. Basic characteristic 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 doi: 10.3969/j.issn.1004-9665.2008.06.002 殷跃平,王猛,李滨,等. 汶川地震大光包滑坡动力响应特征研究[J]. 岩石力学与工程学报,2012,31(10): 1969-1982 doi: 10.3969/j.issn.1000-6915.2012.10.003YIN Yueping, WANG Meng, LI Bin, et al. Dynamic response characteristics of daguangbao landslide triggered by wenchuan earthquake[J]. Chinses Journal of Rock Mechanics and Engineering Geology, 2012, 31(10): 1969-1982 doi: 10.3969/j.issn.1000-6915.2012.10.003 ZHANG Y, CHEN G, ZHENG L, et al. Effects of near-fault seismic loadings on run-out of large-scale landslide:a case study[J]. Engineering Geology, 2013, 166(8): 216-236 裴向军, 黄润秋, 崔圣华, 等. 大光包滑坡岩体碎裂特征及其工程地质意义[J]. 岩石力学与工程学报, 2015(增刊1): 3106-3115PEI Xiangjun, HUANG Runqiu, CUI Shenghua, et al. The rock mass cataclastic characteristic of Daguangbao landslide and its engineering geological significance[J]. Chinese Journal of Rock Mechanism an Engineering, 2015(S1): 3106-3115 YANG C M, CHENG H Y, TSAO C C, et al. The kinematics and initiation mechanisms of the earthquake-triggered Daguangbao landslide[J]. Geophysical Research Abstracts, 2015, 17: 13014-1-13014-11 doi: 10.1007/978-4-431-56205-4_19 张伟锋,黄润秋,裴向军. 大光包滑坡运动特征及其过程分析[J]. 工程地质学报,2015,23(5): 866-885ZHANG Weifeng, HUANG Runqiu, PEI Xiangjun. Analysis on kinematics and movement process of Daguangbao landslide[J]. Journal of Engineering Geology, 2015, 23(5): 866-885 CHIGIRA M, WU X Y, INOKUCHI T, et al. Landslides induced by the 2008 Wenchuan earthquake,Sichuan,China[J]. Geomorphology, 2010, 118(3/4): 225-238 崔圣华. 强震巨型滑坡滑带碎裂岩体微细观分析及静动力破损机制研究 [D]. 成都: 成都理工大学, 2014 WANG W N, CHIGIRA M, FURUYA T. Geological and geomorphological precursors of the Chiu-fen-erh-shan landslide triggered by the Chi-chi earthquake in central Taiwan[J]. Engineering Geology, 2003, 69(2): 1-13 CHANG K J, TABOADA A, CHAN Y C. Geological and morphological study of the Jiufengershan landslide triggered by the Chi-Chi Taiwan earthquake[J]. Geomorphology, 2005, 71(3): 293-309 裴向军,郝永峰,张军新,等. 老鹰岩滑坡成因机制与运动特征研究[J]. 地质灾害与环境保护,2010,21(4): 28-32 doi: 10.3969/j.issn.1006-4362.2010.04.005PEI Xiangjun, HAO Yongfeng, ZHANG Junxin, et al. Movemwnt characteristics and mechanism of Laoyingyan landslide triggered by Wenchuan earthquake[J]. Journal of Geological Hazards and Environment Preservation, 2010, 21(4): 28-32 doi: 10.3969/j.issn.1006-4362.2010.04.005 张培震. 青藏高原东缘川西地区的现今构造变形、应变分配与深部动力过程[J]. 中国科学,2008,38(9): 1041-1056 doi: 10.3321/j.issn:1006-9267.2008.09.001ZHANG Peizhen. The tectonic deformation,strain distribution and deep dynamic process eastern in the margin of the Tibetan Plateau[J]. Science China, 2008, 38(9): 1041-1056 doi: 10.3321/j.issn:1006-9267.2008.09.001 四川省地质矿产勘查开发局化探队. 中华人民共和国区域地质图说明书[R]. 德阳: [s.n.], 1995 HUMAIR F, PEDRAZZINI A, EPARD J L, et al. Structural characterization of Turtle Mountain anticline (Alberta,Canada) and impact on rock slope failure[J]. Tectonophysics, 2013, 605(3): 133-148 CUI S H, PEI X J, HUANG R Q. Effects of geological and tectonic characteristics on the earthquake-triggered Daguangbao landslide, China[J]. Landslides, 2018, 15(4): 649-667 doi: 10.1007/s10346-017-0899-3 李勇,黄润秋,周荣军,等. 龙门山地震带的地质背景与汶川地震的地表破裂[J]. 工程地质学报,2009(1): 3-18 doi: 10.3969/j.issn.1004-9665.2009.01.002LI Yong, HUANG Runqiu, ZHOU Rongjun, et al. Geological background of longmen shan seismic belt and surface ruptures in Wenchuan earthquake[J]. Journal of Engineering Geology, 2009(1): 3-18 doi: 10.3969/j.issn.1004-9665.2009.01.002 ZHU S, SHI Y, LU M, et al. Dynamic mechanisms of earthquake-triggered landslides[J]. Science China, 2013, 56(10): 1769-1779 doi: 10.1007/s11430-013-4582-9 OUIMET W B. Landslides associated with the May 12,2008 Wenchuan earthquake:Implications for the erosion and tectonic evolution of the Longmen Shan[J]. Tectonophysics, 2010, 491(1): 244-252 罗璟,裴向军,黄润秋,等. 强震作用下滑坡岩体震裂损伤程度影响因素研究[J]. 岩土工程学报,2015,37(6): 1105-1114LUO Jing, PEI Xiangjun, HUANG Runqiu, et al. Influencing factors for damage degree of shattered landslide rock mass under high seismic action[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(6): 1105-1114 黄润秋,王运生,裴向军,等. 4.20芦山Ms7.0级地震地质灾害特征[J]. 西南交通大学学报,2013,48(4): 581-589 doi: 10.3969/j.issn.0258-2724.2013.04.001HUANG Runqiu, WANG Yunsheng, PEI Xiangjun, et al. Characteristics of co-seismic landslides triggered by the Lushan Ms7.0 Earthquake on the 20th of April,Sichuan Province,China[J]. Journal of Southwest Jiaotong University, 2013, 48(4): 581-589 doi: 10.3969/j.issn.0258-2724.2013.04.001 DAVIS G H, REYNOLDS S. Structural geology of rocks and regions[M]. New York: Wiley, 1984: 438-447 刘小丽,周德培. 有软弱夹层岩体边坡的稳定性评价[J]. 西南交通大学学报,2002,37(4): 382-386 doi: 10.3969/j.issn.0258-2724.2002.04.007LIU Xiaoli, ZHOU Depei. Stability evaluation of rock mass slope with weak intercalated layers[J]. Journal of Southwest Jiaotong University, 2002, 37(4): 382-386 doi: 10.3969/j.issn.0258-2724.2002.04.007 CHIGIRA M, WANG W, FURUYA T, et al. Geological causes and geomorphological precursors of the Tsaoling landslide triggered by the 1999 Chi-Chi earthquake,Taiwan[J]. Engineering Geology, 2003, 68(S3/4): 259-273 -