Formation mechanism and motion processes of the Aizigou giant paleolandslide，Jinshajiang river
-
摘要:
通过详细的野外调查，并结合遥感解译、室内试验以及数值模拟等手段，对矮子沟巨型古滑坡的基本特征、形成机制及运动演化过程进行了深入研究。矮子沟古滑坡的形成条件为：滑坡剪出口与坡脚之间存在巨大的高差，为滑坡的形成创造了良好的临空条件；顺向岸坡结构以及坡体内发育的多组控制性结构面是滑坡发生的结构基础；玄武岩系中的凝灰岩软弱夹层削弱了岩体的完整性，地表水及地下水长期入渗，水的软化作用降低了软弱夹层的抗剪强度；地震作用是造成岩体最终滑动失稳的关键因素。该滑坡的动力学过程可划分为四个阶段：① 启程活动阶段。斜坡地形效应使得地震波在斜坡上部表现出异常放大现象，当短时间内积聚的振动能量超过岩土体的强度时，易形成高位滑坡，滑坡的变形破坏机制为拉裂-滑移；② 近程活动阶段。近3.82×108 m3的滑坡物质高位高速下滑，与矮子沟右岸坡体发生猛烈碰撞后进一步碎裂解体；③ 高速远程碎屑流阶段。碎屑流继续沿矮子沟高速运动约3 km；④ 堆积堵江阶段。滑坡物质最终形成体积为2.73×108 m3的巨型堰塞坝，堵塞金沙江并形成堰塞湖，其回水区域向上游延续75 km至格勒一带。对金沙江矮子沟口上游河段采集到的堰塞湖沉积物进行光释光测年，结果指示该古滑坡发生于距今2.5万年左右。
Abstract:The Emeishan basalt group widely distributed in southwestern China is often considered as an ideal place for building large hydropower dams due to its deep canyon landforms. Historically, the large-scale high-position and long-distance landslides in Emeishan basalt have caused a large number of casualties, property damage and far-reaching environmental effects. Currently, there is no deep and systematic research for the development process of such landslides over the world, and it is still difficult to objectively access the risk of large-scale high-position landslides in southwestern China. Therefore, the study on the formation mechanism of large-scale high-position and long-distance landslides in Emeishan basalt has significant scientific and practical implication. The Aizigou paleolandslide, which represents a typical large basalt high-speed remote landslide, is located in the Liucheng section of the Lower Jinshajiang river, southwestern Sichuan Province, China. Geological investigation and interpretation of remote sensing imagery, in addition to experiments, numerical simulation, and geochronological dating were used to determine the formation mechanism and motion process of the landslide. Four factors were critical: a substantial height difference between the landslide shearing surface and the foot of the slope, multiple structural defects within the rock mass, a tuff intercalation in the basaltic series with long-term softening due to surface water and groundwater, and seismic activity. The dynamic behavior of the landslide is divided into four stages as follows. ① Seismic waves in the upper part of the slope were significantly amplified by the topography of the mountain. When the accumulation of vibrational energy over a short time period exceeded the strength of the rock mass, the landslide was formed easily. Accordingly, the deformation and failure mechanism was sliding and fracturing. ② Landslide materials with a volume of approximately 3.82×108 m3 slid down from a high position at a high velocity, after which they disintegrated into a debris flow after colliding with the mountainside on the right bank of Aizigou valley. ③ The high-speed debris flow moved approximately 3 km downstream within Aizigou valley. ④ Landslide materials blocked the Jinshajiang river, thereby forming a large barrier dam with a volume of 2.73×108 m3 and establishing a dammed lake behind the barrier dam; its backwater region reached 75 km upstream. An analysis of the landslide-dammed lake deposit samples via optically stimulated luminescence reveals that the landslide occurred approximately 25 thousand years ago.
-
引言
遥感技术因其快速获取大范围区域影像数据的优势，早在二十世纪八十年代就用于地震热异常的分析研究，并且取得了诸多成果（Tramutoli et al，2001；Choudhury et al，2006；Yao，Qiang，2012），震前热异常也在诸多震例中得到验证（吕月琳等，2009）。强祖基等（1990）通过多个震例探讨了临震卫星热红外异常是否与地表增温异常相关，其结果显示临震出现异常增温是普遍现象。吴立新等（2004）对岩石破裂的热红外辐射规律开展了试验研究，在探索地震热异常的形成机理方面取得了进展。
地震热异常具有强干扰、弱信号的特点，即与温度有关的气象信息容易受到多种复杂因素的干扰，而与地震有关的热异常较为微弱，因此需要排除强烈干扰背景以有效地提取热异常信息（李继业等，2011）。吴立新等（2017）引入热红外背景场，以热异常定量指标为判断依据，排除气象干扰构建较为稳定的背景场，对卫星影像观测值与背景值的偏差程度予以分析，发展了多种方法用于热异常提取。近年来，差值法、涡度法、小波功率谱法等多种方法广泛应用于地震热异常的提取分析，取得了诸多成果（陈顺云等，2006；郭晓等，2010；张元生等，2010；康春丽等，2011；刘军等，2014；魏从信等，2018）。陈顺云等（2006）基于NOAA气象卫星的甚高分辨率扫描辐射计（Advanced Very High Resolution Radiometer，缩写为AVHRR）陆地数据集，运用小波理论将地表热红外辐射分解为三个频率成分，结合气象资料和地震活动资料进行综合分析，其结果表明低频成分的热异常更接近于断层活动所引起的热异常信息，低频成分的高温异常区与构造特征、地震活动密切相关，因此低频成分热异常可为断层或地震活动提供一定的指示信息。马未宇等（2006）关于美国国家环境预报中心（National Centers for Environmental Prediction，缩写为NCEP ）提供的再分析数据温度异常图像与天体引潮力附加构造应力的联合分析结果显示，西藏的三次地震均发生于天体引力潮附加构造应力变化达到相对稳定值的末端。康春丽等（2011）基于NOAA-18提供的分辨率为1°×1°网格化的夜间长波辐射（outgoing longwave radiation，缩写为OLR）数据，使用涡度算法分析了四次大地震震前的热异常，结果显示震前一到两天热辐射异常明显增强，震后逐渐消失。
卫星数据鲁棒分析技术（robust satellite techniques，缩写为RST）自Tramutoli （1998）提出以来，被用于大量的震例研究，相关结果显示多个震例前均有热红外异常现象出现（Filizzola et al，2004；Tramutoli et al，2013；Genzano et al，2015；宋冬梅等，2016）。Ma等（2012）基于欧洲第二代气象卫星（Meteosat Second Generation，缩写为MSG）红外波段数据运用RST算法验证并提取了意大利阿布鲁佐大区MS6.3地震前的热异常。宋冬梅等（2016）基于地表温度（land surface temperature，缩写为LST）数据使用RST算法再次对汶川地震进行分析，结果显示所提取的热异常与地质断裂更加吻合。Eleftheriou等（2015）也使用RST算法提取了2004—2013年希腊地区的长时序地震热异常，并对其特点予以分析。
上述研究用于热异常提取的遥感数据通常为夜间某一时刻或夜间均值数据，但是某时刻数据并不能准确地反映夜间热环境的整体情况，且均值数据也有可能导致热信息被湮没。为更加准确地判断热异常是否与地震相关，本文拟综合静止气象卫星的高时间分辨率特性，改进热异常与地震相关性的判别方法，运用RST算法提取夜间均值数据热异常后，进一步提取夜间多时刻热异常，以期通过分析多时刻热异常的时空分布来判断热异常是否可能与地震相关，进而为地震热异常判定提供参考。
1. 数据源选取及热异常提取方法
1.1 数据源
本文采用MSG搭载的旋转增强型可见光及红外成像光谱仪（Spining Enhanced Visible and Infrared Imager，缩写为SEVIRI） 10.8 μm波段数据。该型成像仪具有4个可见光—近红外通道和8个红外通道，其数据具有较高的时空分辨率，在空间上，中红外通道分辨率为3 km，比第一代的5 km有较大提高，且图像量化等级提高到10 bit （徐清，2013），极大地提升了数据质量。
1.2 去云
在热异常提取过程中，由于经常受到云层的干扰，可能无法提取出热异常。为去除遥感数据中云层带来的干扰，本文采用一种云检测方法（Cuomo et al，2004），即基于历年同期数据，通过设置迭代相应阈值，识别出云像元，具体公式如下：
${S_{{\rm{IR}}}}\!\!\!\!{\text{（}}\!{r{\text{，}}\!\!\!t} \!{\text{）}}\!\!\!\! {\text{＝}} \frac{{IR(r{\text{，}}\!\!\!t) {\text{－}} {\mu _{{\rm{IR}}}}\!\!\!\!{\text{（}}\!r \!{\text{）}}\!\!\!\!}}{{{\sigma _{{\rm{IR}}}}\!\!\!\!{\text{（}}\!r \!{\text{）}}\!\!\!\!}}{\text{，}}\!\!\!$
(1) $IR(r{\text{，}}\!\!\!t) {\text{－}} {\mu _{{\rm{IR}}}}\!\!\!\!{\text{（}}\!r \!{\text{）}}\!\!\!\! {\text{≤}} N{\text{，}}\!\!\!$
(2) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
1.3 RST算法
RST算法是一种基于多时相数据的分析方法，其优势在于突出时空域上相对空间域的异常，最初被用于云检测、火山、森林火灾监测等方面（Tramutoli，1998；Genzano et al，2009；Tramutoli et al，2013），后经不断改进应用于地震热红外监测。该方法将地震引起的热红外异常看作信号，其它因素引起的异常看作噪声，引入地震热异常信噪比的概念（孙珂，2017）。热异常稳健估计（robust estimator of thermal infrared anomalies，缩写为RETIRA）值R （Tramutoli et al，2001）是以多年同期的遥感数据为基础，运用RST算法构建去除噪声的背景场，提取相对于背景场的偏离值，表达式如下：
$R {\text{＝}} \frac{{\Delta T\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y{\text{，}}\!\!\!t\!{\text{）}}\!\!\!\! {\text{－}} {\mu _{\Delta T}}\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y\!{\text{）}}\!\!\!\!}}{{{\sigma _{\Delta T}}\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y\!{\text{）}}\!\!\!\!}}{\text{，}}\!\!\!$
(3) $\Delta T\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y{\text{，}}\!\!\!t\!{\text{）}}\!\!\!\! {\text{＝}} T\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y{\text{，}}\!\!\!t\!{\text{）}}\!\!\!\! {\text{－}} T\!\!\!\!{\text{（}}\!t\!{\text{）}}\!\!\!\!{\text{，}}$
(4) 式中，T （x，y，t）为（x，y）处的像元值，T （t）为影像区域像元均值，ΔT （x，y，t）为热异常待提取点位的像元值相对于区域均值的差值，μΔT（x，y）为ΔT （x，y，t）多年相同位置、相同时刻的平均值，σΔT（x，y）为对应的标准差。
由于在构建背景场所用到的数据存在干扰噪声，需进一步去除这些噪声。通过
$A\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y{\text{，}}\!\!\!t\!{\text{）}}\!\!\!\! {\text{＝}} \left\{ {\begin{array}{*{20}{c}} \!\!\!1&{{\rm{|}}T\!\!\!\!{\text{（}}\! {x{\text{，}}\!\!\!y{\text{，}}\!\!\!t} \!{\text{）}}\!\!\!\!{\rm{ }} {\text{－}} {\mu _T}\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y{\text{，}}\!\!\!t\!{\text{）}}\!\!\!\!| {\text{＜}} k{\delta _T}\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y{\text{，}}\!\!\!t\!{\text{）}}\!\!\!\!}{\text{，}}\\ \!\!\!0&{\text{其它}}{\text{，}} \end{array}} \right.$
(5) ${\mu _{\Delta T}}\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y\!{\text{）}}\!\!\!\! {\text{＝}} \frac{{{{\displaystyle\sum\limits_{\forall t \in T} {[\Delta T^2\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y{\text{，}}\!\!\!t\!{\text{）}}\!\!\!\! \cdot A^2\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y{\text{，}}\!\!\!t\!{\text{）}}\!\!\!\!]} }}}}{{\displaystyle\sum\limits_{\forall t \in T} {A\!\!\!\!{\text{（}}\!x{\text{，}}\!\!\!y{\text{，}}\!\!\!t\!{\text{）}}\!\!\!\!} }}{\text{，}}$
(6) 对遥感影像进行处理，构建去除噪声的差值背景场。式（5）和（6）中，A（x，y，t）用来剔除构建背景场的异常值，μΔT（x，y）为剔除背景噪声后的背景场（Zhang，Meng，2019）。构建稳定背景场的具体要求为：对于该点位的时间序列，与均值之差的绝对值大于k （k≥2）倍标准差的A值不可用于构建背景场。将A（x，y，t）赋值为0，迭代直到数值与均值之差均小于k倍中误差（Genzano et al，2009）。由式（6）计算相应的背景场μΔT（x，y），进而得到R值。
2. 多时相地震热异常判定
人类活动及极端天气、大范围洪水、森林火灾等多种因素的影响均可能导致热异常，这使得所提取的地震热异常可能具有假阳性（Filizzola et al，2004），因此需通过筛选来减少非地震因素所造成的干扰。Eleftheriou等（2015）提出的热异常判定条件具体如下：
1） 选取相对偏离阈值为4，根据数据特性选取R≥4可有效检测出热异常；
2） 在时间分布上，第一次出现热异常后七日内应至少再出现一次，且两次热异常的分布位置相同或相邻；
3） 在空间分布上，热异常分布应当是聚集的，检测出的热异常在1°×1°的范围内应当至少覆盖150 km2的面积。
Eleftheriou等（2015）认为满足上述三条即可视为可能与地震相关的热异常。虽然基于影像均值或某时刻数据所获得的热异常提取结果有一定的代表性，可反映夜间一段时间内热辐射的基本情况，但是由于云层及其它干扰因素的不断变化，影像均值可能丢失部分热信息，难以反映夜间热异常的时空变化特征。因此，本文结合Eleftheriou等（2015）的三条热异常判定条件，增加判定条件4）：进一步提取分析夜间多时相数据的热异常，在排除云层干扰的情况下，由多时相数据所提取热异常的时空分布连续。这样，热异常在满足前三条热异常筛选条件的同时，仍需满足条件4），以排除短期天气因素所导致的干扰，满足该条件即可有效地排除异常时空分布强度不足、不连续及空间分布零散等不符合地震热异常的情况。
本文首先选取夜间0时、1时、2时的时刻数据，求均值后提取热异常，对于存在热异常的情况，再选取当日夜间20:30至次日05:00时段每隔30分钟数据，分析夜间18个时刻的遥感数据热异常，以探究其在夜间多个时段的时空演化特征，进而判定热异常与地震的相关性。
3. 热异常提取结果分析
3.1 数据均值热异常
本文基于RST算法探究意大利地区震例的震前热异常，使用2009—2018年SEVIRI遥感10.8 μm数据构建背景场。首先处理夜间0时、1时、2时的均值数据，提取三次热异常，分别为2015年5月6日和5月7日出现在意大利半岛南部和西西里岛北部的热异常、2016年1月9日和12日出现在意大利半岛中部的热异常以及2017年3月21日和23日出现在意大利半岛中部的热异常，其均值热异常如图1所示。根据研究区地形以及热异常序列的分布，本文认为未来地震可能发生于热异常中心附近。事实上，2015年5月6日和5月7日热异常出现后的第三天附近区域发生了MS4.7地震，2016年1月9日和12日热异常出现后的第七天附近区域发生了MS4.4地震，2017年3月21日和23日热异常出现后30日内附近区域未发生地震。
![]()
图 1 意大利三次夜间均值热异常分布（a） 2015年5月6日和7日；（b） 2016年1月9日和12日；（c） 2017年3月21日和23日Figure 1. Three nightly mean thermal anomalies in Italy（a） On May 6 and 7，2015；（b） On January 9 and 12，2016；（c） On March 21 and 23，2017如图1a所示：2015年5月6日在意大利半岛南部和西西里岛北部出现了热异常，7日意大利南部半岛部分区域出现了云层，未能检测出热异常；5月6日热异常分布较7日的分布范围更广。根据热异常的判定条件，西西里岛北部等区域出现的单日热异常不属于有效热异常。由图1b可见：2016年1月9日和1月12日在意大利半岛中部出现了热异常，相关区域的热异常空间分布较为连续；1月12日在意大利半岛中北部和南部出现了热异常。这几处热异常不连续，因此不认为其为有效热异常序列。由图1c可见：2017年3月21日和23日意大利半岛中部区域出现了热异常，两日热异常序列空间分布较为连续，也存在部分区域出现零散热异常的情况。
3.2 多时相数据热异常提取
基于RST算法进一步提取夜间均值多时刻影像热异常，本文选取从当日20:30至次日05:00 （时间均为UTC）每隔30分钟的遥感影像作为基础数据，分析不同时刻的热异常分布特征，并根据文中的判定准则判断该热异常是否可能与地震相关。
通过分析意大利2015年5月6日和5月7日多时刻数据，得到热异常的时空分布序列，如图2所示。可见：5月6日各时刻均出现了热异常，本次热异常序列限于意大利半岛南部和西西里岛部分区域；5月6日凌晨前后多时刻热异常随着时间推移呈现逐渐增多的扩展趋势，其中意大利半岛南部和西西里岛北部的热异常均不断扩展，在震中附近的热异常分布呈现不断增多的趋势；5月7日凌晨前后多时刻热异常时空分布较为连续；由于云层的干扰，5月6日20时30分至5月7日0时0分这期间的7个时刻未出现显著热异常，然而随着时间推移，云层逐渐消散，相关区域热异常逐步显现、增多、扩展，因此认为此次热异常存在的不连续情况是由于云层干扰所致，根据文中的判定方法，判定此次热异常可能与地震相关。
![]()
图 2 意大利地区2015年5月6日和7日多时刻热异常时空分布Figure 2. Spatio-temporal distribution of multi-time thermal anomalies in Italy on May 6 and 7，2015基于RST算法提取意大利2016年1月9日和12日的夜间多时相热异常，得到热异常的时空分布序列，如图3所示。可见：1月9日各时刻相关区域均出现了显著热异常，时空连续性较好，云层的变化情况清晰可见，相关区域附近的云量逐步增多，部分出现热异常的区域逐步被云层覆盖，导致无法连续提取出热异常；1月12日意大利中北部出现了较大区域的热异常，根据热异常的分布特性，此部分热异常与9日热异常的空间距离较远且不连续，故认为此部分热异常不属于有效热异常；1月9日和12日夜间意大利半岛中部均出现了连续热异常，时空分布具有较好的连续性，根据文中的判定方法，判定此次热异常可能与地震相关。
![]()
图 3 意大利地区2016年1月9日和12日多时刻热异常时空分布Figure 3. Spatio-temporal distribution of the multi-time thermal anomalies in Italy on 9 and 12 January 2016通过分析意大利2017年3月21日和3月23日多时刻数据得到夜间多时相热异常的提取结果，如图4所示。3月20日和22日在21:30和22:00时刻均出现数据缺失，导致这两个时刻无法提取热异常，根据两个时刻前后热异常的提取结果，推断相关区域可能未出现热异常。从图4可见：3月21日凌晨前后多数时刻未出现热异常，仅01:00—04:00出现热异常，而后热异常再次消失，热异常未出现时刻未有云层出现，可排除存在云层的干扰，推测热异常可能由短期天气变化导致；3月23日凌晨前后各时刻数据多数时刻未提取出热异常，仅02:30—05:00少数时刻出现热异常，加上云层在相关区域分布较少，在一定程度上可排除云层对于热异常的干扰，但相关区域仍未出现连续热异常，根据文中的判定方法，认为此次热异常与地震无关。
![]()
图 4 意大利地区2017年3月21日和23日多时刻热异常时空分布Figure 4. Spatio-temporal distribution of thermal anomalies in Italy at multiple times on March 21 and 23，20173.3 多时相热异常分析与统计
根据三次夜间多时相热异常的时空分布可知，云层对于热异常提取的影响较大，云层的扩散与消失可能导致无法提取热异常或引起热异常时空分布变化较大。2015年5月7日和5月9日、2016年1月9和1月12日的热异常提取结果均受到云层干扰，云层的移动可能会引起热异常的空间分布产生变化。由2016年1月9日多时刻热异常时空序列的提取结果（图3），可明显看出1月9日1时起震中北部方向逐渐有云层出现，导致北部区域的热异常逐渐消失，而震中西北部地区逐渐出现了热异常；在2015年5月6日和5月9日、1月12日相关区域热异常（图2）通常因云层边缘的移动而出现变化，故认为云层在相关区域的时空变化对于热异常提取具有重要影响。
对夜间多时相云层面积占成像区域面积的百分比进行统计，结果如图5所示，可见：2015年5月6日凌晨前后云层缓慢增多，5月7日凌晨前后云层所占比例逐渐下降；2016年1月9日凌晨前后多时刻云层呈现逐步扩展增多的趋势，12日凌晨前后多时刻云层分布相对平稳；2017年3月21日和23日凌晨前后云层均呈现逐步增多，但基本保持稳定。
![]()
图 5 意大利三次多时刻云层及热异常面积统计（a） 2015年5月6日和7日；（b） 2016年1月9日和12日；（c） 2017年3月21日和23日Figure 5. Statistics of three multi-temporal cloud cover and thermal anomalies in Italy（a） On May 6 and May 7，2015；（b） On January 9 and January 12，2016；（c） On March 21 and 23，2017通过统计意大利三次地震热异常的空间分布面积，可进一步了解夜间热异常的分布变化，探究其分布形态。如图5a所示，自5月6日凌晨起各时刻热异常基本呈现逐步增多趋势，其中5月6日3时热异常分布范围最大，面积达15 030 km2；5月7日凌晨前后多时刻热异常面积整体较前一日明显减小；5月6日20:30至24:00时未出现显著热异常。由图2可知，5月6日凌晨前后，因有云层出现无法提取热异常，之后云层逐渐消散，热异常面积逐渐增大。如图5b所示，2016年1月9日凌晨前后各时刻均出现显著热异常，多时刻热异常呈振荡增多趋势；1月8日23时热异常分布范围最大，面积达2 567 km2，随后分布范围逐步缩小，至1月9日1时开始逐步增大，02:30后热异常面积减小，并趋于稳定；1月12日凌晨前后多时刻数据均出现了显著热异常，热异常面积整体呈现先增大后逐渐减小的趋势，在1月12日00:30降至最小值，为486 km2，之后逐渐增大，在05:00达到最大值，为6 228 km2。如图5c所示，3月21日和23日凌晨前后多数时刻未出现热异常，两日在01:30前均未出现显著热异常，其中3月21日03:00热异常分布范围为当日最大，面积为558 km2；3月23日04:00面积达到当日最大值，为459 km2。
4. 讨论与结论
针对热异常与地震相关性的准确判断，本文对Eleftheriou等（2015）提出的热异常与地震相关性的判定条件予以改进，深入地探究基于夜间多时相数据的热异常与地震的相关性，并基于RST算法提取了意大利地区2015—2017年期间的三次热异常，分析了这些热异常的分布及其变化特征。
1） 根据影像均值或某一时刻的影像值判断热异常是否与地震相关可能存在误判情形，具有一定的局限性，本文通过改进热异常与地震相关性的判定条件，探究研究区夜间多时刻连续热环境的变化，可清晰地观察到热异常在夜间的时空变化特征。
2） 根据改进后的热异常与地震相关性判定条件，探究夜间多时相数据的热异常，结果显示：两次夜间多时相热异常序列较为连续，认为可能与地震相关；一次夜间多时相热异常在排除不存在云层干扰后，夜间多数时刻热异常不具有连续性，仅少数时刻存在热异常，认为该次热异常可能与地震无关。本文根据三次热异常后是否发生地震，验证了改进判定条件，提升了热异常与地震相关性判定的准确性。
3） 三次夜间多时刻热异常的时空变化表明，云层对于热异常提取的影响较大。热异常可能因云层的干扰而无法被提取，云层的扩散与消散都可能对热异常提取造成影响，同时云层分布比例与热异常分布面积并无显著相关性。
综上所述，本文基于夜间多时相数据探究了热异常与地震的相关性，进行了意大利地区三次热异常与地震的相关性分析，初步证实了该方法的有效性。为进一步探究夜间多时相数据提取的热异常与地震相关性判定方法的区域适用性，仍需大量震例作进一步研究。
-
![]()
图 1 研究区地质图及构造背景图
F1：则木河断裂；F2：越西断裂；F3：四开—交际河断裂；F4：布拖断裂；F5：茂租断裂；F6：莲峰—巧家断裂；F7：小江断裂；F8：普渡河—大桥河断裂；F9：宁南—会理断裂。右下插图中：① 龙门山断裂带；② 则木河断裂带；③ 小江断裂带；④ 红河断裂带；⑤ 小金江断裂带
Figure 1. The geological map and tectonic settings of the study region
F1：Zemuhe fault；F2：Yuexi fault；F3：Sikai-Jiaojihe fault；F4：Butuo fault；F5：Maozu fault；F6：Lianfeng-Qiaojia fault；F7：Xiaojiang fault；F8：Puduhe-Daqiaohe fault；F9：Ningnan-Huili fault. Illustration in the lower right corner：① Longmenshan fault zone；② Zemuhe fault zone；③ Xiaojiang fault zone；④ Honghe fault zone；⑤ Xiaojinjiang fault zone
![]()
图 3 矮子沟古滑坡平面图（上）和剖面图（下）
Ⅰ ：滑源区；Ⅱ ：高位高速下滑区；Ⅲ ：撞击碎裂区；Ⅳ ：流通区；Ⅴ ：堆积区
Figure 3. Geological map （upper） and cross section profile AB （lower） of the Aizigou paleolandslide
Ⅰ ：Landslide source area；Ⅱ ：High-position and high-speed sliding area；Ⅲ ：Impact fragmentation area；Ⅳ ：Transportation area；Ⅴ ：Accumulation area
![]()
图 4 撞击碎裂区内及其相邻地区出露的玄武岩基岩
B1区和B2区为观测点#24至#25段出露的碎裂状玄武岩基岩；B3区和B4区为观测点#26至#27段出露的碎裂状玄武岩基岩；B5区为撞击碎裂区相邻地区出露的相对完整的玄武岩基岩
Figure 4. Basaltic bedrocks exposed in the impact fragmentation area and its adjacent areas
The subareas B1 and B2 are the cataclastic basalt bedrock outcropped at the section from the observation point #24 to point #25，the subareas B3 and B4 are those at the section from the point #26 to point #27，and the subarea B5 is the relatively intact basaltic bedrock adjacent to the impact fragmentation area
![]()
图 5 滑源区对岸坡体上的滑坡堆积物
C1区和C2区为玄武岩大块石、孤石分布区，C3区和C4区为撞击碎裂区中下部分布的滑坡堆积物
Figure 5. Parts of the landslide deposits distributed on the bank opposite to the landslide source area
The subareas C1 and C2 are an area with large rocks and boulders，and the subareas C3 and C4 are landslide deposits exposed in the middle and lower parts of the impact fragmentation area
![]()
图 7 #11调查点处的滑坡堆积物
Figure 7. Landslide deposits in excavation at the observation point #11
![]()
图 8 矮子沟出口段北侧坡体上广泛分布的滑坡堆积物
（a−i） 调查点#12，#13和#14处的滑坡堆积物；（j，k） 调查点#18处的滑坡堆积物；（l） 山脊处分布的玄武岩大块石
Figure 8. Massive landslide deposits widely distributed across the northern slope of the mouth of Aizigou valley
（a−i） Landslide deposits at the observation points #12，#13，and #14；（j，k） Landslide deposits at the observation point #18；（l） Large basaltic boulders distributed on the ridge
![]()
图 9 金沙江右岸坡体上广泛分布的滑坡堆积物
（a，b） 调查点#2处的滑坡堆积物；（d，e） 调查点#3处的滑坡堆积物；（c，f） 调查点6#处已出现泥化现象的部分玄武岩块碎石；（g，h，j，k） 调查点#7处的滑坡堆积物；（i，l） 金沙江右岸坡体平台处分布的玄武岩大块石
Figure 9. Massive landslide deposits widely distributed on the right bank of Jinshajiang river
（a，b） Landslide deposits at observation point #2；（d，e） Landslide deposits at observation point #3；（c，f） The basaltic blocks showing severe argillitization at observation point #6；（g，h，j， k） Landslide deposits at observation point #7；（i，l） Large basaltic boulders distributed on a platform on the right bank of Jinshajiang river
![]()
图 10 金沙江两岸出露的古堰塞湖沉积物
（a） 古堰塞湖沉积物分布；（b） 新场出露的粉土层；（c） 葫芦口出露的完整的堰塞湖沉积物；（d） 金堂出露的粗砂层；（e） 格勒出露的砂砾石层
Figure 10. The ancient dammed lake sediments exposed along both banks of Jinshajiang river
（a） Distribution of the ancient dammed lake sediments；（b） Silt layers at Xinchang；（c） Dammed lake sediments at Hulukou；（d） Coarse sand layers at Jintang；（e） Sand gravel layers at Gele
![]()
图 12 矮子沟古滑坡的3D离散元数值模型
Figure 12. The 3D discrete element model of the Aizigou paleolandslide
![]()
图 13 边界条件及监测点A1−A5和B1−B6示意图
Figure 13. Boundary conditions and distribution of monitoring points A1−A5 and B1−B6 of 3DEC model
![]()
图 14 地震波水平向（a）、竖直向（b）加速度时程曲线
Figure 14. Acceleration time-history curves in horizontal （a） and vertical （b） directions
![]()
图 15 施加的地震动力荷载
（a） 水平向剪应力时程曲线；（b） 竖直向正应力时程曲线
Figure 15. Applied seismic dynamic loads
（a） Time history curve of horizontal shear stress；（b） Time history curve of vertical normal stress
![]()
图 17 监测点A1和A5的水平向（a）、竖直向（b）加速度时程曲线
Figure 17. Horizontal （a） and vertical （b） acceleration time-history curves of monitoring points A1 and A5
![]()
图 18 滑坡在不同阶段的运动状态
（a） 在启程活动阶段，滑坡体发生明显的变形和错动；（b） 锁固段岩体被剪断，滑坡体启动；（c） 在近程活动阶段，滑坡体高位高速下滑；（d） 滑坡体高速下滑过程中，与矮子沟右岸坡体发生猛烈碰撞后进一步碎裂解体转化为碎屑流；（e，f） 滑坡碎屑流继续沿矮子沟高速滑动，进入高速远程滑动阶段；（g） 滑坡最终堵塞金沙江河道，形成堰塞湖
Figure 18. Movement states of the landslide at different running periods in the numerical model
（a） In the set-out stage，the head scarp of the landslide ruptured with an obvious deformation；（b） The rock mass of the locked segment was sheared off and the body started sliding；（c） The landslide maintained a high speed after its violent start at a high altitude，and slid down with increasing velocity in the short-range stage；（d） Landslide materials disintegrated into a debris flow after colliding with the mountainside along the right bank of Aizigou valley；（e，f） The debris flow continued to move down along Aizigou valley at a high speed and entered the rapid and long-runout debris flow stage；（g） The landslide blocked the Jinshajiang river，and formed a lake
表 1 金沙江下游六城段河谷演化阶段（黄典等，2010；韩刚等，2011）
Table 1 Valley evolution stages in the Liucheng section of the lower Jinshajiang river （Huang et al，2010；Han et al，2011）
地貌类型 河拔高度/m 高程/m 地质年代 形成时代/ka Ⅳ级阶地 250 830 ${\rm{Q}}_{{\rm{p}}}^{2}$ 180 Ⅲ级阶地 140 720 ${\rm{Q} }_{ {\rm{p} } }^{2} $ 150 Ⅱ级阶地 100 680 ${\rm{Q} }_{ {\rm{p} } }^{3} $ 39 Ⅰ级阶地 20 600 ${\rm{Q} }_{ {\rm{h} } } $ 11 
下载: 导出CSV
表 2 滑坡堰塞湖沉积物样品的光释光测年数据
Table 2 Optically stimulated luminescence dating of the landslide-dammed lake deposit samples
实验室编号 采样地点 K含量 Th含量/10−6 U含量/10−6 质量含水量 年龄/ka 年龄误差/ka 17 124-1 六城镇 1.76% 8.66 1.78 2.22% 6.8 0.7 17 124-2 新场 2.92% 12.20 2.26 2.03% 10.0 0.1 17 124-3 棉纱村 2.85% 12.90 3.02 5.72% 15.7 1.6 17 124-4 葫芦口 0.25% 11.83 0.54 22.20% 20.3 1.5 17 124-5 葫芦口 1.74% 12.30 1.85 6.27% 20.5 1.7 17 124-6 葫芦口 1.73% 12.4 1.89 5.41% 25.7 2.9 17 124-7 大崇 1.50% 8.28 1.41 5.64% 20.3 2.2 17 124-8 金堂 1.43% 9.85 1.78 16.53% 12.7 1.4 17 124-9 蒙姑 1.63% 11.80 2.02 14.75% 13.3 1.2 17 124-10 格勒村 1.74% 7.28 1.68 8.33% 19.6 2.3
下载: 导出CSV
岩性 天然密度
/（kg·m−3）黏聚力
/MPa内摩擦角
/°体积模量
/GPa剪切模量
/GPa微风化玄武岩 2950 12.9 55 45.8 21.2 弱风化玄武岩 2860 9.8 45 31.6 17.3
下载: 导出CSV
结构面 法向刚度
/（GPa·m−1）切向刚度
/（GPa·m−1）内摩擦角
/°黏聚力
/MPa抗拉强度
/MPa节理面 26 19 20 1.2 0.03 滑面 9 6 15 0.5 0.01
下载: 导出CSV
表 5 监测点A1，A2，···，A5的峰值加速度及放大系数
Table 5 PGAs and their amplification coefficients of the monitoring points A1，A2，···，A5
监测点号 PGA/（m·s−2） PGA放大系数 水平向 竖直向 水平向 竖直向 A1 3.24 2.77 1.08 1.31 A2 4.68 3.95 1.56 1.87 A3 −7.02 −4.46 2.34 2.11 A4 −8.37 6.75 2.79 3.20 A5 −17.46 −14.24 5.82 6.75
下载: 导出CSV
-
程谦恭,胡厚田,胡广韬,彭建兵. 2000. 高速岩质滑坡临床弹冲与峰残强降复合启程加速动力学机理[J]. 岩石力学与工程学报,19(2):173–176. doi: 10.3321/j.issn:1000-6915.2000.02.010 Cheng Q G,Hu H T,Hu G T,Peng J B. 2000. A study of complex accelerated dynamics mechanism of highspeed landslide by elastic rocky impulse and peak-residual strength drop[J]. Chinese Journal of Rock Mechanics and Engineering,19(2):173–176 (in Chinese).
方健,尹小涛,周磊,严飞. 2017. 金沙江中游凝灰岩强度参数综合识别研究[J]. 应用力学学报,34(3):507–513. Fang J,Yin X T,Zhou L,Yan F. 2017. Comprehensive study on strength parameters of the tuff at Jinshajiang grand bridge site on Huali freeway[J]. Chinese Journal of Applied Mechanics,34(3):507–513 (in Chinese).
中国地震局. 2020. GB/T 17742—2020中国地震烈度表[S]. 北京: 中国标准出版社: 5–7. China Earthquake Administration. 2020. GB/T 17742−2020 The Chinese Seismic Intensity Scale[S]. Beijing: Standards Press of China: 5–7 （in Chinese）.
韩刚,赵其华,彭社琴. 2011. 白鹤滩水电站坝区深部破裂岩体地应力演化特征[J]. 岩土力学,32(增刊1):583–589. doi: 10.16285/j.rsm.2011.s1.107 Han G,Zhao Q H,Peng S Q. 2011. In-situ stress field evolution of deep fracture rock mass at dam area of Baihetan hydropower station[J]. Rock and Soil Mechanics,32(S1):583–589 (in Chinese).
韩刚. 2015. 一类不对称发育的深部变形破裂成因机理: 以白鹤滩水电站坝址区深部变形破裂为例[D]. 成都: 成都理工大学: 57–60. Han G. 2015. A Type of Asymmetric Distributed Deep-Seated Crack Formation Mechanism: A Case Study of Baihetan Hydropower Station[D]. Chengdu: Chengdu University of Technology: 57–60 （in Chinese）.
胡广韬. 1995. 滑坡动力学[M]. 北京: 地质出版社: 75–76. Hu G T. 1995. Landslide Dynamics[M]. Beijing: Geological Publishing House: 75–76 （in Chinese）.
胡卸文,黄润秋,施裕兵,吕小平,朱海勇,汪雪瑞. 2009. 唐家山滑坡堵江机制及堰塞坝溃坝模式分析[J]. 岩石力学与工程学报,28(1):181–189. doi: 10.3321/j.issn:1000-6915.2009.01.024 Hu X W,Huang R Q,Shi Y B,Lü X P,Zhu H Y,Wang X R. 2009. Analysis of blocking river mechanism of Tangjiashan landslide and dam-breaking mode of its barrier dam[J]. Chinese Journal of Rock Mechanics and Engineering,28(1):181–189 (in Chinese).
黄典,杨达源,李郎平,石安池,陆飞,葛兆帅,胥勤勉. 2010. 金沙江白鹤滩河段下切速率初步研究[J]. 第四纪研究,30(5):872–876. doi: 10.3969/j.issn.1001-7410.2010.05.03 Huang D,Yang D Y,Li L P,Shi A C,Lu F,Ge Z S,Xu Q M. 2010. Incision rates of the Baihetan section of the Jinsha river[J]. Quaternary Sciences,30(5):872–876 (in Chinese).
黄润秋. 2009. 汶川地震地质灾害研究[M]. 北京: 科学出版社: 233–235 Huang R Q. 2009. Geohazard Assessment of the Wenchuan Earthquake[M]. Beijing: Science Press: 233–235 （in Chinese）.
李国和,王思敬,孙承志. 2001. 金沙江水电开发区域工程地质环境综合评价[J]. 地球科学:中国地质大学学报,26(3):309–313. Li G H,Wang S J,Sun C Z. 2001. Comprehensive assessment of engineering geological environment in hydroelectric development zone,Jinshajiang River[J]. Earth Science: Journal of China University of Geosciences,26(3):309–313 (in Chinese).
李坪. 1993. 鲜水河—小江断裂带[M]. 北京: 地震出版社: 23–25. Li P. 1993. Xianshuihe and Xiaojiang Fault Zones[M]. Beijing: Seismological Press: 23–25 （in Chinese）.
李忠生. 2003. 国内外地震滑坡灾害研究综述[J]. 灾害学,18(4):64–70. doi: 10.3969/j.issn.1000-811X.2003.04.013 Li Z S. 2003. The state of the art of the research on seismic landslide hazard at home and abroad[J]. Journal of Catastrophology,18(4):64–70 (in Chinese).
王刚,王二七. 2005. 挤压造山带中的伸展构造及其成因:以滇中地区晚新生代构造为例[J]. 地震地质,27(2):188–199. doi: 10.3969/j.issn.0253-4967.2005.02.002 Wang G,Wang E Q. 2005. Extensional structures within the compressional orogenic belt and its mechanism:A case study for the Late Cenozoic deformation in Central Yunnan[J]. Seismology and Geology,27(2):188–199 (in Chinese).
文宝萍,王凡. 2021. 1965年烂泥沟滑坡前兆、高速远程运动及后期演化特征[J]. 水文地质工程地质,48(6):72–80. doi: 10.16030/j.cnki.issn.1000-3665.202108065 Wen B P,Wang F. 2021. Precursors and motion characteristics of the 1965 Lannigou rockslides and the subsequent evolution[J]. Hydrogeology and Engineering Geology,48(6):72–80 (in Chinese).
吴昊,裴向军,崔圣华,黄润秋,宋词. 2021. 强震山区滑坡发育分布的地形地质控制作用研究[J]. 岩石力学与工程学报,40(5):972–986. doi: 10.13722/j.cnki.jrme.2020.0944 Wu H,Pei X J,Cui S H,Hung R Q,Song C. 2021. Study of topographic and geological controls on landslide development and distribution within mountainous regions influenced by strong earthquakes[J]. Chinese Journal of Rock Mechanics and Engineering,40(5):972–986 (in Chinese).
邢爱国,胡厚田,杨明. 2002. 大型高速滑坡滑动过程中摩擦特性的试验研究[J]. 岩石力学与工程学报,21(4):522–525. doi: 10.3321/j.issn:1000-6915.2002.04.014 Xing A G,Hu H T,Yang M. 2002. Testing study on frictional characteristic of large-scale and high-speed landslide during sliding[J]. Chinese Journal of Rock Mechanics and Engineering,21(4):522–525 (in Chinese).
邢爱国,陈龙珠,陈明中. 2004. 岩石力学特性分析与高速滑坡启动速度预测[J]. 岩石力学与工程学报,23(10):1654–1657. doi: 10.3321/j.issn:1000-6915.2004.10.012 Xing A G,Chen L Z,Chen M Z. 2004. Analysis on mechanical properties of rock and prediction of muzzle velocity of high-speed landslide[J]. Chinese Journal of Rock Mechanics and Engineering,23(10):1654–1657 (in Chinese).
胥良. 2004. 金沙江白鹤滩水电站金江滑坡成因机制及稳定性研究[D]. 成都: 成都理工大学: 72–75. Xu L. 2004. A Study on Genetic Mechanism and Stability of Jinjiang Landslide at the Baihetan Hydropower Station in the Jinsha River[D]. Chengdu: Chengdu University of Technology: 72–75 （in Chinese）.
徐湘涛. 2012. 金沙江白鹤滩水电站高边坡岩体力学特性及其稳定性研究[D]. 成都: 成都理工大学: 35–36. Xu X T. 2012. Research on Rock Mass Mechanical Characteristics and Stability of High Slope of Baihetan Hydropower Station[D]. Chengdu: Chengdu University of Technology: 35–36 （in Chinese）.
殷跃平,王猛,李滨,冯振. 2012. 汶川地震大光包滑坡动力响应特征研究[J]. 岩石力学与工程学报,31(10):1969–1982. doi: 10.3969/j.issn.1000-6915.2012.10.003 Yin Y P,Wang M,Li B,Feng Z. 2012. Dynamic response characteristics of Daguangbao landslide triggered by Wenchuan earthquake[J]. Chinese Journal of Rock Mechanics and Engineering,31(10):1969–1982 (in Chinese).
云南省巧家县志编纂委员会. 1997. 巧家县志[M]. 昆明: 云南人民出版社: 61–63. Chorography Compilation Committee of Qiaojia. 1997. History of Qiaojia[M]. Kunming: People’s Press of Yunnan: 61–63 （in Chinese）.
张明,王正波,孙琳. 2016. 滑坡碎屑流高速远程机制环剪试验研究[J]. 岩石力学与工程学报,35(增刊1):2673–2681. doi: 10.13722/j.cnki.jrme.2014.1609 Zhang M,Wang Z B,Sun L. 2016. Research on rapid and long-runout mechanisms of rockslide debris using ring shear tests[J]. Chinese Journal of Rock Mechanics and Engineering,35(S1):2673–2681 (in Chinese).
张曙光. 2007. 金沙江白鹤滩水电站高拱坝建设工程地质适宜性研究[D]. 成都: 成都理工大学: 31–32. Zhang S G. 2007. Engineering Geological Feasibility Study of High Arch Dam Construction at Baihetan Hydropower Station, Jinsha River[D]. Chengdu: Chengdu University of Technology: 31–32 （in Chinese）.
张欣,王运生. 2017. 白鹤滩水电站库区小江断裂带活动性研究[J]. 工程地质学报,25(2):531–540. doi: 10.13544/j.cnki.jeg.2017.02.033 Zhang X,Wang Y S. 2017. Activities of Xiaojiang fault zone in Baihetan hydropower station reservoir[J]. Journal of Engineering Geology,25(2):531–540 (in Chinese).
周洪福, 方甜, 韦玉婷. 2023. 国内外地震滑坡研究: 现状、问题与展望[J]. 沉积与特提斯地质, 43（3）: 615−628. Zhou H F, Fang T, Wei Y T. 2023. Research situations and suggestions on earthquake-induced landslides[J]. Sedimentary Geology and Tethyan Geology, 43（3）: 615−628 （in Chinese）.
Cho S H,Kaneko K. 2004. Influence of the applied pressure waveform on the dynamic fracture processes in rock[J]. Int J Rock Mech Min Sci,41(5):771–784. doi: 10.1016/j.ijrmms.2004.02.006
Cui S H,Yang Q W,Pei X J,Huang R Q,Guo B,Zhang W F. 2020. Geological and morphological study of the Daguangbao landslide triggered by the MS8.0 Wenchuan earthquake,China[J]. Geomorphology,370:107394. doi: 10.1016/j.geomorph.2020.107394
Fan X M,Scaringi G,Korup O,West A J,van Westen C J,Tanyas H,Hovius N,Hales T C,Jibson R W,Allstadt K E,Zhang L M,Evans S G,Xu C,Li G,Pei X J,Xu Q,Huang R Q. 2019. Earthquake-induced chains of geologic hazards:Patterns,mechanisms,and impacts[J]. Rev Geophys,57(2):421–503. doi: 10.1029/2018RG000626
Hungr O,Leroueil S,Picarelli L. 2014. The Varnes classification of landslide types,an update[J]. Landslides,11(2):167–194. doi: 10.1007/s10346-013-0436-y
Keefer D K. 1984. Landslides caused by earthquakes[J]. GSA Bull,95(4):406–421. doi: 10.1130/0016-7606(1984)95<406:LCBE>2.0.CO;2
Liu W,He S M. 2016. A two-layer model for simulating landslide dam over mobile river beds[J]. Landslides,13(3):565–576. doi: 10.1007/s10346-015-0585-2
Richter C F. 1958. Elementary Seismology[M]. San Francisco: W. H. Freeman and Co: 21–23.
Valagussa A,Marc O,Frattini P,Crosta G B. 2019. Seismic and geological controls on earthquake-induced landslide size[J]. Earth Planet Sci Lett,506:268–281. doi: 10.1016/j.jpgl.2018.11.005
-
期刊类型引用(3)
1. 董建明，张志萍，黄纲领，刘静伟. 黄土-泥岩接触面滑坡的启动及运动特征. 铁道勘察. 2025(02): 65-71 . 
百度学术
2. 张宇，徐宗恒，查玲珑，陈云英，陶真鹏. 不同计算方法的云南省永胜县下院滑坡堰塞湖沉积物粒度特征及沉积历史重建. 中国地质. 2024(05): 1761-1775 . 百度学术
3. 曾京，黄飞鹏，吴中海. 基于构造地貌指数的金沙江干流涛源古湖的成因探讨. 地震科学进展. 2024(10): 742-756 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 134
- HTML全文浏览量: 44
- PDF下载量: 39
- 被引次数: 4
