《教学课件：第三章-煤层瓦斯流动理论.ppt》由会员分享，可在线阅读，更多相关《教学课件：第三章-煤层瓦斯流动理论.ppt（36页珍藏版）》请在三一办公上搜索。
1、安全工程专业,中 国 矿 业 大 学安 全 工 程 学 院煤矿瓦斯治理国家工程研究中心,矿 井 瓦 斯 防 治,程远平 王亮,第三章 煤矿瓦斯流动,一、煤层瓦斯流场分类 二、瓦斯流动理论及其模型 三、瓦斯流动模型的应用,煤层瓦斯流场分类,1,煤层瓦斯流动方式,孔隙直径为10-510-4，瓦斯流动表现为自由扩散或慢速的层流渗透,孔隙直径为10-4210-4 cm，瓦斯流动为层流渗透,孔径大于 210-4 cm，瓦斯流动表现为层流渗透或层流与紊流的混合过渡流,孔隙直径小于10-5 cm，瓦斯流动属于分子扩散,孔径小于310-7cm，出现瓦斯表面扩散和固体中的扩散,一、煤层瓦斯流场分类,煤层瓦斯流动
2、,单向流动示意图1-流线；2-等压线；3-巷道,单向流动,煤层瓦斯沿单一方向流动，流线相互平行,一、煤层瓦斯流场分类,煤层瓦斯流动形态,径向流动示意图1-流线；2-等压线；3-钻孔,径向流动,等瓦斯压力线为一组同心圆，瓦斯流线沿圆的径向发展,一、煤层瓦斯流场分类,煤层瓦斯流动形态,球向流动示意图1-揭开煤层的掘进工作面；2-等压线；3-流线,球向流动,等瓦斯压力线为一组同心球状，瓦斯流线沿球的径向发展,一、煤层瓦斯流场分类,煤层瓦斯流动形态,上述三种流动是按照流场的空间流向分类的基本形式。在实际煤矿中，由于煤层的非均质性、煤层顶底板岩性的多变性等自然条件的不同，实际井巷和钻孔中的瓦斯流动是复杂
3、的，有时可能是几种基本流动的综合。,一、煤层瓦斯流场分类,煤矿常见的流场分类,按流场的稳定性分类,一、煤层瓦斯流场分类,2成面后4d,1成面后几小时,3成面后10d,4成面后15d,5面后55d,6成面后150d（稳定）,煤层内瓦斯流动的稳定过程就是煤壁瓦斯涌出速率的稳定过程，这一过程受到多种因素的影响，主要有煤层瓦斯压力梯度、煤层的透气性、煤壁外的大气压力等，其中最关键的一个是煤层的透气性。,一、煤层瓦斯流场分类,按流场的稳定性分类,瓦斯流动理论及其模型,2,二、瓦斯流动理论及其模型,2.1 瓦斯流动理论,煤层瓦斯流动理论是专门研究煤层内瓦斯压力分布及瓦斯流动变化规律的理论，根据应用范围和使
4、用条件的不同，煤层瓦斯流动理论有以下几种。,2.1.1 线性瓦斯流动理论,线性瓦斯渗流理论认为，煤层内瓦斯运移基本符合线性渗透定律达西定律（Darcys law）,达西定律:,v流速，m/s瓦斯动力粘度系数，Pa sK 煤层的渗透率，m2dx 与流体流动方向一致的极小长度，mdp 在dx长度内的压差，Pa,二、瓦斯流动理论及其模型,2.1.2 瓦斯扩散理论,煤是一种典型的多孔介质，根据气体在多孔介质中的扩散机理的研究，可以用表示孔隙直径和分子运动平均自由程相对大小的诺森数Kn：,d孔隙平均直径，m气体分子的平均自由程，m,扩散分类,二、瓦斯流动理论及其模型,对于含瓦斯煤体来说，一般Kn 10，
5、由于孔隙直径远大于瓦斯气体分子的平均自由程，因此扩散是由于瓦斯气体分子之间的无规则运动引起的，可以用菲克扩散定律去描述。,J瓦斯气体通过单位面积的扩散速度，kg/(s m2)C/X沿扩散方向的浓度梯度Df 菲克扩散系数，m2/sC 瓦斯气体的浓度，kg/m2,菲克定律:,瓦斯气体分子在煤层内有效扩散系数：扩散路径因孔隙通道的曲折而增长，孔截面收缩使扩散流动阻力增大，实际扩散通量减少。,J瓦斯气体在煤层内的有效扩散系数，m2/s有效表面孔隙率 曲折因子，为修正扩散路径变化而引入的。,二、瓦斯流动理论及其模型,2.1.3 瓦斯渗透-扩散理论,瓦斯渗透与扩散理论认为，煤层内瓦斯运动是包含了渗透和扩散
6、的混合流动过程。煤层中存在相互沟通的裂隙网络，沿着这些裂隙网络，游离瓦斯流向低压工作面，而煤体的渗透率与该裂隙网络密切相关。与此同时，块煤内部的瓦斯解吸，向裂隙扩散，因此煤层中瓦斯的渗透率和介质的扩散性共同决定了瓦斯的流动状况。,二、瓦斯流动理论及其模型,2.1.4 非线性瓦斯流动理论,低雷诺数区，Re100，为紊流，惯性力占优势，流动阻力和流速的平方成正比。,多孔介质中流动状态及其规律的变化,二、瓦斯流动理论及其模型,2.1.5 地物场效应的煤层瓦斯流动理论,随着煤层瓦斯流动机理研究的深化，许多学者用流体岩石相互作用机制认识煤层内瓦斯运移过程，充分发展和考虑地应力场、地温场以及地电场等地球物
7、理场作用下的煤层瓦斯运移耦合模型及数值方法，使理论模型更能反映客观事实，以及进一步完善理论模型及测试技术和手段，成为当今推动煤层瓦斯渗流力学向前发展的主流方向。,二、瓦斯流动理论及其模型,2.2 瓦斯流动方程,2.1.1 基本假设,煤层中的原始瓦斯压力，温度和含量相同；煤层的透气性系和孔隙率不受煤层中瓦斯压力变化的影响，且处处相等；煤层瓦斯流场内，温度保持不变，并服从理想气体状态方程；瓦斯为理想气体，瓦斯在煤层中的流动为层流渗透，且服从达西定律；煤层瓦斯含量有游离瓦斯和吸附瓦斯组成，可以用抛物线方程近似计算；瓦斯的吸附、解吸是动态平衡，瓦斯的吸附、解吸是瞬间完成的；由于煤层顶底板透气性与煤层相
8、比要小得多，可将其视为不透气层，且其中不含瓦斯。,二、瓦斯流动理论及其模型,2.1.2 数学模型,二、瓦斯流动理论及其模型,2.1.2.1 煤层瓦斯流动的连续性方程,根据质量守恒原理：,微分单元体各面上进出流量示意图vx、vy、vz分别为v在各坐标轴的分量,考虑x方向在dt时间内的质量变化,同理：,dt时间内微元总的质量变化为：,二、瓦斯流动理论及其模型,dt时间内微元总的质量变化又可以表示为：,可以得到：,化简为：,瓦斯压力为p时的密度，kg/m3v 瓦斯渗流速度，m/sM煤层瓦斯含量，kg/m3,二、瓦斯流动理论及其模型,2.1.2.2 煤层瓦斯状态方程,假设煤层瓦斯为理想气体，由理想气体
9、状态方程可得：,瓦斯压力为p时的密度，kg/m3n瓦斯压力为pn时的密度，kg/m3pn标准状态下的瓦斯压力，m/sp瓦斯压力，MPa,二、瓦斯流动理论及其模型,2.1.2.3 煤层瓦斯含量方程,为了计算方便，采用抛物线方程近似取代煤层瓦斯含量：,M煤层瓦斯含量，kg/m3瓦斯含量系数，m3/(m3Pa0.5)p瓦斯压力，MPa,二、瓦斯流动理论及其模型,2.1.2.4 煤层瓦斯运动方程,基于达西定律经过数学推导得出煤层瓦斯线性流动的运动方程：,+,+,+,代入连续性方程,煤层透气性系数，m2/(MP0.5d)瓦斯含量系数，m3/(m3Pa0.5)p瓦斯压力，Mpa；pn 标准状态下的瓦斯压力
10、，Mpak 煤层的渗透率，m2 瓦斯的动力粘度，Pas，对于甲烷气体=1.08106 Pas,二、瓦斯流动理论及其模型,2.1.2.5 初始条件和边界条件,初始条件,煤层中的瓦斯压力分布或煤层中的瓦斯浓度分布,边界条件,第一类边界条件,待求的函数p(x,y,z,t)在边界上为已知函数，即:,第二类边界条件,待求的函数p(x,y,z,t)在边界上为未知函数，但边界上的流速是已知的，即:,第三类边界条件,待求的函数及其导数在边界上均未知，但其关系是己知的，相应的表达式为：,二、瓦斯流动理论及其模型,2.1.3 瓦斯流动模型的简化,2.1.3.1 一维瓦斯流动方程,煤层瓦斯沿垂直于煤壁的方向流动，可
11、以将其视为一维的平行瓦斯流动。,瓦斯流动控制体示意图,为了便于计算进行如下转换,初始条件：,边界条件：,二、瓦斯流动理论及其模型,2.1.3.2 二维瓦斯流动方程,无论是煤层穿层孔还是顺层孔，钻孔周围的煤层瓦斯均沿径向流入钻孔。,a-钻孔半径；b-钻孔瓦斯流动的影响半径；r-瓦斯流动的半径；dr-研究单元的厚度；H-煤层厚度图3-9 瓦斯径向流场示意图,直角坐标系下二维瓦斯流动方程,进行极坐标转换,瓦斯流动模型的应用,3,三、瓦斯流动模型的应用,3.1 瓦斯一微流动方程的求解,对式（3-20）应用拉氏变换，并将边界条件代入，整理得：,其中,上式即为煤层瓦斯一维平行流场的瓦斯压力分布规律，可以看
12、出瓦斯压力的变化主要取决于离煤壁的距离和煤壁暴露时间，由该公式可以推导出许多有用的关系式。将上式进行微分，则,对于巷道煤壁处，而瓦斯流动规律符合达西定律：,此式为均质无限煤层巷道煤壁单位面积不稳定瓦斯涌出量的方程式。,三、瓦斯流动模型的应用,3.1 钻孔周围瓦斯流动的解算（有限差分法）,简介基本思想,按时间步长和空间步长将时间和空间区域剖分成若干网格，把表示变量连续变化关系的偏微分方程离散为有限个代数方程，然后解此线性代数方程组，以求出溶质在各网格结(节)点上不同时刻的浓度。,基本步骤,剖分渗流区，确定离散点，按某种几何形状剖分成网络系统；建立瓦斯渗流问题的差分方程组；求解差分方程组，采用各种
13、迭代法。,三、瓦斯流动模型的应用,几种导数的差分近似,三、瓦斯流动模型的应用,3.3 钻孔径向流方程的有限差分法格式,对于微分方程(3-22)，采用既稳定、精度又高的六点显示差分格式，可较为准确地反映流动的最初阶段和钻孔附近的情况，时间和空间均采用由小到大的变步长差分网格，见图,三、瓦斯流动模型的应用,3.3 钻孔径向流方程的有限差分法格式,三、瓦斯流动模型的应用,钻孔瓦斯流动方程解算程序的应用,根据前面推导的有限差分方程式（3-38），在matlab平台上编制了求解钻孔瓦斯流动的计算机程序。该程序可全屏幕输入各已知参数，包括煤层的原始瓦斯压力、透气性系数、吸附常数，并模拟抽采钻孔周围煤层中瓦斯流动状况，计算一段时间范围内钻孔周围煤层瓦斯涌出参数、最佳抽采参数及煤层瓦斯压力与孔壁距离、抽采时间的关系。,三、瓦斯流动模型的应用,计算结果,三、瓦斯流动模型的应用,计算结果,
