Processing math: 60%
  • 中文核心期刊
  • 中国科技核心期刊
  • 中国科学引文数据库核心期刊

四川泸定昔格达组滑坡灾害运动过程模拟分析

宋德光, 吴瑞安, 马德芹, 郭长宝, 王炀, 倪嘉伟, 李祥

宋德光, 吴瑞安, 马德芹, 郭长宝, 王炀, 倪嘉伟, 李祥. 2023: 四川泸定昔格达组滑坡灾害运动过程模拟分析. 地质通报, 42(12): 2185-2197. DOI: 10.12097/j.issn.1671-2552.2023.12.014
引用本文: 宋德光, 吴瑞安, 马德芹, 郭长宝, 王炀, 倪嘉伟, 李祥. 2023: 四川泸定昔格达组滑坡灾害运动过程模拟分析. 地质通报, 42(12): 2185-2197. DOI: 10.12097/j.issn.1671-2552.2023.12.014
SONG Deguang, WU Ruian, MA Deqin, GUO Changbao, WANG Yang, NI Jiawei, LI Xiang. 2023: Simulation analysis of landslide disaster movement process in Xigeda Formation, Luding County, Sichuan Province. Geological Bulletin of China, 42(12): 2185-2197. DOI: 10.12097/j.issn.1671-2552.2023.12.014
Citation: SONG Deguang, WU Ruian, MA Deqin, GUO Changbao, WANG Yang, NI Jiawei, LI Xiang. 2023: Simulation analysis of landslide disaster movement process in Xigeda Formation, Luding County, Sichuan Province. Geological Bulletin of China, 42(12): 2185-2197. DOI: 10.12097/j.issn.1671-2552.2023.12.014

四川泸定昔格达组滑坡灾害运动过程模拟分析

基金项目: 

国家重点研发计划项目专题《复杂山区典型滑坡物理模拟与危险性评价》 2021YFC3000505-02

中国地质调查局项目《川藏铁路高位远程地质灾害与高边坡稳定性调查》 DD20221816-3

中央科研院所基本科研业务费项目《川藏交通廊道水岩作用下大型深层滑坡灾变》 DZLXJK202009

详细信息
    作者简介:

    宋德光(1999-), 男, 硕士, 地质资源与地质工程专业, 从事工程地质与地质灾害研究。E-mail: sdgsong@163.com

    通讯作者:

    吴瑞安(1991-), 男, 博士, 副研究员, 从事工程地质与地质灾害方面的研究。E-mail: wuruian1991@126.com

  • 中图分类号: P642.22;P642.23

Simulation analysis of landslide disaster movement process in Xigeda Formation, Luding County, Sichuan Province

  • 摘要:

    四川泸定昔格达组以半成岩为主,工程地质特性复杂,在高陡斜坡中常发生浅层蠕滑变形,在强降雨作用下失稳后可转化为泥石流。本文以四川省泸定县海子坪环环村滑坡为例,基于遥感解译、地面调查、数值模拟等方法,对滑坡发育特征、潜在失稳模式和滑坡-泥石流运动过程进行分析。结果表明,环环村滑坡主要发育于昔格达组内,以深度3~5 m的浅层变形为主,整体处于蠕滑变形阶段。滑坡平面上分为强变形区(A区)和弱变形区(B区),体积分别约为5.5×104 m3、5.8×104 m3,在不同降雨条件下,存在仅有A区下滑和A区牵引B区一起下滑并转化为沟道泥石流2种致灾模式。当仅有A区失稳下滑时,最远运动距离可达1325 m,最大堆积厚度为5.2 m,最大运动速度41.6 m/s,滑坡破坏沟口居民区及道路。当A区和B区同时失稳下滑时,最远运动距离可达1345 m,最大堆积厚度为7.7 m,最大运动速度为44.3 m/s,滑体最远能够冲至河流对岸,形成高约3 m的滑坡坝。研究结果对于深化浅层滑坡-泥石流远程致灾效应的认识和防灾减灾具有一定的指导意义。

    Abstract:

    The Xigeda Formation developing in Luding County, Sichuan Province, is mainly composed of semi-diagenetic rocks, which has complex engineering geology properties.Shallow creep deformation of the Xigeda Formation often occurs in the slope with high and steep terrain.It can transform into debris flow after destabilization under the action of heavy rainfall.In this paper, the Huanhuancun landslide in Haiziping developing in the Xigeda Formation was selected as a study case.Based on remote sensing interpretation, field investigation and numerical simulations, we analyzed and calculated the development characteristics of landslides, potential instability modes and the process of landslide-debris flow movement.This study has shown that the Huanhuancun landslide is mainly developed in the Xigeda Formation, mainly in the shallow deformation with a depth of 3~5 m, and is in the creep deformation stage as a whole.The landslide in plane can be divided into strong deformation zone(A) and weak deformation zone(B) with the volume of approximately 5.5×104m3 and 5.8×104m3.Under different rainfall conditions, there are only two forming patterns of disasters: zone A sliding and zone A pulling zone B sliding together and transforming into gully debris flow.When only the zone A slides downslope, the longest movement distance can reach 1325 m, the maximum thickness of deposits is 5.2 m, the maximum velocity is 41.6 m/s, and the sliding body will damage the residential areas and roads.When the zone A and B slides together, the maximum movement distance can reach 1345 m, the maximum thickness of deposits is 7.7 m, the maximum velocity is 44.3 m/s, and the sliding body can rush as far as the other side of the river, forming a landslide dam with a height of about 3 m.The research results have certain guiding significance for deepening the understanding of remote disaster effect of shallow landslide-debris flow and disaster prevention and mitigation.

  • 中国昔格达组主要呈断续条带状、片状分布于大渡河、金沙江、雅砻江等河谷区域(刘慧军等,2004吴焕恒,2010铁永波等,2022阳清青等,2023),主要由粘土、粉砂质粘土、粉砂组成,是一套形成于上新世中晚期的半成岩地层。有学者根据工程特性将其称为沉积型硬土-软岩(蒋复初等,1999张永双等,2000)。该地层具有极强的水敏感性,当含水率超过敏感界线时,力学性质陡然下降,易发生变形破坏(黄绍槟等,2005周平等,2020)。李后强等(1997)通过研究攀西地区昔格达组滑坡特征,提出昔格达组斜坡易滑角度为28°~30°;Zhu et al.(2011)对喇嘛溪沿线17处浅层滑坡进行了调查研究,认为昔格达组饱水是滑坡发生的主要原因。昔格达组不良的工程地质特性,为地质灾害提供了良好的孕生条件。例如,成昆铁路在施工时遇到昔格达组拱顶沉降、围岩变形、边坡滑塌等问题(丁文富等,2017王志杰等,2017);发育于该地层中的四川省米易县林场滑坡,导致成昆铁路长约860 m的线路严重变形,造成行车中断及严重的经济损失(孟祥磊,2015)。有关昔格达组的研究已取得丰硕的成果,多侧重于工程地质特性及滑坡变形破坏机理,而对昔格达组滑坡运动特征研究较少。

    国内外学者采用多种数值模拟软件开展滑坡运动过程与危险性评价。如刘春等(2017, 2019)基于能量守恒定律,利用离散元颗粒接触模型建立了MatDEM软件,为滑坡运动过程研究提供了一种新方法;Poisel et al.(2008)利用PFC3D软件对秘鲁滑坡的运动过程进行了研究,验证了该软件模拟结果的准确性;Roslan et al.(2020)通过分析认为,UDEC软件可以模拟岩石从断裂到抛掷的过程,是研究边坡稳定性的重要方法;吴瑞安等(2018)分析了不同降雨条件下的滑坡失稳概率,并通过DAN3D数值软件对古滑坡复活后的运动过程进行了研究。此外,Ouyang et al.(2013)基于深度连续积分理论,建立了Massflow数值软件,并应用于深圳光明新区滑坡、茂县新磨村滑坡、金沙江白格滑坡等(Ouyang et al.,2017Scaringi et al.,2018Fan et al.,2020),该软件在滑坡运动过程、危险性评价、风险评估等方面均取得良好的效果。相较于其他软件,Massflow具有运算速度快、效率高、支持二次开发等优点,被广泛应用于山地灾害的全过程动力演化研究。

    在开展川西地区地质灾害调查工作时,笔者发现四川泸定海子坪昔格达组内滑坡灾害发育,饱水状态时容易形成浅层溜滑、泥石流,并可能发生远程致灾效应。本文选取海子坪环环村滑坡为研究对象,结合遥感解译、野外地质调查等方法,在查明滑坡空间结构特征和潜在失稳模式的基础上,采用Massflow数值模拟软件对其运动过程进行分析,并评价其危险性,以期对类似可能远程致灾的滑坡科学防控提供参考。

    海子坪位于四川省泸定县,受龙门山断裂带、鲜水河断裂带和大渡河断裂带控制,形成了典型的高山峡谷地貌(邓建辉等,2007郭长宝等,2021梁馨月等,2021秦宇龙等,2021)。受亚热带季风气候及青藏高原冷空气影响,泸定县年平均气温为15.4℃,年平均降雨量约为1080 mm,最大日降雨量72.3 mm,降雨多集中于5—10月,占全年的90%以上(何军,2009杨智勇等,2011)。受河谷深切影响,海子坪昔格达组出露厚度逾400 m,为中国西南地区昔格达组厚度最大的剖面之一(图 1)。其自上而下可分为3段:上段厚约116 m,为粉砂、粉砂质粘土、粘土互层;中段厚约206 m,为厚层状中细砂、中砂、中粗砂,局部含透镜体及细粒层;下段厚约89 m,为薄层状粉砂质粘土与粉砂互层(吴焕恒,2010施云云,2020)。海子坪环环村出露的昔格达组主要为浅黄色粉砂与浅灰色粘土互层(图 2-a),强降雨作用下饱水的昔格达组常见小规模浅层溜滑和坡面泥石流(图 2-b)。昔格达组具有吸水软化、脱水崩解等特点,在含水率较高时极易形成滑坡和滑坡-泥石流灾害,如四川泸定县兴隆镇群发性滑坡-泥石流、干海子滑坡、汉源县喇嘛溪沟多处浅层滑坡-泥石流灾害、攀西地区昔格达组滑坡与泥石流灾害等(刘恒一,1987张永治,1995柴春阳,2008李文杰等,2016)。不同条件下昔格达组抗剪强度差别较大,干燥条件下内摩擦角最高可达45°,粘聚力达105 kPa,但饱和条件下内摩擦角为6°~26°,粘聚力为16~40.3 kPa(周云金等,2000文丽娜等,2005周罕等,2014)。此外,钟成等(2012)通过原位试验与室内试验测得,昔格达组粉砂渗透系数为1.04×10-4~5.01×10-4 cm/s,粘土渗透系数为1.7×10-5~5.29×10-5 cm/s。

    图  1  四川泸定县海子坪工程地质平面图
    Figure  1.  Engineering geological plan of Haiziping area in Luding County, Sichuan Province
    图  2  四川泸定海子坪昔格达组特征
    a—海子坪昔格达组(镜向160°);b—昔格达组饱水后浅层溜滑(镜向80°)
    Figure  2.  Characteristics of the Xigeda Formation in Haiziping, Luding County, Sichuan Province

    环环村滑坡位于泸定县海子坪环环村马桑林沟上游汇水区,距下游沟口居民区水平距离约1150 m。马桑林沟为季节性冲沟,沟道最高点高程约为2550 m,沟口高程约为1570 m,平均纵坡降约为42.5‰,两侧斜坡坡度为30°~50°(图 3),从沟口冲洪积物的堆积特征分析可知,历史上曾多次发生小型泥石流灾害。

    图  3  四川泸定环环村滑坡平面图与发育特征
    a—环环村滑坡平面图;b—环环村滑坡遥感影像(据Google Earth);c—强变形区发育特征(镜向150°);d—弱变形区发育特征(镜向350°)。δ2—闪长岩;Qpgl—冰水堆积物;Qhpal—冲洪积物;NQx—昔格达组;γo2—斜长花岗岩
    Figure  3.  Plan and development characteristics of the Huanhuancun landslide in Luding County, Sichuan Province

    环环村滑坡在平面上呈圈椅状,后缘以环环村公路为界,高程约2060 m,前缘剪出口高程约1990 m,平均坡度约30°,剪出口高出马桑林沟下游沟口约420 m。滑坡向马桑林沟道方向滑动,主滑方向约为280°(图 3-a)。滑坡在平面上可分为2个亚区(图 3-b~d),两者以冲沟为界,左侧强变形区(A区)平均纵长约140 m,平均横宽约130 m,滑体平均厚度约3 m,体积约5.5×104 m3;右侧弱变形区(B区)平均纵长约125 m,平均横宽约155 m,滑体平均厚度约3 m,体积约5.8×104 m3

    环环村滑坡主要发育于昔格达组内,滑体物质为厚8~30 cm的薄层、中厚层状浅黄色粉砂与厚1~3 cm的薄层状浅灰色粘土互层(图 4),地层倾角为15°~30°,属于昔格达组易滑倾角(李后强,1997Hürlimann et al.,2014)。下伏基岩为晋宁期斜长花岗岩。据调查统计,海子坪地区昔格达组滑坡深度多在5 m以内,由于该区昔格达组浅表部节理发育,越往深部,岩土体完整性越好,一定埋深处的浅灰色粘土结构较致密,能起到相对隔水作用。

    图  4  四川泸定环环村滑坡地质剖面
    Figure  4.  Geological profile of Huanhuancun landslide in Luding County, Sichuan Province

    现场访问和调查发现,环环村滑坡的强变形区(A区)于2020年7月发生过滑动,滑动区域纵长约130 m,水平滑动距离约3 m,造成滑坡后缘公路错断,最大垂直错距约2.2 m(图 5-a),居民门前道路破坏(图 5-b),地表可见长85~90 m、宽15~30 cm的拉张裂缝(图 5-c)。该滑坡上目前已安装裂缝计、位移计等监测仪器,但未采取工程治理措施。在降雨条件下,裂缝附近易发生局部汇水,并成为地表水优势入渗通道,加速雨水集中下渗,可能引发滑坡进一步变形破坏。相对A区而言,弱变形区(B区)变形程度较小,地面以小型拉张裂隙为主,环环村公路上可见多处由坡体变形引起的拉张裂缝,裂缝长4~5 m,宽2~5 cm(图 5-d)。

    图  5  四川泸定环环村滑坡破坏特征
    a—错断的公路(镜向80°);b—居民院子下错陡坎(镜向240°);c—滑坡后缘拉张裂缝(镜向260°);d—环环村公路变形破坏(镜向240°)
    Figure  5.  Damage characteristics of the Huanhuancun landslide in Luding County, Sichuan Province

    环环村滑坡发育于凹型汇水区斜坡地带,目前正处于蠕滑变形阶段。昔格达组节理和微裂隙发育,雨水可沿节理或微裂隙入渗。该地层中粉土与粘土富含的粘土矿物亲水性强,遇水具有膨胀性,滑体重度随雨水入渗增大,土体抗剪强度降低,浅表层岩土体易滑动(黄绍槟等,2005)。在强降雨条件下,地形高差和下部沟道可以为滑坡运动提供有利的势能条件和良好通道,近饱水状态的昔格达组一旦发生滑坡,极可能沿沟道向下滑动解体,并形成滑坡-泥石流灾害(刘恒一,1987李树德等,1988张永治,1995彭建兵等,2020)。这类地质灾害通常具有突发性强、历时短、运动速度快、链式成灾等特点(王士革,1999Iverson et al.,2016Nishiguchi et al.,2022),往往能够形成远程致灾效应。

    根据滑坡空间结构和变形破坏特征分析,环环村滑坡主要存在2种潜在破坏失稳模式:强变形区(A区)局部失稳:在强降雨条件下,A区失稳并转化为泥石流;强变形区(A区)和弱变形区(B区)同时失稳:若A区失稳破坏,B区受到的侧向约束力减小,可能受牵引而发生渐进破坏并转化为泥石流。

    在野外地质调查和滑坡潜在失稳模式分析的基础上,采用Massflow数值软件对滑坡不同工况下的滑坡运动过程进行分析。工况一仅为A区失稳下滑并形成泥石流,工况二为A区与B区同时失稳下滑并形成泥石流。

    Massflow软件是采用基于深度积分的连续介质力学原理及改进的MacCormack-TVD有限差分法,满足质量守恒方程及动量守恒方程(公式1~3),同时考虑地形地貌的数值模拟软件,具有精度高、计算速度快等优点(Ouyang et al.,20132015),可以有效计算滑坡、泥石流、堰塞坝溃决等山地灾害的动力演化过程(图 6)。

    图  6  滑坡运动过程数值模拟技术流程图
    Figure  6.  Numerical simulation flow chart of landslide movement process
    ht+hux+hvy=0 (1)
    hut+(hu2+kapgh2/2)x+huvy=kapghzbx(τzx)bρ (2)
    hvt+huvx+(hv2+kapgh2/2)x+huvy=kapghzbx(τzy)bρ (3)

    式中,ρ为流体密度(kg/m3);h为流体高度(m);t为时间(s);uv为分别为x方向和y方向上的流体速度(m/s);kap为土压力系数;g为重力加速度(m/s2);(τzx)b、(τzy)b分别为底部zx方向和zy方向的剪应力(Pa)。

    Massflow软件提供了Coulomb、Manning、Voellmy等不同的基底摩擦模型用于模拟崩滑流灾害运动过程(周琪等,2019殷邦民,2020)。Coulomb模型主要适用于滑坡及碎屑流,Manning模型主要适用于洪水溃坝等(金文祥,2019),Voellmy模型主要适用于泥流或泥石流。结合上述分析,本文采用Voellmy模型对环环村滑坡-泥石流运动过程进行模拟分析,计算公式如下:

    \tau_{\mathrm{b}}=\sigma \mathtt{μ}+\frac{\rho g v^2}{\xi} (4)

    式中,τb为底部剪切应力(kPa);σ为正应力(kPa);μ为摩擦系数;ρ为泥石流容重(kg/m3);ξ为湍流系数(m/s2)。结合前人研究成果(Buser et al.,1980黄绍槟等,2005段学良等,2019林喜珊,2020熊坤勇等,2020)和室内岩土体测试结果,选取了本次数值模拟所需的相关参数(表 1)。

    表  1  数值模拟参数取值
    Table  1.  Numerical simulation parameter values
    岩土类型 含水率(w)/% 泥石流容重(ρ)/(kg·m-3) 摩擦系数(μ) 湍流系数(ξ)/(m·s-2)
    粉质粘土 26 2020 0.1 1000
    下载: 导出CSV 
    | 显示表格

    可按时间间隔将Massflow软件模拟得到的结果输出为ASCII格式文件,并导入ArcGIS软件进行可视化,即可获得不同时刻下堆积体的速度变化情况(图 7图 8)。在数值模拟计算过程中,通过在不同位置布置监测点可以实现滑坡运动监测,但具体网格点的监测难以反映滑坡运动过程中的实际情况。周琪等(2019)提出了基于网格速度贡献值的滑坡整体速度计算公式(公式5),用于判断滑坡运动情况。

    图  7  工况一不同时刻运动速度变化
    Figure  7.  Variation of motion speed at different time of condition 1
    图  8  工况二不同时刻运动速度变化
    Figure  8.  Variation of motion speed at different time of condition 2
    v_a=\sum\limits_{i=1}^n \frac{\bar{\mu}_i V_i^{\prime}}{V}=\sum\limits_{i=1}^n \frac{\overline{\mu_i} \iint s_i d h}{\iint s_i d s d h} (5)

    式中,va为第i秒对应的滑坡整体速度(m/s);μi为第i秒对应的网格平均速度(m/s);V'iV表示第i秒对应坡体的体积与滑坡总体积;si为第i秒对应的网格面积(m2);h为滑坡在网格上方的堆积厚度(m)。

    在数值模拟计算过程中,利用Massflow软件按5 s间隔输出ASCII格式的xy方向的速度。通过ArcGIS软件将ASCII数据转化为栅格数据,利用栅格计算器,结合公式v=SquareRoot(Square(ui)+Square(vi))获得合速度,并通过空间分析工具获得平均速度。根据公式(5)得出滑坡的整体速度(图 9),可知2种工况下滑坡-泥石流整体运动时间约85 s。

    图  9  不同时刻监测点速度变化情况
    a—工况一监测点速度变化;b—工况二监测点速度变化
    Figure  9.  Velocity changes of monitoring points at different time

    环环村滑坡失稳破坏特征和数值模拟监测数据分析表明,2种工况下滑坡转化为泥石流运动过程相似,大致可分为3个阶段:失稳滑动阶段、滑动解体阶段、减速堆积阶段(陈晓清,2006张明等,2010)。0~10 s为滑坡失稳滑动阶段,降雨沿节理、裂隙入渗,粘土矿物吸水具有膨胀性,使孔隙度减小,产生超孔隙水压力,土体抗剪强度降低或消失,产生部分液化,滑坡失稳启动,滑坡前缘剪出口与堆积区较大的高差使重力势能迅速转化为动能,运动速度迅速增加;10~25 s为滑动解体阶段,滑坡启动后沿沟道滑动过程中逐渐解体转化为泥石流,并在此过程中达到最大速度;25~85 s为减速堆积阶段,泥石流流至地形较缓区域逐渐减速,最终堆积于木角沟内。滑坡减速堆积阶段用时占比大,显示滑坡-泥石流灾害链较强的突发性。

    仅A区失稳下滑时,滑坡最大速度出现于监测点A2、A3之间,为41.6 m/s,整体速度最大为14.8 m/s,当滑坡运动到沟口居民区时速度约为26 m/s;A区与B区共同下滑时,最大速度出现于监测点B2、B3之间,为44.3 m/s,整体速度最大为17.5 m/s,当滑坡体运动到沟口居民区时,速度约为35 m/s。

    环环村滑坡下方为马桑林沟主沟,沟道狭窄,滑坡堆积体沿沟道运动时受地形和滑体体积影响,不同时刻的滑体厚度和面积有所区别(图 10图 11)。

    图  10  工况一不同时刻滑坡堆积范围
    Figure  10.  Accumulation range at different time of condition 1
    图  11  工况二不同时刻滑体堆积范围
    Figure  11.  Accumulation range at different time of condition 2

    当仅A区失稳下滑时,滑坡-泥石流运动抵达沟口居民区用时约30 s,运动过程中最大厚度位于堆积区监测点A4处,厚约5.2 m,最远运动距离为1325 m,最终堆积面积约为3.4×104 m2(图 10)。当A区与B区同时失稳时,滑坡转化而成的泥石流到达沟口居民区用时约26 s,最大厚度位于堆积区监测点B1下方,厚约9.5 m,堆积区最大堆积厚度为7.7 m,最远运动距离为1345 m,最终堆积面积约为4.5×104 m2(图 11)。滑坡运动到沟口时速度较大,会对马桑林沟口与木角沟交汇处附近的房屋及道路造成破坏;工况一堆积范围较小,不会造成堵河,工况二可能会堵河,形成高度约3 m的堰塞坝,但长时间堵塞并回水淹没上游的可能性不大。

    结合滑坡运动过程模拟与危险性分析,提出以下几点治理措施:①对滑坡变形区,及时充填已出现的地表裂缝,防止形成降雨入渗优势通道;②在坡表修建截排水渠,减少雨水沿坡面侵蚀和入渗;③开展相应工程勘察,在沟道内适当设置拦挡坝,有效消减洪峰,并修建排导槽,减小对下游建筑和道路造成损害。

    (1) 四川省泸定县环环村滑坡为昔格达组发育的浅层滑坡,位于汇水斜坡地带,滑体平均厚度3~5 m,局部已经发生明显滑动,目前整体处于蠕滑变形阶段,存在进一步失稳下滑的可能。

    (2) 滑坡平面上可以分为强变形区(A区)和弱变形区(B区),A区体积约5.5×104 m3,B区体积约5.8×104 m3,在不同降雨条件下,存在仅有A区下滑和A区牵引B区一起下滑并转化为沟道泥石流2种成灾模式。

    (3) 仅A区失稳下滑时,最远运动距离为1325 m,最大堆积厚度为5.2 m,最大运动速度41.6 m/s,滑坡到达沟口居民区时速度为26 m/s,沟口堆积面积为3.4×104 m2,可直接危害沟口居民点和道路。

    (4) 当A区与B区同时失稳下滑时,最远运动距离为1345 m,最大堆积厚度为7.7 m,最大运动速度44.3 m/s,沟口堆积面积为4.5×104 m2,滑坡会对沟口处居民区及道路造成破坏,并形成高达3 m的滑坡坝,但长时间堵河可能性较小。

    致谢: 中国地质科学院地质力学研究所杨志华、钟宁副研究员,博士研究生张献兵、李彩虹参与部分野外调查工作,西藏大学吴鹏程硕士参与部分数据处理工作,审稿专家提出了宝贵意见和建议,在此一并表示感谢。
  • 图  1   四川泸定县海子坪工程地质平面图

    Figure  1.   Engineering geological plan of Haiziping area in Luding County, Sichuan Province

    图  2   四川泸定海子坪昔格达组特征

    a—海子坪昔格达组(镜向160°);b—昔格达组饱水后浅层溜滑(镜向80°)

    Figure  2.   Characteristics of the Xigeda Formation in Haiziping, Luding County, Sichuan Province

    图  3   四川泸定环环村滑坡平面图与发育特征

    a—环环村滑坡平面图;b—环环村滑坡遥感影像(据Google Earth);c—强变形区发育特征(镜向150°);d—弱变形区发育特征(镜向350°)。δ2—闪长岩;Qpgl—冰水堆积物;Qhpal—冲洪积物;NQx—昔格达组;γo2—斜长花岗岩

    Figure  3.   Plan and development characteristics of the Huanhuancun landslide in Luding County, Sichuan Province

    图  4   四川泸定环环村滑坡地质剖面

    Figure  4.   Geological profile of Huanhuancun landslide in Luding County, Sichuan Province

    图  5   四川泸定环环村滑坡破坏特征

    a—错断的公路(镜向80°);b—居民院子下错陡坎(镜向240°);c—滑坡后缘拉张裂缝(镜向260°);d—环环村公路变形破坏(镜向240°)

    Figure  5.   Damage characteristics of the Huanhuancun landslide in Luding County, Sichuan Province

    图  6   滑坡运动过程数值模拟技术流程图

    Figure  6.   Numerical simulation flow chart of landslide movement process

    图  7   工况一不同时刻运动速度变化

    Figure  7.   Variation of motion speed at different time of condition 1

    图  8   工况二不同时刻运动速度变化

    Figure  8.   Variation of motion speed at different time of condition 2

    图  9   不同时刻监测点速度变化情况

    a—工况一监测点速度变化;b—工况二监测点速度变化

    Figure  9.   Velocity changes of monitoring points at different time

    图  10   工况一不同时刻滑坡堆积范围

    Figure  10.   Accumulation range at different time of condition 1

    图  11   工况二不同时刻滑体堆积范围

    Figure  11.   Accumulation range at different time of condition 2

    表  1   数值模拟参数取值

    Table  1   Numerical simulation parameter values

    岩土类型 含水率(w)/% 泥石流容重(ρ)/(kg·m-3) 摩擦系数(μ) 湍流系数(ξ)/(m·s-2)
    粉质粘土 26 2020 0.1 1000
    下载: 导出CSV
  • Buser O, Frutiger H. Observed maximum run-out distance of snow avalanches and the determ ination of the friction coefficients μ and ξ[J]. Journal of Glaciology, 1980, 26(94): 121-130. doi: 10.3189/S0022143000010662

    Fan X M, Yang F, Subramanian S S, et al. Prediction of a multi-hazard chain by an integrated numerical simulation approach: the Baige landslide, Jinsha River, China[J]. Landslides, 2020, 17(1): 147-164. doi: 10.1007/s10346-019-01313-5

    Hürlimann M, McArdell B W, Rickli C. Field and laboratory analysis of the runout characteristics of hillslope debris flows in Switzerland[J]. Geomorphology, 2014, 232: 20-32.

    Iverson R M, George D L. Modelling landslide liquefaction, mobility bifurcation and the dynamics of the 2014 Oso disaster[J]. Géotechnique, 2016, 66(3): 1-13.

    Nishiguchi Y, Uchida T. Long-runout-landslide-induced debris flow: the role of fine sediment deposition processes in debris flow propagation[J]. Journal of Geophysical Research: Earth Surface, 2022, 127(2).

    Ouyang C J, He S M, Tang C. Numerical analysis of dynamics of debris flow over erodible beds in Wenchuan earthquake-induced area[J]. Engineering Geology, 2015, 194: 62-72. doi: 10.1016/j.enggeo.2014.07.012

    Ouyang C J, He S M, Xu Q, et al. A MacCormack-TVD finite difference method to simulate the mass flow in mountainous terrain with variable computational domain[J]. Computers & Geosciences, 2013, 52: 1-10.

    Ouyang C J, Zhou K Q, Xu Q, et al. Dynamic analysis and numerical modeling of the 2015 catastrophic landslide of the construction waste landfill at Guangming, Shenzhen, China[J]. Landslides, 2017, 14(2): 705-718. doi: 10.1007/s10346-016-0764-9

    Poisel R, Preh A, Hungr O. Run out of landslides-continuum mechanics versus discontinuum mechanics models[J]. Geomechanics and Tunnelling, 2008, 1(5): 358-366. doi: 10.1002/geot.200800036

    Roslan R, Omar R C, Putri R F, et al. Slope stability analysis using Universal Distinct Element Code(UDEC)method[J]. IOP Conference Series: Earth and Environmental Science, 2020, 451: 012081. doi: 10.1088/1755-1315/451/1/012081

    Scaringi G, Fan X M, Xu Q, et al. Some considerations on the use of numerical methods to simulate past landslides and possible new failures: the case of the recent Xinmo landslide(Sichuan, China)[J]. Landslides, 2018, 15(7): 1359-1375. doi: 10.1007/s10346-018-0953-9

    Zhu B L, Wu X Y. An analysis of rainfall-Induced landslide in colluvial and eluvial soils overlying Xigeda strata, southwestern Sichuan, China[J]. Advanced Materials Research, 2011, 250/253: 2682-2688. doi: 10.4028/www.scientific.net/AMR.250-253.2682

    柴春阳. 雅泸高速文武坡喇嘛溪沟重力侵蚀及其对工程影响研究[D]. 西南交通大学硕士学位论文, 2008.
    陈晓清. 滑坡转化泥石流起动机理试验研究[D]. 西南交通大学博士学位论文, 2006.
    邓建辉, 陈菲, 尹虎, 等. 泸定县四湾村滑坡的地质成因与稳定评价[J]. 岩石力学与工程学报, 2007, 26(10): 1945-1950. doi: 10.3321/j.issn:1000-6915.2007.10.001
    丁文富, 张广泽, 宋章. 成昆铁路昔格达地层工程地质特性及对策研究[J]. 铁道工程学报, 2017, 34(4): 1-5. doi: 10.3969/j.issn.1006-2106.2017.04.001
    段学良, 马凤山, 郭捷, 等. 基于Massflow模型的西藏仁布杰仲沟泥石流运动特征分析[J]. 中国地质灾害与防治学报, 2019, 30(6): 25-33.
    郭长宝, 倪嘉伟, 杨志华, 等. 川西大渡河泸定段大型古滑坡发育特征与稳定性评价[J]. 地质通报, 2021, 40(12): 1981-1991. doi: 10.12097/j.issn.1671-2552.2021.12.001
    何军. 大渡河干海子堰塞堆积体成因机制及稳定性研究[D]. 成都理工大学硕士学位论文, 2009.
    黄绍槟, 吉随旺, 朱学雷, 等. 西攀路昔格达地层滑坡分析[J]. 公路交通科技, 2005, 6(S1): 41-44.
    蒋复初, 吴锡浩, 肖华国, 等. 四川泸定昔格达组时代及其新构造意义[J]. 地质学报, 1999, 73(1): 1-6. doi: 10.3321/j.issn:0001-5717.1999.01.001
    金文祥. 西藏定日县达仓沟泥石流危险性评价[D]. 成都理工大学硕士学位论文, 2019.
    李后强, 丁晶, 艾南山. 攀西地区滑坡发育的理论研究[J]. 中国地质灾害与防治学报, 1997, 8(1): 13-17.
    李树德, 曾思伟. 论泥石流的另一种类型——滑坡型泥石流——以甘肃刘家堡泥石流为例[J]. 水土保持学报, 1988, 2(4): 66-71.
    李文杰, 冯文凯, 魏昌利, 等. 地形对泸定县兴隆镇群发性滑坡型泥石流的影响分析[J]. 水利与建筑工程学报, 2016, 14(1): 149-154.
    梁馨月, 曾璐, 葛永刚, 等. 川西高原鲜水河断裂带炉霍—道孚段泥石流分布特征[J]. 地质通报, 2021, 40(12): 2061-2070. doi: 10.12097/j.issn.1671-2552.2021.12.009
    林喜珊. 流动型黄土滑坡数值模拟研究[D]. 兰州大学硕士学位论文, 2020.
    刘春, 范宣梅, 朱晨光, 等. 三维大规模滑坡离散元建模与模拟研究——以茂县新磨村滑坡为例[J]. 工程地质学报, 2019, 27(6): 1362-1370.
    刘春, 张晓宇, 许强, 等. 三维离散元模型的滑坡能量守恒模拟研究[J]. 地下空间与工程学报, 2017, 13(S2): 698-704.
    刘恒一. 攀西地区的滑坡和泥石流的关系[J]. 水土保持通报, 1987, 7(1): 21-26.
    刘惠军, 聂德新. 昔格达地层研究综述[J]. 地球科学进展, 2004, 19(S1): 80-82.
    孟祥磊. 昔格达地层软弱围岩隧道变形机理及支护技术研究[D]. 西南交通大学硕士学位论文, 2015.
    彭建兵, 王启耀, 庄建琦, 等. 黄土高原滑坡灾害形成动力学机制[J]. 地质力学学报, 2020, 26(5): 714-730.
    秦宇龙, 吴建亮, 詹涵钰, 等. 川西甘孜地区活动断裂与地质灾害分布相关性探讨[J]. 地质力学学报, 2021, 27(3): 463-474.
    施云云. 大渡河泸定段海子坪昔格达组的宇生核素等时线埋藏测年及地貌意义[D]. 南京师范大学硕士学位论文, 2020.
    铁永波, 张宪政, 卢佳燕, 等. 四川省泸定县Ms6.8级地震地质灾害发育规律与减灾对策[J]. 水文地质工程地质, 2022, 49(6): 1-12.
    王士革. 山坡型泥石流的危害与防治[J]. 中国地质灾害与防治学报, 1999, 10(3): 46-51.
    王志杰, 周平, 徐海岩, 等. 昔格达地层隧道围岩灾变特征及致灾因子研究[J]. 铁道工程学报, 2017, 34(11): 67-75.
    文丽娜, 朱学雷, 白志勇, 等. 西攀高速公路新九地区昔格达地层岩土特性[J]. 公路, 2005, 7(7): 145-148.
    吴焕恒. 西昌某边坡昔格达组地层蠕变试验及其边坡蠕变变形分析[D]. 西南交通大学硕士学位论文, 2010.
    吴瑞安, 张永双, 郭长宝, 等. 川西松潘上窑沟古滑坡复活特征及危险性预测研究[J]. 岩土工程学报, 2018, 40(9): 1659-1667.
    熊坤勇, 范宣梅, 杨帆, 等. 九寨沟九寨天堂泥石流沟特征分析与模拟预测[J]. 科学技术与工程, 2020, 20(27): 10989-10995.
    阳清青, 余秋兵, 张廷斌, 等. 基于GDIV模型的大渡河中游地区滑坡危险性评价与区划[J]. 中国地质灾害与防治学报, 2023, 34(5): 130-140.
    杨智勇, 吉锋. 大渡河干海子滑坡稳定性研究[J]. 防灾减灾工程学报, 2011, 31(2): 201-206.
    殷邦民. 尼续村高速远程滑坡流态化堆积与运动学特征研究[D]. 西南交通大学硕士学位论文, 2020.
    张明, 胡瑞林, 殷跃平, 等. 滑坡型泥石流转化机制环剪试验研究[J]. 岩石力学与工程学报, 2010, 29(4): 822-832.
    张永双, 曲永新. 硬土-软岩的厘定及其判别分类[J]. 地质科技情报, 2000, 19(1): 77-80.
    张永治. 攀西地区昔格土综述[J]. 攀枝花大学学报, 1995, 12(2): 75-78.
    钟成, 范德平. 川南昔格达岩组工程地质特性研究[J]. 四川水力发电, 2012, 31(1): 97-99.
    周罕, 曹平, 张科. 昔格达组黏土岩和粉砂岩现场直剪试验研究[J]. 中南大学学报(自然科学版), 2014, 45(10): 3544-3550.
    周平, 王志杰, 侯伟名, 等. 昔格达地层隧道局部浸湿失稳特征及突变预测研究[J]. 岩土工程学报, 2020, 42(3): 503-512.
    周琪, 许强, 周书, 等. 基于数值模拟的突发型黄土滑坡运动过程研究——以黑方台陈家8#滑坡为例[J]. 山地学报, 2019, 37(4): 528-537.
    周云金, 曾联明. 红格提水工程二级泵站昔格达地层特性及坡体变形成因分析[J]. 水电站设计, 2000, 16(2): 61-65.
  • 期刊类型引用(8)

    1. 李祥,吴瑞安,郭长宝,倪嘉伟,王炀,李彩虹,宋德光. 川西泸定昔格达地层黏土岩工程地质特性研究. 西北地质. 2024(02): 262-274 . 百度学术
    2. 陈锡锐,刘虹强,杨剑红,艾其凯,钟波,蔡国军. 四川宝兴张家沟危岩体稳定性及运动学分析. 中国地质灾害与防治学报. 2024(02): 81-89 . 百度学术
    3. 李金秋,张永双,任三绍,冉丽娜. 金沙江上游扎马古滑坡复活特征及堵河危险性分析. 水利学报. 2024(04): 481-492 . 百度学术
    4. 张友谊,杨经纬,谭亦函. 汉源昔格达组半成岩微观结构与宏观力学的相关性. 科学技术与工程. 2024(18): 7542-7551 . 百度学术
    5. 纪佑军,熊军,蒋国斌,王泽根. 考虑应变软化的鸡场镇降雨型滑坡数值分析. 水文地质工程地质. 2024(04): 178-188 . 百度学术
    6. 向兰兰,徐恒志. 达木斯乡英巴格村滑坡风险性评价. 科技创新与应用. 2024(23): 105-109 . 百度学术
    7. 刘东升,吴越,李珂,王艳磊. 基于失稳概率的单体滑坡灾害风险定量评价及其工程应用. 中国地质灾害与防治学报. 2024(04): 67-74 . 百度学术
    8. 陈锡锐,陈思尧,杨剑红,刘虹强,袁兆平,朱国宝,谢晓文,蔡国军. 四川什邡市地质灾害基本特征与易发性评价. 中国地质灾害与防治学报. 2024(06): 153-163 . 百度学术

    其他类型引用(3)

图(11)  /  表(1)
计量
  • 文章访问数:  963
  • HTML全文浏览量:  290
  • PDF下载量:  112
  • 被引次数: 11
出版历程
  • 收稿日期:  2022-03-28
  • 修回日期:  2022-05-01
  • 网络出版日期:  2024-01-03
  • 刊出日期:  2023-12-14

目录

/

返回文章
返回