导 读
巷道是地下金属矿山开拓及生产运输的主要工程,地下巷道围岩稳定性将直接影响矿山的安全生产。通过开展巷道岩体结构分析,获取围岩结构面空间特征,能够解析出结构面所切割的块体特征及失稳情况,识别出巷道块体失稳的危险区域,从而为巷道稳定提供安全保障。现阶段主要采用全人工识别、半人工识别和全自动识别等方法进行岩体结构面识别,不同方法对岩体结构面识别的精度也不同,但普遍存在地下巷道的围岩缺陷与结构面难以准确区分的问题,总体上识别出的结构面数据精度较差。因此,有必要研发一种安全、精确的岩体结构面测量识别技术,为巷道结构面的精准解构和稳定性分析提供数据支持。近年来,利用三维激光扫描技术识别并提取巷道结构面信息已成为巷道稳定性研究热点之一,国内外学者从不同角度相继开展诸多研究并取得一定成果,但是关于提高地下巷道结构面识别精度的技术仍有待深入研究。为了实现对地下巷道围岩结构面点云数据高质量、高精度的识别和优化,华北有色工程勘察院有限公司赵晓明高级工程师和中南大学资源与安全工程学院李杰林副教授联合科研团队采用多种巷道结构面识别方法,以彝良驰宏矿业有限公司毛坪铅锌矿为研究对象,开展了三维激光扫描与结构面信息识别与几何特征提取工作,以期为巷道稳定性分析提供基础数据。
1 基于点云数据的结构面识别方法
通过三维激光扫描获取地下巷道的点云数据,经过去噪、配准和坐标匹配处理后,即可获得完整还原的巷道点云模型,相当于将真实的地下巷道空间场景搬运至计算机内。获取结构面信息首先需要从点云数据中识别结构面,在点云数据中通过全自动、半自动和人工提取的方式来识别巷道围岩的结构面信息。人工识别结构面方法是将现场采用罗盘、直尺等手动测量方式转化至计算机内进行,通过人工手动选取地下巷道出露的结构面点云数据。在点云数据中,结构面可以通过其所出露的点样本进行平面拟合,确定出结构面的综合产状。由于结构面出露的不规则性,部分区域结构面形态变化大,通过计算机算法难以自动识别,而利用人工选取方式往往更精确。相较现场的人工测量,在点云数据中通过人工识别结构面的方法,可使作业人员的工作环境更舒适且平面拟合结果更准确。然而,人工结构面识别方法耗时较长,且不同测量人员获得的岩体结构面数量和产状可能不同,人为因素的主观影响较大。全自动结构面识别方法基于计算机算法实现点云数据的自动识别,利用三维激光扫描设备所获取的地下巷道点云数据,快速、智能地获取点云模型中的结构面信息。全自动识别方式能够极大地减少结构面识别中人为因素的主观影响,减少作业人员的工作时间,提高作业效率。但是,全自动识别方式对于地下巷道不规则结构面出露的识别精度不高,且当点云数据量过大时,往往存在程序运行缓慢和冗余数据量大等缺点。半自动结构面识别方法是通过人工干预方式有针对性地选取识别范围,进而指导计算机进行自动识别,可以减少自动识别程序运行的时间,提高识别精度。对于地下金属矿山而言,由于巷道围岩硬度较高,因此普遍采用钻爆法进行掘进,在爆破过程中容易产生岩体损伤,造成次生结构面发育。采用全自动算法难以在点云数据中准确区分构造性结构面与非构造性结构面,而采用人工识别方法存在测量周期较长以及人工测量造成主观偏差等问题。因此,通常选取“人工+半自动”识别相结合的方法进行地下巷道岩体结构面的识别。
毛坪铅锌矿位于云南省昭通市彝良县,巷道内岩体结构类型为层状碎裂结构。该岩组结构面主要为节理、断层及层间软弱泥灰岩夹层。由于受破碎带和软弱夹层的影响,局部地段岩体发生塑性变形,其变形随时间推移具有明显的流变性。该岩组稳定性差,硐室易形成局部冒落和整体变形,在泥质页岩的下盘部位,巷道侧帮发生片状剥落,整体稳定性差。由于调查区域还可能存在巷内结构面出露区域不连续,如采用传统方式进行采样需要多次布点,测量工作繁琐,且测量结果不精确。因此,选择具有数据获取快、精度高和非接触测量等优势的三维激光扫描技术,对毛坪铅锌矿910 m中段运输巷道进行岩体结构面工程调查。选取毛坪铅锌矿南部910 m中段的120 m运输巷道进行现场扫描测量,如图1所示。结合现场扫描区域,采用Maptek SR3三维激光扫描仪进行测量,以5 m间距为一个点位确定架设站点,如图2所示。三维激光能够全面、精确地捕捉到巷道表面的起伏,尽可能减少扫描盲区,从而准确识别巷道内部的结构面特征。在扫描过程中,需选取4X扫描精度,该精度可达到每秒百万点,但是单个站点扫描所需时间较长,扫描完成一个站点所需时间为7 min 24 s。图1 测量区域选取Fig.1 Selection of measurement area图2 测量站点架设Fig.2 Erection of measuring station将现场扫描所获取的初始点云数据,导入Pointstudio点云处理软件中,经过过滤、去噪和拼接等操作,获取所扫描巷道点云数据,如图3所示。本次测量的巷道长度为120 m,所扫描点数为4 430万个,平均点间距为6 mm,距离误差为±3 mm,符合测量的要求。图3 巷道点云数据模型Fig.3 Roadway point cloud data model点云模型的坐标和方位需要对照真实坐标点进行转换,若巷道坐标和方位与实际矿山巷道不一致,将会影响所提取结构面倾向和倾角的准确性。将扫描作业过程中所提取的对照真实坐标进行输入匹配,并与带真实坐标的巷道平面图进行校准(图4),即可得到真实的点云坐标。图4 点云数据坐标匹配Fig.4 Coordinate matching of point cloud data测量总耗时4 h,相较传统的人工工程地质调查,Maptek SR3三维激光扫描仪的扫描速度较快。从扫描数据可以看出,在地下巷道扫描过程中,受扫描区域管线、积水和喷锚支护等遮挡物的影响,所获取的巷道点云数据不可避免地存在部分空白,但点云数据的可视性总体较好,获得的围岩结构面的整体扫描效果较好,如图5所示。图5 围岩实测照片和点云数据扫描结构面特征对比Fig.5 Comparison of structural plane characteristics between measured photos of surrounding rock and point cloud data scanning为了验证基于三维激光扫描点云数据的巷道围岩结构面特征识别方法的准确性、适用性和高效性,矿山现场技术人员使用罗盘对结构面进行了人工测量。将采用人工识别方法获取的结果与人工测量获取的结果进行对比验证,结果如图6所示。图6 结构面特征识别的对比验证Fig.6 Comparative verification of structural plane characteristics identification选取本次实测巷道中的某一个围岩壁面点云数据进行验证。通过点云数据识别测量的单个结构面产状与现场实际同一点位使用罗盘测量的产状进行对比,如图6所示。从获取的点云数据中识别得到岩体结构面的平均倾向为2°,平均倾角为60°。现场人工测量结构面的平均倾向为4°,平均倾角为64°。综合来看,基于点云数据的巷道围岩结构面特征人工识别方法与现场使用罗盘测量获得的产状误差在5°以内,说明基于点云数据的巷道围岩结构面特征人工识别方法能够满足测量精度要求。选取实测巷道中的某一个围岩壁面点云数据进行验证。通过在点云数据中布置4个测量点,识别点云数据中的单个结构面产状并与现场使用罗盘测量获取的产状进行对比(图7),二者误差结果如表1所示。由表1可知,测量点1、测量点2和测量点3的产状与现场罗盘测量结果的误差较大,倾向最大误差为24°,倾角最大误差为14°,而测量点4与现场罗盘测量的产状误差在5°以内。图7 结构面拟合验证Fig.7 Fitting verification of structural plane表1 人工选取结构面识别与传统罗盘测量结果误差Table 1 Error of manually selected structural plane identification and traditional compass measurement results
出现上述现象的主要原因是岩体结构面存在一定程度的起伏,因此在同一个结构面中的不同测点所获取的结构面产状会存在差异。同时,该结果也间接反映了传统的罗盘测量方式容易产生人为选点的客观误差,且测点数量有限。然而,在基于点云数据的人工识别方法中,可以方便地选取全结构面区域进行结构面产状测量,实现多点测量,减少了不同测量点及结构面起伏度带来的产状测量误差,测量结果更加准确。选取现场扫描巷道的某一段围岩壁面进行结构面半自动识别与人工识别验证,如图8所示。通过比较半自动识别与人工识别结果可知:半自动识别方法花费时间为1 min,共识别出295条结构面数据,而人工识别花费时间为10 min,仅识别出75条件结构面数据。其原因是,半自动识别是根据计算机程序的查找功能自动拟合识别结构面,因此相比人工识别,半自动识别方法所获得的岩体结构面数量更多且效率更高,所得到的岩体结构特征更全面。同时,结构面半自动识别分组结果获取的3组结构面数据(222°∠67°、88°∠75°和55°∠69°)也比人工识别结果获取的1组结构面产状(223°∠64°)更加全面。图8 点云数据结构面产状划分Fig.8 Structural plane occurrence division of point cloud data从结构面识别的点云数据特征来看,人工识别所选取的结构面特征更明确,对于结构面边缘的把控更为精细,半自动识别不可避免地出现点云数据的过分割或欠分割现象,因此针对点云数据结构面分割的参数优化十分重要。在Pointstudio软件基础上,针对点云数据结构面特征识别算法进行改进研究,主要利用“人工+半自动”识别方法进行结构面识别。对于巷道围岩中出露形态不规整、出露面积较小(线状出露)以及形态复杂的结构面,采用人工识别方法进行识别,而对于局部出露较规整和出露面积较大的结构面则采用半自动识别方法进行识别。按照制定的地下巷道岩体结构面特征识别技术规程,考虑到该段巷道地质结构发育的差异性以及现场管线、喷浆支护等因素对扫描数据质量的影响,将长度为120 m的巷道划分成3段研究区域,每段研究区域的长度为40 m,如图9所示。图9 研究区域划分Fig.9 Division of the study area以第一段巷道区域为例,在进行岩体结构面识别之前,将巷道点云数据划分为巷道顶板、巷道南侧边帮和巷道北侧边帮,如图10所示。由图10可见,在巷道北侧区域布满了管线,可视性较差,围岩壁面出露区域较少,为此巷道北侧区域采用Pointstudio软件进行岩体结构面识别。巷道南侧边帮靠近底板区域,围岩壁面的可视性较好,但靠近顶板区域布满电缆,影响了巷道围岩的结构面出露视线,为此巷道南侧边帮靠近底板部分可采用半自动识别程序进行结构面识别,靠近顶板部分采用Pointstudio软件进行人工结构面识别。巷道顶板区域结构面出露特征较为明显,且管线等明显遮挡物体较少,为此巷道南侧边帮靠近底板部分采用半自动识别程序进行结构面识别,同时采用Pointstudio软件进行辅助识别。图10 巷道区域划分Fig.10 Division of roadway area参照第一段巷道的划分,对三段巷道的巷道顶板、北部边帮和南部边帮进行岩体结构面“人工+半自动”识别,采用Dips软件进行结构面优势组划分,如图11所示。从识别分组结果来看,第一段巷道共识别出2组优势结构面,产状分别为240°∠15°和110°∠10°;第二段巷道共识别出2组优势结构面,产状分别为87°∠75°和260°∠66°;第三段巷道共识别出1组优势结构面,产状为189°∠83°。图11 岩体结构面识别结果Fig.11 Identification results of rock mass structural plane由图11可知,云南彝良驰宏矿业有限公司毛坪铅锌矿南部巷道的总体结构面发育程度不高,但第二段巷道所识别出的2组结构面优势组的倾向与巷道轴线方向相近,且结构面倾角较大,容易切割岩体形成可动岩块。显然,第二区段的点云数据质量最高,相比之下,第一区段和第三区段的点云数据质量较低,因而识别出来的结构面产状信息少,甚至缺失。巷道的三段区域所识别出来的优势结构面产状和数量不同,这是由每段区域的三维激光扫描点云数据质量决定的。巷道的三段区域虽然相连,但是受围岩喷浆支护和管线架设的影响,采用三维激光扫描所获取的点云数据差别较大,部分区段的点云数据质量较差,甚至空白,严重影响了岩体结构面识别作业。因此,开展井下巷道三维激光扫描作业,必须在裸露围岩上进行,以实现高质量地质编录和工程地质调查工作。在岩体工程地质调查中,结构面迹长也是一个很重要的参数,基于点云数据可以方便、快捷地开展结构面迹长估计分析工作。选取南部区域巷道扫描点云数据中长度为13.81 m的南侧巷道壁面区域作为研究区域(图12),采用“半自动+人工”点云数据结构面识别方法进行结构面提取分析。图12 岩体结构面识别区域Fig.12 Identification area of rock mass structural plane采用基于点云数据的结构面识别方法,开展巷道研究区域结构面识别工作,在所选取的研究区域巷道中共识别出28条结构面数据(表2),识别结果如图13所示。其中,表2中的“距离”为距东侧起始点的直线距离,可以看出利用本文方法可以得到结构面产状(倾角、倾向)和长度等几何特征(表中未列出结构面填充物特性、渗流情况和岩性特征等信息)。表2 结构面识别结果Table 2 Identification results of structural plane
图13 岩体结构面识别区域识别结果Fig.13 Identification results of rock mass structural plane identification region事实上,结构面的填充物厚度、张开度、粗糙度和渗流特性均可从点云数据中获取,例如:结构面的填充物厚度和张开度可通过点云数据进行尺寸长度的测量;结构面的粗糙度可通过曲率检测算法来计算面状结构面特征的起伏度;结构面的渗流特性可结合三位激光扫描仪器的实时摄像功能,通过点云—图像匹配来确认现场的渗流特性。此次扫描巷道现场管线和巷道喷浆支护对点云数据造成了影响。由图14可以看出,受扫描工况的影响,存在部分点云数据缺失,因此部分参数未在此次测量中体现。因此,在实际工程地质调查中,工程人员可利用三维激光扫描技术将结构面几何特征通过点云数据整合传输到计算机来减少繁琐的罗盘及现场测量过程。同时,在激光扫描过程中,重点对区域节理裂隙的填充物硬度和岩石脆性等具有代表性的物理力学特征进行分析。图14 不同区域岩体结构面识别结果Fig.14 Identification results of structural plane of rock mass in different regions表2只展现了结构面的出露特征,结构面的几何特征还需借助其他软件和算法进行描述。根据表中的结构面数据,利用Dips软件对结构面数据产状进行分析,可以看出该区域出露结构面产状为309°∠66°,Fisher系数为14.2,如图15所示。图15 结构面产状Fig.15 Occurrence of structural plane(1)利用三维激光扫描技术能够快速、精确地开展矿山巷道中岩体结构面识别工作,其扫描的高密度点云数据能够满足地下巷道所需要的精度要求,但在现场测量过程中应注意避开巷道管线、喷浆支护和积水等遮挡物。(2)基于三维激光扫描技术开展地下巷道岩体结构面测量及识别研究。针对巷道围岩中出露形态复杂的结构面,采用“人工+半自动”识别相结合的方法进行识别,所获取的结构面数据更加精确、全面,相比单一的测量方法,能够有效提高作业人员的工作效率。(3)由于三维激光扫描具有高精准度和便捷等优点,“人工+半自动”识别相结合的方法在地下巷道岩体结构面测量和识别中得到了进一步优化。但考虑到采用全自动识别方法进行巷道岩体结构面识别,其结果更为准确,在复杂巷道岩体中实现结构面全自动识别需要进一步研究。An P J,Kun F,Jiang Q Q,et al,2021.Measurement of rock joint surfaces by using smartphone structure from motion(SfM)photogrammetry[J].Sensors,56(1):167-189.
Chen Na,2018.Study on Structure Information Extraction and Deformation Monitoring of Slope Rock Mass Based on Three-Dimensional Laser Scanning [D].Wuhan:Wuhan University.
Chen S,Walske M,Davies I,2018.Rapid mapping and analysing rock mass discontinuities with 3D terrestrial laser scanning in the underground excavation[J].International Journal of Mining Science and Technology,110(3):28-35.
Dong Xiujun,Huang Runqiu,2006.Application of 3D laser scanning technology in geological survey of high and steep slope [J].Chinese Journal of Rock Mechanics and Engineering,25(Supp.2):3629-3635.
Ge Yunfeng,Xia Ding,Tang Huiming,et al,2017.Intelligent identification and information extraction of rock mass disjoint based on three-dimensional laser scanning technology [J].Chinese Journal of Rock Mechanics and Engineering,36(12):3050-3061.
Guo Dengshang,Xie Mowen,Liu Weinan,et al,2022.Extraction method and application of slope rock mass structural plane based on three-dimensional point cloud [J].Mining Research and Development,42(5):197-202.
He Bingshun,Ding Liuqian,Sun Ping,et al,2007.Application of three-dimensional laser scanning system in rock mass discontinuance identification [J].Journal of China Institute of Water Resources and Hydropower Research,5(1):45-47.
He Fulian,Zhang Guangchao,2015.Analysis and control of stability of the fractured soft rock surrounding a deep roadway[J].Rock and Soil Mechanics,36(5):1397-1406.
Hou Chaojiong,Wang Xiangyu,Bai Jianbiao,et al,2021.Basic theory and technology study of stability control for surrounding rock in deep roadway[J].Journal of China University of Mining and Technology,50(1):1-12.
Hu Wuting,Liu Chang,Dong Xiujun,et al,2023.Automatic identification method of rock slope structural surface based on 3D laser scanning [J].Journal of Chengdu University of Technology (Science & Technology Edition),50(3):351-360.
Jing Hongdi,Li Yuanhui,Zhang Zhong,et al,2015.Extraction of rock mass structural plane information based on three-dimensional laser scanning [J].Journal of Northeastern University (Natural Science Edition),36(2):280-283.
Kong D,Sarogloub C,Wu F,et al,2021.Development and application of UAV-SfM photogrammetry for quantitative characterization of rock mass discontinuities[J].Internati-onal Journal of Rock Mechanics and Mining Sciences,141(1):23-45.
Li H,Li X,Li W,et al,2019.Quantitative assessment for the rockfall hazard in a post-earthquake high rock slope using terrestrial laser scanning[J].Engineering Geology,248(2):1-9.
Li Jielin,Bai Dewei,Yang Chengye,et al,2022.Recognition and stability analysis of underground tunnel rock mass structural plane based on 3D laser scanning point cloud data[J].Gold Science and Technology,30(3):343-351.
Li Jielin,Yang Chengye,Hu Yuan,et al,2020.Application of UAV 3D laser scanning technology in underground goaf detection [J].Metal Mine,49(12):168-172.
Li Jielin,Yang Chengye,Peng Zhaozhi,et al,2021.Application of three-dimensional laser scanning technology in identification of structural plane of underground roadway rock mass [J].Gold Science and Technology,29(2):236-244.
Li Pengyu,Jiang Yue,Zong Qi,et al,2018.Measurement and calculation of stope volume of a gold mine based on three-dimensional laser scanning technology [J].Metal Mine,47(4):145-149.
Li Zebang,2018.Research on Application of Three-Dimensional Laser Scanning Technology in Mine Goaf Slope Deformation Monitoring[D].Kunming:Kunming University of Science and Technology.
Liang Yufei,Pei Xiangjun,Cui Shenghua,et al,2021.Structural characteristics of rock mass at landslide failure boundary based on ground three-dimensional laser point cloud[J].Chinese Journal of Rock Mechanics and Engineering,40(6):1209-1225.
Qiu Junling,Xia Qinglin,Yao Lingqing,et al,2012.Mine geological modeling and application based on three-dimensional laser scanning technology[J].Earth Science,37(6):1209-1216.
Xu Du,Feng Xiating,Li Shaojun,et al,2018.Technology and application of laser scanning to test tunnel deformation and rock mass discontinuity[J].Chinese Journal of Geotechnical Engineering,40(7):1337-1342.
陈娜,2018.基于三维激光扫描的边坡岩体结构信息提取和变形监测研究[D].武汉:武汉大学.
董秀军,黄润秋,2006.三维激光扫描技术在高陡边坡地质调查中的应用[J].岩石力学与工程学报,25(增2):3629-3635.
葛云峰,夏丁,唐辉明,等,2017.基于三维激光扫描技术的岩体结构面智能识别与信息提取[J].岩石力学与工程学报,36(12):3050-3061.
郭登上,谢谟文,刘卫南,等,2022.基于三维点云的边坡岩体结构面提取方法与应用[J].矿业研究与开发,42(5):197-202.
何秉顺,丁留谦,孙平,等,2007.三维激光扫描系统在岩体结构面识别中的应用[J].中国水利水电科学研究院学报,5(1):45-47.
何富连,张广超,2015.深部破碎软岩巷道围岩稳定性分析及控制[J].岩土力学,36(5):1397-1406.
侯朝炯,王襄禹,柏建彪,等,2021.深部巷道围岩稳定性控制的基本理论与技术研究[J].中国矿业大学学报,50(1):1-12.
胡武婷,刘昶,董秀军,等,2023.基于三维激光扫描的岩质边坡结构面自动识别方法[J].成都理工大学学报(自然科学版),50(3):351-360.
荆洪迪,李元辉,张忠,等,2015.基于三维激光扫描的岩体结构面信息提取[J].东北大学学报(自然科学版),36(2):280-283.
李杰林,白德威,杨承业,等,2022.基于三维激光扫描点云数据的地下巷道岩体结构面识别及稳定性分析[J].黄金科学技术,30(3):343-351.
李杰林,杨承业,胡远,等,2020.无人机三维激光扫描技术在地下采空区探测中的应用研究[J].金属矿山,49(12):168-172.
李杰林,杨承业,彭朝智,等,2021.三维激光扫描技术在地下巷道岩体结构面识别的应用[J].黄金科学技术,29(2):236-244.
李鹏宇,姜岳,宗琪,等,2018.基于三维激光扫描技术的某金矿采场体积测量与计算[J].金属矿山,47(4):145-149.
李泽邦,2018.基于三维激光扫描技术在矿山采空区边坡变形监测中的应用研究[D].昆明:昆明理工大学.
梁玉飞,裴向军,崔圣华,等,2021.基于地面三维激光点云的滑坡破坏边界岩体结构特征分析[J].岩石力学与工程学报,40(6):1209-1225.
邱俊玲,夏庆霖,姚凌青,等,2012.基于三维激光扫描技术的矿山地质建模与应用[J].地球科学,37(6):1209-1216.
许度,冯夏庭,李邵军,等,2018.激光扫描隧洞变形与岩体结构面测试技术及应用[J].岩土工程学报,40(7):1337-1342.
引用本文: 赵晓明,李想,李杰林,等.巷道岩体结构面“人工+半自动”识别方法及应用研究[J].黄金科学技术,2023,31(5):773-784.
ZHAO Xiaoming,LI Xiang,LI Jielin,et al.Research on “Manual+Semi-automatic” Identification Method and Application of Roadway Rock Mass Structural Plane[J].Gold Science and Technology,2023,31(5):773-784.
作者简介:赵晓明(1979-),男,河北石家庄人,高级工程师,从事矿山防治水工作。707473672@qq.com
通信作者:李杰林(1982-),男,湖南宁远人,副教授,从事金属矿山开采、采空区处理及矿山岩石力学等研究工作。lijielin@163.com
基金项目: 中南大学研究生校企联合创新项目“采空区三维激光扫描精细探测及安全评价研究”(编号:2022XQLH104)和金属矿山安全与健康国家重点实验室开放课题“深部高应力巷道围岩结构面与危险块体自动识别方法研究”(编号:2020-JSKS-SYS-06)联合资助