基于无人机摄影测量的露天矿边坡非连续变形分析.pdf

Series No.567September 2023 金 属 矿 山METAL MINE 总 第567 期2023 年第 9 期收稿日期 2020-06-20基金项目 国家自然科学基金项目(编号:51474050)。作者简介 李金威(1983),男,讲师,硕士。基于无人机摄影测量的露天矿边坡非连续变形分析李金威(黄河交通学院继续教育学院,河南 焦作 454950)摘 要 岩体结构面在很大程度上影响了露天矿岩质边坡的稳定性。为提高岩质边坡结构面信息获取的精度及获取动态破坏过程的能量耗散问题,首先借助无人机非接触摄影测量的优势,提出无人机多层次全方位摄影测量方法获取露天矿边坡结构面信息;其次对动态破坏过程施加黏性阻尼,实现能量耗散,并与试验结果进行对比;最后应用所提方法对某露天矿岩质边坡进行非连续变形分析(Discontinuous Deformation Analysis,DDA)。研究结果表明:多层次全方位无人机摄影测量方法对于边坡结构面信息的采集精度良好,施加 2.45%2.53%黏性阻尼的非连续变形分析结果与试验结果基本一致。对某露天矿岩质边坡实际破坏过程的模拟效果良好,且 DDA 分析得出的滑落面位置与现场位移监测分析结果一致;滑落面安全系数为 1.167,处于失稳临界状态,故应提前剥离上覆预失稳岩体,从源头上解决该边坡的滑坡灾害问题。所提方法也适用于其他结构面发育的边坡体分析。关键词 露天矿岩质边坡 无人机 结构面 DDA 能量耗散 多层次全方位 中图分类号TD325,P237 文献标志码A 文章编号1001-1250(2023)-09-134-06DOI 10.19614/ki.jsks.202309018Discontinuous Deformation Analysis of Open-pit Slope Based on Photogrammetry of Unmanned Aerial VehicleLI Jinwei(College of Continuing Education,Huanghe Jiaotong University,Jiaozuo 454950,China)Abstract The rock mass structural plane affects the stability of rock slope in open-pit mine to a great extent.In order to improve the accuracy of rock slope structural plane information acquisition and energy dissipation in dynamic failure process,a multi-layer and omni-directional photogrammetry method is proposed to obtain surface information of open-pit slope with the ad-vantage of unmanned aerial vehicle(UAV)non-contact photogrammetry.Secondly,viscous damping is applied to the dynamic failure process to achieve energy dissipation,and the results are compared with the test results.Finally,the proposed method was applied to Discontinuous Deformation Analysis(DDA)of a rock slope in an open-pit mine.The study results show that the multilayer and omnidirectional UAV photogrammetry method has a good accuracy in collecting the information of slope structur-al surface,and the analysis results of discontinuous deformation with 2.45%2.53%viscous damping are basically consistent with the experimental results.The simulation results of the actual failure process of a rock slope in an open-pit mine are good,and the position of the sliding surface obtained by DDA analysis is consistent with the results of in-situ displacement monito-ring.The safety factor of sliding surface is 1.167,which is in the critical state of instability.Therefore,the overlying pre-unsta-ble rock mass should be stripped in advance to solve the problem of landslide hazard from the source.The proposed method can also be applied to the slope body analysis with other structural plane development.Keywords open-pit rock slope,UAV,structural plane,DDA,dissipation of energy,multi-level and omni-directional 露天矿岩质边坡是由众多节理面、裂隙面、断层面等地质构造切割岩体所形成的复杂结构体系,与土质边坡构造具有较大差异,其稳定性不仅受制于岩石本身的强度,而且在很大程度上依赖于岩体结构面的力学性质,故岩质边坡结构面信息获取的准确性在其稳定性分析中尤为重要1-2。传统结构面信息测量方法一般通过罗盘等工具直接对边坡进行接触测量,该方法工作效率较低、工作量偏大,具有一定的危险性,且对于高陡边坡,测量人员难以进行接触实测3-4。工程领域中三维激光和431近景摄影测量技术逐渐被广泛应用,LI 等5利用激光扫描技术完成边坡三维建模,并根据获取的结构面信息进行了深度学习分组。王述红等6将近景摄影测量技术应用于岩体结构面信息采集,并据此建立岩体三维模型,实现了岩体的稳定性分析与关键块体治理。非连续性是岩体变形过程固有的属性。DDA 方法7具有严密的理论及高效的计算方法,在模拟大变形、大位移等不连续变形方面具有独特优势。DDA方法在本质上与有限元方法(FEM)类似,它使用有限元类型的网格,但不同的是其所有块体都是真正的相互独立块,并由预先设置的不连续性构造限制其运动。DDA 块体可以是任意凸形或凹形,当块体间发生接触时,利用摩尔-库伦准则对接触界面进行处理,并且通过加载函数及步长设定平衡方程并求解。在FEM 方法中,未知数是所有节点的自由度之和,在DDA 方法中,未知数是所有块体的自由度之和,因此从理论上讲,DDA 方法是 FEM 方法的推广8。基于 DDA 独具的优势,近年来许多研究人员对其初始代码进行了一些修改,以提高其性能并解决DDA 中的一些问题。AMIR 等9将 DDA 中的现有块体作为可变形盘引入系统中进行求解,提高了求解效率。JIANG 等10提出了一种有效方法来解决由于大尺度转动而导致的体积增加误差问题,取得了良好的效果。MORGAN 等11提出了一种利用非连续变形分析进行水力压裂的地质力学模型,与 Griffith 分析模型能够很好地匹配。FU 等12利用非连续变形分析方法开发了一种用于地下洞穴挖掘的系统。WU等13研究了大角度转动条件下的弹性扰动问题。YU 等14使用高阶位移函数描述块体的运动,提高了算法精度。ZHENG 等15引入了离散体结构和组件的双重形式解决了由于开闭迭代和应用硬弹簧而引起的问题。FAN 等16通过虚功原理提高了变形计算的准确性,降低了大角度转动的误差。在 DDA 能量耗散研究中,HATZOR 等17-18、TSESARSKY 等19深入探讨了 DDA 动力系数与其运算结果的关系,证明了在 DDA 中添加适当阻尼会使得其模拟与实际更加吻合。KOO 等21通过在 DDA 计算步骤中添加校正函数来减少由于大块体旋转引起的误差,并通过增加质量比例阻尼来考虑能量损失,以满足大块体平移要求。姜清辉等21在模拟准静态过程时,运用变分机理来考虑阻尼力,以此表征其能量耗散形式。刘永茜等22采用自动调整时间步长的方法,在不同速度运动阶段选用不同的阻尼运算规则,取得了较好的运算结果。露天矿边坡稳定性分析的关键在于求得安全系数及其滑落面,而这一求解过程的关键又在于岩体结构面信息的精准获取;再者岩质边坡破坏是一个动态过程,其中的能量耗散有必要充分考虑。以往研究鲜有针对结构面信息精准获取并将其用于非连续变形分析计算中,多数默认其精度满足计算要求,但在实际情况下结构面形态各异,须重点关注信息采集过程及其精度,并考虑其动态破坏过程的能量耗散,才可保证后续非连续变形分析的准确性。为了对上述问题提供解决方案,本研究借助无人机非接触摄影测量优势,提出多层次全方位无人机摄影测量方法获取边坡结构面信息,保证其结构面信息的采集精度,同时施加适当的黏性阻尼,确保非连续变形分析过程的能量耗散,并将改进后的方法应用到某露天矿岩质边坡的非连续变形分析中。1 结构面信息采集 1.1 多层次全方位摄影测量以某露天矿岩质边坡为研究对象,该边坡基岩裸露,节理发育,且时常发生块体滑落或局部滑塌,失稳趋势明显,对该区域造成了极大安全隐患。坡体主要由条带状砂页岩和砂岩组成,上覆 0.6 m 左右岩体强风化严重,岩体总体呈东西走向,坡高约 48 m,边坡角为 7686。采用无人机多层次全方位摄影测量采集边坡结构面信息,设置合理航道,通过三维摄影获取结构面信息,如图 1 所示。图 1 多层次全方位摄影测量Fig.1 Multilevel and omnidirectional photogrammetric measurement其中多层次指在距离边坡面 1、3、5 m 处分别构建无人机飞行平面;全方位指在各个方位对边坡体进行摄影测量,区别于以往地面定点摄影测量,用覆盖边坡体的百分比来表征该变量。相较于传统测量方法,多层次全方位摄影测量可以最大程度地还原带有坐标的实物模型,并运用计算机算法获取结构面信息,精度高,对于复杂边坡体结构面信息采集有很强的鲁棒性。多层次全方位摄影测量精度见表 1。由表 1 可知:无人机多层次全方位摄影测量精度与多层次量及全方位率成正比,当多层次量为 3、全方位率为100%时,与全站仪测量的相对误差小于531 李金威:基于无人机摄影测量的露天矿边坡非连续变形分析 2023 年第 9 期表 1 多层次全方位摄影测量精度Table 1 Precision of multilevel omnidirectional photogrammetry%多层次不同全方位率下与全站仪测量的相对误差5075100多层次量112.817.925.65210.626.343.5738.223.390.951%,满足计算要求。1.2 结构面信息计算根据结构面产状倾向、倾角概念可得,其产状由从属面的法向量控制,故任意结构面可表示为Z=AX+BY+C,(1)式中,A、B、C 分别指代结构面参数,且从属面的法向量 n 可由(-A,-B,1)表示。根据三点成面原则,将不在同一直线的 n 个点(n3)提取可得一结构面,结构面可表示为X1Y11X2Y21X3Y31XnYn1ABC =Z1Z2Z3Zn,(2)则A,B,CT可表示为ABC =X1Y11X2Y21X3Y31XnYn1TX1Y11X2Y21X3Y31XnYn1-1X1Y11X2Y21X3Y31XnYn1TZ1Z2Z3Zn,(3)进而确定其产状信息倾向、倾角分别为=arctan B/A,(4)=arctanA2+B2().(5)边坡实际点云坐标输出结果及出露迹线提取结果如图 2、图 3 所示。图 2 点云坐标输出结果Fig.2 Output results of point cloud coordinates图 3 出露迹线提取结果Fig.3 Extraction results of exposed trace lines根据该计算过程所得结构面产状信息及其与人工测量所得信息对比结果见表 2,对比结果显示两者误差均在 0.70%以内,反映出本研究算法具有良好的精度。表 2 结构面产状信息的无人机摄影测量与人工测量结果对比Table 2 Comparison results of structural plane information between UAV photogrammetry and manual measurement结构面编号产状信息无人机摄影测量人工测量倾向/()倾角/()迹线/m倾向/()倾角/()迹线/m相对误差/%倾向 倾角 迹线133.535.46.5433.8 35.8 6.54 0.59 0.84 0.002109.178.86.68108.6 79.1 6.67 0.28 0.25 0.303263.183.48.52263.1 83.8 8.61 0.00 0.36 0.104291.783.14.25291.4 83.4 4.24 0.10 0.24 0.265282.282.82.87282.3 82.7 2.88 0.04 0.24 0.356265.882.81.21266.0 82.9 1.20 0.04 0.00 0.00796.573.15.8796.8 73.1 5.91 0.31 0.14 0.348250.277.81.46251.0 78.0 1.45 0.28 0.13 0.689284.689.92.05285.0 89.2 2.06 0.11 0.22 0.481026.4089.51.9426.4 89.2 1.94 0.00 0.45 0.002 考虑能量耗散的 DDA 法2.1 DDA 基本原理任意块体 i 中任意点(x,y)处的位移(u,v)在两个维度上与 6 个位移变量相关,即:D=u0,v0,r0,x,y,xyT,(6)式中,(u0,v0)为块体内特定点(x0,y0)的刚性位移;r0为块体绕特定点(x0,y0)的转动角度;x,y为该块体的正应变;xy为该块体的切应变。对于块体内任意点(x,y),其位移(u,v)可以表示为631总第 567 期 金 属 矿 山 2023 年第 9 期uv=10-(y-y0)(x-x0)0(y-y0)/201(x-x0)0(y-y0)(x-x0)/2u0v0r0 xyxy.(7)以最小势能原理为基础,组建的块体系统整体平衡方程为K11K12K1nK21K22K2nKn1Kn2KnnD1D2Dn=F1F2Fn,(8)式中,Kij是 6 6 子矩阵;Kii为块体的几何与物理参数;Kij(ij)指块体间接触参数;Di和 Fi为6 1 子矩阵;Di为块体 i 的变形参数(d1i,d2i,d3i,d4i,d5i,d6i);Fi为块体 i 上分拨给 6 个变形参数的应力矩阵。2.2 能量耗散岩体破坏过程是不可逆过程,该过程存在岩体变形损耗的能量及结构面相对滑动所损失的能量等,即过程中伴随着能量耗散。非连续变形分析过程可施加适当的黏性阻尼来完成能量耗散机制的设定,并通过经典斜坡案例测试动量耗散阈值,岩体参数取值见表 3,边坡模型如图 4 所示。与工程实验室试验动态位移结果对比,结果(图 5)显示:在设置 2.45%2.53%黏性阻尼时与室内试验结果最为匹配,故选取2.5%黏性阻尼作为能量耗散阈值,动力系数 k 取0.985,用以完成能量耗散设定,使得边坡动态破坏过程与实际破坏过程更加吻合。表 3 边坡模型岩体参数Table 3 Mechanical parameters of the slope rock参数容重/(kN/m3)弹性模量/MPa黏聚力/kPa泊松比取值27 0001 000400.35图 4 滑块沿斜坡滑动模型Fig.4 Model of sliding block along with slope通过上述计算分析,考虑能量耗散的 DDA 程序运算流程如图 6 所示。3 露天矿岩质边坡非连续变形分析3.1 计算模型及参数将无人机采集的结构面信息输入至 DDA,边坡体图 5 不同黏性阻尼下 DDA 与试验位移结果对比Fig.5 Comparison between DDA and experimental displacement under different viscous damping图 6 DDA 程序运算流程Fig.6 Calculating flow of the DDA procedures由结构面分割成块体,作为 DDA 计算的基本单元,DDA 具体数值计算分析模型如图 7 所示,其中 1#、2#孔为工程现场测斜孔,内置测斜仪用以获取边坡位移,边坡模型岩体参数取值见表 4。3.2 现场位移监测及边坡非连续变形分析测斜孔布置参考图 7,测点间隔为 1 m,1#、2#测斜孔内不同深度所发生的累积位移如图 8 所示。由图 8 可知:1#孔附近深度为 012 m 的岩体位移随时间变化较大,深度 12 m 附近以下的岩体变形量较小,且趋于稳定状态,说明 1#孔内 012 m 深度附731 李金威:基于无人机摄影测量的露天矿边坡非连续变形分析 2023 年第 9 期图 7 露天矿岩质边坡 DDA 数值计算分析模型Fig.7 Numerical calculation and analysis model of DDA for open pit rock slope表 4 边坡模型岩体参数Table 4 Mechanical parameters of the slope rock材料名称容重/(kN/m3)弹性模量 E/MPa内聚力 c/MPa内摩擦角/()砂页岩21 00016 0002.019砂岩22 00014 2003.320近岩体已发生相对滑动,随时有发生失稳的可能,且滑动面就在该孔深度为 12 m 附近位置。随着时间的变化,2#孔内前 8 m 深度左右测点总位移呈现规律式上升,且越靠近上坡面的岩体位移量越大,深度 8 m以下岩体位移趋于稳定,说明滑动面位置在深度 8 m附近,且根据位移上升的规律推断该位置岩体存在倾倒破坏可能。图 8 测斜孔测点累积总位移值Fig.8 Cumulative displacements of monitoring points of the inclination hole为预测失稳面确切位置及相应安全系数,并采取相应措施治理该露天矿边坡,本研究运用考虑能量损耗的非连续变形分析方法进行计算,具体动态失稳过程如图 9 所示。图 9 公路岩质边坡非连续变形分析过程Fig.9 Discontinuous deformation analysis of highway rock slope 由图 9 可知:在到达第 1.67106时步前,边坡整体处于稳定状态,当到达第 1.67106时步时,即以折减系数为 1.167 参与程序运算时,边坡触发变形失稳过程,各块体发生不同程度位移,且后续变形愈演愈烈,最终发生图 9 中所示的大型滑坡,该边坡安全系数为 1.167,危险性较高。分析可得,本研究方法计算出的滑落面位置与根据 1#、2#测斜孔所得数据预测的滑落面位置基本一致,从现场监测数据角度印证了该方法的可行性及精度。该露天矿岩质边坡案例可归为牵引式滑动破坏模式,滑坡体基本为强风化砂页岩,存在极大失稳可能性,预失稳岩体均为坡体表层 4 m 以内岩体。考虑上述因素及实际情况,建议按照本研究方法计算得出的滑落带位置,剥离上覆预失稳岩体,从源头上解决831总第 567 期 金 属 矿 山 2023 年第 9 期该露天矿边坡的滑坡灾害问题。4 结 论(1)提出了一种多层次全方位无人机摄影测量获取边坡结构面信息的方法,为边坡结构面信息采集提供了一种新手段。(2)考虑了露天矿边坡非连续变形分析过程中的能量耗散,当施加 2.45%2.53%黏性阻尼时模拟边坡动态破坏过程与实际破坏过程吻合度较好。(3)基于多层次全方位无人机摄影测量及考虑能量耗散的非连续变形分析方法,分析了某露天矿岩质边坡的稳定性,讨论了其最危险滑落面位置,并考虑到其安全系数及预失稳岩体体积等因素。建议按照本研究方法计算得出的滑落带位置,剥离上覆预失稳岩体,从源头上解决该边坡的滑坡灾害问题。参 考 文 献1 黄笃学,董鑫,王建曦,等.露天矿节理参数敏感性研究及边坡稳定性分析J.金属矿山,2020(3):64-70.HUANG Duxue,DONG Xin,WANG Jianxi,et al.Study on the sensi-tivity of joint parameters and slope stability analysis of open-pit mineJ.Metal Mine,2020(3):64-70.2 陈欢欢,雍睿,刘子航.基于 Barton-Bandis 非线性破坏准则的边坡破坏概率分析J.金属矿山,2020(11):44-52.CHEN Huanhuan,YONG Rui,LIU Zihang.Analysis of slope failure probability based on Barton-Bandis nonlinear failure criterionJ.Metal Mine,2020(11):44-52.3 李德鑫,闫志刚,孙久运.基于无人机视觉的河道漂浮垃圾分类检测技术研究J.金属矿山,2021(9):199-205.LI Dexin,YAN Zhigang,SUN Jiuyun.Study on classification and de-tection technology of river floating garbage based on UAV visionJ.Metal Mine,2021(9):199-205.4 廉旭刚,蔡音飞,胡海峰.我国矿山测量领域三维激光扫描技术的应用现状及存在问题J.金属矿山,2019(3):35-40.LIAN Xugang,CAI Yinfei,HU Haifeng.Application status and exist-ing problems of 3D laser scanning technique in mine surveying in ChinaJ.Metal Mine,2019(3):35-40.5 LI H B,LI X W,LI W Z,et al.Quantitative assessment for the rock-fall hazard in a post-earthquake high rock slope using terrestrial laser scanningJ.Engineering Geology,2019,248(8):1-13.6 王述红,朱承金,张紫杉,等.基于动态强度折减 DDA 法的边坡多滑面稳定性分析J.煤炭学报,2019,44(4):1084-1091.WANG Shuhong,ZHU Chengjin,ZHANG Zishan,et al.Stability a-nalysis of multi-slip surface of slope based on dynamic strength re-duction DDA methodJ.Journal of China Coal Society,2019,44(4):1084-1091.7 SHI G H,GOODMAN R E.Two dimensional discontinuous deforma-tion analysisJ.International Journal for Numerical and Analytical Methods in Geomechanics,1985,9:541-556.8 SHI G H.Discontinuous Deformation Analysis:a New Numerical Model for the Statics and Dynamics of Block SystemsD.Berke-ley:University of California,1988.9 AMIR R B S,BAGTZOGLOU A C.Non-rigid disk-based DDA with a new contact modelJ.Computers and Geotechnics,2013,49:25-35.10 JIANG W,ZHENG H.An efficient remedy for the false volume ex-pansion of DDA when simulating large rotationJ.Computers and Geotechnics,2015,70:18-23.11 MORGAN W E,ARAL M M.An implicitly coupled hydro-geome-chanical model for hydraulic fracture simulation with the discontin-uous deformation analysisJ.International Journal of Rock Me-chanics and Mining Sciences,2015,73:82-94.12 FU X,SHENG Q,ZHANG Y,et al.Investigations of the sequential excavation and reinforcement of an underground cavern complex u-sing the discontinuous deformation analysis methodJ.Tunnelling and Underground Space Technology,2015,50:79-93.13 WU J H.The elastic distortion problem with large rotation in dis-continuous deformation analysisJ.Computers and Geotechnics,2015,69:352-364.14YU Y,YIN J.Some modifications to the process of discontinuous deformation analysisJ.Journal of Rock Mechanics and Geotech-nical Engineering,2015,7(1):95-100.15 ZHENG H,ZHANG P,DU X.Dual form of discontinuous deforma-tion analysisJ.Computer Methods in Applied Mechanics and En-gineering,2016,305:196-216.16FAN H,ZHENG H,ZHAO J.Discontinuous deformation analysis based on strain-rotation decompositionJ.International Journal of Rock Mechanics and Mining Sciences,2017,92:19-29.17 HATZOR Y H,FEINTUCH A.The validity of dynamic block dis-placement prediction using DDAJ.International Journal of Rock Mechanics and Mining Sciences,2001,38(4):599-606.18 HATZOR Y H,ARZI A A,ZASLAVSKY Y,et al.Dynamic stability analysis of jointed rock slopes using the DDA method:King Herod”s Palace,Masada,IsraelJ.International Journal of Rock Mechan-ics and Mining Sciences,2004,41(5):813-832.19 TSESARSKY M,HATZOR Y H,SITAR N.Dynamic displacement of a block on an inclined plane:analytical,experimental and DDA resultsJ.Rock Mechanics and Rock Engineering,2005,38(2):153-167.20 KOO C Y,CHERN.Modification of the DDA method for rigid block problemsJ.International Journal of Rock Mechanics and Mining Sciences,1998,35(6):683-693.21 姜清辉,周创兵.岩土工程不连续变形分析计算中的若干问题J.岩石力学与工程学报,2007,26(10):2014-2026.JIANG Qinghui,ZHOU Chuangbing.Some issues in discontinuous deformation analysis(DDA)for geotechnical engineeringJ.Chi-nese Journal of Rock Mechanics and Engineering,2007,26(10):2014-2026.22 刘永茜,杨军.一种改进步长自调的非连续变形分析法J.岩土力学,2011,32(8):2544-2548.LIU Yongqian,YANG Jun.An improved DDA method with adaptive step sizeJ.Rock and Soil Mechanics,2011,32(8):2544-2548.931 李金威:基于无人机摄影测量的露天矿边坡非连续变形分析 2023 年第 9 期