-
在LANGMUIR（兰格苗尔）等温累积吸附量经验方程[1] 发表100周年之际，陈元千等于2018年发表了等温吸附量方程和等温解吸方程的推导[2]。应当指出，无论是等温吸附量曲线的数据点，或是等温累积吸附量曲线的数据点，都不是由吸附仪直接测试得到的，而是利用吸附仪测试的数据，由相关的计算公式求得的。根据吸附仪测试的机理和特点，可将目前的两类吸附仪分别称为测压吸附仪和称重吸附仪。对于测压吸附仪，MAVOR等[3-4] 未经推导直接给出的计算方法是不正确的。文献[5]已经提出了计算等温吸附量的方法。到目前为止，在国外文献中尚看不到有关称重吸附仪计算等温吸附量方法的报道。自2012年以来，中国有关单位以几十万到上百万美元不同的价格，从美国、英国、德国、荷兰和澳大利亚等国家购进了许多台称重吸附仪。然而，令人费解的是，仪器销售公司为何不直接有效地提供等温吸附量的计算方法及其软件系统，致使中国页岩气和煤层气重点实验室的专家和教授，在苦忙于进行大量实验工作的同时，还要开展计算方法的研究工作。在对中国专家、教授发表的文章[6-14] 分析后认为，他们所提出的所谓过剩吸附量和绝对吸附量，以及两者之间的关系式都是不正确的。应当指出，人们只知道吸附量和不再吸附时的饱和吸附量，哪里来的过剩吸附量和绝对吸附量。因此，过剩吸附量是一个子虚乌有的伪命词； 绝对吸附量是一个牵强附会的借用词；两者之间的关系式也是一个缺少严格理论推导以人为定义直接给出的计算公式。
-
1 计算等温吸附量方法的推导
-
当利用称重吸附仪进行注入甲烷气吸附实验时，在某一稳定压力点测试的总质量为 m t，等于测试桶的质量 m b、测试桶内装入岩样的质量 m c、岩样吸附甲烷气的质量 m ag和测试桶空隙体积内自由甲烷气的质量m fg的总和，其质量平衡关系式为：
-
由于测试桶的质量和测试桶内岩样的质量均为常数，因此，到某一压力测试点的累积注入甲烷气的质量m i为：
-
在升压注入甲烷气测试过程中，根据（2）式可得两个相邻压力测试点的质量增量为：
-
由（3）式可得岩样吸附甲烷气的质量增量为：
-
已知气体的摩尔量等于气体的质量除以气体的分子量，因此，由（4）式可得在两个相邻升压测试点之间吸附甲烷气的摩尔量增量为：
-
假设1克岩样吸附甲烷气的摩尔量为q m，那么， m c克岩样吸附甲烷气的摩尔量增量为：
-
根据Avogadro定律[15]，1摩尔的任何气体，在273.15K和1atm（物理大气压）下，占有的气体体积为22 414cm3。那么，在293K和0.101MPa（1atm）的地面标准条件下，1摩尔的任何气体，占有的气体体积应为24 056cm3。因此，Δn ag摩尔的甲烷气，在293K和0.101MPa（1atm）条件下，占有的气体体积应为24 056Δn ag cm3。因此，由（6）式可得1克岩样对甲烷气的吸附量为：
-
将（5）式代入到（7）式得：
-
已知甲烷气的分子量为16.04g/mol，将其代入 （8）式得：
-
应当指出，（7）式中q的单位为cm3/g，由于1m3=106cm3 和1t=106g，因此由（7）式计算 q 的单位等同于m3/t。
-
在测试桶内的空隙体积中自由甲烷气的质量增量为：
-
将（10）式代入（9）式得：
-
则单位岩样的等温累积吸附量，由累加法表示为：
-
2 判断吸附甲烷气程度的方法
-
由（9）式看出，当 时，q=0。因此，可以写出甲烷气吸附程度的判断因子为：
-
由（9）式和（13）式看出，当 η> 0时，q> 0；当 η=0时，q=0；当 η< 0时，q< 0（q为负值），负值的情况是不存在的。
-
3 等温累积吸附方程及求解方法
-
在文献[2]中，陈元千等建立的等温累积吸附量方程为：
-
为了求得方程常数A和B值，将（14）式改写为：
-
利用称重吸附仪的计算方法，在求得不同压力 p下的累积吸附量Q值之后，可利用不同压力下的Q 值，给定不同的A值，利用（15）式的线性迭代试差求解，可以得到相关系数最高的直线关系，经线性回归后，求得正确的B值。
-
4 气体密度和气体偏差系数的计算方法
-
由文献[16]可以写出不同压力和温度下甲烷气的密度为：
-
其中
-
甲烷气的偏差系数可利用相关经验公式[16] 计算，其表达式为：
-
其中
-
5 应用实例
-
已知称重吸附仪实验的有关基础数据为：T=323.5K，m b=124.947g，页岩气的 m c=118.269g，v tv=68.787cm3；煤层气的m c=79.6g，v tv=62.557cm3，页岩气和煤层气吸附的实验数据和有关计算结果分别列于表1和表2。
-
利用本文方法计算得到页岩气和煤层气的等温累积吸附量数据，分别列于表1和表2，绘制其与压力的关系曲线（图1，图2）。根据表1和表2中的 Q值，利用（15）式的线性迭代试差法，求得页岩气等温吸附方程常数 A 和 B 值分别为2.66和0.304 9，相关系数为0.994 3；求得煤层气等温吸附方程常数 A 和 B 值，分别为20.20和0.626 5，相关系数为0.996 5。将页岩气和煤层气的 A 和 B 值，分别代入 （14）式，可以得到预测页岩气和煤层气的等温累积吸附量与吸附压力的关系式为：
-
图1 页岩气的等温累积吸附量与压力的关系
-
Fig.1 Plot of isothermal cumulative adsorption capacity vs p of shale gas
-
图2 煤层气的等温累积吸附量与压力的关系
-
Fig.2 Plot of isothermal cumulative adsorption capacity vs p of coalbed methane
-
利用（25）式和（26）式，预测得到页岩气和煤层气的累积吸附量曲线（图1，图2）。由图1和图2可以看出，页岩气和煤层气累积吸附量的预测曲线与与实际测试资料计算的数据点吻合得很好。利用交会法（图3，图4）确定页岩气和煤层气的饱和吸附压力（p s）分别为10.6和9.2MPa。将其分别代入 （25）式和（26）式，求得页岩气和煤层气的饱和累积吸附量（Q s）分别为2.56和20.14m3/t。
-
图3 页岩气的Δm i与Δm fg交会图
-
Fig.3 Plot of Δm i and Δm fg of shale gas
-
图4 煤层气的Δm i与Δm fg交会图
-
Fig.4 Plot of Δm i and Δm fg of coalbed methane
-
6 结论
-
应用质量平衡原理和Avogadro定律，经过理论推导得到了称重吸附仪计算等温吸附量的方法和判断吸附甲烷气程度的方法，以及确定饱和吸附压力和饱和累积吸附量的方法。实例应用结果表明，所建立的计算方法是正确有效的。但应当指出，由文献[7-14]没有任何物理定义所提出和应用的所谓过剩吸附量和绝对吸附量，及其两者之间的关系式，都是不正确的。他们提出的所谓过剩吸附量和绝对吸附量，是两个子虚乌有和牵强附会的伪命词。而两者之间的关系式，也是一个未经严格推导和没有引用文献，人为直接给出的结果。
-
致谢：对为本文提供技术咨询和实验测试数据的中国石油勘探开发研究院邓泽高级工程师表示诚挚的感谢！
-
符号解释：
-
m t ——某一测试压力点称重吸附仪测试的总质量，g； m b——测试桶的质量，g；m c——实验岩样的质量，g；m ag—— 岩样吸附甲烷气的质量，g；m fg——测试桶的空隙体积内自由甲烷气的质量，g；m i ——某一测试压力点累积注入甲烷气质量，g；Δm i ——某一测试压力点注入甲烷气的质量增量，g；Δm ag ——某一测试压力点吸附甲烷气的质量增量，g； Δm fg——某一测试压力点测试桶空隙体积内自由甲烷气的质量增量，g；Δn ag——某一测试压力点吸附甲烷气的摩尔量增量，mol；M ——甲烷气的分子量，g/mol；q m——1克岩样质量吸附甲烷气的摩尔量，mol/g；q ——1克岩样质量对甲烷气的吸附量，cm3/g；ρg——自由甲烷气的密度，g/cm3；v tv——由氦气标定的测试桶内的空隙体积，cm3；Q ——等温累积吸附量，cm3/g；η——判断因子，dim；A，B ——陈氏等温累积吸附方程常数；p ——称重吸附仪测试的压力，MPa；γg——甲烷气的相对密度，dim；Z ——甲烷气的偏差系数，dim；T ——称重吸附仪测试的恒定温度，K；M air——空气的分子量，g/mol； p r ——甲烷气的对比压力，dim；T r ——甲烷气的对比温度， dim；p c——甲烷气的临界压力，MPa；T c——甲烷气的临界温度，K；p s ——饱和吸附压力，MPa；Q s ——饱和累积吸附量， cm3/g。
-
参考文献
-
[1] LANGMUIR I.The adsorption of gases on plane surfaces of glass，mica and platinum［J］.Journal of the American Chemical Society，1918，40（9）：1 361-1 402.
-
[2] 陈元千，傅礼兵，郝明强.气体吸附方程和解吸方程的推导及应用［J］.中国海上油气，2018，30（2）：85-89.CHEN Yuanqian，FU Libing，HAO Mingqiang.Derivation and ap⁃ plication of gas adsorption equation and desorption equation［J］.China Offshore Oil and Gas，2018，30（2）：85-89.
-
[3] MAVOR M J，OWEN L B，PRATT T J.Measurement and evalua⁃ tion of coal sorption isotherm data［C］.SPE 20728，1990.
-
[4] MAVOR M J，HARTMAN C P T，TIMOTHY J.Uncertainty in sorption isotherm measurements［C］.University of Alabama，Tus⁃ caloosa：Proceedings 2004 International Coalbed Methane Sympo⁃ sium，2004.
-
[5] 陈元千，汤晨阳，陈奇.等温吸附量计算方法的推导及应用［J］.油气地质与采收率，2018，25（6）：56-62.CHEN Yuanqian，TANG Chenyang，CHEN Qi.Derivation and ap⁃ plication of isothermal adsorption rate calculation method［J］.Pe⁃ troleum Geology and Recovery Efficiency，2018，25（6）：56-62.
-
[6] 周理，李明，周亚平.超临界甲烷在高表面活性炭上的吸附测量及其理论分析［J］.中国科学：B辑，2000，30（1）：49-56.ZHOU Li，LI Ming，ZHOU Yaping.Adsorption measurement and theoretical analysis of supercritical methane on high surface area activated carbon［J］.Science China：Series B，2000，30（1）：49-56.
-
[7] 俞凌杰，范明，陈红宇，等.富有机质页岩高温高压重量法等温吸附实验［J］.石油学报，2015，36（5）：557-563.YU Lingjie，FAN Ming，CHEN Hongyu，et al.Isothermal adsorp⁃ tion experiment of organic-rich shale under high temperature and pressure using gravimetric method［J］.Acta Petrolei Sinica，2015，36（5）：557-563.
-
[8] 周尚文，王红岩，薛华庆，等.页岩过剩吸附量与绝对吸附量的差异及页岩气储量计算新方法［J］.天然气工业，2016，36（1）：12-20.ZHOU Shangwen，WANG Hongyan，XUE Huaqing，et al.Differ⁃ ence between excess and absolute adsorption capacity of shale and a new shale gas reserve calculation method［J］.Natural Gas Industry，2016，36（1）：12-20.
-
[9] 周尚文，薛华庆，郭伟，等.基于重量法的页岩气超临界吸附特征实验研究［J］.煤炭学报，2016，41（11）：2 806-2 812.ZHOU Shangwen，XUE Huaqing，GUO Wei，et al.Supercritical isothermal adsorption characteristics of shale gas based on gravi⁃ metric method［J］.Journal of China Coal Society，2016，41（11）：2 806-2 812.
-
[10] 周尚文，王红岩，薛华庆，等.页岩气超临界吸附机理及模型 ［J］.科学通报，2017，62（35）：4 189-4 200.ZHOU Shangwen，WANG Hongyan，XUE Huaqing，et al.Super⁃ critical methane adsorption on shale gas：Mechanism and model（in Chinese）［J］.Chinese Science Bulletin，2017，62（35）：4 189-4 200.
-
[11] 周尚文，李奇，薛华庆，等.页岩容量法和重量法等温吸附实验对比研究［J］.化工进展，2017，36（5）：1 690-1 697.ZHOU Shangwen，LI Qi，XUE Huaqing，et al.Comparative study on the volumetric and gravimetric method for isothermal adsorp⁃ tion experiment of shale［J］.Chemical Industry and Engineering Progress，2017，36（5）：1 690-1 697.
-
[12] 朱汉卿，贾爱林，位云生，等.蜀南地区富有机质页岩孔隙结构及超临界甲烷吸附能力［J］.石油学报，2018，39（4）：391-401.ZHU Hanqing，JIA Ailin，WEI Yunsheng，et al.Pore structure and supercritical methane sorption capacity of organir-rich shales in southern Sichuan Basin［J］.Acta Petrolei Sinica，2018，39（4）：391-401.
-
[13] 田文广，邓泽，王红岩，等.低吸附量煤岩和页岩等温吸附实验的负吸附原因［J］.天然气工业，2018，38（6）：19-26.TIAN Wenguang，DENG Ze，WANG Hongyan，et al.Negative ad⁃ sorption in the isotherm adsorption experiments of low-adsorption coal and shale［J］.Natural Gas Industry，2018，38（6）：19-26.
-
[14] 端祥刚，胡志明，高树生，等.页岩高压等温吸附曲线及气井生产动态特征实验［J］.石油勘探与开发，2018，45（1）：119-127.DUAN Xianggang，HU Zhiming，GAO Shusheng，et al.Shale high pressure isothermal adsorption curve and production dynamic ex⁃ periments of gas well［J］.Petroleum Exploration and Develop⁃ ment，2018，45（1）：119-127.
-
[15] 萨韦利耶夫 И В.普通物理学（第一卷，力学与分子物理学）［M］.钟金城，何伯珩，译.北京：高等教育出版社,1992.SAVELYEV И В.General physics(Volume I，Mechanics and mo⁃ lecular physics)［M］.ZHONG Jincheng，HE Boheng，trans.Bei⁃ jing：Higher Education Press，1992.
-
[16] 陈元千，李璗.现代油气藏工程［M］.北京：石油工业出版社，2001：1-11.CHEN Yuanqian，LI Dang.Modern petroleum engineering［M］.Beijing：Petroleum Industry Press，2001：1-11.
-
摘要
当利用称重吸附仪进行吸附实验时，某一测试压力点测试的总质量m t ，是由测试桶的质量m b、测试桶内装入岩样的质量m c、岩样吸附甲烷气的质量m ag和测试桶空隙体积内自由甲烷气的质量m fg等四部分组成。而注入甲烷气的质量m i仅包括m ag和m fg两部分。由于测试桶的质量和装入测试桶内的岩样质量都是常数，因此为建立利用称重吸附仪计算等温吸附量的方法提供了方便。根据质量平衡原理和 Avogadro定律，经理论推导，得到利用称重吸附仪测试数据计算等温吸附量的方法。同时，还提出了判断甲烷气吸附程度的方法，以及确定饱和吸附压力和饱和累积吸附量的方法。页岩气和煤层气2个实例的应用结果表明，本文提供的方法是正确有效的。
Abstract
When the adsorption experiment is carried out by using weighing method，the total mass m t at a test pressure point is composed by four parts，which are the mass of the test barrel m b，the mass of the core sample m c，the adsorbed meth- ane gas m ag and the free methane gas m fg in the void volume of the test barrel. And the mass of the injected methane gas m i only includes two parts，m ag and m fg. Since the mass of the sample barrel and the mass of the sample loaded into the sample barrel are constant，this provides convenience for establishing a method for calculating the isothermal adsorption capacity by the weighing adsorption apparatus. Based on the principle of mass balance and Avogadro’s law，the method of calculat- ing isothermal adsorption capacity using the data tested by the weighing adsorption apparatus was theoretically derived. Meanwhile，a criteria for judging adsorption degree，determining the saturated adsorption pressure and the saturated cumu- lative adsorption capacity was proposed. The application of shale gas and coalbed methane adsorption showed that the said method was correct and effective.
