张嘉旺 陈忠海 刘龙飞 马志远

(河北建筑工程学院，河北 张家口 075000)

0 引 言

随着社会对矿物的需求不断加大，矿井的开采速度和深度随之增加，地温随之升高，用于开采综掘面的机械化程度也不断提高，机械散热也越来越大.在热害问题越来越严重的同时，生产现场的产尘量也急剧增加，这两个因素不仅危害作业人员的健康，还会大大降低劳动生产率，严重影响着矿井的安全生产.肖俊峰[1]等人通过建立单风筒长压短抽通风方式的DPM模型，研究了风口相对位置以及抽压比对巷道内粉尘运移的影响，得到了布置风筒的最佳位置和抽压比.王汉青[2]等人利用Fluent模拟了两种风筒的送风效果，总结出0.7m的送风口、6m/s的送风速度以及20℃的送风温度能提供更好的人体舒适条件.解彬[3]等人通过改变送风风速、风温等因素对单风筒和双风筒两种送风方式对巷道的降温进行模拟，总结出双风筒降温效果比单风筒降温效果好.周智勇[4]等人选取了四种通风布置参数，分别对应三种模拟方案对巷道内的粉尘分布进行了研究，得到了该巷道最理想排尘的通风参数.下面以鲍店煤矿7309综掘面为例，将降温与粉尘分布结合讨论，利用Fluent模拟探讨单风筒与双风筒两种通风方式下掘进工作面的降温效果与粉尘分布.

1 实验地点概况

图1 掘进巷道断面

鲍店煤矿7309综掘面采用单风筒压入式通风，防尘措施为掘进机内外喷雾，但除尘效果较差.瓦斯绝对涌出量：0.23 m3/min，二氧化碳绝对涌出量：0.54 m3/min，为满足通风需要，压入风量为425.6 m3/min，除风量确定为250-450 m3/min，工作面防尘措施主要为掘进机内外喷雾[5]，风筒布置如图1所示.

2 数学模型建立

由于实际施工现场环境更为复杂，风流属于湍流运动，所以通过数值模拟进行研究时需要做出如下假设：(1)模拟中的气体均为不可压缩气体，忽略重力，不考虑黏性力做功，巷道壁面和风筒的物性均视为常数；(2)气流在掘进巷道的流动遵守质量守恒定律、动量守恒定律和能量守恒定律.

(1)连续性方程

(1)

(2)动量守恒方程

(2)

(3)

(4)

(3)能量守恒方程

(5)

(4)K方程

(6)

(5)ε方程

(7)

(8)

式中：k为湍流动能，m2/s2；ε为湍流动能耗散率，m2/s3；t为时间，s；v为层流粘度系数，Pa,s；p为修正时均压力，Pa；ui为速度分量，m/s；xi,j为坐标分量，m.

3 建立物理模型

3.1 建立几何模型

根据文献[5]中数据，使用Gambit绘制横断面宽4.6 m、高为3 m、长度为40 m的长方体形巷道.压入风筒与除尘风筒均为圆柱体，压入风筒吊装于巷道一侧，风筒中心距巷道顶板与侧壁均为0.6 m，风筒直径为0.8 m，进风口距掘进工作面10 m；除尘风筒中心高1.5 m，直径0.6 m，出风口距掘进断面5 m.

3.2 边界条件及计算方法

模型入口边界条件(Velocity Inlet)：入口温度T=20 ℃，湍流强度为6，水力直径为0.8；单进风筒入口风速V=14.2 m/s，双进风筒的入口风速均为7.1 m/s.出口边界条件(Velocity Inlet)：风速V=-10.8 m/s，湍流强度为6，水力直径为0.7.壁面温度为40 ℃.

巷道内空气的密度为1.225 kg/m·s，温度为300 k，鉴于粉尘颗粒直径较小，随动性良好，DPM模型采用离散相的非耦合法计算，颗粒沿X方向扩散速度为-1 m/s,最小直径为1×10-5，最大直径为5×10-5，密度为1200 kg/m3，颗粒分布采用Rosin-Rammler分布.紊流模型使用标准k-ε双方程模型，同时加强壁面处理(Thermal Effects),采用隐式分离三维稳定流求解器.

4 模拟结果与讨论

根据上文简化的矿道模型与模拟参数的设定，在环境参数不变，送风量一定的工况下，利用Fluent模拟高温掘进面在压入式通风方式下的三维流场,分析对比单风筒与双风筒两种压风方式下风速、温度及颗粒分布的不同.

4.1 通风方式与降温效果分析

图2图3分别为压抽比1.7、距离掘进面0-15米内两种不同通风方式的速度和温度的切片云图.受限于巷道的封闭环境，由进风筒喷射出的射流属于有限贴壁射流，主流区的射流容易受到壁面的影响.对比两种通风方式可以发现，射流流速大的卷吸作用更强，对该有限空间内空气的影响也更大，更容易在局部巷道形成涡流，双风筒通风方式风流分布更均匀，射流到达掘进面的风速只有0.3 m/s.由于单风筒通风射流初始速度大，射流与空气进行充分热湿交换的时间不足，而且到达掘进面时的风速仍然可以达到2.4 m/s，直到射流冲击掘进面后折返于巷道左侧一定距离后，射流在速度逐渐衰减的过程中与空气进行热湿交换.

图2 单风筒(上)、双风筒(下)速度切片云图 图3 单风筒(上)、双风筒(下)温度切片云图

图4为两种通风方式在某一直线上的压力(a)、风速(b)、温度(c)、涡流强度(d)对比，x轴为距巷道口的距离(巷道长40 m)，该直线选取掘进机操作员的大概位置：据左侧墙壁1.5 m、高1.5 m，范围选取距离掘进面0-15 m区域.结合图三中(a)(b)(d)图可知，单风筒通风方式下，在距离掘进面6.5 m处的气流不仅速度回升而且局部产生涡流，不利于整个工作区域的整体降温.(c)图也非常直观的显示出双风筒通风方式下，工作区域的温度平均要低2℃.相比之下，双风筒通风方式在这一区域的气流稳定，对操作员的干扰小.

图4 压力、速度、温度、涡流强度曲线图

4.2 通风方式与粉尘分布规律分析

图5图6为单风筒和双风筒通风方式下巷道内的粉尘分布情况.从该三维云图的侧面可以看出，单风筒通风方式下，巷道内的粉尘分布均匀，原因是在较小的空间内射流受巷道壁面的影响大，过大的风速对掘进面散发的粉尘造成冲击的同时又将粉尘卷起，使粉尘几乎均匀的分布在工作区，不利于人员作业及除尘.相比之下双风筒通风的射流虽然也存在涡流，但因其射流速度小，所以涡流强度较低，该通风方式不仅对掘进面冲击小还引导粉尘集中在巷道下方排出，有利于除尘的同时还为工作人员创造了良好的工作环境.

图5 单风筒通风粉尘分布图

图6 双风筒通风粉尘分布图

5 结 论

模拟基于长压短抽式通风系统，在确定掘进面通风量与进出口风筒直径的条件下，对比分析了单风筒通风与双风筒通风两种方式对高温掘进面的降温效果.结合DPM模型模拟粉尘质量浓度运移的结果，得出几点结论：

1)为满足通风量需求，单风筒通风射流速度过大，射流冲击掘进面后回流，在测定的范围内距巷道6 m左右的位置，压力、风速、涡流强度变化幅度大，使空气不能充分进行热湿交换.在其他因素不改变的情况下，双风筒通风适合应用在该巷道，其射流不会强烈冲击掘进面，所以在测定范围内气流风速平稳，温度合适，压力变化稳定，回风流畅.巷道温度满足《煤矿安全规程》规定掘进面风速大于0.25 m/s小于4 m/s、温度不超过26 ℃的要求.

2)当压抽比固定时，单风筒气流到达掘进面仍有一定的速度，紊乱的气流容易在狭小的空间形成涡流，使粉尘在工作区扩散.采用双风筒通风时，粉尘主要沉降区位于掘进面1.5 m内，气流回流能够将粉尘控制在巷道下方.分析结果表明，在该掘进面采用双风筒通风能够有效降低巷道温度，并能大大降低掘进机司机工作区域的粉尘浓度.