一种地层压力下致密岩石孔渗联测装置及方法

1.本发明涉及对岩石孔隙度、渗透率测量的技术领域，特别涉及一种地层压力下致密岩石孔渗联测装置及方法。


背景技术：

2.储层岩石的孔隙性和渗透性是认识与研究储层的重要性质。孔隙度是定量描述岩石储集能力大小的重要参数，它是指岩石中的孔隙体积(或岩石中未被固体物质充填的空间体积)与岩石总体积的比值。渗透率表示在一定压差下的岩石允许流体通过的性质，储层渗透性的大小用渗透率来表示。随着地下油气的开采，油气流体将从地下采出，储层岩石中的孔隙流体压力就会降低，储层岩石就会变形，这导致储层岩石的孔隙度和渗透率不再是某一固定的常数，而是会发生改变；储层岩石孔隙度和渗透率在地层条件下随孔隙流体压力变化关系的定量表征是油气有效开采的重要基础。
3.目前国内外学者主要通过实验来确定地层条件下岩石孔隙度和渗透率(后面简称“孔渗”)随孔隙流体压力的变化关系式。一般地，首先测定当孔隙流体压力等于大气压时，在不同围压下的孔渗；然后，计算获取孔渗随有效压力的关系式；最后，结合terzaghi有效应力理论计算获取等效地层条件下的孔渗值。然而，现有的孔隙度和渗透率测试仍存在以下几点不足：
4.(1)现有方法没有直接测定地层条件下储层岩石的孔隙度和渗透率，属于储层岩石孔渗确定的间接方法，已有学者证明现有方法计算获取的等效孔渗值不具有地层条件下代表性(肖文联等，2019)。
5.(2)目前已有的岩样孔隙度、渗透率测试设备多数需要将孔隙度及渗透率分开进行测试，或者是同一设备中分步进行测试，并没有实现同时、准确的测定；针对地下取出岩石中可能会存在微裂缝等不稳定的因素，分开测试或分步测试难以保证测试时岩石性质的一致性，从而导致测试的岩石孔渗不在相同条件。
6.(3)实验过程中，测试方案的不同、设备维修等都可能引起管线及阀门体积变化，这将会影响测试结果，往往需要及时进行死体积标定和校正。
7.(4)致密岩石的渗透率很低，一般不能用常规的稳态法(何更生等，2011)来测定其渗透率，通常会选择压力脉冲法来测定其渗透率；目前压力脉冲法在使用时，测试仪器中的压力室体积一般多为固定体积或者依靠经验进行选择，然而实际可变压力室体积的大小会影响脉冲压力衰减曲线，造成测试结果误差。
8.因此，针对现有测试装置及方法的不足，有必要设计一种地层压力下致密岩石孔渗联测装置及方法。


技术实现要素：

9.针对上述问题，本发明旨在提供一种地层压力下致密岩石孔渗联测装置及方法，通过标定管线死体积与优选可变体积压力室体积，绘制脉冲压力衰减曲线，优化孔隙度与
渗透率计算方法，结合压力脉冲实验实现对致密岩石在地层压力下孔隙度及渗透率的同时、有效测定。
10.本发明的技术方案如下：
11.一方面，提供一种地层压力下致密岩石孔渗联测装置，包括高压气源、可变体积压力室、岩心夹持器、围压泵、真空泵、恒温箱、数据采集装置；
12.所述高压气源与所述可变体积压力室的入口端相连，且相连的管路上依次设有阀门一、压力控制器一；
13.所述可变体积压力室与压力监测装置一相连，所述可变体积压力室的出口端与所述岩心夹持器的入口端相连，且相连的管路上设有阀门二；
14.所述岩心夹持器的侧壁与所述围压泵相连，且相连的管路上依次设有压力监测装置二、压力控制器二、阀门三；所述岩心夹持器与温度监测装置相连；所述岩心夹持器的入口端与出口端之间设有压差监测装置；所述岩心夹持器的出口端与所述真空泵相连，且相连的管路上依次设有压力监测装置三和三通阀，所述真空泵上设有压力监测装置四，所述三通阀的另一个出口与放空管相连，所述放空管上设有阀门四；
15.所述可变体积压力室与所述岩心夹持器设置在所述恒温箱内；所述压力控制器一、压力监测装置一、压力监测装置二、压力控制器二、温度监测装置、压差监测装置、压力监测装置三、压力监测装置四分别与所述数据采集装置相连。
16.另一方面，还提供一种地层压力下致密岩石孔渗联测方法，采用上述一种地层压力下致密岩石孔渗联测装置进行测试，包括以下步骤：
17.s1：标定所述联测装置的死体积，优选所述可变体积压力室的体积；
18.s2：将所述可变体积压力室的体积调整为优选后的体积，然后将目标岩心装入所述岩心夹持器内进行地层压力下的脉冲压力衰减测试，记录测试过程中可变体积压力室的压力随时间的变化，并绘制压力衰减曲线；
19.s3：根据所述死体积以及所述压力衰减曲线，确定对应地层压力下所述目标岩心的孔隙度和渗透率。
20.作为优选，步骤s1具体包括以下子步骤：
21.s11：将已知孔隙体积的标准岩心装入所述岩心夹持器内，确定所述可变体积压力室的当前体积；
22.s12：打开围压泵、阀门三，通过压力控制器二将所述岩心夹持器的围压升至预设压力一；
23.s13：打开阀门一、阀门二、三通阀，开启真空泵，对整个所述联测装置进行抽真空；
24.s14：关闭阀门二，打开高压气源、阀门一，通过压力控制器一将所述可变体积压力室的压力升至预设压力二，关闭高压气源、阀门一，待所述压力监测装置一监测的压力稳定后，记录所述可变体积压力室压力；
25.s15：打开阀门二，让气体等温膨胀进入所述岩心夹持器内，待所述压力监测装置一监测的压力稳定后，记录所述可变体积压力室压力；
26.s16：打开三通阀、阀门四，排出气体后卸载围压；
27.s17：更换i块不同孔隙体积的标准岩心，重复步骤s12
‑
s16，记录岩心夹持器内装入各块不同孔隙体积的标准岩心时，所述可变体积压力室在气体等温膨胀前后的压力；
28.s18：绘制标准岩心孔隙体积与气体等温膨胀后可变体积压力室压力倒数的关系曲线，获得所述关系曲线的斜率及纵截距，根据下式计算所述可变体积压力室的优选体积以及所述联测装置的死体积：
[0029][0030]
v
d
＝
‑
(b+v
k
)
ꢀꢀꢀ
(2)
[0031]
式中：v
k
为可变体积压力室的优选体积，cm3；k为v与1/p之间关系曲线的斜率，无量纲；v为标准岩心的孔隙体积，cm3；p为气体等温膨胀后可变体积压力室压力，mpa；p
k
为气体等温膨胀前可变体积压力室压力，mpa；v
d
为联测装置的死体积，cm3；b为v与1/p之间关系曲线的纵截距，无量纲。
[0032]
作为优选，还包括步骤s19：将所述可变体积压力室的体积调节至步骤s18计算获得的优选体积，重复步骤s11
‑
s19，直至第n次调节计算获得的联测装置死体积与第n
‑
1次调节计算获得的联测装置死体积之间的差值在阈值范围内；此时，第n次调节计算获得的联测装置死体积即为所述联测装置的最终死体积；第n次调节计算获得的可变体积压力室优选体积即为所述可变体积压力室的最终优选体积。
[0033]
作为优选，所述阈值范围为|v
dn
‑
v
d(n
‑
1)
|≤2％，其中v
dn
为第n次调节计算获得的联测装置死体积，cm3；v
d(n
‑
1)
为第n
‑
1次调节计算获得的联测装置死体积，cm3。
[0034]
作为优选，步骤s17中，i为小于10的正整数。
[0035]
作为优选，所述标准岩心为钢制岩心。
[0036]
作为优选，步骤s2具体包括以下子步骤：
[0037]
s21：将所述可变体积压力室的体积调整为优选后的体积，然后将目标岩心装入所述岩心夹持器内；
[0038]
s22：打开围压泵、阀门三，通过压力控制器二将所述岩心夹持器的围压升至模拟地层条件下的上覆岩石压力；打开阀门一、阀门二、三通阀，开启真空泵，对整个所述联测装置进行抽真空；
[0039]
s23：关闭三通阀，打开高压气源、阀门一、阀门二，通过压力控制器一将所述可变体积压力室的压力升至模拟地层条件下的孔隙流体压力，关闭高压气源、阀门一和阀门二；
[0040]
s24：打开高压气源、阀门一，通过压力控制器一将所述可变体积压力室的压力升至预设压力三，关闭高压气源、阀门一，待所述压力监测装置一监测的压力稳定后，记录所述可变体积压力室压力；
[0041]
s25：打开阀门二，让气体等温膨胀进入所述岩心夹持器内，通过数据采集装置采集气体等温膨胀过程中，所述可变体积压力室的压力衰减情况，直到整个系统压力达到平衡，同时绘制压力衰减曲线。
[0042]
作为优选，绘制所述压力衰减曲线时，去除所述压力衰减曲线中的波动点，直至所述衰减曲线的对数曲线为相关程度大于等于0.99的直线。
[0043]
作为优选，步骤s3中，所述目标岩心在地层条件下的孔隙度通过下式进行计算：
[0044]
[0045][0046]
式中：φ为目标岩心的孔隙度，％；v
pi
为目标岩心的孔隙体积，cm3；v
b
为目标岩心的体积，cm3；p
ik
为气体等温膨胀前可变体积压力室压力，mpa；p
ib
为气体等温膨胀后可变体积压力室压力，mpa；p
i
为孔隙流体压力，mpa，v
k
为可变体积压力室的优选体积，cm3；v
d
为联测装置的死体积，cm3；
[0047]
所述目标岩心的渗透率通过以下子步骤获得：
[0048]
绘制压力监测装置一监测的压力随时间变化的半对数压力衰减曲线，测得所述半对数压力衰减曲线的斜率，通过下式计算所述目标岩心的渗透率：
[0049][0050]
式中：k
gi
为目标岩心的渗透率，md；α
i
为半对数压力衰减曲线的斜率，s
‑1；μ为气体粘度，mpa
·
s；l为目标岩心的长度，cm；f
z
为与气体性质相关系数，无量纲；a为目标岩心的横截面积，cm2；v
d
为联测装置死体积，cm3。
[0051]
本发明的有益效果是：
[0052]
本发明能够直接同时测定目标岩样在地层压力下的孔隙度和渗透率，能够避免分开测试或岩石存在微裂缝等不稳定因素而导致的测试条件的不一致性，并提高致密岩石孔隙度和渗透率测试效率；本发明能够准确的计算获得联测装置的死体积，优选出可变体积压力室体积，使得目标岩样的孔隙度和渗透率测试结果更加准确。
附图说明
[0053]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0054]
图1为本发明地层压力下致密岩石孔渗联测装置的结构示意图；
[0055]
图2为本发明测试方法中v与1/p之间的关系曲线示意图；
[0056]
图3为本发明测试方法中压力衰减曲线的示意图；
[0057]
图4为本发明测试方法中压力随时间变化的半对数压力衰减曲线示意图。
[0058]
图中标号：
[0059]1‑
高压气源、2
‑
阀门一、3
‑
压力控制器一、4
‑
可变体积压力室、5
‑
压力监测装置一、6
‑
阀门二、7
‑
岩心夹持器、8
‑
压差监测装置、9
‑
温度监测装置、10
‑
压力监测装置二、11
‑
压力监测装置三、12
‑
三通阀、13
‑
阀门四、14
‑
压力监测装置四、15
‑
真空泵、16
‑
围压泵、17
‑
阀门三、18
‑
压力控制器二、19
‑
恒温箱、20
‑
数据采集装置。
具体实施方式
[0060]
下面结合附图和实施例对本发明进一步说明。
[0061]
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的技术特征可
以相互结合。
[0062]
需要指出的是，除非另有指明，本申请使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
[0063]
在本发明中，在未作相反说明的情况下，使用的术语“第一”、“第二”等是用于区别类似的对象，而不是用于描述特定的顺序或先后次序。应该理解这样使用的术语；使用的术语中“上”、“下”、“左”、“右”等通常是针对附图所示的方向而言，或者是针对部件本身在竖直、垂直或重力方向上而言；同样地，为便于理解和描述，“内”、“外”等是指相对于各部件本身的轮廓的内、外。但上述方位词并不用于限制本发明。
[0064]
如图1所示，本发明提供一种地层压力下致密岩石孔渗联测装置，包括高压气源1、可变体积压力室4、岩心夹持器7、围压泵16、真空泵15、恒温箱19、数据采集装置20；
[0065]
所述高压气源1与所述可变体积压力室4的入口端相连，且相连的管路上依次设有阀门一2、压力控制器一3；
[0066]
所述可变体积压力室4与压力监测装置一5相连，所述可变体积压力室4的出口端与所述岩心夹持器7的入口端相连，且相连的管路上设有阀门二6；
[0067]
所述岩心夹持器7的侧壁与所述围压泵16相连，且相连的管路上依次设有压力监测装置二10、压力控制器二18、阀门三17；所述岩心夹持器7与温度监测装置9相连；所述岩心夹持器7的入口端与出口端之间设有压差监测装置8；所述岩心夹持器7的出口端与所述真空泵15相连，且相连的管路上依次设有压力监测装置三11和三通阀12，所述真空泵15上设有压力监测装置四14，所述三通阀12的另一个出口与放空管相连，所述放空管上设有阀门四13；
[0068]
所述可变体积压力室4与所述岩心夹持器7设置在所述恒温箱19内；所述压力控制器一3、压力监测装置一5、压力监测装置二10、压力控制器二18、温度监测装置9、压差监测装置8、压力监测装置三11、压力监测装置四14分别与所述数据采集装置20相连。
[0069]
在一个具体的实施例中，所述压力监测装置一5、压力监测装置二10、压力监测装置三11、压力监测装置四14均采用压力传感器；所述温度监测装置9采用温度传感器；所述数据采集装置20采用计算机。需要说明的是，所述压力传感器、温度传感器、计算机等均为现有技术，具体结构在此不再赘述。需要说明的是，除了所本实施例中采用的压力监测装置、温度监测装置、数据采集装置外，还可采用现有技术中其他的压力监测装置、温度监测装置、数据采集装置。需要说明的是，压力控制器、岩心夹持器等也均为现有技术，其具体结构在此不再赘述。
[0070]
在一个具体的实施例中，所述可变体积压力室4内设有可移动隔板或活塞，通过调整所述隔板或活塞的位置，从而改变所述可变体积压力室4的体积。
[0071]
在一个具体的实施例中，所述阀门二6具有控制流体单向流动的功能，且单向流动的方向为由所述可变体积压力室4向所述岩心夹持器7流动的方向，如此能够保证可变体积压力室内的气体等温膨胀进入所述岩心夹持器7，并防止岩心夹持器7中的气体倒流进入所述可变体积压力室4。可选地，所述阀门二6由普通阀门与单向阀组成，且所述普通阀门靠近所述可变体积压力室4，所述单向阀靠近所述岩心夹持器7。可选地，所述阀门二6直接采用具有单向流动功能的阀门，例如国威阀门q641f型号的气动球阀。
[0072]
在一个具体的实施例中，所述高压气源1为高压氮气气源。需要说明的是，除了本
实施例中采用的高压气源外，也可采用现有技术中其他用于测试岩样孔隙度或渗透率的高压气源。
[0073]
另一方面，本发明还提供一种地层压力下致密岩石孔渗联测方法，采用上述地层压力下致密岩石孔渗联测装置进行测试，包括以下步骤：
[0074]
s1：标定所述联测装置的死体积，优选所述可变体积压力室4的体积；具体包括以下子步骤：
[0075]
s11：将已知孔隙体积的标准岩心装入所述岩心夹持器7内，确定所述可变体积压力室4的当前体积；
[0076]
s12：打开围压泵16、阀门三17，通过压力控制器二18将所述岩心夹持器7的围压升至预设压力一(围压，例如3mpa)；
[0077]
s13：打开阀门一2、阀门二6、三通阀12，开启真空泵15，对整个所述联测装置进行抽真空；
[0078]
s14：关闭阀门二6，打开高压气源1、阀门一2，通过压力控制器一3将所述可变体积压力室4的压力升至预设压力二(脉冲压力衰减测试时的脉冲衰减压力，例如0.6mpa)，关闭高压气源1、阀门一2，待所述压力监测装置一5监测的压力稳定后，记录所述可变体积压力室4压力；
[0079]
s15：打开阀门二6，让气体等温膨胀进入所述岩心夹持器7内，待所述压力监测装置一5监测的压力稳定后，记录所述可变体积压力室4压力；
[0080]
s16：打开三通阀12、阀门四13，排出气体后卸载围压；
[0081]
s17：更换i块不同孔隙体积的标准岩心，重复步骤s12
‑
s16，记录岩心夹持器7内装入各块不同孔隙体积的标准岩心时，所述可变体积压力室4在气体等温膨胀前后的压力；
[0082]
s18：绘制标准岩心孔隙体积与气体等温膨胀后可变体积压力室4压力倒数的关系曲线，获得所述关系曲线的斜率及纵截距，根据下式计算所述可变体积压力室4的优选体积以及所述联测装置的死体积：
[0083][0084]
v
d
＝
‑
(b+v
k
)
ꢀꢀꢀ
(2)
[0085]
式中：v
k
为可变体积压力室的优选体积，cm3；k为v与1/p之间关系曲线的斜率，无量纲；v为标准岩心的孔隙体积，cm3；p为气体等温膨胀后可变体积压力室压力，mpa；p
k
为气体等温膨胀前可变体积压力室压力，mpa；v
d
为联测装置的死体积，cm3；b为v与1/p之间关系曲线的纵截距，无量纲。
[0086]
其中式(1)和式(2)根据波义耳定律化简而来，具体的，根据波义耳定律可知：
[0087]
p
k
v
k
＝pv+pv
d
+v
k
ꢀꢀꢀ
(6)
[0088]
化简可得：
[0089][0090]
绘制如图2所示的v与1/p之间关系曲线，斜率k即为p
k
v
k
，化简即可得到式(1)所示的优选体积计算公式；纵截距b即为v
d
与v
k
之和的相反数，化简即可得到式(2)所示的死体积计算公式。
[0091]
在一个具体的实施例中，还包括步骤s19：将所述可变体积压力室4的体积调节至步骤s18计算获得的优选体积，重复步骤s11
‑
s19，直至第n次调节计算获得的联测装置死体积与第n
‑
1次调节计算获得的联测装置死体积之间的差值在阈值范围内；此时，第n次调节计算获得的联测装置死体积即为所述联测装置的最终死体积；第n次调节计算获得的可变体积压力室4优选体积即为所述可变体积压力室4的最终优选体积。需要说明的是，步骤s19能够获得更加准确的联测装置死体积和可变体积压力室4的优选体积，但是该步骤并非必要步骤，通过步骤s11
‑
s18即可获得联测装置的死体积和可变体积压力室4的优选体积，具体应用时要不要进行步骤s19需要根据用户需求的计算精度而定。
[0092]
可选地，所述阈值范围为|v
dn
‑
v
d(n
‑
1)
|≤2％，其中v
dn
为第n次调节计算获得的联测装置死体积，cm3；v
d(n
‑
1)
为第n
‑
1次调节计算获得的联测装置死体积，cm3。需要说明的是，本发明的阈值范围仅为优选的阈值范围，阈值范围越小，最终结果越准确，但是根据不同的精度要求，阈值范围也可以设置在2％以上或2％以下，例如|v
dn
‑
v
d(n
‑
1)
|≤3％、|v
dn
‑
v
d(n
‑
1)
|≤4％、|v
dn
‑
v
d(n
‑
1)
|≤5％、|v
dn
‑
v
d(n
‑
1)
|≤1％等。
[0093]
在一个具体的实施例中，步骤s17中，i为小于10的正整数。需要说明的是，更换标准岩心的块数越多，测试结果越准确，在块数10以上时的测试结果变化较小，为节约成本，标准岩心的块数在10以内即可。
[0094]
在一个具体的实施例中，所述标准岩心为钢制岩心。需要说明的是，钢制岩心的孔隙度和渗透率不会应为加载、卸载应力而发生改变，用其进行测试能够使测试结果更加准确，但是，本发明也可采用已知孔隙体积的真实岩样或者其他合金制成的岩心进行测试，具体选择根据用户需求的精度等来定，钢制岩心并非本发明的唯一选择，本发明的测试方法中，重点在于已知该标准岩心的孔隙体积。
[0095]
在一个具体的实施例中，所述标准岩心为钢制岩心，所述钢制岩心的轴心挖空模拟裂缝，挖空部分为圆柱形、方体型、楔形等规则的形状，方便计算孔隙体积。
[0096]
可选地，制作多块相同的标准岩心，采用本实施例的标准岩心进行标定联测装置的死体积以及确定可变体积压力室4的优选体积时，步骤s11中，将其中一块标准岩心放入所述岩心夹持器7中，步骤s17第一次更换标准岩心时，向所述岩心夹持器7中放入两块所述标准岩心，步骤s17第二次更换标准岩心时，向所述岩心夹持器7中放入三块所述标准岩心，依次类推，进行实验。
[0097]
可选地，与上述实施例不一样的是，制作多块不同的标准岩心，各标准岩心的长度之间和上述实施例一样，即，步骤s11中，放入的标准岩心长度为l，步骤s17第一次更换的标准岩心长度为2l，第二次更换的标准岩心长度为3l，依次类推。需要说明的是，本实施例中，各标准岩心的模拟裂缝形状相同，即当模拟裂缝为圆柱形时，各标准岩心的模拟裂缝内径相同，当模拟裂缝为方体型时，各标准岩心的模拟裂缝横截面积相同，依次类推。
[0098]
可选地，与上述实施例不一样的是，步骤s17更换的各标准岩心其长度与步骤s11中放入的标准岩心的长度不成等比关系，可以等差，也可以是任意其他关系。
[0099]
在另一个具体的实施例中，与上述实施例相同，所述标准岩心为轴心挖空的钢制岩心，不同的是，该标准岩心还设有与挖空部分形状匹配的钢制内芯，该内芯与所述标准岩心间隙配合，设置多个不同长度(等比、等差、或其他任意关系)或不同外径(等比、等差、或其他任意关系)的内芯，在步骤s11和步骤s17中的标准岩心中放入不同的所述内芯。
[0100]
s2：将所述可变体积压力室4的体积调整为优选后的体积，然后将目标岩心装入所述岩心夹持器7内进行地层压力下的脉冲压力衰减测试，记录测试过程中可变体积压力室4的压力随时间的变化，并绘制压力衰减曲线；具体包括以下子步骤：
[0101]
s21：将所述可变体积压力室4的体积调整为优选后的体积，然后将目标岩心装入所述岩心夹持器7内，目标岩心装入前根据《gb/t 29172
‑
2012岩心分析方法》进行洗盐和烘干；
[0102]
s22：打开围压泵16、阀门三17，通过压力控制器二18将所述岩心夹持器7的围压升至模拟地层条件下的上覆岩石压力；打开阀门一2、阀门二6、三通阀12，开启真空泵15，对整个所述联测装置进行抽真空；
[0103]
s23：关闭三通阀12，打开高压气源1、阀门一2、阀门二6，通过压力控制器一3将所述可变体积压力室4的压力升至模拟地层条件下的孔隙流体压力，关闭高压气源1、阀门一2和阀门二6；
[0104]
s24：打开高压气源1、阀门一2，通过压力控制器一3将所述可变体积压力室4的压力升至预设压力三(预设压力三为孔隙流体压力与脉冲衰减压力之和，脉冲衰减压力与预设压力二保持一致，脉冲衰减压力保持一致能够保证步骤s1计算得到的死体积与优选出来的可变体积压力室体积可用)，关闭高压气源1、阀门一2，待所述压力监测装置一5监测的压力稳定后，记录所述可变体积压力室4压力；
[0105]
s25：打开阀门二6，让气体等温膨胀进入所述岩心夹持器7内，通过数据采集装置20采集气体等温膨胀过程中，所述可变体积压力室4的压力衰减情况，直到整个系统压力达到平衡，同时绘制如图3所示的压力衰减曲线。
[0106]
在一个具体的实施例中，绘制所述压力衰减曲线时，去除所述压力衰减曲线中的波动点，直至所述衰减曲线的对数曲线为相关程度大于等于0.99的直线。需要说明的是，本发明并不限制该相关程度一定大于等0.99，相关程度越接近1，后续计算结果的准确性越高，0.99仅为优选的相关程度，小于0.99相关程度的直线例如0.98、0.95等时也行，这仅能说明波动点去除的不够，而波动点的去除是为了使后续结果更加准确，因此可根据计算精度要求去除波动点。
[0107]
s3：根据所述死体积以及所述压力衰减曲线，确定对应地层压力下所述目标岩心的孔隙度和渗透率。所述目标岩心的孔隙度通过下式进行计算：
[0108][0109][0110]
式中：φ为目标岩心的孔隙度，％；v
pi
为目标岩心的孔隙体积，cm3；v
b
为目标岩心的体积，cm3；p
ik
为气体等温膨胀前可变体积压力室压力，mpa；p
ib
为气体等温膨胀后可变体积压力室压力，mpa；p
i
为孔隙流体压力，mpa，v
k
为可变体积压力室的优选体积，cm3；v
d
为联测装置的死体积，cm3；
[0111]
所述目标岩心的渗透率通过以下子步骤获得：
[0112]
绘制如图4所示的压力监测装置一5监测的压力随时间变化的半对数压力衰减曲线，测得所述半对数压力衰减曲线的斜率，通过下式计算所述目标岩心的渗透率：
[0113][0114]
式中：k
gi
为目标岩心的渗透率，md；α
i
为半对数压力衰减曲线的斜率，s
‑1；μ为气体粘度，mpa
·
s；l为目标岩心的长度，cm；f
z
为与气体性质相关系数，无量纲；a为目标岩心的横截面积，cm2；v
d
为联测装置死体积(即岩心夹持器出口端至三通阀入口端之间的管线体积)，cm3。
[0115]
在一个具体的实施例中，若已知某储层岩石在地层压力下的的孔隙流体压力随时间下降的曲线，则在采用本发明所述的地层压力下致密岩石孔渗联测方法进行孔渗联测时，步骤s2进行脉冲压力衰减测试过程中，步骤s23将可变体积压力室的压力升至模拟地层条件下的孔隙流体压力时，该压力为孔隙流体压力随时间下降的曲线上的p1，此时步骤s25绘制的压力衰减曲线为p1孔隙流体压力下的压力衰减曲线。步骤s25之后还包括步骤s26：打开三通阀12和阀门13，使可变体积压力室的压力下降至p2，p2为孔隙流体压力随时间下降的曲线上的另一个孔隙流体压力，且p2小于p1，重复步骤s24
‑
s25，获得p2孔隙流体压力下的压力衰减曲线，重复步骤s26，每次重复时将可变体积压力室4的压力依次下降至孔隙流体压力随时间下降的曲线上的p3、p4……
p
i
，从而获得各孔隙流体压力下的压力衰减曲线。后续步骤s3可计算得到p1、p2……
p
i
等各孔隙流体压力下对应的孔隙度和渗透率，从而获得该地层孔隙度和渗透率随时间变化的曲线。
[0116]
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。