从正常通风到异常爆炸:矿井瓦斯与煤尘物理过程的仿真技术解剖

摘要:矿井里最关键的物理过程,平时看起来只是通风、稀释和排放,一旦进入异常状态,就会迅速转化为瓦斯积聚、点火燃烧、冲击波传播、煤尘卷吸和混合爆炸的连续链条。真正要把这件事讲清楚,不能只靠经验口号,必须回到流体力学、燃烧学和数值仿真的共同语言。

矿井里最关键的物理过程,平时看起来只是通风、稀释和排放,一旦进入异常状态,就会迅速转化为瓦斯积聚、点火燃烧、冲击波传播、煤尘卷吸和混合爆炸的连续链条。真正要把这件事讲清楚,不能只靠经验口号,必须回到流体力学、燃烧学和数值仿真的共同语言。

从技术角度看,井下环境并不是一个“静止空间”,而是由主通风系统、局部通风、采空区渗流、巷道阻力、设备扰动和煤层释气共同构成的动态流动系统。正常情况下,通风的任务是把瓦斯浓度持续压在爆炸下限以下,并把局部释放出来的可燃气体及时带走。NIOSH 的资料给出的工程常识很明确:主通风系统的基本目标,就是把甲烷稀释在远低于爆炸区间的位置;在美国煤矿实践中,机载监测通常在甲烷浓度达到 1% 时触发停机,而爆炸区间大致位于 5% 到 15% 或 16% 之间。也就是说,正常通风不是“有风就行”,而是一个必须长期保持浓度裕度的控制问题。

真正的难点在于,矿井通风从来不是均匀场。长壁工作面、尾巷拐角、回风死区、采空区边界、密闭区附近,以及大气压快速变化时的异常渗流,都会让局部瓦斯浓度分布出现显著非均匀性。NIOSH 关于长壁通风的研究反复强调,尾巷角和采空区相邻区域,是甲烷积聚最敏感的位置之一。平时人们看到的是风量数字,真正决定风险的却是空间里的速度场、浓度场和渗流路径。

这就是为什么矿井通风研究必须依赖数值模拟。仅靠单点传感器,你只能知道某个位置此刻的浓度;而靠 CFD 或通风网络模型,你才有机会看到“气体是怎么从煤体释放出来、怎么沿着巷道和采空区迁移、怎么在局部形成高浓度带”的全过程。对正常工况而言,核心物理量是速度、压力、湍流强度、甲烷源项、巷道阻力和边界渗透性;对异常工况而言,问题会立刻升级为燃烧波、冲击波、颗粒卷吸和二次传播。

矿井瓦斯燃烧与煤尘传播过程的典型示意

先说异常链条的第一步:瓦斯积聚和点火。井下甲烷爆炸要满足三个基本条件,文献里通常归纳为:甲烷浓度进入爆炸区间、氧含量足够、并存在足够能量和持续时间的点火源。技术上最值得注意的不是这个三要素本身,而是它们如何在时空上重合。正常工况下,释气、扩散和通风稀释是同时发生的;一旦局部风流紊乱、回流增强、障碍物后形成滞止区,或者采空区释气在短时间内压过通风稀释能力,浓度就可能从“可接受波动”跃迁到“可燃区间”。这个阶段的关键不是宏观平均值,而是局部团簇。

点火之后,第二步不是立刻“大爆炸”,而通常是先形成甲烷-空气混合物的火焰传播和超压上升,也就是典型的 deflagration 过程。实验和数值研究都说明,火焰前沿的传播速度、压力上升速率、障碍物诱导的湍流增强,以及巷道几何变化,会决定这一步的强弱。对工程仿真来说,这一阶段的核心是可压缩流、燃烧放热和湍流耦合,而不只是“点着了没有”。

更关键的第三步,是甲烷初级爆炸如何把煤尘带进来。经典矿山爆炸研究早就指出,很多严重的煤尘爆炸不是凭空开始,而是由较小的甲烷爆炸触发。原因很物理:甲烷燃烧产生的超压和风浪会掠过巷道底板、帮壁和顶板附着层,把沉积煤尘重新卷吸到空气中。美国矿山局早期关于煤尘爆炸抑制的模型写得很直白:气压扰动把煤尘带入膨胀燃烧区,颗粒受热后释放挥发分和焦油,可燃产物与氧气反应继续放热,而放热又反过来维持更强的气动卷吸,于是爆炸从“气体火焰”转成“煤尘自维持传播”。这一步,是矿井异常过程从局部燃烧升级成长距离传播灾变的分水岭。

煤尘一旦参与,问题就不再是单一甲烷燃烧,而变成甲烷-煤尘混合爆炸。近年的管道和容器试验给出过很有价值的量化结果:在一定粒径和浓度范围内,煤尘与甲烷叠加后,最大爆炸压力和最大压力上升速率会明显高于煤尘单独爆炸;更细的颗粒由于更容易分散和挥发,危险性还会进一步提高。这里的物理本质是,甲烷提供更容易点燃的气相燃料,煤尘则提供后续更持久的热释放和空间延展能力。

从点火、冲击波到卷尘参与传播,异常过程本质上是多物理场耦合问题

因此,从“正常通风”到“异常爆炸”的完整物理过程,可以概括成一条连续链:煤层释气与通风稀释之间的平衡被打破,局部甲烷进入爆炸区间;点火源触发甲烷火焰传播,形成超压和高速流动;冲击波和风浪把沉积煤尘卷吸进入空气;煤尘受热、挥发、燃烧并参与火焰前沿传播;混合爆炸继续沿巷道传播,同时生成更多高温、高压和有毒气体。

这也是为什么仿真技术会分层展开,而不是只靠一种模型。第一类是通风网络和低维流动模型,用来算正常工况下风量分配、阻力变化、采空区与回风系统之间的连通关系,以及异常气体释放时的整体通路。近年的研究里,1D FDM 或 Flowmaster 一类网络模型也被用于模拟甲烷爆炸冲击波在复杂巷道网络中的传播。

第二类是三维 CFD 与专用爆炸求解器,例如 FLACS、GASFLOW-MPI、以及基于可压缩 Navier-Stokes 方程的燃烧模型。这一类模型擅长处理局部高浓度区、障碍物诱导湍流、火焰加速、压力波反射与干涉等问题。已有研究表明,针对实验矿井网络的甲烷爆炸传播,GASFLOW-MPI 这类工具与全尺度实验结果可以做到较小误差量级。

第三类是颗粒-气体耦合和冲击动力学模型,用来处理煤尘卷吸、颗粒尺寸效应、爆炸超压对巷道结构和人体的冲击效应。这一层模型往往会引入颗粒输运、热解挥发和显式动力学分析。

真正成熟的技术路线,应该是把这三层联合起来:平时用通风网络和甲烷迁移模型看正常工况的浓度安全裕度;对关键区域用三维 CFD 做局部流场和积聚风险精细化分析;对异常工况再用爆炸传播和颗粒耦合模型研究火焰、冲击波与煤尘参与机制。

从学术角度看,矿井瓦斯与煤尘问题的本质,是一个“多尺度、多相、多物理场”的耦合系统。正常阶段它是通风与释气的流动问题,异常早期它是可燃混合物形成与点火的问题,传播阶段它又变成可压缩燃烧、冲击波传播和颗粒卷吸的问题。看清这条链条,本身就是理解和控制风险的前提。

参考来源

  1. NIOSH, Mining Feature: Coal Mine Explosion Prevention. https://archive.cdc.gov/www_cdc_gov/niosh/mining/Features/coalmineexplosion.html
  2. NIOSH, Mining Topic: Ventilation Overview. https://archive.cdc.gov/www_cdc_gov/niosh/mining/topics/ventilation.html
  3. Chi, D. N. & Perlee, H. E., Numerical Simulation of Flame-Induced Aerodynamics in a Coal Mine Passageway, Semi-Empirical Model, U.S. Bureau of Mines, RI 8016, 1975. https://stacks.cdc.gov/view/cdc/237402/cdc_237402_DS1.pdf
  4. Cashdollar et al., Model of Coal Dust Explosion Suppression by Rock Dust Entrainment. https://stacks.cdc.gov/view/cdc/10608/cdc_10608_DS1.pdf
  5. Maier, P., Hartlieb, P., Brune, J. F., Laboratory Scaled Coal Dust Explosions and Physical Test Results for CFD Explosion Models. https://link.springer.com/article/10.1007/s00501-020-00985-0
  6. Tan, B. et al., Comparative study of the explosion pressure characteristics of micro- and nano-sized coal dust and methane–coal dust mixtures in a pipe. https://link.springer.com/article/10.1007/s40789-019-00289-w
  7. Gao et al., Numerical Investigation of Network-Based Shock Wave Propagation of Designated Methane Explosion Source in Subsurface Mine Ventilation System Using 1D FDM Code. https://www.mdpi.com/2071-1050/16/22/9935
  8. Zhang et al., Damage Effect and Injury Range of Shock Waves in Mine Methane Explosion. https://www.mdpi.com/2674-0389/3/4/33
分享到