防灭火注氮钻孔周围氮气运移规律数值模拟研究

2015-05-08吴伯庆赵先科李海涛

江西煤炭科技 2015年1期

吴伯庆，赵先科，李海涛

（1.潞安新疆煤化工（集团）有限公司，新疆哈密839003；2.西安科技大学能源学院，陕西西安710054）

近年来我国对能源的需求已经进入高峰期，煤炭开采量日益增多，采空区面积不断扩大，留下大量浮煤，采空区自然发火事故频发，形势依然严峻。目前，常见的预防采空区自然发火的技术措施有注水、阻化剂、堵漏风等，但在地质条件复杂、煤层状况不稳定的地区效果并不理想。资料表明〔1〕，注氮防灭火技术能够有效降低采空区内的氧气浓度，对于地势错综复杂的矿井具有不可替代的优势。

大量学者对采空区注氮技术进行了研究。徐精彩、文虎〔2〕等根据传热学和化学动力学原理，提出了采空区下限氧浓度和氧化带宽度的计算方法，为采空区注氮后自燃火灾的预测及预防提供了理论依据；李庆源、张峰等〔3〕通过埋管测试方法，研究采空区注氮前后标志性气体在采空区内的分布状态，为优化注氮参数，提高注氮效果提供了现实依据；何宗礼〔4〕利用FLUENT软件对U型通风采空区注氮防灭火技术参数进行了研究，为日后研究采空区最佳注氮量提供了借鉴。纪忠、梁栋〔5〕等根据相似理论，采用实验模拟的方法研究了氮气在采空区三维空间中的运移规律，对于采空区注氮参数优化具有现实意义。何勇〔6〕通过分析采空区注氮技术在淮南矿业集团新庄孜矿的运用过程，总结了注氮防灭火技术的具体实施方案，为今后采空区注氮技术的运用提供了参考。

对于防灭火注氮设备而言，氮气是在注氮压力的作用下冲出注氮钻孔，具有一定的运移规律。本文从压缩气体膨胀的原理出发，采用MATLAB软件对氮气在钻孔周围的运移规律进行了模拟，对注氮钻孔周围氮气的运移规律进行了研究。研究方法和结论为今后氮气在多孔介质内扩散规律的研究以及防灭火注氮设备参数优化提供了理论基础。

1 注氮钻孔周围氮气流动数学模型

氮气在采空区的运移可看做是在多孔介质流场的均匀流动〔1〕，满足连续性方程、能量守恒定律和达西定律。

（1）连续性方程

注氮钻孔周围氮气流动满足质量守恒定律：

式中：ux、uy、uz分别为氮气流动x、y、z三个方向的速度分量，t为流动时间，ρ为氮气的密度。

（2）能量守恒

氮气从注氮钻孔喷出，钻孔周围空间足够大，钻孔内流体阻力损失可以忽略，假设钻孔内气体密度恒定，根据流体力学可知，注氮钻孔周围氮气流动过程满足可压缩空气的绝热变化能量方程：

式中：P、P2分别为注氮压力和大气压力；γ1、γ2分别为注氮钻孔内压缩气体密度和空气密度；k为空气的绝热常数；v1、v2分别为注氮速度和钻孔的出口速度；g为重力加速度。

（3）达西定律

由多孔介质中气体流动理论知，注氮钻孔周围的采空区氮气流动符合达西定律，其径向不稳定流动时的偏微分方程为：

式中：P为采空区氮气压力的平方；λ为采空区氮气渗透率；α为采空区氮气含量系数；r为注氮钻孔半径。

2 物理模型与边界条件

2.1 物理模型

选用直径0.1m的注氮钻孔，其二维物理模型如图1所示，利用PDE工具箱画出计算区域，采用布尔运算方法划分与加密网格，划分结果见图2。

图1 二维物理模型

图2 模型网格划分

2.2 边界条件

在没有外界影响的条件下，采空区内各点氮气压力保持起始值P0，注氮时，钻孔周围径向上的氮气压力开始增加，影响范围不断扩大。氮气流场边界条件分为有限和无限两种，无限流场即假设钻孔无法影响到钻孔边界，因此其边界压力始终保持原有压力。当钻孔影响到期边界时，边界条件为：

式中：P0为采空区氮气起始压力；P1为注氮时钻孔内压力。

初始压力：P＝P0＝P02（t＝0）

3 数值计算及模拟结果分析

3.1 数值计算

为求解方程，对达西定律采用拉普拉斯变换，得：

将（5）式化为标准化的贝塞尔方程，并解之得：

由（4）式得：

将（7）式带入（6）式，得：

方程（8）的求解异常复杂，常规方法难以对其进行求解，目前利用计算机求解。本文采用PDETOOL进行求解。

3.2 模拟结果分析

从图3（a）可以看出，氮气的运移图形呈喇叭形，越往下氮气的扩张面积越大，这是由于氮气从注氮管道喷出以后，在注氮口附近形成局部补给源，在注氮压力和采空区负压作用下氮气沿着注氮方向发生运移，与此同时，由于注氮压力纵向分力大于采空区的环境压力，氮气不断地沿着垂直于注氮管道方向发生横向蔓延，并形成了以钻孔为中心的近似圆形的影响半径，占据越来越多的区域。在氮气流经的区域，原有采空区内气体不断被驱替，区域内各气体组分情况发生变化，新形成的混合气体在注氮推力作用下向前流动，改变了采空区原有的风流运动，形成新的渗流场。从图3（b）可以清晰地看出同一截面上离钻孔不同距离的氮气压力分布情况。由于氮气气体的粘性力和管道表面作用力，靠近管道中心一侧的气体压强大，而管道壁面一侧的压强小，因而同一注氮管道内同一横断面上不同点处氮气的压强不同，呈向外依次减小趋势。矿井采空区内充满冒落的矸石和遗煤，可看做是多孔介质流场，氮气在采空区的运移近似均匀流动。