煤系页岩瓦斯吸附−解吸迟滞效应核磁共振谱实验研究1)

2021-11-10唐巨鹏田虎楠潘一山

力学学报 2021年8期

唐巨鹏田虎楠,2) 潘一山

* (辽宁工程技术大学力学与工程学院,辽宁阜新 123000)

† (辽宁大学物理学院,沈阳 110036)

引言

煤系页岩气俗称瓦斯,是陆相页岩气的重要组成部分之一,主要以吸附态和游离态形式存在.页岩气单井产能并非恒定不变,而是随着累计产量增加而不断递减[1],初期单井产量高递减快,以游离气为主;后期产量稳定,递减慢,以吸附气为主[2],随着煤系页岩瓦斯开采深度不断增加,这一现象愈发明显,而煤系页岩瓦斯吸附−解吸过程差异即吸附−解吸迟滞效应在这一现象中影响不可忽略,因此通过实验和理论从微细观角度探究煤系页岩瓦斯吸附−解吸迟滞规律就显得格外重要.

煤系页岩瓦斯吸附过程中随着孔隙瓦斯压力不断增加,瓦斯(甲烷)分子在煤系页岩基质和孔喉表面不断吸附直到平衡如图1(a)所示;裂隙游离态瓦斯采出引起的裂隙和基质压力差使得基质表面吸附态瓦斯解吸[3],直到解吸平衡此时仍有部分瓦斯由于孔喉、吸附质以及纳米孔隙变形等诸多因素影响而被阻塞或封闭在微孔隙和基质中如图1(b)所示,与吸附平衡状态相比,这些被阻塞和封闭瓦斯导致明显滞后现象[4].

图1 吸附−解吸迟滞效应示意图Fig.1 Schematic diagram of adsorption-desorption hysteresis

国内外学者对页岩瓦斯迟滞效应进行大量理论和实验研究.Ekundayo 等[5-6]研究发现高温高压条件下页岩甲烷吸附−解吸存在明显滞后现象且迟滞程度与黏土和有机质碳含量有关;Zapata 和Sakhaee-Pour[7]认为页岩气吸附−解吸迟滞效应与孔隙毛细管冷凝以及吸附层厚度有关;Hazraa 等[8]研究发现贫瘠黑页岩甲烷吸附−解吸迟滞曲线形状与高岭石含量相关;Elizabeth 等[9]认为吸附−解吸过程纳米孔隙阻塞效应和孔隙弹性流体性质是产生迟滞效应原因;Mehmani 和Prodanović[10]则认为页岩气吸附−解吸滞后效应与样品微孔隙结构喉道空间分布有关;Zhao 等[11]发现甲烷吸附−解吸迟滞效应可能与页岩样品孔隙改变有关;Chen 等[12]从分子动力学角度研究发现常温条件下,页岩气吸附−解吸迟滞效应主要与孔喉尺寸变化有关,毛细管冷凝影响较小;Xu 等[13]通过采用玻尔兹曼法对页岩纳米孔气体吸附−解吸模拟发现解吸迟滞曲线形状高度依赖于孔隙形状和孔隙空间凸曲率;关富佳等[14]研究发现页岩气解吸滞后程度与孔裂隙发育和页岩物性参数有关;唐巨鹏等[15]通过核磁共振技术研究发现煤系页岩瓦斯吸附−解吸存在迟滞现象,且存在临界滞后瓦斯压力;周银波等[16]和陆壮等[17]分别从热效应和变质程度方面对煤粉吸附−解吸迟滞效应进行了定量研究.

通过国内外学者大量理论和实验研究,可以确定页岩中甲烷吸附−解吸迟滞效应的确存在,但是相关研究还存在以下问题:对页岩气吸附−解吸迟滞效应发生机理,不同学者研究结论相去甚远,且缺乏适用于页岩气吸附−解吸迟滞效应微细观定量评价指标;研究对象多为海相或海陆过渡相页岩气而对储量丰富煤系页岩瓦斯吸附−解吸迟滞效应研究甚少.

对于以上研究不足,以双鸭山盆地东保卫煤矿三采区36# 煤层底板煤系页岩为研究对象,提出了一种基于核磁共振谱理论煤系页岩瓦斯迟滞效应微细观定量评价指标;结合煤系页岩瓦斯吸附−解吸迟滞效应核磁共振谱实验进一步研究了吸附态、游离态以及微细观方法测定的宏观瓦斯迟滞规律,以期为双鸭山盆地深部煤系页岩开采提供理论基础和实验参考.

1 煤系页岩瓦斯吸附−解吸迟滞效应微细观定量评价指标

1.1 低场核磁共振谱实验原理

低场核磁共振作为一种针对含H1流体(油/气/水)等新型无损检测技术,在能源勘探领域应用日趋广泛.由低场核磁共振弛豫原理[18-19]可知横向弛豫时间T2为

式中,FS为几何形状因子,无量纲;r平均孔隙半径,μm;A为表面积,μm2;V为多孔介质体积,μm3;ρ2为表面弛豫强度,μm/ms;T2为瓦斯横向弛豫时间,ms.

低场核磁共振谱技术探测的是煤系页岩孔裂隙甲烷中氢原子核(H1),其核磁共振谱信号幅度与探测范围内的甲烷质量成正比[20],单一孔裂隙煤系页岩瓦斯质量与T2谱纵坐标幅值相对应,可以采用T2谱曲线纵坐标幅值积分表征煤系页岩孔裂隙瓦斯量.煤系页岩多孔介质结构的存在使其孔裂隙甲烷核磁共振横向弛豫时间比自由态甲烷小得多,多孔介质中不同孔径对应不同横向弛豫时间,结合自由态瓦斯核磁共振谱特性试验,可以在核磁共振T2谱定量划分微孔隙、中−大孔−裂隙以及不受限自由空间范围,进而对T2谱纵坐标积分求得吸附态和游离态瓦斯[15].

1.2 迟滞效应微细观评价指标

煤系页岩瓦斯主要以吸附态和游离态形式存在,有煤系页岩(80～100 目页岩粉)瓦斯吸附−解吸核磁共振谱实验研究发现,吸附态和游离态瓦斯解吸过程均有迟滞效应[15].国内外学者在高分子聚合物、煤岩和土壤等领域研究从宏观角度给出了吸附−解吸迟滞效应定量评价指标,如表1 所示,HI(hysteresis index)为宏观定量评价指标,但它们对于煤系页岩迟滞效应微细观评价并不完全适用.在文献[25]基础上,本文提出一种适用于煤系页岩瓦斯迟滞效应微细观定量评价指标,评价模型如图2 所示.S为煤系页岩瓦斯核磁共振T2谱幅值积分(煤系页岩瓦斯量).

图2 迟滞评价模型示意图Fig.2 Schematic diagram of hysteresis model

表1 迟滞效应宏观定量评价指标Table 1 Macro quantitative evaluation index of hysteresis

选取常用吸附−解吸模型分别拟合煤系页岩瓦斯核磁共振谱数据,采用相关系数R2作为最优等温吸附−解吸表征模型选取参数.在寻找到最优表征模型后,根据微细观迟滞评价指标MHI (microscopic hysteresis index) 定量计算煤系页岩瓦斯迟滞效应,如式(2)～式(4)所示

式中 σmin为积分起始边界压力,MPa;σmin为积分终止边界压力,MPa;Ahf为完全不可逆迟滞区域面积;Ahy为实验迟滞区域面积;fad(σ)为最优等温吸附模型;fde(σ) 为最优等温解吸模型;f(σmax)为压力为σmax最优等温吸附/解吸模型对应煤系页岩瓦斯核磁共振T2谱幅值积分;M(σmax,f(σmax))为最优等温吸附和解吸模型曲线交点.MHI准确反映了煤系页岩瓦斯微细观迟滞程度,当吸附−解吸过程完全可逆时,吸附和解吸曲线重合,实验迟滞区域面积趋于零,由式(4)可得MHI= 0;当吸附−解吸过程完全不可逆时等温解吸曲线不存在,此时Ahf=Ahy.由式(4)可得MHI= 1.

2 煤系页岩瓦斯吸附−解吸迟滞效应核磁共振谱实验

2.1 试样采集和实验设备

2.1.1 试样采集

双鸭山矿区煤系页岩瓦斯储量大且具有较高开采价值.试样取自双鸭山盆地东保卫煤矿三采区36#煤层底板煤系页岩,埋深578 m,物理参数如表2所示.

表2 页岩特征参数Table 2 Shale characteristic parameters

取未风化煤系页岩,实验室粉碎研磨制取不同粒径页岩粉,并将足量岩粉在烘干箱中105 °C 烘干24 h.冷却后将18.70 g (40～60 目)粗岩粉、37.45 g(80～110 目)细岩粉和3.15 g (30～40 目)松香,三者充分混合后在200 t 压力实验机上压制成型,试样直径为24.87 mm,高度为50.12 mm,型岩试样迅速用保鲜膜包裹,胶带缠绕后,进行蜡封.

2.1.2 实验设备

煤系页岩瓦斯吸附−解吸迟滞效应核磁共振谱现场实验如图3 所示,实验装置连接如图4 所示:主要有动力加载系统、压力调节系统、温度控制系统、MacroMR12-150H-1 低场核磁共振测试分析仪和数据采集系统,煤系页岩吸附气采用高纯度甲烷,实验温度25 ± 0.5 °C.

图3 煤系页岩瓦斯吸附−解吸迟滞效应现场实验Fig.3 Field experiment on hysteresis effect of gas adsorption-desorption in coal shale

图4 煤系页岩瓦斯吸附−解吸迟滞效应实验装置连接示意图Fig.4 Connection diagram of gas adsorption-desorption hysteresis effect experiment device for coal shale

2.2 实验步骤

为进一步研究煤系页岩瓦斯吸附−解吸迟滞效应,对同一试样首次进行三维应力状态吸附−解吸迟滞效应全过程低场核磁共振谱实验,力求能从微细观尺度对原位应力状态下煤系页岩瓦斯迟滞规律进行更为精细的实验表征.

试样和核磁共振夹持器上下压头之间存在微小空间,实验过程中此空间内会充满来自甲烷气瓶的自由态瓦斯.为了能够采用核磁共振T2谱方法计算吸附态和游离态瓦斯,有必要先对夹持器空隙自由态瓦斯进行核磁共振谱特性实验,再进行煤系页岩瓦斯吸附−解吸迟滞效应核磁共振谱实验.

2.2.1 自由态瓦斯核磁共振谱特性实验

将外径25.10 mm、内径4.2 mm、高49.60 mm的聚四氟乙烯试样放入核磁共振夹持器中,如图4所示连接好实验装置进行自由态瓦斯核磁共振谱特性实验.

首先用氦气驱除管线及核磁共振夹持器内其它干扰性吸附气体,检查管线气密性,用真空泵负压10 kPa 连续抽真空30 min,启动低场核磁共振分析仪进行数据采集.

氮气通过减压阀、压力表和恒温水浴箱进入轴压活塞腔内给试样施加轴压;启动高压氟油循环温控系统(循环氟油温度25 ± 0.5 °C),将高压氟油由核磁共振夹持器围压入口注入,围压出口循环流出给试样施加围压,同时给试样提供稳定温度场;接通高压瓦斯瓶,通过减压阀、水浴箱和压力表将瓦斯注入瓦斯参考缸中预热至25 °C,预热后瓦斯经六通阀流入核磁共振夹持器实验腔体中,给试样施加孔隙压力.自由态瓦斯核磁共振T2谱标定实验方案如表3所示,按顺序依次完成实验并采集和记录相关数据.

表3 自由态瓦斯核磁共振谱特性实验方案Table 3 Experiment of free gas by NMR spectrums

2.2.2 煤系页岩瓦斯吸附−解吸迟滞效应核磁共振谱实验

页岩埋深578 m,综合考虑实验安全和非磁性核磁共振夹持器极限承载能力,在这里采用地应力折减、瓦斯压力不折减的方法对模型施加地应力和瓦斯压力[25],选取煤系页岩瓦斯压力比例尺Cρ= 1,垂直地应力比例尺Cσ= 4,侧压系数为1.2,地层平均容重为25 kN/m3[26],可得实验工况轴压为3.61 MPa,围压为2.97 MPa,实验甲烷压力为2.24 MPa.同时对相邻和深部煤系页岩瓦斯进行实验室相似工况模拟,模拟埋深分别为350 m,550 m,750 m,850 m,1050 m,1200 m 和1500 m.自由态瓦斯标定实验结束后,将直径24.87 mm、长度50.12 mm 型页岩试样放入核磁共振夹持器中,按图4 所示重新连接实验装置,进行煤系页岩瓦斯吸附−解吸迟滞效应核磁共振谱实验.

(1)吸附实验

首先用氦气对管线和夹持器腔体进行冲扫,负压10 kPa 连续抽真空120 min.启动低场核磁共振分析仪进行数据采集,压力加载顺序为轴压−围压−孔隙压力,压力加载方式和自由态瓦斯标定实验相同,具体实验方案如表4 所示,按顺序依次完成实验并采集和记录相关数据.

表4 吸附过程核磁共振谱实验方案Table 4 Experimental adsorption of gas by NMR spectrums

(2) 解吸实验

在煤系页岩瓦斯吸附实验结束后,立即进行煤系页岩瓦斯解吸实验,具体实验方案如表5 所示,按顺序依次完成实验并采集和记录相关实验数据.

表5 解吸过程核磁共振谱实验方案Table 5 Experimental desorption of gas by NMR spectrums

3 实验结果

煤系页岩瓦斯吸附−解吸实验过程中,核磁共振夹持器内测试样品同时受到轴压−围压−孔隙压力共同作用,采用平均有效应力σe[27]表示三者对试样应力场综合作用效果,相关实验结果如下.

3.1 自由态瓦斯核磁共振谱特性实验结果

自由态瓦斯核磁共振谱特性实验结果如图5 所示,随着平均有效应力增加,自由态瓦斯横向弛豫时间T2逐渐增大,当平均有效应力在1.76～8.98 MPa,对应自由态瓦斯T2谱范围为93.00～2710.63 ms,且T2谱曲线只有一个特征峰,这与前人研究结论一致[28].

图5 自由态瓦斯T2 谱Fig.5 T2 spectrum of free gas

3.2 煤系页岩瓦斯吸附−解吸迟滞效应核磁共振谱实验结果

煤系页岩瓦斯迟滞效应核磁共振谱实验结果如图6 和图7 所示.

图6 吸附过程T2 谱Fig.6 T2 spectrum of shale gas adsorption

图7 解吸过程T2 谱Fig.7 T2 spectrum of shale gas desorption

由图5 可以确定图6 中横向弛豫时间T2大于103.69 ms 部分为自由态瓦斯谱峰(>T2C2= 103.69 ms).由式(1)可知孔隙越小T2时间越短,因此微孔隙横向弛豫时间要比中−大孔隙−裂隙横向弛豫时间短,所以图6 中横向弛豫时间T2小于2.86 ms 部分为吸附态瓦斯谱峰(T2C4= 115.61 ms)、吸附态瓦斯T2谱峰(

4 煤系页岩瓦斯吸附−解吸迟滞效应分析

由图6 和图7 可得吸附态和游离态瓦斯迟滞效应评价模型如图8 和图9 所示.

图8 吸附态瓦斯迟滞效应评价模型Fig.8 Evaluation model for hysteresis effect of adsorbed gas

图9 游离态瓦斯迟滞效应评价模型Fig.9 Evaluation model for hysteresis effect of porous mediumconfined gas

4.1 吸附态瓦斯迟滞效应分析

4.1.1 吸附态瓦斯迟滞现象

由图8 可知吸附态瓦斯量(吸附态瓦斯T2谱幅值积分S)与平均有效应力 σe之间具有明显非线性特征,使用D-A,D-R 和BET 等温模型对吸附态瓦斯吸附过程模拟时均需饱和蒸汽压力P0,而实验温度25 ± 0.5 °C(甲烷临界温度−82.6 °C),因此采用虚拟饱和蒸汽压力替代饱和蒸汽压力,采用Dubinin 算法[29]可得实验条件虚拟饱和蒸汽压力为11.24 MPa.

选取具有代表性的吸附和解吸模型,分别模拟吸附态瓦斯吸附−解吸核磁共振谱实验数据,模型公式及其相关系数R2如表6 和表7 所示.由表6 可得采用D-R 模型拟合吸附态瓦斯吸附实验数据,对应相关系数R2值为0.99475 拟合程度最高;由表7 可得采用Weibull 函数模型拟合吸附态瓦斯解吸实验数据,对应相关系数R2值为0.996 6 拟合程度最高;因此选用D-R 模型和Weibull 函数模型分别作为吸附态瓦斯吸附和解吸过程最佳表征模型.

表6 常用吸附模型Table 6 Common adsorption models

表7 常用解吸模型Table 7 Common desorption models

由图8 可得吸附过程拟合曲线

解吸过程拟合曲线

吸附过程S1(T2)与解吸过程S2(T2)在高应力阶段相交于M3点,且存在临界滞后压力7.70 MPa.当平均有效应力大于7.70 MPa,吸附曲线S1(T2)比解吸曲线S2(T2)高,无明显迟滞现象;当平均有效应力在1.88～7.70 MPa,解吸曲线S2(T2) 在吸附曲线S1(T2)之上,解吸过程具有明显滞后环.

4.1.2 吸附态瓦斯迟滞区域面积

由图8 可知,吸附过程发生在平均有效应力1.88～8.94 MPa 范围内,而解吸过程平均有效应力为1.10～8.94 MPa,因此以 σe1= 1.88 MPa 作为吸附态瓦斯迟滞起始边界,并将其代入解吸曲线S2(T2) =求得S2(T2) =175.10,即解吸实验曲线模拟终点M1(1.88,175.10);吸附−解吸曲线交点横坐标为7.70 MPa,对应吸附态瓦斯幅值积分为2133.91 即M3(7.70,2133.91);吸附−解吸迟滞区域起始边界和上边界交点为M2(1.88,2133.91).在此边界条件结合式(2)～式(3)可得吸附态瓦斯完全不可逆迟滞区域面积Ahf−ad和实验迟滞区域面积Ahy−ad如式(5)～式(6)所示, σe−min—积分起始边界压力,MPa

4.2 游离态瓦斯迟滞效应分析

4.2.1 游离态瓦斯吸附−解吸迟滞现象

由图9 可得游离态瓦斯量(游离态瓦斯T2谱幅值积分S)与平均有效应力 σe之间具有明显线性特征.设游离态瓦斯吸附和解吸最佳表征模型均为S=Kσe+C,K为斜率,C为截距.

吸附过程拟合曲线

解吸过程拟合曲线

当平均有效应力大于8.47 MPa,游离态瓦斯解吸曲线S4(T2)低于吸附曲线S3(T2),无明显迟滞现象;当平均有效应力在1.88～8.47 MPa 范围内,游离态瓦斯解吸曲线S4(T2)高于吸附曲线S3(T2),游离态瓦斯解吸过程具有明显迟滞效应,且存在临界迟滞压力8.47 MPa.

4.2.2 游离态瓦斯迟滞区域面积

由图9 可知游离态瓦斯吸附过程发生在平均有效应力1.88～8.94 MPa 范围内,而解吸过程平均有效应力为1.10～8.94 MPa,因此以 σe2= 1.88 MPa 作为游离态瓦斯迟滞起始边界,并将其代入解吸曲线S4(T2) = 39.472σe+ 240.124,求得S4(T2) = 314.33,即得游离态瓦斯解吸实验曲线模拟终点M1(1.88,314.33);游离态瓦斯吸附−解吸曲线交点横坐标为8.47 MPa,对应幅值积分为574.48 即M3(8.47,574.48);在此边界条件结合式(2)～式(3)可得游离态瓦斯完全不可逆迟滞区域面积Ahf−pm和实验迟滞区域面积Ahy−pm如式(7)～式(8)所示, σe−min—积分起始边界压力,MPa

4.3 煤系页岩瓦斯微细观迟滞规律

煤系页岩瓦斯主要以吸附态和游离态形式存在(忽略极少量基质固溶态瓦斯),因此可以用吸附态和游离态瓦斯幅值积分之和,近似表征吸附−解吸过程中煤系页岩瓦斯宏观吸附量和解吸量,从微细观角度探究煤系页岩瓦斯宏观吸附−解吸迟滞规律.结合式(2)～式(8)可得煤系页岩瓦斯宏观迟滞系数MHImac,以及吸附态和游离态瓦斯迟滞系数MHIad−mac和MHIpm−mac如式(9)～式(11)所示,并由此可得瓦斯迟滞系数与平均有效应力关系,如图10所示

由图10 可得煤系页岩瓦斯宏观迟滞系数与平均有效应力具有明显幂函数关系.宏观迟滞效应中由吸附态或游离态瓦斯引起的迟滞系数与平均有效应力均可采用二次多项式拟合,且两曲线交于点M(6.30,15.13).

图10 瓦斯迟滞系数与平均有效应力关系Fig.10 Relationship between hysteresis index of gas and average effective stress

当平均有效应力 σe在高应力状态(6.30～7.14 MPa)游离态瓦斯迟滞系数曲线(b)在吸附态瓦斯迟滞系数曲线(c)之上且两曲线变化趋势相反,游离态瓦斯对宏观瓦斯迟滞效应贡献大于吸附态瓦斯;当平均有效应力 σe在(1.88～6.30 MPa)游离态瓦斯迟滞系数曲线(b)在吸附态瓦斯迟滞系数曲线(c)之下,此阶段宏观瓦斯迟滞效应受控于吸附态瓦斯迟滞效应.

4.3.1 煤系页岩瓦斯宏观迟滞规律

由核磁共振谱方法测定的煤系页岩瓦斯宏观迟滞系数MHImac= 0.23σe2.021+ 20.407,R2= 0.97708,随着平均有效应力减小,瓦斯宏观迟滞系数呈幂函数减小趋势如图10 中(a)曲线所示,当平均有效应力从7.14 MPa 降低到1.88 MPa,对应瓦斯宏观迟滞系数由34.02%降低到20.54%,瓦斯宏观迟滞效应逐渐减小.随着宏观迟滞系数减小,煤系页岩瓦斯解吸量逐渐增大.

4.3.2 吸附态瓦斯迟滞规律

宏观迟滞效应中由吸附态瓦斯引起的迟滞系数MHIad−mac= −0.664σ2e+ 6.038σe+ 3.481,R2=0.82395.

随着平均有效应力降低,吸附态瓦斯迟滞系数呈现出先增大后减小趋势如图10 中(b)曲线所示.当平均有效应力从7.14 MPa 降低到4.55 MPa,对应迟滞系数反而从11.35%增加到17.21%,然而当平均有效应力从4.55 MPa 降低到1.88 MPa,对应迟滞系数则由17.21%逐渐降低到13.12%.

这一变化规律与煤系页岩原生微孔隙重新吸附和新生孔隙吸附有关.煤系页岩微孔隙吸附态瓦斯处于吸附−解吸可逆动平衡状态,卸压初始阶段(4.55～7.14 MPa)虽然平均有效应力有所降低,但此时煤系页岩仍处于高应力状态,当平均有效应力大于部分微孔隙端部应力强度因子,微孔隙裂纹会继续扩展,形成新微孔隙[30],同时新生微孔隙也会使吸附态瓦斯扩散通道迂曲度变大,吸附态瓦斯向外扩散受限.此时解吸出的吸附态瓦斯一部分在高瓦斯浓度条件下重新在原生微孔隙表面吸附,一部分则会在新生微孔隙吸附,高平均有效应力阶段,虽然平均有效应力降低36.27%,但由吸附态瓦斯引起宏观迟滞系数反而增加51.63%.

随着继续卸压(1.88～4.55 MPa),平均有效应力由4.55 MPa 逐步减小,新微孔隙产生速度减慢,孔隙变形损伤减小.同时随着游离态瓦斯不断采出,孔裂隙瓦斯浓度梯度逐渐变大,在瓦斯浓度梯度作用下吸附态瓦斯加速向中−大孔隙−裂隙结构扩散,吸附态瓦斯引起的宏观迟滞效应逐渐变小,吸附态瓦斯开始迅速大量解吸.

4.3.3 游离态瓦斯迟滞系数

宏观迟滞效应中由游离态瓦斯引起的迟滞系数MHIpm−mac= 0.895σ2e− 5.856σe+ 16.499,R2=0.91806.

随着平均有效应力降低,迟滞系数呈现先减小后增大趋势如图10 中(c)曲线所示.当平均有效应力从7.14 MPa 降低到3.27 MPa,对应迟滞系数由22.67%减小到6.92%,但是当平均有效应力从3.27 MPa 继续降低到1.88 MPa,对应迟滞系数反而从6.92%增加到8.65%.

煤系页岩初始卸压阶段(3.27～7.14 MPa),随着平均有效应力逐步减小中−大孔隙−裂隙损伤逐渐变大,游离态瓦斯散失通道迂曲度增加.但此时游离态瓦斯流动近似为达西流动且大部分处于超临界状态,其渗透性远超出常温常压状态[31],即使孔隙裂隙迂曲度增加游离态瓦斯仍然快速向井壁渗透,游离态瓦斯引起的迟滞系数急速下降.

继续卸压进入低应力阶段(1.88～3.27 MPa),中−大孔隙−微裂隙部分继续闭合,游离态瓦斯向井壁散失通道进一步受限,且此时游离态瓦斯压力在其临界压力4.59 MPa 以下,渗透性明显减弱,游离态瓦斯渗透能力降低,迟滞系数随着平均有效应力减小反而增加25.04%.

5 迟滞效应对煤系页岩瓦斯开采影响

随着煤系页岩瓦斯开采深度不断增加,受高渗透压、高地应力和高地温综合影响,煤岩体中瓦斯高度聚集,瓦斯压力急剧增大[32],煤系页岩瓦斯深部原位赋存环境变得更为复杂.由本文研究结论可知:随着埋深增加平均有效应力逐渐变大,煤系页岩瓦斯吸附−解吸迟滞系数呈幂函数增加趋势.但深部煤系页岩在高地应力和高孔裂隙瓦斯压力作用下,内部孔裂隙损伤会进一步加大,进而加剧微孔隙吸附态瓦斯扩散受限和中−大孔隙−裂隙游离态瓦斯渗透阻力变大,渗流路径迂曲度增加[33],而这些会进一步加剧深部煤系页岩瓦斯迟滞效应,使深部煤系页岩瓦斯解吸−扩散−渗流机理变得更为复杂.因此建议使用三维应力状态煤系页岩瓦斯等温解吸模型去评估深部煤系页岩瓦斯可采储量,替代常用一维(仅考虑瓦斯压力)等温吸附模型.

深部煤系页岩储层孔隙度和渗透率都很低,同时受煤系页岩瓦斯吸附−解吸迟滞效应影响,煤系页岩瓦斯自然采收率普遍偏低,因此煤系页岩瓦斯商业化开采必须依赖于有效的水力压裂改造措施.由上述研究可知煤系页岩瓦斯宏观迟滞效应中,由吸附态或游离态瓦斯引起的迟滞系数与平均有效应力均具有明显二次多项式函数关系,因此在煤系页岩瓦斯水力压裂开采过程中,应该科学把控支撑剂充填时间节点和不同粒度充填剂合理配比,同时尝试在高应力卸压阶段向微孔隙结构中输送和铺置粒径更小,强度更高支撑剂,以增加对煤系页岩微孔隙结构改造,进一步减小吸附−解吸过程中由吸附态瓦斯引起的迟滞效应,加速吸附态瓦斯解吸扩散,提高煤系页岩瓦斯产量.