考虑滑脱效应的煤层气渗流数学模型及数值模拟
1煤层气藏三维数值模拟_冯文光
卷(Volume)19,期(Number)1,总(Total)75矿物岩石 页(Pages )43-48,1999,3,(M ar ,1999)J M INE RAL PET ROL 煤层气藏三维数值模拟冯文光1 梅世昕2 侯鸿斌2 羊裔常1(1 成都理工学院石油系,成都 610059)(2 华北石油地质局,郑州 450007)【摘 要】 本文根据煤的储气特征、煤层气的吸附特征、排水降压解吸产气机理建立了煤层气藏数学模型和三维全隐式差分模型,采用最佳变松驰求解,收敛度快,稳定性好。
【关键词】 煤层气 解吸 数值模拟 最佳变松驰法 全隐式中图法分类号:T E 319ISSN 1001-6872(1999)01-0043-48; C ODEN:KUYAE2收稿日期(M anuscript received):1998-12-10 改回日期(Accepted for publication):1999-02-05第一作者简介:冯文光 男 51岁 教授(博士导师) 油藏工程专业 研究方向:油气田开发、数值模拟,试井分析及煤层气藏0 引 言在地下煤层形成过程中,生成大量的天然气,一部分以游离状态通过运移、聚集形成常规天然气藏;另一部分以吸附状态为主储存在煤层的孔隙中,这种储存在煤层气中以甲烷为主的天然气称煤层气,或煤层瓦斯。
60年代和70年代初,美国矿业局为改善煤矿安全,对煤层气做了大量的研究工作。
70年代末,煤层气的开发利用还没有被引起重视,甚至把以煤层气作为商品的少数企业家视为“疯了”。
美国1977年开采煤气以来,开采井数成倍增加,年产煤层气产量倍增。
1989年美国已有400口井生产,年产量逾20×108m 3。
1992年美国煤层气年产量达200×108m 3。
我国煤层气蕴藏量十分丰富,与常规天然气蕴藏量相近,华北石油地质局自1986年以来,进行了多项攻关研究。
然而煤层气这一巨大而洁净的资源在我国至今尚未被开发和利用。
大倾角煤层浅部露出区煤层气溢出问题数值仿真模拟
大倾角煤层浅部露出区煤层气溢出问题数值仿真模拟王洪利;张遂安;黄红星;张潇【期刊名称】《煤炭科学技术》【年(卷),期】2022(50)10【摘要】中国新疆阜康白杨河矿区42号煤层倾角大,浅部煤系地层露出,为研究煤层气从该区域溢出问题,防止煤层气资源浪费,基于大倾角煤层中气-水分异现象,修正了大倾角煤层三维气-水两相渗流模型,结合阜康白杨河矿区42号煤层实际地质资料与历史拟合储层再描述,展开了大倾角煤层浅部露出区煤层气溢出问题数值模拟仿真研究。
研究结果表明:大倾角煤层浅部露出区煤层气溢出问题的存在;埋深300、500和700 m直井生产15 a累计溢出量分别为1.02×10^(5)、1.08×10^(7)以及1.33×10^(7) m^(3),和产量比值分别是5.06%、668.02%以及335.77%,直井生产时会造成大量的煤层气资源浪费;溢出现象并非伴随整个生产周期,溢出量先增后降直到溢出现象结束,使用顺层井生产时可以完整观测这一过程,溢出现象持续了1918 d,日溢出量上升和下降过程近似对称,累计溢出量为7.18×10^(5) m^(3),和累计产气量比值为1.54%。
补充400 m和600 m埋深直井模拟方案后,发现煤层气最大溢出量和累计溢出量与直井埋深呈正相关线性关系,而溢出现象开始时间随着埋深增加逐渐提前,500 m以深溢出现象开始时间保持不变。
综合分析煤层气溢出现象是造成该区直井产能较差的原因之一,使用顺层井排采能有效提高经济效益和降低煤层气溢出现象。
【总页数】8页(P143-150)【作者】王洪利;张遂安;黄红星;张潇【作者单位】中国石油大学(北京)石油工程学院;中联煤层气国家工程研究中心有限责任公司;煤层气开发利用国家工程研究中心;中国石油大学(北京)理学院【正文语种】中文【中图分类】TD713【相关文献】1.煤层气羽状水平井数值模拟技术在大宁煤层气试验区的应用2.基于Matlab的大倾角煤层飞矸冲击后果数值模拟3.大倾角煤层开采覆岩移动及应力变化数值模拟分析4.煤层倾角对煤层气水平井产能影响的数值模拟5.厚煤层大倾角煤顶下掘进围岩控制数值模拟研究因版权原因，仅展示原文概要，查看原文内容请购买。
温度作用下煤层瓦斯解吸渗流规律数值模拟
温度作用下煤层瓦斯解吸渗流规律数值模拟
张永利;张乐乐;马玉林;王金铜
【期刊名称】《防灾减灾工程学报》
【年(卷),期】2014(34)6
【摘要】如何提高煤层气渗透率是目前煤层气开采研究中的重要课题。
基于煤层瓦斯渗流规律数学模型,利用COMSOL Multiphysics软件,对流-固-热耦合条件下的非等温煤层气解吸、渗流变化规律进行了数值模拟。
结果表明,在注热条件下,煤层气渗流压力随着温度的增加而下降,且下降速度加剧,压力差越大,气体从高压区域流向低压区域的渗流速度越快。
气体在煤层中径向流向井口,井口附近压力的梯度增大,气体渗流速度较快;在未受到加热影响的区域,煤层气不受外加热量影响,煤层气解吸速率保持不变;注热后煤层温度升高,可以加快煤层气渗流速度、提高渗透率、增加煤层气产量。
研究成果可为煤层中注热开采煤层气的工程实践提供相应的理论依据。
【总页数】7页(P671-677)
【作者】张永利;张乐乐;马玉林;王金铜
【作者单位】辽宁工程技术大学力学与工程学院
【正文语种】中文
【中图分类】TD8
【相关文献】
1.超声热效应促进煤层瓦斯解吸扩散数值模拟
2.厚关键层下突出煤层煤岩瓦斯卸压规律数值模拟研究
3.变温条件下深部煤层瓦斯分布规律的数值模拟
4.粒径-温度耦合作用下煤中瓦斯解吸规律试验研究
5.倾斜高瓦斯煤层抽采条件下采空区漏风规律数值模拟
因版权原因，仅展示原文概要，查看原文内容请购买。
考虑Klinkenberg效应的多孔介质气体渗流模型
抽采效率明显提高。
为了从本质上更加清楚地认识 Ｋｌｉｎｋｅｎｂｅｒｇ效
应对气体渗流的影响，在前人研究［９１１］的基础上，本 文基于渗流力 学 理 论［１２］，推 导 出 了 考 虑 气 体 可 压 缩
并将其应用 于 单 一 钻 孔 的 瓦 斯 抽 采 数 值 模 型，对 比
分析了 Ｋｌｉｎｋｅｎｂｅｒｇ效应对瓦斯压力、渗透率和孔隙 率动态变 化 的 影 响 规 律；刘 清 泉 等［８］建 立 了 煤 体 孔
隙裂隙二重介质 ＰＭ 渗透模型，模拟分析了深部首 采层顺层钻孔预抽消突过程中煤层瓦斯压力及渗透
第１３卷 第３期 ２０１８年２月
中国科技论文 ＣＨＩＮＡＳＣＩＥＮＣＥＰＡＰＥＲ
Ｖｏｌ．１３Ｎｏ．３ Ｆｅｂ．２０１８
考虑 犓犾犻狀犽犲狀犫犲狉犵效应的多孔介质气体渗流模型
茹忠亮１，２，简 阔１，２，马国胜１
（１．山西省晋城市质量技术监督检验测试所，山西晋城 ０４８０００；２．太原科技大学资源与能源学院，山西晋城 ０４８０００）
多孔介 质 中 气 体 渗 流 机 理 是 天 然 气、煤 层 气 抽 采工程中最基础的科学问题。页岩气、煤层气、致密
砂岩气等非常规天然气的储层均为致密多孔介质，
以微米及纳米级孔隙为主，储层渗透率低著［１３］。Ｋｌｉｎｋｅｎｂｅｒｇ 效应 也 称 为 气 体 渗 流 滑 脱 效 应，１９４１ 年 Ｋｌｉｎｋｅｎ ｂｅｒｇ［４］提出多孔介质中气体的滑脱效应会影响其渗 透率，并给出了渗透率与气体压力之间的关系，很多 学者对此现象进行了研究。刘圣鑫等［５］通过计算分
性及 Ｋｌｉｎｋｅｎｂｅｒｇ效应的多孔介质气体渗流数学模 型，给出了压 力 边 界 条 件 和 流 量 边 界 条 件 下 多 孔 介 质气体渗流 方 程 的 解 析 解，并 结 合 数 值 算 法 对 所 得 结果进行验证，分析了 Ｋｌｉｎｋｅｎｂｅｒｇ效应对气体渗流 的影响。
渗流_水平井流耦合数学模型和数值模拟
( 2)
图 1 有侧向流入的水平井流模 型简图 Fig 1 Schematic of horizontal well model with side flow
A 为水平井断面面积。根据动量定理 , 可得到 x 长
度水平井流体的运动方程 ( P 1 - P 2) A 式中 S x =
2 A ( U2 2 - U1 ) +
图3
层流压强水头变化
图4
层流流速变化
Fig 3 Pressure head variation of side flow
表 1 模型参数 Table 1 Parameter values of model
参数 模型长度 模型宽度 模型高度 模型顶部水头 孔隙渗透系数 取值 400 cm 50 cm 105 cm 110 cm 4 103cm/ d 参数 管道长度 管道直径 管道当量粗糙度 管道糙率 迭代精度 取值 400 cm 10 cm 0 6 cm 0 015 10 7
[ 22 23]
( 8) 为修正后的摩擦系数。另外, 在层流时也 ( 9)
常用到 Ouyang Liang Biao 提出的摩擦系数经验修正公式 0 6142 = ( 1 + 0 043 04Rew ) 式中 Rew 为管壁雷诺数
[ 23]
。
第 6期
董贵明 , 等 : 渗流 水平 井流耦合数学模型和数值模拟
摘要 : 在有侧向流入 ( 流出 ) 的水平井流运动和连续性方程的基础上 , 建立了渗流与水平井流耦合数学 模型 , 主要针 对混合水头损失问题 , 进行了渗流与水平井流耦合数值计算 , 分 析了常用的 摩擦系数 修正方法对 计算结果的 影响。 数值模拟计算表明 : 层流或者光滑紊流流态时 , 不同的摩擦系数修正方法对结果影响较 小 ; 粗糙 紊流流态时 , 不同 的摩擦系数修正方法对结果影响较大。当水平井流是层流 光滑紊流流 态时 , 使用了 一种新 的流态 分界点识 别方法 对水平井流态进行判断 , 计算结果证实了方法的可靠性、合理性。 关 键 词 : 渗流 水平井流耦合模型 ; 数值模拟 ; 流态 ; 水头损失 ; 摩擦系数 文献标识码 : A 文章编号 : 1001 6791( 2009) 06 0830 08
煤矿采空区瓦斯渗流规律及其数值模拟研究
煤矿采空区瓦斯渗流规律及其数值模拟研究煤矿采空区上覆岩层结构和移动规律分析综放工艺在开采高含量瓦斯厚煤层的推广应用中之所以遇到困难，往往是由于综放面上隅角瓦斯易超限，从而被迫断电撤人、中断生产所导致的。
上隅角瓦斯的主要来源一是工作面煤壁释放出的瓦斯，二是采煤工作面新采落下来的煤炭中散发出来的瓦斯，三是从采空区涌出的瓦斯，其中采空区涌出瓦斯是主要的来源。
由于采动影响在采动断裂带形成的破断裂隙和离层裂隙，采动裂隙网络与采空区相连通形成采动断裂带，由于瓦斯的升浮、扩散和渗透作用，在采动断裂带形成瓦斯富集区，这是瓦斯抽采的重点区域。
因此，要研究采空区内瓦斯的渗流规律，有必要先研究采空区岩体的垮落特征，按照采场覆岩横向采动特征，将采空区按照自然堆积区、载荷影响区和压实稳定区在横向进行划分，弄清各区碎胀系数、空隙率的分布特点；研究采空区上覆岩层采动断裂带的高度、碎胀系数及空隙率等特征，以便较全面地分析和研究采空区内空气—瓦斯混合气体在冒落带和采动断裂带内的渗流规律。
采空区瓦斯流场数学模型研究煤矿采空区内的瓦斯流动情况，建立起瓦斯流场的数学模型，对于认识采空区内瓦斯的真实流动状况以及对于进行数值模拟都有重要的基础意义。
垮落带之上的采动断裂带，在存在破断裂隙和离层裂隙相互贯通的同时，煤岩体内的裂隙还会与综放采场和采空区连通。
研究瓦斯在采动断裂带内的渗流、升浮和扩散原理，可以为解释采动断裂带是瓦斯聚集带，为其内布置钻孔抽采、巷道排放等瓦斯治理技术提供科学依据。
求解方法的选择FLUENT提供三种不同的解格式：分离解；隐式耦合解；显式耦合解。
三种解法都可以在很大流动范围内提供准确的结果，但是它们也各有优缺点。
分离解和耦合解方法的区别在于，连续性方程、动量方程、能量方程以及组分方程的解的步骤不同，分离解是按顺序解，耦合解是同时解。
两种解法都是最后解附加的标量方程（比如：湍流或辐射）。
隐式解法和显式解法的区别在于线化耦合方程的方式不同。
煤层气开发数值模拟-CMG
CMG煤层气数值模拟软件介绍加拿大计算机模拟软件集团煤层气模拟的基本概念 层气模拟 本概2煤层构造和煤层气流动机理 原生孔隙：基质 次生孔隙：割理（裂缝）3煤层中的流动状态CH4 CO2 N2 煤层气 次开采 煤层气一次开采 CO2 提高采收率 (CO2-ECBM) N2 提高采收率 (N2-ECBM) 烟道气提高采收率4煤层气一次开采机理 通过排水降低割理压力 煤层气从基质中解吸附出来 扩散到节理/裂缝当中 煤层气从基质中解吸附出来，扩散到节理 煤层气和水流动到井筒 在裂缝中为达西流动 • 面割理和端割理 • 面割理和端割理的垂直连接部分 在大型裂缝中的达西流动或者管流 • 大型节理 大型 • 次生裂缝z1 在生产油管和井筒中的管流 裂缝渗透率受基质影响CH4Coal Matrix煤层割理和裂缝CH4 CH4 H2O CH45CH4CH4Slide 5 z1zll, 2012-3-14提高煤层气采收率机理ECBM烟道气CO2分离 注入N2N2 CH4 出售煤 深部煤层绿色电厂CO2CH4CH4CH4• 提高煤层气采收率 • 温室气体( (GHG) )封存6煤层属性：多重孔隙度系统 原生孔隙度系统（煤层基质） 微孔隙度 (< 2 nm) 中孔隙度 (2 – 50 nm)+ 非常低的流通能力：渗透率在微达西范围 只有扩散流动 次生孔隙度系统（煤层节理） 宏观孔隙度 (> 50 nm) 天然裂缝 更强的流通能力：渗透率在毫达西范围 达西流动7在GEM中煤作为多重孔隙系统 需要多重孔隙度模型例如 DUALPOR SHAPE GK在裂缝（割理）系统中为标准的达西流动裂缝间距, I,J,K IJK例如 DIFRAC CON 0.2 或 DIFRAC ALL array MATRIX 表示基质系统 FRACTURE 表示裂缝系统 基质中允许非达西流动 煤层气从基质扩散到裂缝（注意：如果基质渗透率被指定为0，那么基质到裂缝之间没有流动。
地球物理场中煤层气渗流控制方程及其数值解
(1)
采用实验模拟了加载和卸载
[7]
三维应力下煤层气的渗流运移过程；徐剑良等 研 究了煤层气开采过程中的煤岩变形与渗流运移的关 系。但以往研究中的渗流方程大多未能将地应力和 地温列入方程进行计算，缺乏实验中应力与原岩地 应力的联系过程。 目前以重庆大学为代表，地球物理场中煤层气
式中： 为流动甲烷密度(g/cm3)； vx ， v y ， vz 分别 为直角坐标系中 x ，y ，z 方向上甲烷流动速度(cm/s)；
， 9]
，但均未进行实验
1
引
言
测试和计算检验，同时部分参数的计算方法上尚存 在一些不足，如王宏图等[9]得到的应力和温度理论 上同时为 0 时的初始渗透率与实验存在较大误差。 温度与渗透率的关系至今仍然没有统一的认识，甚 至有相反的结论， 如贺玉龙等[10 度升高而降低，而梁 冰等 出统一的合理解释。 针对目前地球物理场中煤层气渗流运移理论及 实验中存在的不足，本研究将地应力和地温列入计 算，建立了新的煤层气渗流控制方程，并进行了考 虑温度和应力的各项实验测试。以重庆沥鼻峡矿区 为研究目标区，根据测试参数，结合实际地质资料， 进行了地球物理场中煤层气渗流方程的一维、二维 数值模拟。尝试考虑地温和地应力因素，动态反映 煤层气渗流运移能力对煤层气储集和富集能力的影 响，希望寻找到提高抽采率的有效途径。
收稿日期：2008–03–22；修回日期：2008–07–20 基金项目：国家自然科学基金创新研究群体科学基金项目(50621403)；重庆市院士基金专项资助项目(CSTC2005AB007)；中国博士后科学基金资助项 目(2005037020) 作者简介：李志强(1975–)，男，1999 年毕业于焦作工学院采矿工程系矿山通风与安全专业，现为博士研究生，主要从事煤层气开发及矿井瓦斯防治 等方面的研究工作。E-mail：maximos592@