道路开挖位置对边坡稳定性影响的数值模拟
陈晓利1), 王明明2), 张凌3),*
1)中国地震局地质研究所, 活动构造与火山实验室, 北京 100029
2)四川省地震局, 成都 610041
3)中国地震应急搜救中心, 北京 100049
*通讯作者: 张凌, 1961年生, 高级工程师, 研究方向为灾害应急与救援, 电话: 13901043361, E-mail: zhangling903@163.com

〔作者简介〕陈晓利, 女, 1969年生, 研究员, 主要研究方向为地震地质灾害与工程地质, 电话: 010-62009056, E-mail: chenxl@ies.ac.cn

摘要

公路边的滑坡崩塌在中国西南地区是常见的一类地质灾害。道路的开挖一方面改变了边坡原有的地表形态, 形成了新的较陡的边坡, 如果缺乏保护措施, 新的陡坡很容易失稳, 从而形成以滚石或浅层滑坡为主的边坡破坏。另一方面, 道路开挖有可能会改变原有坡体的结构, 降低坡体的安全系数, 引起较大规模的滑坡, 造成更为严重的灾害。因此, 坡体开挖的位置对其稳定程度会产生一定的影响。文中以岩质边坡为例, 在建立边坡模型的基础上, 采用二维极限平衡数值模拟方法, 分析了道路开挖位置对整体边坡稳定性的影响; 同时, 对加载地震作用的边坡稳定性也进行了模拟分析。研究结果表明, 对于有潜在滑动面的坡体, 在不同位置进行开挖对坡体稳定性的影响是不同的: 在坡脚处或接近滑动面滑出位置开挖道路, 将减小阻滑力而导致整体稳定性降低; 在接近坡顶处开挖道路则会因卸载坡体物质而提高整体稳定性。合理的开挖位置与坡体坡度有密切关系, 在不降低整体安全系数的情况下, 坡度较陡的坡体, 其开挖位置相应要高一些。无论坡体角度大小, 地震作用能够显著降低边坡的稳定性。

关键词: 边坡稳定性; 道路开挖; 极限平衡法; 安全系数; 地震加载
中图分类号:P315.2 文献标志码:A 文章编号:0253-4967(2018)06-1390-12
SIMULATION STUDY OF ROAD-CUT EFFECTS ON SLOPE STABILITY
CHEN Xiao-li1), WANG Ming-ming2), ZHANG Ling3)
1)Key Laboratory of Active Tectonics and Volcano, Institute of Geology, China Earthquake Administration, Beijing 100029, China
2)Sichuan Earthquake Agency, Chengdu 610041, China
3)National Earthquake Response Support Service, Beijing 100049, China
Abstract

Landslides and rock falls along the highway are common geological hazards in Southwest China. As an influencing factor on potential landslides behavior, roads or distance to roads have been successfully used in landslide susceptibility assessments in mountainous area. However, the relationship between the road-cut and the slope stability is not clear. Therefore, we performed two-dimensional slope stability calculation using the general limit equilibrium(GLE)method incorporated in the software SLOPE/W of GeoStudio for stability analysis of slopes. Our studies show that the man-made roads influence on the slope stability mainly exists in two ways: One is to create a new steep slope, which will result in rock falls and shallow landslides along the roads; the other is to influence the stability of the original slope, which will result in comparatively huge landslides. For the latter, our simulation study reveals that the road location, namely at which part of a natural slope to construct a road is important for the slope stability. For a natural slope with a potential slip surface, if a road is constructed at or near the slope toe where the potential slip surface surpasses, it will greatly degrade the slope’s factor of safety(Fs) and make the slope unstable; however, if a rode-cut is near the top of the slope, it will increase the slope’s Fs and make the slope more stable. The safety location is different for different slope angle, steeper slope needs a higher location for a safety road-cut in comparison with gentle slopes. Moreover, the slope stability decreases when loading a seismic force and it varies with the slope angle. Firstly, the Fs decreases when the slope angle increasing, and when the slope angle reaches 45°, the Fs then becomes greater with the slope angle increasing.

Keyword: slope stability; road-cut; limit equilibrium method; factor of safety(Fs); seismic loading
0 引言

滑坡崩塌灾害在山岳地区分布较为广泛, 在阻塞交通、 河流、 毁坏建筑设施等的同时, 也严重威胁到当地居民的生命财产安全。山岳地区迅速发展的道路交通建设在促进当地的运输、 贸易、 旅游业等经济产业发展的同时, 对周围地质地貌环境的侵蚀、 沉积过程等也会产生重要的影响, 而沿公路分布的滑坡、 崩塌等地质灾害则是这种影响作用的表现之一, 如川藏公路时有发生滑坡崩塌阻塞道路的现象。在山岳地区修筑公路需要翻山越岭, 道路往往需要盘山而上, 尽管在道路的建设过程中也会采取一些边坡加固措施, 但是小规模的滑坡、 崩塌还是难以避免, 降雨或地震作用等则加剧了这类灾害的发生。滑坡崩塌沿公路密集分布已成为山岳地区强震诱发滑坡的特点之一(黄润秋等, 2008; 许冲等, 2013)。例如, 2013年MS7.0芦山地震在四川省宝兴县附近S210省道5km长的路段上触发了79个滑坡、 崩塌(Chen et al., 2014)。而强震中由于滑坡崩塌造成的交通阻塞问题将严重影响震后救援工作的顺利进行, 因而震后的道路交通安全保障受到了更多的关注。

滑坡的影响因素多种多样, 与边坡所处的地质背景、 坡体结构类型、 物质组合、 地形地貌条件、 水文地质条件(内部条件)以及发震时的降雨、 融雪(外部条件)等因素密切相关(Guzzetti et al., 1999; 丁彦慧等, 2000; Bommer et al., 2002; 祁生文等, 2004; Wang et al., 2007; 黄润秋等, 2008; 殷跃平, 2008; Su et al., 2017)。随着人类对自然的改造, 更多的人为因素参与到滑坡演化的进程中, 道路交通的建设就是其中之一。作为滑坡的外部影响因素之一, 道路的开挖一方面改变了边坡原有地表的形态, 形成了新的较陡的边坡, 如果缺乏必要的防护措施, 则新形成的边坡很容易发生失稳, 形成规模不等的滚石或浅层滑坡; 另一方面, 一些学者认为道路的修建影响到整个原有坡体的排水体系, 并且路面交通带来的震动也会影响边坡的应力状况, 进而影响到边坡的稳定状态(Devkota et al., 2013)。因此, 目前在很多的滑坡敏感性、 稳定性评价中也把与道路的距离(Distance to roads, RDS)作为评价参数之一。

然而, 具体到某一特定边坡, 虽然开挖施工造成的新的较陡凌空坡面会成为滑坡崩塌隐患, 但道路开挖对原始边坡本身稳定性的影响还没有得到充分的认识。与较小规模的滑坡崩塌相比, 具有深大埋深滑动面的滑坡具有更强的致灾性, 造成的危害也更为严重(许强等, 2010; Chen et al., 2012; 黄润秋等, 2014)。2014年MS6.5鲁甸地震在震区触发了大量的滑坡崩塌灾害, 在震后进行的野外勘查工作中发现, 大量的滑坡崩塌沿公路发生, 但是灾害体的规模大小及危害程度都不尽相同。较大规模的滑坡崩塌一般都发生在山脚附近, 相比之下, 半山腰的滑坡崩塌规模较小。图1为鲁甸地震在震中附近照壁山触发的滑坡, 其中, 图1a中的滑坡发生在半山腰, 规模较小, 没有对公路造成严重损坏, 清理路面的堆积物后即可通行; 而图1b中的滑坡则发生在山脚, 规模巨大, 原有的公路被滑坡体切断, 需要在堆积体上开挖出1条临时通道来保证交通的通畅, 修复工程极为浩大。

图1 2014年鲁甸地震中的公路滑坡
a 公路位于半山腰; b 公路位于山脚
Fig. 1 Landslides along the road during the 2014 Ludian earthquake.

总体而言, 修筑公路在开挖或填筑路基的过程中, 无论采取什么样的施工步骤和方式, 都在不同程度上和范围内改变了斜坡岩上体的原有形状和自然平衡(杨有辉等, 2006)。为了深入研究公路开挖对原有坡体稳定性的影响, 本文通过建立1个岩质边坡二维模型, 分析在边坡的不同位置进行道路开挖对整体边坡稳定性的影响; 同时, 对地震作用对于边坡稳定性的影响也进行了分析, 并探讨了不同坡体角度对边坡稳定性的影响。

1 方法
1.1 极限平衡法简介

一般地, 人们把斜坡岩体顺坡向下的一切运动现象称为滑坡。无论是自然边坡还是人工边坡, 其稳定性分析都是研究的核心问题。如何确定具有最小安全系数(Factor of safety, Fs)的滑动面是边坡稳定研究的关键问题。安全系数是通过比较阻止坡体滑动的力量与促使坡体滑动的力量的大小来获得的, 其计算涉及坡体的物质组成及孔隙水压力等参数, 可简单地表达为: Fs=SresistanceSMobilized, 分子表示阻止坡体滑动的力的总和, 分母表示驱使坡体滑动的力的总和。如果阻止滑动的力小于驱动滑动的力, 即Fs< 1, 则坡体处于失稳状态; Fs> 1 , 则坡体处于稳定状态; Fs=1则是临界状态。

在过去的几十a中, 涌现出很多种计算安全系数的方法, 除有限元方法使用有限元计算应力外, 其他的方法都是基于极限平衡原理。这些以条分方法为基础的极限平衡方法都很相似, 不同之处在于满足的平衡方程以及对条间力和条间力关系的假设不同(Fellenius, 1936; Bishop, 1955; Morgenstern et al., 1965)。通用极限平衡法(General limit equilibrium method)由Fredlund 教授于20世纪70年代提出, 其建立在2个安全系数方程的基础上, 一个是关于滑坡体条块整体的力平衡方程, 另一个是关于力矩的平衡方程(Fredlund et al., 1977)。其允许条间的剪切力-法向力在一定范围内变化, 同时也不受滑动面形状的限制, 适用于运动学上允许的任意滑面形状, 因而在岩土边坡的稳定性评价中获得了广泛的应用。

近年来, 计算机技术的迅速发展使得应用极限平衡方法进行边坡分析的可操作性更强, 其中, GeoStudio中Slope/W模块是采用该方法进行边坡稳定性评价的有力工具之一, 本文的数值模拟计算正是基于此软件, 具体计算时选用的是数学意义更为严格的Morgenstern-Price(M-P)方法。

1.2 边坡模型的建立

本文利用Slope/W岩土结构稳定性分析软件建立了1个二维模型, 用于分析道路开挖对整个坡体的稳定性影响。模型及坡高、 长度、 水平距离、 坡体坡角等参数定义如图2所示。其中, 坡体高度H=80m, 开挖道路的宽度w=10m。

图2 边坡模型及参数定义示意图
α 坡角; L 距坡顶水平距离; h 开挖高度; H 坡体高度
Fig. 2 Slope model and its parameters definition.

1.3 模型的输入参数

考虑到本文的研究旨在讨论道路开挖对边坡安全系数的影响, 因此, 为便于分析, 在此设定坡体组成为成分单一的岩质边坡。计算中使用到的边坡岩石物理力学参数如表1所示。

表1 边坡岩石物理力学参数 Table1 Rock parameters in the study area

结合实际道路开挖所选择的地形特点, 在具体分析中, 以坡角α =30° 和45° 为例, 保持坡体高度(H=80m)不变, 计算在坡体不同高度位置(h)进行道路开挖后对整体稳定性的影响。为了研究地震作用的影响, 加载地震作用时取水平地震作用系数k=0.3。

2 计算结果及分析

首先对原始坡体的安全系数进行计算, 之后, 从坡脚处开始筑路, 计算坡脚开挖之后坡体的安全系数; 然后, 每次路面抬升h=10m, 依次计算坡体在不同高度处开挖后的安全系数。同时, 加载地震作用, 对每一不同高度分别计算地震作用下的安全系数。

2.1 坡角为30° 的坡体

对于坡角为30° 的坡体, 计算所得坡体未施工前的安全系数为4.567(图3a); 在加载水平地震作用系数0.3时, 坡体的安全系数降为2.876(图3b)。坡脚筑路开挖后坡体的安全系数为4.490(图4a), 加载地震作用后安全系数为2.850(图4d)。在坡体不同位置开挖道路后的安全系数变化如图5所示。

图3 坡角为30° 的坡体最小安全系数及潜在滑动面
a 初始状态时; b 加载地震作用系数k=0.3时
Fig. 3 Minimum Fs and its corresponding slip plane of a 30° slope.

图4 坡角为30° 的坡体在坡脚处开挖时的最小安全系数及潜在滑动面
a 初始状态时; b 加载地震作用系数k=0.3时
Fig. 4 Minimum Fs and its corresponding slip plane of a 30° slope with road-cut at the slope bottom.

图5 坡角为30° 的坡体的最小安全系数随开挖道路高度(h)的变化
a 没有地震作用; b 地震作用系数k=0.3
Fig. 5 Minimum Fs of a 30° slope with road-cut at the different slope height.

2.2 坡角为45° 的坡体

对于坡角为45° 的坡体, 计算所得坡体的安全系数为4.051(图6a); 在加载地震作用系数0.3时, 坡体的安全系数降为2.699(图6b)。在坡脚筑路开挖后坡体的安全系数为3.882(图7a), 加载地震作用后安全系数为2.586(图7b)。在不同位置开挖道路后的安全系数变化如图8所示。

图6 坡角为45° 的坡体最小安全系数及潜在滑动面
a 初始状态时; b 加载地震作用系数k=0.3时
Fig. 6 Minimum Fs and its corresponding slip plane of a 45° slope.

图7 45° 的坡体在坡脚处开挖时的最小安全系数及潜在滑动面
a 初始状态时; b 加载地震作用系数k=0.3时
Fig. 7 Minimum Fs and its corresponding slip plane of a 45° slope with road-cut at the slope bottom.

图8 45° 的坡体的最小安全系数随开挖道路高度(h)的变化
a 没有地震作用; b 地震作用系数k=0.3时
Fig. 8 Minimum Fs of a 45° slope with road-cut at the different slope height.

3 讨论

通过本文的模拟计算结果可以看出, 道路开挖在改变边坡形态的同时, 对坡体的安全系数存在着不同程度的影响。无论是30° 坡体还是45° 坡体, 其安全系数计算结果显示, 与原有坡体的安全系数相比, 在坡脚附近的道路开挖会较大幅度地降低原始坡体的安全系数。而随着路面开挖高度的改变, 安全系数反而会增高, 甚至超过原始坡体的安全系数。对于30° 的坡体而言, 开挖路面高度h≥ 30m时原有坡体的稳定性安全系数大于原始坡体自身安全系数(图5), 表明在该高度以上建设道路不仅没有破坏原始坡体的稳定性, 反而对于原有坡体能够起到稳定作用。而对于具有45° 的坡体而言, 要不降低原有坡体安全系数, 道路开挖的高度应> 50m(图8)。由此可见, 对于具有不同坡度的坡体, 在不同高度处进行的道路开挖会对坡体原始的稳定性产生不同程度的影响。

那么坡体的坡角对安全系数的影响又是如何呢?本文在前述模型的基础上, 通过改变距坡顶水平距离L的大小(60m≤ L≤ 210m, 改变步长10m, 相当于坡角α 的变化范围为22° ≤ α ≤ 63° ), 计算了坡角变化对坡体稳定性、 破坏形式的影响。

计算结果如图9所示。对于滑坡, 普遍存在的认识是坡度越高, 滑坡发生的可能性就越大(Keefer, 1984; 王余庆等, 2001; Wen et al., 2004; 黄润秋等, 2008)。但是, 通过计算发现, 对于没有先存滑动面的坡体, 随着坡度的增大, 安全系数不断降低, 坡度达到45° 左右时安全系数值最小, 之后, 安全系数则是随坡度增大而增大(图9), 坡度在67° 左右时, 安全系数达到最大。

图9 边坡最小安全系数随坡角(水平距离L)改变产生的变化Fig. 9 Slope minimum Fs vs slope angle(horizontal distance L).

值得注意的是, 本研究采用的边坡稳定性分析软件GeoStudio中的Slope/W模块提供的极限平衡分析方法是建立在静力学原理上的, 该方法的平衡公式中缺少应力应变关系, 不能满足位移的兼容性, 因此, 对一些特殊状态缺乏合理的解释, 如坡体角度 ≥ 70° 时, 计算过程出现不收敛的情况。此外, 在进行地震动加载时, 采用拟静力分析法来表示地震动动力载荷的效果, 因而对地震作用下的边坡反应分析有一定的局限性。并且, 具体计算地震动载荷时, 该软件提供的分析方法认为垂直方向的地震作用对安全系数几乎没有影响。因为垂直方向的作用力在改变分析对象的重量时, 对其底部的正应力和剪切力也会产生相应的改变, 这2种作用的效果可以相互抵消, 因此整体上, 垂直方向的地震作用对安全系数几乎没有影响。因此, 应用该软件模块进行的静力和拟静力分析有一定的应用局限。

总之, 滑坡是一个复杂的地质力学过程。具体斜坡的演化过程具有随机性、 复杂性和不确定性的特点。任何一个物理模型只是对现实世界中相应认识状态的一种反映, 本文模型的设置及物理参数的设定是对边坡问题的概化描述, 不能全面反映道路开挖对坡体稳定性的影响。具体到山岳地区的道路施工建设, 应依据坡体的外形、 物质组成以及地震动时程输入等, 建立合理的物理模型进行更加细致的分析研究。

4 结论

道路作为一种滑坡影响因素在边坡稳定性分析中越来越受到重视。山区的道路建设不可避免地会改变沿线的工程地质环境平衡, 或轻或重地引发滑坡崩塌等地质灾害, 选择合理的开挖位置则会减少或减轻这类灾害。本文采用极限平衡法研究了道路开挖对边坡稳定性的影响, 获得的认识有助于山岳地区建设道路时选择合理的开挖位置, 并采取针对性的防治措施。

(1)道路开挖施工会影响坡体的安全系数, 但并不是所有的道路开挖都会降低坡体的安全系数。在坡脚处的道路开挖会最大程度地导致坡体安全系数的降低, 而接近坡顶处的道路施工则会提高坡体安全系数。

(2)对于具有不同坡角的坡体, 在不同高度处进行的道路开挖会对坡体原始的稳定性产生不同程度的影响。本模型中, 对于30° 的坡体而言, 开挖路面高度h≥ 30m时形成的安全系数大于原始坡体自身安全系数, 而对于坡角> 45° 的坡体而言, 不影响原始坡体安全系数的道路开挖高度应≥ 50m。

(3)地震作用会较大程度地降低边坡的安全系数, 加载地震作用后, 无论是否有道路施工, 坡体安全系数都有较为明显的下降。

致谢 感谢审稿人的宝贵建议, 使得本文得以完善!

The authors have declared that no competing interests exist.

参考文献
1 丁彦慧, 王余庆, 孙进忠, . 2000. 地震崩滑预测方法及其工程研究应用[J]. 工程地质学报, 8(4): 475-480.
DING Yan-hui, WANG Yu-qing, SUN Jin-zhong et al. SUN Jin-zhong 2000. Research on the method for prediction of earthquake-induced land slides and its application to engineering projects[J]. Journal of Engineering Geology, 8(4): 475-480(in Chinese). [本文引用:1]
2 黄润秋, 李为乐. 2008. “5·12”汶川大地震触发地质灾害的发育规律研究[J]. 岩石力学与工程学报, 27(12): 2585-2592.
HUANG Run-qiu, LI Wei-le. 2008. A study on the development and distribution rules of geohazards triggered by “5·12” Wenchuan earthquake[J]. Chinese Journal of Rock Mechanics and Engineering, 27(12): 2585-2592(in Chinese). [本文引用:3]
3 黄润秋, 张伟峰, 裴向军. 2014. 大光包滑坡工程地质研究[J]. 工程地质学报, 22(4): 557-585.
HUANG Run-qiu, ZHANG Wei-feng, PEI Xiang-jun. 2014. Engineering geological study on Daguangbao land slide[J]. Journal of Engineering Geology, 22(4), 557-585(in Chinese). [本文引用:1]
4 祁生文, 伍法权, 刘春玲, . 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). [本文引用:1]
5 王余庆, 辛鸿博, 高艳平, . 2001. 预测岩土边坡地震崩滑的综合指标法研究[J]. 岩土工程学报, 23(3): 311-314.
WANG Yu-qing, XIN Hong-bo, GAO Yan-ping, et al. 2001. Study on comprehensive index method for predicating earthquake-induced land slide[J]. Chinese Journal of Geotechnical Engineering, 23(3): 311-314(in Chinese). [本文引用:1]
6 许冲, 徐锡伟, 郑文俊, . 2013. 2013年四川省芦山 “4·20” 7. 0级强烈地震触发滑坡[J]. 地震地质, 35(3): 641-660. doi: DOI: 103969/j. issn. 0253-4967. 2013. 03. 018.
XU Chong, XU Xi-wei, ZHENG Wen-jun, et al. 2013. Land slides 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). [本文引用:1]
7 许强, 李为乐. 2010. 汶川地震诱发大型滑坡分布规律研究[J]. 工程地质学报, 18(6): 818-826.
XU Qiang, LI Wei-le. 2010. Distribution of large scale land slides induced by the Wenchuan earthquake[J]. Journal of Engineering Geology, 18(6): 818-826(in Chinese). [本文引用:1]
8 杨有辉, 郑治. 2006. 公路滑坡预测预报研究现状及前景展望[J]. 公路交通技术, (6): 18-22.
YANG You-hui, ZHENG Zhi. 2006. Current situation and prospect of expressway land slide forecast[J]. Technology of Highway and Transport, (6): 18-22(in Chinese). [本文引用:1]
9 殷跃平. 2008. 汶川八级地震地质灾害研究[J]. 工程地质学报, 16(4): 433-444.
YIN Yue-ping. 2008. Research on the geo-hazards triggered by Wenchuan earthquake, Sichuan[J]. Journal of Engineering Geology, 16(4): 433-444(in Chinese). [本文引用:1]
10 Bishop A W. 1955. The use of the slip circle in the stability analysis of slopes[J]. Geotechnique, 5: 7-17. [本文引用:1]
11 Bommer J J, Carlos E, Rodríguez C R. 2002. Earthquake-induced land slides in Central America[J]. Engineering Geology, 63(3): 189-220. [本文引用:1]
12 Chen X L, Ran H L, Yang W T. 2012. Evaluation of factors controlling large earthquake-induced land slides by the Wenchuan earthquake[J]. Natural Hazards and Earth System Sciences, 12(12): 3645-3657. [本文引用:1]
13 Chen X L, Yu L, Wang M M, et al. 2014. Brief communication: Land slides triggered by the MS=7. 0 Lushan earthquake, China[J]. Natural Hazards and Earth System Sciences, 14(5): 1257-1267. [本文引用:1]
14 Devkota K C, Regmi A D, Pourghasemi H R, et al. 2013. Land slide susceptibility mapping using certainty factor, index of entropy and logistic regression models in GIS and their comparison at Mugling-Narayanghat road section in Nepal Himalaya[J]. Natural Hazards, 65(1): 135-165. [本文引用:1]
15 Fellenius W. 1936. Calculation of the stability of earth dams[C]. Proceedings of the Second Congress of Large Dams, 4: 445-463. [本文引用:1]
16 Fredlund D G, Krahn J. 1977. Comparison of slope stability methods of analysis[J]. Canadian Geotechnical Journal, 14(3): 429-439. [本文引用:1]
17 Guzzetti F, Carrara A, Cardinali M, et al. 1999. Land slide hazard evaluation: A review of current techniques and their application in a multi-scale study, Central Italy[J]. Geomorphology, 31(1-4): 181-216. [本文引用:1]
18 Keefer D K. 1984. Land slides caused by earthquakes[J]. Geological Society of America Bulletin, 95(4): 406-421. [本文引用:1]
19 Morgenstern N R, Price V E. 1965. The analysis of the stability of general slip surfaces[J]. Geotechnique, 15(1): 79-93. [本文引用:1]
20 Su L J, Hu K H, Zhang W F, et al. 2017. Characteristics and triggering mechanism of Xinmo land slide on 24 June 2017 in Sichuan, China[J]. Journal of Mountain Science, 14(9): 1689-1700. [本文引用:1]
21 Wang H B, Sassa K, Xu W Y. 2007. Analysis of a spatial distribution of land slides triggered by the 2004 Chuetsu earthquakes of Niigata Prefecture, Japan[J]. Natural Hazards, 41: 43-60. [本文引用:1]
22 Wen B P, Wang S J, Wang E Z, et al. 2004. Characteristics of rapid giant land slides in China[J]. Land slides, 1(4): 247-261. [本文引用:1]