Research on intelligent control of air volume of mine ventilation network
-
摘要: 现有矿井通风网络风量智能优化算法在求解调风参数时普遍存在模型复杂、收敛速度慢、易陷入局部最优等缺陷，同时也缺乏与调风分支优化选择相结合的研究。针对上述问题，提出了一种基于改进天牛须搜索（BAS）算法的矿井通风网络风量智能调控方法。首先，以用风分支的风量需求为优化目标，构建风量优化调节数学模型，针对该模型中的风量调节约束条件，采用不可微精确罚函数并结合模拟退火算法优化惩罚项，实现模型的去约束化。然后，通过求解灵敏度矩阵，结合风量灵敏度和分支支配度理论选择最优的调节分支集，确定其风阻调节范围，并作为模型的初始解集。最后，基于改进BAS算法求解出最优调风参数，进而控制对应的调风设施，实现风量调控。基于矿井通风实验平台对该方法的可靠性进行实验验证，结果表明：相比于标准BAS算法和粒子群优化（PSO）算法，改进BAS算法综合寻优性能更优越，解得的风量平均值和最优解均高于PSO算法和标准BAS算法，平均运行时间虽略长于标准BAS算法，但远短于PSO算法，平均收敛代数最多，精度最高，容易跳出局部循环得到最优解；在设定风量调节目标后，基于改进BAS算法的矿井通风网络风量智能调控方法可快速精准求解出待调分支的风量最优值，调节后的分支风量满足矿井安全生产的调风要求，风量上调高达 46.5%。Abstract: The existing intelligent optimization algorithm of air volume of mine ventilation network has the defects of complex model, slow convergence speed, easy falling into local optimum when solving the air adjustment parameters. There is a lack of research on the combination of optimal selection of air adjustment branches. To solve the above problems, an intelligent control method of air volume of mine ventilation network based on improved beetle antennae search (BAS) algorithm is proposed. Firstly, the mathematical model of air volume optimal adjustment is established by taking the air volume demand of the air consumption branch as the optimization objective. In view of the air volume adjustment constraint conditions in the model, the non-differentiable exact penalty function and the simulated annealing algorithm are adopted to optimize the penalty term, so that the model is unconstrained. Secondly, by solving the sensitivity matrix and combining the theory of air volume sensitivity and branch dominance, the optimal adjustment branch set is selected. The air resistance adjustment range is determined as the initial solution set of the model. Finally, based on the improved BAS algorithm, the optimal air adjustment parameters are solved. The corresponding air adjustment facilities are controlled to realize air volume adjustment. The reliability of the method is verified by experiments based on the mine ventilation experimental platform. The results show that compared with the standard BAS algorithm and particle swarm optimization (PSO) algorithm, the improved BAS algorithm has superior comprehensive optimization performance. The average value and optimal solution of air volume are higher than those of the PSO algorithm and standard BAS algorithm. Although the average running time is slightly longer than the standard BAS algorithm, it is far shorter than the PSO algorithm. The average convergence algebra is the most, the precision is the highest, and it is easy to jump out of the local loop to get the optimal solution. After setting the air volume adjustment target, the intelligent control method of air volume of the mine ventilation network based on the improved BAS algorithm can quickly and accurately solve the optimal value of the air volume of the branch to be adjusted. The adjusted branch air volume meets the air volume adjustment requirements of mine safety production, and the air volume is increased up to 46.5%.
-
0. 引言
煤矿巷道的形变和破坏是煤矿生产实践中需要着重考虑解决的问题。煤矿巷道形变是指巷道岩层在受到外力因素作用下产生的形状和尺寸的变化，巷道形变会造成井下交通中断、生产设备损坏、人员伤亡等事故，严重威胁煤矿正常生产。因此，在煤矿开采过程中，及时、科学、准确地监测巷道形变，掌握巷道状态及形变趋势，对保障煤矿安全生产具有十分重要的意义[1]。
随着计算机技术的快速发展，多种数值模拟技术被广泛应用于巷道形变的研究中，如有限元法、离散元法、边界元法、有限差分法、十字观测法等。王伟等[2]运用理论计算与数值模拟分析，分析了动压巷道的应力分布规律，为巷道的布设位置选取与支护提供了指导。姜耀东等[3]对开滦矿区赵各庄矿、唐山矿深部开采过程中巷道形变、破坏特征进行了研究，采用FLAC2D数值模拟软件研究地应力与巷道岩层之间的密切关系。余伟健等[4]利用离散裂隙网络在FLAC3D软件中对矿山现场调查得到的裂隙进行了重构还原，研究分析了裂隙发育岩体巷道及典型支护结构体的变形破坏特征。杨景贺[5]采用离散元数值模拟方法，建立了不同支护条件下的高应力软岩巷道模型，研究了高应力软岩巷道变形破坏和支护加固的机理。
目前，针对巷道形变破坏的研究以现场实测为主，通过实测可掌握巷道形变破坏情况，直观了解导致巷道形变破坏的因素及作用效果。叶美图等[6]针对现有矿井巷道表面位移十字交叉测量方法存在的测量不便、耗时长、工作强度大、误差不可控等问题，设计了一种矿井巷道表面位移激光测量装置。该装置基于激光测距原理，利用激光测距仪完成十字定心和变形测量，通过单个激光测距仪可分别测得巷道中每个测试站顶板、底板、左帮和右帮4个方向的表面位移量。井下试验结果表明，该装置实现了无绳测量，测量效率高，结果准确。宋康磊等[7]采用地质雷达测定巷道松动圈，研究了极薄煤层跨采条件下的松动圈大小及巷道形变规律。荣耀等[8]提出了综采工作面三维激光扫描建模总体方案，对煤壁与顶板交线提取、标靶球检测、点云拼接及坐标转换等关键技术进行了研究，实现了三维地质坐标系下煤壁与顶板交线信息的近实时获取，该信息可直接发送给采煤机滚筒，为采煤机下一刀截割提供数据参考。杜江丽等[9]提出了多边形分割法，将模型的重点关注区域划分为多个多边形，对矿区边坡形变监测数据进行研究。亓玉浩等[10]提出了一种综采工作面实时三维建图方法，可满足综采工作面监控的精度需求。刘晓阳等[11]针对岩性变化和构造的分析监测越来越困难的问题，研究了当前巷道顶板稳定性监测方法。王峰[12]提出了以煤层透明化、综采装备透明化、决策及控制透明化为核心的基于透明工作面的智能化开采概念。邹筱瑜等[13]面向移动机器人在煤矿巷道的建图需求，提出了一种基于集成式因子图优化的煤矿巷道移动机器人三维地图构建方法。高晓进等[14]提出了软岩巷道超前爆破预裂顶板防控技术，以解决软岩巷道结构的形变问题。
上述方法为巷道形变监测提供了可行的解决方案，但还存在如下问题：① 单一的形变量指标难以准确反映巷道形变演化规律及围岩状态与发展趋势。② 巷道形变具有区域性、连续性、突发性的特点，现有方法存在监测站点单一、测点有限的问题，难以实现大范围、多站点、区域性同步测量。因此需要选取多个指标建立巷道形变监测指标体系，且需要研究一种具有区域性、连续性、实时性特点的巷道形变监测方法，及时、准确、全面地监测巷道形变以确保煤矿开采安全。
三维激光扫描技术可快速获取被测物体表面三维坐标，点云数据精度高、密度大，包含反射强度等信息[15]，可实现对巷道形变连续监测。因此，针对煤矿巷道形变监测面临的问题，本文提出了一种基于三维激光扫描的煤矿巷道形变监测方法。首先，利用深度学习模型VoxelNet对三维激光扫描数据进行检测去噪。然后，利用alphashape算法对断面离散点进行拟合，并采用基于差值的方法，对数据进行多维差值计算，得到巷道形变的具体数据，实现矿区巷道形变监测的全覆盖。最后，利用中煤大同能源有限责任公司塔山煤矿为期1个月的实际监测数据，对巷道形变进行截面分析和三维整体分析。
1. 煤矿巷道形变监测方法
1.1 三维激光扫描技术
三维激光扫描技术是一种先进的全自动、高精度立体扫描技术，具有操作简单、速度快、灵活性高、精度高等特点[16]，数字化、自动化、实时性强，解决了目前空间信息技术发展的瓶颈。
三维激光扫描技术借助水平方位偏转控制器及对应的高度角偏转控制器科学控制反射棱镜的实际转动，使激光测距系统发射的激光能够顺着不同的坐标轴实时移动并扫描测量。扫描仪自身定义的坐标系统中， X轴在横向扫描面内，Y轴在横向扫描面内与X轴垂直，Z轴与横向扫描面垂直。借助测得的三维扫描仪中心到目标点的斜距S、激光束水平方向偏转角a和竖直方向偏转角b来校准并计算激光脚点的实际三维坐标P（X，Y，Z）[17]。
$$ X = S \cos b \cos a $$ (1) $$ Y = S \cos b \sin a $$ (2) $$ Z = S \sin b $$ (3) 三维激光点坐标计算原理如图1所示，通过测量大量物体表面点位的三维坐标，读取并集成现场物体的反射率和纹理等信息，生成被测物点云模型，获取高精度、高分辨率的三维巷道点云数据。
1.2 点云去噪
对三维巷道点云数据去噪是煤矿巷道形变监测得以有效应用的重要预处理步骤[18]。将噪点设为巷道内的人物和内壁附着的螺丝钉和管道等，通过VoxelNet网络检测出噪点后，再对数据进行分离。
首先使用点云标注工具Semantic Segmentation Editor制造数据集，标注人物和管道等信息后送入VoxelNet训练。然后使用VoxelNet中的Voxel将三维空间中的一定空间划分为一个格子，并使用PointNet网络对格子的数据进行特征提取。最后用提取出来的特征代表这个格子，并放回到三维空间。将无序的点云数据变成高维特征数据，并使用三维卷积抽取有序的三维Voxel数据。
1.3 点云旋转
给定2个来自不同坐标系的三维数据点集，找到2个三维数据点集的空间变换关系，使2个三维数据点集处于同一坐标系中，即可实现三维点云配准[19]。配准目标是在全局坐标框架中找到单独获取视图的相对位置和方向，使2个三维数据点集之间的相交区域完全重叠。通过求点云M和点云N 之间的旋转矩阵和平移向量，使点云M和点云N的坐标处于同一坐标系下。通过欧拉角生成旋转矩阵，将点云数据进行平移和旋转。设新的空间直角坐标轴分别为X'，Y'，Z'，绕X'，Y'，Z' 3个轴旋转的角度分别为α，β，γ，绕X'，Y'，Z'旋转分别得到的旋转矩阵为
$$ {{\boldsymbol{R}}_{X'}}\left( \alpha \right) = \left[ \begin{gathered} 1\quad \quad 0\quad \quad 0 \\ 0\quad \cos \alpha \quad - \sin \alpha \\ 0\quad \sin \alpha \quad \cos \alpha \\ \end{gathered} \right] $$ (4) $$ {{\boldsymbol{R}}_{Y'}}\left( \beta \right) = \left[ \begin{gathered} \cos \beta \quad \;0\quad \sin \beta \\ \quad 0\quad \quad 1\quad \;\;0 \\ - \sin \beta \;\;\;0\quad \cos \beta \\ \end{gathered} \right] $$ (5) $$ {{\boldsymbol{R}}_{Z'}}\left( \gamma \right) = \left[ \begin{gathered} \cos \gamma \quad - \sin \gamma \;\;\;0 \\ \sin \gamma \quad \;\cos \gamma \quad \;0 \\ 0\;\;\;\;\;\;\;\;\;\;\;\;0\quad \,\,\cos \gamma \\ \end{gathered} \right] $$ (6) 对3个旋转矩阵的点云坐标进行逐元素相乘，得到最终的旋转矩阵R。
$$ {\boldsymbol{R}} = {{\boldsymbol{R}}_{Z'}}\left( \gamma \right) {{\boldsymbol{R}}_{Y'}}\left( \beta \right) {{\boldsymbol{R}}_{X'}}\left( \alpha \right) $$ (7) 1.4 点云连续断面
将巷道切分为连续的断面，根据断面对巷道进行形变分析，分别向点云的X'OY'和Y'OZ'面进行数据投影，计算当前断面的中心点$ {P^*} \left( {{x^ * },{y^ * },{{\textit{z}}^ * }} \right) $，其中，$ {x^ * },{y^ * },{{\textit{z}}^ * } $为中心点$ {P^*} $在X'，Y'，Z' 3个坐标轴上对应的数值。
随后，构造与Z'轴平行的平面，平面与Z'轴平行且与Y'轴相切，则该平面的法向量u=[0 1 0]，该平面的方程为$ y - {y^ * } = 0 $。由于落于平面的点比较稀疏，所以赋予平面一定的厚度，将点到平面一定距离d内的点都算作该平面的点。
1.5 形变分析
在实际扫描作业中，由于障碍物遮挡或受巷道空间扫描角度的限制，会出现部分点云空洞，难以保证断面上处处都有激光斑点，所以需要对投影后的断面离散点集进行拟合。alphashape算法可控制点集边界拟合以创建非凸区域，采用该算法拟合投影后的断面离散点集，其计算圆心的原理如图2所示，计算步骤如下：
![]()
图 2 alphashape算法计算圆心原理Figure 2. Diagram of alphashape algorithm calculating the center of a circle1） 选取任意一点$ p \left( {x,y} \right) $，以半径为r的圆进行滚动，在断面点集中搜索所有距离$ p $点2r以内的点的集合，形成点集Q。
2） 选取点集Q中的任意点$ {p_1} \left( {{x_1},{y_1}} \right) $，根据$ p $和$ {p_1} $2个点的坐标和半径r，计算出经过$ p $和$ {p_1} $2个点且半径为r的圆心${p_2} \left( {{x_2},{y_2}} \right) $，${p_3} \left( {{x_3},{y_3}} \right) $。
$$ \left\{ \begin{array}{l} {x_2} = x + 1/2({x_1} - x) - H ({y_1} - y) \\ {y_2} = y + 1/2({y_1} - y) - H (x - {x_1}) \\ {x_3} = x + 1/2({x_1} - x) + H ({y_1} - y) \\ {y_3} = y + 1/2({y_1} - y) + H (x - {x_1}) \\ \end{array} \right. $$ (8) 式中H为参数。
3） 计算点集Q中（除$ {p_1} $点）剩余点分别到$ {p_2},{p_3} $点的距离。若所有点到$ {p_2} $或$ {p_3} $的距离均大于半径r，则证明$ p $点为边界点。
4） 若剩余的点到$ {p_2} $或$ {p_3} $点的距离不全都大于半径r，则遍历点集Q内的所有点，将其轮换作为$ {p_1} $点，重复上述过程。若存在某点满足2）和3），则该点为边界点，终止该点的判断，继续判断下一点。若点集Q中所有临近点中均不存在$ {p_1} $这样的点，则$ p $点为非边界点。
根据alphashape算法拟合巷道断面轮廓图，构建不同监测时期同一位置的巷道断面，将其叠加后，采用对比不同时期同一位置巷道断面的高度和宽度等信息的方法，可以更直观地分析巷道形变趋势。具体分析过程：首先，由alphashapes算法计算第一期巷道断面的轮廓点，分别计算相邻2个点$p'\left( {{x_0},{y_0},{{\textit{z}}_0}} \right), q\left( {{x_1},{y_1},{{\textit{z}}_1}} \right)$所构成的线段$ {l_1} $的方向向量，默认$p',q$处于同一平面内，故只保留X轴和Z轴值，记为$p'\left( {{x_0},{{\textit{z}}_0}} \right), q\left( {{x_1},{{\textit{z}}_1}} \right)$。然后，计算过$p '$点与$ {l_1} $垂直的直线$ {l_2} $，$ {l_2} $即法向直线。最后，计算第二期与第一期相同法向直线的点的坐标，计算2个点间的距离，将其记为形变值。
2. 工程实践
2.1 工程概况
中煤大同能源有限责任公司塔山煤矿的开采工作面为30507工作面，工作面采用单巷布置方式，即在工作面分别布置胶带巷和回风巷。工作面西侧布置回风巷，东侧布置胶带巷，胶带巷东侧为实煤区，回风巷西侧为30509工作面采空区（2018年6月停采），30507工作面与30509工作面之间区段煤柱为8 m，因此将30507工作面回风巷作为监测地点，进行为期1个月的实际形变监测。
通过对30507工作面进行实地勘查，其具有以下特点：
1） 靠近工作面的回风巷区域形变量大，已有控制点坐标经过一定时间会发生变化，无法直接利用。
2） 对于井下空间地理坐标系，除用全站仪进行测点定位外，没有更好的方式将点云数据和矿井坐标系进行联系。
3） 矿区井下环境相对复杂，煤尘、空气湿度、振动等工况极易影响工作效率和数据精度。
2.2 实践方案
塔山煤矿正在回采的工作面为30507工作面，开始监测时采煤机所处位置在里程1 000 m处，监测范围为里程600～900 m。在测绘区域起始位置即里程600 m处，前后上下交错架设多个360°棱镜。根据30507回风巷实际情况，考虑到工作面附近棱镜位置会产生变化，因此以风门处于终采线以内的基准点为基准，测量各个360°棱镜的三维坐标，使棱镜与井下基准点处于同一坐标系下。采用后方交会的方式，在每次进行巷道形变监测数据采集时，以相同基准点作为初始测点的定向点，将矿井坐标系赋予点云数据，以保证每次进行形变监测工作所采集的点云数据都处于同一个空间坐标系下。
受巷道结构影响，只能选取接近于直线的布站路线进行三维激光扫描。为避免由于线性布站及井下煤尘、湿气等因素造成的拼接误差，相邻测站之间采用棱镜定位的方式进行拼接定位，选择应用可识别棱镜并记录棱镜信息的三维激光扫描仪。相邻测站间需保证有2个及以上共同的棱镜作为拼接基准。
将周期性采集的点云数据置于同一坐标系下，对多次的点云数据进行对比，得出巷道顶底板和煤柱侧、回采侧两侧壁的形变量。
2.3 数据采集
1） 控制测量。由于三维激光扫描仪具有独立的坐标系统，为了将点云数据坐标系转换为实际地理坐标系，需结合全站仪进行控制测量。在里程600 m处做3个能够固定在锚杆上的棱镜放置装置，并在其上安置3个360°棱镜。以风门处7号、9号、11号基准点为基准，采用全站仪测出3个棱镜在矿区地测坐标系中的坐标。
2） 在里程600～900 m段，按30 m的间距，在巷道两帮粉刷并做标记，作为扫描的架站位置；在里程600 m处开始架设三维激光扫描仪，并进行对中整平等操作。
3） 定位棱镜测量。扫描仪架设完毕后，使用棱镜扫描功能采集设置的360°棱镜点位数据，将定位棱镜信息同步至点云模型中。
4） 根据项目要求选择合适的扫描参数，建立测站文件夹，使用棱镜的扫描功能，识别测站间用于点云拼接所架设的棱镜，识别完成后开始三维激光扫描，获取该测站全景点云数据。
5） 检查数据质量，若数据质量合格，即可将扫描仪迁移至下一预设测站，并以同样步骤进行该测站的三维激光扫描。依次循环，直至待测区域全部扫描完成，即数据采集工作完成。
2.4 数据处理
数据处理过程只需将相邻测站间相同的棱镜目标进行匹配，即可达到点云拼接目的，过程如下：
1） 控制点坐标录入。在数据处理软件中，使用起始测站中采集的定位棱镜点创建目标，将全站仪测量的定位棱镜坐标赋予对应目标，控制点坐标录入，如图3所示。
2） 后方交会拼接。起始测站点云利用其中的定位棱镜点与已赋予现场坐标系坐标的目标点位匹配。第二测站点云以起始测站为基准，利用与起始测站共视的2个或2个以上棱镜进行后方交会，将其与起始测站拼接[20]。后续测站依次进行，得到定位于现场坐标系中的整体测区点云模型。
3） 数据去噪。为采集完整的巷道顶板、底板及煤柱侧、回采侧两侧壁数据，采用360°全景扫描仪采集现场所有可视数据。由于噪点数据[21]（图4）会对实施效率及结果精度造成影响，所以需对噪点数据进行去噪，去噪后数据如图5所示。
4） 数据输出。对数据进行拼接、去噪后，在确保点云定位信息固定的前提下，导出相应格式的点云模型，如图6和图7所示。
![]()
图 6 煤矿巷道点云数据内部视角Figure 6. Internal perspective of point cloud data in coal mine roadway![]()
图 7 煤矿巷道点云数据外部视角Figure 7. External perspective of point cloud data of coal mine roadway2.5 数据分析
在数据处理过程中，将所有三维点云模型置于同一坐标系下，因此在将同一测区的2期三维点云数据导入对比分析软件后，2期三维点云模型会自动放置在同一位置。
1） 三维整体分析。将置于同一位置的2期点云数据沿其法线方向进行差值分析，得到整体形变量，以色谱的形式将差值进行显示，不同颜色的分布表示不同区域2期点云模型偏差的分布，颜色的深浅表示该区域偏差大小。通过对模型上颜色的分布和深浅分析，得到该测段巷道形变的整体趋势，如图8所示。如果巷道受到向内的挤压，会变成红色，红色越深，表明挤压程度越大；如果巷道趋近于无变化，则会变成绿色；如果巷道向外扩张，会变成蓝色，蓝色越深，表明向外扩张程度越大。
2） 截面分析。以三维整体分析为基础，在整体巷道点云模型的竖直方向上按一定距离布置截面，截取2期点云在截面上的偏差，以色谱的形式将差值进行显示，截面分析结果如图9所示。可看出离工作面距离越近，附色模型的颜色越偏向红色和蓝色，且颜色越深，代表巷道形变越大。
两期巷道截面对比如图10所示，通过2个巷道轮廓的相对偏差，可看出该区域主要存在的形变：如果上轮廓向内偏离，则顶板发生塌陷；如果下轮廓向外偏离，则底板发生鼓起。
3） 预警记录。点选预警记录可查看预警位置、里程位置和最大形变值等信息，数据精确，可随时观察巷道形变情况，以便采取相应维护措施，巷道形变预警记录如图11所示。
3. 结论
1） 将三维激光扫描技术应用于塔山煤矿30507工作面的形变监测，将人工智能与传统方法相结合，提出了基于三维激光扫描的煤矿巷道形变监测方法。利用深度学习VoxelNet模型对点云数据进行去噪，利用alphashape算法对断面离散点进行拟合，并采用基于差值的方法，将数据进行多维差值计算以得到巷道形变数据。
2） 利用塔山煤矿的实际监测数据，对巷道三维点云数据进行了截面分析和三维整体分析。截面分析主要基于巷道三维点云按固定间距选取截面，比对2期点云截面数据，通过截面轮廓差异判断巷道形变情况；三维整体分析利用2期巷道三维点云数据进行整体比对，通过色谱显示比对结果，分析巷道形变区域和变化趋势。随着测点距工作面的距离越近，附色模型越偏向红色和蓝色，且颜色越深，表明巷道形变越大。
-
![]()
图 5 灵敏度
$ {d}_{\mathrm{4,18}} $ 随风阻$ {R}_{18} $ 的变化Figure 5. Variation of sensitivity
$ {d}_{\mathrm{4,18}} $ with air resistance$ {R}_{18} $ ![]()
图 6 不同算法所得分支4风量适应度曲线
Figure 6. Air volume fitness curves of branch 4 obtained by different algorithms
![]()
图 7 分支4瓦斯体积分数随风量变化曲线
Figure 7. Change curves of gas volume fraction in branch 4 with air volume
表 1 通风网络初始参数
Table 1 Initial parameters of ventilation network
分支
编号始节点 末节点 风阻/
(N·s2·m−8)初始风量/
(m3 ·s−1)最小需风量/
(m3 ·s−1)分支
编号始节点 末节点 风阻/
(N·s2·m−8)初始风量/
(m3 ·s−1)最小需风量/
(m3 ·s−1)1 ① ② 0.455 43.37 41.92 12 ⑥ ⑩ 1.376 5.48 2.59 2 ② ③ 0.208 20.60 16.85 13 ⑦ ⑨ 1.206 5.43 5.46 3 ② ④ 0.124 22.78 19.43 14 ⑧ ⑩ 0.336 5.62 4.26 4 ③ ⑪ 1.156 9.57 6.65 15 ⑨ ⑪ 0.209 14.79 9.94 5 ③ ⑤ 0.197 11.02 7.29 16 ⑩ ⑬ 0.137 11.11 6.71 6 ④ ⑥ 0.040 14.88 11.26 17 ⑪ ⑫ 0.074 24.37 20.38 7 ④ ⑬ 1.076 7.90 6.14 18 ⑬ ⑫ 0.296 19.01 12.42 8 ⑤ ⑨ 0.415 9.36 6.07 19 ⑫ ⑭ 0.129 43.37 41.92 9 ⑦ ⑤ 0.327 1.66 1.24 20 ⑭ ⑮ 0.727 43.37 41.92 10 ⑧ ⑦ 0.647 3.77 2.12 21 ⑮ ① 0 43.37 41.92 11 ⑥ ⑧ 0.349 9.39 7.34 
下载: 导出CSV
表 2 分支4风量的灵敏度和支配度
Table 2 Sensitivity and dominance of the air volume of branch 4
分支编号 灵敏度 支配度 分支编号 灵敏度 支配度 1 3.2169 124.3317 12 0.0573 5.9822 2 3.5545 119.3313 13 0.1969 6.2548 3 3.4561 146.8510 14 0.1522 11.5349 4 2.8358 15.6208 15 4.2291 56.1725 5 2.7374 40.1354 16 0.8291 41.8228 6 0.7791 78.3699 17 6.8767 152.2138 7 0.1961 13.1203 18 3.5631 90.9727 8 1.1088 23.9805 19 3.2169 124.3317 9 0.0271 1.1894 20 3.2169 124.3317 10 0.0455 5.2998 21 3.2169 124.3317 11 0.1424 29.8293
下载: 导出CSV
表 3 灵敏度
$ {d}_{4,j} $ 随$ {R}_{j} $ 的变化Table 3 Variation of sensitivity
$ {d}_{4,j} $ with$ {R}_{j} $ ${\mathit{R} }_{15}/(\rm{N}\cdot {\rm{s} }^{2}\cdot {\rm{m} }^{-8})$ $ {\mathit{d}}_{4,15} $ ${\mathit{R} }_{18}/(\rm{N}\cdot {\rm{s} }^{2}\cdot {\rm{m} }^{-8})$ $ {\mathit{d}}_{4,18} $ ${\mathit{R} }_{5}/(\rm{N}\cdot {\rm{s} }^{2}\cdot {\rm{m} }^{-8})$ $ {\mathit{d}}_{4,5} $ ${\mathit{R} }_{8}/(\rm{N}\cdot {\rm{s} }^{2}\cdot {\rm{m} }^{-8})$ $ {\mathit{d}}_{4,8} $ 0.209 4.2291 0.296 3.5631 0.197 2.7374 0.415 1.1088 0.315 3.3323 0.425 2.7233 0.265 2.3562 0.525 0.9211 0.425 2.7041 0.525 2.2890 0.425 1.7563 0.615 0.8072 0.625 1.9785 0.625 1.9658 0.615 1.3303 0.765 0.6664 0.855 1.4831 0.875 1.4339 0.825 1.0362 0.875 0.5891 1.215 1.0379 1.125 1.1134 1.125 0.7755 1.225 0.4253 1.755 0.6906 1.875 0.6391 1.875 0.4565 1.765 0.2910 2.215 0.5254 2.125 0.5537 2.125 0.3974 2.125 0.2376 3.115 0.3454 3.145 0.3474 3.125 0.2543 3.125 0.1531 4.075 0.2448 4.125 0.2484 4.125 0.1816 4.125 0.1101
下载: 导出CSV
表 4 不同算法优化结果
Table 4 Optimization results of different algorithms
算法 风量优化平
均值/(m3∙s−1)风量优化最优
解/(m3∙s−1)平均收敛
代次数平均运行
时间/sPSO算法 13.2015 13.3148 56 16.71 标准BAS算法 13.4163 13.5067 87 9.46 改进BAS
算法13.9817 14.0185 93 10.13
下载: 导出CSV
表 5 优化调节后各分支风量分配
Table 5 Air volume distribution of each branch after optimal adjustment
m3/s 分支
编号最小
需风量调节后
风量分支
编号最小
需风量调节后
风量1 41.92 42.41 12 2.59 4.99 2 16.85 21.33 13 5.46 5.98 3 19.43 21.08 14 4.26 4.26 4 6.65 14.02 15 9.94 12.05 5 7.29 7.31 16 6.71 9.26 6 11.26 13.99 17 20.38 26.07 7 6.14 7.08 18 12.42 16.34 8 6.07 6.07 19 41.92 42.41 9 1.24 1.24 20 41.92 42.41 10 2.12 4.74 21 41.92 42.41 11 7.34 8.99
下载: 导出CSV
-
[1] ZHOU Lihong,YUAN Liming,THOMAS R,et al. Determination of velocity correction factors for real-time air velocity monitoring in underground mines[J]. International Journal of Coal Science & Technology,2017,4(4):322-332.
[2] 周福宝,魏连江,夏同强,等. 矿井智能通风原理、关键技术及其初步实现[J]. 煤炭学报,2020,45(6):2225-2235. DOI: 10.13225/j.cnki.jccs.zn20.0338 ZHOU Fubao,WEI Lianjiang,XIA Tongqiang,et al. Principle,key technology and preliminary realization of mine intelligent ventilation[J]. Journal of China Coal Society,2020,45(6):2225-2235. DOI: 10.13225/j.cnki.jccs.zn20.0338
[3] 钟德云,王李管,毕林,等. 基于回路风量法的复杂矿井通风网络解算算法[J]. 煤炭学报,2015,40(2):365-370. DOI: 10.13225/j.cnki.jccs.2014.0286 ZHONG Deyun,WANG Liguan,BI Lin,et al. Algorithm of complex ventilation network solution based on circuit air-quantity method[J]. Journal of China Coal Society,2015,40(2):365-370. DOI: 10.13225/j.cnki.jccs.2014.0286
[4] KARACAN C O. Development and application of reservoir models and artificial neural networks for optimizing ventilation air requirements in development mining of coal seams[J]. International Journal of Coal Geology,2007,72(3/4):221-239.
[5] BABU V R, MAITY T, BURMAN S. Energy saving possibilities of mine ventilation fan using particle swarm optimization[C]. IEEE International Conference on Electrical, Electronics, and Optimization Techniques （ICEEOT）, Chennai, 2016: 676-681.
[6] 吴新忠,张兆龙,程健维,等. 矿井通风网络的多种群自适应粒子群算法优化研究[J]. 煤炭工程,2019,51(2):75-81. WU Xinzhong,ZHANG Zhaolong,CHENG Jianwei,et al. Optimization of multi-group self-adaptive particle swarm algorithm for mine ventilation network[J]. Coal Engineering,2019,51(2):75-81.
[7] 崔传波,蒋曙光,王凯,等. 基于风量可调度的矿井风量调节[J]. 工矿自动化,2016,42(2):39-43. CUI Chuanbo,JIANG Shuguang,WANG Kai,et al. Adjustment of mine air volume based on air volume dispatchable model[J]. Industry and Mine Automation,2016,42(2):39-43.
[8] 卢新明,尹红. 矿井通风智能化理论与技术[J]. 煤炭学报,2020,45(6):2236-2247. DOI: 10.13225/j.cnki.jccs.ZN20.0365 LU Xinming,YIN Hong. The intelligent theory and technology of mine ventilation[J]. Journal of China Coal Society,2020,45(6):2236-2247. DOI: 10.13225/j.cnki.jccs.ZN20.0365
[9] 王凯,蒋曙光,马小平,等. 瓦斯爆炸致灾通风系统实验及应急救援方法[J]. 中国矿业大学学报,2015,44(4):617-622,643. WANG Kai,JIANG Shuguang,MA Xiaoping,et al. Experimental of the ventilation systems hazard by gas explosion and the methods of emergency rescue[J]. Journal of China University of Mining & Technology,2015,44(4):617-622,643.
[10] 尹晓玉,谢贤平,李建功. 矿井通风网路的线性解算法及其程序设计[J]. 安全与环境学报,2015,15(4):69-73. DOI: 10.13637/j.issn.1009-6094.2015.04.015 YIN Xiaoyu,XIE Xianping,LI Jiangong. Linear algorithm for the mining ventilation network and program[J]. Journal of Safety and Environment,2015,15(4):69-73. DOI: 10.13637/j.issn.1009-6094.2015.04.015
[11] 丰胜成,付华. 矿井通风网络风量调节分支的优化选择[J]. 辽宁工程技术大学学报(自然科学版),2019,38(6):513-516. FENG Shengcheng,FU Hua. The optimal selection of air flow adjustment branch of mine ventilation network[J]. Journal of Liaoning Technical University (Natural Science),2019,38(6):513-516.
[12] ZHOU Hongbing,ZENG Xianlin,HONG Yiguang. Adaptive exact penaly design for constrained distributed optimization[J]. IEEE Transactions on Automatic Control,2019,64(11):4661-4667. DOI: 10.1109/TAC.2019.2902612
[13] KHAN A H,CAO X W,LI S,et al. BAS-ADAM:an ADAM based approach to improve the performance of beetle antennae search optimizer[J]. IEEE/CAA Journal of Automatica Sinica,2020,7(2):461-471. DOI: 10.1109/JAS.2020.1003048
[14] JIANG Xiangyuan, LI Shuai. BAS: beetle antennae search algorithm for optimization problems[J]. International Journal of Robotics and Control, 2018, 1（1）. DOI: 10.5430/ijrc.v1n1p1.
[15] 戴英彪. 基于拉丁超立方试验设计的事故再现结果不确定性分析[D]. 长沙: 长沙理工大学, 2011. DAI Yingbiao. Uncertainty analysis of accident recurrence results based on Latin hypercube experimental design[D]. Changsha: Changsha University of Science and Technology, 2011.
[16] FAN Chaojun,LI Sheng,LUO Mingkun,et al. Coal and gas outburst dynamic system[J]. International Journal of Mining Science and Technology,2017,27(1):49-55. DOI: 10.1016/j.ijmst.2016.11.003
-
期刊类型引用(9)
1. 刘登坤. 三维激光扫描技术在机场飞行区的应用与价值分析. 中国信息界. 2025(02): 31-33 . 
百度学术
2. 张全平，任助理，郝英豪，邓浩坤，苏士杰，袁瑞甫，方程. 基于移动式三维激光扫描技术的巷道围岩变形监测研究. 煤矿机械. 2025(05): 205-208 . 百度学术
3. 连忠文，任助理，郝英豪，杨帆，白刚，方程，袁瑞甫. 基于区域生长的非结构巷道点云去噪方法. 工矿自动化. 2024(03): 48-55 . 本站查看
4. 谢松君. 三维激光点云技术在智慧城市空间信息测绘中的应用. 中国信息界. 2024(03): 84-86 . 百度学术
5. 刘金生，张连杰. 港口常压储罐检验检测技术方法. 水道港口. 2024(04): 526-535 . 百度学术
6. 朱威霖，杨科，池小楼. 基于光纤光栅传感技术的深部巷道围岩状态精准表征. 中国矿业. 2024(10): 177-184 . 百度学术
7. 于洋，陈少青，程林. 三维激光扫描技术在储煤筒仓耐磨衬板检测中的应用. 港口航道与近海工程. 2024(05): 89-94 . 百度学术
8. 豆曙杰，李健. 长平煤矿Ⅲ53102巷三维激光扫描巷道变形监测研究. 煤. 2024(12): 49-53 . 百度学术
9. 刘敬东，李旭，郑志强，苟丙荣，韩维新，巩泽文. 激光SALM技术在煤矿巷道形变监测与支护中的应用. 矿山机械. 2024(12): 52-57 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 337
- HTML全文浏览量: 72
- PDF下载量: 46
- 被引次数: 10
