第29卷 第1期地质科技情报Vol. 29 No. 1
2010年 1月Geological Science and Technology Information Jan. 2010
煤层气渗流理论及其研究进展
田智威
(中国地质大学资源学院煤及煤层气工程系, 武汉430074)
3
摘 要:综述了煤层气渗流研究的现状及取得的研究进展, 包括线性渗透、线性扩散、扩散—渗透以及非线性渗流等基本理论, 基于各种地球物理场效应影响的渗流理论及其修正模型, 以及多煤层瓦斯越流相关的研究; (格子Boltzmann 方法与分子动力学方法) 在煤层瓦斯渗流模拟中的应用; 关键词:煤层气; 渗流; 数值模拟; 研究进展中图分类号:P618.11 文献标志码:A2) 0120061205 应用的科学。自]律以来, , 并逐渐与其他学科交叉, 在能源资源的开发与利用以及工程建设中得到了非常广泛的应用。渗流力学最先应用在水利工程和地下水资源开发等领域; 随后又成为石油和天然气工业的一项基础理论。
随着煤层气这一新型清洁能源的重视与开发[223], 渗流理论又应用在煤层气的开发与利用中[326]。煤层气渗流力学是研究煤层内瓦斯压力分布及其流动变化规律的理论[427], 是由渗流力学、煤地质学、固体力学及采矿学等学科互相交叉渗透发展形成的。自煤层气渗流力学创立至今深受有关研究人员的关注, 尤其自20世纪80年代以来发展更为迅速, 表现在:应用范围更广; 基本理论不断深化; 研究手段及方法不断现代化。
式中:v 为流速; k 为煤层的渗透率; μ为瓦斯黏度系数; p 为瓦斯压力; 5p/5n 为瓦斯压力在流动方向上的偏导数。
此外还导出了瓦斯流量方程[式(2) ]:
(2) q =-λn
式中:q 为瓦斯流量; λ为煤层透气系数。
20世纪80年代, 多位研究者[9214]在修正和完善数学模型、流动方程方面开展了相应的工作。由于大多数井下瓦斯流动都可简化为一维的平行流动和径向流动的有限流场、无限流场或其组合, 为此, 郭勇义等[9]针对一维流动, 结合相似理论提出了修正的流动方程。谭学术等[10]利用渗流理论和热力学理论, 提出了应用真实气体状态方程(而不是理想气体状态方程) 的煤层气渗流方程, 并做了线性化简化处理。余楚新等[11]提出参与渗流的瓦斯量只是可解吸的那一部分, 在瓦斯吸附与解吸过程完全可逆的假设条件下建立了相应的控制方程, 做出了有一定新意的探索。随后, 孙培德[12214]又进一步修正和完善了线性流动数学模型。
随着电子计算机和数值模拟方法的发展, 针对煤层气流场内的压力分布及其流动规律进行数值模拟已经成为研究的主要手段。针对压力是煤层气流动的主要动力因素, 李英俊等[15]重点模拟研究了煤层瓦斯压力的分布规律, 为预测瓦斯涌出量及抽放效果提供了参考依据; 魏晓林[16]采用有限差分法对瓦斯流场中的压力分布以及流量变化进行了数值模拟; Yu 等[17218]分别采用有限元法(FEM ) 和边界单元法(B EM ) 实现了对瓦斯渗流流动的数值模拟。
1 线性渗流理论
1. 1线性渗透理论
为了适应采矿采煤业的大力发展, 控制瓦斯技术已成为当时研究的关键技术之一, 早在20世纪40年代末, 前苏联学者就已经建立起考虑吸附煤层瓦斯作用的瓦斯控制方程[式(1) ]。在我国, 周世宁院士等[8]首先进行了将达西定律应用于煤层瓦斯流动理论的开拓性研究, 认为煤层瓦斯的流动基本符合线性渗流规律, 其观点对煤层瓦斯渗流的应用和瓦斯动力学研究具有相当重要的指导意义。
(1) v =-μn
3收稿日期:2009202216 编辑:禹华珍
基金项目:中国地质大学(武汉) 优秀青年教师资助计划项目(CU GQNL0904) ; 煤燃烧国家重点实验室开放基金项目(FSKLCC0815)
作者简介:田智威(1980—) , 男, 主要从事煤层气渗流理论与数值模拟的教学与科研工作。
62
地质科技情报
2010年
1. 2线性扩散理论
瓦斯扩散理论在欧美国家有较多研究, 而在我
国这方面的研究较少[4], 其扩散规律的理论依据是Fick 扩散定律[式(3) ]:
(3) J =-D
n
式中:J 为单位面积的扩散速度; D 为瓦斯扩散系数; c 为煤体瓦斯浓度。
然而瓦斯的涌出过程是一个很复杂的过程, 杨其銮等[19]将煤层看作微孔隙裂隙体来讨论瓦斯的扩散过程, 认为其属于气体在多孔介质中的扩散, 其涌出规律符合扩散定律。依据扩散运动与热传导类似的原理, 可以借助热传导的相关研究方法来讨论瓦斯的扩散问题, 探讨和实测对比分析。1. 3扩散—瓦斯扩散—, 煤层内瓦斯的运动是包含了渗透与扩散的混合流动过程。煤层中存在相互沟通的裂隙网络, 游离瓦斯沿着这些裂隙流向低压工作面; 同时煤块内部的瓦斯解吸附而向裂隙扩散, 因此, 煤层中瓦斯的渗透率及其扩散性共同决定了其流动状况。
随着瓦斯运移规律研究的深入, 国内外许多学者都赞同这一理论, 周世宁院士[20]指出在煤层裂隙及大孔隙中的流动属于渗透, 符合达西定律; 而在微孔隙结构中的流动主要是扩散, 符合Fick 扩散定律; 并按扩散—渗透规律建立了煤层瓦斯流动方程式, 与按单纯达西定律导出的方程式进行了对比, 明确了达西定律的适用范围。
Saghafi 等[21]从扩散力学出发, 依据Fick 扩散定律建立了煤块瓦斯扩散方程; 又从渗流力学出发, 依据达西定律建立了瓦斯渗流方程; 最后耦合成瓦斯渗透-扩散的流动方程, 结合边值条件建立了瓦斯渗透-扩散的动力模型, 并进行了数值模拟。
吴世跃[22]借鉴石油天然气相关的渗流理论, 考虑了煤的结构特性与压缩性、吸附瓦斯与游离瓦斯流动的差异及其关系, 在均质扩散模型和渗透模型的基础上, 建立了扩散与渗流的物理数学模型, 该模型在一定条件下又可转化为纯扩散或纯渗透模型, 具有一定的普适性。
流动的幂定律[式(4) ]:
m
(4) v N =-A
d x
式中:v N 为标准状态下的瓦斯流速; A 为无量纲的瓦斯渗透率系数; m 为状态常数; d p /d x 为沿x 轴向的压力梯度。
该理论主要适用于中Re 区(雷诺数Re =10～100) , 此时流动区为非线性层流区域。
孙培德[23224]根据非均质多孔介质水动力学理论, , 建。尽, 但仍[20]。[, 提出了考虑Klinkenberg , 并建立了相应的瓦斯流动数学模型, 指出了达西定律的适用范围。
尽管关于幂定律是否优于线性达西定律存在不同的观点[26], 但研究和发展非线性瓦斯流动理论仍具有一定的意义。
3 地球物理场效应的渗流理论
一般将地层中的地应力场、地温场和地电场统称为地球物理场。随着对瓦斯流动机理研究的深化, 许多学者[27236]认识到了地球物理场等对瓦斯流动的影响; 基于对达西定律的各种修正, 建立和发展了相应的流动模型及数值方法。
文献[27231]是国外在该领域的代表性研究成果。Somerton [27]研究了裂隙煤体在三轴应力作用下氮气及甲烷气体的渗透性, 得出煤样的渗透性依赖于作用的应力, 且与应力史有关; 并指出随着地应力的增大, 煤层透气率按指数关系减小。文献[28231]从煤层赋存的地质条件和地球物理场环境出发, 开创性地考察和研究了含气煤样的力学性质, 以及瓦斯渗流和煤岩体之间的固气力学效应。
在国内, 刘保县等[32]首先研究了含瓦斯煤层所处的各地球物理场对瓦斯吸附特性的影响; 随后王宏图等[33234]根据地应力场、地温场以及地电场中瓦斯的渗流特性, 确定了煤层瓦斯渗透率k 与有效应力σ、温度T 及电场强度E 的关系k (σ, T , E ) , 建立了地球物理场中煤层瓦斯渗流方程, 该渗流方程可反映地应力场、地温场和地电场对煤层瓦斯渗流的共同影响, 为分析研究地球物理场作用下瓦斯的运移规律提供了一定的理论基础。易俊等[35]针对应力场与温度场、王宏图等[36]针对地电场作用分别完成了在相应地球物理场影响下的煤层甲烷气体渗流特性的实验研究。
2 非线性渗流理论
达西定律偏离的主要原因可归结为:①流量过大; ②分子效应; ③离子效应; ④流体非牛顿态势。在大量试验研究的基础上得到了更能符合煤层瓦斯
第1期 田智威:煤层气渗流理论及其研究进展 63
总的来说, 创建和发展地球物理场效应的瓦斯流动理论是国内外学者竞相研究的热点。用流体—岩石的相互作用去认识煤层内瓦斯运移的机制, 充分发展考虑地应力场、地温场以及地电场等地球物理场作用下的瓦斯流动模型及其数值方法, 尤其要注重发展可变形的孔隙—裂隙介质的多相耦合模型及其数值方法[37], 使物理模型更能反映客观事实, 进一步完善理论模型及测试技术。
4 多煤层瓦斯越流理论
在多层系统煤层孔隙压力梯度的驱动下, 邻近层的煤层气越过弱透气性的夹层向开采煤层的采场及采空区或采气钻孔迁移并涌出, 系统的煤层气越流场[38](图) 孙培德等[39], 根, 提出了双层系统越流的固气耦合数学模型, 为解决煤层气越流实际工程问题提供了一定的理论依据
。
测和控制提供了可行的研究方法。
根据孙培德[38,41]对越流场的定义, 关于煤层群开采中瓦斯涌出、保护层开采的有效保护范围的确定、井下邻近层瓦斯抽放的合理布孔设计、地面钻孔抽放多气层瓦斯工程的合理设计以及地下多气层之间瓦斯运移规律的预测和评估等问题, 均可归结为瓦斯越流问题。我国对于地下多气层之间瓦斯运移规律以及地面钻孔抽放多气层瓦斯工程的效果预估及合理设计等问题的研究甚少, 也未得到应有的重视, 。
1格子Boltzmann 方法
图1 多煤层瓦斯越流场示意图
Fig. 1 Sketch of coal gas leakage field in multiple coal seams 实际的越流大都是三维问题, 虽然在求解方法上并没有技术难度, 但其计算量与二维问题相比要大很多。孙培德等[39]根据实际情况, 对耦合模型进行了一定的简化处理:其中煤岩固体变形为平面应变模型, 煤层气流动为包含越流的一维平行流动模型, 这样既达到了数值求解的目的, 又大大简化了数学模型, 节省了大量的计算工作量。而在数值求解的算法方面, 孙培德[40]采用的是有限差分的强隐式法(SIP ) , 对耦合模型进行了离散化, 并作了详细的数值分析。随后, 孙培德等[41242]还应用三维模拟仿真技术实现了煤层气越流的三维可视化, 直观再现了煤层气越流场与煤岩体变形场耦合效应的宏观变化特征, 为地下煤层气越流场中邻近层瓦斯涌出预
格子Boltzmann 方法[43246](LBM ) 作为一种不同于传统数值方法的新兴流体计算和建模方法, 具有简单高效、适合大规模并行等优势, 自其诞生之日起就受到众多研究者的广泛关注, 除了在一般简单流动问题中得到了成功的应用外,LBM 还可应用于多相流、化学反应扩散、微尺度流、粒子悬浮流、多孔介质流及磁流体力学等领域。
滕桂荣等[47]采用格子Boltzmann 方法, 通过建立瓦斯渗流的动力学模型, 对瓦斯在裂隙煤体中的渗流运移规律进行了初步模拟, 研究了采动压力差对瓦斯流动速度及流动稳定性的影响规律, 模拟结果与用其他数值方法得到的瓦斯渗流规律比较吻合, 表明LBM 能够比较有效地模拟瓦斯在裂隙煤体中的运移规律, 为探讨煤与瓦斯两相耦合机制、瓦斯抽放方案的设计提供了新的思路。
滕桂荣等[48]又针对裂隙煤体内瓦斯的抽放, 采用LBM 对流场压力和流动速度的分布进行了模拟, 并针对不同的抽放压力与钻孔位置2个关键参数对瓦斯抽放效果的影响规律进行了分析, 为深入探讨瓦斯抽放机理提供了新的技术手段。
由于低渗透气藏的开发利用在我国有着极为重要的作用, 因此, 研究低渗透气藏中气体的渗流规律具有理论和应用价值。朱益华等[48]针对低渗透气藏中的Klinkenberg 效应(滑脱效应) , 利用LBM 易处理多孔介质孔隙复杂边界条件的优势, 很好地再现了滑脱效应, 不仅可以克服Klinkenberg 渗透率测量过程中存在的问题与困难, 而且可以对滑脱因子的影响因素进行定量分析, 为今后更深入的研究奠定了一定的基础。5. 2分子动力学方法
分子动力学方法[40251](MD ) 是一种基于微观领域的模拟方法, 通过模拟物质体系中微观粒子(原
64
地质科技情报
2010年
子、分子) 之间的相互作用及运动, 得到体系中粒子
的运动轨迹, 再按统计物理的方法计算得出宏观运动属性等。随着计算机技术的飞速发展, 该方法逐渐成为模拟微观系统、预测宏观性能的主要方法之一。
李希建等[52]首先从理论上分析了分子动力学模拟研究煤层瓦斯流动的可行性, 认为MD 模拟在流体微观性质方面的研究已经比较成熟, 尤其是在微通道中气体流动方面, 完全有可能根据MD 的模拟过程对煤层瓦斯流动的特性与规律进行模拟。接着, 他们又针对分子动力学模拟的各基本步骤进行了逐一分析, 指出了模拟中所需解决的几个关键问题:合适的势能函数的选择; 参量的确定; , , 是可以实现MD 模拟的, 。参考文献:
[1] Darcy H. L es Fontaines Publique de la V ille de Di j on [M ].
Paris :Victor Dalmont , 1856.
[2] 武法东, 王生维, 陈钟惠, 等. 新兴能源领域———煤层气研究开
重庆大学学报:自然科学版,1994,17(4) :58-63.
[18]Yu C , Xian X. A boundary element met hod for inhomoge 2
neous medium problems [C ]∥Anon. Proceedings of t he 2nd
W orl d Congress on Com p utational Mechanics .
Stuttgart :
FR G ,1990:308-311.
[19]杨其銮, 王佑安. 煤屑瓦斯扩散理论及其应用[J].煤炭学报,
1986, 11(3) :87-94.
[20]周世宁. 瓦斯在煤层中流动的机理[J].煤炭学报, 1990, 15
(1) :61-67.
[21]Saghafi A , J eger C , Tauziede C , et al. A new computer simu 2
lation of in seam gas flow and it s application to gas emission prediction and gas drainage [Dai Guoquan. Proceedings
of t nd I of S af et y in Mines Re 2search nstit Coal Industry Publishing
, -]. [J].山西矿业学
, 1994, 12(3) :259-263.
[23]孙培德. 煤层瓦斯流场流动规律的研究[J].煤炭学报, 1987,
12(4) :74-82.
[24]孙培德. 煤层瓦斯流动理论及其应用[C ]∥佚名. 中国煤炭学
会1988年学术年会论文集. 北京:煤炭工业出版社,1988:5.
[25]罗新荣. 煤层瓦斯运移物理模型与理论分析[J].中国矿业大学
学报, 1991, 20(3) :36-42.
[26]姚宇平. 煤层瓦斯流动的达西定律与幂定律[J].山西矿业学院
学报, 1992, 10(1) :32-37.
[27]Somerton W H. Effect of st ress on permeability of coal [J].
I nt. J. Rock Mech. Min. Sci.
&Geomech. A bst r . , 1975,
发的历史和现状[J].地质科技情报, 1992, 11(2) :69-74.
[3] 蔚远江, 杨起, 刘大锰, 等. 我国煤层气储层研究现状及发展趋
12(2) :151-158.
[28]Harpalaini S. Gas flow t hrough stressed coal [D ].Ph. D.
Thesis. Berkeley :Univ.of California Berkeley , 1985. [29]Harpalani S , Mopherson M J. The effect of gas evacatiom on
coal permeability test specimens [J ].I nt. J. Rock Mech.
Min. Sci. &Geomech. A bst r . , 1984, 21(3) :361-364.
势[J].地质科技情报, 2001, 20(1) :56-60.
[4] 梁冰, 孙可明. 低渗透煤层气开采理论及其应用[M ].北京:科
学出版社, 2006.
[5] 孙培德, 鲜学福. 煤层瓦斯渗流力学的研究进展[J].焦作工学
院学报, 2001, 20(3) :161-167.
[6] 李宏艳. 采动应力场与瓦斯渗流场耦合理论研究现状及趋势
[J].煤矿开采, 2008, 13(3) :4-7.
[7] 杨清岭, 杨鹏. 煤层瓦斯流动理论研究[J].煤炭技术, 2007, 26
(2) :67-69.
[8] 周世宁, 孙辑正. 煤层瓦斯流动理论及其应用[J].煤炭学报,
1965, 2(1) :24-36.
[9] 郭勇义, 周世宁. 煤层瓦斯一维流场流动规律的完全解[J].中
[30]Gawuga J. Folw of gas t hrough stressed carboniferous st rata
[D ].Ph. D. Thesis. Nottingham :Univ. of Nottingham , 1979.
[31]Khodot V V. Role of Met hane in t he stress state of a coal
seam[J].S oviet Mining S ci nce , 1981, 17(5) :460-466. [32]刘保县, 鲜学福, 徐龙君, 等. 地球物理场对煤吸附瓦斯特性的
影响[J].重庆大学学报, 2000, 23(5) :78-81.
[33]王宏图, 杜云贵, 鲜学福, 等. 地球物理场中的煤层瓦斯渗流方
国矿业学院学报, 1984, 13(2) :19-28.
[10]谭学术, 袁静. 矿井煤层真实瓦斯渗流方程的研究[J].重庆建
程[J].岩石力学与工程学报, 2002, 21(5) :644-646.
[34]王宏图, 杜云贵, 鲜学福, 等. 受地应力、地温和地电效应影响的
筑工程学院学报, 1986, 8(1) :106-112.
[11]余楚新, 鲜学福, 谭学术. 煤层瓦斯流动理论及渗流控制方程的
煤层瓦斯渗流方程[J].重庆大学学报, 2000, 23(增刊) :47-49.
[35]易俊, 姜永东, 鲜学福. 应力场、温度场瓦斯渗流特性实验研究
[J].中国矿业, 2007, 16(5) :113-116.
[36]王宏图, 李晓红, 鲜学福, 等. 地电场作用下煤中甲烷气体渗流
研究[J].重庆大学学报, 1989, 12(5) :1-10.
[12]Sun P. Coal gas dynamics and it s applications [J].Scientia
Geologica S inica , 1994, 3(1) :83-90.
[13]孙培德. 瓦斯动力学模型的研究[J].煤田地质与勘探, 1993,
21(1) :32-40.
[14]孙培德. 煤层瓦斯流动方程补正[J].煤田地质与勘探, 1993,
21(5) :34-35.
[15]李英俊, 赵均. 煤层瓦斯压力分布的研究[J].煤矿安全, 1980
(5) :6-11.
[16]魏晓林. 煤层瓦斯流动规律的实验和数值方法的研究[J].粤煤
性质的实验研究[J].岩石力学与工程学报, 2004, 23(2) :303
-306.
[37]赵阳升, 段康廉, 胡耀青, 等. 块裂介质岩石流体力学研究新进
展[J].辽宁工程技术大学学报:自然科学版,1999,18(5) :459
-462.
[38]孙培德. 煤层气越流的固气耦合理论及其计算机模拟研究
[D ].重庆:重庆大学博士学位论文, 1998.
[39]孙培德, 鲜学福. 煤层气越流的固气耦合理论及其应用[J].煤
科技, 1981, (2) :35-41.
[17]余楚新, 鲜学福. 煤层瓦斯渗流有限元分析中的几个问题[J].
炭学报, 1999, 24(1) :60-64.
第1期 田智威:煤层气渗流理论及其研究进展 65
[40]孙培德. 煤层气越流固气耦合数学模型的SIP 分析[J].煤炭学
[47]滕桂荣, 谭云亮, 高明. 基于Lattice Boltzmann 方法对裂隙煤
报, 2002, 27(5) :494-498.
[41]孙培德. SUN 模型及其应用———煤层气越流固气耦合模型及
体中瓦斯运移规律的模拟研究[J].岩石力学与工程学报,
2007, 26(增刊1) :3503-3508.
[48]滕桂荣, 谭云亮, 高明, 等. 基于LBM 方法的裂隙煤体内瓦斯
可视化模拟[M ].杭州:浙江大学出版社, 2002.
[42]孙培德, 万华根. 煤层气越流固-气耦合模型及可视化模拟研
抽放的模拟分析[J].煤炭学报, 2008, 33(8) :914-919.
[49]朱益华, 陶果, 方伟, 等. 低渗气藏中气体渗流Klinkenberg 效
究[J].岩石力学与工程学报, 2004, 23(7) :1179-1185.
[43]Chen S , Doolen G D. Lattice Boltzmann met hod for fluid flows
[J].A nnu. Rev. Fl ui d Mech. , 1998, 30:329-364. [44]郭照立, 郑楚光, 李青, 等. 流体动力学的格子Boltzmann 方法
[M ].武汉:湖北科学技术出版社, 2002.
[45]Succi S. T he L attice B oltz mann equation f or f l ui d d y namics
and beyond [M ].Oxford :Clarendon Press , 2001.
应研究进展[J ].地球物理学进展, 2007, 22(5) :1591-1596.
[50]Allen M P , Tildesley D J. Com puter si mul ation of liqui ds
[M ].Oxford :Clarendon Press , 1994.
[51]K oplik J , Baanavar J R. Continuum deductions from molecular
hydrodynamics[J].A nnu. ui d Mech. , 1995,27:257-292.
[李希建[46]Sukop M C , Thorne D T. L attice B oltz mann modeling :A n
int roduction f or geoscientists and engineers [M ].Berlin , Hei 2
delberg ,New Y ork :Springer , 2006. [J, 2008, 39(7) :84-87.
R Methane Seepage Mechanics
TIAN Zhi 2Wei
(Dep artment of Coal &Coalbed Met hane Engi neeri ng , Facult y of Eart h Resources ,
Chi na U ni versit y of Geosciences , W uhan 430074, Chi na )
Abstract :This paper present s an overview of t he research progress of coal gas seepage mechanics , bot h in China and abroad. It first int roduces some elementary seepage laws , including linear infiltration , linear dif 2f usion , diff usion 2infilt ration and nonlinear seepage t heories. Also it analyzes t he flow rules concerning t he effect s of different geop hysical fields and multiple coal seam leakflows. Finally , several novel numerical met hods are mentioned in t he application of coalbed met hane seepage flows.
K ey w ords :coalbed met hane ; seepage mechanics ; numerical simulation ; research advance
第29卷 第1期地质科技情报Vol. 29 No. 1
2010年 1月Geological Science and Technology Information Jan. 2010
煤层气渗流理论及其研究进展
田智威
(中国地质大学资源学院煤及煤层气工程系, 武汉430074)
3
摘 要:综述了煤层气渗流研究的现状及取得的研究进展, 包括线性渗透、线性扩散、扩散—渗透以及非线性渗流等基本理论, 基于各种地球物理场效应影响的渗流理论及其修正模型, 以及多煤层瓦斯越流相关的研究; (格子Boltzmann 方法与分子动力学方法) 在煤层瓦斯渗流模拟中的应用; 关键词:煤层气; 渗流; 数值模拟; 研究进展中图分类号:P618.11 文献标志码:A2) 0120061205 应用的科学。自]律以来, , 并逐渐与其他学科交叉, 在能源资源的开发与利用以及工程建设中得到了非常广泛的应用。渗流力学最先应用在水利工程和地下水资源开发等领域; 随后又成为石油和天然气工业的一项基础理论。
随着煤层气这一新型清洁能源的重视与开发[223], 渗流理论又应用在煤层气的开发与利用中[326]。煤层气渗流力学是研究煤层内瓦斯压力分布及其流动变化规律的理论[427], 是由渗流力学、煤地质学、固体力学及采矿学等学科互相交叉渗透发展形成的。自煤层气渗流力学创立至今深受有关研究人员的关注, 尤其自20世纪80年代以来发展更为迅速, 表现在:应用范围更广; 基本理论不断深化; 研究手段及方法不断现代化。
式中:v 为流速; k 为煤层的渗透率; μ为瓦斯黏度系数; p 为瓦斯压力; 5p/5n 为瓦斯压力在流动方向上的偏导数。
此外还导出了瓦斯流量方程[式(2) ]:
(2) q =-λn
式中:q 为瓦斯流量; λ为煤层透气系数。
20世纪80年代, 多位研究者[9214]在修正和完善数学模型、流动方程方面开展了相应的工作。由于大多数井下瓦斯流动都可简化为一维的平行流动和径向流动的有限流场、无限流场或其组合, 为此, 郭勇义等[9]针对一维流动, 结合相似理论提出了修正的流动方程。谭学术等[10]利用渗流理论和热力学理论, 提出了应用真实气体状态方程(而不是理想气体状态方程) 的煤层气渗流方程, 并做了线性化简化处理。余楚新等[11]提出参与渗流的瓦斯量只是可解吸的那一部分, 在瓦斯吸附与解吸过程完全可逆的假设条件下建立了相应的控制方程, 做出了有一定新意的探索。随后, 孙培德[12214]又进一步修正和完善了线性流动数学模型。
随着电子计算机和数值模拟方法的发展, 针对煤层气流场内的压力分布及其流动规律进行数值模拟已经成为研究的主要手段。针对压力是煤层气流动的主要动力因素, 李英俊等[15]重点模拟研究了煤层瓦斯压力的分布规律, 为预测瓦斯涌出量及抽放效果提供了参考依据; 魏晓林[16]采用有限差分法对瓦斯流场中的压力分布以及流量变化进行了数值模拟; Yu 等[17218]分别采用有限元法(FEM ) 和边界单元法(B EM ) 实现了对瓦斯渗流流动的数值模拟。
1 线性渗流理论
1. 1线性渗透理论
为了适应采矿采煤业的大力发展, 控制瓦斯技术已成为当时研究的关键技术之一, 早在20世纪40年代末, 前苏联学者就已经建立起考虑吸附煤层瓦斯作用的瓦斯控制方程[式(1) ]。在我国, 周世宁院士等[8]首先进行了将达西定律应用于煤层瓦斯流动理论的开拓性研究, 认为煤层瓦斯的流动基本符合线性渗流规律, 其观点对煤层瓦斯渗流的应用和瓦斯动力学研究具有相当重要的指导意义。
(1) v =-μn
3收稿日期:2009202216 编辑:禹华珍
基金项目:中国地质大学(武汉) 优秀青年教师资助计划项目(CU GQNL0904) ; 煤燃烧国家重点实验室开放基金项目(FSKLCC0815)
作者简介:田智威(1980—) , 男, 主要从事煤层气渗流理论与数值模拟的教学与科研工作。
62
地质科技情报
2010年
1. 2线性扩散理论
瓦斯扩散理论在欧美国家有较多研究, 而在我
国这方面的研究较少[4], 其扩散规律的理论依据是Fick 扩散定律[式(3) ]:
(3) J =-D
n
式中:J 为单位面积的扩散速度; D 为瓦斯扩散系数; c 为煤体瓦斯浓度。
然而瓦斯的涌出过程是一个很复杂的过程, 杨其銮等[19]将煤层看作微孔隙裂隙体来讨论瓦斯的扩散过程, 认为其属于气体在多孔介质中的扩散, 其涌出规律符合扩散定律。依据扩散运动与热传导类似的原理, 可以借助热传导的相关研究方法来讨论瓦斯的扩散问题, 探讨和实测对比分析。1. 3扩散—瓦斯扩散—, 煤层内瓦斯的运动是包含了渗透与扩散的混合流动过程。煤层中存在相互沟通的裂隙网络, 游离瓦斯沿着这些裂隙流向低压工作面; 同时煤块内部的瓦斯解吸附而向裂隙扩散, 因此, 煤层中瓦斯的渗透率及其扩散性共同决定了其流动状况。
随着瓦斯运移规律研究的深入, 国内外许多学者都赞同这一理论, 周世宁院士[20]指出在煤层裂隙及大孔隙中的流动属于渗透, 符合达西定律; 而在微孔隙结构中的流动主要是扩散, 符合Fick 扩散定律; 并按扩散—渗透规律建立了煤层瓦斯流动方程式, 与按单纯达西定律导出的方程式进行了对比, 明确了达西定律的适用范围。
Saghafi 等[21]从扩散力学出发, 依据Fick 扩散定律建立了煤块瓦斯扩散方程; 又从渗流力学出发, 依据达西定律建立了瓦斯渗流方程; 最后耦合成瓦斯渗透-扩散的流动方程, 结合边值条件建立了瓦斯渗透-扩散的动力模型, 并进行了数值模拟。
吴世跃[22]借鉴石油天然气相关的渗流理论, 考虑了煤的结构特性与压缩性、吸附瓦斯与游离瓦斯流动的差异及其关系, 在均质扩散模型和渗透模型的基础上, 建立了扩散与渗流的物理数学模型, 该模型在一定条件下又可转化为纯扩散或纯渗透模型, 具有一定的普适性。
流动的幂定律[式(4) ]:
m
(4) v N =-A
d x
式中:v N 为标准状态下的瓦斯流速; A 为无量纲的瓦斯渗透率系数; m 为状态常数; d p /d x 为沿x 轴向的压力梯度。
该理论主要适用于中Re 区(雷诺数Re =10～100) , 此时流动区为非线性层流区域。
孙培德[23224]根据非均质多孔介质水动力学理论, , 建。尽, 但仍[20]。[, 提出了考虑Klinkenberg , 并建立了相应的瓦斯流动数学模型, 指出了达西定律的适用范围。
尽管关于幂定律是否优于线性达西定律存在不同的观点[26], 但研究和发展非线性瓦斯流动理论仍具有一定的意义。
3 地球物理场效应的渗流理论
一般将地层中的地应力场、地温场和地电场统称为地球物理场。随着对瓦斯流动机理研究的深化, 许多学者[27236]认识到了地球物理场等对瓦斯流动的影响; 基于对达西定律的各种修正, 建立和发展了相应的流动模型及数值方法。
文献[27231]是国外在该领域的代表性研究成果。Somerton [27]研究了裂隙煤体在三轴应力作用下氮气及甲烷气体的渗透性, 得出煤样的渗透性依赖于作用的应力, 且与应力史有关; 并指出随着地应力的增大, 煤层透气率按指数关系减小。文献[28231]从煤层赋存的地质条件和地球物理场环境出发, 开创性地考察和研究了含气煤样的力学性质, 以及瓦斯渗流和煤岩体之间的固气力学效应。
在国内, 刘保县等[32]首先研究了含瓦斯煤层所处的各地球物理场对瓦斯吸附特性的影响; 随后王宏图等[33234]根据地应力场、地温场以及地电场中瓦斯的渗流特性, 确定了煤层瓦斯渗透率k 与有效应力σ、温度T 及电场强度E 的关系k (σ, T , E ) , 建立了地球物理场中煤层瓦斯渗流方程, 该渗流方程可反映地应力场、地温场和地电场对煤层瓦斯渗流的共同影响, 为分析研究地球物理场作用下瓦斯的运移规律提供了一定的理论基础。易俊等[35]针对应力场与温度场、王宏图等[36]针对地电场作用分别完成了在相应地球物理场影响下的煤层甲烷气体渗流特性的实验研究。
2 非线性渗流理论
达西定律偏离的主要原因可归结为:①流量过大; ②分子效应; ③离子效应; ④流体非牛顿态势。在大量试验研究的基础上得到了更能符合煤层瓦斯
第1期 田智威:煤层气渗流理论及其研究进展 63
总的来说, 创建和发展地球物理场效应的瓦斯流动理论是国内外学者竞相研究的热点。用流体—岩石的相互作用去认识煤层内瓦斯运移的机制, 充分发展考虑地应力场、地温场以及地电场等地球物理场作用下的瓦斯流动模型及其数值方法, 尤其要注重发展可变形的孔隙—裂隙介质的多相耦合模型及其数值方法[37], 使物理模型更能反映客观事实, 进一步完善理论模型及测试技术。
4 多煤层瓦斯越流理论
在多层系统煤层孔隙压力梯度的驱动下, 邻近层的煤层气越过弱透气性的夹层向开采煤层的采场及采空区或采气钻孔迁移并涌出, 系统的煤层气越流场[38](图) 孙培德等[39], 根, 提出了双层系统越流的固气耦合数学模型, 为解决煤层气越流实际工程问题提供了一定的理论依据
。
测和控制提供了可行的研究方法。
根据孙培德[38,41]对越流场的定义, 关于煤层群开采中瓦斯涌出、保护层开采的有效保护范围的确定、井下邻近层瓦斯抽放的合理布孔设计、地面钻孔抽放多气层瓦斯工程的合理设计以及地下多气层之间瓦斯运移规律的预测和评估等问题, 均可归结为瓦斯越流问题。我国对于地下多气层之间瓦斯运移规律以及地面钻孔抽放多气层瓦斯工程的效果预估及合理设计等问题的研究甚少, 也未得到应有的重视, 。
1格子Boltzmann 方法
图1 多煤层瓦斯越流场示意图
Fig. 1 Sketch of coal gas leakage field in multiple coal seams 实际的越流大都是三维问题, 虽然在求解方法上并没有技术难度, 但其计算量与二维问题相比要大很多。孙培德等[39]根据实际情况, 对耦合模型进行了一定的简化处理:其中煤岩固体变形为平面应变模型, 煤层气流动为包含越流的一维平行流动模型, 这样既达到了数值求解的目的, 又大大简化了数学模型, 节省了大量的计算工作量。而在数值求解的算法方面, 孙培德[40]采用的是有限差分的强隐式法(SIP ) , 对耦合模型进行了离散化, 并作了详细的数值分析。随后, 孙培德等[41242]还应用三维模拟仿真技术实现了煤层气越流的三维可视化, 直观再现了煤层气越流场与煤岩体变形场耦合效应的宏观变化特征, 为地下煤层气越流场中邻近层瓦斯涌出预
格子Boltzmann 方法[43246](LBM ) 作为一种不同于传统数值方法的新兴流体计算和建模方法, 具有简单高效、适合大规模并行等优势, 自其诞生之日起就受到众多研究者的广泛关注, 除了在一般简单流动问题中得到了成功的应用外,LBM 还可应用于多相流、化学反应扩散、微尺度流、粒子悬浮流、多孔介质流及磁流体力学等领域。
滕桂荣等[47]采用格子Boltzmann 方法, 通过建立瓦斯渗流的动力学模型, 对瓦斯在裂隙煤体中的渗流运移规律进行了初步模拟, 研究了采动压力差对瓦斯流动速度及流动稳定性的影响规律, 模拟结果与用其他数值方法得到的瓦斯渗流规律比较吻合, 表明LBM 能够比较有效地模拟瓦斯在裂隙煤体中的运移规律, 为探讨煤与瓦斯两相耦合机制、瓦斯抽放方案的设计提供了新的思路。
滕桂荣等[48]又针对裂隙煤体内瓦斯的抽放, 采用LBM 对流场压力和流动速度的分布进行了模拟, 并针对不同的抽放压力与钻孔位置2个关键参数对瓦斯抽放效果的影响规律进行了分析, 为深入探讨瓦斯抽放机理提供了新的技术手段。
由于低渗透气藏的开发利用在我国有着极为重要的作用, 因此, 研究低渗透气藏中气体的渗流规律具有理论和应用价值。朱益华等[48]针对低渗透气藏中的Klinkenberg 效应(滑脱效应) , 利用LBM 易处理多孔介质孔隙复杂边界条件的优势, 很好地再现了滑脱效应, 不仅可以克服Klinkenberg 渗透率测量过程中存在的问题与困难, 而且可以对滑脱因子的影响因素进行定量分析, 为今后更深入的研究奠定了一定的基础。5. 2分子动力学方法
分子动力学方法[40251](MD ) 是一种基于微观领域的模拟方法, 通过模拟物质体系中微观粒子(原
64
地质科技情报
2010年
子、分子) 之间的相互作用及运动, 得到体系中粒子
的运动轨迹, 再按统计物理的方法计算得出宏观运动属性等。随着计算机技术的飞速发展, 该方法逐渐成为模拟微观系统、预测宏观性能的主要方法之一。
李希建等[52]首先从理论上分析了分子动力学模拟研究煤层瓦斯流动的可行性, 认为MD 模拟在流体微观性质方面的研究已经比较成熟, 尤其是在微通道中气体流动方面, 完全有可能根据MD 的模拟过程对煤层瓦斯流动的特性与规律进行模拟。接着, 他们又针对分子动力学模拟的各基本步骤进行了逐一分析, 指出了模拟中所需解决的几个关键问题:合适的势能函数的选择; 参量的确定; , , 是可以实现MD 模拟的, 。参考文献:
[1] Darcy H. L es Fontaines Publique de la V ille de Di j on [M ].
Paris :Victor Dalmont , 1856.
[2] 武法东, 王生维, 陈钟惠, 等. 新兴能源领域———煤层气研究开
重庆大学学报:自然科学版,1994,17(4) :58-63.
[18]Yu C , Xian X. A boundary element met hod for inhomoge 2
neous medium problems [C ]∥Anon. Proceedings of t he 2nd
W orl d Congress on Com p utational Mechanics .
Stuttgart :
FR G ,1990:308-311.
[19]杨其銮, 王佑安. 煤屑瓦斯扩散理论及其应用[J].煤炭学报,
1986, 11(3) :87-94.
[20]周世宁. 瓦斯在煤层中流动的机理[J].煤炭学报, 1990, 15
(1) :61-67.
[21]Saghafi A , J eger C , Tauziede C , et al. A new computer simu 2
lation of in seam gas flow and it s application to gas emission prediction and gas drainage [Dai Guoquan. Proceedings
of t nd I of S af et y in Mines Re 2search nstit Coal Industry Publishing
, -]. [J].山西矿业学
, 1994, 12(3) :259-263.
[23]孙培德. 煤层瓦斯流场流动规律的研究[J].煤炭学报, 1987,
12(4) :74-82.
[24]孙培德. 煤层瓦斯流动理论及其应用[C ]∥佚名. 中国煤炭学
会1988年学术年会论文集. 北京:煤炭工业出版社,1988:5.
[25]罗新荣. 煤层瓦斯运移物理模型与理论分析[J].中国矿业大学
学报, 1991, 20(3) :36-42.
[26]姚宇平. 煤层瓦斯流动的达西定律与幂定律[J].山西矿业学院
学报, 1992, 10(1) :32-37.
[27]Somerton W H. Effect of st ress on permeability of coal [J].
I nt. J. Rock Mech. Min. Sci.
&Geomech. A bst r . , 1975,
发的历史和现状[J].地质科技情报, 1992, 11(2) :69-74.
[3] 蔚远江, 杨起, 刘大锰, 等. 我国煤层气储层研究现状及发展趋
12(2) :151-158.
[28]Harpalaini S. Gas flow t hrough stressed coal [D ].Ph. D.
Thesis. Berkeley :Univ.of California Berkeley , 1985. [29]Harpalani S , Mopherson M J. The effect of gas evacatiom on
coal permeability test specimens [J ].I nt. J. Rock Mech.
Min. Sci. &Geomech. A bst r . , 1984, 21(3) :361-364.
势[J].地质科技情报, 2001, 20(1) :56-60.
[4] 梁冰, 孙可明. 低渗透煤层气开采理论及其应用[M ].北京:科
学出版社, 2006.
[5] 孙培德, 鲜学福. 煤层瓦斯渗流力学的研究进展[J].焦作工学
院学报, 2001, 20(3) :161-167.
[6] 李宏艳. 采动应力场与瓦斯渗流场耦合理论研究现状及趋势
[J].煤矿开采, 2008, 13(3) :4-7.
[7] 杨清岭, 杨鹏. 煤层瓦斯流动理论研究[J].煤炭技术, 2007, 26
(2) :67-69.
[8] 周世宁, 孙辑正. 煤层瓦斯流动理论及其应用[J].煤炭学报,
1965, 2(1) :24-36.
[9] 郭勇义, 周世宁. 煤层瓦斯一维流场流动规律的完全解[J].中
[30]Gawuga J. Folw of gas t hrough stressed carboniferous st rata
[D ].Ph. D. Thesis. Nottingham :Univ. of Nottingham , 1979.
[31]Khodot V V. Role of Met hane in t he stress state of a coal
seam[J].S oviet Mining S ci nce , 1981, 17(5) :460-466. [32]刘保县, 鲜学福, 徐龙君, 等. 地球物理场对煤吸附瓦斯特性的
影响[J].重庆大学学报, 2000, 23(5) :78-81.
[33]王宏图, 杜云贵, 鲜学福, 等. 地球物理场中的煤层瓦斯渗流方
国矿业学院学报, 1984, 13(2) :19-28.
[10]谭学术, 袁静. 矿井煤层真实瓦斯渗流方程的研究[J].重庆建
程[J].岩石力学与工程学报, 2002, 21(5) :644-646.
[34]王宏图, 杜云贵, 鲜学福, 等. 受地应力、地温和地电效应影响的
筑工程学院学报, 1986, 8(1) :106-112.
[11]余楚新, 鲜学福, 谭学术. 煤层瓦斯流动理论及渗流控制方程的
煤层瓦斯渗流方程[J].重庆大学学报, 2000, 23(增刊) :47-49.
[35]易俊, 姜永东, 鲜学福. 应力场、温度场瓦斯渗流特性实验研究
[J].中国矿业, 2007, 16(5) :113-116.
[36]王宏图, 李晓红, 鲜学福, 等. 地电场作用下煤中甲烷气体渗流
研究[J].重庆大学学报, 1989, 12(5) :1-10.
[12]Sun P. Coal gas dynamics and it s applications [J].Scientia
Geologica S inica , 1994, 3(1) :83-90.
[13]孙培德. 瓦斯动力学模型的研究[J].煤田地质与勘探, 1993,
21(1) :32-40.
[14]孙培德. 煤层瓦斯流动方程补正[J].煤田地质与勘探, 1993,
21(5) :34-35.
[15]李英俊, 赵均. 煤层瓦斯压力分布的研究[J].煤矿安全, 1980
(5) :6-11.
[16]魏晓林. 煤层瓦斯流动规律的实验和数值方法的研究[J].粤煤
性质的实验研究[J].岩石力学与工程学报, 2004, 23(2) :303
-306.
[37]赵阳升, 段康廉, 胡耀青, 等. 块裂介质岩石流体力学研究新进
展[J].辽宁工程技术大学学报:自然科学版,1999,18(5) :459
-462.
[38]孙培德. 煤层气越流的固气耦合理论及其计算机模拟研究
[D ].重庆:重庆大学博士学位论文, 1998.
[39]孙培德, 鲜学福. 煤层气越流的固气耦合理论及其应用[J].煤
科技, 1981, (2) :35-41.
[17]余楚新, 鲜学福. 煤层瓦斯渗流有限元分析中的几个问题[J].
炭学报, 1999, 24(1) :60-64.
第1期 田智威:煤层气渗流理论及其研究进展 65
[40]孙培德. 煤层气越流固气耦合数学模型的SIP 分析[J].煤炭学
[47]滕桂荣, 谭云亮, 高明. 基于Lattice Boltzmann 方法对裂隙煤
报, 2002, 27(5) :494-498.
[41]孙培德. SUN 模型及其应用———煤层气越流固气耦合模型及
体中瓦斯运移规律的模拟研究[J].岩石力学与工程学报,
2007, 26(增刊1) :3503-3508.
[48]滕桂荣, 谭云亮, 高明, 等. 基于LBM 方法的裂隙煤体内瓦斯
可视化模拟[M ].杭州:浙江大学出版社, 2002.
[42]孙培德, 万华根. 煤层气越流固-气耦合模型及可视化模拟研
抽放的模拟分析[J].煤炭学报, 2008, 33(8) :914-919.
[49]朱益华, 陶果, 方伟, 等. 低渗气藏中气体渗流Klinkenberg 效
究[J].岩石力学与工程学报, 2004, 23(7) :1179-1185.
[43]Chen S , Doolen G D. Lattice Boltzmann met hod for fluid flows
[J].A nnu. Rev. Fl ui d Mech. , 1998, 30:329-364. [44]郭照立, 郑楚光, 李青, 等. 流体动力学的格子Boltzmann 方法
[M ].武汉:湖北科学技术出版社, 2002.
[45]Succi S. T he L attice B oltz mann equation f or f l ui d d y namics
and beyond [M ].Oxford :Clarendon Press , 2001.
应研究进展[J ].地球物理学进展, 2007, 22(5) :1591-1596.
[50]Allen M P , Tildesley D J. Com puter si mul ation of liqui ds
[M ].Oxford :Clarendon Press , 1994.
[51]K oplik J , Baanavar J R. Continuum deductions from molecular
hydrodynamics[J].A nnu. ui d Mech. , 1995,27:257-292.
[李希建[46]Sukop M C , Thorne D T. L attice B oltz mann modeling :A n
int roduction f or geoscientists and engineers [M ].Berlin , Hei 2
delberg ,New Y ork :Springer , 2006. [J, 2008, 39(7) :84-87.
R Methane Seepage Mechanics
TIAN Zhi 2Wei
(Dep artment of Coal &Coalbed Met hane Engi neeri ng , Facult y of Eart h Resources ,
Chi na U ni versit y of Geosciences , W uhan 430074, Chi na )
Abstract :This paper present s an overview of t he research progress of coal gas seepage mechanics , bot h in China and abroad. It first int roduces some elementary seepage laws , including linear infiltration , linear dif 2f usion , diff usion 2infilt ration and nonlinear seepage t heories. Also it analyzes t he flow rules concerning t he effect s of different geop hysical fields and multiple coal seam leakflows. Finally , several novel numerical met hods are mentioned in t he application of coalbed met hane seepage flows.
K ey w ords :coalbed met hane ; seepage mechanics ; numerical simulation ; research advance