采空区氧化自燃带宽度与工作面推进速度关系的数值模拟.pdfVIP

2011年lO月 矿 业安全 与环 保 第38卷第5期 采空区氧化 自燃带宽度与工作面 推进速度关系的数值模拟 李尚国 ，陆 伟 ，梁林秀。 (1．瓦斯灾害监控与应急技术 国家重点实验室，重庆400037；2．中煤科工集团重庆研究院，重庆400037； 3．神华乌海能源有限责任公司路天煤矿，内蒙古乌海016000) 摘 要：基于采空区渗流方程和多孔介质动量方程 ，应用FLUENT软件模拟了采空区漏风风速变化 情况，结合采空区自燃 “三带”划分标准，得出采空区氧化 自燃带宽度与工作面推进速度的关系：随着工 作面推进速度提高，进风侧、回风侧、工作面中部氧化 自燃带与工作面距离都将增大；工作面中部氧化 自 燃带与工作面距离增幅较大，进风侧、回风侧氧化 自燃带与工作面距离增幅较小。 关键词：采空区；漏风风速；氧化 自燃带；推进速度 中图分类号：TD75 2．2 文献标志码 ：B 文章编号：1008—4495(2011)05—0010—03 采空区浮煤 自燃是矿井实际生产过程 中常常发 — — 多孑L介质粒子形状系数； 生的灾害现象。采空区靠近工作面附近区域，由于 — — 压力坡度 ； 漏风风速较大热量难 以积聚，采空区深部漏风风速 — — 采动裂隙椭抛带的风速，m／s； 较小且风流经过长距离的输送，O，大量损耗，浮煤 — — 运动黏性系数，m ／s； 常常难以自燃 ¨I4』。因此 ，在采空区内形成了散热 D ——平均调和粒径 ； 带、氧化 自燃带和窒息带。 “三带”的划分是从宏观 u——渗流速度 ，m／s。 上分析采空区 自然发火的主要依据 ，而 以漏风风速 I．1．2 采空区多孔介质动量方程 划分 “三带”在工作面 回采过程 中很难实际测量获 多孔介质模型是在所谓的多孔介质区域中加上 取。为克服上述困难，笔者用计算机数值模拟方法， 由经验公式确定的流动阻力 ，实质是在流体流动控 借助于图形显示技术 ，给出更详细直观的氧化 自燃 制方程的动量方程中加上额外的动量源(汇)项。源 带形状的数值解，以及当工作面推进速度发生变化 项由两部分组成，一部分是黏性损失项 (Darcy)，另 引起氧化 自燃带变化的结果。 一 个是惯性损失项 ： 2 2 1 l 模型及参数确定 s=∑Ddzvj+∑C如 (2) J 1 J= 1 1．1 采空区多孔介质模型 式中 Ls—— 向(，Y)动量源项 ； 1．1．1 采空区渗流控制方程 D，c—— 规定的矩阵； 将综放采空区视为由煤岩混合体组成的多孔介 — — 流体的动力黏度 ，Pa·S； 质空间，松散煤体中的漏风流场十分复杂，将采空区 _ 向流体的渗透速度 ，rn／s。 内风流流动看作不可压缩气体在二维空间的非线性 在多孔介质单元 中，动量损失对于压力梯度有 渗流，包括紊流、层流、过渡流，所服从的基本方程为