一种基于图像识别的井巷通风阻力系数计算方法技术

本发明专利技术公开了一种基于图像识别的井巷通风阻力系数计算方法，属于矿山智能通风领域。通过对待测风阻巷道的围岩拍照获得多张包含粗糙度信息照片，图片采用高斯差分算法识别出图片中粗糙度，并提取照片中粗糙度的键点。结合相机内外参数、坐标系变换将二维信息转换在世界坐标系并构件点云三维模型。将点云模型分割为几个区域，通过各区域内的散点拟合三维平面，计算各个点云到平面距离的均值作为该区域的粗糙度。同样的方法计算其他分割区域的粗糙度。对所有分割区域的粗糙度求平均值作为这张照片的粗糙度；同样方法识别其他照片中的围岩粗糙度，将所有照片识别的粗糙度计算均值作为待测巷道的粗糙度。最后联立Colebrook数学模型和摩

【技术实现步骤摘要】
一种基于图像识别的井巷通风阻力系数计算方法


[0001]本专利技术专利涉及矿井通风防灾
的通风阻力系数智能预测技术，尤其涉及图像 识别的井巷通风阻力系数计算方法。

技术介绍

[0002]矿井通风作为矿山的“呼吸系统”对采掘生产、防治瓦斯积聚和井巷内遗煤自燃、营 造良好的作业环境具有重要意义。空气沿井巷中流动时，由于风流的黏滞性、惯性和井巷壁 面等对井巷的阻滞、扰动作用而形成通风阻力造成巷道的需风地点供风不足。通风摩擦阻力 系数作为衡量通风难以程度的重要参数之一，如何合理的预测长期以来是研究的热点问题。
[0003]目前，预测巷道内通风摩擦阻力系数的方法可归结于三类：(1)经验取值法，即根据以往的经验数据，在某一范围内随机取值；(2)现场测试法，通过现场的实测获得数据；(3)数值计算法，通过建立数值计算模型，通过计算机计算而得。
[0004]经验法易受到主观因素的影响，预测结果和实际结果误差较大，难以客观反映真实结 果；现场测试易受到井下复杂环境的影响，其次开采年限较长的矿井巷道分支有数千条之多， 无法对全部巷道进行实际测量；数值计算法需要求解偏微分方程，采用不同数学模型结果差 别较大，并且迭代次数对计算机的性能要求较高。除上述原因外，风流本身脉动影响导致测 量结果“测不准、总在变”，难以获得令人满意的测量结果。针对通风摩擦阻力系数，部分学 者们的研究结果具有代表性。
[0005]张庆华等提出一种基于大数据的通风巷道摩擦阻力系数快速确定方法，通过建立巷道 通风参数数据与巷道摩擦阻力系数的分类数据母几何和子几何，新城各行综合、动态更新、 精确的摩擦阻力系数系统大数据集，通过巷道顺序，在其中搜索类似巷道的数据集，然后在 相似巷道的数据集中搜索相似的影响因子，快速确定摩擦阻力系数。
[0006]杨卓明提出一种通风摩擦阻力系数测量方法和通风阻力系数测量装置。在巷道中设置 可调节的阻断程度的阻断物；分N此调节所述阻断物对所述巷道的阻断程度，测量N个风量 数值及所述N个风量数值对应的N个巷道压力差值；根据所述N个风量数值及其对应的N 个巷道压力差值确定所述巷道的通风阻力系数。
[0007]高菊茹等提出一种测试曲线隧道通风沿程摩擦阻力系数的实验装置。根据相似比选用 的与原型巷道壁面粗糙度相应的贴合在曲线模型隧道内表面的砂纸、静压测试管、带测试小 孔的整流圆管、微压计、通过变频器控制的实验风机所构成，曲线模型隧道的一段与直线模 型隧道连接、另一端与整流圆管连接，整流圆管与实验风机连接，静压测试管安装在曲线模 型隧道的待测长度的两端点位置，整流圆管用L型皮托管与微压计相连。由风机提供测试风 量，将测得的静压差值和动压差值带入公式计算即可得到曲线模型隧道内的摩擦阻力系数。
[0008]刘承东等提出一种综合管廊通风系统当量摩擦阻力系数确定方法。分别通过理论
公式 和数值计算获得忽略管线和支架铺设情况时综合管廊舱室内的摩擦阻力系数，对比二者获得 第一修正系数；分别通过理论公式和数值计算获得铺设管线和支架情况时综合管廊舱室内的 摩擦阻力系数，对比二者获得第二修正系数；根据第一修正系数和第二修正系数获得管廊通 风系统的当量摩擦阻力系数。
[0009]高科等提出一种矿井通风摩擦阻力系数测定的方法。将壁面粗糙度作为通风摩擦阻力 系数的固有属性，基于激光三维扫描仪获得围岩表面点云，遍历每段小巷道点云得到待测巷 道的散点集合并提取所有点云的坐标点，根据散点坐标求得巷道壁面粗糙度；通过点云计算 小段巷道的体积和实验比例系数；对所有的小巷道求出巷道实验比例系数平均值，进而求出 整个巷道的通风摩擦阻力系数。这种算法对于粗糙度的处理较为简化，测量易受尺度效应的 影响。

技术实现思路

[0010]本专利技术要解决的技术问题是提高巷道内通风摩擦阻力系数测试的智能化程度，减少井 下测阻投入的人力和财力，为此提出一种基于图像识别的通风摩擦阻力系数测定新方法，本 方法不受环境因素的干扰，直接从图片中重构巷道围岩的三维点云预测粗糙度，进一步通过 粗糙度预测通风摩擦阻力系数。该测阻方法为矿井通风系统的优化设计提供设计参数。
[0011]为解决上述技术问题，本专利技术所采取的技术方案是：提供的一种基于图像识别的井巷通 风阻力系数计算方法，流程如图1所示，包括如下步骤。
[0012]步骤1：获取围岩表面图片，在巷道内架设相机，并对巷道的围岩拍照，在此过程中不 断调整拍摄角度，照片需要覆盖预测通风摩擦阻力系数巷道区段围岩的所有区域，得到的照 片集合P，遍历该区域所有的照片集合P＝{p1,p2,
…
,p
m
，
…
，p
n
}。
[0013]步骤2：识别照片中的围岩粗糙度，如图2所示采用高斯函数差分法处理照片，本步骤 以集合P中的任意一张照片为例说明粗糙度识别的计算过程，集合内的其他照片都采用相同 的计算过程处理，如果照片的长度为x m，宽度为y m，滤波参数分别为σ1和σ2，那么两次 滤波后的图片表示如下：将两幅图像g1和g2在不同参数下的高斯滤波结果相减，g3是得到高斯差分图，表示如下。
[0014]步骤3：提取照片中识别出的边缘点作为三维重构的关键点，建立世界坐标系(world)、 相机坐标系(camera)、图像坐标系(picture)、像素坐标系(pile)间的映射关系，涉及外参矩阵和 内参矩阵，相机的外参是矩阵的旋转R和平移T，相机的内参是焦距的参数f
x
和f
y
，主点的 偏移量x0和y0；
①
首先需要根据外参矩阵将世界坐标系O
world
转换为以相机
为中心的坐标系 O
camera
，如图3所示，转换涉及坐标系的旋转和平移，此时相机坐标系和世界坐标系方向一致；
ꢀ②
根据相机成像原理(如图4所示)，借助相机的内参将像素点映射在相机坐标系中，为了统 一标准和以后方便计算，归一化处理相机坐标系中的位置坐标，使其变为与成像平面平行且 距离相机原点距离相等的平面，此时也是建立三维空间点和二维空间点的关系，然而，物理 平面坐标系中的图像是以图像的左上角为原点，而我们之前的归一化相平面坐标系中是成像 中心O
imagin
(u0,v0)点为原点(如图5所示)，因此需要平移变换。
[0015]对上述过程做个总结，即从像素坐标系
→
图像坐标系
→
相机坐标系
→
世界坐标系；计 算过程如公式(3)，等式最右端矩阵包含世界坐标系中的点的信息，将点的集合命名为S， 遍历集合内的所有点S＝(s1，s2，...，s
m
，...，s
n
)，其中任意元素为)，其中任意元素为
[0016]步骤4：通过方差减小规则(SDR)将点的集合S分割为多个子集合(如图6所示)，遍历 所有的子集合S＝{S1,S2,...,S
m
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于图像识别粗糙度的通风摩擦阻力系数预测方法，其特征在于包括如下步骤：步骤1：获取围岩表面图片；在巷道内架设相机，并对巷道的围岩拍照，在此过程中不断调整拍摄角度，照片需要覆盖预测通风摩擦阻力系数巷道区段围岩的所有区域，得到的照片集合P，遍历该区域所有的照片集合P＝{p1,p2,
…
,p
m
，
…
，p
n
}；步骤2：识别照片中的围岩粗糙度，如图2所示采用高斯函数差分法处理照片，本步骤以集合P中的任意一张照片为例说明粗糙度识别的计算过程，集合内的其他照片都采用相同的计算过程处理，如果照片的长度为x m，宽度为y m，滤波参数分别为σ1和σ2，那么两次滤波后的图片表示如下：将两幅图像g1和g2在不同参数下的高斯滤波结果相减，g3是得到高斯差分图表示如下：步骤3：提取照片中识别出的边缘点作为三维重构的关键点，建立世界坐标系(world)、相机坐标系(camera)、图像坐标系(picture)、像素坐标系(pile)间的映射关系，涉及外参矩阵和内参矩阵，相机的外参是矩阵的旋转R和平移T，相机的内参是焦距的参数f
x
和f
y
，主点的偏移量x0和y0；
①
首先需要根据外参矩阵将世界坐标系O
world
转换为以相机为中心的坐标系O
camera
如图3所示，转换涉及坐标系的旋转和平移，此时相机坐标系和世界坐标系方向一致；
②
根据相机成像原理(如图4所示)，借助相机的内参将像素点映射在相机坐标系中，为了统一标准和以后方便计算，归一化处理相机坐标系中的位置坐标，使其变为与成像平面平行且距离相机原点距离相等的平面，此时也是建立三维空间点和二维空间点的关系，然而物理平面坐标系中的图像是以图像的左上角为原点，而我们之前的归一化相平面坐标系中是成像中心O
imagin
(u0,v0)点为原点(如图5所示)，因此需要平移变换；步骤4：通过方差减小规则(SDR)将点的集合S分割为多个子集合(如图6所示)，遍历所有的子集合S＝{S1,S2,...,S
m
,...,S
n
}，分割准则如公式(4)，其中，sd是集合的方差，|E|是E集合里数目；步骤5：通过偏最小二乘法分别拟合S子集合中的点，平面的集合A，遍历集合A＝{A0,A1,A2,
…
,A
...

【专利技术属性】
技术研发人员：高科，戚志鹏，刘玉姣，
申请(专利权)人：辽宁工程技术大学，
类型：发明
国别省市：