基于fluent综采工作面采空区内瓦斯浓度场研究.docVIP

基于fluent综采工作面采空区内瓦斯浓度场研究 　　摘要： 本文通过数值模拟采空区内瓦斯浓度场，从理论上确定了瓦斯浓度的分布情况。进而在对综采工作面瓦斯治理提供了科学的理论依据。 　　Abstract： Through the numerical simulation of methane concentration field in mined-out area， this paper determined the distribution of methane concentration in theoretical aspect. Then the scientific theoretical basis was provided in the methane control of fully mechanized working face. 　　关键词： 瓦斯；浓度场；数值模拟 　　Key words： methane；concentration field；numerical simulation 　　中图分类号：TD712 文献标识码：A 文章编号：1006-4311（2014）21-0163-02 　　1 182403工作面概况 　　井上182403工作面位于水鱼岗村西北方向，地势北高南低，地表为丘陵地形，地表大部分已被当地居民开垦为农田，跃峰渠流经该处。井下位于四采区左翼，本采区182408工作面和182402工作面已结束，附近无其它采掘情况。 　　403工作面煤层厚度为3.01m～3.19m，平均煤厚3.1m，局部有夹矸，结构简单。老顶为细砂岩，厚6.10m；直接顶为砂页岩，厚3.24m；直接底为砂页岩，厚5.50m；老底为中砂岩，厚4.75m。 　　2 采空区瓦斯浓度场基础理论 　　为了进行插管抽采设计，需要确定抽采泵的流量，因为它决定着设备选型、管径等参数的计算等。根据煤矿瓦斯抽采工程设计规范，标准状态下抽采泵的流量与最大的设计瓦斯抽采量以及抽采泵入口处预计的瓦斯浓度等参数有关，目前主要采用类比法根据经验确定这些参数，急需理论上的指导。 　　在此次研究过程中采用解析法研究U型通风的漏风场和由插管抽采产生的流场，用数值模拟说明了上隅角封堵以及插管抽采在治理上隅角瓦斯浓度超限方面的作用，给出了详细的数值模拟方法以及有关参数的取法，插管抽采采空区瓦斯是治理U型通风上隅角瓦斯浓度超限的根本措施；数值模拟工作在插管抽采设计方面具有一定作用，所用方法可供采空区瓦斯抽采设计参考。 　　2.1 U型通风的漏风场 为了研究采空区内的复杂流场（包括氧浓度场、瓦斯浓度场、温度场等），首先研究采空区内的漏风场。 　　图1是U型通风采空区漏风场二维示意图。记工作面长度为L，采空区沿深度方向的瞬时长度为H，下隅角处的气体静压为P1，上隅角处为P2。根据无粘不可压缩无旋流理论，流场压力p满足拉普拉斯方程， 　　■+■=0（1） 　　边界条件为，■■ y=0=■■ y=L =■■ x=H =0，（x0）（2） 　　p■ x=0=P1-■y，（0yL）（3） 　　定解问题（1）～（3）是椭圆型线性偏微分方程，可采用傅里叶变换法求解（数学软件MATHEMATICA 8暂不能求解该类方程）。压力场为， 　　p（x，y）=■+■■■Exp 　　[-（2n-1）？仔x/L]（4） 　　我们主要关心x方向的流动。 　　根据达西定律，vx=-■■，由此得， 　　vx=■■■Exp[-（2n-1）？仔x/L] 　　（5） 　　在x=0的边界上，vx=■■■（6） 　　式中，r为采空区阻力系数（见后）。 　　式（6）在y=0和y=L处发散，是由无粘流假设造成的，实际上，由于粘性作用，在y=0和y=L处vx=0；我们关心的是0yL的情况，此时式（6）收敛。 　　定义？浊=y/L，Kx=vxrL/（P1-P2），则在x=0的边界上， 　　Kx=■■■（7） 　　Kx～？浊关系曲线如图2所示。由图2可见，在？浊=0和？浊=1处，Kx=±∞（级数发散）；在？浊=0.5处，Kx=0（沿回采工作面倾向采空区一半进风、一半出风）。 　　关于r，根据达西定律的常见形式vx=-■可知，r=■，这里，？滋为空气的动力粘性系数，主要与温度有关，标准状态下约为1.8×10-5Pa?s；K为采空区的渗透率，单位为m2，越往采空区内部K值越小，理论计算可采用乌沙阔夫公式，有关参数可通过试验获得。文献给出的某矿试验结果平均值为r=11.8Pa?s/m2（单一开采）和39.2Pa?s/m2（分层开采），换算为渗透率分别为K=1.525×10-6m2和4.592×10-7m2，可供参考。 　　关于P1、P2，可通过井下实测