CN106198297B - 一种精确计算甲烷在页岩上真实吸附量的方法 - Google Patents
一种精确计算甲烷在页岩上真实吸附量的方法 Download PDFInfo
- Publication number
- CN106198297B CN106198297B CN201610482070.4A CN201610482070A CN106198297B CN 106198297 B CN106198297 B CN 106198297B CN 201610482070 A CN201610482070 A CN 201610482070A CN 106198297 B CN106198297 B CN 106198297B
- Authority
- CN
- China
- Prior art keywords
- sample
- total
- cylinder
- absorption
- methane
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N5/00—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
- G01N5/02—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by absorbing or adsorbing components of a material and determining change of weight of the adsorbent, e.g. determining moisture content
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/44—Sample treatment involving radiation, e.g. heat
Abstract
本发明提出了一种精确计算甲烷在页岩上真实吸附量的方法,包括:将泥页岩样品粉碎;置于恒温箱中烘干;得到泥页岩样品的质量随气体压力的变化数据,计算出样品和样品缸的初始体积和质量;计算样品在不同压力下的过剩吸附量;将高压段实测样品质量和气体密度在笛卡尔直角坐标系中作图;拟合得到解析式,要求拟合R2=1;计算样品的饱和吸附质量和饱和吸附相体积,二者的比值为吸附相的密度值;对在绝对直线上的饱和吸附点用吸附相体积校正其过剩吸附量,对前面的其他点用吸附相密度校正其过剩吸附量,从而得到甲烷在页岩上的真实吸附量。
Description
技术领域
本发明涉及地质勘探技术领域,尤其是涉及一种精确计算甲烷在页岩上真实吸附量的方法。
背景技术
国内外学者在对煤系泥岩和页岩的吸附性能进行研究时,大多采用容量法等温吸附仪,近几年随着页岩气勘探开发的蓬勃发展,重量法等温吸附仪也正逐渐被国内外学者所接受,后者以磁悬浮天平为主要检测工具,测试结果更精准。
国内外学者在利用容量法对煤系泥岩和页岩等温吸附特征进行研究时发现,在利用容量法对页岩吸附甲烷进行测试时,吸附相体积的存在将减小样品池中自由气体的体积,进而对页岩吸附能力的测试结果产生影响。而在此之前的相关研究中均忽略了吸附相体积存在,这必将导致对于吸附能力的评价存在较大偏差。针对这样的问题,一部分学者试图从不同的角度校正吸附相体积对吸附气能力测试所产生的影响。然而无论以何种方式对吸附相体积所产生的影响进行校正,终将遵循物质守恒原理,因此吸附相密度的确定是无法回避的问题,且吸附相密度的准确性将直接影响校正结果。
由于目前对于等温吸附实验的研究大多采用容量法等温吸附仪,重量法等温吸附仪的使用相对较少,测试手段的单一化导致无法从实验角度对吸附相密度进行精确测试,因此,目前国内外对于吸附相密度的研究仍处于理论假设和数学优化阶段。Yee D.(1994)认为当温度高于临界温度时(通常情况下),吸附相密度等于气体压缩的极限密度,以此推论甲烷的吸附相密度为0.375g/cm3。Arri(1992)假设吸附相密度为常压沸点下的液体密度,以此推论甲烷的吸附相密度为0.421g/cm3。Katsuyuki Murayta(2001)利用经验公式得到甲烷吸附相密度约为0.59g/cm3。以上研究均将吸附相密度作为定值,未考虑温度、压力等条件产生的影响,为此,德国学者Berhard M.Krooss及其团队(2008,2010,2012,2013,2014)将Gibbs关于吸附的定义与兰格缪尔方程结合,将吸附相密度作为待定参数,从数学优化角度对甲烷和二氧化碳在煤系泥岩和页岩上的吸附相密度进行评价,得到不同吸附剂和吸附质组合在不同温度条件下的吸附相密度。类似的,国内学者周理(2000,2001)将Gibbs对吸附的定义式与D-A方程相结合,推导出绝对吸附量模型,并优化出不同温度条件下甲烷在AX-21活性炭上的吸附相密度,温度在233K-333K时甲烷吸附相密度为0.15g/cm3-0.350g/cm3。这种与等温吸附模型相结合进而对吸附相密度进行优化的方法,其对吸附相密度评价结果的准确性将很大程度上依赖于等温吸附模型的适用性,而页岩孔隙结构复杂,等温吸附模型的适用性难以确定,因此利用该方法所得到的评价结果尚待进一步研究。
综上,前人对吸附相密度的研究仍处于理论假设或数学优化阶段,尚未从实验角度对其进行精确测试,而不同学者研究结果相差很大,终无定论,这必将对客观评价岩吸附气能力产生较大影响。
发明内容
本发明的目的是提供一种基于实验方法且适用性强的精确计算甲烷在页岩上真实吸附量的方法。
本发明的技术方案是这样实现的:一种精确计算甲烷在页岩上真实吸附量的方法,包括以下步骤:
S1、将泥页岩样品粉碎至60~80目;
S2、将上述样品置于恒温箱中,设定温度80℃烘干,每隔24h称重一次,直到相邻两次质量变化小于2%为止;
S3、使用氦气对上述处理后泥页岩样品进行浮力测试,得到泥页岩样品的测量质量m测(g)随测试压力P(MPa)变化数据,利用测量装置自带软件分别计算出样品以及样品和样品缸的初始体积(V0样、V0总)和质量(m0样、m0总);
S4、使用甲烷对泥页岩样品进行吸附测试,得到泥页岩样品的测量质量m测(g)随测试压力P(MPa)变化数据。
S5、利用NIST化学数据库查出每个测点在对应温度压力条件下甲烷的自由气密度值ρ气(g/cm3);
S6、计算样品在不同压力点的过剩吸附量A过剩(cm3/g);
S7、将不同压力点对应的实测样品质量m测和甲烷密度ρ气在笛卡尔直角坐标系中作图,其中实测样品质量m测为纵轴,甲烷密度ρ气为横轴;
S8、找到高压段中最后几个在一绝对直线上的点,拟合得到其线性方程,要求拟合R2=1;
S9、上述线性方程的截距减去样品加样品缸的初始质量m0总即得到样品的饱和吸附质量m饱吸(g);上述线性方程的负斜率减去样品加样品缸的初始体积V0总即得到样品的饱和吸附相体积V饱吸(cm3),二者的比值即为吸附相的密度值ρ吸(g/cm3);
S10、对在绝对直线上的饱和吸附点用吸附相体积校正其过剩吸附量,对前面的其他点用吸附相密度校正其过剩吸附量,从而精确得到甲烷在页岩上的真实吸附量A真实。
作为一种优选的技术方案,所述步骤S4中,通过等温吸附实验测量过剩吸附量原理及步骤如下:
①采用磁悬浮等温吸附仪,该吸附仪的顶部为一磁悬浮天平,精度为10μg,天平下部为一密闭的壳体,在壳体内的悬挂端挂一样品缸;当开关置于零点档位时,天平不称量;当置于称量档位时,天平称量样品和样品缸的质量,天平示数等于样品加样品缸的质量;由于磁悬浮天平称量精度高,且不直接与称量的样品部分接触,所以该装置可精确测量出样品缸内的过剩吸附量;
②进行空白测试:将空样品缸放入壳体中,对壳体内空间依次进行加热和抽真空处理,要求加热至样品缸中温度恒定,向壳体内充入氦气,分别在不同压力下进行测量,测量结果满足如下关系式:
m测=m缸-ρhe·V缸
上式中,m测为磁悬浮天平读数,g;m缸为样品缸的质量,g;ρhe为氦气的密度,g/cm3;V缸为样品缸的体积,cm3;在笛卡尔坐标系中将上述测点描点,纵轴为测试质量,横轴为相应压力下氦气的密度,再将点连线进行拟合得到直线关系式,其纵截距即为样品缸的质量m缸,负斜率即为样品缸的体积V缸;
③进行浮力测试:将泥页岩样品放入样品缸,对壳体内空间依次进行加热和抽真空处理,要求加热至样品缸中温度恒定,然后向壳体内充入氦气,并分别在不同压力下进行浮力测试,测量过程满足关系式:
m测=m0总-ρhe·V0总
上式中,m测为磁悬浮天平读数,g;m0总为样品加样品缸的初始质量,g;ρhe为氦气的密度,g/cm3;V0总为样品加样品缸的初始体积,cm3;在笛卡尔坐标系中将上述测点描点,纵轴为天平测试质量,横轴为相应压力下氦气的密度,再将点连线进行拟合得到直线关系式,其纵截距即为样品加样品缸的初始质量m0总,负斜率即为样品加样品缸的初始体积V0总;
④进行吸附测试:对壳体内的进行抽真空处理,待壳体达到测试温度并稳定后,向内充入甲烷至目标压力进行吸附测试,等温吸附实验程序自动采集系统内的温度、压力和天平的读数,待吸附达到稳定后,增加压力继续测量,如此依次加压测试直到设计的压力点测完为止,可得到各个测点天平称量质量与对应的压力数据。
作为一种优选的技术方案,所述步骤S5中,利用NIST化学数据库,直接输入每个测点温度压力条件,即可查询得到甲烷准确的自由气密度值ρ气(g/cm3)。
作为一种优选的技术方案,对所述步骤S6中计算过剩吸附量,对吸附系统受力分析可得到如下关系式:
m测=m0总+m吸附-ρ气·V0总
上式中,m测为磁悬浮天平测量得到的质量,g;m0总为样品和样品缸的初始质量,g;m吸附为样品吸附甲烷的质量,g;ρ气为甲烷在当前气压下的密度,(g/cm3);V0总为样品和样品缸的初始体积,cm3;
由上式可得页岩吸附甲烷质量:
m吸附=m测+ρ气·V0总-m0总
从而可得地面标准条件下过剩吸附量为:
A过剩=(m吸附/M甲烷·22.4·1000)/m0样
其中M甲烷为甲烷相对分子质量,可取16.04276。
作为一种优选的技术方案,所述步骤S8中,找到高压阶段最后几个在一绝对直线上的点,拟合得到其计算式,原则上要求拟合R2=1;现对其数学原理推导如下:
当我们不考虑吸附在页岩表面的吸附相体积时:
m测=m0总+m吸附-ρ气·V0总
但是实际上吸附相的体积是真实存在的,并且会严重影响实验数据,考虑吸附相体积时,有如下关系式:
m测=m0总+m吸附-ρ气·(V0总+V吸附)
上式中V吸附为吸附相的体积,cm3;
当达到吸附饱和后,可以得到:
m测=m0总+m饱吸-ρ气·(V0总+V饱吸)
令a=V0总+V饱吸,b=m0总+m饱吸,y=m测,x=ρ气,既有:
y=ax+b
步骤7)笛卡尔坐标系中所作的点中,可以得到后面几个点应该是能连成一绝对直线的,可选取最后几个直线点拟合,并且从理论上其R2=1。
作为一种优选的技术方案,所述步骤S9中,利用步骤S8拟合计算式的截距减去样品加样品缸的初始质量m0总即得到样品的饱和吸附相质量m饱吸=b-m0总,上述计算式的负斜率减去样品加样品缸的初始体积V0总即得到样品的饱和吸附相体积V饱吸=-a-V0总,二者的比值即为页岩表面吸附相的密度值ρ吸=m饱吸/V饱吸;
在吸附达到饱和前,同样考虑页岩表面吸附相体积的增加,应该有:
m实吸=(m测-m0总)+ρ气·(V0总+V实吸)
上式中,m实吸为样品实际吸附气体的质量,g;V实吸为页岩表面吸附气体的吸附相体积,cm3;
此时,将上式变换可得:
从而有:
作为一种优选的技术方案,所述步骤S10中,对过剩吸附量校正为绝对吸附量,方法如下:
①对在绝对直线上的饱和吸附点用吸附相体积校正其过剩吸附量:
②对前面的其他点用吸附相密度校正其过剩吸附量:
从而精确得到甲烷在页岩上的真实吸附量A真实。
采用了上述技术方案,本发明的有益效果为:本发明能够精确计算甲烷在页岩上真实吸附量的方法,适用性强,克服了困扰行业多年的技术难题,在地质勘探技术领域具有深远的影响。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种精确计算甲烷在页岩上真实吸附量的方法,包括以下步骤:
步骤一、将泥页岩样品粉碎至60~80目。
步骤二、将上述样品置于恒温箱中,设定温度80℃烘干,每隔24h称重一次,直到相邻两次质量变化小于2%为止。
步骤三、对上述处理后泥页岩样品进行重量法等温吸附实验,得到泥页岩样品的测量质量m测(g)随测试压力P(MPa)变化数据,以及利用测量装置自带软件分别计算出样品以及样品和样品缸的初始体积(V0样、V0总)和质量(m0样、m0总)。
通过重量法等温吸附实验测量过剩吸附量原理及步骤如下:
①采用磁悬浮等温吸附仪,该吸附仪的顶部为一磁悬浮天平,精度为10μg,天平下部为一密闭的壳体,在壳体内的悬挂端挂一样品缸;当开关置于零点档位时,天平不称量;当置于称量档位时,天平称量样品和样品缸的质量,天平示数等于样品加样品缸的质量;由于磁悬浮天平称量精度高,且不直接与称量的样品部分接触,所以该装置可精确测量出样品缸内的过剩吸附量;
②进行空白测试:将空样品缸放入壳体中,对壳体内空间依次进行加热和抽真空处理,要求加热至样品缸中温度恒定,向壳体内充入氦气,分别在不同压力下进行测量,测量结果满足如下关系式:
m测=m缸-ρhe·V缸
上式中,m测为磁悬浮天平读数,g;m缸为样品缸的质量,g;ρhe为氦气的密度,g/cm3;V缸为样品缸的体积,cm3;在笛卡尔坐标系中将上述测点描点,纵轴为测试质量,横轴为相应压力下氦气的密度,再将点连线进行拟合得到直线关系式,其纵截距即为样品缸的质量m缸,负斜率即为样品缸的体积V缸;
③进行浮力测试:将泥页岩样品放入样品缸,对壳体内空间依次进行加热和抽真空处理,要求加热至样品缸中温度恒定,然后向壳体内充入氦气,并分别在不同压力下进行浮力测试,测量过程满足关系式:
m测=m0总-ρhe·V0总
上式中,m测为磁悬浮天平读数,g;m0总为样品加样品缸的初始质量,g;ρhe为氦气的密度,g/cm3;V0总为样品加样品缸的初始体积,cm3;在笛卡尔坐标系中将上述测点描点,纵轴为天平测试质量,横轴为相应压力下氦气的密度,再将点连线进行拟合得到直线关系式,其纵截距即为样品加样品缸的初始质量m0总,负斜率即为样品加样品缸的初始体积V0总;
④进行吸附测试:对壳体内的进行抽真空处理,待壳体达到测试温度并稳定后,向内充入甲烷至目标压力进行吸附测试,等温吸附实验程序自动采集系统内的温度、压力和天平的读数,待吸附达到稳定后,增加压力继续测量,如此依次加压测试直到设计的压力点测完为止,可得到各个测点天平称量质量与对应的压力数据。
步骤四、利用NIST化学数据库查出每个测点温度压力条件下甲烷的自由气密度值ρ气(g/cm3)。利用NIST(美国国家标准与技术局)化学数据库,直接输入每个测点温度压力条件,即可查询得到甲烷准确的自由气密度值ρ气(g/cm3)。
步骤五、计算出样品的过剩吸附量A过剩(g/cm3)。对吸附系统受力分析可得到如下关系式:
m测=m0总+m吸附-ρ气·V0总
上式中,m测为磁悬浮天平测量得到的质量,g;m0总为样品和样品缸的初始质量,g;m吸附为样品吸附甲烷的质量,g;ρ气为甲烷在当前气压下的密度,g/cm3;V0总为样品和样品缸的初始体积,cm3;
由上式可得页岩吸附甲烷质量:
m吸附=m测+ρ气·V0总-m0总
从而可得地面标准条件下过剩吸附量为:
A过剩=(m吸附/M甲烷·22.4·1000)/m0样
其中M甲烷为甲烷相对分子质量,可取16.04276。
步骤六、将各个测试点对应的实测样品质量m测和甲烷密度ρ气在笛卡尔直角坐标系中作点,其中实测样品质量m测为纵轴,样品密度ρ气为横轴。
步骤七、找到最后几个在一绝对直线上的点,拟合得到其计算式,要求拟合R2=1。
现对其数学原理推导如下:
当我们不考虑吸附在页岩表面的吸附相体积时:
m测=m0总+m吸附-ρ气·V0总
但是实际上吸附相的体积是真实存在的,并且会严重影响实验数据,考虑吸附相体积时,有如下关系式:
m测=m0总+m吸附-ρ气·(V0总+V吸附)
上式中V吸附为吸附相的体积,cm3;
当达到吸附饱和后,可以得到:
m测=m0总+m饱吸-ρ气·(V0总+V饱吸)
令a=V0总+V饱吸,b=m0总+m饱吸,y=m测,x=ρ气,既有:
y=ax+b
步骤六笛卡尔坐标系中所作的点中,可以得到后面几个点应该是能连成一绝对直线的,可选取最后几个直线点拟合,并且从理论上其R2=1。
步骤八、上述计算式的截距减去样品加样品缸的初始质量m0总即得到样品的饱和吸附相质量m饱吸(g),上述计算式的负斜率减去样品加样品缸的初始体积V0总即得到样品的饱和吸附相体积V饱吸(cm3),二者的比值即为吸附相的密度值ρ吸(g/cm3)。
在吸附达到饱和前,同样考虑页岩表面吸附相体积的增加,应该有:
m实吸=(m测-m0总)+ρ气·(V0总+V实吸)
上式中,m实吸为样品实际吸附气体的质量,g;V实吸为页岩表面吸附气体的吸附相体积,cm3;
此时,将上式变换可得:
从而有:
步骤九、对在绝对直线上的饱和吸附点用吸附相体积校正其过剩吸附量,对前面的其他点用吸附相密度校正其过剩吸附量,从而精确得到甲烷在页岩上的真实吸附量A真实。
对过剩吸附量校正为绝对吸附量,方法如下:
①对在绝对直线上的饱和吸附点用吸附相体积校正其过剩吸附量:
②对前面的其他点用吸附相密度校正其过剩吸附量:
从而精确得到甲烷在页岩上的真实吸附量A真实。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种精确计算甲烷在页岩上真实吸附量的方法,包括以下步骤:S1、将泥页岩样品粉碎至60~80目;其特征在于,还包括:
S2、将上述样品置于恒温箱中,设定温度80℃烘干,每隔24h称重一次,直到相邻两次质量变化小于2%为止;
S3、使用氦气对上述处理后泥页岩样品进行浮力测试,得到泥页岩样品的测量质量m测随测试压力P变化数据,利用测量装置自带软件分别计算出样品以及样品和样品缸的初始体积和质量;
S4、使用甲烷对泥页岩样品进行吸附测试,得到泥页岩样品的测量质量m测随测试压力P变化数据;
S5、利用NIST化学数据库查出每个测点在对应温度压力条件下甲烷的自由气密度值ρ气;
S6、计算样品在不同压力点的过剩吸附量A过剩;
S7、将不同压力点对应的实测样品质量m测和甲烷密度ρ气在笛卡尔直角坐标系中作图,其中实测样品质量m测为纵轴,甲烷密度ρ气为横轴;
S8、找到高压段中最后几个在一绝对直线上的点,拟合得到其线性方程,要求拟合R2=1;
S9、上述线性方程的截距减去样品加样品缸的初始质量m0总即得到样品的饱和吸附质量m饱吸;上述线性方程的负斜率减去样品加样品缸的初始体积V0总即得到样品的饱和吸附相体积V饱吸,二者的比值即为吸附相的密度值ρ吸;
S10、对在绝对直线上的饱和吸附点用吸附相体积校正其过剩吸附量,对前面的其他点用吸附相密度校正其过剩吸附量,从而精确得到甲烷在页岩上的真实吸附量A真实。
2.如权利要求1所述的一种精确计算甲烷在页岩上真实吸附量的方法,其特征在于,所述步骤S4中,通过等温吸附实验测量过剩吸附量原理及步骤如下:
①采用磁悬浮等温吸附仪,该吸附仪的顶部为一磁悬浮天平,精度为10μg,天平下部为一密闭的壳体,在壳体内的悬挂端挂一样品缸;当开关置于零点档位时,天平不称量;当置于称量档位时,天平称量样品和样品缸的质量,天平示数等于样品加样品缸的质量;由于磁悬浮天平称量精度高,且不直接与称量的样品部分接触,所以该装置可精确测量出样品缸内的过剩吸附量;
②进行空白测试:将空样品缸放入壳体中,对壳体内空间依次进行加热和抽真空处理,要求加热至样品缸中温度恒定,向壳体内充入氦气,分别在不同压力下进行测量,测量结果满足如下关系式:
m测=m缸-ρhe·V缸
上式中,m测为磁悬浮天平读数,g;m缸为样品缸的质量,g;ρhe为氦气的密度,g/cm3;V缸为样品缸的体积,cm3;在笛卡尔坐标系中将上述测点描点,纵轴为测试质量,横轴为相应压力下氦气的密度,再将点连线进行拟合得到直线关系式,其纵截距即为样品缸的质量m缸,负斜率即为样品缸的体积V缸;
③进行浮力测试:将泥页岩样品放入样品缸,对壳体内空间依次进行加热和抽真空处理,要求加热至样品缸中温度恒定,然后向壳体内充入氦气,并分别在不同压力下进行浮力测试,测量过程满足关系式:
m测=m0总-ρhe·V0总
上式中,m测为磁悬浮天平读数,g;m0总为样品加样品缸的初始质量,g;ρhe为氦气的密度,g/cm3;V0总为样品加样品缸的初始体积,cm3;在笛卡尔坐标系中将上述测点描点,纵轴为天平测试质量,横轴为相应压力下氦气的密度,再将点连线进行拟合得到直线关系式,其纵截距即为样品加样品缸的初始质量m0总,负斜率即为样品加样品缸的初始体积V0总;
④进行吸附测试:对壳体内的进行抽真空处理,待壳体达到测试温度并稳定后,向内充入甲烷至目标压力进行吸附测试,等温吸附实验程序自动采集系统内的温度、压力和天平的读数,待吸附达到稳定后,增加压力继续测量,如此依次加压测试直到设计的压力点测完为止,可得到各个测点天平称量质量与对应的压力数据。
3.如权利要求1所述的一种精确计算甲烷在页岩上真实吸附量的方法,其特征在于,所述步骤S5中,利用NIST化学数据库,直接输入每个测点温度压力条件,即可查询得到甲烷准确的自由气密度值ρ气。
4.如权利要求1所述的一种精确计算甲烷在页岩上真实吸附量的方法,其特征在于,对所述步骤S6中计算过剩吸附量,对吸附系统受力分析可得到如下关系式:
m测=m0总+m吸附-ρ气·V0总
上式中,m测为磁悬浮天平测量得到的质量,g;m0总为样品和样品缸的初始质量,g;m吸附为样品吸附甲烷的质量,g;ρ气为甲烷在当前气压下的密度;V0总为样品和样品缸的初始体积;
由上式可得页岩吸附甲烷质量:
m吸附=m测+ρ气·V0总-m0总
从而可得地面标准条件下过剩吸附量为:
A过剩=(m吸附/M甲烷·22.4·1000)/m0样
其中M甲烷为甲烷相对分子质量,可取16.04276。
5.如权利要求1所述的一种精确计算甲烷在页岩上真实吸附量的方法,其特征在于,所述步骤S8中,找到高压阶段最后几个在一绝对直线上的点,拟合得到其计算式,原则上要求拟合R2=1;现对其数学原理推导如下:
当不考虑吸附在页岩表面的吸附相体积时:
m测=m0总+m吸附-ρ气·V0总
但是实际上吸附相的体积是真实存在的,并且会严重影响实验数据,考虑吸附相体积时,有如下关系式:
m测=m0总+m吸附-ρ气·(V0总+V吸附)
上式中V吸附为吸附相的体积,cm3;
当达到吸附饱和后,可到:
m测=m0总+m饱吸-ρ气·(V0总+V饱吸)
令a=V0总+V饱吸,b=m0总+m饱吸,y=m测,x=ρ气,即有:
y=ax+b
步骤7)笛卡尔坐标系中所作的点中,可以得到高压段的后几个点应该是在一直线上,选取最后几个直线点拟合,并且从理论上其R2=1。
6.如权利要求5所述的一种精确计算甲烷在页岩上真实吸附量的方法,其特征在于,所述步骤S9中,利用步骤S8拟合计算式的截距减去样品加样品缸的初始质量m0总即得到样品的饱和吸附相质量m饱吸=b-m0总,上述计算式的负斜率减去样品加样品缸的初始体积V0总即得到样品的饱和吸附相体积V饱吸=-a-V0总,二者的比值即为页岩表面吸附相的密度值ρ吸=m饱吸/V饱吸;
在吸附达到饱和前,同样考虑页岩表面吸附相体积的增加,应该有:
m实吸=(m测-m0总)+ρ气·(V0总+V实吸)
上式中,m实吸为样品实际吸附气体的质量,g;V实吸为页岩表面吸附气体的吸附相体积,cm3;
此时,将上式变换可得:
从而有:
7.如权利要求4所述的一种精确计算甲烷在页岩上真实吸附量的方法,其特征在于,所述步骤S10中,对过剩吸附量校正为绝对吸附量,方法如下:
①对在绝对直线上的饱和吸附点用吸附相体积校正其过剩吸附量:
②对前面的其他点用吸附相密度校正其过剩吸附量:
从而精确得到甲烷在页岩上的真实吸附量A真实。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610482070.4A CN106198297B (zh) | 2016-06-24 | 2016-06-24 | 一种精确计算甲烷在页岩上真实吸附量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610482070.4A CN106198297B (zh) | 2016-06-24 | 2016-06-24 | 一种精确计算甲烷在页岩上真实吸附量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106198297A CN106198297A (zh) | 2016-12-07 |
CN106198297B true CN106198297B (zh) | 2019-08-02 |
Family
ID=57462477
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610482070.4A Active CN106198297B (zh) | 2016-06-24 | 2016-06-24 | 一种精确计算甲烷在页岩上真实吸附量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106198297B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106840536B (zh) * | 2017-01-25 | 2018-12-21 | 中国地质大学(北京) | 一种页岩岩样密封性检测装置及其密封性检测方法 |
CN106840957B (zh) * | 2017-04-05 | 2019-08-16 | 中国石油大学(华东) | 一种评价页岩中可动油饱和量的装置及方法 |
CN107271314A (zh) * | 2017-05-24 | 2017-10-20 | 中国矿业大学 | 一种测量煤岩吸附膨胀系数的方法 |
CN108181200B (zh) * | 2017-12-27 | 2021-03-23 | 中国石油化工股份有限公司 | 泥页岩等温吸附实验曲线校正方法及装置 |
CN108446526B (zh) * | 2018-05-22 | 2020-06-26 | 中国地质大学(武汉) | 基于吸附特征曲线与特征方程的页岩吸附量的预测方法 |
CN108593497A (zh) * | 2018-05-30 | 2018-09-28 | 中国工程物理研究院激光聚变研究中心 | 一种精确测量自支撑铍薄膜密度的方法 |
CN109269938A (zh) * | 2018-09-26 | 2019-01-25 | 西安石油大学 | 一种煤和岩石真密度的测试方法 |
CN110220817A (zh) * | 2019-07-08 | 2019-09-10 | 河南理工大学 | 容量法等温吸附实验中吸附量的自由体积校正方法 |
CN111175170A (zh) * | 2020-01-17 | 2020-05-19 | 西安石油大学 | 一种基于重量法的等温吸附曲线校正方法 |
CN111896421B (zh) * | 2020-07-29 | 2021-05-28 | 西南石油大学 | 一种基于吸附势理论计算甲烷在页岩中真实吸附量的方法 |
CN112113881B (zh) * | 2020-09-01 | 2021-05-07 | 中国矿业大学(北京) | 煤中吸附性甲烷的密度测量方法 |
CN112051182B (zh) * | 2020-09-07 | 2022-01-28 | 西南石油大学 | 一种快速预测不同深度页岩储层吸附甲烷能力的方法 |
CN112051179B (zh) * | 2020-09-11 | 2022-08-19 | 中国地质科学院地质力学研究所 | 含水页岩吸附气量测试方法 |
CN112229754B (zh) * | 2020-09-29 | 2022-01-28 | 西南石油大学 | 一种迭代平均法计算甲烷在页岩中吸附量的方法 |
CN112697632B (zh) * | 2020-12-04 | 2023-02-24 | 中国石油天然气股份有限公司 | 一种煤岩、页岩大样量重量法等温吸附测量装置及方法 |
CN112504904B (zh) * | 2020-12-08 | 2021-09-03 | 中国科学院力学研究所 | 一种岩石不同含水条件下气体吸附等温线测试分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104677771A (zh) * | 2015-03-26 | 2015-06-03 | 成都理工大学 | 一种基于磁悬浮重量法确定页岩孔隙度的方法 |
CN104698152A (zh) * | 2015-02-27 | 2015-06-10 | 中国石油大学(华东) | 一种研究泥页岩吸附天然气主控因素的方法 |
CN104713803A (zh) * | 2015-03-16 | 2015-06-17 | 中国石油大学(华东) | 一种精确测量甲烷在页岩上吸附相密度的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8738295B2 (en) * | 2010-05-05 | 2014-05-27 | Conocophillips Company | Shale analysis methods |
-
2016
- 2016-06-24 CN CN201610482070.4A patent/CN106198297B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104698152A (zh) * | 2015-02-27 | 2015-06-10 | 中国石油大学(华东) | 一种研究泥页岩吸附天然气主控因素的方法 |
CN104713803A (zh) * | 2015-03-16 | 2015-06-17 | 中国石油大学(华东) | 一种精确测量甲烷在页岩上吸附相密度的方法 |
CN104677771A (zh) * | 2015-03-26 | 2015-06-03 | 成都理工大学 | 一种基于磁悬浮重量法确定页岩孔隙度的方法 |
Non-Patent Citations (6)
Title |
---|
First international inter-laboratory comparison of high-pressure CH 4 , CO 2 and C 2 H 6 sorption isotherms on carbonaceous shales;M Gasparik等;《International Journal of Coal Geology》;20141231;第1-44页 |
Highest Pressure Adsorption Equilibria Data: Measurement with Magnetic Suspension Balance and Analysis with a New Adsorbent/Adsorbate-Volume;F Dreisbach等;《Adsorption-journal of the International Adsorption society》;20021231;第95-109页 |
富有机质页岩高温高压重量法等温吸附实验;俞凌杰等;《石油学报》;20151231;第557-563页 |
考虑多因素的页岩气吸附能力计算模型;梁彬等;《特种油气藏》;20151231;第121-124页 |
超临界条件下煤层甲烷视吸附量、真实吸附量的差异及其地质意义;杨兆彪等;《天然气工业》;20111231;第13-17页 |
页岩容量法等温吸附实验中异常现象分析;张庆玲;《煤田地质与勘探》;20151231;第31-34页 |
Also Published As
Publication number | Publication date |
---|---|
CN106198297A (zh) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106198297B (zh) | 一种精确计算甲烷在页岩上真实吸附量的方法 | |
CN205138940U (zh) | 一种测试煤体瓦斯吸附解吸热效应实验系统 | |
CN104713803B (zh) | 一种精确测量甲烷在页岩上吸附相密度的方法 | |
CN103994943B (zh) | 一种煤/页岩等温吸附实验装置 | |
CN102901803B (zh) | 受载含瓦斯煤不同含水率吸附-解吸实验的方法 | |
CN103115844B (zh) | 一种煤页岩等温吸附/解吸曲线的测定方法 | |
CN106290061B (zh) | 测定吸附质气体在吸附质上的吸附相密度的方法 | |
CN103018132B (zh) | 煤体吸附解吸过程中动态变形特征测试方法 | |
CN103424421B (zh) | 一种采用低场核磁共振进行煤样甲烷吸附量测量的方法 | |
CN103994960B (zh) | 一种煤/页岩等温吸附实验方法 | |
CN106198306B (zh) | 一种测定吸附质气体在吸附质上的吸附相密度的方法 | |
CN209745750U (zh) | 煤体吸附瓦斯气体过程的变形量-吸附量同步测试装置 | |
CN103033125B (zh) | 煤体吸附解吸瓦斯气体过程变形测试装置 | |
CN108872045A (zh) | 一种页岩碎样总孔隙度的测量方法 | |
CN113011116B (zh) | 一种煤粒微孔扩散系数及无因次数值反演方法 | |
CN103063544A (zh) | 一种模块化煤岩吸附解吸性能评价装置及方法 | |
CN104897514A (zh) | 一种测定煤页岩表面气体吸附及煤页岩解吸曲线的装置 | |
CN106644840A (zh) | 注气驱替煤层瓦斯过程中双向扩散系数的测定装置及方法 | |
CN104729974B (zh) | 一种考虑温度效应的气测孔隙度测量方法 | |
CN201749071U (zh) | 天然气水合物储气量测定的实验装置 | |
CN104267133A (zh) | 一种测定气固表面吸附等温线的方法 | |
Ye et al. | Abnormal adsorption and desorption of nitrogen at 77 K on coals: Study of causes and improved experimental method | |
CN103335927B (zh) | 一种谷物孔隙率测量装置及其测量方法 | |
CN102042942A (zh) | 天然气水合物储气量测定的实验装置 | |
CN106323807A (zh) | 一种通用的密度测定方法及其装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |