ISSN 1003-8035 CN 11-2852/P
  • 中国科技核心期刊
  • CSCD收录期刊
  • Caj-cd规范获奖期刊
  • Scopus 收录期刊
  • DOAJ 收录期刊
  • GeoRef收录期刊
欢迎扫码关注“i环境微平台”

抚顺西露天矿弱层强度衰减特性及边坡滑移大变形规律

陶泽, 孙闯, 金淳哲, 刘晓龙, 李东

陶泽,孙闯,金淳哲,等. 抚顺西露天矿弱层强度衰减特性及边坡滑移大变形规律[J]. 中国地质灾害与防治学报,2023,34(3): 31-39. DOI: 10.16031/j.cnki.issn.1003-8035.202205016
引用本文: 陶泽,孙闯,金淳哲,等. 抚顺西露天矿弱层强度衰减特性及边坡滑移大变形规律[J]. 中国地质灾害与防治学报,2023,34(3): 31-39. DOI: 10.16031/j.cnki.issn.1003-8035.202205016
TAO Ze,SUN Chuang,JIN Chunzhe,et al. Characteristics of strength reduction in the weak layer and large-slip displacement of the cut slope on the Fushun west open-pit mining area[J]. The Chinese Journal of Geological Hazard and Control,2023,34(3): 31-39. DOI: 10.16031/j.cnki.issn.1003-8035.202205016
Citation: TAO Ze,SUN Chuang,JIN Chunzhe,et al. Characteristics of strength reduction in the weak layer and large-slip displacement of the cut slope on the Fushun west open-pit mining area[J]. The Chinese Journal of Geological Hazard and Control,2023,34(3): 31-39. DOI: 10.16031/j.cnki.issn.1003-8035.202205016

抚顺西露天矿弱层强度衰减特性及边坡滑移大变形规律

基金项目: 国家重点研发计划项目(2017YFC1503101)
详细信息
    作者简介:

    陶 泽(1996-),男,硕士,主要从事岩土力学数值模拟方面的研究工作。E-mail:2627137567@qq.com

    通讯作者:

    孙 闯(1983-),男,博士,教授,主要从事岩石力学试验与数值计算方面的研究工作。E-mail:sunchuang88@163.com

  • 中图分类号: P642.22

Characteristics of strength reduction in the weak layer and large-slip displacement of the cut slope on the Fushun west open-pit mining area

  • 摘要: 软弱夹层的强度衰减特性是顺层岩质边坡失稳的关键因素之一。选取抚顺西露天矿南帮边坡为研究对象,通过UDEC数值模拟方法,结合离散元理论,建立南帮顺层边坡二维离散元模型,开展南帮边坡弱层强度衰减特性及滑坡大变形规律研究。通过不同围压下的三轴试验,拟合得到弱层残余强度随时间衰减曲线,并通过FISH语言实现数值计算过程中弱层强度随拟合方程进行衰减。数值模拟结果表明:坡脚处最先产生变形,随着开挖的不断进行,弱层强度逐渐衰减,坡顶后缘被拉裂,原有应力平衡被破坏,滑坡体沿弱层产生滑动。通过建立回填数值模型进行边坡稳定性得知,采取压脚回填100 m以后,边坡稳定系数得到了明显的提高,边坡处于基本稳定状态。
    Abstract: The strength attenuation characteristics of the weak interlayers are one of the key factors leading to instability of the bedding rock slopes. In this study, the Nanbang slope of Fushun west open-pit mine was selected as the research object. A two-dimensional discrete element model of the Nanbang bedding slope was established using the UDEC numerical simulation method in combination with the discrete element theory, and the strength attenuation characteristics and large deformation law of the weak layer were studied. The residual strength decay curve of the weak layer over time was obtained by fitting the results of triaxial tests conducted under different confining pressures, and the FISH language was used to implement the attenuation of weak layer strength during numerical calculations according to the fitting equation. The numerical simulation results showed that deformation first occurred at the toe of the slope. As excavation continued, the strength of the weak layer gradually decreased, resulting in the trailing edge of the slope top being pulled apart, and the original stress balance being destroyed. Eventually, the landslide body slides along the weak layer. By establishing a backfilling numerical model for slope stability, it was found that after backfilling with a presser foot for 100 m, the slope stability coefficient had significantly improved, and the slope was in a basically stable state.
  • 滑坡是极具破坏性的地质灾害之一,给人类社会造成了巨大的经济损失和人员伤亡,严重威胁人类生命财产安全,制约人类社会的可持续发展[12]。开展区域滑坡危险性评价,识别滑坡高危险区和阐明地质灾害的潜在分布对地质灾害防治具有重要意义,同时也是决策者预防和减少地质灾害损失的重要手段[3]

    滑坡危险性评价经历了从定性模型到半定量模型再到定量模型的发展过程[4]。定性模型主要是依据专家的知识和野外经验,受主观性限制较多;定量模型主要使用统计或数学建模技术,评价结果具有客观性[57]。随着科学技术的发展和定量风险管理的需要,滑坡定量评价模型成为主要的研究方向[810]。常用的定量模型主要包括信息量模型[1112]、确定系数法[13]、逻辑回归模型[1416]、地理加权回归模型[17]、支持向量机模型[1819] 、人工神经网络[20]和随机森林[2122]等模型。上述模型为区域滑坡危险性评价提供了重要参考,已在省、市和县域尺度的地质灾害危险性评价中得到广泛应用[2324]。但在评价过程中,多数研究在定性确定评价因子后,直接应用于危险性评价[25],很少开展评价因子的定量优选工作;各评价因子分级量化主观性较强,评价因子权重赋值多依靠专家经验或先验知识,缺少客观依据[26];运用常权栅格叠加各评价因子时,因子间存在抵消问题和极值区域埋没问题[27]。因此,结合区域滑坡成灾特征建立科学的评价指标体系,构建有效的评价模型是区域滑坡危险性评价的关键问题[28]

    大渡河中游地区位于青藏高原东缘的川西地区,区域地形地貌和地质条件复杂,生态环境较脆弱,气候复杂多变,是我国滑坡地质灾害高发地区[29]。由于地理环境恶劣,地质灾害调查和监测难度大,地质灾害综合防治能力较薄弱。综上,本文以大渡河中游地区为研究区,构建一种基于地理探测器、信息量法和变权栅格叠加的评价模型(geographic detector , information quantity and variable weight raster overlay model,GDIV)。通过地理探测器定量描述各因子对滑坡发生的影响力,基于影响力的显著程度实现评价因子的优选,同时以因子探测结果作为因子定权基础,增强权重的客观性,解决了定性评价因子选择主观性强的问题。采用变权栅格叠加评价因子信息量值,解决了评价因子叠加存在的抵消问题。基于GDIV模型开展区内滑坡危险性评价,评价模型为区域滑坡地质灾害评价提供新参考,评价结果为研究区地质灾害防治提供决策支持。

    大渡河中游地区地处四川省阿坝州和甘孜州交界处,区域面积约8 179.6 km2。该区地质构造复杂,包括滇藏“歹”字形构造体系和金汤弧形构造体系,出露的地层主要有三叠系、志留系和泥盆系。地处青藏高原东缘地势一级阶梯向二级阶梯过渡地带,区内海拔1 787~5 782 m,相对高差大。气候以高原型季风气候为主,年平均气温18.4℃,年均降雨量788 mm。区内水系发达,以大渡河为主流,包括革什扎河和东谷河等支流,河流汇集于丹巴县。在复杂的地质地貌条件、降雨和地震等因素影响下,地质灾害发生频繁。大渡河中游地区地质灾害类型主要包括滑坡、崩塌和泥石流,地质灾害详查数据显示,截至2021年12月,区内分布滑坡地质灾害隐患点313处(图1)。滑坡发育规模以中小型为主,大型规模的地质灾害6处,巨型规模的滑坡1处。

    图  1  大渡河中游地区滑坡分布图和地质条件背景图
    Figure  1.  Map of landslide distribution and geological conditions in the middle reach area of Dadu River

    本研究所需的基础数据包括大渡河中游地区滑坡点、地层岩性、土壤类型、数字高程模型和地震动参数等。其中滑坡隐患点数据来源于四川省国土空间生态修复与地质灾害防治研究院1∶50 000地质灾害详查数据;地层岩性和地质构造数据来源于1∶25万炉霍县幅、康定幅和马尔康幅地质图和相关地质报告;降雨量数据来自国家气象科学数据中心—中国气象数据网(http://data.cma.cn/);土壤类型数据来源于中国科学院资源环境科学数据中心(http://www.resdc.cn),基于2010—2020年降雨量插值生成年平均降雨量;数字高程模型(digital elevation model,DEM)和2020年8月两景Landsat8 OLI_TIRS数据来源于地理空间数据云(http://www.gscloud.cn/),空间分辨率30 m。利用DEM生成坡度、坡向、地形湿度指数(TWI)和径流强度指数(SPI),基于Landsat8影像计算归一化植被指数(NDVI)。地震动参数源自中国地震动参数区划图GB 18306—2015(http://www.gb18306.net/)。同时将年平均降雨量、地震动参数等栅格数据统一空间分辨率为30 m,地理参考为WGS_1984_UTM_zone_47N。

    滑坡地质灾害是由内部因素和外部因素共同作用的结果[30],本研究基于区域滑坡发育规律,从地质特征、地震、地形地貌、地表水系、降雨和人类活动等6个方面, 选取13个影响因子构建大渡河中游地区滑坡危险性初级评价指标体系(图2)。其中工程地质岩组分为4类,分别为坚硬岩、较坚硬岩、较软岩、松散土。坚硬岩主要为花岗岩、闪长岩、石英砂岩和石英岩等;较坚硬岩主要为板岩、石英岩状砂岩和硅质岩等;较软岩主要为页岩、千枚岩和粉砂岩等;松散土主要是第四系河漫滩等。

    图  2  大渡河中游地区滑坡危险性初级评价指标体系分级图
    Figure  2.  Grading chart of the primary hazard assessment index system for landslides in the middle reach area of Dadu River Basin

    将滑坡点密度作为因变量,高程、坡度、坡向、地震动参数、土壤类型、与道路距离和与断层距离等13个因子作为自变量,运用地理探测器的因子探测器和交互探测器对因子进行计算分析。为消除因子类型差异的影响和优化各级因子对地质灾害影响程度的体现,先利用信息量法对自变量进行处理[31]

    I(Xi,Y)=lnP(B/Xi)P(B)=lnNi/NSi/S (1)

    式中:Xi——成灾因素X中的第i区间;

    Y——成灾因素X中第i区间滑坡发生的信息量值;

    B——滑坡事件;

    Ni——研究区内包含评价因素Xi的单元数且存在滑坡的单元数;

    Si——研究区内包含评价因素Xi的单元数;

    N——研究区存在滑坡单元总数;

    S——研究区中评价单元的总数。

    因子探测器的计算如下[3233]

    q=1h=1LNhσh2Nσ2=1SSWSST (2)
    SSW=h=1lNhσh2,SST=Nσ2 (3)

    式中:h=1, 2, ···, L——变量Y或因子X的分层;

    NhN——层h和全区的单元数;

    σ2hσ2——层h和全区的Y值的方差;

    SSWSST——层内方差之和和全区总方差;

    q ——因子对滑坡发育的解释力,值域为[0,1],值越大说明对滑坡的影响力越大,因子对滑坡发育的解释力越强。

    交互探测器是评估因子X1X2共同作用时是否会增加或减弱对因变量Y的解释力,或这些因子对Y的影响是相互独立的。因子间若存在交互作用便能通过对比q值体现,交互作用方式的判断依据如表1所示。

    表  1  交互作用探测器因子关系
    Table  1.  Factor relationships of interaction detectors
    因子关系交互作用
    q(X1X2)<Min(q(X1), q(X2))非线性减弱
    Min(q(X1), q(X2))< q(X1X2)< Max (q(X1), q(X2))单因子非线性减弱
    q(X1X2)> Max (q(X1), q(X2))双因子增强
    q(X1X2)= q(X1)+q(X2)独立
    q(X1X2)> q(X1)+q(X2)非线性增强
    下载: 导出CSV 
    | 显示表格

    LRI(logistic regression and informative)模型是逻辑回归模型和信息量模型的组合模型,逻辑回归模型是一种对二元分类变量进行回归分析的统计方法,定量描述二元因变量和一系列自变量之间的关系。在危险性评价时,选取评价因子为自变量,以信息量模型计算出的各个指标等级的信息量值表征各评价因子,以区内是否存在滑坡为因变量(1代表存在,0代表不存在)。构建逻辑回归函数:

    P=11+e(α+β1X1++βnXn) (4)

    式中:P——滑坡发生的概率,范围为0~1;

    α——截距;

    β——回归系数;

    i——评价因子种类数目。

    将式(5)两边取自然对数ln[P/1−P],将评价因子Xi(i=1, 2, ···, n)作为自变量,建立线性回归方程:

    ln(p1p)=α+β1X1+β2X2++βnXn=α+βX (5)

    变权栅格叠加是变权理论在栅格叠加分析中的应用,该方法根据各叠加图层的初始权重及叠加栅格的像元值,通过一定规则对极值像元进行惩罚型、激励型或混合型变权,可以凸显个别极值栅格的效应[28],实现层内因子重要区间的最大化应用。变权栅格叠加的实现主要包括状态变权函数确定和归一化变权权重计算,其过程如图3所示。图3中,W为各评价因子的初始权重,S(X)是根据状态变权函数获得各图层每个像元的状态变权值,W(X)是各个栅格像元的状态变权值和原始权重值乘积的归一化值,M是所有叠加图层各个栅格像元值和变权权重乘积的综合值。本文的状态变权函数采用局部激励型变权,当信息量值>0.6时,相应权重翻倍,即调权幅度为200%;当信息量值≤0.6时,权重保持不变,即调权幅度为100%。

    图  3  变权栅格叠加过程
    Figure  3.  The variational raster overlay process

    GDIV模型是本文提出的地理探测器、信息量法和变权栅格叠加的组合模型,该模型定量描述各因子对滑坡发生的影响力,并以定量结果作为因子定权基础,增强了权重的客观性;同时以变权栅格叠加因子值,消除了因子叠加存在的抵消问题。该模型的实现主要有四步:(1)基于信息量法,对初选因子分级和量化,增强因子各级的差异性;(2)基于地理探测器的因子探测和交互探测结果,选取影响力较强的因子作为评价因子;(3)构建权重方程,见式(6),归一化各评价因子q值以获得权重;(4)确定状态变权函数,生成各因子的变权权重,利用变权栅格叠加实现危险性评价。采用自然间断法将评价结果划分为地质灾害极高、高、中、低和极低5个等级的危险区,实现滑坡危险性评价区划。具体流程如图4所示。

    图  4  GDIV模型计算流程图
    Figure  4.  The flowchart of GDIV model calculation process
    ωi=qii=1nqi (6)

    式中:ωi——归一化各因子q值所得的权重值;

    qi——各因子因子探测q值。

    因子探测结果表明,地形地貌、地质特征、人类活动、地表水系、地震和降雨6个方面的因子对地质灾害发生的影响力较大,但各因子对地质灾害发生的解释力存在差异(表2)。高程、土壤类型、年平均降雨量、地震动参数、与河流距离、工程地质岩组和与道路距离对滑坡发生的影响力较显著,这7个因子是滑坡发育的重要因子,其q值均大于0.1,其中高程的平均解释力度为58.3%,土壤类型的平均解释力为41.5%,则高程和土壤类型是研究区滑坡的重要背景因子。与断层距离、归一化植被指数、坡向、径流强度指数和坡度对滑坡发生的解释力度相对一般,q值仅介于0~0.1。地形湿度指数p值大于0.05,未通过显著性检验,该因子对大渡河中游地区滑坡发生的影响程度较小,在本研究中为次要因子。

    表  2  滑坡初级评价指标q值统计
    Table  2.  Statistical analysis of primary evaluation index q-values for landslides
    类别指标qp
    地质特征工程地质岩组(X10.1560.000
    与断层距离(X20.0870.000
    地震地震动参数(X30.1640.000
    地形地貌高程(X40.5830.000
    坡度(X50.0210.023
    坡向(X60.0380.003
    地形湿度指数(X70.0170.297
    归一化植被指数(X80.0720.000
    土壤类型(X90.4150.000
    地表水系与河流距离(X100.1580.000
    径流强度指数(X110.0320.015
    降雨年平均降雨量(X120.1820.000
    人类活动与道路距离(X130.1150.000
    下载: 导出CSV 
    | 显示表格

    基于因子探测器中对滑坡发生解释力较强的7个因子,采用交互探测器分析两因子交互作用后对滑坡发育的解释力(表3)。结果显示,21对交互因子中,双因子增强16对,非线性增强5对,不存在相互独立或者非线性减弱的关系,与单因子对滑坡发育的解释力相比,两因子交互作用明显增强了对滑坡发育的解释力。交互作用解释力最高的是高程与工程地质岩组(q=0.736),其次是高程与地震动参数(q=0.676),解释力均高于65% ,两者均是滑坡发生的重要因子组合。高程与土壤类型、工程地质岩组、与道路距离、地震动参数、年平均降雨量、与河流距离均呈现双因子增强作用类型,解释力均在59.6%以上,说明高程增强了土壤类型、工程地质岩组、与道路距离、与河流距离、地震动参数、年平均降雨量对滑坡发生的影响力。土壤类型与地震动参数、工程地质岩组的解释力均大于53%,是滑坡发育的较重要因子组合。从各因子的交互作用结果来看,高程、土壤类型和工程地质岩组等7个因子的协同作用是大渡河中游地区滑坡发生的主导交互作用方式,也是滑坡形成的重要影响因子。

    表  3  部分滑坡初级评价指标交互作用
    Table  3.  Interactions of primary evaluation indicators for landslides
    Xi∩Xjq(Xi)q(Xj)q(Xi∩Xj)q(Xi)+q(Xj)交互类型
    X4∩X10.5830.1560.7360.739双因子增强
    X3∩X40.1640.5830.6760.747双因子增强
    X9∩X40.4150.5830.5960.998双因子增强
    X10∩X40.1580.5830.6030.741双因子增强
    X13∩X40.1150.5830.5970.698双因子增强
    X12∩X40.1820.5830.6720.765双因子增强
    X9∩X30.4150.1640.5370.579双因子增强
    X9∩X10.4150.1560.5550.571双因子增强
    X9∩X100.4150.1580.4340.573双因子增强
    X9∩X130.4150.1150.4280.53双因子增强
    X9∩X120.4150.1820.5270.597双因子增强
    X10∩X30.1580.1640.3120.322双因子增强
    X10∩X10.1580.1560.3440.314非线性增强
    X13∩X30.1150.1640.2760.279双因子增强
    X13∩X10.1150.1560.2780.271非线性增强
    X3∩X10.1640.1560.3290.320非线性增强
    X13∩X100.1150.1580.2260.273双因子增强
    X10∩X120.1580.1820.3430.340非线性增强
    X13∩X120.1150.1820.2920.297双因子增强
    X3∩X120.1640.1820.2690.346双因子增强
    X12∩X10.1820.1560.3480.338非线性增强
    下载: 导出CSV 
    | 显示表格

    通过地理探测器分析了研究区内高程、土壤类型和工程地质岩组等13个因子与滑坡发生频率变化关系, 结合因子探测器和交互探测器结果,最终选定影响力较强和交互作用显著的7个因子作为大渡河中游地区滑坡危险性评价因子,即高程、年平均降雨量、土壤类型、地震动参数、与道路距离、工程地质岩组和与河流距离。

    通过式(1)得到各评价因子信息量值(表4),其中高程值<2700 m的区间信息量值大于2,表明区内滑坡主要发生在高程值<2700 m的区间;土壤类型为淋溶土和人为土的区间信息量值大于1,表明区内滑坡主要发生在这两类土壤中。工程地质岩组通过基岩和堆积体的物理力学性质来影响滑坡发生,区内较软岩的信息量值大于1.8,意味着较软岩分布区更有利于滑坡的发生。归一化植被指数、与河流距离、地震动参数和与道路距离是滑坡形成的重要因子,年平均降雨量在750~775时,信息量值最大;道路在建设中因剥离地表土壤和松散岩石对周边环境产生影响,区内滑坡发生概率与道路距离成负相关关系。区内与道路距离小于100 m的缓冲区信息量值最高,随着与道路距离的增加信息量值依次递减。

    表  4  危险性评价因子分级与信息量值
    Table  4.  Grading and information value of hazard evaluation factors
    评价因子分级信息量值评价因子分级信息量值
    高程/m<2 7002.058年平均
    降雨量/mm
    <750−0.557
    2 700~3 2001.308750~7750.438
    3 200~3 600−1.37775~800−1.014
    3 600~4 000−2.445800~840−0.055
    4 000~4 400−3.76840~880−0.404
    > 4400>880−0.231
    土壤类型淋溶土1.685地震动
    参数
    <0.10.151
    半淋溶土0.1~0.150.464
    初育土−3.9210.15~0.2−1.059
    高山土0.1070.2~0.3
    人为土1.429与道路
    距离/m
    <1001.500
    铁铝土0.890100~2001.227
    与河流
    距离/m
    <400−1.204200~3001.148
    400~800−0.826300~4001.053
    800~1 200−0.025400~5000.789
    1 200~1 6000.004>500−0.335
    1 600~2 0000.577
    >2 0001.038
    工程地质
    岩组
    坚硬岩0.023
    较坚硬岩0.443
    较软岩1.878
    松散土类−1.086
    下载: 导出CSV 
    | 显示表格

    基于大渡河中游地区各评价因子的信息量值,运用LRI模型分析得到各因子的逻辑回归系数值(表5)。结果表明,各评价因子的显著性水平值均小于0.05,通过显著性检验。各评价因子的回归系数均为正值,表明各评价因子对滑坡发生均具有促进作用。高程的回归系数最大(B=4.992),区内地质灾害的发生与高程相关程度最大。其次是土壤类型(B=3.001)、工程地质岩组(B=1.606)、年平均降雨量(B=1.103)、与道路距离(B=0.995)和地震动参数(B=0.802),回归系数均大于0.8,对滑坡的发生具有较强的影响。基于各因子的回归系数构建线性回归方程,实现区内滑坡危险性评价,见图5(a)。评价结果表明,随着危险等级的降低对应区内的面积逐步增加,其中极高和高危险的面积占比为18.1%,分布295处滑坡,随着危险等级的降低对应区内的滑坡点也逐渐减少。

    表  5  滑坡危险性评价因子逻辑回归分析结果
    Table  5.  Results of logistic regression analysis for landslide hazard evaluation factors
    评价因子BSEWalddfsigExp(B)
    高程4.9920.55182.21010.000147.24
    土壤类型3.0010.55029.78510.00020.110
    工程地质岩组1.6060.8373.38710.0004.666
    年平均降雨量1.1030.3798.46810.0003.013
    与道路距离0.9950.3962.57310.0002.435
    地震动参数0.8020.4691.65710.0001.830
    与河流距离0.1480.3985.25910.0010.739
    常数−7.1320.696104.81510.0000.001
      注:B为模型中各变量的回归系数、SE是标准差、Wald是卡方统计、Sig为显著性水平,dfExp(B)为逻辑回归的结果参数。
    下载: 导出CSV 
    | 显示表格
    图  5  滑坡危险性区划图
    Figure  5.  Landslide hazard zoning map

    利用式(6)获得各评价因子权重(表6),经变权栅格叠加实现基于GDIV模型的滑坡危险性评价(图5b)。结果表明,各评价因子的权重大小顺序为:高程(0.329)>土壤类型(0.234)>年平均降雨量(0.103)>地震动参数(0.092)>与河流距离(0.089)>工程地质岩组(0.088)>与道路距离(0.065)。研究区极高和高危险区分布滑坡共302处,占总地质灾害的96.5%。其中极高危险区滑坡数量占比为65.5%,高危险区滑坡数量占比为31%,中等及以下危险区滑坡发育较少。随滑坡危险等级的降低,对应区内的滑坡点数量也依次减少,评价结果与实际情况相符。滑坡危险性等级从低到高的各区面积占比分别是25.2%(极低)、40%(低)、13.2%(中)、11.2%(高)和10.5%(极高)。GDIV模型充分考虑因子间的抵消效应和极值埋没问题,变权栅格叠加评价因子的综合信息量值在极值区间占比较小,极低危险区面积占比小于低危险区。极高危险区和高危险区主要沿东谷河、大渡河和革什扎河两侧分布,区内海拔小于3200 m,以淋溶土为主,人类活动频繁。低危险区和极低危险区主要分布在水系发育程度较低、海拔相对较高和距交通干线相对较远的地区,集中分布在大渡河中游地区的阿科里乡、毛日乡、丹东乡、银恩乡、七美乡和协德乡等地。

    表  6  滑坡危险性评价因子权重值
    Table  6.  Weight values of landslide hazard assessment factors
    因子q权重
    高程0.5830.329
    土壤类型0.4150.234
    年平均降雨量0.1820.103
    地震动参数0.1640.092
    与河流距离0.1580.089
    工程地质岩组0.1560.088
    与道路距离0.1150.065
    下载: 导出CSV 
    | 显示表格

    基于ROC曲线对GDIV模型和LRI模型滑坡危险性评价结果进行精度验证和对比分析。研究区共分布滑坡点313处,在研究区生成313处随机点作为非灾害样本点合并绘制ROC曲线(图6)。LRI模型的AUC值为0.903,标准误差为0.013,在满足显著性条件下AUC值最高可大0.929。GDIV模型的AUC值为0.917,标准误差为0.012,在满足显著性条件下AUC值最高可大0.941。两种模型的评价效果均较好,但GDIV模型的评价精度相对较高。由此表明基于GDIV模型的滑坡危险性评价方法在大渡河中游地区的实用性更好,能较好地评估滑坡危险情况。

    图  6  滑坡危险性评价结果ROC曲线
    Figure  6.  ROC curve of landslide hazard evaluation results

    本文以大渡河中游地区为研究区,初选高程、土壤类型和与河流距离等13个因子构建初级评价指标体系,利用地理探测器和信息量法实现因子的筛选和优化,构建GDIV模型开展研究区滑坡危险性评价,为验证模型的精度引入LRI模型,基于ROC曲线对比分析和验证两种模型的评价精度,得出以下结论。

    (1)地理探测器结果表明,高程(q=0.583)和土壤类型(q=0.415)是区内滑坡危险性的重要背景因素,高程与工程地质岩组(q=0.736)、高程与地震动参数(q=0.676)、高程与年平均降雨量(q=0.672)的协同作用对滑坡的发生具有重要影响。

    (2)基于GDIV模型的大渡河中游地区滑坡危险性评价结果显示,区内大部分地区为中度及以下危险区,面积占比为78.3%;而极高和高危险区面积占比为21.7%。极高危险区主要集中分布在河谷地区,工程地质岩组为松散土类均分布在该区内,土壤类型为淋溶土;高危险区要分布在东谷河、大渡河和革什扎河两侧极高危险区以外的区域;中危险区主要集中在海拔3200~4000 m范围内,极低和低危险区主要集中在人口稀少的高海拔区域。

    (3)GDIV模型(AUC=0.917)评价精度高于LRI模型(AUC=0.903),在极高危险区和高危险区内GDIV模型所含滑坡点总量较LRI模型多7处。GDIV模型评价结果更符合区内滑坡分布情况,该模型适用于大渡河中游地区滑坡危险性评价。

  • 图  1   抚顺西露天矿全景图

    Figure  1.   Aerial panoramic view of Fushun west open-pit mine

    图  3   GDS高精度软岩流变仪

    Figure  3.   GDS high precision soft rock rheometer

    图  2   千台山滑坡裂缝

    Figure  2.   Development of landslide cracks in Qiantai Mountain

    图  4   岩石试样

    Figure  4.   Rock samples

    图  5   常规三轴压缩试验破坏岩样

    Figure  5.   Damaged rock samples in conventional triaxial compression test

    图  6   不同围压下应力应变曲线

    Figure  6.   Stress-strain relationship curves under different confining pressures

    图  7   峰值强度与残余强度拟合曲线

    Figure  7.   Fitting curve of peak and residual intensities

    图  8   岩石峰后变形简化曲线

    Figure  8.   Simplified curve of post-peak deformation of rocks

    图  9   E400剖面数值模型

    Figure  9.   Numerical Model of cross-section diagram for E400 Section

    图  10   各级开挖变形图

    Figure  10.   Deformation diagram for each level of excavation

    图  11   南帮回填数值模型

    Figure  11.   Numerical model of cross-section diagram for Nanbang backfilling

    图  12   南帮回填位移图

    Figure  12.   Displacement diagram of Nanbang backfill

    表  1   弱层强度拟合参数

    Table  1   Fitting parameters for weak layer strength

    目标参数abcR2
    残余强度0.157360.814460.69310.99969
    峰值强度0.170030.347144.895170.99155
    下载: 导出CSV

    表  2   岩体力学参数

    Table  2   Summary table of rock mechanical parameters

    岩性密度
    /(kg·m−3
    弹性模量
    /GPa
    泊松比
    抗剪断强度
    内摩擦角/(°)黏聚力/MPa
    玄武岩280070.14420.300
    花岗片麻岩280080.2450.400
    弱层23000.10.416.90.075
    煤层13001.20.24350.140
    油母页岩23003.40.26350.200
    断层23000.10.416.90.075
    杂填土18000.10.4180.100
    凝灰岩260080.25390.250
    下载: 导出CSV

    表  3   节理力学参数

    Table  3   Summary table of joint mechanical parameters

    岩性刚度参数抗剪断强度
    法向刚度
    /(GPa·m−1
    切向刚度
    /(GPa·m−1
    节理面
    摩擦角/(°)
    节理面黏聚力
    /kPa
    玄武岩3.333.332590
    花岗片麻岩5.935.9345150
    弱层1.811.812010
    煤层2.432.432540
    油母页岩2.772.772540
    断层1.811.812010
    凝灰岩3.673.6728100
    下载: 导出CSV

    表  4   滑坡稳定状态划分

    Table  4   Division of landslide stable state

    滑坡稳定系数F<1.001.00≤F<1.051.05≤F<1.15F≥1.15
    滑坡稳定状态不稳定欠稳定基本稳定稳定
    下载: 导出CSV
  • [1] 邹宗兴,唐辉明,熊承仁,等. 大型顺层岩质滑坡渐进破坏地质力学模型与稳定性分析[J]. 岩石力学与工程学报,2012,31(11):2222 − 2231. [ZOU Zongxing,TANG Huiming,XIONG Chengren,et al. Geomechanical model of progressive failure for large consequent bedding rockslide and its stability analysis[J]. Chinese Journal of Rock Mechanics and Engineering,2012,31(11):2222 − 2231. (in Chinese with English abstract)

    ZOU Zongxing, TANG Huiming, XIONG Chengren, et al. Geomechanical model of progressive failure for large consequent bedding rockslide and its stability analysis[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(11): 2222-2231. (in Chinese with English abstract)

    [2] 张泽林,王涛,吴树仁,等. 泥岩中软弱夹层的剪切力学特性研究[J]. 岩石力学与工程学报,2021,40(4):713 − 724. [ZHANG Zelin,WANG Tao,WU Shuren,et al. Study on shear mechanical properties of mudstone with weak intercalation[J]. Chinese Journal of Rock Mechanics and Engineering,2021,40(4):713 − 724. (in Chinese with English abstract)

    ZHANG Zelin, WANG Tao, WU Shuren, et al. Study on shear mechanical properties of mudstone with weak intercalation[J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(4): 713-724. (in Chinese with English abstract)

    [3]

    SLOAN S W. Geotechnical stability analysis[J]. Géotechnique,2013,63(7):531 − 571.

    [4] 孔祥曌,李滨,贺凯,等. 柱状岩体崩塌动力特征与破碎规律—以重庆甑子岩崩塌为例[J]. 中国地质灾害与防治学报,2022,33(5):1 − 10. [KONG Xiangzhao, LI Bin, HE Kai, et al. Dynamic characteristics and fragmentation evolution of columnar rockfall: A case study of the Zengziyan rockfall in Chongqing, China[J]. The Chinese Journal of Geological Hazard and Control,2022,33(5):1 − 10. (in Chinese with English abstract)

    KONG Xiangzhao, LI Bin, HE Kai, et al. Dynamic characteristics and fragmentation evolution of columnar rockfall: a case study of the Zengziyan rockfall in Chongqing, China[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33(5)1-10(in Chinese with English abstract)

    [5] 张家勇,邹银先,杨大山. 基于PFC3D的鱼鳅坡滑坡运动过程分析[J]. 中国地质灾害与防治学报,2021,32(4):33 − 39. [ZHANG Jiayong,ZOU Yinxian,YANG Dashan. Analysis of Yuqiupo landslide motion process based on PFC3D[J]. The Chinese Journal of Geological Hazard and Control,2021,32(4):33 − 39. (in Chinese with English abstract)

    ZHANG Jiayong, ZOU Yinxian, YANG Dashan. Analysis of Yuqiupo landslide motion process based on PFC3D[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32(4)33-39(in Chinese with English abstract)

    [6] 成浩,韩培锋,苏有文. 含石量和坡度对土石混合体崩塌运动规律的影响[J]. 中国地质灾害与防治学报,2020,31(4):11 − 23. [CHENG Hao,HAN Peifeng,SU Youwen. Influence of rock content and slope on collapse in soil-rock accumulation[J]. The Chinese Journal of Geological Hazard and Control,2020,31(4):11 − 23. (in Chinese with English abstract)

    CHENG Hao, HAN Peifeng, SU Youwen. Influence of rock content and slope on collapse in soil-rock accumulation[J]. The Chinese Journal of Geological Hazard and Control, 2020, 31(4)11-23(in Chinese with English abstract)

    [7] 宋子岭,杨添,赵立春. 含多层软弱夹层的顺向岩质边坡稳定性评价方法对比分析[J]. 中国地质灾害与防治学报,2016,27(2):20 − 25. [SONG Ziling,YANG Tian,ZHAO Lichun. Applicatyion of combination algorithm in stability analysis of rock slope with contains multilayer weak bedrock layers[J]. The Chinese Journal of Geological Hazard and Control,2016,27(2):20 − 25. (in Chinese with English abstract)

    SONG Ziling, YANG Tian, ZHAO Lichun. Applicatyion of combination algorithm in stability analysis of rock slope with contains multilayer weak bedrock layers[J]. The Chinese Journal of Geological Hazard and Control, 2016, 27(2): 20-25. (in Chinese with English abstract)

    [8]

    XUE Demin,LI Tianbin,ZHANG Shuai,et al. Failure mechanism and stabilization of a basalt rock slide with weak layers[J]. Engineering Geology,2018,233:213 − 224. DOI: 10.1016/j.enggeo.2017.12.005

    [9] 张社荣,谭尧升,王超,等. 多层软弱夹层边坡岩体破坏机制与稳定性研究[J]. 岩土力学,2014,35(6):1695 − 1702. [ZHANG Sherong,TAN Yaosheng,WANG Chao,et al. Research on deformation failure mechanism and stability of slope rock mass containing multi-weak interlayers[J]. Rock and Soil Mechanics,2014,35(6):1695 − 1702. (in Chinese with English abstract)

    ZHANG Sherong, TAN Yaosheng, WANG Chao, et al. Research on deformation failure mechanism and stability of slope rock mass containing multi-weak interlayers[J]. Rock and Soil Mechanics, 2014, 35(6): 1695-1702. (in Chinese with English abstract)

    [10]

    LI Jinglong,ZHANG Bo,SUI Bin. Stability analysis of rock slope with multilayer weak interlayer[J]. Advances in Civil Engineering,2021,2021:1 − 9.

    [11]

    TAN Xin,REN Yakun,LI Tenglong,et al. In-situ direct shear test and numerical simulation of slate structural planes with thick muddy interlayer along bedding slope[J]. International Journal of Rock Mechanics and Mining Sciences,2021,143:104791. DOI: 10.1016/j.ijrmms.2021.104791

    [12]

    ZHENG Qingsong,LIU Enlong,YU Di,et al. Fatigue and damage properties of non-consecutive jointed mudstone samples subjected to cyclic triaxial loading[J]. Bulletin of Engineering Geology and the Environment,2020,79(5):2467 − 2481. DOI: 10.1007/s10064-019-01693-x

    [13]

    YANG Shengqi,TIAN Wenling,JING Hongwen,et al. Deformation and damage failure behavior of mudstone specimens under single-stage and multi-stage triaxial compression[J]. Rock Mechanics and Rock Engineering,2019,52(3):673 − 689. DOI: 10.1007/s00603-018-1622-y

    [14]

    YANG Bing,HOU Jiangrong,LIU Yifei,et al. Dynamic response and failure characteristics of slope with weak interlayer under action of near-fault ground motion[J]. Shock and Vibration,2021,2021:1 − 18.

    [15] 王来贵,向丽,赵娜,等. 地震作用下顺倾多弱层岩质边坡动力响应[J]. 中国地质灾害与防治学报,2021,32(6):18 − 25. [WANG Laigui,XIANG Li,ZHAO Na,et al. Dynamic response of down-dip multi-weak-layer rock slope under earthquake[J]. The Chinese Journal of Geological Hazard and Control,2021,32(6):18 − 25. (in Chinese with English abstract)

    WANG Laigui, XIANG Li, ZHAO Na, et al. Dynamic response of down-dip multi-weak-layer rock slope under earthquake[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32(6): 18-25. (in Chinese with English abstract)

    [16] 李鹏, 苏生瑞, 王闫超, 等. 含软弱层岩质边坡的动力响应研究[J]. 岩土力学, 2013, 34(增刊1): 365-370

    LI Peng, SU Shengrui, WANG Yanchao, et al. Research on dynamic response of rock slope with weak layer[J]. Rock and Soil Mechanics, 2013, 34(Sup 1): 365-370. (in Chinese with English abstract)

    [17] 皮晓清,李亮,唐高朋,等. 基于有限元极限上限法的含软弱夹层边坡稳定性分析[J]. 铁道科学与工程学报,2019,16(2):351 − 358. [PI Xiaoqing,LI Liang,TANG Gaopeng,et al. Stability analysis for soil slopes with weak interlayers using the finite element upper bound limit analysis[J]. Journal of Railway Science and Engineering,2019,16(2):351 − 358. (in Chinese with English abstract)

    PI Xiaoqing, LI Liang, TANG Gaopeng, et al. Stability analysis for soil slopes with weak interlayers using the finite element upper bound limit analysis[J]. Journal of Railway Science and Engineering, 2019, 16(2): 351-358. (in Chinese with English abstract)

    [18] 杨令强,练继建,陈祖坪. 软弱夹层的接触问题与强度问题探讨[J]. 水利水电技术,2002,33(6):17 − 19. [YANG Lingqiang,LIAN Jijian,CHEN Zuping. Discussion on contact problem and strength problem of weak interlayer[J]. Water Resources and Hydropower Engineering,2002,33(6):17 − 19. (in Chinese)

    YANG Lingqiang, LIAN Jijian, CHEN Zuping. Discussion on contact problem and strength problem of weak interlayer[J]. Water Resources and Hydropower Engineering, 2002, 33(6): 17-19. (in Chinese)

    [19] 吕布, 杨之俊, 韦秀东, 等. 坝基不同倾角软弱夹层的破坏模式及D-C模型参数研究[J]. 岩土工程学报, 2019, 41(增刊1): 189 − 192

    LU Bu, YANG Zhijun, WEI Xiudong, et al. Failure modes and constitutive model for weak interlayer of dam foundation with different inclination angles[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(Sup 1): 189 − 192. (in Chinese with English abstract)

    [20] 郑青松,刘恩龙,刘明星. 三轴试验下结构面倾角对制备岩样力学特性的影响[J]. 岩土力学,2019,40(5):1854 − 1861. [ZHENG Qingsong,LIU Enlong,LIU Mingxing. Influence of dip angle of structural planes on mechanical properties of artificial rock samples under triaxial test conditions[J]. Rock and Soil Mechanics,2019,40(5):1854 − 1861. (in Chinese with English abstract)

    ZHENG Qingsong, LIU Enlong, LIU Mingxing. Influence of dip angle of structural planes on mechanical properties of artificial rock samples under triaxial test conditions[J]. Rock and Soil Mechanics, 2019, 40(5): 1854-1861. (in Chinese with English abstract)

    [21] 李东. 抚顺西露天矿南帮高陡弱层边坡变形规律研究[D]. 阜新: 辽宁工程技术大学, 2021

    LI Dong. Study on deformation law of high, steep and weak layer slope in the south side of Fushun west open pit mine[D]. Fuxin: Liaoning Technical University, 2021. (in Chinese with English abstract)

    [22] 刘传正,崔原,陈春利,等. 辽宁抚顺西露天矿南帮滑坡成因[J]. 地质通报,2022,41(5):713 − 726. [LIU Chuanzheng,CUI Yuan,CHEN Chunli,et al. Research on the south side landslide at west open-pit coal mine in Fushun City, Liaoning Province of China[J]. Geological Bulletin of China,2022,41(5):713 − 726. (in Chinese with English abstract)

    LIU Chuanzheng, CUI Yuan, CHEN Chunli, et al. Research on the south side landslide at west open-pit coal mine in Fushun City, Liaoning Province of China[J]. Geological Bulletin of China, 2022, 41(5)713-726(in Chinese with English abstract)

    [23] 孙闯, 张向东, 张涛, 等. 深部大跨度泥质顶板剪切冒落失稳区预测研究[J]. 中国安全生产科学技术, 2016, 12(1): 23 − 27

    SUN Chuang, ZHANG Xiangdong, ZHANG Tao, et al. Prediction on shear caving instability region of deep and large span muddy roof[J]. Journal of Safety Science and Technplogy, 2016, 12(1): 23 − 27. (in Chinese with English abstract)

    [24] 孙闯,张树光,贾宝新,等. 花岗岩峰后力学特性试验与模型研究[J]. 岩土工程学报,2015,37(5):847 − 852. [SUN Chuang,ZHANG Shuguang,JIA Baoxin,et al. Physical and numerical model tests on post-peak mechanical properties of granite[J]. Chinese Journal of Geotechnical Engineering,2015,37(5):847 − 852. (in Chinese with English abstract)

    SUN Chuang, ZHANG Shuguang, JIA Baoxin, et al. Physical and numerical model tests on post-peak mechanical properties of granite[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(5): 847-852. (in Chinese with English abstract)

    [25] 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会. 滑坡防治工程勘查规范GB/T 32864-2016[S], 2016.

    General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China, Standardization Administration of China.Specification of geological investigation for landslide stabilization GB/T 32864-2016[S], 2016. (in Chinese)

    [26] 汪美华,李勇,裴叶青. 甘肃临夏盆地韩集北山滑坡群致灾特征与稳定性评价[J]. 地质通报,2023,42(2/3):460 − 468. [WANG Meihua, LI Yong, PEI Yeqing. Disaster characteristics and stability evaluation of the Hanjin Beishan landslide group in Linxia Basin, Gansu Province[J]. Geological Bulletin of China,2023,42(2/3):460 − 468. (in Chinese with English abstract)

    [WANG Meihua, LI Yong, PEI Yeqing. Disaster characteristics and stability evaluation of the Hanjin Beishan landslide group in Linxia Basin, Gansu Province[J]. Geological Bulletin of China, 2023, 42(Sup 1): 460-468.(in Chinese with English abstract)

  • 期刊类型引用(11)

    1. 杨磊,张耀磊,宋光啸,覃剑文,潘光明. 覆盖型岩溶区地面塌陷灾变规律与影响因素. 应用基础与工程科学学报. 2025(01): 192-204 . 百度学术
    2. 刘动,昝子卉,杨贝贝,黄文彬. 岩溶塌陷精细化探查与治理关键技术. 广东土木与建筑. 2024(08): 53-58 . 百度学术
    3. 姜明映,陈娥梅,刘永成,朱亚东,蒋文奎,谢瑞庭. 岩溶区钻孔桩泥浆漏失引起土体变形机制与破坏模式试验研究. 中外公路. 2024(05): 59-64 . 百度学术
    4. 张晓龙. 宝浪苏木分水枢纽闸门建设地质分析. 海河水利. 2023(02): 100-103 . 百度学术
    5. 任涛,田国亮,宁志杰,周爱红,李宽,陈石. 基于地理探测器和随机森林的岩溶塌陷易发性评价. 灾害学. 2023(03): 227-234 . 百度学术
    6. 孙伟. 广东省英德市城南社区岩溶塌陷发育特征及成因分析. 中国地质灾害与防治学报. 2023(05): 74-80 . 本站查看
    7. 张赞萍. 基于综合物探方法的岩溶塌陷地质灾害勘查及地质材料属性研究. 现代应用物理. 2023(04): 207-214 . 百度学术
    8. 贾龙,雷明堂,程小杰. 基于井中超声波成像的岩溶特征高精度探测和评价. 地质通报. 2022(Z1): 453-460 . 百度学术
    9. 蔡洋洋. 高密度电法在岩溶塌陷勘察中的应用研究. 世界有色金属. 2022(08): 184-187 . 百度学术
    10. 赵建军. 广东英德沙口镇岩溶地面塌陷发育特征及形成机理分析. 中国地质灾害与防治学报. 2022(05): 59-65 . 本站查看
    11. 韩晓鹏. 岩溶地区地下工程溶洞处理分析设计. 城市建筑. 2021(30): 134-136+139 . 百度学术

    其他类型引用(4)

图(12)  /  表(4)
计量
  • 文章访问数:  2006
  • HTML全文浏览量:  613
  • PDF下载量:  680
  • 被引次数: 15
出版历程
  • 收稿日期:  2022-05-12
  • 修回日期:  2022-07-22
  • 网络出版日期:  2023-04-25
  • 刊出日期:  2023-06-24

目录

/

返回文章
返回