基于fluent的采空区瓦斯分布规律研究

摘要

近年来随着矿井开采的不断延深，瓦斯梯度也不断增大，进而造成工作面瓦斯涌出量很大，而采空区瓦斯涌出在工作面瓦斯涌出中占有很大的比例，因此只有了解采空区瓦斯的赋存、运移、涌出规律，才能合理地调整采场通风系统，采取行之有效的瓦斯治理措施防止瓦斯集聚。

因此本文以采空区瓦斯、空气混合气体中的瓦斯（CH4）为主要研究对象，通过理论分析的方法较系统分析了综采工作面及其采空区瓦斯运移的基本特征，以及影响采空区瓦斯影响涌出的主要因素，并以此为理论基础建立了采空区瓦斯运移的数值模型。紧接着以顾桥矿1123(1)工作面为原形在FLUENT软件中搭建模型，通过数值模拟对比分析了仅有漏风和布置有高抽巷情况下的采空区内的瓦斯分布规律，得出高抽巷的抽排作用使得综放面采空区内的高浓度瓦斯向抽采,口运移，同时扩大了工作面漏风向采空区内深入的距离，减少了采空区内瓦斯向工作面的涌出。最后，根据数值模拟结果结合现场测量数据对1123（1）工作面的瓦斯抽放效果进行分析评价。

通过本文的分析验证数值模拟与现场实测结果基本吻合，由此可推知采用FLUENT进行数值模拟对采空区瓦斯抽采措施的制定具有一定的指导意义。

关键词：采空区；瓦斯渗流；FLUENT软件；数值模拟；瓦斯抽采。

Abstract

Recent years,with the increase of the depth of coal seams mined, methane gradient increased in a huge speed, causing a large amount of working face gas emission .And a large proportion of working face gas comes from gas emission of gob. So only if we get a good knowledge of the laws of gas storage, gas movement, gas emission, then we can adjust the ventilation system reasonably and take effective methane control measures to prevent the gas gathering.

This dissertation analysis the basic features of fully Mechanized mining face, goaf gas migration rules and the main factors which has the effects on goaf gas emission through theoretical analysis, taking the CH4 constituents in goaf gas-air mixed gas as the research object. And the , on the basis of analysis ,the third chapter 2 established the equations of the gas seepage in gob. The fluid mechanics calculation software FLUENT is used to establish the model of 1123(1) working face of the Guqiao colliery and simulate two different 3D models of gas distribution in gob of fully mechanized mining face between the model with air leakage only and the model with air leakage and high-located gas drainage roadway. The conclusion is high-located gas drainage roadway makes the high-concentration gas in the gob deliver to the drainage port, the distance of the air leakage go deep in the gob is increased too, so the result is the gas emission from the gob to the workface reduced. Finally the chapter 4 makes an evaluation about the effects of the methane control according to the simulated results and the field metrical data.

The simulated results are identical with the field metrical data, which is proved by this dissertation. So using FLUENT software to simulate the gas seepage in gob has the fixed guidance meaning to methane control measures.

Key words: goaf, methane seepage flow, FLUENT software, numerical simulation, methane extraction.

目录

第一章绪论 (1)

1.1 研究意义 (1)

1.2 国内外的研究现状 (3)

1.2.1 采空区混合气体流动场研究现状 (3)

1.2.2 瓦斯渗流理论的研究历程 (3)

1.3 研究思路 (5)

1.3.1 研究内容 (5)

1.3.2 研究方法 (6)

1.3.3 技术路线 (6)

第二章综采工作面瓦斯运移理论研究与数值模型 (8)

2.1 综采工作面基本特征 (8)

2.2 综采面采空区瓦斯来源及涌出影响因素 (9)

2.2.1 综采面采空区瓦斯来源 (9)

2.2.2 瓦斯涌出量计算 (10)

2.2.3 影响采空区瓦斯涌出的基本因素 (12)

2.3 采空区瓦斯运移特征 (14)

2.3.1 多孔介质的特性 (14)

2.3.2 采空区瓦斯的升浮 (16)

2.3.3 采空区瓦斯的扩散 (17)

2.3.4 采空区瓦斯的聚集 (17)

2.4综放面采空区瓦斯运移基本方程 (18)

2.4.1 采空区渗流基本方程 (18)

2.4.2 扩散运动方程 (21)

2.5 控制微分方程的建立 (22)

2.5.1 基本假设 (22)

2.5.2采空区气体流动方程 (22)

2.5.3 采空区瓦斯的动力弥散方程 (23)

第三章 FLUENT 软件介绍以及综采面采空区瓦斯运移规律数值模拟 (25)

3.1 FLUENT软件介绍 (25)

3.1.1 FLUENT软件组成 (25)

3.1.2 FLUENT求解 (26)

3.2 物理模型的构建与边界条件设定 (27)

3.2.1 物理模型的规定和假设 (27)

3.2.2 计算模型的建立与网格生成 (27)

3.2.3采空区瓦斯运移数值模拟参数及边界条件的设定 (28)

3.3 数值模拟结果的分析 (31)

3.3.1 只考虑漏风时采空区瓦斯浓度的分布模拟与结果分析 (31)

3.3.2设置底抽巷时的瓦斯浓度的分布模拟与结果分析 (37)

第四章工程实践 (44)

4.1 工作面概况 (44)

4.2 1123（1）工作面瓦斯涌出量预测 (44)

4.3 1123（1）工作面瓦斯综合治理技术 (45)

4.3.1 风排瓦斯 (45)

4.3.2 1123（3）底抽巷抽采 (46)

4.3.3 采空区埋管抽采 (47)

4.3.4 轨运顺顺层钻孔抽采 (48)

4.3.5 初期抽采方式 (48)

4.4 1123（1）工作面瓦斯治理效果分析 (48)

4.4.1 1123（1）工作面风量与风排瓦斯的关系 (48)

4.4.2 瓦斯抽采量与抽采率随时间的变化关系 (50)

4.4.3 抽采瓦斯量与抽采浓度随工作面推进速度的关系 (55)

4.4.4 1123（1）工作面在回采期间生产以及瓦斯综合治理效果 (55)

第五章结论与展望 (58)

5.1 论文主要结论 (58)

5.2 论文后续工作的展望 (59)

第一章 绪 论

1.1 研究意义

我国是世界上的产煤大国，同时是煤矿事故发生大国。我国的煤矿事故一般

来说有以下几个特点：一是事故总量过高，这几年来，全国煤矿每年死亡人数呈下降趋势，但总数仍约为其他产煤大国死亡人数总和的3倍；二是特大事故尚未得到有效遏制。每年一次死亡人数在10人以上的特大事故约有十数起，尤其是近几年来一次死亡百人以上的特别重大事故连续发生数次，后果相当严重;三是瓦斯事故比例高，人员伤亡大，每年发生的特大事故中有80%以上是瓦斯事故，而因瓦斯事故造成的死亡人数约2000人，占煤矿事故总死亡人数的1/3。从我国煤矿每年事故统计资料来看，瓦斯爆炸事故造成的死亡人数在全部伤亡人数中所占比例呈上升趋势。因此，瓦斯事故的预防是我们煤矿安全工作的重中之重。

根据国家煤监局公布的数据，“十一五”也即2006年~2010年我国部分国

有重点煤矿事故的次数与伤亡人数的统计绘制成表1-1。

表1-1 十一五期间部分国有重点煤矿各类事故统计表

2006 2007 2008 2009 2010 起数 人数 起数 人数 起数 人数 起数 人数 起数 人数

顶板 1633 1902 1299 1518 1032 1222 805

939 702 829 瓦斯 327 1319 272 1084 182 778 157

755 145 623 运输 467 517 409 453 348 400 285

319 246 281 水害 100 417 63 255 59 263 47

166 38 224 机电 94 93 90 86 106 109 97

97 71 78 火灾 7 26 10 72 11 111 4

31 12 168 放炮 78 90 69 77 39 55 45

75 34 37 其它 239 382 209 241 177 277

176 249 155 193 由表1-1及图1-1可以得出瓦斯事故不仅发生的次数多，而且平均每起事故

的死亡人数高达4.3人。这说明瓦斯灾害，特别是瓦斯煤尘爆炸和煤与瓦斯突出，一般出现大都是重大事故，甚至是特大事故，对井下人员的生命造成巨大的威胁，并且摧毁矿井设施，造成的直接经济损失和恶劣的社会影响也是十分巨大的，同时还迫使矿井停产，投入大量的人力物力进行抢险救灾，其造成的的间接经济损

年

份 类 型

失也是难以估量的。

顶板瓦斯运输水害机电火灾放炮其他0

1000

2000

3000

4000

5000

60007000

事故起数 死亡人数

图1-1 十一五期间部分国有煤矿各类事故与死亡人数当今随着采煤机械化程度越来越高，开采深度越来越深，瓦斯方面的各种问题越来越突出，尤其是综放开采，由于其产量集中，瓦斯散发面与采空区高度相比其他采煤方式加大许多。采空区空间的加大和覆岩活动范围的增大，使综放工作面瓦斯涌出出现了一些新的特点:①在采煤机割煤移架过程中顶煤的应力状态和煤体的结构相应发生改变，并产生大量的裂隙，从而使煤体中本来的存在的大量吸附性瓦斯解析成游离状态释放到采空区，同时采空区的遗煤也会释放出一部分瓦斯，造成在采空大大增加，贯通区的瓦斯大量积存。②在工作面后方，原先在采空区上部积聚的部分瓦斯由于顶煤的冒落而涌入工作面，使采空区的瓦斯涌出量比例占综放工作面总涌出量的50%~60%，最高可达70%~80%。③采空区的瓦斯涌出是工作面和矿井瓦斯涌出的主要来源，其大量涌出经常造成工作面瓦斯频繁超限并被迫停产，甚至会酿成严重事故，造成大量人员伤亡和巨大经济损失。

近年来计算机技术发展迅速，这同时为我们研究采空区瓦斯的运移及进行

数值模拟提供了新的工具。由于采空区的环境非常复杂，不易测，因此研究采空区流场及瓦斯分布规律过程中，计算机技术提供了相当有力的条件，利用计算机对采空区进行数值模拟，模拟采空区流场及瓦斯的分布规律，为煤矿的安全生产、防治瓦斯灾害等工作提供了可靠的理论依据，对防灾技术及安全理论的研究具有重要的推动意义，同时对深入认识采空区瓦斯分布、运移、涌出规律也有一定的理论价值。数值模拟方法是场流研究的基本手段。在现实中，采空区的研究当中能够得到的现场试验研究的条件是极其有限的，而模型的理论分析和理论实验功能却能帮助我们较为真实的认识与研究采空区瓦斯分布等安全问题的内在机制

和基本规律。

1.2国内外的研究现状

1.2.1 采空区混合气体流动场研究现状

早期的采空区气体流动规律研究主要是应用并联风路模型，即把采空区和回采工作面看作是两条并联的风路来分析采空区内的漏风量。该模型对分析回采工作面后无漏风源和汇的U型通风方式具有一定的指导意义，但当漏风源汇数量较多或采空区几何形状比较复杂时不适用。

采空区管网模型和渗流模型是较早的具有代表性的采空区气体流动模型，其把采空区介质看做是纵横交错相互连通的管网，认为气体在采空区的流动可以近似认为是在众多管网中的流动。冯小平(1995)、徐军(1999)和张国枢(2002）等人以采空区空气流动属于层流和湍流之间的过渡流为基础，对采空区进行了管网化划分，建立采空区渗流场数学模型，并用网络解算的方法对采空区风流状态进行了数值模拟。他们还以能量方程和质量方程为基础对采空区的温度场和浓度场进行了模拟。刘泽功(1996)也将采空区冒落带的气体流动假设为通过许多纵横交错的短管流动，将采空区多孔介质简化为网络系统，从而根据网络模拟技术，对采空区漏风场进行模拟解算。丁文通等(2003)通过建立采空区滤流场数学模型，对采空区进行了管网化划分，并以实例用计算机模拟采空区的流场分布，应用极限风速法对采空区”三带”划分进行了研究。

1.2.2 瓦斯渗流理论的研究历程

（1）线性瓦斯渗流理论

早期在等温和低速条件下的煤层瓦斯渗流是线性渗流。研究历程大致如下图所示：

1965年周世宁院士提出采空区瓦斯渗流符合达西定律

章梦涛利用多孔介质流体动力学理论建立了采空区气体流动和瓦斯运

移的二维数学模型

1984年郭勇义和周世宁1984通过相似理论方法研究了四种一维瓦斯流

场的完全解，提出了修正的瓦斯流动方程

图1-2 线形瓦斯渗流理论研究历程图 (2) 非线性瓦斯渗流理论

1990年余楚新、鲜学福和谭学术(1990) 建立了煤层瓦斯渗流的控制方程

印度的Patrick 利用气体扩散方程和Darcy 定律建立了三维的矿区气体扩

散模型，并用有限差分法进行解算

1993年黄云飞和孙广忠讨论了提出“煤一瓦斯介质力学”的观点，并对

煤一瓦斯介质的变形、渗透率、强度等力学特性进行了系统研究

郑哲敏院士课题组阐述了煤与瓦斯突出的孕育、启动与停止过程的机

理，指出煤与瓦斯的突出机理缘于煤的破碎起动与瓦斯渗流的耦合建立

了煤层瓦斯线性渗流数学模型。

章梦涛、梁冰和赵阳升等总结以上几种理论也较系统地研究了非等温、

煤体变形及瓦斯渗流耦合力学模型及其有效应力规律。

齐庆杰和冯圣根据渗流理论，建立了三维的采空区气体流动和瓦斯运移

的数学模型及数值解算方法。

孙培德基于幂定律的推广形式得出幂定律更符合煤层内瓦斯流动的基本规律

Khristanovi 和siemek 等认为在压力梯度较大时，二项式非线性渗流定律比线性的达西定律更符合实际 E.M.Alle 证明了

瓦斯在煤岩体

中的渗流存在

不遵守线性的

达西定律，存在

非线性关系 罗新荣基于实验研究，提出了基于克林伯格效应的修正达西定律，建立

了非线性瓦斯渗流数学模型。

丁广骧等引用了大雷诺数渗流力学的Bachmat 方程，得出了采空区非线

性的渗流力学数学模型及数值解算方法

陈全利用渗流理论和质量守恒方程，建立了采空区气体三维不可压缩非

线性渗流模型，并用有限元法进行解算

辽宁工程技术大学的单亚飞、李宗翔等利用有限元方法对采空区自然发火的“三带”进行了划分，并对采空区自然发火的早期过程、综放工作面的自然发火、综放面煤柱内的漏风和耗氧过程、综放沿空巷周围煤体自燃升温过程、以及采空区温度场温度分布问题等煤自然发火各个方面进行了数

值模拟

图1-3 非线性瓦斯渗流理论的研究与早期数值模拟情况

近来，更多研究利用FLUENT软件包来求解采空区气体分布规律。肖肠和朱毅通过实验和现场观测，拟合出孔隙率等参数的分段函数，建立了采空区漏风松散体空隙分段分布方程、采空区不同深度的分段渗透系数方程等自然发火预测模型，并利用FLUENT解算。Yuan Liming和A.c.smith根据小尺寸自燃实验数据，数值模拟了遗煤自燃下的考虑老空区影响的采空区的气体流动和传热情况，研究了其温度场和氧气浓度场的分布规律。王俊峰和邬剑明等在示踪气体试验的基础上，建立了模拟综采工作面采空区空气、瓦斯流动的数学物理模型，利用FLUENT解算了采空区氧气和瓦斯的分布规律。得出CFD模型模拟结果与现场长期监测结果相吻合，高瓦斯涌出可以加速采空区的自我惰化，有利于防治采空区自燃。

1.3研究思路

1.3.1 研究内容

在查阅和分析过内外研究资料和文献的基础上，总结采空区气体流动的相关研究成果和存在的问题，运用多孔介质流体力学、渗流力学等学科理论，结合顾桥煤矿综采工作面采空区的实际环境，利用FLUENT软件对采空区瓦斯的渗流情况进行数值模拟。

具体内容有：

1. 综放工作面采空区瓦斯来源构成及其影响因素进行研究。

2. 建立综放工作面采空区的物理及数值模型。

3. 以碎胀特性推导渗流特性，将采空区空间视为非均匀孔隙介质空间，用FLUENT软件建立采空区瓦斯运移的三维数值模拟，研究采空区混合气体的分布以及沿工作面走向、倾向和竖向三个不同方向上的浓度的分布规律。

4. 以顾桥矿1123（1）工作面的通风与瓦斯抽放数据基础，并结合模拟结果，对该工作面已经采取的瓦斯治理措施的效果进行分析。

1.3.2 研究方法

综合应用理论分析和数值模拟技术和现场分析等研究方法，系统研究采空区瓦斯——空气混合气体不同组分的分布情况，以总结出采空区气体三维渗流场、浓度场和温度场的多场耦合的普遍规律。

(1) 理论分析

充分吸收、利用和发展了多孔介质渗流力学、计算流体力学等科学理论和方法，并结合现场资料和现场观测数据及学科前沿对采空取瓦斯的渗流规律与运移理论进行研究。

(2) 数值模拟

利用数值模拟技术，建立1123(1)工作面采空区混合气体的三维流动模型。通过CFD解算和后处理，进一步形象地揭示采空区气体三维渗流的浓度分布规律。

(3) 现场实测

顾桥煤矿1123（1）工作面可采走向长度2754.6m，面长265m，2012年6月初开始试采，剩余可采走向1630m。工作面1440m以外处于-830m以下。预计2013年底回采结束。搜集该工作面从2012年7月份到2013年4月份的风量及瓦斯抽采数据。为理论分析和数值模拟提供一手资料，也用于验证瓦斯综合治理效果和数值模拟结果。

1.3.3 技术路线

以采空区上覆岩层区带划分为依据，以各区带内岩体的碎胀系数为特征，建立采空区瓦斯渗流数学模型，利用FLUENT 建立瓦斯渗流数值模型，通过数学模型给出的边界条件确定数值模拟的边界条件，进行迭代计算，得到综放采空区瓦斯渗流的数值解，本文的技术路线见下图。

FLUENT 模拟得出的结果现场实际得到的瓦斯抽采数据

两者基本吻合

对顾桥1123（1）工作面已经采取的瓦斯治理措施进行效果分析给出评价

论文

图1-4 技术路线图

第二章综采工作面瓦斯运移理论研究与数值模型

2.1 综采工作面基本特征

采空区是人类进行地下矿层开采活动而形成的空间。我国有很多学者长期以来致力于采空区基本特征的研究工作，付出了大量艰辛劳动，同时也取得了颇丰的成绩。其中最为著名的是钱鸣高院士提出的采场上覆岩层的“砌体梁”结构模型，钱院士在总结预成裂隙假说及铰接岩块假说的基础上，又经过长期的大量的对岩层内部移动特征的现场观测和生产实践活动，在20世纪70年代末80年代初期提出了岩体结构“砌体梁”力学模型，从而发展了以上学说。近年来随着对岩层控制科学研究的不断深入以及为了解决岩层活动中更为广泛的问题，如岩层移动过程中的卸压瓦斯抽放与突水防治、离层区充填与地表沉陷控制等，在砌体梁理论的基础上，进一步提出了“岩层控制关键层理论”。“砌体梁模型”如下图所示，该图表示回采面前后的岩体形态其中A为煤壁支撑影响区，B为离层区，C为重新压实区；I为冒落带，Ⅱ为裂隙带，Ⅲ为弯曲下沉带。“横三区”和“竖三带”的提出为研究瓦斯在采空区的运移提供了不同的运移介质和较为科学的理论依据。

A-煤壁支撑影响区；B-离层区；C-重新压实区；

I -冒落带；Ⅱ-裂隙带；Ⅲ-弯曲下沉带

图2-1 砌体梁结构模型中的“横三区”与“竖三带”

2.2 综采面采空区瓦斯来源及涌出影响因素

2.2.1 综采面采空区瓦斯来源

在人类对煤层进行开采活动之前，原始的瓦斯气体与煤岩体受到地层压力、瓦斯压力等作用下，所组成的系统处于一个动态平衡状态。在煤层进行开采时，由于受采掘活动的影响原始的平衡被打破。随着工作面的不断推进，采空区上覆煤岩层发生形变、断裂，形成裂隙，裂隙在相互沟通后，构成了一个动态变化的采动裂隙带，这就为瓦斯的运移、储集提供了便利的运输通道。采空区内由于漏风的存在、再加上压差的驱动，贯通的孔隙空间将充满流动瓦斯或瓦斯一空气混合气体。同时也为开采煤层上下邻近层及围岩中的瓦斯涌出提供了通道，大量的瓦斯在浓度压力差的驱动下通过该贯穿的通道以扩散或渗透的形式源源不断的涌向开采煤层回采工作面采空区，以至达到新的平衡。

综采工作面采场范围内涌出瓦斯的地点称为瓦斯涌出源。瓦斯涌出源的多少，各源涌出瓦斯量的大小直接影响着采场的瓦斯涌出量。工作面的主要瓦斯来源有三个方面主要是煤壁瓦斯涌出、采空区瓦斯涌出、落煤瓦斯涌出。而采空区瓦斯涌出源可以分为遗煤瓦斯涌出、围岩瓦斯涌出、未采分层瓦斯涌出和邻近已采老空去瓦斯涌出。如图2-2所示。各部分涌出瓦斯随着采场内煤岩层的变形或跨落而卸压，涌入采空区并混合在一起形成采空区涌出瓦斯。

图2-2 工作面瓦斯来源构成示意图

(1) 煤壁瓦斯涌出

煤层开采后，由于在矿山压力的作用下，老顶的初次来压以及基本顶的周期来压，支架上方顶煤的应力平衡状态受到破坏，同时出现了卸压带，煤的透气性系数大大的增加，使得回采工作面暴露的煤壁及煤壁位移卸压带的瓦斯解吸进入风流当中。当综放开采时，瓦斯通过卸压带的裂隙从顶煤壁涌入顶板的断裂带。

瓦斯涌出的强度随着煤壁暴露时间的增加而逐步的降低，煤壁的瓦斯释放与割煤速度之间的关系成近似线形关系。煤壁瓦斯涌出量主要取决于煤层原始瓦斯压力、透气性及工作面推进速度等因素。

(2) 采空区遗煤瓦斯涌出

采空区遗煤的瓦斯涌出主要是残留在采空区的放落煤炭继续释放的瓦斯，其

主要由煤层的采出率所控制并随时间的推移逐渐减少。采空区气体渗流规律数值模拟。

(3) 采落煤瓦斯涌出

随着回采工作面不断向前推进，工作面煤壁上的煤被不断采落、运出，其所

含有的瓦斯也源源不断地释放到回采空间。采落煤呈块粒状，使得煤体的暴露面积成倍增加，从而提高了瓦斯解吸强度和速度，并导致瓦斯涌出量的增加。不同块粒度采落煤的瓦斯放散速度是不同的，粒度越小，瓦斯放散速度越快；粒度越大，瓦斯放散速度越慢。采落煤瓦斯涌出量主要取决于采落煤量、原始瓦斯压力、透气性、煤体的坚固系数等因素，其涌出强度与煤壁相似，也随着时间的增加而降低。

2.2.2 瓦斯涌出量计算

瓦斯涌出量是指在矿井建设和生产过程中从煤与岩石内涌出的瓦斯量。其表达方式有两种包括绝对瓦斯涌出量和相对瓦斯涌出量。绝对瓦斯涌出量指在单位时间内涌出的瓦斯量，单位为m 3/min 年相对瓦斯涌出量指平均日产一吨煤同期所涌出的瓦斯量，单位为m 3/t 。

在采煤工作面，瓦斯的涌出主要包括两部分:本煤层瓦斯的涌出与邻近层瓦

斯的涌出，即Q=Q B +Q L 。当采煤工作面匀速向前推进时本煤层瓦斯涌出量可以按下列两种方法来计算：

① 瓦斯在煤层中的含量来计算

01B Q 1440

H （x -x ）mV(l-l ) (2-1) 在式中：Q B —— 瓦斯从本煤层涌入工作面的量，m 3

/min ； X 0 —— 本煤层瓦斯原有的含量，m 3/t ；

X l —— 残余瓦斯的含量，m 3/m 3 ；

m —— 煤层厚度，m ；

V —— 工作面平均推进度，m/d ；

l —— 工作面长度，m ；

l H 一— 进风巷和回风巷瓦斯排放带的长度，m 。

② 按涌出的规律来计算：

12B B B Q Q Q =+ (2-2)

110(1)()1440

bt H B x e mV l l Q ---= (2-3) 1202[1(1)]()1440

bt n H B x e t mV l l Q ---+-= (2-4) 在式中： Q B1 ——工作面煤壁瓦斯的涌出量，m 3

/min Q B2 ——落煤瓦斯涌出量，m3 /min

B ——系数 d-1

t 1 ——煤壁寿命，d ， t 1=B/V

S ——采煤循环的一个进度，m 。

n ——系数。

.将（2-2）（2-3）带入（2-4）中得

102H [1(1)]()1440

bt n B x e t mV l l Q ---+-= (2-5) 根据上式(2-3)、（2-4)和(2-5)，煤壁和采落煤的瓦斯涌出量会随着回采工作

面向前推进速度的加快而相应的增加，但是相对于煤壁瓦斯的涌出量，采落煤炭瓦斯涌出量增加的速度更快，说明在开采层瓦斯的涌出量中，采落煤炭的绝对涌出量所占的比例增大，而煤壁瓦斯的绝对涌出量占的比例会减少。

在邻近层瓦斯涌出量的计算公式是：

L 01

Q =1440n

k k k k lV x n m =∑ (2-6) 在上式中: Q L —— 上、下邻近层涌入开采层的瓦斯量，m 3/min;

x 0k —— 第k 邻近层瓦斯原始的含量，m 3/m 3

m k —— 第k 邻近层的厚度，m

n k —— 第k 邻近层瓦斯涌出率， 00k k k k

x x n x -= x k —— 第k 邻近层残余瓦斯含量，m 3/m 3。

由(2-6)可得，邻近层瓦斯的涌出量与工作面推进的速度成正比，即工作面推

进速度加快，邻近层瓦斯涌出量增大。

2.2.3 影响采空区瓦斯涌出的基本因素

(1) 自然因素

①在煤层与围岩中瓦斯的含量

在这里瓦斯含量的多少直接会影响瓦斯的涌出量。当单一的薄煤层和中厚煤层开采时，煤层暴露面和采落的煤炭是瓦斯的主要来源，开采时的瓦斯涌出量随着煤层随着煤层瓦斯含量的增大而增大。在煤层开采的过程中，如果附近存在煤层或岩层的瓦斯含量比较大时，因为受到煤层回采的影响，在采空区就会形成许多的裂隙，这些在煤层中的瓦斯，就会不停地向开采煤层的采空区流去，在流进生产空间，使得矿井的瓦斯涌出量增大。在这种情况下，煤层开采时的瓦斯涌出量要比它的瓦斯含量大的多。

②地面大气压的变化

在一年内的夏冬两季，地面大气压的差值可达到5.3~8kPa，在个别的情况下，一天之内的气压差值可达到2~2.7kPa，而地面大气压的改变会引起井下大气压发生相应的变化。它对瓦斯在采空区或坍冒处的涌出的影响比较明显。如果地面的大气压突然下降时，风流的压力将低于瓦斯积存区内气体的压力，引起了涌入风流中的瓦斯增多，造成工作面瓦斯的涌出量增大。反之，工作面瓦斯的涌出量将减少。

(2)开采技术因素

①开采顺序与回采方法

首先开采的煤层(或分层)，其相对瓦斯涌出量增大，而后开采的煤层(或分层)，其涌出量减少。回收率低的回采方法，相对瓦斯涌出量增大。陷落式顶板管理方法比充填式造成更大范围的围岩破坏与卸压，邻近层瓦斯涌出的分量增大。因此，陷落式比充填式的相对涌出量高。水采水运的采煤方法，比旱采相对瓦斯涌出量降低，这是因为湿煤残余瓦斯含量增大的缘故。

②配风量

随着风量的增加，会使采空区漏风量增大，这样漏入的风量会带出更多的瓦斯，使采面瓦斯涌出量增大，造成上隅角瓦斯超限。而当配风量小时，稀释瓦斯的能力有限，容易造成回风流及上隅角瓦斯超限。可见，合理配风对控制采空区瓦斯涌出具有重要的作用。

③生产工艺

从采落的煤炭与煤层的暴露面内涌出瓦斯的特征是，瓦斯在初期涌出的强

度会很大，随后大致的方向是依照指数函数的曲线逐渐减少，因此在落煤时瓦斯的涌出量总是比其它工序要大。当落煤时，瓦斯涌出量的增加值与落煤量、新显露煤面的大小和煤块的破碎程度有关。比如：在风镐落煤时，涌出瓦斯的数量可增大1.1~1.3倍；在放炮时可以增大1.4~2.0倍。由于综合机械化工作面推进速度较快，产量较高，当在瓦斯含量大的煤层内工作时，瓦斯的涌出量也会很大。 ④ 工作面推进速度的影响

当工作面推进速度越快，单位时间内工作面卸压的范围越大，绝对瓦斯涌出

量就越大，反之则越小。此外，推进速度的改变意味着产量也发生变化，则相对瓦斯涌出量也发生变化。推进速度快的相对瓦斯涌出量小，推进速度慢的相对瓦斯涌出量大。

以顾桥煤矿1123（1）工作面为例，图2-3为日进尺与瓦斯涌出的关系。

2012/6/262012/7/32012/7/102012/7/172012/7/242012/7/31

20

40

6080100

120瓦斯涌出量 (m 3/m i n )推 进 距 离 (m )时间(yy-mm-dd)

3040

50

6070

图2-3 工作面推进速度与瓦斯涌出量的关系 当工作面沿11-2煤层从开切眼向采区边界推进时，在工作面老顶第一次冒

落以前的时间内，采空区的瓦斯涌出量仍然是以留在采空区内的遗煤涌出的瓦斯为主。当工作面老顶第一次冒落时，就会从卸压后的邻近煤层和岩层向采空区涌出大量的瓦斯，采空区的瓦斯涌出量将显著增加。并且，在工作面继续推进到老顶下一次冒落前。瓦斯涌出量一般情况下将逐渐减少，但仍大于老顶第一次冒落前的瓦斯涌出量。再以后发生顶板冒落时，将会重复出现上面的过程，于是采空区中瓦斯涌出量将逐渐增大。随着工作面推进到一定程度，采空区深度也推进到

一定程度，采空区深部的瓦斯涌出量达到某一极限位并趋于稳定。

⑤通风压力与采空区封闭质量

通风压力与采空区密闭质量都对老空区的瓦斯涌出有一定的影响。通风压力小，采空区密闭质量好，可减少老采区瓦斯涌出不均匀系数及涌出量，这对老矿井具有很大的意义。

⑥采空区内部与采面间的压差

气体总是从压力高处向压力低处流动，采空区深部瓦斯浓度较高，所以瓦斯从采空区深部不断向采面运移，从支架后涌入工作面，甚至造成采面或上隅角瓦斯超限。控制采空区内瓦斯或有毒气体涌出，一般采用均压法，即增加工作面的绝对压力，使工作面和采空区间的压力达到平衡，这样把将涌入到工作面的瓦斯或有害气体“顶”回去，其在控制采空区瓦斯涌出和均压抽放采空区瓦斯等方面有较广泛的应用。

⑦瓦斯抽放的影响

随着综放工作面开采采高的扩大，使采空区的垮落高度也随之加大，扩大了上邻近层瓦斯的排放范围，从煤体中解吸出来的瓦斯不断的涌入综放工作面。通过引巷抽放、埋管抽放、高位钻孔抽放、高抽巷抽放、注浆道抽放等多种抽放方法，使工作面的瓦斯涌出量减少。因此，通过采取瓦斯抽放方法，能有效的降低工作面的瓦斯涌出量。

2.3 采空区瓦斯运移特征

2.3.1 多孔介质的特性

简单的说，多孔介质是指含有大量孔隙的固体。一般认为富含孔隙和裂隙的煤体是多孔介质，具有如下特点：

(1) 多相性。指可以同时存在固相、液相和气相或同时存在固相和液相、固相和气相。固相部分称为固体骨架，固体骨架以外的部分成为多孔介质的孔隙空间。

(2) 在多孔介质所占据的范围内，固体骨架遍布于整个多孔介质中。相比于固体骨架具有很大的比面积，孔隙却比较狭窄。

(3) 孔隙空间中只有少数的空隙是封闭的，绝大多数空隙是相互连通的，流体在连通的孔隙中流动。

1、孔隙性

孔隙的大小和多少是多孔介质的重要特征，通常用孔隙率这一指标来衡量多

孔介质中孔隙的多少。孔隙率是孔隙体积与多孔介质总体积之比，其计算公式为：

n 1S S P V V V V V V V

-===- (2-7) 式中：n ──孔隙率；

V ──多孔介质总体积

V P ──孔隙体积，

V S ──固体骨架体积，

孔隙率n 是量纲为1的量，与颗粒大小无关，只和颗粒的排列方式有关。采

空区孔隙率目前都是按照顶板岩石的碎胀系数计算，即：

11p

n K =- (2-8) 式中：K p ──采空区岩石碎胀系数。

2、渗透性

多孔介质的渗透性既可以用渗透率来衡量，即渗透率是多孔介质对流体的

渗透能力，由达西定律给出；又可以用渗透系数来衡量。渗透率仅与多孔介质的组成结构(颗粒大小、颗架分布、颗粒充填、比面、弯曲度和空隙率等)有关，与其中的流体性质(流体的密度、粘滞性等)无关；而渗透系数则不但与多孔介质的组成结构有关，而且还与其中的流体性质有关。

现有文献中关于渗透率的计算公式，多为经验性的和半经验性的公式，下

面将介绍一种有代表性的计算公式。由多孔介质的Carman 公式,渗透系数的计算公式为：

23

__

2180(1)m D n K n =- (2-9) 式中：式中：m D ──多孔介质骨架的平均粒径，m

n ──多孔介质的孔隙率。

3、比面

比面__M 定义为单位体积多孔介质中的总表面积，它的量纲是长度的倒数。可表示为：

___A M =

V

(2-10) 式中：M ──比面，1cm -

V ──多孔介质的体积, 3cm

A ──体积V 中的总表面积, 2cm

如果多孔介质是由若干半径为r 的圆球堆积排列构成，则其比面可以表示为：

2344/3/3

M r r r ππ== (2-11) 对于有m 种不同半径球形颗粒所组成的多孔介质，其中半径为r i 的圆球个数为个，则颗粒的总面积A 和总体积V 分别表示为：

21

4m i

i A r N π==∑ (2-12)

3143m i i i V r N π==

∑ (2-13) 其比面___M 表示为：___221144/3

m m

i i i i i i M r N r N ππ===∑∑ (2-14)

4、压缩性

多孔介质的压缩性用体积压缩系数来衡量。体积压缩系数指，当温度一定

时，作用于多孔介质上的压强变化一个单位而引起的整个体积的变化率，其计算公式为：

1dV

V dP ?=- (2-15)

式中： ?──体积压缩系数，1MPa -；

V ──多孔介质的总体积，3cm

P ──作用于多孔介质上的压强，MPa 。

2.3.2 采空区瓦斯的升浮

根据环境流体力学理论，气体升浮产生的条件有两种：一是气体由于受热体积膨胀，密度减小，从而产生密度差；二是气体对象中的含有物浓度相对周围气体中含有物浓度存在差异，从而产生密度差。采空区瓦斯升浮的根本条件是局部对象相对于它周围的气体存在着小于周围气体密度的密度差。

瓦斯的主要成分是甲烷，它是无色、无味、无臭、可以燃烧或爆炸的气体。

甲烷的密度为0.716kg/m 3 (标准状况下)，为空气密度的0.554倍，所以矿井内空气中如果有瓦斯聚集体的存在，则会因其本身与周围气体的密度差而会上升并漂