火源位置对井巷外因火灾参数影响的数值模拟研究

2023-10-07张军亮朱云飞张学亮

煤矿安全 2023年9期

张军亮，朱云飞，张学亮

（1.山西工程技术学院矿业工程系，山西阳泉 045000；2.中煤华晋集团晋城能源有限公司，山西晋城 048200）

矿井火灾是煤矿常见的事故灾害之一，由于矿井巷道属于受限空间，矿井火灾尤其是外因火灾具有火势发展快、烟流蔓延迅速的特征，往往会导致人员受困，严重者甚至造成人员死亡，因此许多学者对矿井巷道中的火灾参数规律进行了深入研究。在巷道火灾温度分布规律上，Kurioka等[1]对巷道内顶板下的最大温度的变化规律进行了研究，得出了最大温度与火源功率曾指数关系，LI 等[2]和YAO 等[3]基于羽流理论对其经验公式进行了进一步修正；李晴等[4]基于理论分析，建立了火灾烟流温度分布的预测模型；刘雨晴等[5]分析了纵向通风下不同巷道倾角火灾烟气的蔓延特性和温度分布。在巷道火灾烟气分布规律上，薛彦平[6]对比了不同通风方式下工作面巷道烟流逆退距离规律；郝海清等[7]基于Froude 公式，提出了烟气逆流临界风速计算公式，并对数值模拟分析了羊场湾煤矿通风网络的排烟风量效果；贾静等[8]通过数值模拟分析了火灾烟气的逆退现象。在火源位置方面，张洪杰等[9]分析了火源与掘进面不同距离下的火灾参数分布规律；吴佳平等[10]对比不同封闭墙与火源距离分析了火区内瓦斯气体的爆炸危险性。上述研究文献均假设火源位于巷道底板的情况，在实际矿井中，巷道在发生冒顶等事故时，矿井往往采取复合材料或木背板填充[11]，因此，一些底部火灾可能引燃顶部的背板，同样顶部背板等可燃物燃烧后也极可能掉落至地面继续燃烧，造成双火源现象，这势必会造成火烟参数分布的变化，从而影响火烟控制措施的制定，然而现有文献鲜有研究。为此，为了重点研究巷道垂向断面不同火源位置对火灾参数的影响，选取里必煤矿某巷道实际生产通风工况，分别对顶板单一火源、底板单一火源和同时发生在顶、底板的双火源3 种工况下巷道内烟雾、温度、CO 参数分布进行了探讨分析，为井巷相应火灾应急救援措施的制定提供参考。

1 模型构建

1.1 模型及参数设定

研究选用FDS（Fire Dynamic Simulation）火灾模拟软件，其适用于模拟低马赫数火灾烟流和热量传输情况，因此常被用来分析隧道及巷道火灾。数值模型依照里必煤矿某段直巷道进行建模，选取的巷道长、宽、高依次为100、4、3 m。火源位于巷道轴向和横向中央，尺寸为1 m×1 m。巷道内传感器设置如图1。

图1 巷道模型及传感器布置示意图Fig.1 Tunnel model and sensor layout position

为了获取巷道内的参数，在火源上方及巷道轴向上2 m 高处每隔10 m 依次设立温度传感器和CO 气体体积分数传感器，为了进一步分析巷道顶部的温度分布，在火源上方的顶板沿巷道轴向每间隔10 m 也布设了温度传感器。巷道四周边壁由于采用了喷浆，选取软件材料库中混凝土的密度、热导率等参数。巷道边壁及初始环境温度为20 ℃，环境压力为1.0×105Pa。

1.2 工况参数设定

为研究正常通风情况下不同起火位置火灾参数的分布情况，设定巷道左边为风流速度入口，速度值为1 m/s，另一边设为开口边界，模拟时间为900 s。由于研究主要针对巷道顶板冒落后填充的有机物或锚网、支架后方的背板被引燃后的火灾，因此火灾功率不大，参考傅培舫等[12]在巷道内进行的全尺寸可燃物火灾参数试验，结合其他文献综合考虑[13-14]，将火灾总功率设置为1 MW，即火源单独位于顶板和底板时，火源功率为1 MW，火源同时位于顶底板情况时，顶板火源和底部火源分别为0.5 MW。

1.3 网格尺寸划分

在有限元分析软件中，网格尺寸越密集，模拟结果越精确。依据FDS 软件用户手册中火源特征直径及网格划分标准，计算得出模型网格尺寸应位于0.061～0.241 m 范围内。为了进一步检查网格独立性，对所采用的4 种网格方案模拟了巷道顶板处火源处的温度分布，不同网格尺寸下巷道纵向温度分布如图2。

图2 不同网格尺寸下巷道纵向温度分布Fig.2 Longitudinal temperature distribution in the tunnel under different grid sizes

由图2 可得：当网格尺寸采用0.1 m 和0.2 m混合方案时，和全部采用0.1 m 网格及0.08 m 的网格方案模拟结果趋于一致。为节省计算时间，采用混合网格方案，即在火源上下游附近20 m 内采用0.1 m 网格，其它区域采用0.2 m 网格，且该方案已在HU 等[15]的研究中得到验证。网格尺寸划分示意图如图3。

图3 网格尺寸划分示意图Fig.3 Schematic diagram of grid size division

2 结果分析

2.1 烟流分布

不同火源位置状态下烟流达到稳定状态后的分布情况如图4。

由图4 可知：火源位于底板和顶底板时，下方火源会产生火焰区、浮力区和羽流区，位于下部的烟流在上升撞击顶板后沿上下游运动；而火源位于顶板时，烟流产生后直接沿顶棚向上下游运动，3 种模拟工况在该风速下烟流均有逆退现象，且烟流在运行至一定距离后受巷道机械通风和摩擦阻力影响而停止逆退，转而向下风向运动直至出口处。

虽然3 种状态下火源功率总和一致，但3 种状态下的烟流分布情况和烟流逆退距离却表现不同。

在烟流分布上，尤其是火源下游位置，火源位于顶板时，烟流整体分布于巷道顶部，和底部的新鲜空气分界面清晰，而火源位于底板和同时位于顶底板时，在下游处烟流流动出现一定程度的紊流现象，烟流沉降明显，整个巷道几乎被烟流充斥。

不同工况下火源附近烟流速度矢量分布如图5。

图5 不同工况下火源附近烟流速度矢量分布图Fig.5 Velocity vector distribution of smoke flow near the fire source with different conditions

由图5 可知：当火源位于顶板时，火源附近烟流基本保持沿水平面层流运动，顶板附近烟流受巷道壁面黏性阻力影响，速度小于下部风流；当火源位于底板和同时位于顶底板时，在火焰下角出现旋涡等紊流现象，这可能进一步加剧了下游方向烟流的卷吸和沉降；此外，下部火源产生的烟流受浮力和纵向风流综合作用，向下游顶板方向偏移，但火源同时位于顶底板时，由于顶板火源产生的烟流与底板火源产生的烟流在火源上方靠近顶板处相互作用，底板火源的火焰偏移角度要略大于单一火源位于底板的工况。

在烟流逆退距离方面，当火源位于底板时烟流逆退的距离最短，而火源位于顶板时次之，火源同时位于顶板和底板时烟流逆退距离最大。

根据烟流逆退分析理论，烟流逆退是热浮力和惯性力作用的结果，即烟流前锋处静压差和巷道通风的动压相等时，烟流停止逆退[16]，如式（1）：

式中： Δp为静压差，Pa；hv为动压，Pa。

其中静压差可表示为：

式中： Δρ为烟流密度与空气密度差，kg/m3；g为重力加速度，m/s2；h为烟气厚度，m。

巷道纵向通风的动压可表示为：

式中： ρ为纵向风流空气密度，kg/m3；v为巷道速度，m/s。

将式（2）和式（3）代入式（1）可得：

由理想状态方程可得：

式中：T为纵向新鲜风流温度，K； ΔT为烟流与新鲜风流温度之差，K。

前人研究[17]表明，纵向烟流温度随火源距离呈指数分布，如式（6）：

式中：ΔTmax为巷道顶部烟流与新鲜风流最大温差，K；k1、k2为其他因子影响系数；H为巷道宽度，m；x为与火源距离，m。

将式（6）和式（5）代入式（4），进一步分析可得逆退距离L为：

由式（7）可知，烟流逆退受顶板烟流最大温差、风速、烟流厚度、巷道宽度等参数影响。在本文模拟工况下，风速、巷道宽度等参数均一致，因此逆退距离仅受烟流厚度h和巷道顶部烟流最大温差ΔTmax影响，烟流逆退距离函数关系图如图6。但由图4 可知3 种模拟工况下上游烟流的厚度差别不大，火烟逆流的动力更多决定于顶部烟流的最大温差，可在对巷道顶部烟流最大温度的分析可以进一步得到验证。

2.2 温度分布

不同火源位置状态下烟流稳定时巷道顶板附近的温度分布如图7。

图7 不同火源位置下巷道顶板处纵向温度分布Fig.7 Longitudinal temperature distribution at the roof with different fire source locations

由图7 可得：3 种状态下温度分布均表现出火源附近高，上下游温度随与火源距离增加而降低的规律；但火源位于底板时，最大温度出现在火源下游的底板处，这是因为火源产生的烟流受纵向风流影响，火烟上升过程中向下游偏移，当上升至顶板处时已位于火源上方顶板的下游。

对比3 种情况下的火源附近顶板处的最大温度，可得火源同时位于顶底板时略大于火源位于顶板，火源位于底板时最小，且明显低于其它2种工况。分析其原因：当火源位于底板时，由于距离顶板较远，且烟流上升过程中与巷道纵向风流对流交换影响，导致热量损失，而火源位于顶板时，火源距顶板近，燃烧产生的热量绝大多数用来加热顶部空气，受浮力作用卷吸下方的冷空气较少，因此顶板烟流最大温度明显没有火源位于顶板时高；而当顶底板均存在火源时，结合图5分析可知，底板火源产生的高温烟流撞击顶板时已在顶板火源的下风向处，下部火源产生的烟流上升中卷吸空气产生的烟流加上撞击顶板后向上游的射流进一步阻碍了顶板火源产生的烟流向下风向流动，使得更多的热量集中在顶板火源处，温度更高。

在火源上游方向，随着与火源距离的增加，顶底火灾工况下的温度衰减趋势要大于火灾位于顶板的工况，但随着离火源距离进一步增大，反而发生反转，表现出火源位于顶底板时的温度要大于火源位于顶板时，而火源位于底板的情况虽然随距离增大下降趋势较缓，但始终温度最低。在火源下游方向，离火源较近处，温度随距离增加同上游一样急剧衰减，在距离火源远端，3 种工况下的顶板处烟流温度分布则与上游完全相反，火灾位于底板时温度最大，但3 种工况下的温差并不明显，表明火源位于底板时，产生的高温烟流更多地流向下游方向，从而使下游方向顶板处的烟流温度更高。

分析煤矿巷道火灾时温度对人员逃生的影响，常考虑人员头顶高度处火烟的温度，鉴于巷道的底板铺设有轨道、巷道发生底鼓和劳保鞋的高度，绘制的不同火源位置状态下900 s 时巷道2 m 高处的温度分布如图8。

图8 不同火源位置下巷道2 m 高处纵向温度分布Fig.8 Longitudinal temperature distribution at 2 m height with different fire source locations

由图8 可得：与顶部温度分布规律一样，2 m高处的温度也呈现出火源附近高，随火源距离增加温度逐渐降低的规律。但3 种工况下，无论火源附近的最大温度值还是火源上下游的空气温度，火源位于顶板时温度值均要显著低于其他两者情况，且在火源两侧的变化不大，其中火源位于底板时，各处烟流温度值最大。这是因为顶板火灾烟流分布较为稳定且主要集中于上部，而底板火源产生的烟流上升过程中卷吸了更多的空气，烟气质量大，且新鲜风流及上部逆退的烟流再次流经火源后被进一步加热，使得下游底部空气的温度较高，与上部烟流温差进一步缩小，降低了烟流浮力，进一步使得烟流下沉，造成烟流厚度增加，从而影响了低处的空气温度。参考对火灾逃生临界温度60 ℃的规定[18]，可知火源位于顶板时整个巷道满足逃生条件，而火源位于底板时和火源位于顶底板时，在火源下游一定范围内，温度则不满足逃生条件，因此矿井救援需考虑火源发生的位置。

2.3 CO 体积分数分布

不同火源位置状态下2 m 高处CO 体积分数沿巷道轴向的分布情况如图9。

图9 不同火源位置下巷道2 m 高处纵向CO 体积分数分布Fig.9 Longitudinal CO concentration distribution at 2 m height with different fire source locations

由图9 可得：对于顶板火灾和顶底板火灾，2种工况下CO 峰值均出现在逆退烟流的前锋附近且两者峰值差别不大，这是由于顶板火灾产生的烟雾在热动力作用下直接沿顶板向上游蔓延，在烟流前锋处才向下卷吸新鲜空气；从烟流最大逆退处到火源位置，顶板火灾和顶底板火灾上游CO 体积分数明显高于火源位于底板时，随着空气对烟气的不断稀释，CO 体积分数迅速降低；而在火源下游，由于2 m 高处烟流和空气已充分混合，随着距离增加，各位置CO 体积分数略有下降，但差值不大；对比于底板火灾，CO 体积分数峰值位置则出现在火源附近的下游，且峰值明显低于前2种工况，这是由于底板火源产生的火烟受浮力作用向上运动过程中就会卷吸新鲜空气从而被稀释，且在纵向风流作用下直接向下游运动，因此峰值体积分数较低。由于底板火源产生的火羽流撞击顶板后方沿顶板向上下游方向蔓延，因此CO 体积分数随着与火源距离的增加逐渐降低。由于不同火源位置会造成CO 体积分数分布差异，因此，矿井应急救援工作同样要考虑火源位置。