CN114894671B - 基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法 - Google Patents
基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法 Download PDFInfo
- Publication number
- CN114894671B CN114894671B CN202210312854.8A CN202210312854A CN114894671B CN 114894671 B CN114894671 B CN 114894671B CN 202210312854 A CN202210312854 A CN 202210312854A CN 114894671 B CN114894671 B CN 114894671B
- Authority
- CN
- China
- Prior art keywords
- adsorption
- gas
- coal
- desorption
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N7/00—Analysing materials by measuring the pressure or volume of a gas or vapour
- G01N7/02—Analysing materials by measuring the pressure or volume of a gas or vapour by absorption, adsorption, or combustion of components and measurement of the change in pressure or volume of the remainder
- G01N7/04—Analysing materials by measuring the pressure or volume of a gas or vapour by absorption, adsorption, or combustion of components and measurement of the change in pressure or volume of the remainder by absorption or adsorption alone
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N2013/003—Diffusion; diffusivity between liquids
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Combustion & Propulsion (AREA)
- Engineering & Computer Science (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
本发明公开了基于扩散‑渗流解耦实验的煤体甲烷吸附时间测定方法,包括以下步骤:(1)选取块煤样品,分析获得煤样的灰分、水分,并测定煤样的吸附常数的a、b值;(2)对步骤(1)选取的块煤样品进行吸附解吸实验;(3)对步骤(2)实验后的煤样进行渗流通道改造;(4)对改造完毕后的样品进行吸附解吸实验;(5)对步骤(4)实验完成后的煤样继续进行渗流通道改造,继续进行吸附解吸实验;(6)分析实验后获得的多组吸附解吸数据,确定差异消失时间t1,由解吸速率计算t1时刻窜流量;(7)通过双重孔隙介质拟稳态窜流方程计算吸附时间;能真实反映煤体在煤层中解耦扩散行为和渗流行为。
Description
技术领域
本发明涉及煤体瓦斯运移特性研究技术领域,具体涉及基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法。
背景技术
瓦斯是一种赋存在煤层中的清洁能源。瓦斯抽采不仅可以提高煤炭开采回采率,还可以预防煤与瓦斯突出事故的发生。在煤炭瓦斯抽采过程中,瓦斯要经过脱附-扩散-渗流等环节的耦合过程才能被抽采出来。瓦斯的扩散特性是影响瓦斯解吸速率的重要因素,也直接影响瓦斯的抽采效果。
吸附时间是一个反映了瓦斯在煤岩基质中的扩散能力的重要参数,体现了扩散系数和形状因子对扩散通量的控制作用。在中国《煤层气测试方法》标准中(GB/T 19559-2008),吸附时间是指样品所含气体被解吸出来63.2%所用的时间。然而,使用块煤测试吸附时间时,由于测得的所含气体经脱附-扩散-渗流耦合过程,获得的扩散相关参数受渗流影响。为准确测量吸附时间等扩散参数、缺少一种扩散-渗流解耦实验方法,用于单独测试扩散行为或渗流行为。
发明内容
针对上述存在的技术不足,本发明的目的是提供基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法,本发明的测试方法操作简单,且能有效避免块煤样品的渗流特性对扩散特性测试导致的误差,能真实反映煤体在煤层中解耦扩散行为和渗流行为。
为解决上述问题,本发明采用如下技术方案:
本发明提供基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法,包括以下步骤:
(1)选取块煤样品,记录煤样重量,对煤样进行工业分析获得煤样的灰分、水分,测试煤样密度,并采用高压容量法测定煤样的吸附常数的a、b值;
(2)抽真空后,对步骤(1)选取的块煤样品进行吸附解吸实验,实验过程中用压力传感器实时记录压力变化,实验时间持续至吸附解吸平衡为止;
(3)通过中心钻孔的方式对步骤(2)实验后的煤样进行渗流通道改造;
(4)对改造完毕后的样品进行吸附解吸实验,吸附解吸条件与步骤(2)相同,实验过程中用压力传感器实时记录压力变化,实验时间持续至吸附解吸平衡为止;
(5)对步骤(4)实验完成后的煤样继续进行渗流通道改造,扩大煤样的中心钻孔,对扩大了中心钻孔的煤样进行吸附解吸实验,吸附解吸条件与步骤(2)相同,实验过程中用压力传感器实时记录压力变化,实验时间持续至吸附解吸平衡为止;
(6)分析步骤(2)、步骤(4)、步骤(5)实验后获得的多组吸附解吸数据,根据吸附速率曲线的差异确定差异消失时间t1,从而确定扩散-渗流行为的主控阶段,在解耦扩散-渗流行为后,由解吸速率以及步骤(1)测得的煤样密度计算t1时刻窜流量;
(7)将步骤(1)中得到的吸附常数a、b值和通过步骤(2)、步骤(4)、步骤(5)实验后获得的块煤在步骤(6)确定的t1时刻的瓦斯吸附量带入朗格缪尔方程,通过朗格缪尔方程反演计算不同渗流通道对应的t1时刻的基质中的气相瓦斯压力;将得到的基质中的气相瓦斯压力以及步骤(6)求得的t1时刻的窜流量代入到双重孔隙介质拟稳态窜流方程计算吸附时间。
优选地,步骤(2)、步骤(4)、步骤(5)实验后获得的多组吸附解吸数据包括瓦斯吸附量、吸附速率、瓦斯解吸量、解吸速率游离瓦斯含量。
优选地,步骤(6)中,所述窜流量采用以下公式计算:
式中:Qs为“窜流量”,即单位体积煤基质与裂隙系统的质量交换率,kg/(m3·s);为t时刻瓦斯解吸率,mL/(g·min);/>为标准状态下甲烷的密度,0.71428kg/m3;ρcoal为煤样密度,g/mL;当t=t1时,式(1-1)求得t1时刻的窜流量。
优选地,步骤(7)中,所述的双重孔隙介质拟稳态窜流方程如下所示:
式中:Qs为“窜流量”,即单位体积煤基质与裂隙系统的质量交换率,kg/(m3·s);D为瓦斯在基质孔隙中的扩散系数,m2/s;σc为基质形状因子,m-2;Cm、Cf分别为基质、裂隙中的气相瓦斯浓度,kg/m3;Pm为基质裂隙中的气相瓦斯压力,Pa;Pf为裂隙中的气相瓦斯压力,Pa;M为甲烷的摩尔质量;R为理想气体常数,J/(mol·K);T为温度,K;τ为吸附时间。
优选地,所述基质中的气相瓦斯浓度Cm、裂隙中的气相瓦斯浓度Cf分别采用以下公式计算:
式中:M为甲烷的摩尔质量,kg/mol;R为理想气体常数,J/(mol·K);T为温度,K;Pm、Pf分别为基质、裂隙中的气相瓦斯压力,Pa。
优选地,所述朗格缪尔方程如下所示:
式中:V为瓦斯吸附量,m3/t;P为吸附平衡时瓦斯压力,MPa;a为吸附常数,m3/t;b为吸附常数,MPa-1。
将步骤(1)中得到的吸附常数a、b值和通过步骤(2)、步骤(4)、步骤(5)得到的在步骤(6)确定的t1时刻的块煤基质瓦斯吸附量带入朗格缪尔方程,采用如下公式计算t1时刻的基质中气相瓦斯压力Pm:
式中:Pm为基质中气相瓦斯压力,Pa;Vm为基质中的瓦斯吸附量,mL/g。优选地,所述基质中的瓦斯吸附量Vm采用如下公式计算:
Vm=Vtotal-Vde-Vfree (0-11)
式中:Vtotal——为解吸实验开始前煤样罐内吸附量,mL/g;
Vde——为t时刻瓦斯解吸量,mL/g;
Vfree——为常压下游离瓦斯含量,mL/g。
本发明的有益效果在于:
本发明基于扩散-渗流解耦实验,可以准确测试窜流量并计算吸附时间,有利于独立评价扩散行为和渗流行为,能够为煤层瓦斯运移研究提供参考。
本发明的测试方法操作简单,且能有效避免块煤样品的渗流特性对扩散特性测试导致的误差,准确测试吸附时间等扩散特性参数,真实反映煤体在原位地应力条件下的扩散行为和渗流行为,对进一步推动煤层瓦斯运移的相关研究有较好指导意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法应用流程图;
图2为实验样品示意图,图2(a)为原始样品,图2(b)为第一次渗流通道改造后的样品,图2(c)为第二次渗流通道改造后的样品;
图3为样品解吸速率图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1至图3所示,本实施例提供一种基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法,具体步骤如下:
(1)选取块煤样品,记录煤样重量,根据《煤的工业分析测试方法》(GBT/T 212-2008)对煤样进行工业分析获得煤的灰分、水分,测试煤样密度,并根据《煤的甲烷吸附量测定方法(高压容量法)》(MT/T 752-1997)获得对煤进行吸附常数测试吸附常数a、b值;
(2)抽真空后,对步骤(1)选取的块煤样品进行吸附解吸实验,实验过程中用压力传感器实时记录压力变化,实验时间持续至吸附解吸平衡为止;
(3)步骤(2)实验结束后取出经过吸附解吸实验的块煤样品,对块煤样品中心钻孔以改变样品渗流通道,所得实验样品如图2(b)所示;
(4)改造完毕后,对改造好的样品进行吸附解吸实验,实验中吸附解吸条件与未改造前即步骤(2)的实验条件相同,实验过程中用压力传感器实时记录压力变化,实验时间持续至吸附解吸平衡为止;;
(5)继续对步骤(4)实验后的煤块样品进行渗流通道改造,所得实验样品如图2(c)所示,改造完毕后,对改造好的样品进行吸附解吸实验,实验中保证吸附解吸条件与步骤(2)的实验条件相同,实验过程中用压力传感器实时记录压力变化,实验时间持续至吸附解吸平衡为止;
(6)分析步骤(2)、(4)、(5)中获得的样品解吸数据如图3所示,可以确定扩散-渗流行为的主控阶段:在0-t1阶段为渗流控制阶段,在t1-t阶段为扩散控制阶段;在扩散控制阶段,由t1时刻解吸速率可以计算此时的窜流量Qs;
所述的窜流量Qs计算公式如下所示:
式中Qs——为“窜流量”,即单位体积煤基质与裂隙系统的质量交换率,kg/(m3·s);——为t时刻瓦斯解吸率,mL/(g·min);/>——为标准状态下甲烷的密度,0.71428kg/m3;ρcoal——为煤样密度,g/mL;
式中:V为瓦斯吸附量,m3/t;P为吸附平衡时瓦斯压力,MPa;a为吸附常数,m3/t;b为吸附常数,MPa-1;
以及通过步骤(2)、(4)、(5)的实验获得的t1时刻的煤柱瓦斯吸附量Vm反演计算t1时刻的基质的气相瓦斯压力Pm;根据双重孔隙介质拟稳态窜流方程,求得吸附时间τ;
所述的双重孔隙介质拟稳态窜流方程如下所示:
式中:D——为瓦斯在基质孔隙中的扩散系数,m2/s;σc——为基质形状因子,m-2;Cm、Cf——分别为基质、裂隙中的气相瓦斯浓度,kg/m3,一般定义为:
式(0-14)和式(0-15)中M——为甲烷的摩尔质量,kg/mol;R——为理想气体常数,J/(mol·K);T——为温度,K;Pm、Pf——分别为基质、裂隙中的气相瓦斯压力,Pa;
所述的基质孔隙压力Pm计算公式如下所示:
式中:Vm为基质中的瓦斯吸附量,mL/g。
Vm=Vtotal-Vde-Vfree (0-17)
式中:Vtotal为解吸实验开始前煤样罐内吸附量,mL/g;
Vde为t时刻瓦斯解吸量,mL/g;
Vfree为常压下游离瓦斯含量,mL/g;
在扩散控制阶段,裂隙中的瓦斯几近释放完毕,该阶段的解吸量主要来自于从基质扩散出来的瓦斯,裂隙中气相瓦斯压力Pf等于大气压力P0。将式(1-5)带入式(1-2)可知,所述的计算吸附时间计算公式如下所述:
表1块煤甲烷吸附时间计算结果
本实施例的方法操作简单,且能够通过吸附解吸数据分析有效避免块煤样品的渗流特性对扩散特性测试导致的误差,准确测试吸附时间等扩散特性参数,能够真实反映煤体在原位地应力条件下的扩散行为和渗流行为,对进一步推动煤层瓦斯运移的相关研究有较好指导意义。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (6)
1.基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法,其特征在于,包括以下步骤:
(1)选取块煤样品,记录煤样重量,对煤样进行工业分析获得煤样的灰分、水分,测试煤样密度,并采用高压容量法测定煤样的吸附常数a、b值;
(2)抽真空后,对步骤(1)选取的块煤样品进行吸附解吸实验,实验过程中用压力传感器实时记录压力变化,实验时间持续至吸附解吸平衡为止;
(3)通过中心钻孔的方式对步骤(2)实验后的煤样进行渗流通道改造;
(4)对改造完毕后的样品进行吸附解吸实验,吸附解吸条件与步骤(2)相同,实验过程中用压力传感器实时记录压力变化,实验时间持续至吸附解吸平衡为止;
(5)对步骤(4)实验完成后的煤样继续进行渗流通道改造,扩大煤样的中心钻孔,对扩大了中心钻孔的煤样进行吸附解吸实验,吸附解吸条件与步骤(2)相同,实验过程中用压力传感器实时记录压力变化,实验时间持续至吸附解吸平衡为止;
(6)分析步骤(2)、步骤(4)、步骤(5)实验后获得的多组吸附解吸数据,根据吸附速率曲线的差异确定差异消失时间t1,从而确定扩散-渗流行为的主控阶段,在解耦扩散-渗流行为后,由解吸速率以及步骤(1)测得的煤样密度计算t1时刻窜流量;
所述窜流量采用以下公式计算:
式中:Qs为“窜流量”,即单位体积煤基质与裂隙系统的质量交换率,kg/(m3·s);为t时刻瓦斯解吸率,mL/(g·min);/>为标准状态下甲烷的密度,0.71428kg/m3;ρcoal为煤样密度,g/mL;
(7)将步骤(1)中得到的吸附常数a、b值和通过步骤(2)、步骤(4)、步骤(5)实验后获得的块煤在步骤(6)确定的t1时刻的瓦斯吸附量带入朗格缪尔方程,通过朗格缪尔方程反演计算不同渗流通道对应的t1时刻的基质中的气相瓦斯压力;将得到的基质中的气相瓦斯压力以及步骤(6)求得的t1时刻的窜流量代入到双重孔隙介质拟稳态窜流方程计算吸附时间。
2.如权利要求1所述的基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法,其特征在于,步骤(2)、步骤(4)、步骤(5)实验后获得的多组吸附解吸数据包括瓦斯吸附量、吸附速率、瓦斯解吸量、解吸速率游离瓦斯含量。
6.如权利要求5所述的基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法,其特征在于,所述基质中的瓦斯吸附量Vm采用如下公式计算:
Vm=Vtotal-Vde-Vfree (0-6)
式中:Vtotal——为解吸实验开始前煤样罐内吸附量,mL/g;
Vde——为t1时刻瓦斯解吸量,mL/g;
Vfree——为常压下游离瓦斯含量,mL/g。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210312854.8A CN114894671B (zh) | 2022-03-28 | 2022-03-28 | 基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210312854.8A CN114894671B (zh) | 2022-03-28 | 2022-03-28 | 基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114894671A CN114894671A (zh) | 2022-08-12 |
CN114894671B true CN114894671B (zh) | 2023-06-02 |
Family
ID=82714777
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210312854.8A Active CN114894671B (zh) | 2022-03-28 | 2022-03-28 | 基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114894671B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115828784B (zh) * | 2022-12-07 | 2023-09-26 | 重庆中环建设有限公司 | 穿越煤系地层隧道施工瓦斯涌出特征的预测方法及系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102830213B (zh) * | 2012-08-10 | 2015-07-29 | 河南理工大学 | 变温条件下受载含瓦斯煤吸附-解吸-渗流实验系统 |
KR101782294B1 (ko) * | 2017-04-20 | 2017-10-23 | 한국가스공사 | 탄층 메탄가스 저류층의 가스 생산 모델링방법 |
CN108458951B (zh) * | 2018-04-02 | 2020-06-26 | 河南理工大学 | 通过解吸率快速确定煤层瓦斯含量的方法 |
CN109142669A (zh) * | 2018-09-27 | 2019-01-04 | 西安科技大学 | 一种基于数据相关的煤层瓦斯参数快速测定方法 |
CN111337388A (zh) * | 2018-12-19 | 2020-06-26 | 河南理工大学 | 一种基于扩散率曲线差减法快速测定煤层瓦斯含量方法 |
CN113899654A (zh) * | 2021-10-11 | 2022-01-07 | 重庆大学 | 一种煤样瓦斯逸散损失量测定的实验系统和实验方法 |
-
2022
- 2022-03-28 CN CN202210312854.8A patent/CN114894671B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114894671A (zh) | 2022-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106525690B (zh) | 致密砂岩稳态法气水相对渗透率曲线测定方法 | |
CN111257202A (zh) | 一种含吸附气条件下页岩压裂液强制渗吸及返排实验方法 | |
CN103792334B (zh) | 一种预测缝洞型碳酸盐岩酸蚀蚓孔宽度的方法 | |
Harpalani et al. | Impact of CO 2 injection on flow behavior of coalbed methane reservoirs | |
Harpalani et al. | Influence of gas production induced volumetric strain on permeability of coal | |
CN108458951B (zh) | 通过解吸率快速确定煤层瓦斯含量的方法 | |
US7406857B2 (en) | Electronic humidity chamber for vapor desorption to determine high capillary pressures | |
CN104196524B (zh) | 一种欠饱和煤储层开发的气水产出动态相渗曲线测量方法 | |
CN105910971A (zh) | 富有机质致密岩心气体渗透率和扩散系数的联测方法 | |
WO2018236637A1 (en) | METHOD FOR CORRECTING LOW PERMEABILITY LABORATORY MEASUREMENTS FOR LEAKS | |
CN114894671B (zh) | 基于扩散-渗流解耦实验的煤体甲烷吸附时间测定方法 | |
CN103226089B (zh) | 一种页岩气体渗透率测定方法 | |
CN113075108B (zh) | 一种考虑束缚水饱和度的岩心多次应力敏感测试方法 | |
Wang et al. | Seepage law and permeability calculation of coal gas based on Klinkenberg effect | |
CN112878987B (zh) | 利用生产数据计算页岩气井控制储量的方法 | |
CN107917868A (zh) | 一种围压下页岩自吸能力的测试装置及其测试方法 | |
CN104268377A (zh) | 一种基于煤岩工业组分的脆性指数确定方法 | |
CN112924354A (zh) | 一种基于气体扩散系数的页岩储层应力敏感实验评价方法 | |
CN111337388A (zh) | 一种基于扩散率曲线差减法快速测定煤层瓦斯含量方法 | |
CN112432889A (zh) | 一种液体饱和法测岩石孔隙度的方法及校正方法 | |
Xu et al. | CO 2 storage capacity of anthracite coal in deep burial depth conditions and its potential uncertainty analysis: a case study of the No. 3 coal seam in the Zhengzhuang Block in Qinshui Basin, China | |
CN114720655A (zh) | 同时测量岩心不同赋存状态气体产出特征的系统及方法 | |
CN110715879B (zh) | 基于气水分布的高演化页岩储层微孔隙吸附气量评价方法 | |
CN108593493A (zh) | 一种确定致密储层含气量的方法和装置 | |
CN111579584A (zh) | 一种基于红外热成像技术识别煤中气体富集区域的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |