用matlab模拟采空区瓦斯涌出量,采煤工作面瓦斯流动模型及COMSOL数值解算

工作面瓦斯涌出量大,易在采空区中发生积聚,特别是上隅角瓦斯超限事故时有发生。因此,掌握工作面的瓦斯压力分布和瓦斯涌出规律,准确预测工作面的瓦斯涌出量,分析瓦斯分布规律,对有效地治理瓦斯具有重要的意义;对工作面瓦斯流动机理、运移规律以及瓦斯压力分布的基础研究显得很有必要。林柏泉从研究煤层瓦斯流动机理出发,提出了“煤层瓦斯运移数学模型”[1],这种方法为研究煤层瓦斯运移、瓦斯压力分布和瓦斯涌出规律提供了理论基础。随着计算机技术的发展,利用CFD(计算流体力学)模型对采空区瓦斯流动规律进行研究越来越受到欢迎。然而采煤工作面瓦斯流动是一个复杂的过程,基本的CFD模型很难准确地对采煤工作面进行描述。本文针对采煤工作面和采空区流动特性的差异以及瓦斯的扩散运移特性创建瓦斯流动模型,利用COMSOLMultiphysics软件耦合求解非线性方程的高效性,对模型进行数值解算[2-4]。1工作面气体流动模型1.1工作面流动方程Navier-Stoke方程能很好地描述管道内流体的流动规律,无论是微风还是湍流,都可以通过Navier-Stoke方程求解,煤矿巷道和工作面通风可以视为管道流动,本文采用Navier-Stoke方程作为工作面空气流动控制方程:-(uns+(uns)T)+unsuns+pns=0(1)uns=0(2)式中,为黏性系数,kg/(ms);u为速度矢量,m/s;为液体密度,kg/m3;p为压力;Navier-Stokes方程的因变量分别为速度(u)和压力(p),下标用“ns”表示。1.2采空区流动方程采空区内充满冒落的矸石,而矸石间又存在大量的空隙和裂缝,气体在采空区内的流动规律是介于Darcy渗流和Navier-Stokes管道流动之间的过渡流动,Brinkman方程描述的孔隙介质中流体快速流动正好适用于采空区,因此常用来描述采空区内的空气流动规律。采空区渗透性系数的选取成为一个关键。李宗翔等曾描述过采空区渗透性系数与岩石垮落碎胀系数的关系[5-6],Brinkman方程的因变量分别为速度矢量和压力。-(ubr+(ubr)T)-ubr+pbr=0(3)式中,为孔隙率;k为渗透系数。1.3空气流动模型通过研究对比发现,Navier-Stokes方程能很好地适用于巷道内空气流动,Brinkman方程针对破碎带兼顾流体压力梯度和运动作用的特点,比较适合破碎带渗流对于描述垮落带的空气流动[2-6]。利用Navier-Stokes方程和Brinkman方程构建出采煤工作面空气流动模型(图1)。将方程(1)和(3)联立求解,建立如图1所示的空气流动模型。求解模型可以得出空气流动平衡状态下压力场分布和速度场分布,即得出工作面和采空区速度和压力分布情况。在N-S模型中,与Brinkman流动区域的速度和压力在界面上是一致的:uns=ubr,pns=pbr。这表明在图1空气流动模型Fig.1Atmosphereflowmodelvn,vb,pn,pb分别为N-S与Brinkman流动区域的速度和压力漏风面上,工作面内的速度和压力与采空区是一致的(即漏风面为连续面)。N-S模型的边界条件约束如下:在巷道壁上速度为0,对应无滑移边界条件,进风巷风量也是已知的,巷道的中心线应用对称边界条件,回风巷道的出口压力已知。2瓦斯扩散运移模型瓦斯在采场的运移遵守流体动力弥散定律,而煤屑内瓦斯运动基本符合线性扩散定律菲克定律[5,7],仅用达西定律并不能完全反映瓦斯扩散和运移规律,溶质在流场中的运移规律可用平移-扩散方程来描述:sct+(-sDLc+uc)=Sc(4)其中,s为流体体