一种页岩多组分吸附实验计算方法
【专利摘要】本发明提供了一种页岩多组分吸附实验计算方法，包括：操作等温吸附实验装置，获取参考釜体积Vref和样品釜体积Vsam；操作所述等温吸附实验装置，利用所述参考釜体积Vref和样品釜体积Vsam，计算获取样品釜中装入页岩样品后的自由空间体积Vvoid和中间区体积Vmid；操作所述等温吸附实验装置，利用所述参考釜体积Vref、所述样品釜体积Vsam、所述样品釜中装入页岩样品后的自由空间体积Vvoid和所述中间区体积Vmid，计算获取每个样品釜的压力点所对应的混合气体吸附量以及混合气体各组分的吸附量。
【专利说明】
一种页岩多组分吸附实验计算方法
技术领域
[0001] 本发明涉及吸附量测量领域，特别涉及一种页岩多组分吸附实验计算方法。
【背景技术】
[0002] 目前开展的页岩气吸附实验研究多集中于页岩对甲烷单组分气体的吸附，而实际 页岩储层中吸附的却是混合气体，吸附气体除甲烷以外还包括乙烷、丙烷等烃类气体以及 二氧化碳、氮气等非烃类气体，这些气体在吸附的过程中会相互干扰，因此，单单测量单组 分吸附数据不满足页岩储层实际情况。单组分气体吸附时，无论是在吸附相还是在游离相 中，成分不发生变化;而多组分气体吸附时，由于气体的吸附能力各不相同，每种气体也不 是独立吸附的，而是存在竞争吸附，吸附能力较强的气体会优先占据页岩表面的吸附位，吸 附能力较弱的气体会被从页岩表面置换下来，最终表现出来的宏观现象是吸附相和游离相 中各组分气体的浓度不断发生变化。
[0003] 目前单组分气体吸附实验方法主要有容积法和重量法两大类，其中容积法实验装 置主要由气瓶、真空栗、增压栗、参考釜、样品釜、压力和温度传感器以及恒温箱所组成，实 验过程中，先向参考釜中注入吸附质气体，压力稳定后记录参考釜压力，然后打开参考釜和 样品釜之间的中间阀，让气体等温膨胀进入样品釜，等待吸附平衡后记录平衡压力，根据质 量守恒原理就能计算出每次的吸附增量，累积求和得到对应压力下的吸附量。重量法是根 据吸附前后试样的重量变化来计算吸附量，其主要装置由电磁悬挂天平、高压密闭釜以及 恒温箱组成。其中，密闭釜被分隔为两个密闭室，待测岩样和高压气体分别置于两个密闭室 内，实验过程中，首先记录样品釜的质量读数，然后逐步向样品釜内充入待测气体，等待天 平读数稳定后记录吸附平衡时的样品质量，前后两次质量之差即为吸附气体的质量。
[0004] 与单组分吸附实验不同的是，多组分吸附实验增加气相色谱仪装置，且在每个吸 附平衡点前后，需从参考釜和样品釜中分别取气，采用气相色谱仪测量取出气体组成，进而 得到参考釜和样品釜中游离相混合气体组成。
[0005] 现有的气体吸附装置详细介绍单组分气体吸附实验流程和计算方法的比较多，介 绍多组分气体吸附的实验过程多是一笔带过。因此，需要一个测量多组分吸附的完整实验 方法。
【发明内容】
[0006] 为了能够页岩多组分混合气体吸附实验困难的现状，本发明提供了一种页岩多组 分吸附实验计算方法，给出了吸附实验过程中具体的计算方法，完善了现有技术中的实验 流程。
[0007] 为实现上述目的，本发明提供了一种页岩多组分吸附实验计算方法，包括：
[0008] 操作等温吸附实验装置，获取参考釜体积Vref和样品釜体积Vsam;
[0009] 操作所述等温吸附实验装置，利用所述参考釜体积Vrrf和样品釜体积Vsam，计算获 取样品釜中装入页岩样品后的自由空间体积V TOld和中间区体积Vmid;
[0010] 操作所述等温吸附实验装置，利用所述参考釜体积Vref、所述样品釜体积Vsam、所述 样品釜中装入页岩样品后的自由空间体积V TO1d和所述中间区体积Vmid，获取每个样品釜的 压力点所对应的混合气体吸附量以及混合气体各组分的吸附量。
[0011] 可选的，在本发明一实施例中，所述获取参考爸体积Vref和样品爸体积Vsam时操作 等温吸附实验装置的步骤包括：
[0012] 每次向等温吸附实验装置的样品釜中装入不同体积的不锈钢珠，实验开始前对等 温吸附实验装置抽真空，设置实验温度恒定；
[0013] 向参考釜中充入氦气，待温度压力稳定后记录不同体积的不锈钢珠对应的参考釜 初始压力；
[0014] 将所述样品釜与所述参考釜相连通，氦气等温膨胀进入样品釜，待平衡后记录不 同体积的不锈钢珠对应的平衡压力；
[0015] 根据不同体积的不锈钢珠所对应的初始压力和平衡压力计算得到参考釜体积Vrrf 和样品圣体积Vsam。
[0016] 可选的，在本发明一实施例中，所述计算得到参考釜体积Vrrf和样品釜体积Vsam的 步骤包括：
[0017] 根据He状态方程，利用不同体积的不锈钢珠所对应的初始压力和平衡压力，得到 氦气在不同参考釜初始压力下的密度以及氦气在不同平衡压力下的密度p saml;
[0018] 利用同一体积的不锈钢珠所对应的密度psaml和密度0^{1进行计算，将密度0 3^减 去密度Prrfl之后除以密度Psaml获得的值作为横坐标值，将对应的不锈钢珠的体积AV作为纵 坐标值，在平面直接坐标系中，由所述横坐标值和对应的纵坐标值构成的点拟合成一条直 线，所述直线的斜率为参考釜体积Vrrf，所述直线的截距为样品釜体积乂_。
[0019] 可选的，在本发明一实施例中，所述获取样品釜中装入页岩样品后的自由空间体 积操作等温吸附实验装置的步骤包括：
[0020] 在等温吸附实验装置的样品釜中装入待测页岩样品，实验开始前对等温吸附实验 装置抽真空，设置实验温度恒定；
[0021] 向参考釜中充入氦气，待压力温度稳定记录第二参考釜压力Pref2 ;
[0022] 将所述样品釜与所述参考釜相连通，氦气等温膨胀进入样品釜，待吸附平衡后记 录第二平衡压力Psam2 ;
[0023] 利用所述参考爸体积Vref和样品爸体积Vsam、第二参考爸压力Pref2、第二平衡压力 psam2,计算得到样品釜内自由空间体积。
[0024]可选的，在本发明一实施例中，所述计算得到样品釜内自由空间体积的步骤包括： [0025]根据He状态方程，利用第二参考釜压力Pref 2、第二平衡压力Psam2，得到氦气在第二 参考釜压力Pref2下的密度p ref2以及氦气在第二平衡压力下的密度psam2;
[0026] 根据质量守恒原理，利用所述参考釜体积Vrrf和样品釜体积Vsam、氦气在第二参考 釜压力P rrf2下的密度Prrf2以及氦气在第二平衡压力下的密度psam2，计算出页岩样品骨架体 积Vf rame ；
[0027] 将样品釜体积1_减去页岩样品骨架体积Vframe3获得样品釜中装入页岩样品后的自 由空间体积V?id。
[0028] 可选的，在本发明一实施例中，获取所述中间区体积Vmid时操作等温吸附实验装置 的步骤包括：
[0029] 在等温吸附实验装置的样品釜中装入待测页岩样品，实验开始前对等温吸附实验 装置抽真空，设置实验温度恒定；
[0030] 向参考釜中充入氦气，待压力温度稳定记录第二参考釜压力pref2;将所述样品釜 与所述参考釜相连通，氦气等温膨胀进入样品釜，待吸附平衡后记录第二平衡压力p sam2; [0031 ] 利用所述参考爸体积Vref和样品爸体积Vsam、第二参考爸压力Pref2、第二平衡压力 Psam2，计算得到样品爸内自由空间体积Vvoid ;
[0032] 将等温吸附实验装置的中间区与所述参考釜不相连通，对样品釜抽真空处理之 后，将所述中间区与所述参考釜连通，将所述中间区与所述样品釜置于不连通状态，向参考 釜中充入氦气，待压力温度稳定后记录第三参考釜压力Pref3 ;
[0033] 将所述中间区与所述样品釜连通，氦气等温膨胀进入样品釜，待吸附平衡后记录 第二平衡压力Psam3 ;
[0034] 利用所述参考釜体积Vref和自由空间体积VTOid、第三参考釜压力Pref 3、第三平衡压 力psam3，计算得到中间区体积。
[0035] 可选的，在本发明一实施例中，所述计算得到中间区体积的步骤包括：
[0036]根据He状态方程，利用第三参考釜压力Pref3、第三平衡压力Psam3,得到氦气在第三 参考釜压力Pref3下的密度pref3以及氦气在第三平衡压力下的密度psam3;
[0037] 根据质量守恒原理，利用所述参考釜体积Vref、自由空间体积VTO1d、在第三参考釜 压力p ref3下的密度pref3以及氦气在第三平衡压力下的密度psam3,计算出中间区体积v mid。
[0038] 可选的，在本发明一实施例中，所述获取每个样品釜的压力点所对应的混合气体 吸附量以及混合气体各组分的吸附量时操作等温吸附实验装置的步骤包括：
[0039] 在等温吸附实验装置的样品釜中装入待测页岩样品，实验开始前对等温吸附实验 装置抽真空，设置实验温度恒定；
[0040] 向参考釜中充入混合气体，压力达到预设值之后，将参考釜与中间区相连通，混合 气体等温膨胀进入中间区，待平衡后记录参考釜充气压力为Ρ ιη;
[0041] 将参考釜与中间区置于不相连通的状态，从中间区中取出混合气体，运用气相色 谱仪分析混合气体气源组成浓度yin(j);其中，j代表混合气体中各个组分气体；
[0042] 经过抽真空处理之后让中间区处于真空状态，将中间区和参考釜相连通，混合气 体等温膨胀进入中间区，待平衡后将参考釜与中间区置于不相连通的状态，把中间区与样 品釜相连通，混合气体等温膨胀从中间区进入样品釜，待吸附平衡之后，将样品釜和中间区 置于不相连通的状态；
[0043] 从中间区中取出混合气体，运用气相色谱仪分析吸附平衡后的混合气体组成浓度 yequ+dj);其中，表示第i次吸附平衡后混合气体中第j个组分的摩尔浓度；
[0044]所述参考釜内压力逐步降低，所述样品釜内压力逐步增加，依次反复操作，记录η 个样品釜压力Psmn个参考釜压力PrrfjW及吸附平衡后的混合气体组成浓度yequ+Kj);其 中，i表示混合气体吸附平衡实验中达到吸附平衡的次数，i为自然数，i小于等于n，j代表混 合气体中各个组分气体；
[0045] 利用参考釜充气压力Pin、混合气体气源组成浓度yin(j)、样品釜压力Psam_i、对应的 参考釜压力prefl以及吸附平衡后的混合气体组成浓度yequl(j)，获取第i次吸附平衡时所 对应的混合气体的吸附量、混合气体中第j个组分的吸附量。
[0046]可选的，在本发明一实施例中，所述第i次吸附平衡时所对应的混合气体的吸附量 的表达式为：
[0048] 式中，为测量第i次吸附平衡时混合气体的吸附量，mmol/g;pin为温度为T、压 力为Pir^混合气体气源组成浓度为yin时的气体密度，kg/m3;pref_i为温度T、压力P ref_i时混 合气体气源组成浓度为的气体密度，kg/m3;psam_i为温度T、压力混合气体组成 浓度为yequ时的气体密度，kg/m 3;Min表示混合气体气源组成浓度为yin时的气体摩尔质量， g/mol ;MeqU_i表示混合气体组成浓度为yeqU_i时的气体摩尔质量，g/mol ;Mj为混合气体中第j 个组分的摩尔质量，g/mol;m表示置于样品釜中页岩样品的质量，测量第i次吸
附平衡时等温吸附实验装置总的气量，mmol/g;其中， i = l时， ，
[0049] 可选的，在本发明一实施例中，所述第i次吸附平衡时所对应的混合气体中第j个 组分的吸附量的表达式为：
[0051] 式中，j)为测量第i次吸附平衡后的混合气体中各个气体组分的吸附量。
[0052] 与传统的方法相比，本技术方案的效果和优点体现在以下几个方面：
[0053] 1、本发明第一次完整的阐述了多组分吸附的实验方法和步骤，且实验过程精简、 耗时较短，解决了多组分吸附实验困难的现状；
[0054] 2、本发明中的中间区有两个作用，一是通过中间区的压力过渡作用，实现参考釜 一次注气，样品釜中样品多次吸附，精简了实验流程，减少了实验时间；二是可通过控制中 间区阀门，从中间区取气进行多组分气体摩尔组成分析，取气过程发生在中间区，对样品釜 内平衡压力不会造成波动，使得整个实验过程比较稳定；
[0055] 3、本发明中，预留的较长的取气管线可将从中间区中取出气体的压力降为大气 压，进而直接对取出气体进行气色谱分析，实现了在线色谱测量，省去了人工取气的繁琐步 骤，且避免了人工取气过程中带入杂质气体。
【附图说明】
[0056] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。
[0057] 图1为本发明提出的一种页岩多组分吸附实验计算方法流程图；
[0058]图2为实施例的实验装置示意图。
[0059] 附图标识：
[0060] 1-混合气体钢瓶，2-氦气钢瓶，3-气体增压栗，4-缓冲容器，5-五通阀，6-排气口， 7-真空栗，8-参考釜，9-样品釜，10-气相色谱仪，11-恒温油浴，12-温度传感器，13-第一压 力传感器，14-第二压力传感器，15-数据转换模块，16-计算机数据采集模块，17-第一手动 阀，18-第二手动阀，19-调压阀，20-第一气动阀，21-第二气动阀，22-第三气动阀，23-第四 气动阀，24-取气阀。
【具体实施方式】
[0061] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0062] 根据【背景技术】部分记载可知，由于完成整个吸附实验所需要的气体量较多，因此 一次性注入足够量的吸附质气体后，参考釜的初始压力很高，而样品釜的初始压力近似为 真空，吸附实验过程中参考釜需要逐步降压，在现有技术的情况下很难通过阀门进行准确 控制，为了实现参考釜的逐步降压过程，本技术方案在实验装置上设置了中间区，利用中间 区来进行逐步降压的方法。
[0063]具体技术思路为:在参考釜和实验釜中间，设置一个体积很小（3ml~5ml)的中间 区，由于中间区体积较参考釜体积小很多，因此，即使在高压情况下，样品釜每次的压力增 加也是一定的，可以通过阀门控制中间区的气体向样品釜内逐步增压。这样，在本技术方案 中，中间区具有以下两个方面的作用：
[0064] 1、利用中间区较小体积来实现参考釜的逐步降压、样品釜的逐步增压过程，该方 法简单方便、压力可控制性强;并且，在吸附实验过程中，进行单组分气体实验时，由于通过 中间区向参考釜一次性注入足够量的吸附质气体，并且通过中间区对参考釜进行逐步降低 压力，对样品釜进行逐步增加压力，得到多个不同吸附平衡后样品釜内压力下对应的吸附 量。传统技术假如要得到5个不同吸附平衡后样品釜内压力下对应的吸附量，就需要先后5 次向参考釜内注入吸附质气体，使得测量结果的实验误差得到累积，所以本发明的技术方 案能够降低测量结果的实验误差。
[0065] 2、在吸附实验过程中，当进行多组分气体实验时，在吸附平衡后，利用中间区取 气，并通过气相色谱仪进行气体组分分析。
[0066] 基于上述分析，本发明提出一种页岩多组分吸附实验计算方法，如图1所示。包括：
[0067] 步骤101):操作等温吸附实验装置，获取参考釜体积Vref和样品釜体积V sam;
[0068] 步骤102):操作所述等温吸附实验装置，利用所述参考釜体积Vref和样品釜体积 Vsam，计算获取样品釜中装入页岩样品后的自由空间体积V TOld和中间区体积Vmid;
[0069] 步骤103):操作所述等温吸附实验装置，利用所述参考釜体积Vref、所述样品釜体 积Vsam、所述样品釜中装入页岩样品后的自由空间体积V TOld和所述中间区体积Vmid，计算获 取每个样品釜的压力点所对应的混合气体吸附量以及混合气体各组分的吸附量。
[0070] 在步骤101中，操作等温吸附实验装置的步骤包括：
[0071]每次向等温吸附实验装置的样品釜中装入不同体积的不锈钢珠，实验开始前对等 温吸附实验装置抽真空，设置实验温度恒定；
[0072] 向参考釜中充入氦气，待温度压力稳定后记录不同体积的不锈钢珠对应的参考釜 初始压力；
[0073] 将所述样品釜与所述参考釜相连通，氦气等温膨胀进入样品釜，待平衡后记录不 同体积的不锈钢珠对应的平衡压力；
[0074] 根据不同体积的不锈钢珠所对应的初始压力和平衡压力计算得到参考釜体积Vrrf 和样品釜体积Vsam。其中，计算得到参考釜体积Vrrf和样品釜体积V_的步骤包括：
[0075] 根据He状态方程，利用不同体积的不锈钢珠所对应的初始压力和平衡压力，得到 氦气在不同参考釜初始压力下的密度以及氦气在不同平衡压力下的密度p saml;
[0076] 利用同一体积的不锈钢珠所对应的密度psaml和密度0^{1进行计算，将密度03^减 去密度Prrfl之后除以密度Psaml获得的值作为横坐标值，将对应的不锈钢珠的体积AV作为纵 坐标值，在平面直接坐标系中，由所述横坐标值和对应的纵坐标值构成的点拟合成一条直 线，所述直线的斜率为参考釜体积Vrrf，所述直线的截距为样品釜体积乂_。
[0077] 在步骤102中，操作等温吸附实验装置的步骤包括：
[0078] 在等温吸附实验装置的样品釜中装入待测页岩样品，实验开始前对等温吸附实验 装置抽真空，设置实验温度恒定；
[0079]向参考釜中充入氦气，待压力温度稳定记录第二参考釜压力Pref2 ;
[0080]将所述样品釜与所述参考釜相连通，氦气等温膨胀进入样品釜，待吸附平衡后记 录第二平衡压力Psam2 ;
[0081 ] 利用所述参考釜体积Vrrf和样品釜体积Vsam、第二参考釜压力Prrf2、第二平衡压力 psam2，计算得到样品釜内自由空间体积。其中，计算得到样品釜内自由空间体积的步骤包 括：
[0082] 根据He状态方程，利用第二参考釜压力Pref 2、第二平衡压力Psam2，得到氦气在第二 参考釜压力Pref2下的密度p ref2以及氦气在第二平衡压力下的密度psam2;
[0083] 根据质量守恒原理，利用所述参考釜体积Vref和样品釜体积V_、氦气在第二参考 釜压力Prrf2下的密度Prrf 2以及氦气在第二平衡压力下的密度Psam2，计算出页岩样品骨架体 积Vf rame ；
[0084] 将样品釜体积1_减去页岩样品骨架体积Vframe3获得样品釜中装入页岩样品后的自 由空间体积V?id。
[0085] 在步骤102中，在获得样品釜中装入页岩样品后的自由空间体积VTOld时操作等温 吸附实验装置的步骤的基础上还包括：
[0086] 将等温吸附实验装置的中间区与所述参考釜不相连通，对样品釜抽真空处理之 后，将所述中间区与所述参考釜连通，将所述中间区与所述样品釜置于不连通，向参考釜中 充入氦气，待压力温度稳定后记录第三参考釜压力Pref3 ;
[0087] 将所述中间区与所述样品釜相连通，氦气等温膨胀进入样品釜，待吸附平衡后记 录第三平衡压力Psam3;
[0088] 利用所述参考釜体积Vref和自由空间体积VTOid、第三参考釜压力Pref 3、第三平衡压 力psam3，计算得到中间区体积。其中，计算得到中间区体积的步骤包括：
[0089] 根据He状态方程，利用第三参考釜压力Pref3、第三平衡压力Psam3,得到氦气在第三 参考釜压力Pref3下的密度pref3以及氦气在第三平衡压力下的密度psam3;
[0090]根据质量守恒原理，利用所述参考釜体积Vref和自由空间体积Vvoid、在第三参考釜 压力pref3下的密度pref3以及氦气在第三平衡压力下的密度psam3,计算出中间区体积v mid。 [0091 ]在步骤103中，操作等温吸附实验装置的步骤包括：
[0092] 在等温吸附实验装置的样品釜中装入待测页岩样品，实验开始前对等温吸附实验 装置抽真空，设置实验温度恒定；
[0093] 向参考釜中充入混合气体，压力达到预设值之后，将参考釜与中间区相连通，混合 气体等温膨胀进入中间区，待平衡后记录参考釜充气压力为Ρ ιη;
[0094] 将参考釜与中间区置于不相连通的状态，从中间区中取出混合气体，运用气相色 谱仪分析混合气体气源组成浓度yin(j);其中，j代表混合气体中各个组分气体；
[0095] 经过抽真空处理之后让中间区处于真空状态，将中间区和参考釜相连通，混合气 体等温膨胀进入中间区，待平衡后将参考釜与中间区置于不相连通的状态，把中间区与样 品釜相连通，混合气体等温膨胀从中间区进入样品釜，待吸附平衡之后，将样品釜和中间区 置于不相连通的状态；
[0096]从中间区中取出混合气体，运用气相色谱仪分析吸附平衡后的混合气体组成浓度 yequ+dj);其中，表示第i次吸附平衡后混合气体中第j个组分的摩尔浓度；
[0097]所述参考釜内压力逐步降低，所述样品釜内压力逐步增加，依次反复操作，记录η 个样品釜压力匕^:及对应的η个参考釜压力Pref_1;其中，i表示混合气体吸附平衡实验中达 到吸附平衡的次数，i为自然数，i小于等于η;
[0098] 利用参考釜充气压力Pin、混合气体气源组成浓度yin(j)、样品釜压力Psam_i及对应 的参考釜压力prefl，获取第i次吸附平衡时所对应的混合气体的吸附量、混合气体中第j个 组分的吸附量。其中，
[0099] 所述第i次吸附平衡时所对应的混合气体的吸附量的表达式为：
[0101] 式中，为测量第i次吸附平衡时混合气体的吸附量，mmol/g;pin为温度为T、压 力为Pir^混合气体气源组成浓度为yin时的气体密度，kg/m3;pref_i为温度T、压力P ref_i时混 合气体气源组成浓度为的气体密度，kg/m3;psam_i为温度T、压力混合气体组成 浓度为yequ时的气体密度，kg/m 3;Min表示混合气体气源组成浓度为yin时的气体摩尔质量， g/mol ;MeqU_i表示混合气体组成浓度为yeqU_i时的气体摩尔质量，g/mol ;Mj为混合气体中第j 个组分的摩尔质量，g/mol;m表示置于样品釜中页岩样品的质量，测量第i次吸 附平衡时等温吸附实验装置总的气量，mmol/g;其中
_i = l时，
[0102]所述第i次吸附平衡时所对应的混合气体中第j个组分的吸附量的表达式为：
[0104] 式中，j)为测量第i次吸附平衡后的混合气体中各个气体组分的吸附量。
[0105] 以80%CH4+20%⑶2混合气体在页岩上的吸附实验，对本发明作详细介绍。但这些 实例不能被理解为本发明可实施范围的限定。
[0106] 在本发明的实施例中，采用如图2所示的实验装置。实验的具体过程如下：
[0107] 首先进行参考釜和样品釜体积的计算。
[0108] 向样品釜9中装入数个不锈钢钢珠，使用真空栗7对整个实验装置抽真空，所有阀 门处于关闭状态，设置油浴11温度为30°C。打开第一气动阀20和第二手动阀18,氦气钢瓶2 中的氦气在气体增压栗3的作用下经过缓冲容器4到达参考釜8,关闭第一气动阀20,通过温 度传感器12和第一压力传感器13监测参考釜8的压力和温度，待压力温度稳定记录参考釜8 压力P refl。打开第二气动阀21和第三气动阀22,氦气等温膨胀进入样品釜9,通过第二压力 传感器14监测样品釜的压力，待吸附平衡后记录平衡压力为P saml。改变样品釜9中的钢珠 数，重复上述步骤，如此循环两次，最终得到三组参考釜压力和平衡压力。根据质量守恒定 律，有以下公式：
[01 09] PreflVref - Psaml ( Vref+Vsam_ A V) ( 1 )
[01 ?Ο] 式中，Prefl、Psaml分别为氦气在压力Prefl和压力Psaml下的密度，kg/m3，通过氦气状态 方程求得;Vrrf、Vsam分别为参考釜8的体积、样品釜9的体积，cm3;AV为装入样品釜9中的样 品体积。其中，当测量样品釜9和参考釜8体积时，该体积为装入样品釜的不锈钢珠的体积， 当测量自由空间体积和中间体积时，该体积为装入样品釜的页岩样品骨架体积，cm 3。
[0111] 将方程(1)整理成以下形式：
(:2)
[0113] 方程(2)为直线方程y = kx+b的形式。其中，参考釜8体积Vref为直线斜率k，样品釜9 体积Vsam为直线截距b。通过改变样品釜9中钢珠数，重复氦气膨胀实验，循环三次，得到三组 压力值，利用He状态方程计算得到三组p refl、psaml值，进而将三组AV、prefl、psaml值按照方程 (2)的形式拟合成一条直线，根据直线的斜率和截距可得参考釜和样品釜体积。
[0114] 由图2可知，在缓冲容器4的气体输入端设置有调压阀19,该阀门主要用于控制增 压栗3作用下的气体压强在40MPa以内。增长实验装置的寿命。混合气体钢瓶1和氦气钢瓶2 输出的气体在气体增压栗3的作用下经过缓冲容器4后才到参考釜8。这样避免输入的气体 压强过大对参考釜8造成损坏。温度传感器12、第一压力传感器13和第二压力传感器14监测 的数据通过数据转换模块15处理，转换成被数据采集模块16所接受格式的数据。
[0115]下面进行自由空间体积和中间区体积标定。
[0116]在样品釜9中装入页岩样品（样品经过干燥箱110°C恒温干燥12小时），使用真空栗 7对系统抽真空，所有阀门处于关闭状态，设置油浴11温度为30°C。打开第二手动阀18和第 一气动阀20,氦气钢瓶2中氦气在气体增加栗的作用下经过缓冲容器4到达参考釜8。关闭第 一气动阀20，待压力温度稳定记录参考釜压力P ref2，打开第二气动阀21和第三气动阀22，氦 气等温膨胀进入样品釜9,待压力温度稳定后记录平衡压力为P sam2，通过以下公式可计算出 页岩样品骨架体积Vfran? : (3)
[01 18] 式中，Vframe为页岩样品骨架体积，〇113;01^2、0331112分别为氦气在压力?1^2(单位]\ /〇：)&) 和卩_2 (单位MPa)下的密度，kg/m3。
[0119] 则样品釜中装入页岩样品后的自由空间体积为：
[0120] Vvoid = Vsam-Vframe (4)
[0121 ] 式中，V?id为自由空间体积，cm3。
[0122]紧接着关闭第二气动阀21，打开第三气动阀22和第四气动阀23,使用真空栗7对样 品釜9抽真空，关闭第三气动阀22和第四气动阀23，打开第二气动阀21，并同时打开第一气 动阀20和第二手动阀18，在0-10MPa的范围内再次向参考釜中充入氦气，在气体增加栗3的 作用下，氦气经过缓冲容器4到达参考釜8,关闭第一气动阀20,待压力温度稳定记录参考釜 压力P ref3,打开第三气动阀22,氦气等温膨胀进入样品釜，待吸附平衡后记录平衡压力为 Psam3，根据质量守恒定律，可得以下等式：
[01 23] Pref 3 ( Vref+Vmid )-Psam3 (Vref+Vvoid) (5)
[0124]整理可得到计算中间体积的表达式：
(6)
[01 26] 上式中，Pref 3、Psam3分别为気气在压力Pref 3 (单位MPa)和Psam3(单位MPa)下的密度， kg/m3; Vmid为中间体积，cm3。
[0127] 下面进行混合气体吸附的测量。
[0128] 使用真空栗7对实验系统抽真空，所有阀门关闭，设置油浴温度为恒定值T。打开第 一气动阀20和第一手动阀17,混合气体钢瓶1中的混合气体在气体增压栗3的作用下经过缓 冲容器4到达参考釜8，压力达到预设值后，关闭第一气动阀20，待压力温度稳定后，打开第 二气动阀21，混合气体进入中间区，待平衡后记录参考釜充气压力为Ρ ιη。该中间区为第二气 动阀21、第三气动阀22和第四气动阀23均关闭之后由装置的管线构成。中间区的体积远小 于参考釜的体积。关闭第二气动阀21，打开第四气动阀23和取气阀24,从中间区中取出混合 气，运用气相色谱仪10分析混合气体气源组成浓度y in(j)(j代表混合气体中各个组分气 体），真空栗7利用五通阀5与缓冲容器4、参考釜8、中间区通过管线相连通，对中间区和取气 管线抽真空后关闭第四气动阀23。打开第二气动阀21，参考釜8中的混合气体进入中间区， 关闭第二气动阀21，打开第三气动阀22,中间区中气体进入样品釜9,待吸附平衡后，记录参 考釜压力为IW」，样品釜压力为P_j。关闭第三气动阀22,再次从中间区中取气，运用气相 色谱仪分析吸附平衡后混合气体组成浓度y equ」(j)(j代表混合气体中各个组分气体）。
[0129] 第1次吸附平衡之后开始测试，实验装置的总气量为： (7)
[0131] 其中，
[0132] Min (8)
[0133] 第1次吸附平衡之后混合气体吸附量： (69:)
[0135] 其中， 2.
[0136] Mequ x^yequ l(j)Mi{W)
[0137] 第1次吸附平衡时混合气体的第j个组分的吸附量为：
(11)
[0139]由于测量第1次吸附平衡时，要取出部分气体进行吸附平衡后混合气组分分析，因 此第二个吸附点开始测试时，装置的总气量要在原始总气量的基础上减去取出的这部分气 体：
(12)
[0141]依次类推，测量第η个吸附平衡点时，装置的总气量为：
(13) L0143」第η次吸附平衡时混合气体吸附量为：
(14)
[0145]第η次吸附平衡时混合气体的第j个组分的吸附量为：
(1:5)
[0147] 以上各个式子中，测量第i次吸附平衡时实验装置总气量， 为测量第i次吸附平衡时混合气体的吸附量，为测量第i次吸附平衡时混 合气体中各个气体组分的吸附量，j代表混合气体中各个气体组分，mm〇l/g; yin( j)为混合气 体中各个气体组分的摩尔浓度，j代表混合气体中各个气体组分;yequl(j)为第i次吸附平衡 后混合气体中各个气体组分的摩尔浓度，j代表混合气体中各个气体组分;ρ ιη为温度T、压力 Pin时混合气体组成浓度为yin时的混合气体密度，kg/m3;Pref_i为温度T、压力Pref_i时混合气 体组成浓度为y in时的混合气体密度，kg/m3;psam_i为温度T、压力Psam_i时混合气体组成浓度 为y equ_i时的混合气体密度，kg/m3; Min是混合气体组成浓度为yin时的混合气体摩尔质量，g/ mol ;MeqU_i是混合气体组成浓度为yeqU_i时的混合气体摩尔质量，g/mol ;Mj为混合气体中第j 个组分的摩尔质量，g/mol。!!!表示置于样品釜中页岩样品的质量，gd表示混合气体吸附平 衡实验中达到吸附平衡的次数，i为自然数，i小于等于η。
[0148] 由图2可知，五通阀5上有一个排气口 6，用于实验完成后，将实验装置内残余的气 体排出，这样就彻底完成了吸附平衡实验。
[0149] 最后应说明的是:上述仅用以说明本发明而并非限制本发明所描述的技术方案； 尽管本说明书对本发明已进行了详细的说明，但是，本领域的技术人员仍然可以对本发明 进行修改或等同替换，一切不脱离本发明的精神和范围的技术方案及其改进，其均应涵盖 在本发明的权利要求范围中。
【主权项】
1. 一种页岩多组分吸附实验计算方法，其特征在于，包括： 操作等温吸附实验装置，获取参考釜体积Vrrf和样品釜体积Vsam; 操作所述等温吸附实验装置，利用所述参考釜体积Vrrf和样品釜体积Vsam，计算获取样 品釜中装入页岩样品后的自由空间体积Vvoid和中间区体积Vmid; 操作所述等温吸附实验装置，利用所述参考釜体积Vrrf、所述样品釜体积Vsam、所述样品 釜中装入页岩样品后的自由空间体积Vvcild和所述中间区体积Vmid，获取每个样品釜的压力 点所对应的混合气体吸附量以及混合气体各组分的吸附量。2. 如权利要求1所述的计算方法，其特征在于，所述获取参考釜体积Vre3f和样品釜体积 Vsam时操作等温吸附实验装置的步骤包括： 每次向等温吸附实验装置的样品釜中装入不同体积的不锈钢珠，实验开始前对等温吸 附实验装置抽真空，设置实验温度恒定； 向参考釜中充入氦气，待温度压力稳定后记录不同体积的不锈钢珠对应的参考釜初始 压力； 将所述样品釜与所述参考釜相连通，氦气等温膨胀进入样品釜，待平衡后记录不同体 积的不锈钢珠对应的平衡压力； 根据不同体积的不锈钢珠所对应的初始压力和平衡压力计算得到参考釜体积Vrrf和样 品圣体积Vsam。3. 如权利要求2所述的计算方法，其特征在于，所述计算得到参考釜体积Vrrf和样品釜 体积Vsam的步骤包括： 根据He状态方程，利用不同体积的不锈钢珠所对应的初始压力和平衡压力，得到氦气 在不同参考爸初始压力下的密度Prefl以及氦气在不同平衡压力下的密度Psaml; 利用同一体积的不锈钢珠所对应的密度Psaml和密度仏^进行计算，将密度Psaml减去密 度Prrfl之后除以密度Psaml获得的值作为横坐标值，将对应的不锈钢珠的体积AV作为纵坐标 值，在平面直接坐标系中，由所述横坐标值和对应的纵坐标值构成的点拟合成一条直线，所 述直线的斜率为参考釜体积Vrrf，所述直线的截距为样品釜体积乂_。4. 如权利要求1所述的计算方法，其特征在于，所述获取样品釜中装入页岩样品后的自 由空间体积VvcilJt操作等温吸附实验装置的步骤包括： 在等温吸附实验装置的样品釜中装入待测页岩样品，实验开始前对等温吸附实验装置 抽真空，设置实验温度恒定； 向参考釜中充入氦气，待压力温度稳定记录第二参考釜压力Pref2 ; 将所述样品釜与所述参考釜相连通，氦气等温膨胀进入样品釜，待吸附平衡后记录第 二平衡压力Psam2 ; 利用所述参考爸体积Vref和样品爸体积Vsam、第二参考爸压力Pref 2、第二平衡压力Psam2， 计算得到样品釜内自由空间体积。5. 如权利要求4所述的计算方法，其特征在于，所述计算得到样品釜内自由空间体积的 步骤包括： 根据He状态方程，利用第二参考釜压力Prrf2、第二平衡压力Psam2，得到氦气在第二参考 釜压力Pref2下的密度Pref2以及氦气在第二平衡压力下的密度Psam2; 根据质量守恒原理，利用所述参考釜体积Vrrf和样品釜体积Vsam、氦气在第二参考釜压 力Prrf2下的密度Prrf2以及氦气在第二平衡压力下的密度Psam2，计算出页岩样品骨架体积 Vframe； 将样品釜体积Vsam减去页岩样品骨架体积Vframe3获得样品釜中装入页岩样品后的自由空 间体积Vvoid。6.如权利要求1所述的计算方法，其特征在于，获取所述中间区体积Vmid时操作等温吸 附实验装置的步骤包括： 在等温吸附实验装置的样品釜中装入待测页岩样品，实验开始前对等温吸附实验装置 抽真空，设置实验温度恒定； 向参考釜中充入氦气，待压力温度稳定记录第二参考釜压力Pref 2;将所述样品釜与所述 参考釜相连通，氦气等温膨胀进入样品釜，待吸附平衡后记录第二平衡压力Psam2; 利用所述参考爸体积Vref和样品爸体积Vsam、第二参考爸压力Pref 2、第二平衡压力Psam2， 计算得到样品釜内自由空间体积VTO1d; 将等温吸附实验装置的中间区与所述参考釜不相连通，对样品釜抽真空处理之后，将 所述中间区与所述参考釜连通，将所述中间区与所述样品釜置于不连通状态，向参考釜中 充入氦气，待压力温度稳定后记录第三参考釜压力Pref 3 ; 将所述中间区与所述样品釜连通，氦气等温膨胀进入样品釜，待吸附平衡后记录第三 平衡压力Psam3 ; 利用所述参考釜体积Vref和自由空间体积Vvoid、第三参考釜压力Pref3、第三平衡压力 P_3，计算得到中间区体积。 7 .如权利要求6所述的计算方法，其特征在于，所述计算得到中间区体积的步骤包括： 根据He状态方程，利用第三参考釜压力P ref 3、第三平衡压力Psam3，得到氦气在第三参考 爸压力Pref3下的密度Pref3以及氦气在第三平衡压力下的密度Psam3 ; 根据质量守恒原理，利用所述参考釜体积Vref、自由空间体积Vvcild、在第三参考釜压力 Pref3下的密度Pref3以及氦气在第三平衡压力下的密度PSam3,计算出中间区体积V mld。8.如权利要求1所述的计算方法，其特征在于，所述获取每个样品釜的压力点所对应的 混合气体吸附量以及混合气体各组分的吸附量时操作等温吸附实验装置的步骤包括： 在等温吸附实验装置的样品釜中装入待测页岩样品，实验开始前对等温吸附实验装置 抽真空，设置实验温度恒定； 向参考釜中充入混合气体，压力达到预设值之后，将参考釜与中间区相连通，混合气体 等温膨胀进入中间区，待平衡后记录参考釜充气压力为Pin; 将参考釜与中间区置于不相连通的状态，从中间区中取出混合气体，运用气相色谱仪 分析混合气体气源组成浓度yin(j);其中，j代表混合气体中各个组分气体； 经过抽真空处理之后让中间区处于真空状态，将中间区和参考釜相连通，混合气体等 温膨胀进入中间区，待平衡后将参考釜与中间区置于不相连通的状态，把中间区与样品釜 相连通，混合气体等温膨胀从中间区进入样品釜，待吸附平衡之后，将样品釜和中间区置于 不相连通的状态； 从中间区中取出混合气体，运用气相色谱仪分析吸附平衡后的混合气体组成浓度yen (j);其中，yequ+dj)表示第i次吸附平衡后混合气体中第j个组分的摩尔浓度； 所述参考釜内压力逐步降低，所述样品釜内压力逐步增加，依次反复操作，记录n个样 品釜压力Psmii个参考釜压力Pm1W及吸附平衡后的混合气体组成浓度y^U);其中，i 表示混合气体吸附平衡实验中达到吸附平衡的次数，i为自然数，i小于等于η，j代表混合气 体中各个组分气体； 利用参考釜充气压力Pin、混合气体气源组成浓度yin(j)、样品釜压力Psam_i、对应的参考 釜压力及吸附平衡后的混合气体组成浓度获取第i次吸附平衡时所对应 的混合气体的吸附量、混合气体中第j个组分的吸附量。9. 如权利要求8所述的计算方法，其特征在于，所述第i次吸附平衡时所对应的混合气 体的吸附量的表达式为：式中，nmix_i为测量第i次吸附平衡时混合气体的吸附量，mm〇l/g;pir^温度为T、压力为P ir^ 混合气体组成浓度为yir^的气体密度，kg/m3;pref_i为温度T、压力Pref_i时混合气体组成浓度为 7』寸的气体密度，kg/m 3;pS￡m_i为温度T、压力Psa^时混合气体组成浓度为yequ_i时的气体密度，kg/ m3;Min表不混合气体组成浓度为yin时的气体摩尔质量，g/m〇l ;Mequ_i表不混合气体组成浓度为 yequ_i时的气体摩尔质量，g/mo I; Mj为混合气体中第j个组分的摩尔质量，g/mo I; m表不置于样品 釜中页岩样品的质量，测量第i次吸附平衡时等温吸附实验装置总的气量，mmol/10. 如权利要求8所述的计算方法，其特征在于，所述第i次吸附平衡时所对应的混合气 体中第j个组分的吸附量的表达式为：式中，nsingie_i (j)为测量苐i次吸附肀衡后的混合气体中谷个气体组分的吸附量。
【文档编号】G01N30/02GK106018165SQ201610602168
【公开日】2016年10月12日
【申请日】2016年7月27日
【发明人】齐荣荣, 王庆, 宁正福, 张宇琪, 杜华明
【申请人】中国石油大学(北京)