Calculation method of slope stability coefficients based on spatial multi-profile slopes
-
摘要:
在边坡稳定性分析中,通常需要对多个剖面的边坡稳定性进行单独计算。为了提高多剖面边坡稳定性系统计算的效率,本研究提出了一种基于空间多剖面边坡稳定性系数计算的方法。首先,将二维剩余推力法扩展至空间剖面,形成空间剖面的剩余推力法;其次,通过在三维边坡工程地质模型上实现剖面图的自动生成和管理,进一步优化了计算流程;为了实现快速高效计算边坡稳定性系数,引入了基于多线程并行计算的方法,并结合空间剖面的剩余推力法,对空间多剖面边坡稳定性系数进行了更为准确和高效的计算。以内蒙古锡林浩特某露天矿为应用实例,建立内排土场边坡三维边坡工程地质模型,在此基础上,自动生成7个空间剖面图,采用多线程并行计算出7个剖面的稳定性系数,最后,通过可视化手段对计算结果进行了直观展示。采用多线程并行计算多剖面边坡稳定性系数能够充分利用计算资源,提高计算效率,同时,在实际工程应用中进一步证实了该方法的可行性和适用性。
Abstract:In slope stability analysis, the conventional practice involves independently calculating slope stability for various cross-sections.
Objective To enhance the efficiency of slope stability system calculations for multiple profiles, this study introduces a method based on spatial slope stability coefficient calculations. Initially, the two-dimensional residual thrust method is expanded to incorporate spatial profiles, resulting in the residual thrust method for spatial profiles.
Methods Subsequently, the calculation process is further optimized by realizing the automatic generation and management of the profile on the 3D slope engineering geological model. To achieve swift and efficient calculation of slope stability coefficients, a method based on multi-threaded parallel computation is introduced. Combining this with the residual thrust method for spatial profiles allows for more precise and efficient computation of stability coefficients for multiple spatial profiles.
Results Using a specific open-pit mine in Xilinhot, Inner Mongolia as a case study, we establish a three-dimensional slope engineering geological model for the internal dump site slope. Seven spatial profiles are automatically generated and stability coefficients are calculated using multi-threaded parallel computation. Finally, we present visual results through visualization techniques.
Conclusion The research findings suggest that adopting multi-threaded parallel computation for computing slope stability coefficients of multiple profiles can fully utilize computing resources and significantly improve computational efficiency. Furthermore, practical implementation of this method in engineering projects affirms its applicability and feasibility.
-
图 1 第i个条块受力状态
Wi . 第i条块的重量;Di . 第i条块的剩余推力;Di-1. 第i-1条块剩余推力;Wp,i ,Wp,i-1. 作用于条块的水压力;αi. 第i条块的底滑面与水平线的夹角;Ti . 第i条块底面的剪切力;Ni . 第i条块的法向力;(tx,i,ty,i),(tx,i-1,ty,i-1),(bx,i,by,i)和(bx,i-1,by,i-1). 分别为条块4个顶点坐标;oy. 高程方向;ox. 垂直于oy;(Wx,i,Wy,i),(Wx,i-1,Wy,i-1). 均为条块中间部分的坐标
Figure 1. Force state of the ith block
表 1 边坡剖切线数据结构
Table 1. Data structure of slope profile lines
字段名 数据类型 说明 SectionName navrchar 剖面名称 SlopeCuttingtangent geometry 边坡剖切线 TimeProperty datatime 边坡剖切线存储的日期 表 2 岩层多边形图元数据结构
Table 2. Rock polygon map metadata structure
字段名 数据类型 说明 SectionName navrchar 剖面名 Rockpolygon geometry 岩层多边形 sectionElementPara navrchar 剖面图元属性参数 ColorProperty string 岩层多边形属性颜色 表 3 岩石物理力学指标
Table 3. Rock physico-mechanical indexes
地层 黏聚力C/MPa 内摩擦角φ/(°) 容重γ/(kN·m−3) 排弃物料 0.0151 14.0 18.0 第四系 0 26.9 20.8 新近系、古近系 0.120 24.0 19.3 煤 0.100 30.0 12.2 泥岩 0.050 25.0 19.6 -
[1] 瞿婧晶,陆燕,吴曙亮,等. 基于不同方法的江苏镇江地区下蜀土边坡稳定性分析与评价[J]. 中国地质灾害与防治学报,2021,32(1):35-42.QU J J,LU Y,WU S L,et al. Evaluation of Xiashu loess slope stability in Zhenjiang area using different methods[J]. The Chinese Journal of Geological Hazard and Control,2021,32(1):35-42. (in Chinese with English abstract [2] 林姗,郭昱葵,孙冠华,等. 边坡稳定性分析的虚单元强度折减法[J]. 岩石力学与工程学报,2019,38(增刊2):3429-3438.LIN S,GUO Y K,SUN G H,et al. Virtual element strength reduction method for slope stability analysis[J]. Chinese Journal of Rock Mechanics and Engineering,2019,38(S2):3429-3438. (in Chinese with English abstract [3] 尉迟衍春. 复杂边坡多失效模式的确定方法及其工程应用[D]. 山东青岛:青岛理工大学,2018.YUCHI Y C. Method fordetermining multiple failure modes of complexslope andits engineering application[D]. QingdaoShandong:Qingdao University of Technology,2018. (in Chinese with English abstract [4] 陈亚烽,陈国庆,严明,等. 基于一阶线性应变软化理论的边坡稳定性研究[J]. 地质科技通报,2022,41(6):180-188.CHEN Y F,CHEN G Q,YAN M,et al. Study on slope stability considering first-order linear strain-softening theory[J]. Bulletin of Geological Science and Technology,2022,41(6):180-188. (in Chinese with English abstract [5] 易富,孟兴涛,赵文华,等. 基于因素空间的露天采坑边坡稳定性评价[J]. 地质科技通报,2023,42(5):1-9.YI F,MENG X T,ZHAO W H,et al. Evaluation of slope stability of open pit based on factor space[J]. Bulletin of Geological Science and Technology,2023,42(5):1-9. (in Chinese with English abstract [6] 王伟,邹艺璇,李尧,等. 基于GIS与Scoops 3D的水动力型滑坡区域稳定性分析[J]. 防灾减灾工程学报,2024,44(1):79-89.WANG W,ZOU Y X,LI Y,et al. Stability analysis of hydrodynamic landslides in the reservoir area based on GIS and Scoops 3D[J]. Journal of Disaster Prevention and Mitigation Engineering,2024,44(1):79-89. (in Chinese with English abstract [7] 魏婷婷,范文,于宁宇. 基于典型剖面分析的区域斜坡稳定性分区方法[J]. 工程地质学报,2017,25(6):1518-1526.WEI T T,FAN W,YU N Y. Method of zoning of regional slope stability based on typical profiles analysis[J]. Journal of Engineering Geology,2017,25(6):1518-1526. (in Chinese with English abstract [8] 张年学,李晓,盛祝平,等. 多剖面合力法分析滑坡稳定性[J]. 工程地质学报,2017,25(5):1190-1204.ZHANG N X,LI X,SHENG Z P,et al. Stability analysis of landslide using joint forces method on multiple profiles[J]. Journal of Engineering Geology,2017,25(5):1190-1204. (in Chinese with English abstract [9] 陈应显,王鹏飞. 露天矿含弱层边坡三维危险滑面构建及稳定性计算[J]. 安全与环境学报,2024,24(3):986-992.CHEN Y X,WANG P F. Three-dimensional dangerous sliding surface of slopes with a weak layer in open pit construction and stability calculation[J]. Journal of Safety and Environment,2024,24(3):986-992. (in Chinese with English abstract [10] LIU X,LI D Q,CAO Z J,et al. Adaptive Monte Carlo simulation method for system reliability analysis of slope stability based on limit equilibrium methods[J]. Engineering Geology,2020,264:105384. [11] 范军富,刘志斌. 露天煤矿边坡稳定性(剩余推力法)计算机辅助设计[J]. 露天采矿技术,2005,20(2):21-23. doi: 10.3969/j.issn.1671-9816.2005.02.009FAN J F,LIU Z B. Computer aided design of slope stability in open-pit coal mine[J]. Opencast Coal Mining Technology,2005,20(2):21-23. (in Chinese with English abstract doi: 10.3969/j.issn.1671-9816.2005.02.009 [12] 韩彦峰,王家序,周广武,等. 滑动轴承混合润滑多线程并行计算数值方法[J]. 华中科技大学学报(自然科学版),2016,44(6):7-12.HAN Y F,WANG J X,ZHOU G W,et al. Mixed lubrication study of multi-thread parallel computing algorithm for journal bearings[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition),2016,44(6):7-12. (in Chinese with English abstract [13] 张宏莉,田耕,胡铭曾. 多线程技术与并行计算[J]. 计算机科学,1999,26(9):48-50.ZHANG H L,TIAN G,HU M Z. Multi-thread and parallel computing[J]. Computer Science,1999,26(9):48-50. (in Chinese with English abstract [14] 王永亮,王娟. 流体驱动岩体裂缝扩展计算新技术:多线程并行计算[C]//北京力学会. 北京力学会第二十八届学术年会论文集(上).S. l. ]:[S. n. ],2022. WANG Y L,WANG J. New technology of fluid-driven rock fracture extension calculation:Multithreaded parallel computation[C]//Beijing Mechanics Society. Proceedings of the Twenty-eighth Annual Academic Conference of the Beijing Mechanics Society (Above). [S. l. ]:[S. n. ],2022. (in Chinese) [15] 王文鼐,陈邦文,韩鹏,等. 采用多线程并行调度的网络仿真加速[J]. 南京邮电大学学报(自然科学版),2015,35(1):33-37.WANG W N,CHEN B W,HAN P,et al. Network simulation speedup by multi-thread and parallel scheduling[J]. Journal of Nanjing University of Posts and Telecommunications (Natural Science Edition),2015,35(1):33-37. (in Chinese with English abstract [16] 葛梦凡,李志,徐南,等. 三维模型自旋图的多线程并行算法[J]. 小型微型计算机系统,2017,38(10):2369-2373. doi: 10.3969/j.issn.1000-1220.2017.10.037GE M F,LI Z,XU N,et al. Multithreaded parallel algorithm of 3D model spin image[J]. Journal of Chinese Computer Systems,2017,38(10):2369-2373. (in Chinese with English abstract doi: 10.3969/j.issn.1000-1220.2017.10.037 [17] ZHANG W J,ZHENG H,JIANG F Y,et al. Stability analysis of soil slope based on a water-soil-coupled and parallelized Smoothed Particle Hydrodynamics model[J]. Computers and Geotechnics,2019,108:212-225. doi: 10.1016/j.compgeo.2018.12.025 [18] 赵宏坚,周翠英. 基于实体建模的三维地层构造[J]. 岩土力学,2010,31(4):1257-1263. doi: 10.3969/j.issn.1000-7598.2010.04.044ZHAO H J,ZHOU C Y. Solid-oriented 3D-stratum construction[J]. Rock and Soil Mechanics,2010,31(4):1257-1263. (in Chinese with English abstract doi: 10.3969/j.issn.1000-7598.2010.04.044 [19] 田宜平,刘维安,张夏林. 基于等角度变比例投影的矿体轮廓线自动匹配方法研究[J]. 地质科技通报,2020,39(1):175-180.TIAN Y P,LIU W A,ZHANG X L. Automatic matching of ore body contour line based on equal-angle and variable proportion projection[J]. Bulletin of Geological Science and Technology,2020,39(1):175-180. (in Chinese with English abstract [20] 梁师俊. 基于虚拟钻孔的工程地质三维剖切的实现[J]. 地理与地理信息科学,2010,26(1):11-14.LIANG S J. 3D sectioning realization of engineering geology based on virtual borehole[J]. Geography and Geo-Information Science,2010,26(1):11-14. (in Chinese with English abstract [21] 李楠,曾俊钢,肖克炎,等. 一种线框矢量模型表示的地质体折面剖切方法[J]. 测绘科学,2011,36(6):28-31.LI N,ZENG J G,XIAO K Y,et al. Cutting and splitting on geological body with wire frame model[J]. Science of Surveying and Mapping,2011,36(6):28-31. (in Chinese with English abstract [22] 杨凯祥,陈鸿,谢明师. STL模型相交性检测及交线追踪算法[J]. 科学技术与工程,2017,17(1):255-259. doi: 10.3969/j.issn.1671-1815.2017.01.046YANG K X,CHEN H,XIE M S. Intersection test and cross line tracking algorithms of the STL models[J]. Science Technology and Engineering,2017,17(1):255-259. (in Chinese with English abstract doi: 10.3969/j.issn.1671-1815.2017.01.046 -