同震滑坡是最具破坏的地震次生灾害之一,可以造成大范围严重的破坏,尤其是山区,了解同震滑坡的易发性控制因素是预测潜在地震区同震滑坡易发性、危险性和风险性的基础,对防震抗灾具有重要意义。然而,控制同震滑坡发生的因素众多,包括气象水文、地形地貌、地质构造、地层岩性、岩体破碎程度、同震断层、地震、场地烈度、人类活动、土地利用等,各个因素之间相互影响,相互依存。为此,前人发展了信息量法、主成分分析法、频率比法、多元统计方法、逻辑回归方法、蒙托卡罗方法、模糊数学方法、神经网络方法等来开展同震滑坡易发性、危险性及风险性的评价。然而由于主控因素不明,现有的研究往往要采集各因素的海量数据,在讨论各因素与同震滑坡空间分布的统计关系上,进而开展潜在震区同震滑坡易发性、危险性及风险性评价,难以高效准确实现同震滑坡空间分布预测。
针对上述重要问题,中国科学院地质与地球物理研究所邹宇工程师、祁生文研究员、郭松峰副研究员等,以2014年8月3日中国云南省鲁甸县发生的Mw6.1级地震为例开展研究。该地震由前期未测绘的隐伏断层(包谷垴-小河断层)引发,同震滑坡分布并未沿发震断层两侧展布,如汶川地震等揭示的情况。基于建立的同震滑坡及影响因子数据库,通过空间分析揭示了同震滑坡空间分布的主控因素。通过对震前、震后山体滑坡数据进行分析,结合图像解译和实地调查,在面积为704.7 km2的研究区内,共识别出1414个滑坡,建立了鲁甸地震滑坡空间数据库(图1)。详细绘制了研究区构造、岩性等地质背景图。利用地理信息系统(GIS)分析滑坡的空间特征,通过GIS空间分析,表征了滑坡发生与坡度、高程、坡向、地震烈度、距发震断层距离(以下简称断层距)、距震中距离、岩性、距非发震断层距离、距河流距离、距公路距离等因素之间的相关性。
图1 鲁甸地震位置图及同震滑坡数据库
利用空间分析方法,分析了各个因素与同震滑坡易发性之间的关系,在此基础上,利用曲线下面积方法,对各个因素之间的相关性展开研究。通过统计分析发现,坡度和断层距是同震滑坡易发性的两个控制性因素,随着断层距的增加,同震滑坡发育密度呈明显的指数下降趋势,但随着坡度的增大,同震滑坡密度表现出良好的韦伯累积分布(图2,图3)。其他因素的影响受到这两个因素的控制,其对同震滑坡易发性的贡献通过坡度和断层距两个因素得以体现。
图2 滑坡密度与断层距的关系
图3 滑坡密度与坡度的关系
以高程为例,分析坡度、断层距与高程之间的关系。在700-1200 m范围内的高程区域,边坡坡度在10°以下的百分比约9%,在20°以下的百分比约28%;而在>2700 m范围内的高程区域,边坡坡度在10°以下的百分比约为30%。在20°以下的百分比约74%。可见,随着高程的增加,边坡坡度逐渐降低(图4)。
与此同时,在700-1200 m范围内的高程区域,断层距在3 km内的百分比约为29%,在6km内的百分比约为52%;而在>2700 m的高程区域,断层距在3 km内的百分比约为0,在6 km内的百分比约为18%。可见,随着高程的增加,断层距逐渐增加(图5)。图4、图5阐述了坡度、断层距与高程的关系,清晰阐明随着高程的增加,滑坡密度逐渐下降的现象是由坡度、断层距控制的结果(图6)。
图4 高程与坡度之间的关系
图5 高程与断层距之间的关系
图6 滑坡密度与高程的关系
该研究表明,对于即使隐伏的同震断层,其同震滑坡易发性和非隐伏同震断层易发性类似,比如汶川地震(Qi et al, 2010),坡度和断层距两因素是其控制性因素,这一研究对于同震滑坡易发性的预测评价具有重要意义。
研究成果发表于国际学术期刊Engineering Geology (邹宇,祁生文*,郭松峰,郑博文,詹志发,何乃武,黄晓林,侯晓坤,刘海洋. Factors controlling the spatial distribution of coseismic landslides triggered by the Mw 6.1 Ludian earthquake in China [J]. Engineering Geology, 2022, 296: 106477. DOI: 10.1016/j.enggeo.2021.106477)。研究得到第二次青藏高原综合科学考察研究等项目资助(QZKK20190904)。