大断面巷道风量精准测量数值模拟研究.docxVIP

? ? 大断面巷道风量精准测量数值模拟研究 ? ? 杜 斌，朱 蕾 （1.中煤科工集团沈阳研究院有限公司，辽宁 抚顺 113122；2.煤矿安全技术国家重点实验室，辽宁 抚顺 113122） 我国新开发建成投产的现代化矿井具有资源赋存条件好、井型大、产量高、自动化智能化程度高等特点。与之相对应，井下巷道的断面往往也比较大。在我国西北地区的一些大型生产矿井中，巷道断面可以达到20.5m2，井下车辆可以实现“双向两车道”行驶。我国现行的风量测定方法中，主要使用“走线法”。但是当测风员在井下大断面巷道测风地点进行测风作业时，由于身高及臂长等因素的限制，测量高度往往偏低，达不到操作规范的要求；同时，空气作为1种流体，由于黏滞力的存在，巷道风流的风速往往是中间轴心部位流速大，巷道壁面（边界层）流速小；导致了大断面巷道风量测量长期存在测量值与真实值相差偏大的问题。但是，通风系统作为矿井5大系统之一，肩负着通风供氧排尘、稀释有毒有害气体的重大任务，同时风量是矿井通风技术管理人员的眼睛，实时掌握矿井通风系统的风量分配情况是通风技术管理的关键。风量充足，就可以稀释瓦斯，降低瓦斯积聚、超限事故发生的概率；还能降低粉尘以及其他有毒有害气体的体积分数等。矿井智能化通风系统[1]建设也对大断面巷道风量精准测量技术提出了更高的要求；大断面巷道风量精准测量的需求，不仅出现在煤炭开采行业，同时亦存在于大型隧道建设施工、水电站水流通道建设施工等行业。因此，大断面巷道的风量精准测量是一个跨行业的共性技术难题，亟需得到妥善解决。 对于大断面巷道风量精准测量的技术难题，国内大量专家学者开展了研究。霍曼[2]根据切贝切夫法规则，采用在拱形巷道的测风断面设置了36个风速传感器的方法对风量进行精准测量，但是由于测试支架的黏滞作用和流体的不可压缩性会导致风速在测风断面处发生变化，从而使得测试结果产生较大的误差；为了避免测试工作给巷道内风速分布带来的扰动，刘明浩[3]采用激光多普勒测速仪（LDA）对巷道风速进行非接触式测量，由于激光汇聚点就是测量探头，所以排除了测量工作对流场的干扰；丁翠等[4-5]借助数值模拟方法对拱形和梯形巷道内风流的分布规律进行研究，同时建立同尺寸巷道通风实验系统对模拟结果进行对比验证，分析巷道断面上平均风速点的分布规律，明确了巷道中平均风速点的位置仅与巷道截面的几何特征有关，而与通风风速无关。也就是说在巷道断面的几何参数保持不变的情况下，不论巷道风量如何变化，巷道中风速与平均风速相等的位置是不会发生改变的，这对精准测量大断面巷道风量具有很好的启发意义。测试时可以在巷道断面的平均风速处使用非接触式高精度风速传感器测量风速[6]，得到的风速就是平均风速，然互配合使用激光断面仪测量巷道断面的面积，二者之积就是流经巷道的风量。但是文献[4]只给出了拱形断面巷道中平均风速点的轨迹方程，对于矩形巷道或者非规则形状巷道中平均风速点的分布规律尚未说明，即没有明确给出风速传感器的安设位置。鹿广利[7]开展了关于巷道风流状态及控制方程的理论研究，推导出了巷道断面风速分布公式： 式中：v为巷道断面的平均风速，m/s；a为巷道摩擦阻力系数，（kg·s2）/m4；r0为巷道半径，m；r为巷道某点到轴心的距离，m。 式（1）表明只要确定巷道断面内某一点的风速值和这个测量点距离巷道断面中心的距离，便可以确定巷道断面平均风速的大小。当巷道断面为规则几何图形时，可以准确快速地确定测点巷道断面中心的距离，但是当巷道断面为不规则几何图形时，难以确定巷道中心的位置。 受文献[4]的启发，测量大断面巷道的风量可以使用巷道断面中的平均风速与巷道断面面积的乘积。巷道断面面积可以使用激光断面仪准确测得，而巷道断面中的平均风速却不易获得。文献[4]只指出大断面巷道中风速等于平均风速的位置不会随风流风速的改变而发生变化，只与巷道断面的几何形状有关，但尚未给出平均风速位置的确定方法，所以无法指导现场人员使用非接触式高精度风速传感器在平均风速位置测量大断面巷道的风量。为此，使用COMSOL Multiphysics研究大断面巷道中风流风速的分布特性，模拟大断面巷道不同入口风速条件下巷道断面中平均风速的位置，验证文献［3］所述的“关键环”，然后利用后处理功能确定大断面巷道中瞬时风速等于平均风速的位置，并将计算模型进行封装，开发出可以独立运行的APP。 1 大断面巷道风量精准测量 1.1 巷道几何模型 以煤矿现场常见的矩形巷道和半圆拱形巷道为例开展研究，矩形巷道高4.6m、宽6.9m、模拟计算长度20m，半圆拱形巷道半径3.45m、高4.8m、模拟计算长度20m。建立的矩形断面巷道几何模型如图1，半圆拱形巷道几何模型如图2。 图1 矩形断面巷道几何模型Fig.1 Geometric model of r