东北大学 资源与土木工程学院, 辽宁 沈阳 110819
收稿日期:2017-01-14
基金项目:国家自然科学基金资助项目(51474050, U1602232);地质灾害防治与地质环境保护国家重点实验室项目(SKLGP2014K011);辽宁省高等学校优秀人才支持计划项目(LN2014006)。
作者简介:张紫杉(1990-),男,河北秦皇岛人,东北大学博士研究生;
王述红(1969-),男,江苏泰州人,东北大学教授,博士生导师。
摘要:将关键块体理论与模糊理论相结合, 对于岩质边坡的稳定性进行了综合评价.选取抚顺露天矿南帮典型坡段, 对岩坡的物理力学参数进行测定, 并运用自主研发的关键块体搜索系统——GeoSMA-3D (Geotechnical Structure and Model Analysis-3D)进行关键块体的搜索与表征, 根据表征特点进行评价参数的选择, 最后运用基于AHP的模糊综合评价技术, 引入总体稳定性系数, 对边坡稳定性进行了客观定量的评价.最后通过VTK技术实现了边坡稳定性可视化平台的构建, 为后续的边坡加固及灾害防治提供了可靠依据.
关键词:露天矿岩质边坡层次分析法块体理论模糊评价
Comprehensive Assessment of Rock Slope Stability Based on Spatial Block Identification
ZHANG Zi-shan1, WANG Shu-hong1, WANG Fei-li1
School of Resources & Civil Engineering, Northeastern University, Shenyang 110819, China
Corresponding author: WANG Shu-hong, E-mail: shwang@mail.neu.edu.cn
Abstract: A comprehensive assessment on rock slope stability was carried out by combining the key block theory with the fuzzy method. Based on typical slope of strip mine in Fushun, the physical and mechanical parameters were measured and a block identification system—GeoSMA-3D (Geotechnical Structure and Model Analysis-3D) was developed to identify the key blocks and to calculate the factors of safety. The assessment factors were detemined based on the characteristics of the spatial blocks. According to the information acquired, a fuzzy comprehensive assessment method based on AHP was applied to assess the stability of the slopes in an objective and quantitative way. Finally, an visualizaion plantform is established to visulize the terrain and total safety factors of the slope, which provides a reliable basis for disaster control and reinforcement on the slopes.
Key Words: strip minerock slopeanalytic hierarchy processblock theoryfuzzy evaluation
近些年, 关于边坡稳定性的研究越来越引起学者们的关注.边坡稳定性就是指在多元因素影响下边坡发生滑落、倾倒等破坏的可能性.目前采用基于大量分析指标的多元多级综合评价方法, 如层次分析法、灰色聚类法、神经网络法等, 都被运用到了边坡稳定性的计算之中, 并且取得了较好的评价效果[1-3].由于涉及的地形复杂, 影响因素多, 层次关系复杂, 单纯运用直接从现场采集的数据进行评价, 容易造成评价结构混乱、主次不清晰等问题, 给后期综合分析带来一定困难.
本研究根据已有的研究成果, 将块体理论与模糊综合评价方法进行结合, 提出了一种全新的基于块体理论的岩质边坡稳定性的评价方法.该方法弥补了传统块体理论方法只考虑块体可动性而没有考虑渗流、风化等水文地质条件及岩坡坡面形态对岩体稳定性影响的缺点, 使得评价可以更加客观准确.
1 研究方法1.1 岩坡现场测量技术为了获取实时有效的现场数据, 本研究运用了先进的非接触测量技术, 为边坡后期稳定性的分析计算提供了有力的技术支持.为了获得岩体几何参数和结构面参数, 运用了来自奥地利的三维非接触测量系统——ShapeMetrix3D, 该系统可以通过高清摄影图片对岩体及其表面形状进行三维建模, 并且可以得到岩体的空间形态和表面节理的相关信息, 包括岩体的外部轮廓参数和岩体结构面的倾向、倾角、中心点坐标等[4].除此之外, 本研究还引进了红外线成像技术, 来对边坡进行分幅扫描观测, 获得边坡的温度场分布情况, 根据温度场变化情况提取出可能发生渗水的地点.
1.2 块体信息处理块体理论是由美籍华人科学家石根华提出的, 该方法运用拓扑学的原理, 解决了岩体稳定性的问题, 目前已经广泛应用到包括三峡工程在内的大型工程之中[5-6].该理论假设节理岩体是由平面切割整块完整岩石后形成的结构面和结构体组合成的.运用几何方法, 通过判断块体空间可动性及有限性, 利用赤平投影法, 寻找出关键块体.由于岩体本身强度比岩石结构面强度要大得多, 在考虑岩石破坏时, 只考虑结构面的剪切破坏和滑移破坏.通过式(1), 计算各个块体在滑动力作用下的安全系数K.
(1) |
为了更加准确搜索到关键块体的相关信息, 本文运用GeoSMA-3D (Geotechnical Structure and Model Analysis-3D)系统进行分析, 确定关键块体的体积、数量及相应的安全系数[7-8].同时, 为了对块体系统的稳定性进行评价, 对搜索出的关键块体信息需要进行进一步的处理.由于边坡的体积并不相同, 为了得到统一的评价标准, 要对搜索的块体体积进行归一化处理, 即计算块体体积占总体边坡的比例v, 以及按块体体积加权的方法得出边坡总体安全系数Fs总:
(2) |
(3) |
1.3 AHP评价理论层次分析模型又称AHP, 是一种基于权重分析的多层次的综合评价模型.该方法同时具备了定性分析及定量分析的优点, 获得了广泛的应用[9].层次分析法的主要思路是把复杂的因素分类汇总, 形成有一定规律的结构层次体系, 再通过比对同一层次中的各要素的重要性, 给出合理的权重, 并在通过了一致性检验后, 得到最终的权重体系.层次分析法的过程主要分为以下3步:
1) 建立递进式层次结构.按照分类分级的思想, 对研究指标进行分析, 构建包含目标层、指标层、准则层的递进式的层次分析结构.根据之前的研究成果, 除空间块体搜索表征的参数外, 选出了两类6项主要影响岩坡稳定性的指标, 即地形特征(坡面形状、坡角大小、边坡高度)和地质特征(岩体强度、岩体渗水性、岩体风化程度), 构建边坡稳定性评价层次分析体系, 如图 1所示.
图 1(Fig. 1)
图 1 抚顺露天矿边坡综合评价体系Fig.1 Analytical hierarchy diagram of rock slope in Fushun strip mine |
2) 建立判断矩阵, 并进行一致性检验.建立了递进式层次分析评价结构后, 针对隶属关系, 对重要性标度进行量化, 建立判断矩阵W, 见式(4).其中wij表示同一层中第i指标对第j指标的重要性标度值.并采用1~9的标度对各层中因子相对于上一目标层的重要性进行比较[9].标准的量化过程既考虑工程经验发挥作用, 又可以有效减少主观臆断, 使评价更加合理.
(4) |
(5) |
(6) |
(7) |
(8) |
1) 确定评价因素集及评价等级集.评价因素集是指由被评价对象组成的集合, 设被评价对象有m个, 分别用u1, u2, …, um表示;评价等级集表示将最终的评价结果划分成几个不同的等级, 设分为n个等级, 分别用v1, v2, …, vn表示:
(9) |
(10) |
(11) |
(12) |
(13) |
(14) |
(15) |
(16) |
(17) |
(18) |
为了更加合理地作出评价, 本文选取了具有典型特征的坡段区域进行分析评价, 并且运用GeoSMA-3D进行关键块体的搜索, 如图 2所示.最后结合测量内容以及已有的科研成果, 建立相应的评价体系和评价指标(表 1), 得到评价需要的参数(表 2).
图 2(Fig. 2)
图 2 +20E600坡段的关键块体搜索Fig.2 Key block identification on+20E600 slope (a)—ShapeMetix-3D结构面采集;(b)—岩体表面切割迹线;(c)—关键块体搜索. |
表 1(Table 1)
表 1 稳定性评价标准Table 1 Stability assessment standard
| 表 1 稳定性评价标准 Table 1 Stability assessment standard |
表 2(Table 2)
表 2 典型边坡的评价参数Table 2 Assessment parameters of the typical rock slope
| 表 2 典型边坡的评价参数 Table 2 Assessment parameters of the typical rock slope |
根据前面提到的理论建立的评价标准, 计算隶属度(表 3)并运用层次分析法计算各因素间的权重向量(表 4), 经过一致性检验, CR值分别为0.015 8, 0.046 2, 0.060 6和0.038 0, 小于上限值0.1, 符合一致性检验要求.根据上述求出的隶属度值以及AHP方法得到的各指标权重系数, 进行模糊综合计算, 得到边坡模糊综合评价结果(表 5).根据稳定性综合评价的结果进行整体稳定性评价系数?的计算:
表 3(Table 3)
表 3 各评价因素隶属度值Table 3 Membership degree values of assessment factors
| 表 3 各评价因素隶属度值 Table 3 Membership degree values of assessment factors |
表 4(Table 4)
表 4 边坡权重系数分配表Table 4 Weight distribution of assessment indexes for the slope
| 表 4 边坡权重系数分配表 Table 4 Weight distribution of assessment indexes for the slope |
表 5(Table 5)
表 5 抚顺露天矿典型边坡评价结果Table 5 Assessment results of typical slopes of Fushunstrip mine
| 表 5 抚顺露天矿典型边坡评价结果 Table 5 Assessment results of typical slopes of Fushunstrip mine |
(19) |
为了更直观反映整个边坡中各坡段区域的稳定性情况, 本文结合VTK技术完成了三维岩坡建模可视化(3DRM)模块的C++程序编写, 通过运用Delaunay三角剖分理论算法及颜色映射算法, 实现岩质边坡工程地形及总体稳定性的可视化显示,如图 3所示.从结果可以看出, 采用上述评估方法得到四处坡段均属于稳定的范围, 但总体稳定性系数不同, 表现出不同坡段间稳定性的差异.根据实际的勘测, 四处边坡都没有明显的失稳破坏现象, 都属于稳定的范围内, 但-20E800坡段临近下部区域有局部剥落和裂隙扩展等现象, -200E900坡段下方有起鼓的现象, 上述两处坡段相对于+20E600和+20E1000两坡段明显处于较不稳定的状态,这可能与整体边坡的中上部附近断层以及边坡底部大量涌水有关.结果表明,稳定性综合评价计算可以反映正确稳定性关系, 计算与工程现场的探测与调查情况基本相同.
图 3(Fig. 3)
图 3 抚顺露天矿E600~E1000坡段可视化Fig.3 Visualization of Fushun strip mine E600~E1000 |
3 结论本研究将块体理论与模糊综合评价技术相结合, 对岩质边坡的稳定性进行了综合评价, 通过非接触测量手段搜集现场资料, 对岩质边坡的关键块体进行搜索, 并且引入模糊理论和层次分析理论, 克服了传统判定方法的主观性.同时, 将块体理论与模糊综合评价相结合, 确定岩坡的稳定性参数, 弥补了块体理论考虑因素单一的缺点, 可以很好反映岩体的稳定性特性.本文提出的评价模型在抚顺露天矿项目中应用效果良好, 基本符合工程实际情况, 可以为后续的矿山边坡加固及灾害防治提供可靠依据.
参考文献
[1] | Feng X T, Zhao H B, Li S J. A new displacement back analysis to identify mechanical geo-material parameters based on hybrid intelligent methodology[J].International Journal for Numerical and Analytical Methods in Geomechanics, 2004, 28(11): 1141–1165.DOI:10.1002/(ISSN)1096-9853 |
[2] | Kang F, Han S X, Salgado R, et al. System probabilistic stability analysis of soil slopes using Gaussian process regression with Latin hypercube sampling[J].Computers and Geotechnics, 2015, 63(9): 13–25. |
[3] | Sahoo S, Dhar A, Kar A. Environmental vulnerability assessment using grey analytic hierarchy process based model[J].Environmental Impact Assessment Review, 2016, 56(1): 145–154. |
[4] | 王述红, 杨勇, 王洋, 等. 基于数字摄像测量的开挖空间模型及不稳块体的快速识别[J].岩石力学与工程学报, 2010, 29(sup 1): 3432–3438. ( Wang Shu-hong, Yang Yong, Wang Yang, et al. Spatial modeling and quick identification of unstable rock blocks based on digital photogrammetry[J].Chinese Journal of Rock Mechanics and Engineering, 2010, 29(sup 1): 3432–3438.) |
[5] | Goodman R E, Shi G H. Block theory and its application to rock engineering[M]. Englewood Cliff: Prentice Hall, 1985. |
[6] | Lu J. Systematic identification of polyhedral rock blocks with arbitrary joints and faults[J]. Computers and Geotechnics, 2002, 29(1): 49–72.DOI:10.1016/S0266-352X(01)00018-0 |
[7] | Wang S H, Ni P P, Guo M D. Spatial characterization of joint planes and stability analysis of tunnel blocks[J].Tunnelling and Underground Space Technology, 2013, 38(1): 357–367. |
[8] | Wang S H, Ni Y, Viete D R, et al. Rock slope fracture process using a combined block theory and discontinuous deformation analysis approach based on the strength reduction method[J].Disaster Advances, 2014, 7(5): 76–83. |
[9] | 夏元友, 朱瑞赓. 不稳定边坡治理方案的多层次模糊群决策[J].自然灾害学报, 1998, 8(1): 88–91. ( Xia Yuan-you, Zhu Rui-geng. Multi-fuzzy decision for unstable slope strengthening[J].Journal of Natural Hazards, 1998, 8(1): 88–91.) |
[10] | 张勇慧, 李红旭, 盛谦, 等. 基于模糊综合评判的公路岩质边坡稳定性分级研究[J].岩土力学, 2010, 31(10): 3152–3156. ( Zhang Yong-hui, Li Hong-xu, Sheng Qian, et al. Study of stability gradation of highway rock slopes based on fuzzy comprehensive evaluation[J].Rock and Soil Mechanics, 2010, 31(10): 3152–3156.) |