岩性油气藏  2021, Vol. 33 Issue (2): 171-179       PDF    
×
基于逸度与压力计算页岩吸附甲烷的等量吸附热差异分析——以延长探区延长组页岩为例
薛培1,2, 张丽霞2, 梁全胜2, 师毅3    
1. 西北政法大学, 西安 710063;
2. 陕西延长石油 (集团) 有限责任公司研究院, 西安 710075;
3. 延长石油集团油气勘探公司, 陕西 延安 716000
摘要: 为了完善吸附热力学参数的计算方法,明确陆相页岩吸附CH4的热力学特征,以延长探区延长组页岩为研究对象,开展了不同温度下页岩吸附CH4的等温吸附实验,并利用绝对吸附量曲线对比分析了基于逸度-绝对吸附量与压力-绝对吸附量的等量吸附热的差异。结果表明:①逸度小于压力,在低压力区间(0.36~2.21 MPa)内,逸度与压力的差值较小,随着温度降低和压力升高,逸度与压力的差值增加。②基于逸度与压力获得的等量吸附热曲线均具有明显的单调递增的线性特征,但采用压力-绝对吸附量获得的等量吸附热曲线斜率大于采用逸度-绝对吸附量获得的等量吸附热曲线,这表明基于压力-绝对吸附量获得的等量吸附热计算结果中吸附质分子间作用力的影响过大。③绝对吸附量为0.01~0.35 mmol/g,基于逸度-绝对吸附量与压力-绝对吸附量获得的等量吸附热的相对偏差为-86.54%~57.01%,页岩气吸附体系的热力学评价应当采用逸度-绝对吸附量数据作为等量吸附热计算的基础数据。
关键词: Clasius-Clayperon方程    逸度    等量吸附热    页岩    延长组    
Difference analysis of isosteric heat of methane adsorption on shale based on fugacity and pressure: a case study of Yanchang Formation in Yanchang exploration area
XUE Pei1,2, ZHANG Lixia2, LIANG Quansheng2, SHI Yi3    
1. Northwest University of Political Science and Law, Xi'an 710063, China;
2. Research Institute, Shaanxi Yanchang Petroleum(Group) Co., Ltd., Xi'an 710075, China;
3. Yanchang Petroleum Group Exploration Company, Yan'an 716000, Shaanxi, China
Abstract: In order to improve the calculation method of isosteric heat of adsorption, clarify the thermodynamic characteristics of adsorption of CH4 by continental shale and reveal its adsorption mechanism, shale samples from Yanchang Formation in Yanchang exploration area were selected to carry out isothermal adsorption experiment of CH4 adsorption on shale at different temperatures, and the difference of isosteric heat of adsorption based on fugacity-absolute adsorption capacity and pressure-absolute adsorption capacity was analyzed by using absolute adsorption capacity curve. The results show that: (1) Fugacity is less than pressure, and in the pressure range of 0.36-2.21 MPa, the difference between fugacity and pressure is small. With the temperature decreasing and the pressure increasing, the difference increases.(2) The isotherm curves based on fugacity and pressure have obvious monotonic increasing linear characteristics, but the slope of isosteric heat of adsorption curve obtained by pressure-absolute adsorption capacity is larger than that obtained by fugacity-absolute adsorption capacity, which indicates that the intermolecular force of adsorbate has a great influence on the calculation results of isosteric heat of adsorption based on pressure-absolute adsorption capacity.(3) In the absolute adsorption capacity range of 0.01-0.35 mmol/g, the relative deviation of isosteric heat of adsorption obtained based on fugacity-absolute adsorption capacity and pressure-absolute adsorption capacity is -86.54%-57.01%. The data of fugacityabsolute adsorption capacity should be used as the basic data for the calculation of isosteric heart of adsorption in the thermodynamic evaluation of shale gas adsorption system.
Key words: Clasius-Clayperon equation    fugacity    isosteric heat of adsorption    shale    Yanchang Formation    
0 引言

页岩气的主要赋存空间为页岩中发育的微纳米级孔隙,其赋存状态具有吸附态、游离态和溶解态等3种形式,以吸附态和游离态为主,其中吸附气占比可达20%~85%[1]。由于吸附过程中,存在着气体分子动能及势能的变化,因此页岩气的吸附过程为放热。分析吸附过程的热力学特征有助于深入认识页岩气的吸附机理。

现有的页岩气及煤层气吸附热力学特征分析普遍采用基于Clasius-Clayperon方程的吸附等量线标绘法来计算等量吸附热[2],区别在于等量吸附热计算所用的等温吸附数据类型存在差异。文献[3-10] 采用基于压力的等温吸附数据作为等量吸附热计算的基础数据,而文献[11-15]采用基于逸度的等温吸附数据计算等量吸附热。由于经典的ClasiusClayperon方程是在假设吸附气为理想气体并忽略吸附相体积的基础上推导而来[16],然而页岩和煤等温吸附CH4过程的温度和压力范围均已超过CH4的临界点(CH4的临界温度为-82.6 ℃,临界压力为4.64 MPa),等温吸附过程属于超临界吸附的范畴[17],因此经典的Clasius-Clayperon方程并不适用于页岩与煤的吸附热力学特征分析,对于非理想气体的吸附热力学分析,需将Clasius-Clayperon方程中的压力替换为逸度[18]。目前尚未见到关于逸度与压力的选择对等量吸附热计算结果影响的量化评价。

笔者以延长探区甘泉地区延长组页岩为研究对象,开展不同温度下(30~70 ℃)页岩吸附CH4的等温吸附实验及绝对吸附量的校正,对比分析逸度与压力之间的差异,以及基于逸度-绝对吸附量和压力-绝对吸附量获得的等量吸附热的差异,以期完善页岩等量吸附热的计算方法,进而明确陆相页岩吸附热力学特征。

1 页岩等温吸附实验 1.1 实验样品与方法

等温吸附实验样品来自延长探区甘泉地区延长组长7段黑色页岩。实验样品地球化学分析表明,该样品总有机碳(TOC)质量分数为5.32%,有机质类型为Ⅱ1型,有机质成熟度(Ro)为0.97%。矿物组成分析表明,该样品黏土矿物含量较高,质量分数为52.4%,石英质量分数为24.5%,长石质量分数为14.1%,并含有少量的黄铁矿和菱铁矿[19]

等温吸附实验共设定30 ℃,40 ℃,50 ℃,60 ℃,70℃等5个温度点,实验压力为0.4~11.0 MPa,单个压力点下的吸附平衡时间不少于12 h,采用纯度99.99% 的甲烷为吸附气。样品预处理及实验操作规范参照了GB/T 19560—2008的相关规定[20]

1.2 过剩吸附量与绝对吸附量

延长组页岩样品等温吸附曲线如图 1所示。由于过剩吸附量表征了实际吸附相密度减去气相密度所剩部分的吸附量,过剩吸附量小于绝对吸附量。同时采用过剩吸附量计算获得的等量吸附热及初始等量吸附热小于采用绝对吸附量获得的等量及初始等量吸附热。因此需要将过剩吸附量校正为绝对吸附量,采用绝对吸附量计算等量吸附热。

下载原图 图 1 延长组页岩等温吸附曲线 Fig. 1 Adsorption isotherm of shale of Yanchang Formation

根据Gibbs吸附量的定义,过剩吸附量与绝对吸附量间的关系[21]

$ n_{\mathrm{ab}}=\frac{n_{\mathrm{ex}}}{1-\frac{\rho_{\mathrm{g}}}{\rho_{\mathrm{a}}}} $ (1)

式中:nab为绝对吸附量,mmol/g;nex为过剩吸附量,mmol/g;ρg为气相密度,g/cm3ρa为吸附相密度,g/cm3

采用真实气体状态方程计算气相密度,气相密度可表示为

$ \rho_{\mathrm{g}}=\frac{M P}{Z \mathrm{R} T} $ (2)

式中:M为摩尔质量,g/mol;R为通用气体常数,8.314 J/(K·mol);Z为真实气体压缩因子;T为绝对温度,K;P为气体压力,MPa。

本文采用Peng-Robinson方程[22]计算压缩因子。P-R(Peng-Robinson)方程表示为压缩因子Z的形式为

$ {Z^3} - (1 - B){Z^2} + \left( {A - 2B - 3{B^2}} \right)Z - \left( {AB - {B^2} - {B^3}} \right) = 0 $ (3)

参数AB的计算方法为

$ {A = \frac{{a(T)p}}{{{{\rm{R}}^2}{T^2}}}} $ (4)

$ {B = \frac{{bp}}{{{\rm{R}}T}}} $ (5)

式(4)、式(5)中a(T) 为内聚力参数,可表示为

$ a(T) = {a_{\rm{c}}} \cdot \alpha (T) = \left( {\frac{{0.452{\kern 1pt} {\kern 1pt} 724{{\rm{R}}^2}{T_{\rm{c}}}^2}}{{{p_{\rm{c}}}}}} \right) \cdot \alpha (T) $ (6)

其中

$ {\alpha (T) = {{\left[ {1 + k\left( {1 - {T_{\rm{r}}}^{0.5}} \right)} \right]}^2}} $ (7)

$ {b = \frac{{0.077{\kern 1pt} {\kern 1pt} 80{\rm{R}}{T_{\rm{c}}}}}{{{p_{\rm{c}}}}}} $ (8)

式中:Tc为临界温度,K;p为压力,Pa;pc为临界压力,Pa;α(T) 为与对比温度Tr和偏心因子ω有关的函数;b为协体积项;Tr为对比温度,Tr = T/Tck为与偏心因子ω有关的函数。

不同的ω取值范围内,k表达式不同,当偏心因子取值为0 < ω < 0.5时,k可表示为

$ k = 0.374{\kern 1pt} {\kern 1pt} 6 + 1.542{\kern 1pt} {\kern 1pt} 26\omega - 0.269{\kern 1pt} {\kern 1pt} 92{\omega ^2} $ (9)

式中:ω为偏心因子,CH4的偏心因子取值0.008。

经计算,得出在温度为30~70 ℃、压力为0~ 15 MPa时甲烷的气相密度(图 2)。

下载原图 图 2 甲烷的密度曲线 Fig. 2 Density curve of CH4

由于Ozawa等[23]经验公式法适用于较宽温压范围内的吸附相密度计算[19],因此本文采用Ozawa经验公式法计算吸附相密度,其计算公式为

$ {\rho _{\rm{a}}} = {\rho _{\rm{b}}} {\rm{exp}} [ - 0.002{\kern 1pt} {\kern 1pt} 5 \times (T - {T_{\rm{b}}})] $ (10)

式中:ρb表示沸点密度,g/cm3(CH4的沸点密度为0.424 g/cm3);Tb表示沸点,K(CH4的沸点为111.65 K,即-161.5℃)。

经计算,得出在温度为30~70 ℃时CH4的吸附相密度(图 3)。

下载原图 图 3 甲烷的吸附相密度曲线 Fig. 3 Adsorbed phase density curve of CH4

基于式(1),结合CH4的气相密度(图 2)和吸附相密度(图 3),即可实现由过剩吸附量(图 1)向绝对吸附量的校正,校正后的绝对吸附量曲线如图 4所示。采用Langmuir方程对绝对吸附量曲线进行了拟合,拟合相关系数在0.98以上,Langmuir方程拟合参数如表 1所列。表 1VL-pPL-p分别表示压力-绝对吸附量数据的Langmuir方程拟合参数。

下载原图 图 4 绝对吸附量曲线 Fig. 4 Absolute adsorption isotherm curves
下载CSV 表 1 压力-绝对吸附量曲线的Langmuir方程拟合参数 Table 1 Fitting parameters of Langmuir equations for the pressure-absolute adsorption isotherm
2 逸度 2.1 逸度的定义及计算方法

逸度及逸度系数是变组成体系热力学计算中的重要基本概念,也是相平衡和化学反应平衡研究的重要工具[24]。逸度是在处理非理想气体化学势时引入的概念。对于1 mol的纯物质的热力学关系式[25]

$ \mathrm{d} \mu=V_{\mathrm{m}} \mathrm{d} p $ (11)

式中:μ为化学势,J/mol;Vm为摩尔体积,m3/mol。

对于理想气体,其状态方程为pVm = RT,则

$ \mathrm{d} \mu_{\mathrm{ideal}}=\frac{\mathrm{R} T}{p} \mathrm{~d} p=\mathrm{R} T \mathrm{~d} \ln p $ (12)

对于真实气体,其状态方程为pVm = ZRT,则

$ \mathrm{d} \mu_{\mathrm{real}}=\frac{Z \mathrm{R} T}{p} \mathrm{~d} p \neq \mathrm{R} T \mathrm{~d} \ln p $ (13)

据文献[25]报道,1901年,Lewis提出将理想气体与真实气体的偏差放到压力项做一次性修正,并定义新函数,使得

$ \mathrm{d} \mu=\mathrm{R} T \mathrm{dln} f $ (14)

式中:f为逸度,Pa。对于理想气体,其逸度等于压力;对于真实气体,可把逸度看作校正的压力,也表征体系的逃逸趋势。

同时定义逸度系数为

$ \varphi = \frac{f}{p} $ (15)

逸度系数可用来衡量真实气体偏离理想状态的程度,φ越接近1,则气体越接近理想气体。纯组分的逸度和逸度系数可采用状态方程法、图解法、对比状态法和近似法计算[26]。本文采用P-R状态方程计算CH4的逸度,其逸度系数[18]可表示为

$ \ln \varphi=Z-1+\ln \frac{\mathrm{R} T}{Z \mathrm{R} T-p b}-\frac{\sqrt{2} a(T) \alpha(T)}{4 b \mathrm{R} T} \ln \frac{Z \mathrm{R} T+(\sqrt{2}+1) p b}{Z \mathrm{R} T-(\sqrt{2}-1) p b} $ (16)

2.2 逸度与逸度系数

利用式(16)及P-R方程即可得出逸度及逸度系数。当实验温度为30 ℃,40 ℃,50 ℃,60 ℃,70 ℃,压力为0.36~11.00 MPa时的逸度及逸度系数曲线如图 5所示。

下载原图 图 5 逸度及压力与逸度的差值曲线 Fig. 5 Fugacity and difference curve of pressure and fugacity

图 5(a)可看出,在低压力区间(0.36~ 2.21 MPa)内,逸度与压力基本相等,随着压力的增加,逸度与压力的差距逐渐增大。为了清晰地呈现二者的差异与温度和压力的关系,绘制了不同温度下压力与逸度差值随压力的变化曲线[图 5(b)]。从图 5(b)可看出,当压力大于逸度,特别是在低温高压条件下,压力与逸度的差值越发明显。具体而言,低压力阶段(0.36~2.21 MPa),压力与逸度的差值较小;处于0.002~0.010 MPa时,温度对二者差值的影响较小;随着压力的升高,二者的差值逐渐增加,以30℃为例,压力为0.36 MPa时,差值为0.002 MPa,当压力升高至11 MPa时,差值增加至2.10 MPa;同时随着压力的升高,温度对差值的影响越发明显,温度越低,二者的差值越大,例如当压力为11 MPa,温度为30 ℃,40 ℃,50 ℃,60 ℃,70 ℃时,二者的差值分别为2.01 MPa,1.78 MPa,1.59 MPa,1.43 MPa,1.28 MPa。

此外,逸度系数也反映了压力与逸度的差别。由于逸度系数定义为φ = f/p,且受限于不同温压条件下逸度与压力的关系,逸度系数小于1。随着压力的升高,逸度系数逐渐降低,同时,低温条件下的逸度系数随压力增加所表现出的下降速率要大于高温条件下的下降速率(图 6)。

下载原图 图 6 逸度系数曲线 Fig. 6 Fugacity coefficient curve

文中采用逸度-绝对吸附量数据作为等量吸附热计算的基础数据,因此需要将图 4中以压力作为横坐标的等温吸附曲线转换为以逸度作为横坐标的等温吸附曲线,以此作为热力学分析的基础数据(图 7)。由于低压力阶段,压力与逸度差值较小,而随着压力的增加,二者差值增大。因此,图 7所示的逸度-绝对吸附量曲线在高压部分数据点相对于压力-绝对吸附量曲线会向左平移,横坐标平移距离即为压力与逸度的差值,曲线形态相较于图 4所示的压力-绝对吸附量曲线发生改变。若将Langmuir方程中的参数p用逸度f替换,方程形式保持不变,对逸度-绝对吸附量数据进行了拟合。拟合结果表明,方程拟合相关系数均在0.98以上。由于在相同的温度和压力条件下,逸度小于压力,且随着压力的升高,二者差值增大,因此相对于图 4而言,逸度-绝对吸附量的Langmuir方程的拟合参数Langmuir体积(VL-f)及Langmuir压力(PL-f)均略大于压力-绝对吸附量的拟合结果(表 1表 2)。

下载原图 图 7 基于逸度的绝对吸附量曲线 Fig. 7 Absolute adsorption curve based on fugacity
下载CSV 表 2 逸度-绝对吸附量曲线的Langmuir方程拟合参数 Table 2 Fitting parameters of Langmuir equations for fugacity-absolute adsorption isotherm
3 等量吸附热 3.1 等量吸附热的计算方法

等量吸附热(qst)是指当压力(p)、温度(T)及吸附剂表面积(A)恒定时,吸附剂表面吸附无限小量的吸附质,吸附体系所释放的热量,其数学表达式[27]

$ q_{\mathrm{st}}=T\left(v_{\mathrm{g}}-v_{\mathrm{a}}\right) \frac{\mathrm{d} p}{\mathrm{~d} T} $ (17)

式中:qst为等量吸附热,kJ/mol;vg为体相气体的摩尔体积,m3/mol;va为吸附相气体的摩尔体积,m3/mol。

若忽略吸附相体积,且视体相气体为理想气体,则式(17)可写作经典的Clasius-Clayperon方程[10],即

$ q_{\mathrm{st}}=\mathrm{R} T^{2} \frac{\mathrm{d} \ln p}{\mathrm{~d} T} $ (18)

由于页岩及煤岩等温吸附CH4实验的温度和压力范围大多分别集中在20~80 ℃及0~20 MPa[28-29],同时CH4的临界温度为-82.6 ℃,临界压力为4.64 MPa,在此温度和压力范围内的等温吸附过程存在着由气相吸附向超临界流体吸附的转变,气体分子行为已不再满足理想气体假设。以CH4压缩因子为例(图 8),不同温度下CH4的压缩因子随压力的增加,先减小后增加,偏离了理想气体状态。因此针对非理想气体的热力学分析,若采用ClasiusClayperon方程计算等量吸附热,则需要将压力p替换为逸度f[18, 30-31],即将式(18)转变为

$ {q_{{\rm{st}}}} = {\rm{R}}{T^2}\frac{{{\rm{d}}\ln f}}{{{\rm{d}}T}} $ (19)

下载原图 图 8 CH4的压缩因子(采用P-R方程计算) Fig. 8 Compression factor of CH4 by P-R equation

进而采用吸附等量线标绘法计算固-气吸附体系的等量吸附热。

基于Clasius-Clayperon方程的吸附等量线标绘法在等量吸附热计算中应用较为广泛[2]。该方法首先需要对式(19)两端变形并积分可得

$ \ln f = - \frac{{{q_{{\rm{st}}}}}}{{{\rm{R}}T}} + {C_2} $ (20)

由式(20)可知,逸度的自然对数lnf与绝对温度的倒数1/T呈线性关系,由斜率得到吸附量n对应的qst。假设斜率为A,则等量吸附热计算式为

$ {q_{{\rm{st}}}} = - {\rm{R}}A $ (21)

3.2 基于逸度与压力的等量吸附热特征对比

笔者分别采用逸度-绝对吸附量和压力-绝对吸附量数据对页岩吸附CH4的等量吸附热进行计算。同时,在后续行文过程中将基于逸度-绝对吸附量数据获得的等量吸附热标记为qst-f,相对应地将基于压力-绝对吸附量数据获得的等量吸附热标记为qst-p,以便区分。

首先,依据lnp-nab及lnf-nab的数据特征分别选用线性函数和幂函数对lnp-nab及lnf-nab进行了拟合(图 9)。拟合结果表明,lnpnab呈明显的线性关系,lnfnab呈明显的幂函数关系,函数拟合相关系数均超过0.99。

下载原图 图 9 lnp-nab及lnf-nab数据的函数拟合 Fig. 9 Fitting curve of lnp-nab and lnf-nab

利用拟合得到的lnp-nab与lnf-nab的线性函数、幂函数分别计算不同绝对吸附量下的lnp及lnf,进而绘制不同绝对吸附量下的lnp-1/T及lnf-1/T曲线,并分别对曲线进行线性函数拟合(图 10),lnp-1/T及lnf-1/T曲线均符合线性相关性,且相关系数均超过0.95。

下载原图 图 10 lnp-1/T(a)及lnf-1/T(b)数据的线性函数拟合 Fig. 10 Linear fitting curves of lnp-1/T(a)and lnf-1/T(b)

分别利用lnp-1/T[图 10(a)]及lnf-1/T[图 10(b)]曲线线性拟合函数的斜率结合式[式(21)]计算基于压力-绝对吸附量及逸度-绝对吸附量的等量吸附热,计算结果如图 11所示。

下载原图 图 11 基于逸度-绝对吸附量及压力-绝对吸附量获得的等量吸附热 Fig. 11 Isosteric heat of adsorption based on fugacity-absolute adsorption isotherm and pressure-absolute adsorption isotherm

基于逸度的等量吸附热及基于压力的等量吸附热均与绝对吸附量满足良好的线性关系,相关系数均在0.99以上(表 3),且随着绝对吸附量的增加,等量吸附热均逐渐升高(图 11)。绝对吸附量由0.01 mmol/g上升至0.35 mmol/g时,基于逸度的等量吸附热为22.31~72.35 kJ/mol,基于压力的等量吸附热为9.59~134.97 kJ/mol。

下载CSV 表 3 等量吸附热曲线的线性函数拟合参数 Table 3 Linear function fitting parameters of isosteric heat of adsorption

等量吸附热的大小与固体表面的强/弱吸附位点分布的均一性及吸附气分子间作用力有关[32]。吸附过程中,气体分子会优先在固体表面的强吸附位点吸附,随着吸附的进行,当强吸附位点吸附逐渐饱和后,气体分子逐步在弱吸附位点吸附,这造成了等量吸附热随着吸附量的增加而降低;同时,随着吸附量的增加,固体表面的吸附态分子密度增加,吸附气分子间的作用力增加,这造成了等量吸附热随吸附量增加而升高。由此可见,CH4在甘泉地区长7段泥岩孔隙表面的吸附过程中,吸附质分子间的作用力是存在的,且相对于吸附剂表面的非均一性,前者吸附质分子间的作用力对等量吸附热的影响是主要的。

在采用等量吸附线标绘法计算等量吸附时,压力或逸度的选择并不会改变等量吸附热曲线的线性单调性特征,二者均反映了等量吸附热具有明显的单调递增的线性特征。区别在于,基于逸度的等量吸附热(qst-f)曲线的斜率小于基于压力的等量吸附热(qst-p)曲线。由于等量吸附热随吸附量增加而上升是受吸附质分子间相互作用的影响,随着吸附的进行,吸附剂表面吸附位点对吸附质分子的作用力逐渐减弱,等量吸附热曲线将更多地反映出吸附质分子间作用力的影响。因此可以认为,若采用压力-绝对吸附量计算等量吸附热,则会造成等量吸附热计算结果中吸附质分子间作用力的影响过大。

此外,在绝对吸附量小于0.065 mmol/g时,基于逸度的等量吸附热(qst-f)大于基于压力的等量吸附热(qst-p);绝对吸附量大于0.065 mmol/g时,前者(qst-f)小于后者(qst-p);二者在绝对吸附量为0.065 mmol/g处曲线相交。为定量评价逸度与压力在等量吸附热计算结果中的误差,引入相对误差(δ)的概念,即

$ \delta = \frac{{{q_{{\rm{st - f}}}} - {q_{{\rm{st - p}}}}}}{{{q_{{\rm{st - f}}}}}} \times 100\% $ (22)

经计算,绝对吸附量为0.01~0.35 mmol/g时,相对误差为-86.54 %~57.01 %(图 12),若以±5% 作为合理的相对偏差范围,那么仅在绝对吸附量为0.060~0.074 mmol/g时,相对偏差处于合理范围。特别在低吸附量,即低压力阶段,基于逸度与压力的等量吸附热仍然具有较高的误差。例如,当低吸附量为nab = 0.01 mmol/g时,相对偏差高达52%。其原因在于,虽然在低压力条件下,压力与逸度的差值较小,但在高压力阶段,逸度与压力的差值较大,致使高压阶段中逸度-绝对吸附量曲线的形态与压力-绝地吸附量曲线存在较大的差异,进而影响到吸附等量线标绘法计算过程中lnp-nab与lnf-nab曲线的拟合函数发生改变,因此在低压阶段,基于压力与逸度的等量吸附热曲线出现了较高的相对偏差。由此可见,采用压力-绝对吸附量作为计算等量吸附热的基础数据,会造成较大的误差。

下载原图 图 12 等量吸附热的相对偏差 Fig. 12 Relative deviation of isosteric heat of adsorption
4 结论

(1)在低压力区(0.36~2.21 MPa)内,逸度与压力基本相等,随着压力的升高,逸度逐渐小于压力,且温度越低,压力越高,压力与逸度的差值越大。

(2)采用吸附等量线标绘法计算等量吸附热时,压力与逸度的选择并不会改变等量吸附热曲线的线性单调性特征,二者均反映了等量吸附热具有明显的单调递增的线性特征,但若采用压力-绝对吸附量计算等量吸附热,则会造成等量吸附热计算结果中吸附质分子间作用力的影响过大。

(3)基于逸度的等量吸附热与基于压力的等量吸附热的相对偏差为-86.54%~57.01%,仅在绝对吸附量为0.060~0.074 mmol/g时,二者的相对偏差处于±5% 的合理范围内。页岩气吸附体系的热力学评价应当采用逸度-绝对吸附量数据作为等量吸附热计算的基础数据。

参考文献
[1]
CURTIS J B. Fractured shale-gas systems. AAPG Bulletin, 2002, 86(11): 1921-1938.
[2]
DONGMIN S, MARTIN B, FLOR S, et al. Comparison of experimental techniques for measuring isosteric heat of adsorption. Adsorption, 2000, 6(4): 275-286. DOI:10.1023/A:1026551213604
[3]
周来, 冯启言, 秦勇. CO2和CH4在煤基质表面竞争吸附的热力学分析. 煤炭学报, 2011, 36(8): 1307-1311.
ZHOU L, FENG Q Y, QIN Y. Thermodynamic analysis of competitive adsorption of CO2 and CH4 on coal matrix. Journal of China Coal society, 2011, 36(8): 1307-1311.
[4]
崔永君, 张庆玲, 杨锡禄. 不同煤的吸附性能及等量吸附热的变化规律. 天然气工业, 2003, 23(7): 130-131.
CUI Y J, ZHANG Q L, YANG X L. The change rules of adsorption capacity and isosteric heat of adsorption of different coals. Natural Gas Industry, 2003, 23(7): 130-131.
[5]
卢守青, 王亮, 秦立明. 不同变质程度煤的吸附能力与吸附热力学特征分析. 煤炭科学技术, 2014(6): 130-135.
LU S Q, WANG L, QIN L M. Analysis on adsorption capacity and adsorption thermodynamic characteristics of different metamorphic degree coals. Coal Science and Technology, 2014(6): 130-135.
[6]
郭为, 熊伟, 高树生, 等. 温度对页岩等温吸附/解吸特征影响. 石油勘探与开发, 2013, 40(4): 481-485.
GUO W, XIONG W, GAO S H, et al. Impact of temperature on the isothermal adsorption/desorption characteristics of shale gas. Petroleum Exploration and Development, 2013, 40(4): 481-485.
[7]
李晓媛, 曹峰, 岳高凡, 等. 柴达木盆地东部石炭系页岩吸附特性实验研究. 地学前缘, 2016, 23(5): 95-102.
LI X Y, CAO F, YUE G F, et al. The experimental study of adsorption characteristics of Carboniferous shale in eastern Qaidam. Earth Science Frontiers, 2016, 23(5): 95-102.
[8]
白建平, 张典坤, 杨建强, 等. 寺河3号煤甲烷吸附解吸热力学特征. 煤炭学报, 2014, 39(9): 1812-1819.
BAI J P, ZHANG D K, YANG J Q, et al. Thermodynamic characteristics of adsorption-desorption of methane in coal seam 3 at Sihe coal mine. Journal of China Coal Society, 2014, 39(9): 1812-1819.
[9]
李树刚, 白杨, 林海飞, 等. CH4, CO2和N2多组分气体在煤分子中吸附热力学特性的分子模拟. 煤炭学报, 2018, 43(9): 114-121.
LI S G, BAI Y, LIN H F, et al. Molecular simulation of adsorption thermodynamics of multicomponent gas in coal. Journal of China Coal Society, 2018, 43(9): 114-121.
[10]
林海飞, 蔚文斌, 李树刚, 等. 煤体吸附CH4及CO2热力学特性试验研究. 中国安全科学学报, 2018, 28(6): 129-134.
LIN H F, WEI W B, LI S G, et al. Experimental study on thermodynamics characteristics of CH4 and CO2 adsorption on coal. China Safety Science Journal, 2018, 28(6): 129-134.
[11]
蔺亚兵, 马东民, 刘钰辉, 等. 温度对煤吸附甲烷的影响实验. 煤田地质与勘探, 2012(6): 24-28.
LIN Y B, MA D M, LIU Y H, et al. Experiment of the influence of temperature on coalbed methane adsorption. Coal Geology & Exploration, 2012(6): 24-28.
[12]
岳高伟, 王兆丰, 康博. 基于吸附热理论的煤-甲烷高低温等温吸附线预测. 天然气地球科学, 2015, 26(1): 148-153.
YUE G W, WANG Z F, KANG B. Prediction for isothermal adsorption curve of coal/CH4 based on adsorption heat theory. Natural Gas Geoscience, 2015, 26(1): 148-153.
[13]
杨峰, 宁正福, 王庆, 等. 甲烷在页岩上的吸附热力学. 中南大学学报(自然科学版), 2014, 45(8): 2871-2877.
YANG F, NING Z F, WANG Q, et al. Thermodynamic analysis of methane adsorption on gas shale. Journal of Central South University(Science and Technology), 2014, 45(8): 2871-2877.
[14]
杨峰, 宁正福, 刘慧卿, 等. 页岩对甲烷的等温吸附特性研究. 特种油气藏, 2013, 20(5): 133-136.
YANG F, NING Z F, LIU H Q, et al. Methane adsorption characteristics of gas shale. Special Oil & Gas Reservoirs, 2013, 20(5): 133-136.
[15]
马东民, 曹石榴, 李萍, 等. 页岩气与煤层气吸附/解吸热力学特征对比. 煤炭科学技术, 2015, 43(2): 64-67.
MA D M, CAO S L, LI P, et al. Comparison on adsorption and desorption thermodynamics features between shale gas and coalbed methane. Coal Science and Technology, 2015, 43(2): 64-67.
[16]
PAN H, RITTER J A, BALBUENA P B. Examination of the approximations used in determining the isosteric heat of adsorption from the Clausius-Clapeyron equation. Langmuir, 1998, 14(21): 6323-6327.
[17]
周理, 吕昌忠, 王怡林, 等. 述评超临界温度气体在多孔固体上的物理吸附. 化学进展, 1999, 11(3): 221-226.
ZHOU L, LYU C Z, WANG Y L, et al. Physisorption of gases on porous solids at above-critical temperatures. Progress in Chemistry, 1999, 11(3): 221-226.
[18]
ZHENG Y N, LI Q Z, YUAN C C, et al. Thermodynamic analysis of high-pressure methane adsorption on coal-based activated carbon. Fuel, 2018, 230: 172-184.
[19]
薛培, 张丽霞, 梁全胜, 等. 陆相页岩吸附CH4的热力学特征. 天然气工业, 2019, 39(11): 64-73.
XUE P, ZHANG L X, LIANG Q S, et al. Thermodynamic characteristics of CH4 adsorption by continental shale. Natural Gas Industry, 2019, 39(11): 64-73.
[20]
中华人民共和国国家质量监督检验检疫总局, 中华人民共和国国家标准化管理委员会. GB/T 19560-2008煤的高压等温吸附试验方法. 北京: 中国标准出版社, 2008.
State Administration of Quality Supervision, Inspection and Quarantine of the PRC & Standardization Administration of the PRC. GB/T 19560-2008 Experimental method of high-pressure isothermal adsorption to coal. Beijing: Standard Press of China, 2008.
[21]
周尚文, 王红岩, 薛华庆, 等. 页岩过剩吸附量与绝对吸附量的差异及页岩气储量计算新方法. 天然气工业, 2016, 36(11): 12-20.
ZHOU S W, WANG H Y, XUE H Q, et al. Difference between excess and absolute adsorption capacity of shale and a new shale gas reserve calculation method. Natural Gas Industry, 2016, 36(11): 12-20.
[22]
ROBINSON D B, PENG D Y, CHUNG Y K. The development of the Peng-Robinson equation and its application to phase equilibrium in a system containing methanol. Fluid Phase Equilibria, 1985, 24(1/2): 25-41.
[23]
OZAWA S, KUSUMI S, OGINO Y. Physical adsorption of gases at high pressure. Ⅳ. An improvement of the Dubinin-Astakhov adsorption equation. Journal of Colloid and Interface Science, 1976, 56(1): 83-91.
[24]
张珏成. 纯物质逸度定义的讨论. 上海工程技术大学学报, 2004, 18(2): 14-17.
ZHANG J C. Discussion about fugacity definition of pure matter. Journal of Shanghai University of Engineering Science, 2004, 18(2): 14-17.
[25]
鄢浩, 陈晋阳. 化工热力学. 北京: 中国石化出版社, 2011.
YAN H, CHEN J Y. Chemical engineering thermodynamics. Beijing: China Petrochemical Press, 2011.
[26]
傅献彩, 沈文霞, 姚天扬, 等. 物理化学. 北京: 高等教育出版社, 2005.
FU X C, SHEN W X, YAO T Y, et al. Physical chemistry. Beijing: Higher Education Press, 2005.
[27]
AZAHAR F H M, MITRA S, YABUSHITA A, et al. Improved model for the isosteric heat of adsorption and impacts on the performance of heat pump cycles. Applied Thermal Engineering, 2018, 143: 688-700.
[28]
MENG M, QIU Z S, ZHONG R Z, et al. Adsorption characteristics of supercritical CO2/CH4 on different types of coal and a machine learning approach. Chemical Engineering Journal, 2019, 368: 847-864.
[29]
KONG S Q, HUANG X, LI K J, et al. Adsorption/desorption isotherms of CH4 and C2 H6 on typical shale samples. Fuel, 2019, 255: 115623.
[30]
ZHENG Q R, JI X W, GAO S, et al. Analysis of adsorption equilibrium of hydrogen on graphene sheets. International Journal of Hydrogen Energy, 2013, 38(25): 10896-10902.
[31]
WANG K, QIAO S Z, HU X J. Study of isosteric heat of adsorption and activation energy for surface diffusion of gases on activated carbon using equilibrium and kinetics information. Separation and Purification Technology, 2004, 34(1/3): 165-176.
[32]
RUTHVEN D M. Principle of adsorption and adsorption process. New York: John Wilev & Sons Inc., 1984.