铁永波,张宪政,卢佳燕,2,梁京涛,王东辉,马志刚,李宗亮,鲁 拓,2,石胜伟,刘民生,巴仁基,何龙江,张新克,甘 伟,陈 凯,高延超,白永健,龚凌枫,曾孝文,徐 伟
(1.中国地质调查局成都地质调查中心,四川 成都 610081;2.中国地质大学(北京)工程技术学院,北京 100083;3.四川省地质调查院,四川 成都 610084;4.四川省国土空间生态修复与地质灾害防治研究院,四川 成都 610084;5.中国地质科学院探矿工艺研究所,四川 成都 610031;6.应用地质研究中心,四川 成都 610036)
据中国地震台网中心发布的数据,2022年9月5日12 时52 分,四川省甘孜州泸定县磨西镇发生Ms6.8 级地震,震中位于29.25°N、102.08°E,震源深度16 km,最高烈度达到Ⅸ度。地震造成大量交通、通信、电力及建筑等基础设施受损,同时诱发大量崩塌、滑坡等地震地质灾害。截至2022年9月14日17时,地震已经造成93 人遇难,25 人失联。
地震地质灾害是指地震发生时因地面强烈震动而引起的崩塌、滑坡等地质灾害,也可称为同震地质灾害,其破坏模式、表现形式及致灾特征一直以来都受到国内外学者的广泛关注[1-5],在地震地质灾害的早期识别、快速解译、综合评估及区划等研究方法和模型方面取得了较多进展[6-10]。已有研究表明,地震地质灾害的数量及规模与距断层距离、距震中距离、地震烈度、地形及岩土体结构等密切相关[11-18],并根据不同地区、不同震级地震开展了大量的规律统计分析与成因机制研究[19-23],取得了较多的研究成果,为地震地质灾害的机理认识提供了重要基础[24-28]。许多学者认为,对强震引发的地质灾害的防范不能只关注短期防治,需要从更长远、更系统的视角进行设计[29-32]。事实证明,一次强震后,震区内地质灾害的活跃将会持续数十年之久,且震后次生地质灾害的群发性、链式效应等特征会更为突出,且规模会明显增大,这种现象在2008年汶川Ms8.0 级地震发生后近15 a 内仍极为典型[33-35]。因此,深入分析地震地质灾害规律及特征,对震区灾后重建与次生地质灾害防治具有重要的指导意义。此次研究聚焦泸定县Ms6.8 级地震诱发的地质灾害发育分布特征,在分析控灾因素和震后地质灾害演化趋势的基础上,提出了震区地质灾害防灾减灾建议,为震区灾害重建规划与地质灾害防治提供科学依据。
泸定县位于四川省甘孜藏族自治州东南部。此次地震震中位于泸定县南部的磨西镇,区内地形高差大(磨西镇境内最高海拔可达7 556 m,与最低处海拔高差超过5 000 m),地震波随山体高度增加而放大的效应极为突出,为地震诱发高位崩塌、滑坡提供了地形条件。同时,区域内水系发育,有雅砻江、大渡河等大江大河,河谷深切、支流发育,加之区内构造应力强[36-37],岩体被河流切割后沿临空方向卸荷并导致裂隙发育,为同震崩塌、滑坡的形成提供了地质条件。此次震中所在的磨西镇属于典型亚热带季风气候,主要受东南季风和来自四川盆地的暖湿气流影响,区内地形高差大导致降水量随海拔差异较大:磨西镇海拔1 600 m 处多年平均降水量为897.8 mm,海拔3 000 m处多年平均降水量为1 941.5 mm,暴雨频发,为震后次生地质灾害的发生提供充足的水动力条件。
研究区内地层岩性复杂,从区域1∶130 万比例尺地质图(图1)上可以看出,从第四系至震旦系地层均有分布,广泛出露震旦系(Z)、志留系(S)、泥盆系(D)、石炭系(C)、二叠系(P)、三叠系(T)、古近系(E)和第四系(Q)地层, 其中,以中生界、古生界三叠系及元古界变质岩分布最广。研究区西部区域主要出露三叠系的砂岩、板岩和少量石灰岩。研究区东部区域主要出露石英岩、大理岩、白云石、变质砂岩、碳酸盐岩、基性和酸性火山岩。侵入岩体从元古代到新生代均有出露,岩性以花岗岩类、闪长岩类和少量超基性侵入岩为主,属坚硬岩组。花岗岩类主要分布在鲜水河断裂带以西,闪长岩类在研究区北部零星分布。第四系沉积物主要分布在鲜水河及其支流的相邻区域和该区的一些小盆地,包括全新世冲积层、砂、黏土、砾石,更新世冲积层、砾石土、黏土。研究区内多时代、多类型的岩性背景奠定了地质灾害孕灾的复杂性,特别是河谷斜坡上部的破碎岩体在地震下为群发性浅表层崩塌的形成提供了有利条件,而坚硬花岗岩体在高陡临空条件下形成卸荷裂隙,为高位大块体岩质崩塌的形成提供了有利条件。
图1 研究区地质概况图Fig.1 Geological overview of the study area
磨西镇所处地质构造部位属于鲜水河断裂带南端,鲜水河断裂是我国西部“Y”字形断裂带的重要组成部分,为一条大型左旋走滑断裂带,经炉霍、道孚、乾宁(八美)、康定延伸至石棉安顺场一带,全长约350 km,走向整体呈NW—SE。鲜水河断裂带活动强,历史上曾发生过有记载的震级Ms7.0 级以上地震8 次,如1786年6月1日四川康定—泸定Ms7.75 级地震,与此次泸定Ms6.8 级地震空间上较为接近。泸定县Ms6.8 级地震发震断裂为鲜水河断裂带南东段的磨西断裂,震中位于磨西镇海螺沟风景区内,地震最高烈度达到Ⅸ度,影响区范围主要涉及3 个市、州共12 个县(市)。由于地震对部分县(市)的影响区面积较小,且未引发明显新增地质灾害隐患点,故此次研究区只选择其中10 个县(市),分别是甘孜州泸定县、康定市、九龙县、雅江县,雅安市石棉县、天全县、荥经县、汉源县,凉山州甘洛县、冕宁县,震区影响范围约4.36 ×104km2。
根据四川省地质灾害隐患点在库统计数据,此次震区10 县(市)境内有地质灾害隐患点共1 797 处(截至2022年9月4日),灾害类型主要为滑坡、泥石流和崩塌,其中滑坡887 处、泥石流677 处、崩塌233处(表1)。从空间分布特征上看,震前研究区地质灾害主要沿河流呈线状分布。根据应急管理部和中国地震局提供的烈度分布数据,泸定县Ms6.8 级地震烈度在Ⅵ度以上影响范围内,已有地质灾害点722 处,沿大渡河两岸呈线状分布(图2),主要分布在鱼同镇—泸定县城—得妥镇,安顺场镇—石棉县城—汉源县城。
表1 震区震前地质灾害点统计表Table 1 Statistics of geological hazard points before the earthquake in the study area
泸定县Ms6.8 级地震影响范围内,震前地质灾害发育,且其分布密度与烈度等级呈正相关(图3)。烈度在Ⅷ度以上影响范围内,已有地质灾害点81 处,主要分布在燕子沟镇、磨西镇和德妥镇,尤以磨西镇最多,少数分布在草科乡和王岗坪乡。在Ⅶ度影响范围内,已有地质灾害点292 处,主要分布在泸定县城—德威镇,少数分布在安顺场镇附近。在Ⅵ度影响范围内,已有地质灾害点共349 处,主要分布在泸定县—磨西镇段和石棉县—汉源县段:泸定县—磨西镇段地质灾害点共计361 处,沿大渡河两岸密集分布;石棉县—汉源县段地质灾害点共计188 处,沿河流两岸呈散点状分布。
图3 震区震前地质灾害点与地震烈度统计图Fig.3 Geological hazard points and seismic intensity before the earthquake in the study area
为突出泸定县Ms6.8 级地震地质灾害的发育特征,本研究聚焦此次地震烈度Ⅷ度、Ⅸ度区主要分布的泸定县和石棉县。为突出“以人为本”的原则,重点对有威胁对象的地质灾害隐患点开展地面排查和遥感解译,主要包括居民区、重大工程区、公路等区域。地面排查工作主要通过地面现场开展,遥感解译数据来源于高分六号、高分二号、北京三号卫星遥感影像,时间为2022年9月9—11日。
截至2022年9月14日22 时,泸定县和石棉县境内共排查有威胁对象的地震地质灾害隐患点273 处,此次遥感解译有威胁对象的地质灾害隐患点共292处。地震还导致泸定县境内53 处已有地质灾害隐患点加剧变形、石棉县境内28 处已有地质灾害隐患点加剧变形。
研究表明,此次地震地质灾害类型主要为中、小型规模高位崩塌滑坡,且集中分布在震中附近地震烈度Ⅸ度区域,包括泸定县磨西镇、得妥乡、得妥乡—德威镇段大渡河两岸及石棉县草科乡、王岗坪乡;地质灾害较为发育的乡镇有泸定县磨西镇、得妥镇和石棉县草科乡、王岗坪乡,其余地区地震地质灾害发育程度相对较低(图4)。磨西镇地震地质灾害主要集中分布在2 个区域:一个区域是磨西台地临燕子沟和雅家埂河一侧临空区域,因磨西台地最大高度可达110 m,为冰水堆积和泥石流堆积体,虽然台地胶结较好,但在强震作用下仍产生大量崩塌,这些区域的崩塌也是导致磨西镇至得妥乡道路中断的主要原因,如图5(a)所示;另一个区域是磨西镇两岸高位山体,主要发育高位崩塌,多分布在海螺沟两岸、燕子沟磨西段右岸、雅家埂河磨西段左岸坡体中上部,能划分出单个崩塌的边界,分布相对较为分散。
图4 研究区地震新增地质灾害分布图Fig.4 Distribution of new seismic geological disasters in the study area
泸定县得妥镇和石棉县草科乡、王岗坪乡境内地震烈度Ⅸ度区域,地质灾害密集发育在坡体和山脊中上部临空侧,多为浅表层崩塌滑坡,且密度极大,多呈片状分布,难以划出崩塌之间的边界,表现为崩塌群形式,如图5(b)所示。
图5 震区典型地质灾害发育特征Fig.5 Characteristics of typical geological hazards in seismic area
得妥乡—德威镇段大渡河两岸地质灾害多发育在坡体上部,崩塌运动表现出2 种差异显著的模式:相对较为坚硬的块石多以落石方式崩落,且运动距离较远,部分直接进入大渡河,多数停留在坡脚下公路、民房处及低洼地带并成灾,这也是导致得妥镇至德威镇段道路遭受影响的主要原因,如图5(c)所示;相对碎裂的岩体崩落后,在运动过程中解体,以碎块石的方式堆积在斜坡中下部地形低缓处,其运动距离相对较短,造成直接破坏的范围有限,如图5(d)所示。
(1)地形地貌对地震地质灾害的控灾机制
此次地震地质灾害集中分布在泸定县磨西镇、得妥镇和石棉县草科乡、王岗坪乡,该区域属于流水侵蚀形成的沟谷地貌,河谷深切而狭窄,为地震地质灾害的诱发提供了极大的势能条件。同时,地震地质灾害多发育在大渡河深切河谷两岸及一级支流流域内,河流侵蚀基准面相对较低,下切作用强烈,多以V 型谷为主,高陡岸坡和狭窄沟道形成的陡坡地形为地震地质灾害的发生提供了有利条件。
以大岗山水电站至王岗坪乡区域为例,该区域南边长约3.2 km,宽约2.8 km,面积约9 km2,基于如图6(a)所示的9月9日震后航空影像解译可知,该区域共发育地质灾害点119 处,其中滑坡91 处,崩塌28 处,平均发育密度为13.2 处/km2。根据灾害规模统计结果可知,大型规模地质灾害发育2 处,中型发育36 处,小型发育81 处,灾害规模以中小型为主。由图6(b)所示的该地区震前光学影像可知,未见明显地质灾害发育,说明该区域地震前后灾害发育密度和数量发生明显变化。
图6 大岗山水电站至王岗坪乡地震前后地质灾害发育情况对比图Fig.6 Comparison of the distributions of geological hazards in the area from Dagangshan hydropower station to Wanggangping Town before and after the earthquake
(2)地震烈度对地震地质灾害的控灾机制
地震烈度是对地表变形影响最为显著的影响因素。本文基于2022年8月14日与2022年9月7日获取的2 期Sentinel-1 升轨数据,使用DInSAR 技术获得了泸定地震同震InSAR 形变场(图7)。其中红色(负值)表示该区域下沉变形,蓝色(正值)表示抬升变形。通过将InSAR 结果和地震烈度叠加分析可知,磨西镇和草科乡是研究区内地表形变最强烈的区域,这与地震烈度的分布是相吻合的。因此这2 个区域内的地震灾害发育概率相对较大。
图7 研究区地震InSAR 形变场(四川省地质调查院提供)Fig.7 Seismic InSAR deformation field in the study area(provided by Sichuan Geological Survey)
通常情况下,地质灾害发育数量、规模及破坏力均与地震烈度有显著相关性。通过研究泸定县Ms6.8 级地震地质灾害分布数量与地震烈度的对应关系,发现此次新增的地震地质灾害主要分布在地震烈度Ⅸ度区和Ⅷ度区,其数量约占总数的68%(表2)。虽然此次地震地质灾害规模多以中、小型为主,在数量上却表现出与地震烈度等级的极好相关性。同时,地震地质灾害空间分布特征表明,随着地震烈度的衰减,地震地质灾害发育数量呈急剧减少的特点,这主要与地震能量释放随着距震中距离远而明显衰减的控灾机制有关。
表2 震区震前和震后地质灾害点与地震烈度统计表Table 2 Statistics of geological hazard points and seismic intensity before and after the earthquake in the earthquake area
(3)特殊结构土体对地质灾害的控灾机制
泸定县Ms6.8 级地震影响区分布有海拔7 556 m的贡嘎山,因该地区强烈的构造隆升和特殊气候环境,使得该地区海洋性冰川极为发育。早期冰川形成的大量冰碛土和冰水堆积体在磨西河流域内广泛分布,且规模巨大,与其他地区第四系土体成因及分布有明显不同。如磨西镇所在的磨西台地就是早期贡嘎山冰川消融和后期河流下切后形成的冰水堆积阶地,其纵向长约11 km,垂直高差最高可达110 m。从此次地震地质灾害发育特征看,磨西镇境内冰碛土边坡发生破坏的区域主要分布在海螺沟景区3 号营地至4 号营地之间的公路边坡、岸坡,以垮塌为主,主要受坡面临空和冻融作用影响,在地震时产生破坏,成为泥石流的潜在物源。磨西台地两侧临空边坡受地震影响后发生多起垮塌,虽然磨西台地的岩土体具有较好的胶结性,但因临空面坡度较陡,岩土体在长期重力卸荷作用下形成纵向拉裂缝,在强震作用下产生垮塌。除此外,其他类型成因土体区土质崩塌的发育数量相对较少,表明磨西镇境内特殊岩土体对地震地质灾害的明显控灾特征。
(1)加强区域性地震地质灾害成因机制与规律研究
目前已有大量关于地震地质灾害规律及成因机制的研究成果,但不同区域内的地形、地质条件和地震属性存在差异,导致不同区域内的地震地质灾害的发育规律存在明显差异。此次泸定县Ms6.8 级地震诱发的地质灾害发育特征是地形地貌、地质背景条件及地震属性等共同作用的结果,加强研究区不同地震属性、不同地质背景下的地震地质灾害成因机制与规律研究,对未来该地区地质灾害演化趋势及防治技术防范具有重要理论与实践意义,特别是要加强西南高山峡谷区震后交通中断条件下快速遥感识别与调查等新技术在地震地质灾害规律认识中的应用实践。
(2)加强地震地质灾害多致灾模式的演变趋势研判
地震发生后形成的大量不稳定斜坡、地裂缝及松散堆积体,他们在余震、降雨等作用下可能发生震后次生灾害,如崩滑体转化成坡面泥石流,震裂岩体转化成崩塌或滑坡等。因此需要根据此次地震地质灾害的发育特征,结合其所处地貌部位及结构特征等进行分析,研究其在不同外动力作用下可能的破坏模式及成灾表现形式,可通过物理模型试验、数值模拟等定量研究手段,结合定性分析进行综合研判,科学划定地质灾害可能的影响区范围与致灾模式,并在此基础上提出针对性的防范措施建议。
(3)加强震区泥石流防治工程设计参数确定的科学性与合理性
震区泥石流防治工程设计要充分汲取汶川震区泥石流破坏特征,通过类比优化震后泥石流动力学参数计算公式,结合未来地质灾害演变趋势提出针对性防治措施。特别是对坡面有大量松散堆积体的崩塌滑坡区,在防治工程设计中除了考虑对崩塌滑坡本身的防治措施外,还要综合考虑松散土体在降雨作用下转化成泥石流的防治方案。此外,对于已有泥石流防治工程的流域,需要进一步根据震后流域内松散物源增加量评估现有工程的防治成效,评估地震对已建工程的影响,并针对性地提出下一步防治措施建议。
(4)加强地震地质灾害密集发育区微小流域综合治理规划
一次强震后地质灾害的表现形式及链式效应极为突出,对这类地质灾害的治理需要提出系统的防治方案,要结合流域上下游、灾害规模及保护对象进行综合考虑。如此次地震地质灾害密集发育的泸定县得妥乡湾东河流域内,在干流和支流中均发育有大量崩滑体,在强降雨作用下,支流发生的泥石流进入干流会形成一定的堵溃。若多条支流同时暴发泥石流,其堵溃效应会更明显。虽然该流域内沟道纵坡比降大,泥石流堵溃后的规模放大有限,仍会形成超常态的流量,其危险区可能会扩大几倍,因此,需要采用流域水土保持、重点河道清淤、泥石流物源稳固等多管齐下的方式,联合自然资源、水利、林草及水电站业主单位共同推进相关工作,从源头防范震后次生灾害风险。
(1)泸定Ms6.8 级地震地质灾害主要以浅表层中、小规模崩塌和滑坡为主,主要集中分布在坡体中上部,在震中附近具有典型的群发性特征,与地震Ⅸ度烈度区具有较好的对应关系。
(2)陡峻的地形地貌、高地震烈度及特殊岩土体是影响泸定县Ms6.8 级地震地质灾害发育数量、空间分布、规模及破坏模式的主要影响因素,今后在地震活跃区重大工程及村镇建设选址中需加以注意。
(3)在地震地质灾害集中发育的泸定县磨西镇、得妥镇,石棉县王岗坪乡、草科乡境内,地震地质灾害形成的大量松散物将成为泥石流的潜在物源,支沟群发性高频黏性泥石流将是该区域未来地质灾害的主要表现形式,需要加强该地区泥石流的趋势研判与风险防控。
(4)此次初步排查和遥感解译的新增地质灾害隐患点主要针对有威胁对象的点,用于分析此次地震地质灾害的宏观规律与趋势。目前震中区所在的泸定县得妥镇、石棉县草科乡部分区域仍尚未完成地面地质灾害排查,随着后续各项工作的推进,震区新增地质灾害隐患点数量会有所增加,特别是无人类工程活动区群发性崩塌的数量会有明显增加,因此关于地震地质灾害新增点的分布规律还需进一步研究。