y型通风技术在煤矿中的应用.docxVIP

y型通风技术在煤矿中的应用 大部分矿工事故发生在采矿巷道，而采空区的屋顶数量约为30%60%。典型的U型通风方式,由于受其漏风流场影响,极易造成工作面上隅角瓦斯超限。尽管可以对U型通风方式进行改进,但效果并不明显。对于瓦斯涌出量大的矿井,采用两进一回Y型通风方式能有效解决上隅角瓦斯积聚问题,因此弄清两进一回Y型通风采空区流场流动对治理上隅角、回风巷瓦斯超限及采空区遗煤自燃和瓦斯爆炸有着非常重要的理论和现实意义。 1 采空区流场数学模型 1.1 采空区流体流动规律 采空区复杂的空隙结构导致其气体流动也很不稳定。为了研究方便,弄清其与主要问题之间的关系,故作如下简化: 1) 将采空区多孔介质视为各向同性; 2) 采空区流体流动服从达西渗流定律; 3) 将采空区气体视为不可压缩气体,忽略有流体粘性力做功所引起的耗散热,即将流体流动近似为稳定、等温的流动过程; 4) 将组分运输中,气体设置为甲烷、空气的混合气体,忽略水蒸气、二氧化氮等小组分气体。 1.2 组分传输方程 根据质量守恒定律和渗流定律理论,对于Y型通风采空区流风流的稳态流动,由以上假设条件,可得控制方程如下(1)～(6): 连续性方程: ?u?x+?v?y=0?u?x+?v?y=0(1) 式中:u、v分别为速度矢量在x、y方向的分量。 动量(Navier-Stokes)方程: ?(nρuu)?x+?(nρuv)?y=??x(nμeff?u?x)+??y(nμeff?u?y)?n?p?x?μeu(2)?(nρuu)?x+?(nρuv)?y=??x(nμeff?u?x)+??y(nμeff?u?y)-n?p?x-μeu(2) ?(nρvu)?x+?(nρvv)?y=??x(nμeff?v?x)+??y(nμeff?v?y)?n?p?y?μev(3)?(nρvu)?x+?(nρvv)?y=??x(nμeff?v?x)+??y(nμeff?v?y)-n?p?y-μev(3) 式中:u、v分别为x、y方向的速度;ρ为采空区气流密度;μeff=μ+μt;μ为动力粘滞系数;μt为湍流动力粘滞系数;n为多孔介质孔隙率。 组分传输方程: ?(ρcsu)?x+?(ρcsv)?y=??x(Ds?ρcs?x)+??y(Ds?vcs?y)+Ss(4)?(ρcsu)?x+?(ρcsv)?y=??x(Ds?ρcs?x)+??y(Ds?vcs?y)+Ss(4) 式中:cs为组分s的体积浓度;Ss为组分生产率;Ds为该组分的扩散系数。 标准k-ε方程: ?(ρku)?x+?(ρkv)?y=??x(μ+μtσk)?k?x+??y(μ+μtσk)?k?y+Gk?ρε(5)?(ρku)?x+?(ρkv)?y=??x(μ+μtσk)?k?x+??y(μ+μtσk)?k?y+Gk-ρε(5) ?(ρεu)?x+?(ρεv)?y=??x(μ+μtσε)?ε?x+??y(μ+μtσε)?ε?y+εk(C1εGk?C2ερε)(6)?(ρεu)?x+?(ρεv)?y=??x(μ+μtσε)?ε?x+??y(μ+μtσε)?ε?y+εk(C1εGk-C2ερε)(6) 式中:μt=ρCμk2/ε;C1ε=1.44;C2ε=1.92;Cμ=0.09;σk=1.0;σε=1.3;u、v分别为x、y方向的速度分量,m/s;k为湍动能,m2/s2;ε为湍动能耗散率,m2/s3;μ为动力粘滞系数,Pa·s;μt为湍流动力粘滞系数,Pa·s;Gk为由于平均速度梯度引起的湍动能k的产生项。 2 y型开口的物理模型和网格的划分 2.1 gampit建立 根据两进一回Y型通风回采工作面采空区工作面采空区的实际情况,利用Fluent前处理软件Gambit建立如图1所示的两进一回Y型通风回采工作面采空区物理模型。采空区长160 m、宽100 m,工作面长100 m、工作面宽5 m,机巷和风巷沿空留巷宽度均为3 m、机巷和风巷长度均为15 m,沿空留巷长度160 m,沿空留巷宽度3 m。 2.2 模型划分网格描述 使用建模工具Gambit在笛卡尔坐标系下建立模型,并在边界层加密的基础上,对模型进行网格划分,本次模型采用结构化网格划分技术,工作面网格间距0.2 m,采空区网格间距为0.5 m。采用控制容积法对控制方程进行离散,离散过程中的对流项和扩散项分别采用Quick格式,Pressure-Velocity Coupling压力速度耦合方式采用SIMPLE算法。 2.3 工作面与采空区的边界 入口边界:两条进风巷设置为速度入口(Velocity-inlet),机巷入口平均风速为2.25 m/s,风巷入口平均风速为0.75 m/s,均匀风流垂直于进风巷入口进入,进风巷瓦斯的浓度为0;出口边界:设置为自由出流(Out-flow);所有壁面为无滑移边界