金沙江白格滑坡失稳机理及影响因素分析_杨连伟.pdf

第 卷 第 期 年 月人 民 长 江 ，收稿日期：作者简介：杨连伟，男，高级工程师，主要从事水文地质、工程地质、环境地质等方面的研究。：文章编号：（）引用本文：杨连伟 金沙江白格滑坡失稳机理及影响因素分析 人民长江，（）：，金沙江白格滑坡失稳机理及影响因素分析杨 连 伟（江西省地质环境监测总站，江西 南昌）摘要：峡谷区高位滑坡发生后堵塞河道，往往造成巨大的损失。年 月 日及 月 日，金沙江白格滑坡两次堵塞金沙江，形成的堰塞湖淹没了上游村镇，堰塞坝溃决后洪水冲击下游造成了巨大损失。为研究白格滑坡首次失稳破坏的机理及其关键影响因素，依据现场调查建立了失稳前边坡的二维模型，并结合地质过程进行了数值模拟计算。结果表明：白格滑坡是深切河谷斜坡岩体在长期卸荷和表生时效作用下产生的岩质滑坡；斜坡破坏时主应力在软硬岩接触面附近发生应力集中，致使下部硬岩剪出口位移超过 最终失稳破坏；对滑坡区工程地质条件、岩体卸荷及表生改造作用等因素的分析认为，“”白格滑坡是随时间推移、多种因素叠加，量变转化为质变形成的。关 键 词：滑坡堰塞坝；失稳机理；白格滑坡；金沙江中图法分类号：文献标志码：引 言滑坡堵江是指河谷斜坡岩土体因崩塌、滑坡及其转化成的碎屑流等造成江河堵塞和回水的现象。河谷发育过程中，随着河流下切，岸坡岩体卸荷回弹，形成一系列的顺坡向裂缝，其后裂缝间的岩桥在岩体重力或降雨、地震等荷载的作用下被突破，造成滑坡或崩塌堵塞主河道，形成堰塞湖；随着堰塞湖水的增高，堰塞体会溃坝并形成洪水事件。这一过程，会对上下游造成巨大的破坏。目前对滑坡堵江的研究主要集中在历史滑坡堵江的识别统计、堰塞体的特征、溃坝最大洪水量估算、滑坡堵江的危害及治理等方面。年 月 日（另据地震台网监测为 日：）和 月 日，西藏自治区江达县波罗乡白格村与四川省白玉县绒盖乡则巴村交界处两次发生大规模高位滑坡灾害（以下分别简称“”和“”白格滑坡），岩土体阻断金沙江干流，形成堰塞坝。这两次滑坡 堵江事件发生后，先形成堰塞湖淹没上游（江达县波罗乡、白玉县金沙乡等），堰塞体溃决之后，形成洪水冲击下游，尤其是第二次堵江坝体溃决后，使下游多座桥梁被冲毁，多个城镇被淹，在建的苏洼龙水电站也受到冲击，造成了巨大的经济损失和广泛的社会影响。白格滑坡 堵江事件具有严重的灾害链效应，事后诸多学者对其进行了详细研究，主要包括白格堰塞湖应急抢险处置 及水文监测预报 的经验总结、滑坡残留体后期稳定性及堵江风险分析、堰塞湖溃坝洪水数值分析 以及白格滑坡的形成演化过程分析。尽管前人从不同角度对白格滑坡进行了调查研究，但由于早期交通条件限制及研究目的不同，对于其失稳机理及关键影响因素仍未形成统一的认知。因此，本文对滑坡区的地质条件进行了详细的现场调查和分析，在分析相关资料及数值的基础上，对白格滑坡首次滑动的形成机理、演化过程和关键影响因素进行了认真研究，并计算了滑坡的运动速度。人 民 长 江 年研究结果可为青藏高原周边高山峡谷地区斜坡灾害的研判和预警工作提供借鉴。金沙江白格滑坡概况 白格滑坡基本特征白格滑坡区处于青藏高原东部，金沙江缝合带上，周边沟谷纵横，山谷相间。两次滑坡的滑源区均位于金沙江“”形河谷右岸的猪背形山脊上（见图），山脊走向约 ，斜坡朝向 ，海拔高程 ，垂直高差达 。图 主滑坡区出露地层 经调查，滑坡区斜坡中上部主要出露华力西期蛇纹岩（）（见图），高程在 之间；受构造影响，蛇纹岩整体成碎粉岩状，呈灰绿 绿白色。滑源区下游侧及 高程以下主要为元古界雄松群（）英长片麻岩，层面内倾，产状，发育有 、两组结构面，片麻岩表面风化后呈土黄色 青灰色。年 月 日第 次发生滑坡时，滑坡后缘高程约 ，滑源区剪出口高程约 。主滑区共有 万 万 岩土体失稳，堵塞金沙江后形成堰塞坝。堰塞湖水于 月 日开始自然漫坝下泄，至 日，泄洪基本完成。月 日在第 次滑坡的后壁上发生第 次滑坡，约 万 岩体失稳形成碎屑流，与第 次堰塞残留体一起形成新的堰塞坝，第 次堵塞金沙江。经人工干预，至 月 日下午堰塞坝被完全冲开。“”滑坡发生后，在滑坡源区（）形成高陡岩质后壁，而主滑区后部为后部变形体（），下游侧变形体（）和上游侧变形体（）分布于主滑区两侧（见图（）。在下部堆积体中，近滑源区侧（侧）地势相对低洼，主要以蛇纹岩碎屑体为主；远滑源区侧（侧）高程相对较高，主要以片麻岩堆积体为主。“”和“”两次滑坡后的泄洪通道均位于地势相对低洼的蛇纹岩堆积体中。因此，研究“”图 两次白格滑坡后影像 白格滑坡的形成机理，需要深入研究滑坡堰塞体的具体特征。“”白格滑坡 堰塞体特征分析堰塞体主要位于主滑坡区下部的流通堆积区（）中。依据坡面特征，将区划分为图 所示 个亚区；区外侧为涌浪影响区（）。图 “”白格滑坡堆积 堰塞体特征（摄于 年 月 日）“”，（）流通铲刮亚区（）和表层铲刮亚区（和）。岩土体失稳后铲刮下部斜坡覆盖层。区位于滑坡源区正下方，受到铲刮最严重，斜坡表面后退约 ，形成一个槽状区，该区域宽约，高约 。两侧的 区和 区位于剪出口侧边，受失稳岩土体的冲击力和铲刮力较小，仅被铲刮掉表面覆盖层，裸露出两处片麻岩基岩陡壁。第 期 杨连伟：金沙江白格滑坡失稳机理及影响因素分析（）堰塞体（主堆积体，）。整个堰塞体平面呈纺锤形（见图），沿河谷方向长约 ，中部垂直河谷方向宽约 ，水平投影面积约 万。图 中三角形所示为整个堰塞体的最高点。堰塞体近滑源区侧（西侧）地势相对低洼（见图 ），坡体呈灰绿 绿白色，以蛇纹岩碎屑为主，夹有粒径不同的墨绿色蛇纹岩晶体。该区域上游侧宽 ，下游侧宽 。远滑源区侧（东侧）地势相对较高（见图 ），坡体整体呈灰黄色夹青灰色，主要以片麻岩碎屑为主，其间分布大量直径约 的片麻岩块石。该区域上游侧宽 ，下游侧宽 。沿河谷方向，区下部的堆积体地势最高，整体呈条垄状，走向 ，垄脊高出河床约 。条垄南侧为缓斜坡，发育 级与垄走向近平行的下错台坎（见图），台坎下部为张性裂缝。台坎 下部以厚度小于 的松散粉砂性土为主，现场调查期间本区饱水，人踩后易湿陷。垄的北侧为陡坡，堆积体向下被堰塞湖水体淹没（见图）。图 主堆积区南侧的横向台坎（镜向）（）图 主堆积体北侧的陡坡（镜向）（）两次泄流的泄洪槽均位于堰塞体西侧相对低洼的蛇纹岩堆积区中（见图）。蛇纹岩与片麻岩分布区的界线特征如图 所示。图 堰塞体上蛇纹岩与片麻岩分布区界线 （）涌浪影响区（）。金沙江水受失稳岩体冲击后形成涌浪冲向对岸，被对岸斜坡阻挡后再落回江底，期间冲刷坡表的覆盖层。经调查，涌浪影响区斜坡上的灌木和植草大面积倒伏，倒伏方向指向江底。经测量，影响区边界高程最高可达 ，高出金沙江河床约 。饱水土体被涌浪裹挟冲上对岸斜坡后又沿坡表滑下，留下深色残留体。涌浪回落反向冲刷已经停积的堆积体，造成堆积体与涌浪影响区接触部位形成一条冲沟。涌浪反向冲刷并在堆积体表面形成与主滑向相反的泥浆冲刷痕迹（见图）。图 泥浆冲刷痕迹（镜向）（）白格滑坡首次滑动的失稳机理及动力特性 “”白格滑坡的失稳机理分析（）堆积区中的岩石都呈碎块 碎粉状，原岩结构不可辨，碎块石平均粒径为 ，局部出露直径超过 的块石。据地质资料，研究区受深层韧性剪切作用，使蛇纹岩碎裂成碎粉岩状，片麻岩体节理裂隙发育。岩体失稳后经过近 的运动和内部相互碰撞，会变得更加破碎，使原岩的结构在堆积体中难以分辨；但堆积区的岩性层序与滑源区层序一致，上部（西侧）为蛇纹岩，下部（东侧）为片麻岩，故白格滑坡为岩质滑坡。人 民 长 江 年（）斜坡失稳区呈纺锤形（见图），其上部蛇纹岩中滑坡后壁的平均坡度约为，后缘局部陡壁坡度可达 左右；而下部剪出口内侧片麻岩中滑面的平均坡度约为。滑坡失稳后，下滑合力 向前（见图），主应力作用面近于垂直，破裂面与主应力作用面的夹角 （为岩体内摩擦角），则破裂面倾角 ，值越大则 值越小，反之则 越大。由于软质蛇纹岩中的滑面倾角较大，而硬质片麻岩中的滑面倾角较小，从而在滑源区中部形成一个巨大的槽状区。注：和 分别为滑坡物质在运动方向上的水平推力和法向自重力，为合力。图 滑体运动铲刮作用细观机制 （）调查发现：堆积体东侧有涌浪回落时的反向冲刷痕迹；堰塞湖水体在堰塞体西侧的低洼区过流。此处蛇纹岩碎屑体质软，易受冲蚀，很快形成泄洪槽，槽壁因受洪水侧蚀、塌岸而逐渐后退，泄洪槽两侧坡体上无冲刷痕迹；坡体上未见到二次涌浪拍打迹象。因此，“”白格滑坡为一次性失稳破坏并冲击江面。（）综上所述，斜坡上部破碎的蛇纹岩在岸坡卸荷及表生时效作用下，内部的裂隙渐进性地发展为贯通性裂缝（滑坡区 世纪 年代即有变形），裂缝逐渐向下延展，受到下部硬质片麻岩的阻挡。在各种外部因素的影响下，片麻岩体最终被剪断，斜坡完全失稳，岩土体堵塞金沙江，形成白格堰塞湖。（）岩体失稳后，斜坡下部的片麻岩被推挤到堆积体前部，停积时因受到挤压而隆起；蛇纹岩碎屑体在滑坡体后部受到滑坡体牵引而松散，地形相对低洼。“”白格滑坡数值模拟结果分析为验证对“”白格滑坡失稳机理的认识，对滑源区进行二维有限元数值模拟分析。表 为计算所采用的主要岩土体参数，图 所示为所建立的失稳前斜坡剖面模型，计算模型考虑区分强风化和新鲜岩体。表 白格滑坡主要物理力学参数 序号岩土体类型重度（）内聚力内摩擦角（）弹性模量泊松比片麻岩蛇纹岩强风化片麻岩强风化蛇纹岩坡表覆盖层图 “”白格滑坡数值计算剖面 “”结合历史遥感影像、现场调查和数值计算结果，将“”白格滑坡的失稳 停积过程描述如下：受构造深层韧性剪切作用，蛇纹岩碎裂成碎粉岩状。金沙江两岸岩体向临空面卸荷回弹变形，蛇纹岩体中形成卸荷裂隙。裂隙逐渐沿软弱面成集束发展，在坡表形成宏观裂缝。坡内裂隙向下扩展，受到硬质片麻岩的阻碍。片麻岩“锁固体”被剪断，变形体失稳。变形岩体在剪断“锁固体”时，破裂面发生偏转，在软硬岩接触部位形成一个巨型凹槽。失稳岩土体停积在河谷中形成堰塞体。远滑源区侧片麻岩分布区因受到后部岩体的推挤而隆起，而堰塞体近滑源区侧蛇纹岩区则相对低洼。岩体高速撞击江面形成涌浪，冲刷对岸坡表壤土和植被，形成涌浪影响区。堰塞湖水位逐渐抬升，从相对低洼的蛇纹岩区漫坝，水流侵蚀软质的蛇纹岩堆积体，形成贯通的泄洪槽。变形体失稳后，后部岩体突然失去支撑，形成了新的变形体。月 日，变形体失稳，再次堵塞金沙江。数值模拟计算结果显示（见图），上部软质的蛇纹岩在重力作用下向斜坡下部蠕变，剪切面沿强弱风化界面向下延展。上部变形体遇到下部硬质岩阻挡后，主应力在软硬岩接触面附近发生应力集中。直到下部硬岩发生大变形，被剪出，如图（）所示，剪出口上部网格隆起，斜坡失稳，剪出口位移超过 。由计算结果可以看到，岩体失稳后滑移面在软硬岩接触带附近变缓，滑面上部倾角约为，下部滑面倾角降低为 左右，局部反倾，在滑源区下部形成一个凹槽。数值模拟计算结果与前文白格滑坡失稳机理的认识相互印证。第 期 杨连伟：金沙江白格滑坡失稳机理及影响因素分析图 “”白格滑坡滑源区位移有限元计算结果 “”白格滑坡失稳后的动力特性分析“”白格滑坡从滑坡源区陡壁后缘至堆积体最前缘的水平距离为 ，垂直高差 。则滑坡的远程运动特性参数，亦即崩落体的等价摩擦系数 ，远小于 （约等于，为国际公认的岩质材料摩擦系数经验值），白格滑坡失稳后进行了远程运动。采用谢德格尔法对滑坡体的运动速度进行估算，计算公式如式（）所示。（）（）式中：为估算崩滑体速度，；为重力加速度，取 ；为崩塌后壁顶点至崩落体沿途计算点高差，；为崩塌后壁顶点至崩落体沿途计算点水平距离，；为崩落体的等价摩擦系数，取 。利用滑坡主剖面（剖面线位置参见图）计算白格滑坡体首次滑动的全程运动速度，结果曲线如图 所示。岩土体在脱离剪出口不远速度就达到 ，此后随着斜坡高程的降低，岩土体的运动速度不断增大，在河谷底部速度达到峰值 。此后，岩土体受到江水和对岸山体的阻挡，在大约 距离内速度急剧降低，直至停积下来，形成堰塞坝。同理，“”白格滑坡岩土体峰值运动速度为 。图 白格滑坡沿途运动速度分布 影响“”金沙江白格滑坡的主要因素白格滑坡位于金沙江峡谷中，山区复杂的地质背景、长期的斜坡卸荷及表生改造作用和