无人机LiDAR在采空区沉陷监测中的应用.pdf

第 卷第 期 年 月西安科技大学学报 郑俊良，姚顽强，蔺小虎，等 无人机 在采空区沉陷监测中的应用 西安科技大学学报，（）：，（）：收稿日期：基金项目：国家自然科学基金项目（）通信作者：郑俊良，男，河南洛阳人，工程师，博士研究生，：无人机 在采空区沉陷监测中的应用郑俊良，姚顽强，蔺小虎，张联队，张咏，张冬（西安科技大学 测绘科学与技术学院，陕西 西安 ；陕西彬长矿业集团有限公司，陕西 咸阳 ；陕西延长石油靖边煤业有限公司，陕西 榆林 ）摘要：煤炭资源高强度开采对地表造成了极大破坏，严重威胁建（构）筑物安全，现有的开采沉陷监测方法缺乏面域形变监测，限制了采空区地表形变规律研究。针对传统观测方法的不足，以西部煤矿某工作面为研究对象，采用无人机 （）测量方法，进行采空区面域沉陷监测研究。针对采空区地形特征，改进渐进三角网点云滤波过程中种子点的选取方法，提高点云滤波精度；对研究区沉陷模型精度进行验证，结果表明：沉陷模型整体误差小于 ，平均误差小于 ；最后，基于沉陷模型进行数据挖掘分析，获取下沉等值线图、下沉面积统计图、下沉分区统计图等，实现采空区全方位沉陷监测。为煤矿煤炭资源开采、地表下沉监测、土地复垦和生态环境保护等提供参考。关键词：无人机 ；点云滤波；沉陷监测；沉陷模型中图分类号：文献标志码：文章编号：（）：开放科学（资源服务）标识码（）：，（，；，；，）：（），（），；，；，：；引言随着能源生产和消费革命的推进，煤炭在能源结构中的占比逐渐减少，但在一段时间内煤炭仍属于中国主体能源，具有重要战略地位 。西部煤矿地处黄土高原中上游，地形复杂多样，属于典型的西部黄土沟壑区，生态环境十分脆弱。机械化综采、垮落式顶板管理的采煤方法造成上覆岩层破裂、地表塌陷，引起地下水流失、土地沙漠化，破坏沉陷区生态环境、威胁地表建（构）筑物安全 。因此，针对西部煤矿所在的黄土沟壑区，开展可靠的矿区开采沉陷监测，准确掌握开采沉陷动态变化过程，对开采沉陷预计、地表建（构）筑物损坏评价、生态修复具有重要意义。目前，开采沉陷监测方法包括传统预埋监测桩观测（水准测量、全站仪测量、测量等）、监测、地面三维激光测量、倾斜摄影测量和机载 测量。传统预埋监测桩观测方法局限于点位观测，具有较高的点位精度（毫米级、厘米级），但不能准确描述采空区地面整体的形变过程，且西部沟壑区地形结构复杂，存在工作量大和人员无法到达的问题 ；测量能够进行面状观测，且取得较高的测量精度（毫米级），但受限于观测周期，对于地表塌陷速度快、沉陷量大的区域，容易造成失相干现象 ；地面三维激光测量能够取得面状厘米级，甚至毫米级观测精度，但由于其视野范围、点云拼接误差等原因，在超长超宽工作面监测时存在工作量大、累积误差大等问题 ；近年来，随着影像密集特征点匹配算法的成熟，无人机倾斜摄影测量技术提高了摄影测量的观测精度，但在地表植被覆盖度高、地形起伏大的区域，易造成监测结果精度低，且精度受限于像控点密度 ；机载 测量具有精度高、观测范围大的优势，但有人机载 测量存在设备昂贵、天气要求高、空域申请困难等问题；随着无人机技术发展和激光扫描设备改进，无人机 测量技术逐渐成熟，在满足厘米级测量精度的同时能够极大地降低成本、提高效率，弥补有人机载 测量的不足 。无人机 通过对目标区域低空多时相扫描，可获取高密度点云数据，经点云滤波、插值和求差计算，可获取目标区域高程变化。但在激光点云滤波、插值和求差计算过程中不可避免地存在分类错误、噪声和误差，一些学者通过改进滤波算法和插值方法提高 精度，进而改善沉陷模型精度 ，但仍存在一定问题：现有成熟点云滤波算法主要针对传统有人机载 点云数据，相比较而言，无人机 点云具有更高的点云密度和更复杂的噪声数据，如在山区应用效果较好的渐进三角网点云滤波算法（，），在无人机 点云滤波时，由于高密度和多噪声的原因，易出现种子点错误问题，导致滤波错误；无人机 获取的点云、数字高程模型（，）和沉陷模型之间的精度缺乏实测数据验证；在无人机 测量系统和环境条件相同情况下，沉陷模型精度是否优于原始点云和 精度缺乏试验对比；传统 监测成果单一，对沉陷模型进一步研究较少，未进行深入的数据分析和挖掘。研究区及数据获取 研究区概况凉水井矿位于陕西省榆林市东北部、榆神府矿区 东 部，处 于 北 纬 和 东 经 之间（图 ）。井田地面标高 ，地势西北高东南低。以凉水井煤矿 工作面为研究对象，开采煤层为 煤层，该煤层平均埋深 ，平均厚度 ，采用长臂式综采、顶板自然塌陷方法开采；其中上覆 煤层已于 年开采完成，地表于 年达 西 安科技大学学报 年第 卷第 期郑俊良，等：无人机 在采空区沉陷监测中的应用到稳定状态。?图 研究对象位置 基于 煤矿测量规程（）设置地面监测点，结合临近矿区开采沉陷角量参数和现场实际情况，共布设条监测线：走向观测线、倾向观测线和道路观测线，如图所示。地面监测点通过埋设预制观测桩，采用四等水准测量，获取监测点多时相高程，实现地表沉陷监测。在工作面停采线附近，选择面积为 的区域作为研究区，进行无人机 试验数据采集和分析。工作面回采期间共进行 次全面观测，根据开采沉陷一般规律，选取位于研究区域内的 个地面监测点作为无人机 参考数据，对无人机 监测成果进行精度分析，如图 中红色三角形所示。?图 研究区与监测线布设 数据获取无人机 型号为 ，集成了 激光测距系统和高精度 （）系统，以及地面 差分系统，扫描测程 。无人机有效飞行时间约 ，试验飞行高度 ，速度 ，激光脉冲发射频率为 ，扫描速度为 线 。为了保证多时相数据精度一致，在每次数据采集时均采用相同的飞行和航测参数。在 年 月和 年 月共进行 次无人机 测量，数据采集范围覆盖该时间段地表下沉范围，航测时间与工作面回采位置、试验区域如图 所示。?图 研究区与采空区位置关系 图 展示了井下巷道、工作面位置、采空区与地表航测范围位置关系。年 月 日和 年 月 日分别进行无人机 航测，其对应的工作面回采位置如图 蓝色竖线所示。同时航测范围覆盖研究期间采煤导致的地表下沉影响范围，其中采空区 表示在首次数据采集时，已存在的采空区；采空区 表示第 次 数据采集时产生的采空区；图中黑色三角形表示地面实测监测点。研究技术路线无人机 开采沉陷监测主要步骤包括外业航测、点云滤波、点云插值计算、差值计算、沉陷模型获取及数据挖掘分析等。同时，以地面监测桩测量数据为依据进行 点云、和沉陷模型精度验证，技术路线如图 所示。?图 研究技术路线 采用渐进三角网点云滤波算法进行点云滤波，该算法一般由局部最低点作为稀疏种子点，创建初始不规则三角网，通过激光点与三角网的距离和角度进行迭代分类。该算法适用性强、滤波效果较佳，但种子点的选择是渐进三角网滤波的关键。因此，结合矿区地表沉陷特征，对种子点的选取方法进行改进，确保种子点的准确性。点云插值计算参考文献 所述，选择最适用于黄土沟壑区的 插值算法进行计算，获取地表 。点云、和沉陷模型精度采用平均误差（，）、平均绝对误差（，）、中 误 差（，）验证，以水准测量结果作为参考值，计算方法如下（）（）（）（）槡（）式中 为基于 提取的下沉值；为地面监测桩观测下沉值；为参考点数量。点云滤波与插值计算 常用滤波方法无人机 能够获取海量地表点云数据，其中包括地面、植被、建筑物等，根据点云的几何特征、色彩、高程变化等信息提取地面点的过程即为点云滤波。点云滤波方法分为：基于坡度滤波（，）、基于渐进三角网滤 波 、形 态 学 滤 波（，）、布料模拟滤波（，）等。滤波算法采用点之间的坡度作为阈值参数进行滤波，该方法对于缓坡地形效果较好，但对于地形复杂区域，由于单一阈值的限定，滤波效果不足；滤波算法借鉴形态学图像处理方法进行点云滤波，剔除非地面点时不能有效保留地形信息，易造成地形平滑；滤波算法采用点云翻转、模拟布料贴合倒置点云的方法，使用高程阈值分离地面点和非地面点，平地、缓坡和陡坡需采用不同参数，参数确定存在困难；滤波方法由 提出，根据种子点建立初始地形模型，对待定点进行距离和角度阈值判断，确定地面点或非地面点，逐步迭代完成点云滤波，该方法对于复杂的山区、森林地形具有较好的适应性，但对种子点的依赖度较高。滤波方法具有较好的适用性，特别是对于山区、沟壑区效果较好，但其种子点的选取是滤波的关键 。特别是对于煤层开采导致地表出现阶梯裂缝区域种子点的选取，如种子点选取在裂缝或塌陷坑底部会造成局部滤波错误。因此，针对开采沉陷区域梯度变形和伪最低噪声点问题，对种子点选取方法进行改进。种子点选取及滤波针对矿区，特别是浅埋煤层、大采高开采造成的阶梯塌陷地表，无人机 渐进三角网点云滤波过程中种子点的选取至关重要。因此，针对矿区阶梯状地表，在渐进三角网滤波算法的基础上，进行种子点选取方法的改进，去除错误的种子点，提高种子点质量，计算过程如图 所示。?图 改进的 滤波流程 格网划分主要是用规则的虚拟格网覆盖全部点云数据，将点云按照平面坐标分配到相应的虚拟网格内。种子点的选取是滤波的关键，采用改进的扩展局部最小值方法进行种子点选取，原理如图 所示。传统种子点选取主要采用格网内最低点作为种子点，图 中采用传统方法选取的种子点应为 。但 相较于格网内其他点，高程明显过于小，属于最低点噪声数据。文中方法将格网内激光点按高程值由小到大进行排序，并计算各自与后续点之间的高程差值（图 ）。将高差平均值作 西 安科技大学学报 年第 卷第 期郑俊良，等：无人机 在采空区沉陷监测中的应用?图 种子点选取计算示意 为阈值。若当前点与后点的高差小于高差阈值，即 珔（，），则当前点为种子点，若大于则计算第 最低点与第 最低点高差，逐步计算，直到获取满足高差阈值要求的数据点，并将该点作为种子点。该方法能够确保选择的最低点为地面点，避免噪声点、异常点的影响，如图中噪声点 。采用文中算法选取种子点，然后采用 算法进行地面点提取，结果如图 所示。?对 年 月和 年 月点云数据分别进行点云滤波，并对滤波结果进行统计，地面点云占比分别为 和 ，由于不同时期地表植被覆盖度不同，导致了地面点云占比不同。点云插值计算插值算法优劣直接影响地面 精度，目前常用的插值算法包括反距离权重（，）、克里金（）、自然邻域（，）、样条函数（）、不规则三角网（，）等 。插值以采样点与插值点间的距离为权重计算加权平均值，离插值点越近的采样点权重越大，超过一定距离时权重忽略不计，适用于采样点均匀且密集的区域 ；插值与 类似，该方法采用半变异函数作为权重进行无偏最优估计，不仅考虑距离权重，还兼顾空间分布关系，适用于空间相关性较好的采样点，但容易出现边缘效应 ；自然邻域插值是基于 结构的插值计算，该方法对局部特征的继承性较好，但局部点缺失时容易造成失真 ；以最小曲率面逼近各采样点，以获取最佳平滑曲面，顾及大范围的采样点，适用于复杂高曲率曲面，但平滑易造成地形失真 ；不规则三角网（）利用采样点构建不规则三角形，能够在地形复杂区域保留地形特征细节，但其边缘梯度可能存在突然变化 。文献 的研究成果表明，在黄土高原沟壑区，插值在无人机 点云插值计算中具有较好效果。因此，采用 插值算法对地面点云进行插值计算，并将插值格网大小对结果精度的影响进行试验。采用 误差分布和误差值大小对比的方法进行结果分析，不同格网大小引起的 误差分布如图 所示。?图 误差分布 同时，计算不同格网尺寸 的 ，结果见表 。由表 可知，和 在 区间内变化不大，当格网尺寸大于 时，误差急剧增大。同时，结合图 综合考虑，最终选择 的格网尺寸进行插值计算，获取较高精度的研究区 。差值计算沉陷模型是由多时相地面 差值计算获得，与 结构和功能类似，能够反映一段时间内地表高程变化，即地表沉陷量，是一种高精度地图 局部点云滤波效果 表沉陷表示模型 。通过不同时相 之间的差值计算，可获取该段时间内地表的下沉变化。以首次观测时获取的 为基础，与开采结束后获取的 做减法操作，可得到该段时间内地表的沉陷模型，其三维效果如图 所示，由于地表人为破坏、建筑施工等，导致沉陷量部分异常值过大或过小。表 误差 格网大小 平均误差 平均绝对误差 中误差?图 沉陷模型三维效果 精度评定基于地面监测桩实测数据对地面点云、和沉陷模型进行精度评定。由于点云由一系列离散点构成，无法进行点对点差值直接计算，因此采用局部平均值的方法进行差值计算，计算监测桩附近局部点云的高程平均值，然后与实测数据对比，结果见表 ，年 月 测量点云中误差为 ，年 月为 。表 点云误差统计 日期最大 最小 平均 中值 中误差 其次，基于监测桩实测数据对 和沉陷模型误差进行计算，误差分布如图 所示。?图 点云、，沉陷模型误差分布 ，从图 可以看出，原始地面点云、和沉陷模型三者中 的精度较差，沉陷模型精度最优。在无人机 航测过程中，多时相数据采集时使用相同航测参数，能够消除一定的系统误差，如激光扫描设备系统误差、飞行高度、地形变化等因素的影响，提高沉陷模型的精度。基于沉陷模型提取地面监测桩的下沉量，与地表监测点监测值进行对比，绘制下沉折线图，如图 所示。在 年 月至 年 月期间，地面水准测量获取的最大下沉值为 ，最小下沉值为 ，如图 中折线所示；基于沉陷模型获取的最大下沉值为 ；除 监测点实测数据下沉异常，实测下沉折线图与基于沉陷模型获取的下沉折线图走势和拐点基本一致。相同位置实测下沉量与沉陷模型下沉量的误差绝对值最大为 ，最小为 ，平均为 ，中位数为 ，其中小于 的占比 ，误差绝对值整体小于 ，误差分布相对集中。西 安科技大学学报 年第 卷?图 下沉折线 第 期郑俊良，等：无人机 在采空区沉陷监测中的应用 基于沉陷模型的数据挖掘与分析沉陷模型包含丰富的地表下沉信息，包括地表平面位置、任意位置下沉量、下沉变化等，基于沉陷模型进行数据挖掘分析，能够获取开采沉陷引起的下沉范围、边界角、下沉面积、不同程度损害区域面积等。下沉量统计为了准确获取地表不同位置的下沉量和不同下沉量的分布情况，对沉陷模型进行下沉量分级分类统计，结果如图 所示。?图 下沉分布 图 表示不同下沉量区域分布情况，下沉区域呈椭圆形分布，最大下沉区位于采空区中心，由中心逐渐向外扩散，下沉值逐渐减小，相应的下沉面积逐渐增大。图中圆形标注区域出现异常下沉，经实地调查发现是由于煤层开采期间，人为施工造成。对下 沉 区 域 面 积 进 行 分 级 计 算，以 大 于 的下沉值作为开采沉陷引起的地表下沉，计算结果见表 。下沉值大于 的下沉总面积为 ，采空区面积为 。为准确确定开采沉陷引起的不同下沉量所占面积，暂定以下沉值大于 的沉陷量作为煤层开采引起的地表下沉。结果表明，下沉值小于 的面积占比超过一半以上，随着下沉量的增大，对应的下沉面积逐渐减小。下沉剖面信息提取传统开采沉陷监测采用预埋监测桩的方式，监测煤层开采过程中地表变化，该方法获取的监测数据少、范围小。试验基于沉陷模型采用虚拟监测点和剖面线的形式，分别对沉陷区域进行走向和倾向方向下沉信息提取，虚拟监测点与剖面线位置如图 所示。表 下沉面积统计 下沉值 像元数 个像元尺寸 面积 总面积比例 开采沉陷面积比例?虚拟监测点以 间隔在走向和倾向方向均匀分布，如图 中黑色三角形所示；剖面线分为走向剖面线和倾向剖面线，位置如图 中灰色实线所示。基于沉陷模型在虚拟监测点位提取地表下沉值，并进行折线图绘制，结果如图 所示。虚拟监测点最大下沉值为 ，走向下沉折线呈半盆地分布，倾向下沉折线图呈全盆地分布，符合开采沉陷下沉规律。根据沉陷模型分辨率大小，下沉剖面线能够提取大量监测点数据，相比离散监测点具有更好的鲁棒性。因此，在沉陷区域设置一条长 的走向剖面线和长 的倾向剖面线，基于剖面线以沉陷模型分辨率 为间隔，提取下沉值。对提取结果进行散点绘制，并提取中心线，结果如图 所示。基于沉陷模型提取的走向下沉曲线呈盆地形状，说明在走向方向达到超充分采动，符合开采沉陷一般规律；倾向方向下沉曲线呈山谷对称分布，图 虚拟监测点与剖面线位置 其中存在部分异常值，经实地探查，造成该异常的主要原因是由于地表人为开挖施工，去除异常值，整体符合开采沉陷地表下沉特征。?图 虚拟监测点下沉折线?沉陷面域分析传统开采沉陷监测主要对固定点进行监测，基于沉陷模型不仅能够进行传统点监测、剖面线监测，获取下沉值和下沉折线图。同时也能够进行面域监测分析，进一步提取下沉等值线、下沉范围、地表损害区域图等。等值线是一种由相同值连接构成的曲线，如等高线、等温线、等深线等。文中采用等值线表示地表下沉值，能够清晰地表示地表下沉情况。等值线的精度一般取决于等值线的间隔值，间隔值越小，图形越精细、等值线越密、表达的效果越好。但是等值线间隔值越小，越容易受异常值和离散数的影响，图形越复杂、异常数据影响越明显。因此，合适的等值线间隔能够完美表达地表下沉情况，同时能够降低异常值的影响。基于工作面下沉值大小和不同下沉值的分布区间，采用多组、分梯度的方法设置等值线间隔，分别选用 ，和 间隔值进行下沉等值线绘制，结果如图 所示。从图 可以看出，不同等值距的下沉等值线所表达的特征详细度不同。较小的等值线区间值对应较详细的沉陷特征，同时也易受噪声和异常值的影 响。通 过 综 合 比 较 发 现，沉 降 间 隔 为 的等值线既可以保留采煤沉陷特征细节，也能够降低噪声和异常值的影响。传统开采沉陷以边界角确定煤层开采影响范围，但由于西部煤矿地形复杂，不同地形下沉特征有所区别，走向和倾向观测线获取的边界角难以准确描述所有地形条件下的影响边界。因此，基于沉陷模型进行开采边界的提取能够弥补传统测线的不足。但由于地表自然变化（如雨水冲刷、风沙侵蚀等）和无人机 测量精度的影响，获取 下沉边界存在一定困难。因此，考虑自然因素影响，以 下沉值作为煤层开采引起的下沉边界，提取地表下沉边界，如图 所示。从图 可以看出，采空区 开采引起的沉陷区域面积 ，小于表 中计算的下沉面积 ，造成该结果的原因是由于研究区内存在人工开挖区域和远离采空区的下沉面积，该面积并不是由于煤层开采造成的地表下沉。下沉区域整体偏向采空区 ，造成该结果的原因是在 年 月无人机 航测前，采空区 已存在，其对地表造成的影响并未结束，因此后续地表下沉范围偏向采空区 。无人机 监测期间，地表下沉面积与采空区面积比例为 ，该 西 安科技大学学报 年第 卷?图 剖面下沉曲线 第 期郑俊良，等：无人机 在采空区沉陷监测中的应用比例能够为后续煤层开采地表影响范围的确定提供参考。?图 不同等值距的下沉等值线?图 下沉边界（）（）为了明晰煤层开采造成地表下沉的范围和严重程度，以不同下沉量所占下沉面积为参数，进行地表损害的分等定级，将地表下沉值进行分组，并确定其范围和对应面积，结果如图 所示，可以看出，不同沉降等级区域由内向外逐渐扩展，采空区从中心向外延伸。沉降可分为 级：，和 以上，相应的沉降区面积为 ，和 。沉降面积之比为 ，能为工作面开采过程中地表建（构）筑物的保护、移民搬迁提供依据。?图 地表下沉量分级 结论）在煤矿地形复杂区域，采用无人机 进行煤层开采沉陷监测，针对沉陷区地形特征，对 滤波算法种子点的选取方法进行改进，剔除了最低点噪声，实现无人机 地面点云的稳健提取。采用 插值算法对地面点云进行内插生成 ，对格网尺寸进行试验，确定最优的格网尺寸大小（），保证 误差分布和精度均处于最优，获取较高精度的地面 和沉陷模型。）无人机 测量的高程精度可达 ，受地面噪声、飞行条件、地形地貌的影响，的精度可达 ，由于相似误差的抵偿，沉陷模型精度可达 。因此，研究区基于无人机 获取的沉陷模型整体精度优于 ，平均误差小于 。）根据沉陷模型下沉特征进行多方位分析，对沉陷模型进行挖掘和分析，获取地表下沉折线图、下沉等值线图、地表损害分布图、地表下沉面积与采空区比例关系等数据，能够为开采沉陷预计和地表建构筑物保护提供依据。参考文献（）：王双明，段中会，马丽，等 西部煤炭绿色开发地质保障技术研究现状与发展趋势 煤炭科学技术，（）：，（）：王冠，刘静 煤炭经济发展的绿色向度与绿色道路 人民论坛，（）：，（）：芦家欣，汤伏全，赵军仪，等 黄土矿区开采沉陷与地表损害研究述评 西安科技大学学报，（）：，（）：余学义，王昭舜，杨云 大采深综放开采地表移动变形规律 西安科技大学学报，（）：，（）：许家林，钱鸣高 绿色开采的理念与技术框架 科技导报，（）：，（）：张海君，吴奕枢，王飞，等 准格尔煤田龙王沟煤矿特厚煤层综放开采地表沉陷规律 西安科技大学学报，（）：，（）：汤伏全，董龙凯，王宗良，等 基于沉陷对称特征的近水平煤层开采 三维位移反演模型 煤炭学报，（）：，（）：李海滨，唐国茜，赵桂娟，等 三维激光扫描技术在路堤沉降观测中的应用 西安科技大学学报，（）：，（）：，（）：，（）：詹容若，段亮，罗晓容，等 无人机多点位航拍高分辨率三维数字露头建模 西安科技大学学报，（）：，（）：何柯璐，汤伏全，韦书平，等 机载激光点云构建西部矿区开采沉陷模型研究 测绘科学，（）：，（）：，西 安科技大学学报 年第 卷第 期郑俊良，等：无人机 在采空区沉陷监测中的应用 ，（）：，（）：王竞雪，张雪洋，洪绍轩，等 结合形态学和 三角网的机载 点云滤波算法 测绘科学，（）：，（）：，甘斌，郑俊良，姚顽强，等 基于数学形态学的机载 采煤区沉陷信息提取 西安科技大学学报，（）：，（）：，（）：李朋飞，张晓晨，严露，等 复杂地形中机载 点云构建 的插值算法对比 农业工程学报，（）：，（）：，（；）：隋立春，张熠斌，张硕，等 基于渐进三角网的机载 点云数据滤波 武汉大学学报（信息科学版），（）：，（）：，（；）：，（）：李卉，李德仁，黄先锋，等 一种渐进加密三角网 点云滤波的改进算法 测绘科学，（）：，（）：，张孟丹，余钟波，谷黄河，等 无定河流域降水量空间插值方法比较研究 人民黄河，（）：，（）：，：，（）：冯波，陈明涛，岳冬冬，等 基于两种插值算法的三维地质建模对比 吉林大学学报（地球科学版），（）：，（），（）：（责任编辑：高佳）