Laser speckle stripes supersaturated problems in coal flow measurement research
-
摘要:
为解决煤流量三维扫描过程中亮度过饱和激光条纹中心线提取精度的问题,提出了一种对应的中心线提取算法。根据斜射式激光三角法测量原理,搭建了视觉测量系统平台,获取亮度过饱和的煤堆表面激光条纹数字图像并进行预处理。探究激光光强的分布规律,深度剖析激光光条亮度产生过饱和的原因。在分析传统结构光中心线提取方法不适用于本研究的基础之上,针对部分过饱和提出了基于高斯曲线拟合法的中心点提取算法,而针对全部过饱和则提出了一种基于邻域块思想的中心点提取算法。最后确定以激光条纹截面上过饱和采样点数量作为判别依据的中心线提取策略,得到亚像素级坐标精度的条纹中心线。通过设置与灰度重心法、曲线拟合法的对照模拟试验,发现该算法绝对偏差均值更小,提取精度要优于传统算法。最终在煤堆激光条纹图像上开展了试验验证,表明该算法在实际工作环境下抗噪声干扰性能力更强,减少了亮度过饱和现象给煤堆表面激光条纹中心线提取所带来的影响。
Abstract:In order to solve the problem of extraction accuracy of brightness supersaturated laser stripe centerline in three-dimensional scanning of coal flow, a centerline extraction algorithm is proposed. According to the principle of oblique laser triangulation measurement, a visual measurement system platform was built to obtain the digital image of laser stripes with oversaturated brightness on the surface of coal piles and preprocess it. The distribution law of laser light intensity is explored, and the reason for the supersaturation of laser stripe brightness is analyzed. It is found that the traditional structured light centerline extraction method is not suitable for this study. Therefore, a center point extraction algorithm based on Gaussian curve fitting method is used for partial supersaturation, and a center point extraction algorithm based on neighborhood block idea is applied for total supersaturation. The centerline extraction strategy based on the number of supersaturated sampling points on the laser stripe cross section is proposed to obtain the fringe centerline with sub-pixel coordinate accuracy. By setting the contrast simulation experiment with the gray center of gravity method and curve fitting method, it is found that the mean absolute deviation of the proposed algorithm is smaller, and the extraction accuracy is better than the traditional algorithm. Finally, the experimental verification is carried out on the laser stripe image of coal heap, and the experimental results indicate that the proposed algorithm has stronger anti-noise interference ability in the actual working environment, which can reduce the influence of brightness supersaturation on the extraction of laser stripe centerline on the surface of coal heap.
-
Keywords:
- coal flow scanning /
- image processing /
- centerline extraction /
- extraction algorithm
-
0. 引 言
随着我国浅部煤炭资源日趋枯竭,大部分煤矿逐步转入深部开采,造成煤层瓦斯含量增大,极易诱发工作面瓦斯超限,从而制约煤矿安全生产[1]。通过利用采动对煤层瓦斯的卸压增透效应,构建顶板瓦斯抽采巷(简称“顶板巷”)抽采卸压瓦斯是治理工作面瓦斯超限的有效途径之一[2]。顶板巷因其大流量和连续抽采的优点被广泛应用,尤其在山西地区[3]。然而,合理确定顶板巷布置位置是保障顶板巷瓦斯抽采效果的关键。因此,开展顶板巷精准布置方法研究,对防治煤矿瓦斯灾害与开发煤层瓦斯资源具有重要意义。
目前,顶板巷布置位置主要根据采动覆岩裂隙发育的空间形态来确定[3-4],常用方法有现场考察和经验公式理论计算,但上述两种方法均需耗费较多时力、成本较高。近年来,基于GA–BP神经网络的回归预测法因具有考虑因素较为全面、低成本、高效率和高精度的特点,在矿业领域的数据预测中得到了广泛应用[5-6],这为提出一种利用GA–BP神经网络模型来实现顶板巷布置位置智能预测的新方法提供了可能。
因此,通过深入分析顶板巷的布置原则及其影响因素,利用山西地区顶板巷的实际布置参数作为数据样本,从数据驱动角度运用GA–BP神经网络模型进行回归预测,形成了一种基于GA–BP神经网络模型的顶板巷布置位置智能预测方法,并设计开发了相应的智能预测系统,可为顶板巷精准布置提供参考。
1. 顶板瓦斯抽采巷布置原则及影响因素
1.1 顶板瓦斯抽采巷的布置原则
在煤层开采过程中,煤层采动引起上覆岩层破断、移动并产生大量裂隙,为煤层瓦斯的运移与持续解吸创造了有利条件。因此,掌握覆岩采动裂隙空间形态的发育特征是确定顶板巷布置位置的关键。
针对上述问题,袁亮等[7]提出了采动覆岩“高位环形裂隙体”理论,并将其应用于顶板巷布置位置的确定。李树刚、林海飞[8-10]等建立了采动裂隙椭抛带动态演化模型。钱鸣高院士[11]揭示了长壁工作面覆岩采动裂隙的两阶段发展规律与“O”形圈分布特性。近年来,笔者研究团队[2-4,11-16]提出了采动覆岩瓦斯卸压运移“三带”理论以及采动裂隙“O”形圈理论,建立了基于关键层运动的顶板巷布置优化设计方法。上述研究成果为确定顶板巷布置位置提供了基本原则。
根据采动覆岩瓦斯卸压运移“三带”理论,采空区上覆煤岩层可划分为导气裂隙带、卸压解吸带和不易解吸带,其中导气裂隙带中煤岩层裂隙发育程度高且数量多,是煤岩层卸压瓦斯的主要贮存空间和运移通道[12];采动裂隙“O”形圈理论指出,随着工作面的不断推进,离层裂隙增加,采空区中部裂隙先增加后减少然后逐渐被压实,采空区两侧离层裂隙仍能保持,从平面看,在采空区四周存在如图1所示的沿层面横向连通的离层发育区,称之为采动裂隙“O”形圈,“O”形圈裂隙区煤岩体透气性较好,是采场卸压瓦斯流动的主要通道[11]。
由此可见,在垂直方向上,为防止抽采巷过早破坏、保证抽采效率,顶板巷应布置在裂隙充分发育且富集瓦斯的导气裂隙带中上部;在水平方向上,回风巷侧受工作面风流影响,瓦斯体积分数较高,顶板巷应靠近回风巷,且远离采空区中部的重新压实区,并布置在“O”型圈裂隙中[11,15]。顶板巷布置原则如图2所示。
1.2 顶板瓦斯抽采巷布置的影响因素
根据上述的顶板巷布置原则,其布置位置主要由采动覆岩的导气裂隙带和“O”形圈的范围所决定。大量研究和现场实践表明[11-23],采动覆岩导气裂隙带和“O”形圈的范围受到多种因素影响,包括工作面的采厚、倾向长度、覆岩结构、埋深、煤层倾角、推进速度等,且与上述影响因素之间存在着非线性关系,具体见表1。
表 1 覆岩采动裂隙空间形态的影响因素Table 1. Factors influencing spatial morphology of mining-induced fractures影响因素 影响规律 采厚 覆岩导气裂隙带高度和“O”形圈裂隙区宽度随工作面采厚的增加而增大,且“O”形圈宽度变化整体近似呈对数分布[12-13] 倾向长度 在一定的覆岩结构条件下,工作面倾向长度越大,覆岩“O”形圈裂隙区宽度越大[11,15] 覆岩结构 在其他条件相同时,覆岩岩性越坚硬,覆岩的采动裂隙发育越显著[12,14,19] 埋深 随着工作面埋深的增加,顶板岩层的地应力增大,工作面覆岩裂隙越发育[6] 煤层倾角 当煤层倾角小于45°时,覆岩裂隙带高度随着煤层倾角增大而增大;当煤层倾角为45°~60°时,
覆岩裂隙带高度随倾角的增大而减小[23]推进速度 较慢的推进速度增加覆岩的变形量和裂隙分布,加剧覆岩的破坏程度;而较快的推进速度则可以减小
覆岩变形量和裂隙发育,有效抑制覆岩裂隙发育[20]2. 基于GA–BP神经网络的顶板瓦斯抽采巷布置位置智能预测模型
2.1 GA–BP神经网络模型
2.1.1 BP神经网络
BP神经网络是一种使用误差反向传播算法训练的多层前馈网络,是应用最广泛的神经网络模型之一,主要由输入层、隐藏层和输出层构成[23],单隐藏层神经网络数学表达式如式(1)和式(2)所示。BP神经网络包括数据的前向传播、误差计算、误差反向传播、权重更新和迭代等5个过程[24],当预测输出与实际输出之间的误差未满足终止条件或者满意的训练效果,则将误差通过网络反向传播,进行权重更新,通过反复迭代,使神经网络的预测输出达到满意的训练效果。
$$ {\textit{z}}_{j}^{h}=f\left(\sum_{i=1}^{m} v_{i j} x_{i}+b_{j}^{h}\right) $$ (1) $$ y_{k}=f\left(\sum_{j=1}^{q} w_{j k} {\textit{z}}_{j}^{h}+b_{k}\right) $$ (2) 式中:$ {\textit{z}}_{j}^{h} $为隐藏层第j个神经元的输出;$ y^{k} $为第k个输出层神经元的输出;$ x_{i} $为第i个输入层神经元的输出;$ v_{ij} $为输入层到隐藏层的权重;$ w_{j k} $为隐藏层到输出层的权重;$ b_{j}^{h} $为隐藏层的偏置项;$ b_{k} $为输出层的偏置项;m为输入层神经元的个数;q为隐藏层神经元的个数。
2.1.2 遗传算法
遗传算法(Genetic Algorithms)是1962年由美国Michigan大学Holand教授提出的模拟自然界遗传机制和生物进化论而形成的一种启发式搜索和优化方法[25]。该算法借鉴了自然生物体中基因的传承与种群的适应性进化过程,以“适者生存”为核心原则,旨在追求问题的最优解。通过定义适应度函数,并依此函数在种群中执行选择、交叉和变异操作,算法能够筛选出表现良好的个体进行保留,同时淘汰那些适应度较低的个体。新一代的种群在继承前一代优秀特征的基础上,会表现出更加优越的适应性。这一过程将不断重复迭代,直至达到既定的停止条件,从而完成优化任务,其基本步骤如图3所示。
2.1.3 GA–BP神经网络
GA–BP神经网络是利用遗传算法优化BP神经网络的初始权重和偏置,以改善BP神经网络的收敛速度慢和可能陷入局部最优解的问题[26-27]。GA–BP神经网络分为初始化神经网络、遗传算法优化初始权重和偏置和训练BP神经网络3个部分。其中,初始化神经网络包括根据输入层和输出层神经元节点个数确定BP神经网络结构,进而确定遗传算法种群的大小即不同的权重和偏置组合的数量;使用遗传算法优化BP神经网络的初始权值和偏置,种群中的每个个体都包含了一个网络所有的权值和偏置,个体通过适应度函数计算个体的适应度值,遗传算法通过选择、交叉和变异操作找到最优适应度值对应个体;BP神经网络用遗传算法得到最优个体对网络初始权值和偏置赋值,使用收集的数据进行训练和预测。
GA–BP神经网络预测流程如图4所示。
2.2 模型输入指标
由1.2节可知,影响顶板巷布置的物理指标有采厚、埋深、推进速度、覆岩结构、倾向长度、煤层倾角。其中覆岩结构参照文献[19]用覆岩硬岩岩性比例系数来量化,计算在工作面顶板导气裂隙带高度范围内各硬岩地层(细砂岩、中砂岩、粗砂岩和石灰岩)的总厚度和导气裂隙带高度的比值。硬岩岩性比例系数计算公式为
$$ Q=\frac{\sum k}{(15 \sim 20) h} $$ (3) 式中:Q为硬岩岩性比例系数:h为采厚;Σk为15~20倍采厚中硬岩岩层的累计厚度。
以表2寺家庄15303工作面覆岩结构为例,覆岩结构量化值Q=$ \dfrac{42.31}{88.17} $=0.48。
表 2 寺家庄15303工作面覆岩结构Table 2. Overlying strata structure of No. 15303 working face in Sijia Village岩层岩性 厚度/m 距15号煤累厚/m 81号煤 0.90 89.61 砂质泥岩 1.00 88.61 82号煤 0.44 88.17 粉砂岩 4.75 83.42 砂质泥岩 1.39 82.03 粉砂岩 5.80 76.23 泥岩 1.40 74.83 84号煤 1.20 73.63 泥岩 2.30 71.33 9号煤 1.30 70.03 砂质泥岩 1.60 68.43 9号下煤 0.60 67.83 粉砂岩 1.50 66.33 中砂岩 3.73 62.60 细砂岩 1.00 61.60 粗砂岩 6.05 55.55 石灰岩 1.50 54.05 砂质泥岩 0.30 53.75 11号煤 0.30 53.45 砂纸泥岩 2.90 50.55 12号煤 0.90 49.65 粉砂岩 3.90 45.75 细砂岩 3.80 41.95 中砂岩 3.00 38.95 泥岩 0.93 38.02 石灰岩 4.00 34.02 13号煤 0.60 33.42 细砂岩 7.15 26.27 砂质泥岩 3.30 22.97 泥岩 2.00 20.97 石灰岩 4.70 16.27 14号煤 0.23 16.04 泥岩 1.00 15.04 粉砂岩 3.20 11.84 14号下煤 0.26 11.58 砂质泥岩 2.20 9.38 中砂岩 3.56 5.82 细砂岩 3.82 2.00 砂质泥岩 2.00 0 15号煤 5.13 — 砂纸泥岩 4.00 0 细砂岩 6.00 4.00 顶板巷布置受上述众多物理指标的影响,为获得主要影响指标,将所述的7个指标作为候选集,结合灰色关联度分析,得到各指标与顶板巷垂直层位和水平位置的关联度$ \xi_{i}(j) $和权重Wi,根据权重大小绝对模型输入指标。
灰色关联度分析[28]是一种用于研究变量之间关联程度的方法,适用于研究变量之间的关联性,特别是在样本数据较少或者数据分布不规则的情况下。
其关联系数$ \xi_{i}(j) $计算公式为
$$ \begin{split} \xi_{i}(j)=&\frac{\min _{i} \min _{j}\left|x_{i}(j)-y(j)\right|}{\left|x_{i}(j)-y(j)\right|+\rho \max _{i} \max _{j}\left|x_{i}(j)-y(j)\right|} +\\ &\frac{\rho \max _{i} \max _{j}\left|x_{i}(j)-y(j)\right|}{\left|x_{i}(j)-y(j)\right|+\rho \max _{i} \max _{j}\left|x_{i}(j)-y(j)\right|} \end{split} $$ (4) 式中:$ x_{i} $为第i个序列;$ y $为参考序列;j为序列中的第j个数值;$ \rho $为分辨系数,本文取0.5。
关联度$ \gamma_{t} $计算公式为
$$ \gamma_{i}=\frac{1}{m} \sum_{j=1}^{m} \xi_{i}(j) $$ (5) 根据关联度可得权重Wi为
$$ W_{i}=\frac{\gamma_{i}}{\displaystyle\sum_{i=1}^{n} \gamma_{i}} $$ (6) 式中,Wi为第i个指标的权重。
为了消除输入数据指标之间的量纲影响,采用线性归一化法对数据进行归一化处理,使各指标数据转化为[0,1]之间的数。转换公式为
$$ Y=\frac{X-X_{\min }}{X_{\max }-X_{\min }} $$ (7) 式中,Y为归一化处理后的数据;X为输入数据; Xmin为输入的样本数据中的最小值;Xmax为输入的样本数据中的最大值。
从大量相关文献和煤矿企业实际工作面的生产数据中收集了山西地区各个工作面的顶板巷布置位置的实际工况数据,组建了GA–BP神经网络模型的训练集,见表3。
表 3 顶板巷布置位置训练集Table 3. Roof gas drainage roadway location training set序号 工作面 采厚/m 覆岩结构 倾向长度/m 倾角/(°) 埋深/m 日推进距离/m 垂直层位/m 伸入工作面
水平距离/m抽采效果 1 五阳煤矿7603[29,30] 6.10 0.55 180 8.5 480 5.6 35 40 抽采量22 m3/min 2 赵庄矿1310[31] 4.88 0.33 219 10 510 3.6 30 20 抽采浓度29% 3 赵庄矿1307[32,33] 4.50 0.21 220 3 568 4.2 47 34 抽采量30 m3/min 4 漳村煤矿2601[34] 6.00 0.56 225 5 550 4.8 22 21 抽采量8 m3/min 5 寺家庄矿15106[35] 5.40 0.32 286.2 4 480 6.4 30 40 抽采量29 m3/min 6 新大地15201[3] 5.20 0.32 180 10 406 3.0 46 54 抽采率92% 7 李村煤矿1303[36] 5.10 0.40 280 5 482 5.2 45 35 抽采浓度18% 8 阳泉三矿K8206[10] 6.80 0.55 252.2 5 547 3.4 60 50 抽采率87% 9 阳煤三矿K8108[37] 6.35 0.33 190.3 7 580 3.0 70 50 抽采量59 m3/min 10 阳煤三矿K8110[38] 6.42 0.48 190.3 5 527 3.5 80 63 抽采量49 m3/min 11 阳煤一矿S8101[39] 6.02 0.28 240 3 434.5 4.2 67 51 抽采率66% 12 石港煤矿15101[40] 7.02 0.18 152 9 470 2.5 54 60 抽采率70% 13 开元9404[10] 4.23 0.59 180 5 365 3.0 54 28 抽采率77% 14 阳泉一矿S8310[41] 6.51 0.74 220 6 602 5.6 70 50 抽采量28 m3/min 15 开元9801[42] 5.58 0.32 180 8 440 3.6 43 40 抽采浓度45% 16 五阳煤矿7607[43] 6.06 0.52 227 5 450 3.6 35 41 抽采率55% 17 白羊岭煤矿15118[44] 4.60 0.26 240 9 390 3.2 52 45 抽采浓度21% 18 阳煤五矿83206[45] 6.40 0.61 191.5 8 520 2.9 57 45 抽采率90% 19 漳村煤矿2603[46] 5.85 0.46 240 5 538 4.8 20 15 抽采量8 m3/min 20 马兰矿18305高瓦斯工作面[47] 4.20 0.71 228 3 468 6.0 40 25 抽采量16 m3/min 将表3中的数据通过式(7)进行线性归一化处理,依次将顶板巷垂直层位和水平位置作为参考序列;在此基础上,由式(4)和(5)计算出影响顶板巷布置的各个物理指标与参考序列的关联度,而后结合式(6)确定各个物理指标对顶板巷垂直层位和水平位置的影响权重,进而根据各个物理指标的权重之和得到各指标对顶板巷布置的影响程度,其排序结果见表4。
表 4 影响指标权重Table 4. Influence weight of indicators影响指标 垂直层位关联度 水平位置关联度 垂直层位影响权重 水平位置影响权重 总权重 综合排名 采厚 0.652 0.744 0.177 0.189 0.366 1 覆岩结构 0.624 0.624 0.170 0.159 0.329 3 倾向长度 0.598 0.647 0.162 0.164 0.326 5 倾角 0.584 0.666 0.159 0.169 0.328 4 埋深 0.651 0.668 0.177 0.170 0.347 2 推进速度 0.572 0.585 0.155 0.149 0.304 6 由此可知,各物理指标对顶板巷布置位置的综合影响程度由大到小依次为:采厚、埋深、覆岩结构、倾角、倾向长度、推进速度。将总权重大于0.32的采厚、埋深、覆岩结构、倾角、倾向长度等物理指标作为山西地区顶板巷布置位置的主要影响指标,并选取上述5个指标作为GA–BP神经网络输入指标。
2.3 GA–BP神经网络参数的确定
GA–BP神经网络参数由神经网络参数和遗传算法参数2部分组成。在神经网络参数取值方面,选择合适的BP神经网络的隐藏层层数和神经元个数可以提高模型预测的准确性;Kolmogorov定理证明,单隐藏层结构能够以任意精度映射非线性输入输出关系[28],因此本文选择隐藏层层数为1,隐藏层节点数通过式(8)计算,进而借助试错法细化,以达到预测误差最小化[26];在激活函数的选择上,隐藏层采用ReLU函数,以引入非线性并提高训练效率,而输出层则使用线性激活函数,以直接输出预测值,模型性能的评估以均方误差(MSE)为标准。同时,在构建GA–BP神经网络模型过程中,遗传算法编码使用“浮点数编码方法”,选择算法使用“轮盘赌选择法”,交叉算法使用“模拟二进制交叉”,变异算法使用“高斯变异”。对于上述遗传算法参数的取值,通过选取各参数不同数值进行组合进行多次训练,比较模型收敛速度和预测准确率,选择收敛速度快和准确率高的参数组合,具体见表5。
表 5 GA–BP神经网络参数取值Table 5. Parameter setting of GA–BP neural network遗传算法参数 取值 神经网络参数 取值 种群大小 60 输入层节点数 5 遗传代数 70 输出层节点数 2 交叉概率 0.7 隐藏层层数 1 变异概率 0.2 隐藏层节点数 8 $$ b=\sqrt{m+n}+a $$ (8) 式中,m为输入层神经元个数;n为输出层神经元个数;a为[1, 10]的常数。
以2.1节选择出的采厚、覆岩结构、倾向长度、倾角和埋深等5个指标作为GA–BP神经网络的输入,垂直层位(离煤层顶板距离)和水平位置(离回风巷的水平投影距离)作为预测模型的输出,用表3的20组数据作为GA–BP神经网络模型的训练样本。
GA–BP神经网络适应度变化如图5所示,由图5中可以看到随着遗传代数的不断增加,适应度不断减小,遗传代数到60代之后,变化幅度很小,基本维持稳定,因适应度选择用MSE(均方误差)表示,所以适应度越小,表明GA–BP神经网络初始权重和偏置效果越好。
2.4 模型预测结果及分析
利用表6中的3组数据对训练之后的GA–BP神经网络模型进行检验,GA–BP神经网络模型的测试结果见表6。
表 6 模型预测结果Table 6. Model prediction results工作面 采厚/m 硬岩岩性
比例系数倾向长度/m 倾角/(°) 埋深/m 垂直层位/m 垂直层位
相对误差/%伸入工作面水平距离/m 伸入工作面
水平距离/%预测 实际 绝对误差 预测 实际 绝对误差 寺家庄15110 5.67 0.35 180.0 5 327 48.6 53 4.4 0.83 46 48 2 4.17 阳煤K8205 6.90 0.62 260.8 8 580 48.8 55 6.2 11.27 47.2 47 0.2 0.43 阳煤五矿8204 6.50 0.58 225.0 8 623 56.7 60 3.3 0.55 41.2 40 0.2 0.50 由表6的预测结果可知,顶板巷垂直层位预测结果的绝对误差最大为6.2、相对误差最大为11.27%,水平距离预测结果的绝对误差最大为2、相对误差最大为4.17%,预测结果相对可靠,可以为山西地区煤矿回采工作面顶板巷布置位置提供一定的参考,后续研究应进一步搜集更多现场实测原始数据样本,解决顶板巷布置位置预测的“小样本、贫信息”问题。
3. 顶板瓦斯抽采巷布置智能预测系统软件的设计及实现
3.1 预测系统设计
顶板瓦斯抽采巷布置预测系统包含模型训练、指标参数输入、预测结果输出等3个模块,其运行原理如图6所示。
首先,收集模型训练所需的相关数据,将初步处理的数据存放到文件中,并采用GA–BP神经网络算法训练预测模型。而后,用户在预测系统交互界面中输入指标参数,并通过加载训练好的预测模型,计算得到预测结果,并显示在结果输出界面上。预测系统的交互界面设计本着便捷和直观的原则,方便用户输入参数和查看预测结果,且在预测系统交互界面中展示顶板巷布置位置示意。
3.2 预测系统的软件实现
采用PyQt5和python语言开发了顶板瓦斯抽采巷布置智能预测系统软件。该系统软件展示了训练模型所用的样本数据,并可将软件界面输入的相应参数保存到文件中(图7);对于覆岩结构参数量化,可通过在系统中输入岩层岩性和层厚,系统即可自动计算得到硬岩岩性比例系数(图8);另外,用户也可查看预测系统调用的GA–BP模型参数(图9)。该系统软件通过识别输入参数栏中的采厚、覆岩结构、倾向长度、倾角和埋深参数的数值,调用训练好的GA–BP神经网络模型,计算得到顶板巷的空间布置位置,并显示在预测结果栏中,软件界面如图10所示。
例如,寺家庄煤矿15110工作面煤层均厚为5.67 m,可采长度1 915 m,倾向长度180 m,硬岩岩性比例系数为0.35,平均倾角为5°,工作面埋深为327 m,将上述相关物理指标参数输入系统界面后,得到顶板巷的预测结果如图10所示。
同时,在15110工作面实际生产过程中,考虑到覆岩地质条件的影响,15110工作面顶板巷布置在9号煤层下部的砂质泥岩中,其距15煤层顶板平均垂距为53.5 m,与回风巷水平距离为47.8 m。顶板巷的抽采效果如图11所示[48]。
由此可知,15110工作面上隅角和回风巷区域内瓦斯体积分数在0.6%左右,未发生瓦斯超限现象,平均瓦斯抽采量约72 m3/min,抽采率为75.5%,且抽采量随推进距的增加呈现不断增加的趋势。而与之相邻的15106工作面顶板巷瓦斯抽采率为75.6%[35],两者的顶板巷瓦斯抽采效果基本持平,这进一步验证了该预测方法的可靠性。
4. 结 论
1) 在深入分析顶板巷布置原则及其布置位置的主控因素的基础上,提出了一种基于GA–BP神经网络模型的顶板巷布置位置智能预测方法,为顶板巷布置精准设计提供了一种新的技术手段。
2) 采用灰色关联分析法确定了采厚、埋深、覆岩结构、煤层倾角、倾向长度等5个物理指标是顶板巷布置位置的主控因素,并将其作为GA–BP神经网络模型的预测指标。
3) 设计开发了基于GA–BP神经网络模型的顶板瓦斯抽采巷布置位置智能预测系统,软件系统的预测结果与实际工况值的相对误差仅为0.43%~11.27%,在可接受的范围内。该研究可为顶板巷布置精准设计提供一定的参考。
-
表 1 样本数据与偏差
Table 1 Sample data and deviation
样点编号 本文算法 灰度重心法 曲线拟合法 测量值$y$/pixel 绝对偏差$\delta $/pixel 测量值$y$/pixel 绝对偏差$\delta $/pixel 测量值$y$/pixel 绝对偏差$\delta $/pixel 1 492.044 0.044 491.553 0.447 491.992 0.008 2 492.036 0.036 492.249 0.249 491.961 0.039 3 491.965 0.035 492.498 0.498 491.900 0.100 4 492.024 0.024 491.684 0.316 492.123 0.123 5 491.965 0.035 492.019 0.019 491.913 0.087 6 492.014 0.014 492.458 0.458 492.031 0.031 7 491.988 0.012 491.538 0.462 492.092 0.092 8 491.963 0.037 492.212 0.212 491.971 0.029 9 492.025 0.025 491.560 0.440 492.169 0.169 10 491.997 0.003 491.957 0.043 492.175 0.175 绝对偏差均值MAE/pixel 0.027 0.314 0.085 表 2 不同中心线提取算法的标准误差
Table 2 Standard errors of different centerline extraction algorithms
提取方法 本文算法 灰度重心法 曲线拟合法 $\sigma $/pixel 0.029 0.326 0.108 -
[1] 王国法,刘 峰,孟祥军,等. 煤矿智能化(初级阶段)研究与实践[J]. 煤炭科学技术,2019,47(8):1−36. WANG Guofa,LIU Feng,MENG Xiangjun,et al. Research and practice on intelligent coal mine construction(primary stage)[J]. Coal Science and Technology,2019,47(8):1−36.
[2] 葛世荣,胡而已,裴文良. 煤矿机器人体系及关键技术[J]. 煤炭学报,2020,45(1):455−463. GE Shirong,HU Eryi,PEI Wenliang. Classification system and key technology of coal mine robot[J]. Journal of China Coal Society,2020,45(1):455−463.
[3] 李东晓. 机器人技术在煤矿自动化中的应用[J]. 煤炭科学技术,2007(5):62−64. LI Dongxiao. Application of robot technology to mine automation[J]. Coal Science and Technology,2007(5):62−64.
[4] 代 伟,赵 杰,杨春雨,等. 基于双目视觉深度感知的带式输送机煤量检测方法[J]. 煤炭学报,2017,42(S2):547−555. DAI Wei,ZHAO Jie,YANG Chunyu,et al. Detection method of coal quantity in belt conveyor based on binocular vision depth perception[J]. Journal of China Coal Society,2017,42(S2):547−555.
[5] 黄世顶,韩 雷,芦 伟,等. 带式输送机巡检机器人研究[J]. 起重运输机械,2021(1):79−83. HUANG Shiding,HAN Lei,LU Wei,et al. Research on patrol robot of belt conveyors[J]. Hoisting and Conveying Machinery,2021(1):79−83.
[6] YULONG H,HAILI S,RUOFEI Z. Three-dimensional Linear Restoration of a Tunnel Based on Measured Track and Uncontrolled Mobile Laser Scanning[J]. Sensors (basel, Switzerland),2021,21(11):3815.
[7] 于希宁,杨跃权,刘红军,等. 露天煤场存煤量自动计量装置的研制及应用[J]. 仪器仪表学报,2000(2):194−195,199. YU Xining,YANG Yuequan,LIU Hongjun,et al. Research and application of the automatic measuring instrument for estimating coal storage volume in open stockyard[J]. Chinese Journal of Scientific Instrument,2000(2):194−195,199.
[8] 袁 洋. 激光盘煤系统的研究[D]. 保定: 华北电力大学, 2009. YUAN Yang. Research on the laser coal stocktaking system[D]. Baoding: North China Electric Power University, 2009.
[9] 李克伟. 基于机器视觉技术的带式输送机输送量检测系统的应用与研究[D]. 沈阳: 东北大学, 2014. LI Kewei. The research and application of the detection system to measure the material transportation on conveyor belt based on machine vision technology[D]. Shenyang: Northeastern University, 2014.
[10] 张延芝. 基于激光扫描技术的带式输送机非接触式输送量检测系统的应用与研究[D]. 沈阳: 东北大学, 2014. ZHANG Yanzhi. The application and research of the system to measure the conveying capacity by non-contact type on conveyor belt[D]. Shenyang: Northeastern University, 2014.
[11] 李 萍. 基于机器视觉的散状物料动态计量系统研究[D]. 北京: 中国矿业大学(北京), 2017. LI Ping. Research on dynamic measurement system of bulk material based on machine vision[D]. Beijing: China University of Mining & Technology(Beijing), 2017.
[12] 关丙火. 基于激光扫描的带式输送机瞬时煤量检测方法[J]. 工矿自动化,2018,44(4):20−24. GUAN Binghuo. Detection method of instantaneous coal quantity of belt conveyor based on laser scanning[J]. Industry and Mine Automation,2018,44(4):20−24.
[13] DENG E C, ZHAO X Q, YUAN Z P. Research on measuring system of bulk cargo[C]//2015 6th International Conference on Manufacturing Science and Engineering. Atlantis Press, 2015: 1882-1887.
[14] 何 雯. 精确提取光饱和与物体边缘线结构光条纹中心的方法研究[D]. 咸阳: 西北农林科技大学, 2018. HE Wen. Research on the extraction accuracy of laser stripe center in light saturation and at the edge of the object[D]. Xianyang: Northwest A&F University, 2018.
[15] LONG N Q,BUCZEK M M,HIEN L P,et al. Accuracy assessment of mine walls’surface models derived from terrestrial laser scanning[J]. International Journal of Coal Science & Technology,2018,5(3):328−338.
[16] 黎 明. 三维激光测量技术及应用研究[D]. 杭州: 浙江大学, 2005. LI Ming. Research on 3d laser measurement technology and its application[D]. Hangzhou: Zhejiang University, 2005.
[17] 胡而已. 基于激光扫描的综放工作面放煤量智能监测技术[J]. 煤炭科学技术,2022,50(2):244−251. HU Eryi. Research on intelligent monitoring technology of coal caving in fully mechanized caving face based on laser scanning[J]. Coal Science and Technology,2022,50(2):244−251.
[18] 徐兆军,倪申健,马辉贤. 基于斜射式激光三角测量锯材宽度的研究[J]. 光学技术,2018,44(3):339−343. XU Zhaojun,NI Shenjian,MA Huixian. Study on the measurement of sawn timber width based on oblique laser triangulation system[J]. Optical Technique,2018,44(3):339−343.
[19] 方雅媚,王红军,黄矿裕,等. 复杂背景下目标三维重构的图像预处理方法[J]. 激光与光电子学进展,2019,56(13):118−126. FANG Yamei,WANG Hongjun,HUANG Kuangyu,et al. Image-preprocessing algorithm for three-dimensional reconstruction of targets in complex backgrounds[J]. Laser & Optoelectronics Progress,2019,56(13):118−126.
[20] JIANG T, YU Z, CAO G, et al. Effect of laser intensity changing on accuracy of size measurement by laser scanning[C]//Proceedings of the Society of Photo-optical Instrumentation Engineers (spie), 2005: 195-198.
[21] 张凤生. 高斯光束的数值模拟新方法[J]. 光子学报,2008(6):1259−1262. ZHANG Fengsheng. A new numerical simulation method for gaussian beam[J]. Acta Photonica Sinica,2008(6):1259−1262.
[22] 朱梦楠. 三维扫描中激光光条中心提取精度的研究[D]. 咸阳: 西北农林科技大学, 2017. ZHU Mengnan. Research on extraction accuracy of laser stripe center in three-dimensional scanning[D]. Xianyang: Northwest A&F University, 2017.
[23] 胡而已. 运动物体三维形貌测量方法与实验系统研究[D]. 武汉: 华中科技大学, 2009. HU Eryi. Methods study and experimental device design for measuring the 3-d profile of a moving object[D]. Wuhan: Huazhong University of Science and Technology, 2009.
[24] 潘 硕,邵新杰,李晓磊,等. 基于块匹配的激光条纹亚像素中心提取[J]. 激光与红外,2021,51(6):814−819. PAN Shuo,SHAO Xinjie,LI Xiaolei,et al. Sub-pixel extraction of laser stripe center based on block-matching[J]. Laser & Infrared,2021,51(6):814−819.
-
期刊类型引用(3)
1. 张百颂. 高突矿井下行揭煤钻孔施工工艺及瓦斯抽采效果研究与应用. 煤炭与化工. 2025(05): 118-123 . 百度学术
2. 兰安畅,袁欣鹏,付智博. 矿井瓦斯智能化抽采发展现状. 科技和产业. 2024(18): 240-245 . 百度学术
3. 姬静伟. 煤层瓦斯抽采有效半径分析. 山东煤炭科技. 2024(12): 60-65 . 百度学术
其他类型引用(1)