摘要
软弱夹层作为一种典型地质结构是影响整个边坡安全稳定性的关键因素。为了更好地分析软弱夹层的性质及对边坡稳定性的影响,针对某含软弱夹层矿山露天岩质边坡,运用正交对照法对边坡软弱夹层的埋深、倾角、数量及夹层间结构面间距等4种夹层赋存状态,进行变量4×4正交模拟,揭示了不同赋存状态下边坡稳定性及变形规律。研究结果表明,1)边坡软弱夹层的赋存状态对边坡稳定性有显著影响,随软弱夹层埋藏深度的增加边坡安全系数随之增大,当埋深达到一定值后,边坡安全系数趋于稳定;2)随着软弱夹层角度的增大,破坏机制表现为层间错动、顺层蠕滑和沿软弱结构面的剪切滑移递进式变化;3)随软弱夹层层数增多,边坡整体性下降致使变形不均匀,边坡安全系数减小,水平方向位移增大;4)对比同等埋深结构面间距较小的软弱夹层边坡,结构面间距对边坡安全系数的影响较小。
矿山开采过程中常遇到边坡含软弱夹层的问题,这种岩质边坡的存在对矿山的稳定性和变形产生了重大影响。软弱夹层具有较低的强度和稳定性,容易出现滑动、断裂和塌陷等不稳定现象。因此,矿山边坡的稳定性和变形特征与软弱夹层的状态密切相
目前,相关学者关于软弱夹层对边坡稳定性及变形的研究已取得一定的进展。例如,王来贵
基于此,由于地质构造环境的不同,软弱夹层岩质边坡影响因素各异,主要有以下几类:1)岩体与软弱夹层的物理力学性质;2)软弱夹层几何参数(埋深、倾角、夹层数量、夹层间结构面间距等);3)初始地应力、水文气象、地震和人工作用等。但是关于露天矿山岩质边坡软弱夹层赋存状态对边坡稳定性及变形的影响的研究深度与系统性还不够,因此,以云南某露天矿山含软弱夹层岩质开采边坡为基础,研究矿山岩质边坡软弱夹层赋存状态对边坡稳定性及变形的影响。基于Bishop法和有限元强度折减法,探讨岩质边坡在不同软弱夹层埋深、倾角、夹层数量、夹层间结构面间距条件下的安全系数与边坡位移的变化规律,从而揭示含软弱夹层岩质边坡的变形及稳定性特征。力求通过分析边坡的稳定性和变形特征,为矿山设计及施工提供科学依据,减少矿山灾害的发
云南某大型磷资源露天开采矿山,边坡岩体结构类型为层状结构,层间有水云母黏土岩软弱夹层发育。其边坡岩体主要为白云岩及砂质白云岩,含有2层水云母黏土岩软弱夹层。此类软弱夹层含泥质、粉砂质物质,因此岩体抗风化能力差,层面黏结性差。岩体常规物理力学参数如
岩性 | 容重/(kg∙ | 弹性模量/GPa | 内摩擦角/(°) | 黏聚力c/MPa | 泊松比 |
---|---|---|---|---|---|
白云岩及砂质白云岩 | 27.6 | 27.997 | 39.2 | 0.180 | 0.31 |
软弱夹层(水云母黏土岩) | 17.3 | 2.230 | 18.7 | 0.059 | 0.32 |
该典型矿山含软弱夹层岩质边坡呈东-西走向。该矿山边坡上部已经开挖成了5个台阶,分别是2190平台、2160平台、2130平台、2100平台和2070平台。边坡角度约为45°,边坡的最大高差为160 m。剖面位置选取边坡整体最大高程处,剖面包含软弱夹层及边坡所有平台且平台较为清晰,因此选取此剖面具有较好的代表性。典型剖面选取如

图1 采场边坡模型图
Fig. 1 3D topographic model of stope slope
剖面模型长230 m,选取高程2 030~2 190 m,总高度H=160 m,厚度25 m,软弱夹层厚度为h=3.61 m,夹层夹角为15°。采用犀牛软件建模与网格划分(网格边长1~5 m),共26 509网格38 797节点,并将边坡模型导入到FLAC 3D中得到,剖面模型图如2所示。边坡在自然工况下设置边界条件为左右两侧施加法向约束,上部为自由边界,下部为底板固定约束。

图2 边坡剖面计算模型图
Fig. 2 Slope profile calculation model diagram

图3 Bishop法力学模型图
Fig. 3 Bishop method mechanics calculation model diagram
Bishop法是一种极限平衡条分法,将圆弧滑面垂直方向分成n条,该法忽略竖向条间剪力,根据整体力矩平衡条件求解安全系数。由于该方法满足土条间的水平推力平衡条件,其精度要高于瑞典圆弧法。Bishop法力学模型如
该法的计算参数:Fs为滑面安全系数、R为圆弧滑面半径、bi为第i条块宽度、Wi为第i条块重力、αi为第i条块圆弧底面倾角、Xi和Xi+1为第i与i+1条块间横向力、Ei和Ei+1为第i与i+1条块间纵向力、Ti为第i条块圆弧底面剪力、Ni为第i条块法向作用力、ui为第i条块孔隙水压力。其计算公式为:
, | (1) |
。 | (2) |
式中:Ei及Xi表示法向及切向条间力;Wi为i条块自重;Qi为水平力;ci为材料黏结力和i为材料内摩擦角;ei为第i条块重心到滑弧圆心的垂直距离。
传统强度折减法是一种常用于边坡稳定性分析的方法,其原理基于材料的强度折减。在应用该方法时,首先建立边坡的力学模型,再根据不同的受力情况(如坡面滑动、坡脚滑移等),计算出边坡在不同情况下的安全系数。然后,按照一定的强度折减系数,将边坡的强度参数进行折减,重新计算安全系数,直到得到符合要求的安全系数,即为边坡安全系
强度折减系数的计算公式如下,
, | (3) |
式中:k为强度折减系数;m为强度折减的比例;n为强度折减的次数。在边坡稳定性分析中,通常将m取为较小的值(如0.1),n取为较大的值(如10),以保证计算结果的精度。
强度折减后的边坡强度参数可以用以下公式计算:
, | (4) |
。 | (5) |
式中:C'和φ'分别为强度折减后的剪切强度和内摩擦角;C和φ为原始的剪切强度和内摩擦角。在计算边坡安全系数时,需要使用强度折减后的强度参数进行计算。
根据对矿山现场勘察结果及其勘察地形图,在采场边坡带有软弱夹层的凌空面选取代表性剖面作为研究对象,如

图4 软弱夹层不同赋存状态方案图
Fig. 4 Scheme diagram of different occurrence states of soft sandwich
方案 | 埋藏深度 | 倾角/(°) | 夹层数量/条 | 夹层间结构面间距 |
---|---|---|---|---|
方案一(埋深) | 0.1H~0.7H | 15 | 2 | 2h |
方案二(倾角) | 0.5H | 5~35 | 2 | 2h |
方案三(数量) | 0.5H | 15 | 2~5 | 2h |
方案四(厚度) | 0.5H | 15 | 2 | 2h~8h |
基于Bishop法与强度折减法分析矿山岩质边坡软弱夹层不同埋深对边坡稳定性的影响。采取固定15°倾角,软弱夹层为2层,夹层间结构距离2h=7.22 m,选取截面右侧上层软弱夹层埋深作为研究基准点,分别取距离上凌空面0.1H=16 m、0.3H=48 m、0.5H=80 m、0.7H=112 m(H=160 m为选取截面的总高程)作为研究对象,如

图5 软弱夹层不同埋藏深度对边坡安全系数及最大水平位移变化趋势图
Fig. 5 Variation trend of slope safety factor and maximum horizontal displacement of soft interlayer with different burial depths
由
结合

图6 软弱夹层不同埋藏深度下边坡x方向位移云图
Fig. 6 x-direction displacement cloud of slope under different buried depths of soft interlayer
当软弱夹层埋藏深度为0.5H=80 m时,软弱夹层对边坡水平方向位移影响较小,位移主要受基岩位移的影响;当软弱夹层埋藏深度为0.7H=112 m时,软弱夹层埋深已经低于边坡整体滑移带,对边坡水平位移基本不产生影响。
在本节中,以
经计算,在5°~35°倾角下,软弱夹层边坡安全系数与最大水平位移如

图7 不同软弱夹层倾角下边坡安全系数及最大水平位移变化趋势图
Fig. 7 Trend chart of slope safety factor and maximum horizontal displacement under different inclination angles of soft sandwich layers
呈现负相关,随着夹层角度增加,安全系数值减小。当倾角大于15°时,边坡安全系数呈加快下降趋势且角度越大安全系数减小越快。由于2种方法计算原理不同,故计算所得安全系数具有一定的差别,当软弱夹层倾角在5°和35°时差距较大。Bishop法所求安全系数大于同倾角下强度折减法计算结果,这主要是由于Bishop法基于极限平衡原理,假定边坡为剪切破坏。而强度折减法在抗剪强度一定的情况下,上部岩体易出现拉裂破坏,因为岩体的抗拉强度较小,强度折减法所得的安全系数相对较
不同软弱夹层倾角条件下边坡的x方向位移与应变区域位置差异较大,如

图8 软弱夹层不同倾角下边坡x方向位移云图
Fig. 8 Slope x-direction displacement cloud under different inclination angles of soft sandwich
对含有软弱夹层的边坡,边坡中含有软弱夹层的数量对其稳定性有较大的影响,以软弱夹层层数2~5层为例,如
控制软弱夹层角度及层间厚度研究不同软弱夹层数量对矿山开采边坡稳定性的影响,如

图9 不同软弱夹层数量下边坡安全系数及最大水平位移变化趋势图
Fig. 9 Trend chart of slope safety factor and maximum horizontal displacement under different number of soft interlayers
由

图10 不同软弱夹层数量下边坡x方向位移云图
Fig. 10 Slope x-direction displacement cloud under different number of soft interlayers
对含软弱夹层的边坡,不同结构面间距条件下边坡的安全系数存在一定差距。以边坡软弱夹层2层,倾角15°为例,模型图如


图11 不同结构面间距下边坡x方向位移云图
Fig. 11 Slope x-displacement cloud under different structural surface spacing

图12 不同结构面间距下边坡安全系数及最大水平位移变化趋势
Fig. 12 The variation trend of slope safety factor and maximum horizontal displacement under different structural surface spacing
由图
1)软弱夹层埋深对含软弱夹层岩质边坡稳定性影响显著,随着软弱夹层埋藏深度的增加,边坡安全系数增加的同时水平位移减小,在埋深较浅时,软弱夹层对边坡稳定性影响较为显著;随埋深的增大,软弱夹层对边坡稳定性影响逐渐减小,当埋深达到一定深度后,边坡安全系数增加缓慢,水平方向位移变化也随之下降。当软弱夹层埋藏深度大于0.5H时,在此深度下,边坡稳定性几乎不受软弱夹层控制。
2)软弱夹层倾角在5°~35°范围内时,随软弱夹层倾角的增加,边坡岩体破坏大致经历层间错动、顺层蠕滑和沿软弱结构面剪切滑移几个阶段;且边坡滑裂带逐渐变浅,较大应变区域面积增大且集中于软弱夹层上部,滑出口逐渐向坡址推进,边坡稳定性下降。
3)软弱夹层层数对边坡水平方向位移影响较大,随软弱夹层层数增多,边坡安全系数减小,水平方向位移增大;且边坡整体刚度下降使得边坡在受到外力作用时更容易发生水平方向的位移,多层软弱夹层的存在会导致边坡的整体变形不均匀;边坡潜在滑移面主要集中在前缘,受最上层软弱夹层影响较大。
4)对比同等埋深较小结构面间距软弱夹层边坡,结构面间距对边坡安全系数的影响较小。在边坡软弱夹层结构面间距不同的情况下,边坡变形破坏机制基本一致,边坡稳定性主要受最上层软弱夹层埋深位置影响。下层软弱夹层对基岩的扰动导致边坡水平位移增加,结构面间距的增大会导致边坡整体的刚度降低从而导致安全系数减小。
参考文献
张家明. 含软弱夹层岩质边坡稳定性研究现状及发展趋势[J]. 工程地质学报, 2020, 28(3): 626-638. [百度学术]
Zhang J M. State of art and trends of rock slope stability with soft interlayer[J]. Journal of Engineering Geology, 2020, 28(3): 626-638. (in Chinese) [百度学术]
王来贵, 马玮航, 赵娜, 等. 含软弱夹层顺层岩质边坡变形-滑动演化过程研究[J]. 辽宁工程技术大学学报(自然科学版), 2023, 42(2): 129-135. [百度学术]
Wang L G, Ma W H, Zhao N, et al. Deformation-sliding evolution process of bedding rock slope with weak interlayer[J]. Journal of Liaoning Technical University (Natural Science), 2023, 42(2): 129-135.(in Chinese) [百度学术]
宋洋, 马旭琪, 赵常青, 等. 隐伏非贯通软弱夹层岩质边坡剪切蠕变特征及稳定性研究[J]. 岩土工程学报, 2024, 46(4):755-763. [百度学术]
Song Y, Ma X Q,Zhao C Q, et al. Shear creep characteristics and stability of rock slopes with concealed noninterpenetrated weak interlayer[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(4):755-763.(in Chinese) [百度学术]
陈权川, 王安礼, 邬忠虎, 等. 基于尖点突变理论的含软弱夹层缓倾顺层岩质边坡稳定性研究[J]. 河南理工大学学报(自然科学版), 2022, 41(6): 188-195. [百度学术]
Chen Q C, Wang A L, Wu Z H, et al. Study on stability of gently inclined bedding rock slope with weak interlayer based on cusp catastrophe theory[J]. Journal of Henan Polytechnic University (Natural Science), 2022, 41(6): 188-195.(in Chinese) [百度学术]
Tien Y M, Kuo M C. A failure criterion for transversely isotropic rocks[J]. International Journal of Rock Mechanics and Mining Sciences, 2001, 38(3):399-412. [百度学术]
刘杨, 彭亮, 祝鑫, 等. 降雨条件下含缓倾软弱夹层矿山高边坡稳定性分析[J]. 矿业研究与开发, 2022, 42(7): 71-79. [百度学术]
Liu Y, Peng L, Zhu X, et al. Stability analysis of high slope with gently inclined weak interlayer under rainfall condition in mine[J]. Mining Research and Development, 2022, 42(7): 71-79.(in Chinese) [百度学术]
Zheng Q S, Liu E L, Yu D, 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. [百度学术]
范昊, 肖诗荣, 魏瑞琦, 等. 软弱夹层对某矿山边坡稳定性的影响分析[J]. 三峡大学学报(自然科学版), 2022, 44(1): 54-59. [百度学术]
Fan H, Xiao S R, Wei R Q, et al. Analysis of influence of weak interlayer on slope stability of a mine[J]. Journal of China Three Gorges University (Natural Sciences), 2022, 44(1): 54-59.(in Chinese) [百度学术]
范刚, 张建经, 付晓, 等. 含软弱夹层顺层岩质边坡动力破坏模式的能量判识方法研究[J]. 岩土工程学报, 2016, 38(5): 959-966. [百度学术]
Fan G, Zhang J J, Fu X, et al. Energy identification method for dynamic failure mode of bedding rock slope with soft strata[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(5): 959-966.(in Chinese) [百度学术]
Chen Z L, Hu X, Xu Q. Experimental study of motion characteristics of rock slopes with weak intercalation under seismic excitation[J]. Journal of Mountain Science, 2016, 13(3): 546-556. [百度学术]
Wang K L, Lin M L. Initiation and displacement of landslide induced by earthquake: a study of shaking table model slope test[J]. Engineering Geology, 2011, 122(1/2): 106-114. [百度学术]
Fan G, Zhang J J, Wu J B, et al. Dynamic response and dynamic failure mode of a weak intercalated rock slope using a shaking table[J]. Rock Mechanics and Rock Engineering, 2016, 49(8): 3243-3256. [百度学术]
Zhu H H, Wang Z Y, Shi B, et al. Feasibility study of strain based stability evaluation of locally loaded slopes: insights from physical and numerical modeling[J]. Engineering Geology, 2016, 208: 39-50. [百度学术]
张社荣, 谭尧升, 王超, 等. 多层软弱夹层边坡岩体破坏机制与稳定性研究[J]. 岩土力学, 2014, 35(6): 1695-1702. [百度学术]
Zhang S R, Tan Y S, Wang C, 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) [百度学术]
尤耿明, 王光进, 刘文连, 等. 矿山边坡软弱夹层赋存状态及蠕变特性对边坡稳定性影响研究[J]. 金属矿山, 2021(12): 200-206. [百度学术]
You G M, Wang G J, Liu W L, et al. Study on the influence of occurrence state and creep characteristics of weak interlayer on mine slope stability[J]. Metal Mine, 2021(12): 200-206.(in Chinese) [百度学术]
Zhu S N, Yin Y P, Li B, et al. Shear creep characteristics of weak carbonaceous shale in thick layered Permian limestone, southwestern China[J]. Journal of Earth System Science, 2019, 128(2): 28. [百度学术]
方庆红. 含缓倾软弱夹层露天高边坡坡态控制参数优化及稳定性研究[D]. 武汉: 武汉科技大学, 2021. [百度学术]
Fang Q H. Optimization of slope state control parameters and stability study of open-pit high slope with gently inclined weak interlayer[D].Wuhan: Wuhan University of Science and Technology, 2021. (in Chinese) [百度学术]
郑川. 软弱夹层条件下岩质边坡稳定性分析与治理方案研究[D]. 重庆: 重庆交通大学, 2020. [百度学术]
Zheng C. Stability analysis and treatment scheme of rock slope under the condition of weak interlayer[D].Chongqing: Chongqing Jiaotong University, 2020. (in Chinese) [百度学术]
戴加祺. 含软弱夹层的岩质边坡稳定性影响因素的定性分析[J]. 西部探矿工程, 2017, 29(12): 13-15, 19. [百度学术]
Dai J Q. Qualitative analysis of influencing factors on stability of rock slope with weak interlayer[J]. West-China Exploration Engineering, 2017, 29(12): 13-15, 19.(in Chinese) [百度学术]
穆成林. 顺层岩质高边坡开挖过程变形失稳演化机制及预测评价研究: 以织金石化场区边坡为例[D]. 成都: 成都理工大学, 2017. [百度学术]
Mu C L. Study on deformation and instability evolution mechanism and prediction evaluation of bedding rock high slope during excavation: a case study of Zhijin petrochemical field slope[D].Chengdu: Chengdu University of Technology, 2017. (in Chinese) [百度学术]