4 初始及边界条件

　　4.1 火源的处理

　　不考虑矿石自燃的具体燃烧过程，只简单地将火源处理成一高温定温区，核心温度为150℃，不断地向周围释放热量和烟气。高温烟流在巷道内流动过程中不再发生化学反应，产生的高温烟气中主要考虑组分SO 2，浓度为77.4×1O -6。

　　4.2边界条件

　　入口处风流速度V in=1.5 m/s，温度20℃，紊流动能

　　动能耗散率

　　系数C D =1.0。

　　出口边界条件采用局部单向化假设，即认为出口截面上的节点对第一个内节点无影响，离散化方程系数aE=0。

　　在巷道壁面，假定烟流不可渗透，气密性好，在壁面处扩散通量为0，粗糙度取0.03 m，壁面温度等于巷道围岩冷却带的温度为20℃。

　　5 工程实例

　　5.1 工程背景

　　华锡集团公司铜坑矿细脉带矿体属高温一中温热液交代型矿床，矿体含硫量高达13.3%，如暴露在通风好的情况下极易氧化而引起自燃，产生的高温以及SO 2等有毒气体，对井下生产及地面环境造成了严重的影响。

　　近年来，随着开采活动的继续、民采的盗采矿柱等，隔火矿柱和密闭墙遭到了严重的破坏，导致火区“活化”，矿岩重新发生自燃而引起火灾，产生的SO 2气体窜到井下作业面，引起井下作业人员的中毒或窒息，严重威胁着矿山的安全生产。因此，对地下火和毒气蔓延的控制已刻不容缓。

　　5.2 数值计算与安全评价

　　选取一中段水平巷道为模拟对象，长度为60 m，高度为2.8 m，自燃发火区位于巷道顶板处，距巷道入口处的距离为15 m。烟气流动的速度矢量图、温度分布图如图1、图2所示。

　　从图1中可以看出，*近火源区域的风流速度要比远离火源区域的风流速度慢，这是因为矿井发生火灾后，巷道内风流受到高温的加热作用，密度发生变化，在重力场的作用下诱发对流，此时的风流流动主要由浮升力和主要通风机的风压所控制。

　　在*近火源的区域，气体受热后膨胀，能对气体的流动产生一定的阻止作用；而由于入口处的风速较大，烟气上浮的趋势明显地被高风速所抑制，所以烟流在巷道中并没有产生逆流现象，而是沿着巷道的轴向方向蔓延。

　　由图2的温度分布图可以看出火源处的温度随着风流不断向巷道轴向方向扩散，并且逐渐地降低，温度曲线呈下滑状。断面烟流温度分布也不均匀，*近巷道顶板的烟流温度高，而*近底板或壁面烟流温度低。温度场分布的情况与文献中的实际巷道火灾实验结果一致，因而用该数值模拟方法所得出的结果是可*的。

　　由图3可知，在入口至火源的区域内的PMV较低，人体感觉适中，并未出现闷热的现象。其原因是入口处的风流速度较大，抑制了高温的扩散。而火源至出口之间的巷道内热舒适度值都很高，说明了烟流气体均被排至该段巷道中去了，而前段巷道中的风流则为新鲜风流，通过图4中的空气龄指标也表明了上述现象。

　　经计算，整个模拟区中的预计热舒适指标值为：min=-2.3148，max=3，mean=1.12。

　　经计算，整个模拟区中的预计不满意百分比值为：min=5.00，max=99.12，mean=84.73。

　　从图5可以看出，PPD先逐渐变小，后又逐渐增大，笔者认为，入口处的风速较大，在离火源区较远处出现了人体感觉“偏冷”的现象，而随着气流不断接近火源区，吸收了火源区所扩散出来的温度，气流得到了升温，人体感觉逐渐地由冷至适中，所以在一定范围内PPD值得到了减小。

　　在火源区附近，由于各处的温度不同，浮升力大小不均匀，驱动着气体上浮流动，即出现巷道下部的烟气流动速度要比上部快，上部的高温烟气得不到有效的扩散，在巷道上部区域出现了“聚集”，从而使PPD值逐渐的增大，并且下部区域的PPD值始终比上部的要低。

　　6 结论

　　1)得出了自燃火灾巷道中高温烟流的速度场、温度场分布规律。为矿山制定出行之有效的安全预案及准确的救援决策信息，包括高温、有毒有害气体大量涌出时的安全应急方案、人员撤退路钱的选择、风机的布置、避险硐室位置的确定等提供了技术依据。

　　2)当巷道入口处的风速Vin=1.5m/s时，火灾巷道中的烟气并没有出现逆流现象，说明当矿石自燃发生火灾时，采用该风速可以满足巷道内通风排烟的要求。

　　3)采用数值模拟的方法，可以帮助研究人员减少盲目性，增强目的性，并且可以节约时间和经费，可有效地预测火灾蔓延趋势及流场分布情况，对矿山的安全生产有很强的指导意义。

　　4)在建立模型中，对边界条件作了一定的假设和初始化，作为进一步研究，可考虑加入更复杂的边界条件，使之与工程实际更接近，以便能得出更为精确的模拟结果。