GEOFLUIDSgydF4y2Ba GeofluidsgydF4y2Ba 1468 - 8123gydF4y2Ba 1468 - 8115gydF4y2Ba HindawigydF4y2Ba 10.1155 / 2021/5510566gydF4y2Ba 5510566gydF4y2Ba 研究文章gydF4y2Ba 瓦斯抽放的案例研究优化基于硅镁层有效钻孔间距的煤矿gydF4y2Ba https://orcid.org/0000 - 0002 - 5549 - 0639gydF4y2Ba 霁gydF4y2Ba 明gydF4y2Ba 1gydF4y2Ba https://orcid.org/0000 - 0002 - 7231 - 5584gydF4y2Ba 太阳gydF4y2Ba Zhong-guanggydF4y2Ba 2gydF4y2Ba 3gydF4y2Ba 4gydF4y2Ba 太阳gydF4y2Ba 魏gydF4y2Ba 3gydF4y2Ba 雪gydF4y2Ba 易gydF4y2Ba 1gydF4y2Ba 深部煤炭资源开采的重点实验室gydF4y2Ba 中国教育部和中国矿业大学&技术gydF4y2Ba 徐州221116gydF4y2Ba 中国gydF4y2Ba 2gydF4y2Ba 国家重点实验室检测气体的灾难gydF4y2Ba 预防和紧急控制gydF4y2Ba 重庆400037gydF4y2Ba 中国gydF4y2Ba 3gydF4y2Ba 中国煤炭技术和工程集团重庆研究所gydF4y2Ba 重庆400039gydF4y2Ba 中国gydF4y2Ba ccteg.cngydF4y2Ba 4gydF4y2Ba 国家重点实验室煤矿灾害动力学与控制gydF4y2Ba 重庆大学gydF4y2Ba 重庆400044gydF4y2Ba 中国gydF4y2Ba cqu.edu.cngydF4y2Ba 2021年gydF4y2Ba 8gydF4y2Ba 9gydF4y2Ba 2021年gydF4y2Ba 2021年gydF4y2Ba 26gydF4y2Ba 2gydF4y2Ba 2021年gydF4y2Ba 15gydF4y2Ba 6gydF4y2Ba 2021年gydF4y2Ba 14gydF4y2Ba 8gydF4y2Ba 2021年gydF4y2Ba 8gydF4y2Ba 9gydF4y2Ba 2021年gydF4y2Ba 2021年gydF4y2Ba 版权©2021明霁et al。gydF4y2Ba 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。gydF4y2Ba

根据煤层的孔隙率和渗透率的动态表达式派生,multiphysical场耦合数值模型的气体迁移交互作用下的应力场和渗流场成立。瓦斯抽放工程3号煤层由硅镁层煤炭工业有限公司有限公司被选为研究对象。在各种不同的排水时间排水洞的位置考虑,结合之前的硅镁层煤矿1207工作面,气体混合布局排水系统方案具有有效的钻孔间距是通过COMSOL多重物理量模拟。此外,一系列的现场工业试验进行了验证研究结果,揭示全面考虑煤的萃取时间和优化开采布局的孔可以有效提高工程经济效益。gydF4y2Ba

 九龙坡区科技计划项目gydF4y2Ba 2020 - 02 - 005 ygydF4y2Ba 重庆市自然科学基金gydF4y2Ba cstc2020jcyj-msxmX0972gydF4y2Ba cstc2019jcyj-msxmX0633gydF4y2Ba 优先级的学术程序开发江苏高等教育机构gydF4y2Ba 1。介绍gydF4y2Ba

煤层瓦斯抽放是我的一个重要措施来控制气体(gydF4y2Ba 1gydF4y2Ba,gydF4y2Ba 2gydF4y2Ba]然而,大多数研究以有效的排水半径为排水洞布局的基础上,俯瞰的影响重叠排水系统和裂缝扩张没有采取不同排水的排水洞的影响在不同煤层的位置考虑在内。此外,大多数研究采用单个钻孔的布局,这常常导致一些常见的问题,比如浪费资源,不均匀的瓦斯抽放,达到相关标准的失败。许多现有的研究没有考虑建立multiphysical场耦合模型的气体迁移的相互作用下应力场和渗流场和忽视的动态特性煤层的渗透率和孔隙度。所有前面提到的缺乏在先前的研究有必要进行深入优化研究煤层瓦斯抽放。gydF4y2Ba

在这篇文章中,3号煤层的瓦斯抽放优化硅镁层煤炭工业有限公司路安集团有限公司进行了研究。基于煤层应力的变化、气体压力、气体吸附和解吸,往往产生重要影响煤层渗透率和孔隙度、煤层渗透率和孔隙度的动态表达式推导出用来建立移民multiphysical场耦合的数值模型的交互作用下应力场和渗流场。此外,在各种不同的排水时间排水洞的位置考虑,结合推进情况的硅镁层煤矿1207工作面,气体混合布局排水系统方案具有有效的钻孔间距由COMSOL获得多重物理量模拟。此外,进行了一系列现场工业试验的目的是验证研究。gydF4y2Ba

 2。气固耦合分析煤层气体迁移gydF4y2Ba 2.1。煤层气的迁移规律gydF4y2Ba

两种类型的气流在煤层已确定包括层流扩散流和针对不同类型的煤层(gydF4y2Ba 3gydF4y2Ba- - - - - -gydF4y2Ba 5gydF4y2Ba]。一般来说,具有微孔结构的煤层或低渗透、扩散流中起着至关重要的作用在微孔结构或高渗透的煤层,层流控制(gydF4y2Ba 6gydF4y2Ba- - - - - -gydF4y2Ba 8gydF4y2Ba]。在这项研究中,气体流动过程中排水作为层流,遵循达西定律。gydF4y2Ba

许多因素可以影响煤层气的运移,存储、地质条件等发生煤层,煤和岩石的力学性能,煤层渗透率和孔隙度,煤层气体的吸附容量,气体压力。在这些因素中,气体压力、气体吸附能力的煤层,煤层孔隙度,渗透率是主要的。气体压力的驱动力是煤层气体迁移而气体吸附能力决定煤层储存气体的能力。与此同时,孔隙度和渗透率是光滑的指标气体迁移在煤层gydF4y2Ba 9gydF4y2Ba- - - - - -gydF4y2Ba 11gydF4y2Ba]。gydF4y2Ba

3号煤层的瓦斯压力测定为0.31 MPa在硅镁层矿作为本研究的基本的气体压力,尽管气体压力测量保持低于0.74 MPa的临界值。气体吸附常数a和b的# 3煤层被确定在17.52米左右gydF4y2Ba3 /gydF4y2Bat和0.75 MPagydF4y2Ba1gydF4y2Ba,分别,这被认为是低,这表明3号煤层的气体吸附能力弱。在现场排水,可以采取几项措施促进排水和提取包括选择适当的负压引流孔,孔的布局的优化,促进气体的解吸,增加气体的扩散和迁移速度。gydF4y2Ba

 2.2。动态表达式的推导煤层的孔隙度和渗透率gydF4y2Ba

在天然气开采过程中,气体往往是眠和迁移,使煤炭骨架应力变化,导致小变形,孔隙度和渗透率相应地改变。根据现有的研究和田间条件下硅镁层我和2 d模型的数值模拟计算,平面应变gydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba 介绍,导致修改动态表达式煤层孔隙度(gydF4y2Ba 12gydF4y2Ba- - - - - -gydF4y2Ba 14gydF4y2Ba]。gydF4y2Ba (1)gydF4y2Ba φgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba ∆gydF4y2Ba pgydF4y2Ba kgydF4y2Ba 年代gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba 煤的初始孔隙度,gydF4y2Ba εgydF4y2Ba vgydF4y2Ba 是体积应变,gydF4y2Ba ∆gydF4y2Ba PgydF4y2Ba 气体压力的变化,gydF4y2Ba ∆gydF4y2Ba pgydF4y2Ba =gydF4y2Ba PgydF4y2Ba −gydF4y2Ba PgydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba PgydF4y2Ba 是真正的气体压力,gydF4y2Ba PgydF4y2Ba 0gydF4y2Ba 是初始气体压力,gydF4y2Ba kgydF4y2Ba 年代gydF4y2Ba 骨架是煤岩体的弹性模量。gydF4y2Ba

根据Kozeny-Carman公式,结合以前的动态表达式孔隙度、渗透率的修改后的动态表达式可以得到如下:gydF4y2Ba (2)gydF4y2Ba kgydF4y2Ba egydF4y2Ba =gydF4y2Ba kgydF4y2Ba 0gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba +gydF4y2Ba ∆gydF4y2Ba pgydF4y2Ba /gydF4y2Ba kgydF4y2Ba 年代gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba 3gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba kgydF4y2Ba 0gydF4y2Ba 是煤炭和初始磁导率gydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba 平面应变。gydF4y2Ba

 2.3。Multiphysical场耦合的煤层气运移模式gydF4y2Ba

本研究关注的是变形的煤层由于强加的力量和气体流量(gydF4y2Ba 15gydF4y2Ba- - - - - -gydF4y2Ba 18gydF4y2Ba]。因此,以下假设了:gydF4y2Ba

瓦斯抽放孔形成的流场在一定的负压下径向流场gydF4y2Ba

气体的吸附内容所描述的是朗缪尔方程gydF4y2Ba

气体流量是由达西定律描述gydF4y2Ba

气体的渗流过程被认为是一种理想的等温过程gydF4y2Ba

煤层是各向同性的线性弹性变形,只有一个小gydF4y2Ba

水的影响在煤层瓦斯抽放是被忽视的gydF4y2Ba

考虑到上述假设,煤层气的渗流运动的理想气体应遵循如下列出的控制方程的气体流(gydF4y2Ba 13gydF4y2Ba,gydF4y2Ba 14gydF4y2Ba,gydF4y2Ba 19gydF4y2Ba- - - - - -gydF4y2Ba 25gydF4y2Ba]。gydF4y2Ba

2.3.1。质量守恒定律gydF4y2Ba

气体在煤层通常需要两种形式包括自由州和吸附状态。尽管状态,源和煤层气体的总量保持不变。换句话说,气体迁移应该首先遵循质量守恒定律:gydF4y2Ba (3)gydF4y2Ba ∂gydF4y2Ba 米gydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba +gydF4y2Ba ∇gydF4y2Ba ∙gydF4y2Ba ρgydF4y2Ba ggydF4y2Ba 问gydF4y2Ba ggydF4y2Ba =gydF4y2Ba 问gydF4y2Ba pgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba 米gydF4y2Ba 气体含量(公斤/米gydF4y2Ba3gydF4y2Ba),gydF4y2Ba ρgydF4y2Ba ggydF4y2Ba 是气体密度(公斤/米gydF4y2Ba3gydF4y2Ba),gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba 是气体的达西渗流速度(米/秒),gydF4y2Ba 问gydF4y2Ba pgydF4y2Ba 源或汇项目(公斤/ (mgydF4y2Ba3gydF4y2Ba·s))gydF4y2Ba tgydF4y2Ba 是时间变量(s)。gydF4y2Ba

 2.3.2。气体的状态方程gydF4y2Ba

气体渗流运动被视为理想气体,满足理想气体的状态方程:gydF4y2Ba (4)gydF4y2Ba ρgydF4y2Ba ggydF4y2Ba =gydF4y2Ba 米gydF4y2Ba ggydF4y2Ba PgydF4y2Ba RgydF4y2Ba TgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba 米gydF4y2Ba ggydF4y2Ba 是气体的分子量(公斤/ kmol),gydF4y2Ba PgydF4y2Ba 是气体压力(Pa),gydF4y2Ba RgydF4y2Ba 理想气体常数(kJ / (kmol·K)),然后呢gydF4y2Ba TgydF4y2Ba 是绝对温度(KgydF4y2BaogydF4y2Ba)。gydF4y2Ba

 2.3.3。朗缪尔方程gydF4y2Ba

煤层是一种多孔介质,有许多孔隙和裂缝,在气体吸附状态和自由州的共存。朗缪尔方程描述气内容定义了瓦斯吸附和解吸瓦斯压力之间的关系:gydF4y2Ba (5)gydF4y2Ba 米gydF4y2Ba =gydF4y2Ba 米gydF4y2Ba ggydF4y2Ba RgydF4y2Ba TgydF4y2Ba φgydF4y2Ba PgydF4y2Ba 0gydF4y2Ba +gydF4y2Ba 一个gydF4y2Ba bgydF4y2Ba ρgydF4y2Ba 年代gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba bgydF4y2Ba PgydF4y2Ba PgydF4y2Ba 2gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba φgydF4y2Ba 孔隙度,gydF4y2Ba PgydF4y2Ba 0gydF4y2Ba 大气压力(Pa),gydF4y2Ba 一个gydF4y2Ba 朗缪尔常数(mgydF4y2Ba3gydF4y2Ba/公斤),gydF4y2Ba bgydF4y2Ba 朗缪尔常数(gydF4y2Ba 巴勒斯坦权力机构gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba ),gydF4y2Ba ρgydF4y2Ba 年代gydF4y2Ba 是煤层密度(公斤/米gydF4y2Ba3gydF4y2Ba)。gydF4y2Ba

 2.3.4。达西定律gydF4y2Ba

当压力梯度分布确定的煤层瓦斯压力,气体迁移会出现,这通常是一个线性渗流和遵循达西定律:gydF4y2Ba (6)gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba =gydF4y2Ba −gydF4y2Ba kgydF4y2Ba egydF4y2Ba μgydF4y2Ba ggydF4y2Ba ∇gydF4y2Ba pgydF4y2Ba +gydF4y2Ba ρgydF4y2Ba ggydF4y2Ba ggydF4y2Ba ∇gydF4y2Ba zgydF4y2Ba 2gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba 是达西流速度(米/秒),gydF4y2Ba kgydF4y2Ba egydF4y2Ba 煤层的渗透率(mgydF4y2Ba2gydF4y2Ba),gydF4y2Ba μgydF4y2Ba ggydF4y2Ba 是气体动态粘度(Pa·s),然后呢gydF4y2Ba ggydF4y2Ba 重力加速度(gydF4y2Ba 米gydF4y2Ba ·gydF4y2Ba 年代gydF4y2Ba −gydF4y2Ba 2gydF4y2Ba )。gydF4y2Ba

天然气开采的过程中有两个物理现象包括液体和气体的渗流运动的微小变形煤。煤体的应力和变形导致孔隙特征变化,然后影响到气体渗流运动。因此,气体渗流往往改变煤的身体内的气体压力,导致不同有效应力。两种现象就像一个连锁反应的相互作用相互约束和影响。gydF4y2Ba

煤层孔隙度的动态表达式代入方程(gydF4y2Ba 5gydF4y2Ba),然后gydF4y2Ba (7)gydF4y2Ba 米gydF4y2Ba =gydF4y2Ba 米gydF4y2Ba ggydF4y2Ba RgydF4y2Ba TgydF4y2Ba 1gydF4y2Ba PgydF4y2Ba 0gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba PgydF4y2Ba 0gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba ∆gydF4y2Ba PgydF4y2Ba kgydF4y2Ba 年代gydF4y2Ba +gydF4y2Ba 一个gydF4y2Ba bgydF4y2Ba ρgydF4y2Ba 年代gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba bgydF4y2Ba PgydF4y2Ba PgydF4y2Ba 2gydF4y2Ba 。gydF4y2Ba

煤层渗透率的动态表达式代入方程(gydF4y2Ba 6gydF4y2Ba),然后gydF4y2Ba (8)gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba =gydF4y2Ba −gydF4y2Ba kgydF4y2Ba 0gydF4y2Ba ∇gydF4y2Ba pgydF4y2Ba +gydF4y2Ba ρgydF4y2Ba ggydF4y2Ba ggydF4y2Ba ∇gydF4y2Ba zgydF4y2Ba 2gydF4y2Ba μgydF4y2Ba ggydF4y2Ba 1gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba +gydF4y2Ba ∆gydF4y2Ba pgydF4y2Ba /gydF4y2Ba kgydF4y2Ba 年代gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba 3gydF4y2Ba 。gydF4y2Ba

用公式(gydF4y2Ba 4gydF4y2Ba)和派生公式(gydF4y2Ba 7gydF4y2Ba)和(gydF4y2Ba 8gydF4y2Ba)到公式(gydF4y2Ba 3gydF4y2Ba),得到以下结果:gydF4y2Ba (9)gydF4y2Ba ∂gydF4y2Ba 米gydF4y2Ba ggydF4y2Ba /gydF4y2Ba RgydF4y2Ba TgydF4y2Ba 1gydF4y2Ba /gydF4y2Ba PgydF4y2Ba 0gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba /gydF4y2Ba PgydF4y2Ba 0gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba ∆gydF4y2Ba PgydF4y2Ba /gydF4y2Ba kgydF4y2Ba 年代gydF4y2Ba +gydF4y2Ba 一个gydF4y2Ba bgydF4y2Ba ρgydF4y2Ba 年代gydF4y2Ba /gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba bgydF4y2Ba PgydF4y2Ba PgydF4y2Ba 2gydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba −gydF4y2Ba ∇gydF4y2Ba ·gydF4y2Ba 米gydF4y2Ba ggydF4y2Ba PgydF4y2Ba kgydF4y2Ba 0gydF4y2Ba ∇gydF4y2Ba pgydF4y2Ba +gydF4y2Ba ρgydF4y2Ba ggydF4y2Ba ggydF4y2Ba ∇gydF4y2Ba zgydF4y2Ba 2gydF4y2Ba RgydF4y2Ba TgydF4y2Ba μgydF4y2Ba ggydF4y2Ba 1gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 年代gydF4y2Ba +gydF4y2Ba ∆gydF4y2Ba pgydF4y2Ba /gydF4y2Ba kgydF4y2Ba 年代gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba 3gydF4y2Ba =gydF4y2Ba 问gydF4y2Ba pgydF4y2Ba 。gydF4y2Ba

公式(gydF4y2Ba 9gydF4y2Ba)是气固耦合模型最终派生。模型嵌入到COMSOL多重物理量瓦斯抽放软件优化。gydF4y2Ba

 3所示。计算模型的引入gydF4y2Ba 3.1。计算模型的参数设置gydF4y2Ba

在建立数值模型之前,相关参数的确定是至关重要的。根据相关文献和地质力学测试,测试和调整系数后COMSOL多重物理量5.0内置模块,得到了模型中所使用的主要参数,如表所示gydF4y2Ba 1gydF4y2Ba。gydF4y2Ba

模型参数。gydF4y2Ba

项目gydF4y2Ba 象征gydF4y2Ba 价值gydF4y2Ba
气体密度(公斤/米gydF4y2Ba3gydF4y2Ba)gydF4y2Ba rhoggydF4y2Ba 0.716gydF4y2Ba
气体动力学粘度(Pa·s)gydF4y2Ba 六世gydF4y2Ba 1.08gydF4y2Ba egydF4y2Ba 5gydF4y2Ba
等效压缩系数矩阵(1 / Pa)gydF4y2Ba χgydF4y2Ba pgydF4y2Ba 2.18gydF4y2Ba egydF4y2Ba 3gydF4y2Ba
流体压缩系数(1 / Pa)gydF4y2Ba chifgydF4y2Ba 1.3942gydF4y2Ba egydF4y2Ba 5gydF4y2Ba
煤的密度(公斤/米gydF4y2Ba3gydF4y2Ba)gydF4y2Ba rhocgydF4y2Ba 1470年gydF4y2Ba
初始孔隙度的煤gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba 0.0318gydF4y2Ba
泊松比的煤gydF4y2Ba μgydF4y2Ba 0.33gydF4y2Ba
初始渗透率煤(mgydF4y2Ba2gydF4y2Ba)gydF4y2Ba kgydF4y2Ba 0gydF4y2Ba 0.907gydF4y2Ba egydF4y2Ba -15年gydF4y2Ba
杨氏模量的煤(MPa)gydF4y2Ba EgydF4y2Ba 3000年gydF4y2Ba
初始煤气压力(MPa)gydF4y2Ba PgydF4y2Ba 0gydF4y2Ba 0.31gydF4y2Ba
标准大气压力(Pa)gydF4y2Ba PgydF4y2Ba 一个gydF4y2Ba 1.01gydF4y2Ba egydF4y2Ba 5gydF4y2Ba
Biot-Willis系数gydF4y2Ba αgydF4y2Ba BgydF4y2Ba 0.801gydF4y2Ba
3.2。模型的建立gydF4y2Ba

# 3煤层原始瓦斯压力的硅镁层矿山为0.31 MPa。由于气体浓度测量误差大,原来的气体压力降低30%被选为标准的模拟,导致气体排水标准压力设定在0.22 MPa。gydF4y2Ba

为了提供直接指导现场施工,介绍了有效的钻孔间距的模拟。整个煤层的瓦斯压力可以减少最大钻孔间距值在标准压力值排水后一段时间。换句话说,通过数值模拟,可以确定最合理的钻孔间距值,这不仅符合减少气体压力的需要,也减少了现场施工的难度和成本。gydF4y2Ba

方法类似于“水桶理论”是用来计算的有效钻孔间距,旨在找出最低的地区钻孔抽水的叠加效应仿真计算模型,确保区域内的压力下降低于标准的价值。gydF4y2Ba

多孔弹性模块和达西定律模块(gydF4y2Ba 13gydF4y2Ba,gydF4y2Ba 14gydF4y2Ba,gydF4y2Ba 26gydF4y2Ba在COMSOL本模拟中使用,所需的变量和参数输入到软件。根据实际钻孔瓦斯抽放过程,一个理想化的二维计算模型gydF4y2Ba 长度gydF4y2Ba ×gydF4y2Ba 高度gydF4y2Ba 的gydF4y2Ba 16gydF4y2Ba 米gydF4y2Ba ×gydF4y2Ba 6.6gydF4y2Ba 米gydF4y2Ba 建立如图gydF4y2Ba 1gydF4y2Ba。排水井被安排在中间模型底部的固定约束。同时,上,左,右边界的边界载荷8.25 MPa,实际的围压条件。各方的模型是密封的初始气体压力为0.31 MPa在模型中。gydF4y2Ba

计算模型。gydF4y2Ba

作为显示在图gydF4y2Ba 1gydF4y2Ba,钻孔是安排在2 d的计算模块。单一变量法后,提取的有效提取半径对应于不同直径钻孔(95毫米,113毫米和133毫米),萃取压力(20 kPa, 25个kPa, 30 kPa, 35 kPa, 40 kPa, 45 kPa),和提取时间(90 d、180 d和360 d)。并给出了详细的计算表gydF4y2Ba 2gydF4y2Ba。gydF4y2Ba

单孔模型的仿真结果。gydF4y2Ba

单身生gydF4y2Ba 有效排水半径(米)gydF4y2Ba
孔直径gydF4y2Ba 负压gydF4y2Ba 90 dgydF4y2Ba 180 dgydF4y2Ba 360 dgydF4y2Ba
95毫米gydF4y2Ba 20 kPagydF4y2Ba 0.45gydF4y2Ba 0.57gydF4y2Ba 0.72gydF4y2Ba
95毫米gydF4y2Ba 25个kPagydF4y2Ba 0.46gydF4y2Ba 0.57gydF4y2Ba 0.73gydF4y2Ba
95毫米gydF4y2Ba 30 kPagydF4y2Ba 0.46gydF4y2Ba 0.58gydF4y2Ba 0.74gydF4y2Ba
95毫米gydF4y2Ba 35个kPagydF4y2Ba 0.47gydF4y2Ba 0.59gydF4y2Ba 0.75gydF4y2Ba
95毫米gydF4y2Ba 40个kPagydF4y2Ba 0.47gydF4y2Ba 0.60gydF4y2Ba 0.76gydF4y2Ba
95毫米gydF4y2Ba 45 kPagydF4y2Ba 0.48gydF4y2Ba 0.61gydF4y2Ba 0.77gydF4y2Ba
113毫米gydF4y2Ba 20 kPagydF4y2Ba 0.48gydF4y2Ba 0.60gydF4y2Ba 0.76gydF4y2Ba
113毫米gydF4y2Ba 25个kPagydF4y2Ba 0.48gydF4y2Ba 0.61gydF4y2Ba 0.77gydF4y2Ba
113毫米gydF4y2Ba 30 kPagydF4y2Ba 0.49gydF4y2Ba 0.62gydF4y2Ba 0.78gydF4y2Ba
113毫米gydF4y2Ba 35个kPagydF4y2Ba 0.49gydF4y2Ba 0.62gydF4y2Ba 0.79gydF4y2Ba
113毫米gydF4y2Ba 40个kPagydF4y2Ba 0.50gydF4y2Ba 0.63gydF4y2Ba 0.80gydF4y2Ba
113毫米gydF4y2Ba 45 kPagydF4y2Ba 0.51gydF4y2Ba 0.64gydF4y2Ba 0.81gydF4y2Ba
133毫米gydF4y2Ba 20 kPagydF4y2Ba 0.50gydF4y2Ba 0.63gydF4y2Ba 0.80gydF4y2Ba
133毫米gydF4y2Ba 25个kPagydF4y2Ba 0.51gydF4y2Ba 0.64gydF4y2Ba 0.81gydF4y2Ba
133毫米gydF4y2Ba 30 kPagydF4y2Ba 0.51gydF4y2Ba 0.65gydF4y2Ba 0.82gydF4y2Ba
133毫米gydF4y2Ba 35个kPagydF4y2Ba 0.52gydF4y2Ba 0.66gydF4y2Ba 0.83gydF4y2Ba
133毫米gydF4y2Ba 40个kPagydF4y2Ba 0.53gydF4y2Ba 0.67gydF4y2Ba 0.84gydF4y2Ba
133毫米gydF4y2Ba 45 kPagydF4y2Ba 0.53gydF4y2Ba 0.67gydF4y2Ba 0.85gydF4y2Ba

根据表gydF4y2Ba 2gydF4y2Ba,随着提取时间、压力和提取钻孔直径增加,单一的排泄孔的有效半径增加,展示一个积极的相关性。gydF4y2Ba

目前,通常采用萃取钻孔直径大小的at113毫米和133毫米直径的25个kPa和40之间kPa萃取压力。提取通常会消耗大量的时间,超过1年甚至2到3年。由于相对较厚煤在硅镁层密封,采用多个提取洞多个层。提取孔通常安排在三种方式包括三花风格,四朵风格,和five-flower风格,显示在图gydF4y2Ba 2gydF4y2Ba。gydF4y2Ba

三安排排水洞。gydF4y2Ba

一些模拟各种条件进行包括提取钻孔直径113毫米和133毫米的抽汽压力下25 kPa, 30 kPa, 35 kPa, 40 kPa。提取时间被设定为360 d。计算的有效提取半径是显示在图gydF4y2Ba 3gydF4y2Ba在下面。gydF4y2Ba

有效孔之间的距离不同吸力下压力。gydF4y2Ba

提取的部分分为巷道。gydF4y2Ba

根据图gydF4y2Ba 3gydF4y2Ba,直径133米的排水洞证明更好的提取。提取的负压孔35 kpa的煤炭开采实践。提取的five-flower布局洞比三花布局而优越的三花布局优越四朵布局。gydF4y2Ba

 4所示。的模拟和分析1207工作面硅镁层煤gydF4y2Ba 4.1。仿真过程分析gydF4y2Ba

罢工的1207工作面长度3号煤层的硅镁层我的是220的倾斜长度1092米。生产过程中,为了确保生产的安全性,进行了天然气开采。然而,这种泵方案单一而不考虑具体情况。gydF4y2Ba

如前所述,排水时间大大影响了瓦斯抽放。不同的排水时间应采取不同的排水洞布局方案,如不同钻孔间距或模式。没有。1207工作面先进、不同提取时间的不同位置中提取矿山巷道钻孔。在这种情况下,不同部分应该划分基于不同的提取时间,允许个人提取条件设计对应于不同的领域。gydF4y2Ba

根据行车速度的综合分析,工作面布局,和矿业采矿巷道的速度。1207工作面,考虑到排水效果,煤层瓦斯抽放的没有。1207工作面在仿真可以分为三个部分,如图gydF4y2Ba 4gydF4y2Ba。三个部分的平均提取时间是110 d, 210 d和310 d,分别。gydF4y2Ba

平均萃取时间、负压和孔直径的三个部分,优化天然气开采的三个部分也详细研究。优化的计算模型的提取变化单列孔布局five-row模式洞布局,直到得到最优结果。根据提取的实验计划,整个的目标层的# 3煤提取0.22 MPa。gydF4y2Ba

 4.2。仿真结果的分析gydF4y2Ba 4.2.1。准备第一个和中间部分的分析gydF4y2Ba

第一部分的萃取时间和中间部分被设定为110天,210天,分别。根据仿真计算,煤层气体压力可以降低标准的价值在一个固定的时间内只有当five-flower布局采用。由于不同的提取时间,有效的钻井两部分不同的空间。仿真过程中的气体压力云图所示的数字gydF4y2Ba 5(一个)gydF4y2Ba和gydF4y2Ba 5 (b)gydF4y2Ba。gydF4y2Ba

气体压力云从模拟图。gydF4y2Ba

第一部分gydF4y2Ba

中间部分gydF4y2Ba

根据图gydF4y2Ba 5gydF4y2Ba更多的天然气开采,煤层中气体的压强有经验的不断变化。气体压力下降最接近排水洞。进一步排泄孔的距离会导致更少的气体压力下降到阈值距离达成的气体压力保持在初始值,表明煤炭以外的边界并不影响提取的影响。为了描述气体压力分布five-flower安排,结构线的值线图如图gydF4y2Ba 6gydF4y2Ba云是插入到模拟气体压力图,可以生动地说明中的气体压力值的分布区域five-flower包裹的洞。值线图如图gydF4y2Ba 7gydF4y2Ba。gydF4y2Ba

的结构线值线图。gydF4y2Ba

值线图。gydF4y2Ba

根据图gydF4y2Ba 5gydF4y2Ba,最大气体压力的包装区域five-flower孔是0.22 MPa,和其他地方的气体压力低于设定值。在这个时候,两个相邻排水洞之间的距离在同一水平得出的结论是有效的钻孔间距。gydF4y2Ba

根据仿真结果,有效的钻孔间距值five-flower安排了两个部分,如表所示gydF4y2Ba 3gydF4y2Ba。gydF4y2Ba

有效的在不同部分钻孔间距值。gydF4y2Ba

道路分割gydF4y2Ba 提取时间gydF4y2Ba 水平有效的钻孔间距值gydF4y2Ba 纵向行间距gydF4y2Ba
第一次分割gydF4y2Ba 110 dgydF4y2Ba 1.15米gydF4y2Ba 2.20米gydF4y2Ba
中间分割gydF4y2Ba 210 dgydF4y2Ba 2.45米gydF4y2Ba 2.25米gydF4y2Ba
4.2.2。分析的结果结束部分gydF4y2Ba

最后部分的提取时间被设定为310天。仿真表明,煤层气体压力可以减少标准价值三花在一个固定的时间安排。仿真过程中的气体压力云图如图gydF4y2Ba 8gydF4y2Ba。gydF4y2Ba

气体压力云图。gydF4y2Ba

同样,从图gydF4y2Ba 8gydF4y2Ba更接近排泄孔导致更高的压力下降。气体压力达到峰值0.22 MPa的三角形重心。为了更好地描述气体压力的分布在三花洞,value-taking的构造线线路图如图gydF4y2Ba 9gydF4y2Ba云是插入到模拟气体压力图,可以生动地说明中的气体压力值的分布区域包装的三花洞,value-taking线图如图gydF4y2Ba 10gydF4y2Ba。gydF4y2Ba

的结构线值线图。gydF4y2Ba

值线图。gydF4y2Ba

根据图gydF4y2Ba 10gydF4y2Ba,最大气体压力在包装区域的三花洞是0.22 MPa,和其他地方的气体压力低于设定值。这时,提取两个相邻孔之间的距离在同一级别是有效的钻孔间距。gydF4y2Ba

仿真表明,水平有效的钻孔间距是2.00米的垂直行间距2.00米,三花图案布局采用最后一节。gydF4y2Ba

 4.3。混合布局的建议方案gydF4y2Ba

根据前面的计算结果,混合布局方案是专门为开发的瓦斯抽放。我1207工作面硅镁层。提出了如下所示的细节。gydF4y2Ba

第一部分gydF4y2Ba:five-flower模式布局采用水平钻孔间距为1.15米。上部排钻距地面5.50米的中间行纵向钻探3.30离开地板,和底部行drilling1.10 m离开地板。gydF4y2Ba

中间部分gydF4y2Ba:采用five-flower模式布局水平钻孔间距为2.45米。上部排钻是5.55米的地板纵向中间行钻井3.30离开地板,和底部行drilling1.05 m离开地板。gydF4y2Ba

结束部分gydF4y2Ba:采用三花模式安排水平钻孔间距为2.00米。上部排钻是4.92米离地板纵向排钻孔底部1.68米的距离。gydF4y2Ba

 5。工业测试gydF4y2Ba 5.1。详细的测试计划gydF4y2Ba

在这个工业实践,回风巷道的没有。1207工作面被选中。而不影响原来的生产和提取,提取练习测试进行了在中间部分的450米~ 485米和750米~ 785米的最后部分回风巷道。根据前面的计算结果,具体方案如下:gydF4y2Ba

在450 ~ 485回风巷道的中间部分gydF4y2Ba

five-flower孔布局采用的负压35 kPa,直径133毫米,平均长度为161 m,水平钻孔间距为2.45米。纵向,上部排孔是5.55米的距离,孔的中间行3.30米的距离,和底部排孔1.05离开地板。gydF4y2Ba

在750米~ 785米的最后部分回风巷道gydF4y2Ba

三花洞布局采用的负压35 kPa,直径133毫米,平均长度为161 m,水平钻孔间距为2.00米。纵向,上层行和地板之间的距离是4.92米,与底行和地板之间的距离为1.68米。gydF4y2Ba

 5.2。现场试验的分析gydF4y2Ba

瓦斯抽放的工业性试验方案显示,排水钻孔煤层有效、均匀。中间的残余瓦斯含量和结束部分的空气巷道降低了31.24%和32.16%,分别分析气体含量减少了46.86%和49.81%,分别和煤层瓦斯压力降低45.03%和48.33%,分别。gydF4y2Ba

 6。结论gydF4y2Ba

随着提取时间,提取的负压,并提取钻孔直径增加,有效提取单一的排泄孔的半径增加,展示一个积极的相关性。gydF4y2Ba

基于仿真的提案,包括一个有效的提取孔之间的间距。提议,提取孔排列基于有效排水洞的距离,采取不同的提取时间考虑巷道对应于不同的位置。gydF4y2Ba

云图,根据气体压力更大的压降是确定接近提取孔。气体压力保持稳定在达到阈值距离。gydF4y2Ba

气体压力随开采时间的增长,最后变得稳定。一些叠加效果已确定两个排水洞更密集的气体压降叠加区域。gydF4y2Ba

 数据可用性gydF4y2Ba

使用的数据来支持本研究的发现可以从相应的作者。gydF4y2Ba

 的利益冲突gydF4y2Ba

作者宣称没有利益冲突有关的出版。gydF4y2Ba

 确认gydF4y2Ba

本文得到了优先级的学术程序开发江苏高等教育机构和重庆市自然科学基金,中国(批准号cstc2019jcyj-msxmX0633也没有。cstc2020jcyj-msxmX0972),九龙坡区科技计划项目,中国(批准号2020 - 02 - 005 - y)。gydF4y2Ba

 杨ydF4y2Ba Z。gydF4y2Ba 王gydF4y2Ba Y。gydF4y2Ba 风扇gydF4y2Ba J。gydF4y2Ba 黄gydF4y2Ba Y。gydF4y2Ba 他gydF4y2Ba Y。gydF4y2Ba 研究关键参数的定向长钻孔布置在汽油的脸gydF4y2Ba 冲击和振动gydF4y2Ba 2021年gydF4y2Ba 2021年gydF4y2Ba 14gydF4y2Ba 5579967gydF4y2Ba 10.1155 / 2021/5579967gydF4y2Ba 郭gydF4y2Ba X。gydF4y2Ba 雪gydF4y2Ba 年代。gydF4y2Ba 李gydF4y2Ba Y。gydF4y2Ba 郑gydF4y2Ba C。gydF4y2Ba 谢gydF4y2Ba lgydF4y2Ba 研究优化关键领域在ultrathick煤层抽放钻孔的密封gydF4y2Ba 冲击和振动gydF4y2Ba 2021年gydF4y2Ba 2021年gydF4y2Ba 8gydF4y2Ba 5536196gydF4y2Ba 10.1155 / 2021/5536196gydF4y2Ba 风扇gydF4y2Ba J。gydF4y2Ba 瓦斯抽放技术的研究和应用在高气煤层工作面gydF4y2Ba 2016年gydF4y2Ba 中国煤炭研究所gydF4y2Ba 锅gydF4y2Ba R。gydF4y2Ba 加载煤的渗透率演化特征及其应用在卸压瓦斯抽放gydF4y2Ba 2014年gydF4y2Ba 中国矿业大学和技术gydF4y2Ba 周gydF4y2Ba j . Z。gydF4y2Ba 魏gydF4y2Ba c F。gydF4y2Ba 李gydF4y2Ba w·T。gydF4y2Ba 程ydF4y2Ba h·L。gydF4y2Ba saolid弹塑性力学分析的基于数学模块COMSOL多重物理量gydF4y2Ba 工程武汉大学学报gydF4y2Ba 2015年gydF4y2Ba 48gydF4y2Ba 2gydF4y2Ba 195年gydF4y2Ba 201年gydF4y2Ba 如果gydF4y2Ba lgydF4y2Ba 张gydF4y2Ba H。gydF4y2Ba 魏gydF4y2Ba J。gydF4y2Ba 李gydF4y2Ba B。gydF4y2Ba 汉gydF4y2Ba H。gydF4y2Ba 建模和气体有效扩散系数的实验在水饱和煤gydF4y2Ba 燃料gydF4y2Ba 2021年gydF4y2Ba 284年gydF4y2Ba 15日,第118887条gydF4y2Ba 10.1016 / j.fuel.2020.118887gydF4y2Ba 如果gydF4y2Ba lgydF4y2Ba 魏gydF4y2Ba J。gydF4y2Ba 习gydF4y2Ba Y。gydF4y2Ba 王gydF4y2Ba H。gydF4y2Ba 温gydF4y2Ba Z。gydF4y2Ba 李gydF4y2Ba B。gydF4y2Ba 张gydF4y2Ba H。gydF4y2Ba 长期水入侵的影响矿物和煤的孔隙结构gydF4y2Ba 燃料gydF4y2Ba 2021年gydF4y2Ba 290年,第119848条gydF4y2Ba 10.1016 / j.fuel.2020.119848gydF4y2Ba 王gydF4y2Ba F。gydF4y2Ba 他gydF4y2Ba J。gydF4y2Ba 梁gydF4y2Ba Y。gydF4y2Ba 罗gydF4y2Ba Y。gydF4y2Ba 廖gydF4y2Ba Z。gydF4y2Ba 李gydF4y2Ba lgydF4y2Ba 研究煤的渗透率特征包含煤层气不同加载路径下gydF4y2Ba 能源科学与工程gydF4y2Ba 2018年gydF4y2Ba 6gydF4y2Ba 5gydF4y2Ba 475年gydF4y2Ba 483年gydF4y2Ba 10.1002 / ese3.221gydF4y2Ba 2 - s2.0 - 85055088315gydF4y2Ba 李gydF4y2Ba B。gydF4y2Ba 雪gydF4y2Ba 年代。gydF4y2Ba 郑gydF4y2Ba C。gydF4y2Ba 提高天然气pre-drainage基于故障树系统性能失效分析对环境安全gydF4y2Ba 费森尤斯公司环境公报gydF4y2Ba 2021年gydF4y2Ba 30.gydF4y2Ba 2gydF4y2Ba 2137年gydF4y2Ba 2146年gydF4y2Ba 程gydF4y2Ba Z。gydF4y2Ba 锅gydF4y2Ba H。gydF4y2Ba 邹gydF4y2Ba Q。gydF4y2Ba 气体流动特性和优化保护性开采煤层的瓦斯抽放钻孔布置:一个案例研究从Shaqu煤矿、山西省,中国gydF4y2Ba 自然资源研究gydF4y2Ba 2021年gydF4y2Ba 30.gydF4y2Ba 2gydF4y2Ba 1481年gydF4y2Ba 1493年gydF4y2Ba 10.1007 / s11053 - 020 - 09775 - 4gydF4y2Ba 张gydF4y2Ba C。gydF4y2Ba 徐gydF4y2Ba J。gydF4y2Ba 王gydF4y2Ba E。gydF4y2Ba 彭gydF4y2Ba 年代。gydF4y2Ba 实验研究气体流动特性和卸压瓦斯抽放效果在不同卸荷应力路径gydF4y2Ba GeofluidsgydF4y2Ba 2020年gydF4y2Ba 2020年gydF4y2Ba 10gydF4y2Ba 8837962gydF4y2Ba 10.1155 / 2020/8837962gydF4y2Ba 李gydF4y2Ba T。gydF4y2Ba 吴gydF4y2Ba B。gydF4y2Ba LeigydF4y2Ba B。gydF4y2Ba 研究瓦斯抽放钻孔排水层的优化基于采矿断裂的演化特征gydF4y2Ba 能量gydF4y2Ba 2019年gydF4y2Ba 12gydF4y2Ba 23日,第4499条gydF4y2Ba 10.3390 / en12234499gydF4y2Ba 李gydF4y2Ba B。gydF4y2Ba 魏gydF4y2Ba J。gydF4y2Ba 彭gydF4y2Ba lgydF4y2Ba 数值模拟在煤层瓦斯抽放钻孔的基于气固耦合模型gydF4y2Ba 计算机模拟和新技术gydF4y2Ba 2014年gydF4y2Ba 18gydF4y2Ba 12个一个gydF4y2Ba 418年gydF4y2Ba 424年gydF4y2Ba 魏gydF4y2Ba J。gydF4y2Ba 李gydF4y2Ba B。gydF4y2Ba 王gydF4y2Ba K。gydF4y2Ba 太阳gydF4y2Ba D。gydF4y2Ba 3 d数值模拟基于pore-fracture双重瓦斯抽放钻孔的媒体gydF4y2Ba 国际矿业科技杂志》上gydF4y2Ba 2016年gydF4y2Ba 26gydF4y2Ba 4gydF4y2Ba 739年gydF4y2Ba 744年gydF4y2Ba 10.1016 / j.ijmst.2016.05.028gydF4y2Ba 2 - s2.0 - 84976519112gydF4y2Ba 李gydF4y2Ba 年代。gydF4y2Ba 相耦合模型和数值模拟煤含有天然气基于COMSOL multiphysicgydF4y2Ba 先进材料的研究gydF4y2Ba 2013年gydF4y2Ba 616 - 618gydF4y2Ba 515年gydF4y2Ba 520年gydF4y2Ba 杨gydF4y2Ba X。gydF4y2Ba 赵gydF4y2Ba X。gydF4y2Ba 研究复杂应力对断裂的影响破坏煤岩体在卸压钻孔gydF4y2Ba 在中国的科学和技术gydF4y2Ba 2015年gydF4y2Ba 9gydF4y2Ba 2gydF4y2Ba 138年gydF4y2Ba 139年gydF4y2Ba 汉gydF4y2Ba lgydF4y2Ba 基于fluid-solid-heat耦合模型的数值模拟瓦斯抽放gydF4y2Ba 辽宁科技大学学报(自然科学版)gydF4y2Ba 2013年gydF4y2Ba 32gydF4y2Ba 12gydF4y2Ba 1605年gydF4y2Ba 1608年gydF4y2Ba 李gydF4y2Ba C。gydF4y2Ba 钻井参数之间的关系,研究煤层瓦斯抽放效果gydF4y2Ba 2014年gydF4y2Ba 河南理工大学gydF4y2Ba 刘gydF4y2Ba C。gydF4y2Ba 阴gydF4y2Ba G。gydF4y2Ba 李gydF4y2Ba M。gydF4y2Ba 邓gydF4y2Ba B。gydF4y2Ba 首歌gydF4y2Ba Z。gydF4y2Ba 刘gydF4y2Ba Y。gydF4y2Ba 阴gydF4y2Ba 年代。gydF4y2Ba 页岩渗透率模型考虑床上用品在真三轴应力条件下的影响gydF4y2Ba 天然气的科学与工程》杂志上gydF4y2Ba 2019年gydF4y2Ba 68年,第102908条gydF4y2Ba 10.1016 / j.jngse.2019.102908gydF4y2Ba 2 - s2.0 - 85067177345gydF4y2Ba 刘gydF4y2Ba C。gydF4y2Ba 阴gydF4y2Ba G。gydF4y2Ba 李gydF4y2Ba M。gydF4y2Ba 商gydF4y2Ba D。gydF4y2Ba 邓gydF4y2Ba B。gydF4y2Ba 首歌gydF4y2Ba Z。gydF4y2Ba 煤岩的变形和渗透率演化考虑箱包的效果gydF4y2Ba 国际岩石力学和采矿科学杂志》上gydF4y2Ba 2019年gydF4y2Ba 117年gydF4y2Ba 49gydF4y2Ba 62年gydF4y2Ba 10.1016 / j.ijrmms.2019.03.016gydF4y2Ba 2 - s2.0 - 85063086295gydF4y2Ba 太阳gydF4y2Ba Z。gydF4y2Ba 李gydF4y2Ba lgydF4y2Ba 王gydF4y2Ba F。gydF4y2Ba 周gydF4y2Ba G。gydF4y2Ba 软、硬煤的解吸特征及其影响突出预测指标gydF4y2Ba 能源,一个部分:复苏,利用率和环境影响gydF4y2Ba 2020年gydF4y2Ba 42gydF4y2Ba 22gydF4y2Ba 2807年gydF4y2Ba 2821年gydF4y2Ba 10.1080 / 15567036.2019.1618991gydF4y2Ba 2 - s2.0 - 85066074179gydF4y2Ba 雪gydF4y2Ba Y。gydF4y2Ba 刘gydF4y2Ba J。gydF4y2Ba RanjithgydF4y2Ba p·G。gydF4y2Ba 梁gydF4y2Ba X。gydF4y2Ba 王gydF4y2Ba 年代。gydF4y2Ba 调查的影响天然气压裂的压裂特征煤炭质量和天然气开采效率基于一个基本领域模型gydF4y2Ba 石油科学与工程》杂志上gydF4y2Ba 2021年gydF4y2Ba 206年,第109018条gydF4y2Ba 10.1016 / j.petrol.2021.109018gydF4y2Ba 胡gydF4y2Ba 年代。gydF4y2Ba 道gydF4y2Ba G。gydF4y2Ba 萧红gydF4y2Ba lgydF4y2Ba 分析和数值模拟的流固耦合瓦斯抽放钻孔gydF4y2Ba 重庆大学学报gydF4y2Ba 2011年gydF4y2Ba 34gydF4y2Ba 11gydF4y2Ba 105年gydF4y2Ba 110年gydF4y2Ba 魏gydF4y2Ba 年代。gydF4y2Ba 马gydF4y2Ba K。gydF4y2Ba 李gydF4y2Ba B。gydF4y2Ba 对煤与瓦斯突出的研究基于COMSOL多重物理量gydF4y2Ba 煤炭gydF4y2Ba 2010年gydF4y2Ba 19gydF4y2Ba 5gydF4y2Ba 14gydF4y2Ba 16gydF4y2Ba 阴gydF4y2Ba G.-Z。gydF4y2Ba 明辉gydF4y2Ba l . I。gydF4y2Ba 神州gydF4y2Ba l . I。gydF4y2Ba 古文谱gydF4y2Ba l . I。gydF4y2Ba 姚gydF4y2Ba J.-W。gydF4y2Ba 张gydF4y2Ba Q.-G。gydF4y2Ba 3 d数值模拟的瓦斯抽放钻孔煤基于固气耦合模型包含气体gydF4y2Ba 中国煤炭学会杂志》上gydF4y2Ba 2003年gydF4y2Ba 38gydF4y2Ba 4gydF4y2Ba 535年gydF4y2Ba 541年gydF4y2Ba 谢gydF4y2Ba h·P。gydF4y2Ba 周gydF4y2Ba h·W。gydF4y2Ba 雪gydF4y2Ba d . J。gydF4y2Ba 高gydF4y2Ba F。gydF4y2Ba 理论、技术和工程的同时在中国开采煤炭和天然气gydF4y2Ba 中国煤炭学会杂志》上gydF4y2Ba 2014年gydF4y2Ba 39gydF4y2Ba 8gydF4y2Ba 1391年gydF4y2Ba 1397年gydF4y2Ba