为什么地震会引发灾难性的山体滑坡？

  每一次地震都能引发地表位移吗？地震每年都会造成人员伤亡、建筑物和基础设施损坏。长期以来，地震引发的山体滑坡一直被视为次生灾害，主要表现为落石、滑坡、泥石流，可造成严重的损失[1]。20世纪，近80次地震引发了10万到100万次地表位移，造成数万人死亡[2]。

1. 每次地震都能引发地表移动吗？

  为了回答这个问题，基弗（Keefer）搜集并分析了1811年至1980年间世界各地地震引发的所有地表移动，从而确定产生地表移动的地震特征与诱发地表移动特征间的相关性：

• 结论1：大多数地表失稳破坏由中到高震级地震触发[3]，通常大于5级（图1）。不过，我们要谨慎对待“最小级”或震级阈值这类概念，因为其并未考虑地震前边坡的稳定状态。边坡状态通常很难准确评估，它是各种“荷载”（地震、暴雨）长期作用的结果。
• 结论2：大多数地表移动不会离震源太远，最远作用距离取决于地震震级。地震震级越高，越有可能引发远离震中的地表移动（图1）[5]。

2. 地震引发的地表移动有何特点？

  地震通常会引发地面表层移动。失稳体积有时可达数百万立方米，1970年5月31日的秘鲁地震就是典型例证：瓦斯卡兰山的崩塌碎屑流时[6]，约5000万立方米的岩石、冰、雪混合物，以210至280 km/h的速度移动，覆盖了超16.5 km的超长距离，导致位于山谷的兰拉希尔卡（Ranrahirca）和云盖（Yungay）村庄约2万人死亡[7]。

  地震引起的滑坡“强度”可由以下特征表示：

• 滑坡体积；
• 滑坡在运动过程中的流动性；
• 滑坡移动停止后的覆盖范围。

  基弗（1984年）的工作表明，地表移动“强度”取决于地震震级：震级越高，可能引发的地表移动“强度”就越大。

  一些地震甚至能使以前发生过位移的古滑坡再次复活，2010年1月12日的海地地震就是如此：这次地震在距震中46公里的范围内引发了约4500次地表移动。至于2010年地震前记录的1273起山体滑坡中，有572处山体滑坡因2010年地震再次发生滑坡（图2）[8]。

3. 地表移动对居民和地区有何影响？

  长期以来，地震引起的地表移动一直被视为与地面震动伴随的次生效应。然而，全球许多案例表明，地表移动的破坏潜力巨大。例如，2001年1月13日萨尔瓦多地震期间，拉斯科林纳斯滑坡（图3）夺走了500多人的生命，接近地震遇难者总数的三分之二。佩特利（Petley）等人（2006）[9]估计，2005年10月8日巴基斯坦地震期间，有30%的受害者（约26500人）死于地震引发的地表移动。2008年5月12日中国的汶川大地震引发了15,000多次地表移动，仅此一项就夺去了近20,000人的生命，占地震遇难者总数的近三分之一[10]。

  地震引起的地表移动可能有两类影响：

• 直接影响：因滑坡体中的碎块埋压建筑物或车辆、损坏/破坏建筑物及交通网络等造成的死亡。
• 间接影响：地表移动产生的碎块毁坏或阻塞道路，妨碍援助受害者或将其疏散至医疗中心。

  有时，地表移动产生的碎屑会阻塞河道，形成堰塞体，进而在上游形成堰塞湖。中国2008年5月12日的地震期间，该现象频发（达数百次）。这段动画[11]展示了2008年地震引起地表移动前后的山坡状况。

  2016年11月14日新西兰地震期间也发生了类似现象，如图4所示。

  根据新西兰地质灾害研究中心监测[12]，2016年11月14日的地震及其余震引发了80,000至100,000次地表移动，形成约150座堰塞坝。由于这些坝寿命长短不一（几分钟至几百年不等），自然灾害预防机构必须明确以下两类主要问题：

• 如何随着时间推移妥善维护这些堰塞坝，避免附近城市和地区遭受洪水灾害？
• 特别是在遭逢暴雨、水位上涨的情况下，堰塞坝上游形成的湖泊将如何演变？
• 迄今为止，鲜有科学研究准确预测这些情况下会发生什么。

4. 我们如何解释这些边坡失稳？

  要确定边坡是否稳定，必须比较：

• 可能促使边坡失稳的驱动力。这些力由地震、暴雨引起。
• 阻碍地表移动的抵抗力。这些力源于材料的固有性质，特别是其黏聚力和内摩擦角，莫尔-库仑定律中，这两个参数能反映土壤抵抗破坏的能力[13]。

  地震可能会伴随一些诱发边坡失稳的现象。

4.1. 场地效应（site effect）

  如图5所示，地震期间，在地球表面某已知点记录到的地表移动取决于三个参数：

• 与震源有关的参数。该参数反映了震源处的破坏机制、断层长度等地震性质。
• 与震波在地下震源到表面移动安置的传感器间的传播相关的参数：源头发出的体积波在传输过程中会随着所穿越环境性质的变化而变化。

• 

  与当地地质和地形条件（即传感器下方和附近的物质类型、地形）影响路径上地震波传播和所记录地表移动的相关参数，此参数称为场地效应。由于场地效应，坡面的地震位移情况随地点而异。与岩质场地（图6中的Campos站）边坡附近测得的地震移动量相比，沉积盆地（图6中的SCT站）或坡顶记录的地震移动通常振幅更大、持续时间更长。现场观测结果[14]和数值模拟结果[15]都表明，场地效应确实会放大地表移动。

4.2. 材料力学性能下降

  反复地震会导致构成边坡的岩土材料变形，可能降低材料的固有特性，诱发边坡整体失稳。

4.3. 孔隙水压力增加

  地震产生的压缩波在松散土体中传播时，往往会引起土体压缩。如果土体处于饱和状态，即水占据了土体颗粒间的所有孔隙（或间隙），那么随着土体压缩，土体孔隙中的水压力增加，水往往会从土体孔隙中排出到地表。但是，若地震快速发生、震级强大或在土壤排水完成前再次发生，水就无法从土体中排出。土体水压越高，土壤的力学强度越低，边坡受破坏的可能性就越大。

5. 如何评估地震期间边坡的稳定性？

  考虑地震期间的边坡稳定性时，需要回答两个问题：

• 在特定的地震中边坡会被破坏吗？
• 如果受破坏，破坏特点如何？特别是滑坡体能运动多远，速度多快？

  由于滑坡体的传播速度和距离极其复杂，迄今为止很少有科学研究从这两个方面来表征不稳定性。

5.1. 边坡特征描述

  无论引发失稳的原因为何，边坡稳定性研究都需包含如下内容：

• 深入地质研究，旨在确定现场的几何结构、地层和地质历史（图7[16]）。
• 描述土壤渗流的水文地质研究，例如用压力计测量钻孔中地下水位的变化，因为水通常会促进边坡失稳。
• 通过力学试验表征构成边坡物质的力学特性，通常在实验室进行，旨在确定抗剪强度，即黏聚力和内摩擦角。

  基于上述研究内容获得的参数，可利用计算软件比较边坡失稳的下滑力和阻止滑动的抗滑力来确定边坡稳定与否，即确定边坡的静态安全系数。该系数大于1时（边坡抗滑力大于下滑力），边坡是稳定的。系数小于1时，边坡不稳定。而该系数越接近1（例如，对过去多次发生地震的边坡而言），即使中等强度的地震也会导致边坡失稳。

  为分析地震期间边坡的稳定性，还需刻画潜在地震的特征：这是地震期间评估边坡稳定性的基础。

5.2. 地震作用下分析边坡稳定性的方法

  目前主要有三类方法。

  最简单的方法是太沙基准静态法[17]，即计算边坡的准静态安全系数。该方法中，地震可以表征为一个恒定的体积力（volume force），其模量与不稳定坡体的重量成正比。该力可拆分为水平分量F_h和垂直分量F_v（图8）。在实践中，因为对边坡稳定性最具破坏性的力是水平作用力或剪切作用力，因而常忽略F_v。该方法可以判断边坡在地震中稳定（准静态安全系数大于1）还是不稳定（准静态安全系数小于1）。

  该方法的缺点是不提供关于边坡破坏概率或地震动强度的信息。

  因此，建议将其使用限制在地震作用下对边坡稳定性的初步分析，例如在人迹罕至或设备损失风险较低的地区。

  第二类方法是纽马克（Newmark）在1965年引入的“永久位移计算法”[18]（图9）。其步骤包括：根据加速度图，估计地震作用下假定刚性块体在支撑地面上滑动时可能发生的位移，以此表征地表移动情况[19]；通过数学计算获得不稳定块体的位移：将准静态安全系数等于1的加速度定义为“临界加速度”值，一段时间内所有超过该值的加速度峰都进行两次积分（图9）。一旦计算出该位移，专家就能确定边坡材料能否承受该位移而不导致边坡完全破坏。

  该方法的缺点：

• 未考虑地震作用不稳定体的内部变形（第4.2.段），因而仅适用于研究地震作用下变形可能较小的刚性材料中相对浅层的滑坡。
• 未考虑地震作用下土壤中可能产生的孔隙水压力（第4.3段）。
• 未考虑场地效应（第 4.1 节）。

  第三类方法为应力–应变法，强烈建议将其用于研究居民较多和（或）设施损失风险较高，且数据质量足够的特定场地。这类方法需要使用二维或三维计算软件评估边坡中给定地震引起的位移，软件集成了场地的所有复杂特质（不同材料性质、地形地貌等）和地震。图10为某个数值模拟案例，模拟了多次地震、反复荷载作用后边坡逐渐破裂的过程。

6. 地震引发的地表移动研究面临的挑战？

  为更好理解相关问题，目前许多课题都在研究地震引起的地表移动，内容包括：

• 水的影响。目前人们对水对地震引起地表移动的作用知之甚少。
• 三维效果（或3D效果）。由于边坡失稳现象极其复杂，通常以滑坡失稳传播方向上的二维（2D）截面作近似分析。这种近似对于人工边坡较为合理，但自然边坡的近似模拟往往与实际相差甚远。需要在计算机软件中输入描述边坡所有维度（3D）的信息，以实现对地震波传播方向及其对边坡稳定性影响的研究。
• 地震引发的地表移动预测。为预防自然灾害，必须预测地震引发的下一次移动将在何时何地发生，这意味着需要知道边坡某点的稳定状态，而边坡的稳定状态要追溯至其整个“荷载”历史，因而难以明悉。我们还需要预测可能影响某一特定地点的地震类型，以及该地震何时发生。最后，还需要量化边坡失稳的“强度”，即描述滑坡体的体积或不稳定体移动时可能覆盖的范围、滑坡体移动的速度和距离。由于地震引起的地表移动很复杂，不同的地震产生的结果可能不同，因此难以估量上述参数。

  现场观测还表明，地震也可能长期影响边坡稳定性：例如，1999年9月21日台湾地震后的几年里，季风引起的地表移动有所增加[22]。这是因为地震后边坡变得脆弱，更容易因强降雨事件失稳。

7. 要点总结

• 5级以上的地震很可能引发地表移动，而这类地表移动会造成人员伤亡和地震相关损失。
• 地震引发的地表移动类型众多，现象通常非常复杂。
• 要降低此类危害带来的风险，一方面要进一步理解引发地表移动的机制，另一方面要进一步认识控制滑坡体移动速度和移动距离的机制。

 

---

参考资料及说明

封面图片：2016年4月发生地震后，日本小木松市发生的一场滑坡。[图片©共同社/路透社，来源网络相簿（Flickr）,公共领域]

[1] Bird J.J. & Bommer J.J. (2004). Earthquake losses due to ground failure. Engineering Geology, 75 (2), p. 147-179. < https://doi.org/10.1016/j.enggeo.2004.05.006 >.

[2] Keefer D.K. (1984). Landslides caused by earthquakes. Bulletin of the seismological society of America, 95, pp. 406-421.

[3] 震级为地震学中表示测量地震释放能量的地震矩。地震释放的能量越多，震级越高：震级增加1相当于地震释放能量增加30倍。

[4] Gishig V., Preisig G. & Eberhardt E. (2016). Numerical Investigation of Seismically Induced Rock Mass Fatigue as a Mechanism Contributing to the Progressive Failure of Deep-Seated Landslides. Rock Mechanics and Rock Engineering. 49: 2457. < https://doi.org/10.1007/s00603-015-0821-z >.

[5] Dome J G., Bourdeau C., Lenti L., Martino S. & Pluta K. (2017). Mean landslide geometries inferred from a global database of earthquake- and non-earthquake-triggered landslides. Italian Journal of Engineering Geology and Environment, 2.https://doi.org/10.4408/IJEGE.2017-02.O-05.

[6] (https://www.ngdc.noaa.gov/hazardimages/picture/show/785)

[7] Lliboutry L. A. (1975). The Yungay disaster (Peru). Snow and Ice-Symposium-Neiges et Glaces. (Proceedings of the Moscow Symposium, August 1971; Proceedings of the Moscow Colloquium, August 1971): IAHS-AISH Publ. No. 104.

[8] Gorum T., van Westen C. J., Korup O., van der Meijde M., Fan X. & van der Meer F. D. (2013). Complex rupture mechanism and topography control symmetry of mass-wasting pattern, 2010 Haiti earthquake Geomorphology184, p. 127-138. < https://doi.org/10.1016/j.geomorph.2012.11.027 >.

[9] Petley D., Dunning S., Rosser N. & Kausar A.B. (2006). Incipient landslides in the Jhelum Valley, Pakistan following the 8th October 2005 earthquake. SAARC Workshop on Landslide Risk Management in South Asia.

[10] Yin, Y., Wang, F. & Sun, P. (2009). Landslide hazards triggered by the 2008 Wenchuan earthquake, Sichuan, China. Landslides 6, 139-151.

[11] Contains modified Copernicus Sentinel data {{{year}}}} [CC BY 3.0 (https://creativecommons.org/licenses/by/3.0) or Attribution], via Wikimedia Commons

[12] GeoNet, Geological hazard information for New Zealand

[13] 摩尔-库仑定律：该数学定律在岩土工程中用于以近似方式描述大多数长期承受剪切应力（称为剪切破坏）的土体行为。该定律有两个参数：黏聚力（即零法向应力下土体的抗剪强度）和内摩擦角（堆积土体相对于水平面自然形成的角度）。

[14] Miller P., Hovius N. & Haines J.A. (2008). Topographic site effects and the location of earthquake induced landslides. Earth and Planetary Science Letters. 275, 221–232.

[15] Bozzano F., Lenti L., Martino M., Paciello A. & Scarascia Mugnozza G. (2008). Self-excitation process due to local seismicamplification responsible for the reactivation of the Salcito landslide (Italy) on 31 October 2002. Journal of Geophysical Research 113, B10312. <https://doi.org/10.1029/2007JB005309>>.

[16] Bourdeau C., Lenti L., Martino S., Oguz O., Yalcinkaya E., Bigarre P. & Coccia A S. (2017). Comprehensive analysis of the local seismic response in the complex Büyükçekmece landslide area (Turkey) by engineering-geological and numerical modelling. Engineering Geology 218, 90-106. < http://dx.doi.org/10.1016/j.enggeo.2017.01.005 >.

[17] Terzaghi K. (1950). Mechanism of Landslides. Geological Society of America, Harvard University Press, 41 p.

[18] Newmark N.M. (1965). Effects of earthquakes on dams and embankments. Geotechnics 15, 139-160.

[19] Accelerogram: A signal giving the acceleration of the ground as a function of time is called an “accelerogram”.

[20] Cundall P.A. & HART R.D. (1992). Numerical modelling of discontinua. Eng. Comp. 9:101-113

[21] (2011). UDEC-Universal Distinct Element Code, Version 5.0. User’s Manual. Itasca Consulting Group, Inc. In Minneapolis

[22] Hovius N., Meunier P., Lin C.-W., Hongey C., Chen Y.-G., Dadson S., Horng M.-J. & Lines M. (2011). Prolonged seismically induced erosion and the mass balance of a large earthquake. Earth and Planetary Science Letters 304, 347-355