工作面采空区瓦斯涌出量与供风能力的关系.docxVIP

工作面采空区瓦斯涌出量与供风能力的关系 1 工作面通风安全分析 据统计，目前大型和大型开采工作场所的实际供气量超过10003m，这主要是由于传统方法中对各行各业的冷却和开采、重量、灰尘和风速的冷却和影响。尤其对高瓦斯及突出矿井的综放工作面,单以传统的工作面风量分配方法已难于适应新工作法的要求,客观上还必须采用新的分析方法。综放开采自身固有的最突出安全问题就是采空区瓦斯大量涌出和采空区大量遗煤的自然发火,这涉及到复杂的场流安全理论,必须用现代通风理论和方法加以解决。因此,寻找适合综合放面的通风参数是一个新的通风安全重要课题。笔者将用数值模拟方法,对风量调整范围及安全评价从理论进行深入分析和探讨。 2 氧浓度耗氧速度内镜c 多数研究者都将采空区看成包括紊流、层流和过渡流的二维多孔介质非线性渗流场,服从Bachmat方程。为便于求解,在工作面和漏风点附近用大渗透性系数的高速达西渗流耦合计算来拟合紊流和过渡流,并实现流态的逐渐过渡,而采空区内部属于线性达西渗流。如图1所示(qL,q′L为工作面向采空区漏入、漏出风量),对于简化的采空区非均质二维流场,漏风流运动定解数学模型为 div(kHgradp)=0, (1) k=b(Κp-1)3ΚΡ?Η=ΚΡΜΚ0Ρ-1k=b(Kp?1)3KP?H=KPMK0P?1 KP=K′P+(K0Ρ0P-K′P)e-ax,(在采空区上) p=r1·Q2·(l-y)(在工作面边界上) -kgradp=0(在采空区固壁边界上) 流函数方程 div(φgradψ)=0(2)div(φgradψ)=0(2) 其中, φ=max{k(x,y)}k(x,y)?ψ=0?(在采空区固壁边界上)φ=max{k(x,y)}k(x,y)?ψ=0?(在采空区固壁边界上) 同一流场中,假设速度场是定常的,且是位势流动,瓦斯浓度变化过程为 cg=cg0(在漏入新风边界上) 氧浓度消耗过程的定解数学模型为 ndc(Ο2)dt=ndiv(D(Ο2)?gradc(Ο2))-div(c(Ο2)V)+W(Ο2)W(Ο2)=-(γ0?Η1c0+n?Η?WCΗ4n?Η+WCΗ4)c(Ο2)Η(4)c(Ο2)=c0(在漏入新风边界上) 以上所有公式中, p——风压,Pa; k——采空区冒落体的渗透性系数,m2/(Pa·s); b——与k同单位的待定系数; M——采放高度,m; H1——开采高度,m; H——采空区流场高度,H=H(x,y),m; n——冒落空隙度,n=n(x,y); KP——采空区冒落碎胀系数; K0Ρ——初始冒落碎胀系数; K′P——压实碎胀系数; a——碎胀系数的衰减率,由矿压观测确定; r1——工作面单位长度的风阻,Ns2/m9; l——工作面长度,m; Q——工作面风量,m3/min; ψ——流函数; φ——非均匀流场修正函数; cg——瓦斯浓度,mol/m3; cg0——工作面入风流瓦斯浓度,mol/m3; WCH4——是瓦斯涌出源项(暂不考虑放煤瓦斯团高峰影响),mol/(m3·h); W0——采空区瓦斯涌出强度,mol/(m2·h); W′0——衰减初始涌出强度,mol/(m2·h); d——W0的衰减率; c(O2)——氧浓度,mol/m3; c0——标准空气中氧气浓度,mol/m3; γ0——松散煤矸混合体在标准空气中的耗氧速度,mol/(m3·h); W(O2)——综合考虑采空区内煤矸氧化和瓦斯涌出稀释作用的耗氧汇项,mol/(m2·h); Dg,D(O2)——分别是瓦斯和氧气的弥散系数张量,m2/h; t——变化时间,h; V——由式(1)解出的渗流速度场函数。 上述模型采用迎风有限元解法,编制了求解计算机程序,首先求速度场,然后求气体浓度分布。计算时进行了简化:采空区瓦斯涌出量和煤氧化产生的气体远小于采空区漏风量,不改变采空区流态;煤尚未开始氧化自热,不考虑热风压的影响。渗透性系数识别系数b通过观测工作面的风流损失量与数值模型进行拟合得到,煤耗氧速度γ0一般由低温吸氧实验测定,这里按氧浓度观测资料反求获得。 以东滩煤矿4302综放工作面为算例,工作面长度l=200 m,r1=0.00 156 Ns2/m9,M=5.8～8.8 m,H1=2.8 m,煤的最短自然发火期为37天。取a=0.0268,Kp=1.15～1.45,此时n=0.13～0.31,反求得b=0.171,渗透性系数k=0.000 5～0.010 6 m2/(Pa·s),气体弥期散度5 m(速度纵向)、1.5 m(速度横向),而气体分子扩散系数0.075 6 m2/h。取采空区煤岩混合体的耗氧速度,其γ0=0.795 mol/(m2·h)。工作面的氧浓度c0=9.375 mol/m3。选取该矿综放面顶煤壁瓦斯涌出强度0.08 mol/