0.   引言

近年来，随着矿井智能化的普及、矿井通风技术的飞速发展^[1]，矿井智能通风系统覆盖面逐渐扩大，而该系统的建立与稳定运行对通风参数的精确度提出了更高的要求。高精度风速传感器在满足智能化系统对各项参数高准确度要求的前提下，一定程度上消除了井下测风人员的安全隐患，逐渐取代了传统的井下测风工具，成为大部分矿井测风工作的首选。井下巷道通风工作复杂且影响因素多变，尤其体现在异常风流对通风系统的影响上，其原因之一是巷道中堆放设备、杂物等，造成巷道风流局部紊乱，使附近风速测点所得数据产生偏差，不利于巷道通风环境的维持与通风网络的解算。在测风巷道中，风速传感器在不同的空间位置所测得的风速值也不尽相同，因此精准感知巷道通风参数，对异常风流巷道的测风仪器的合理空间位置进行研究，确保智能化通风系统基础信息的科学性与可靠性是至关重要的。

目前许多学者在井下巷道测风领域展开了研究。刘剑等^[2]对通风参数在湍流脉动作用下的变化进行了研究，得到可实现精准测风的理论支撑。李雪冰等^[3]通过理论与实验研究得出巷道断面任意点风速与巷道断面平均风速的关系，采用提高数据采集频率的方法，满足风速时均化要求。丁翠等^[4]通过模拟与实验结合的方式对三心拱截面巷道进行研究，提出“关键环”概念，明确了“关键环”分布的关键影响因素。杨宇等^[5]对拱形巷道断面中特征线上的风流分布进行研究，巷道中没有障碍物时，风速分布呈现左右基本对称、上下差别较大的状况。潘竞涛等^[6]对风速传感器的测量数据采用最小二乘法进行关系拟合与分析，将测量数据转换成平均风速并应用于实际生产。李亚俊等^[7]对圆形巷道断面风流分布规律进行研究，得出可通过曲线拟合对巷道断面任一点风速值计算的结论。宋莹等^[8]通过实验与模拟相结合的方法，得出巷道断面平均风速分布规律及多种截面风流的分布特征。张浪^[9]通过现场实测与模拟，得出井下测风站人员与传感器的最优位置。张士岭^[10]对风速变化规律进行研究，提出了巷道断面平均风速所在位置及测量方法。鹿广利等^[11]通过对不同条件下的巷道进行模拟与分析，得出不同转弯角度下传感器最优布设位置。张京兆等^[12]研究了入口形式对巷道定点测风位置的影响，得出巷道不同入口形式的合理测风定点位置，其中矩形巷道风流平衡时，合理测风高度变化区间为2.59～2.79 m。李虎民等^[13]通过现场实测与数值模拟相结合的方式，确定了巷道相对粗糙度对定点测风位置的影响大小。盛典^[14]通过数值模拟与分析，得出全断面风门能够有效对矿井通风系统进行调节与控制的结论。葛启发等^[15]通过数值模拟与分析，对进路式采场通风参数进行了优化，为降低局部通风机能耗提供了理论指导。郭对明等^[16]对深部掘进巷道进行相关研究，改善了井下高温环境。王春龙等^[17]通过数值模拟分析，对矿井通风降温参数优化进行了分析研究，为矿井深部通风降温参数设置提供了参考。上述学者采用数值模拟软件与理论分析结合实验验证的方法研究了井下巷道风流运移规律、优化了矿井通风相关参数，但未涉及巷道中存在不同形状的障碍物时风流运移规律的探索。

本文以中国华电集团有限公司小纪汗煤矿11218回风巷为例，对井下巷道中障碍物不同位置与不同尺寸对风速的影响展开研究，利用Fluent数值模拟软件模拟井下巷道中不同障碍物的大小和障碍物放置位置对风速的影响，并找到受障碍物放置位置影响最小的监测位置。

1.   模型建立与模拟参数设置

1.1   模型建立

在不考虑巷道壁面发生局部形变的情况下，采用小纪汗煤矿矩形截面11218回风巷的测量尺寸与风速测量平均值进行模拟计算，其数据见表1。

表  1  11218回风巷截面参数信息

Table  1.  11218 return air roadway parameter information

高度/m	宽度/m	周长/m	断面积/m2	进风口风速/（m·s−1）
3.16	5.54	17.4	17.51	2

下载: 导出CSV 

| 显示表格

ANSYS Fluent软件可用于模拟各种复杂流场，湍流模型具有研究范围广、精度高的显著特征，适用于研究矿山井下风流场^[18-19]。本文建立6种障碍物模型，其长、宽、高分别用l，b，h表示，障碍物信息见表2。设障碍物与两帮之间的较近距离（以下简称间距）L=0，0.5，1 m；粗糙度常数为0.8。

表  2  障碍物信息

Table  2.  Obstacle information

	l/m	b/m	h/m
障碍物1	1	1	1
障碍物2	2	1	1
障碍物3	1	1	0.5
障碍物4	2	1	0.5
障碍物5	1	0.5	1
障碍物6	2	0.5	1

下载: 导出CSV 

| 显示表格

以表1、表2数据为基础，采用Design Modeler软件建立不同参数的障碍物模型，分别对应6种障碍物模型，建立了6组巷道模型，将每组巷道模型分为3个子模型（L=0，0.5，1 m），共计18个巷道子模型。L=0.5 m时的各模型关键部位如图1所示。

以x轴为巷道走向，y轴为巷道宽度，z轴为高度，当L=0.5 m时，巷道模型的俯视图如图2所示。根据障碍物模型与间距L，对11218回风巷距上游端口10 m处放置障碍物的情况进行分类。由于该巷道截面具有对称结构特点，可由巷道一侧的结果推导出另一侧的相关结论。

图  1  模型关键部位

Figure  1.  Key parts of the model

下载: 全尺寸图片 幻灯片

图  2  L=0.5 m时巷道模型的俯视图

Figure  2.  Top view plan of the roadway model at L=0.5 m

下载: 全尺寸图片 幻灯片

1.2   模拟参数设置

采用ANSYS Fluent软件建立18个巷道子模型，巷道子模型较为规整，采用cartesianMesh方法进行网格划分，网格单元平均尺寸为0.2 m，最大、最小尺寸分别为0.25，0.1 m，经网格划分的独立性考核确定其合理性；以L=0.5 m时模型2为例，该模型划分共计约116万个节点、100万个网格单元，平均单元质量为0.998。其局部网格划分如图3所示。

图  3  L=0.5 m时模型2局部网格划分

Figure  3.  Model 2 local grid division at L = 0.5 m

下载: 全尺寸图片 幻灯片

井下巷道风流类型多为湍流^[20-21]，在巷道放置障碍物后将导致风流更加紊乱，但风流依旧为湍流类型。本文采用只考虑动量传输，不涉及能量交换的Standard $ k-{\varepsilon } $湍流模型作为流体运动模型，采用双方程湍流模型压力耦合的半隐式方法（Semi-Implicit Method for Pressure Linked Equations，SIMPLE）对流体运动进行仿真，收敛标准除连续性方程为1×10⁻⁴外，其他均为1×10⁻⁷。以空气为模拟风流，密度取实测平均密度1.07 kg/m³，通过插值法结合温度数据确定巷道空气黏性系数为1.818 5×10⁻⁵ kg/（m·s），使得模拟结果更符合11218回风巷的实际通风状况。

边界条件设置：进风口采用固定速度进风（ $ 2\;\mathrm{m}/\mathrm{s} $），出风口采用自由出流；表面粗糙度常数为0.5，其他巷道壁面粗糙度常数取0.7，障碍物粗糙度常数取0.8，壁面类型均为标准壁面函数，其余设置及参数采用软件系统的推荐值或默认值。

2.   障碍物对巷道合理风速分布区间的定量分析

为研究障碍物在不同条件下对11218回风巷风流场的影响，选取障碍物种类、间距大小、入口风速3个变量结合控制变量法展开相关分析，所有巷道的模型尺寸、环境条件、网格划分所采用的方法、求解器参数设置等均保持一致。数据取自巷道断面宽度中心位置（y=2.77 m）、高度为2.7 m的所有计算单元，每0.2 m设置1个分析节点，进风口处风流场云图及分析节点如图4所示。

图  4  进风口风流场云图及分析节点

Figure  4.  Inlet airflow field cloud diagram and analysis node

下载: 全尺寸图片 幻灯片

煤矿企业规定人工测风误差不大于5%，为满足井下测风工作精度需求，指导矿井智能测风工作的开展，本文以3%为风速传感器测风精度标准，开展以下分析工作。当进风口风速为2 m/s时，传感器测风误差为1.94～2.06 m/s，称测风精度内的风速数值为合理风速，风速在测风精度内的风流为合理风流，合理风流处即为风速传感器最优布设区域。由文献[12]可知，矩形巷道断面在无障碍物且风流充分发展条件下，当x＞115 m时，在2.59 m＜z＜2.79 m位置布设风速传感器较为合理。因此本文在 $ x=120\;\mathrm{m} $截面的中垂线 $ {\textit{z}}=2.7\;\mathrm{m} $处布设风速传感器；取障碍物迎风侧前1 m、障碍物迎风侧、障碍物背风侧、障碍物背风侧后1 m、风速传感器截面（ $ x= 120\;\mathrm{m} $）为定量分析内容，根据x坐标升序方式分别命名截面，为截面1—截面5。

2.1   L=0时障碍物对巷道风流场的影响

当 $ L=0 $时，取各个模型所对应的风流场云图做对比，结果如图5所示。可看出截面1、截面3、截面4的风流场云图中，合理风流所呈现形状大致为弧状，且与各自对应的障碍物矩形不符。在截面1中，合理风流分布面积随着障碍物宽度减小而减小。

图  5  L=0时所有模型的5种截面位置所对应的风流场云图

Figure  5.  Cloud plot of wind fields corresponding to the five cross-sectional positions of all models at L =0 m

下载: 全尺寸图片 幻灯片

结合图4与图5中的截面5可看出，顶板位置中垂线的合理风流分布于2.59～2.78 m处。

2.2   L=0.5 m时障碍物对巷道风流场的影响

当 $ L=0.5\;\mathrm{m} $时，取各个模型所对应的风流场云图做对比，结果如图6所示。可看出截面1中合理风流分布面积随障碍物体积减小而减小，截面3中当宽度 $ b=0.5 $ m时，模型5、模型6可体现障碍物大致形状。结合图4与图6中的截面5可知，顶板位置中垂线的合理风流分布在2.59～2.80 m处，其中模型1、模型2风流发展较为缓慢，尚未形成圆角矩形的形状。

图  6  L=0.5 m时所有模型的5种截面位置所对应的风流场云图

Figure  6.  Cloud plot of wind fields corresponding to the five cross-sectional positions of all models at L =0.5 m

下载: 全尺寸图片 幻灯片

2.3   L=1 m时障碍物对巷道风流场的影响

当 $ L=1\;\mathrm{m} $时，取各个模型所对应的风流场云图做对比，结果如图7所示。可看出随着障碍物宽度减小，截面1中合理风流横向变化不明显；截面3中，宽度 $ b=0.5 $ m时，模型5、模型6可体现障碍物大致形状。结合图4与图7中的截面5可知，顶板位置中垂线的合理风流分布于2.61～2.78 m处。

图  7  L=1 m时所有模型的5种截面位置所对应的风流场云图

Figure  7.  Cloud plot of wind fields corresponding to the five cross-sectional positions of all models at L=1 m

下载: 全尺寸图片 幻灯片

结合图4−图7可知，在3种间距中，各模型于断面直角处存在微小部分的合理风速区域，其面积在 $ L=0.5\;\mathrm{m} $时最大， $L=1\;\mathrm{m} $时次之， $L=0$时最小，随着间距L的增加，截面1中合理风流分布区域逐渐增加；模型3、模型4的风流场云图中，合理风流横向分布较广， $ L=1\;\mathrm{m} $时较为明显。截面1、截面3、截面4的风流场云图中，合理风流于横向扩散较快；在放置障碍物的矩形巷道中，合理风流在图4中均匀分布，截面1—4的直角处微量分布，截面5为圆角矩形分布，故合理风流随x坐标增加呈均匀分布→直角处微量分布→空心圆角矩形分布的规律。两帮附近风流发展不对称，障碍物体积越大，风流发展至平衡所需时间越长；巷道放置障碍物一侧的合理风流分布区域更广，即传感器最优布设区域更大；结合表2，对比模型1、模型2和模型5、模型6可知，障碍物宽度越大，此现象越明显。

3.   L不同对风流运移规律的影响分析

以同种障碍物、不同间距为模拟条件的计算结果分析11218回风巷风流发展情况，取巷道断面所有计算节点风速数据，各个模型对比数据如图8所示。可看出随着x的增加，各模型平均风速大小均呈增大−减小−增大−平衡的状态。

图  8  不同模型平均风速对比

Figure  8.  Comparison of the average wind speed of different models

下载: 全尺寸图片 幻灯片

为方便对比各模型数据差异，降低分析研究的复杂程度，现以风速变化为依据，将风流发展分为4个阶段。

在x为0～10 m时，各模型风速呈上升状态，x接近10 m区域时，风速急剧增加至最大值，这是由于风流发展受到阻碍，即此时风流与障碍物相接触，致使通风面积急剧减少，各模型风速在x接近10 m时急剧上升，将此区间命名为第1阶段；当x＞10 m后，风速降低幅度逐渐减小，在x为60～90 m内各模型风速先后降低至谷值，这是由于风流经过障碍物后通风面积突然增大，在障碍物后面形成涡流区域，进而使风速急剧降低，随后风速降幅逐渐减小，将此区间命名为第2阶段；风速在降低至最小值后缓慢增加，各模型风流在x为180～220 m内先后发展至平衡，将风速在此增加的区间命名为第3阶段；风流发展平衡时x坐标至巷道出风口区间命名为第4阶段。

在第1阶段，模型1， $ L=0$时风速峰值已超过 $ 2.225\;\mathrm{m}/\mathrm{s} $；模型3与模型5、模型4与模型6峰值均接近，且模型3与模型5、模型4与模型6所对应的障碍物体积相同。由此可见，障碍物体积相同，风速峰值大致相同。

在第2阶段，风速到达谷值时x坐标与障碍物长度呈负相关；在模型1中， $ L=0.5\;\mathrm{m} $时风速谷值最大，模型4中， $L=0$时风速谷值最小。

在第3阶段，L不同时，模型1、模型2风速差别较大；模型3、模型4风速大致相同；障碍物高度 $ h=1\;\mathrm{m} $时， $ L=0$时风速最大， $L=1\;\mathrm{m} $时次之， $L=0.5\;\mathrm{m} $时最小。

在第4阶段，障碍物长度与风流稳定位置呈负相关；风流发展平衡时，模型1、模型2的风速数值差最大，模型5、模型6次之，模型3、模型4的风流变化幅度最小，当h=0.5 m时，风速集中于 $ 2.05\;\mathrm{m}/\mathrm{s} $。结合表2可知障碍物竖放或增加宽度对风流影响较大。

4.   不同模型对风流运移规律的影响分析

采用相同的数据采集方式及阶段式分析方法，对同间距L条件下，不同模型对风流运移规律的影响进行研究。其对比结果如图9所示。

图  9  同间距风速对比

Figure  9.  Wind speed comparison at the same distance

下载: 全尺寸图片 幻灯片

在第1阶段，各模型在风流经过障碍物时，在 $ L=0 $时风速波动较大，其他间距风速波动大致相同。

在第2阶段，各模型风速降低幅度大致相同， $ L=0.5\;\mathrm{m} $时，模型3—模型6相对模型1、模型2提前 $ 20\;\mathrm{m} $降低至风速谷值；模型1、模型2风速谷值远大于其他模型，风速谷值降序为模型1、模型2、模型5、模型6、模型3、模型4；障碍物长度相同时，风速谷值与障碍物高度、宽度呈正相关，障碍物高度与宽度相同时，风速谷值与障碍物长度呈负相关。

在第3阶段，模型1、模型2风速相对其他模型增幅较小，在 $ L=0.5\;\mathrm{m} $时最为明显，即障碍物体积与本阶段风速增幅呈负相关。

在第4阶段，风流发展平衡，模型2风速最接近初始风速； $ L=0 $， $1\;\mathrm{m} $时风速集中于 $ 2.05\;\mathrm{m}/\mathrm{s}$； $L=0.5\;\mathrm{m} $时，风速不集中，但风速相对最小，接近于初始风速；因此风流发展稳定时， $ L=0.5\;\mathrm{m}$时风速可靠性最高， $L=1\;\mathrm{m} $时次之， $L=0$时可靠性相对最低。

5.   风流运移规律的风速普适性分析

煤矿安全规程规定，采区进回风巷风速应处于0.25～6 m/s，本文以1 $ \mathrm{m}/\mathrm{s} $为最低风速，1 $ \mathrm{m}/\mathrm{s} $为速度增量，5 $ \mathrm{m}/\mathrm{s} $为最高风速，进行相关研究与分析。由于 $ L=0.5\;\mathrm{m} $时风速可靠性最高，所以在 $ L=0.5\;\mathrm{m} $时，不同风速条件下，展开风流运移规律的风速普适性研究是具有代表性的。采用相同的数据采集方式，定义初始风速为 $ {v}_{1} $，巷道平均模拟风速为 $ {v}_{2} $，风速变化率 $ \varphi =({{{v}_{2}-v}_{1}})/{{v}_{1}} $，其对比结果如图10所示。

图  10  同模型风速变化率对比

Figure  10.  Comparison of wind speed change rate with the same model

下载: 全尺寸图片 幻灯片

上述多种条件下，风速变化率均处于上升−下降−再上升−平衡的4个阶段，故采用阶段式分析方法，其分析结果如下。

在第Ⅰ阶段，风速变化率增幅与初始风速大小呈正相关，随着初始风速增大，风速变化率波动幅度逐渐减低；风速变化率峰值与障碍物参数呈正相关，其中宽度（b）影响最大，高度（h）影响次之，长度（l）影响最小。

在第Ⅱ阶段，风速变化率谷值与初始风速呈正相关，初始风速越小，风速变化率波动幅度越明显。风速变化率谷值大小与障碍物长度（l）呈负相关，与障碍物其他参数关联不明显。

在第Ⅲ阶段，模型1、模型2风速变化率相对其他模型增幅较为平缓；除 $ 1\;\mathrm{m}/\mathrm{s} $的风速变化率曲线外，其他曲线于巷道走向长度90～120 m内相继交于一点，障碍物宽度（b）与高度（h）越大，该交点坐标值越大；在曲线相交前，初始风速越大，则风速变化率越大，在曲线相交后，风速越小，则风速变化率越大。

在第Ⅳ阶段，风流发展平衡时，风流变化率均小于3%；模型1、模型2风速变化率较其他模型小，其中模型2风速变化率最小，即该模型最接近初始风速。

6.   结论

1） 在放置障碍物的矩形巷道中，在 $ x=120\;\mathrm{m} $，2.59 m＜z＜2.79 m处布设传感器条件下，截面5中传感器最优布设区域在该截面中垂线高度2.59～2.80 m处，两帮附近风流发展不对称，放置障碍物侧合理风流分布区域更广，障碍物宽度越大，该现象越明显；最优布设传感器区域随x坐标增加呈均匀分布−截面直角处微量分布−空心圆角矩形分布的规律。

2） 障碍物体积相同时，风速峰值大致相同；风速谷值与障碍物高度、宽度呈正相关，障碍物高度与宽度相同时，风速谷值与障碍物长度呈负相关。 $ h=0.5\;\mathrm{m} $时风速变化幅度小，平衡时风速集中于 $ 2.05\;\mathrm{m}/\mathrm{s} $。

3） 在多种风速模拟条件下，第Ⅳ阶段中，模型2的风流变化率仍为各个模型的最小值，与单一风速结论形成对应。故在6类障碍物模型中，模型2对应的障碍物在 $ L=0.5\;\mathrm{m} $时对回风巷风流运移规律影响较小的结论具有风速普适性，在此类条件下可确保在巷道中存在障碍物时，控制其对风速监测的影响最小。