一种基于毛管压力的岩石电阻率频谱测试方法、联测装置

1.本文涉及油气田开发领域，尤其是一种基于毛管压力的岩石电阻率频谱测试方法、联测装置。


背景技术：

2.在油气储层勘探开发过程中，地球物理测井是评价储层的一种重要手段，电法测井技术是储层饱和度定量评价的主要依据。
3.现有技术中通常利用普通电阻率测井测量地层电阻率，结合archie公式评价储层饱和度信息。其中，普通电阻率测井可以看作单一频率点测量，其无法对低孔、低渗、低阻储层以及水淹层等杂储层进行准确评价。渗透率是岩石的另一项重要物性参数。对于常规储层，孔隙度与渗透率之间通常存在正相关关系，即孔隙度越大，岩石的渗透率也就越大。但对于探勘开发中的杂储层，对应的孔隙度与渗透率则相关性较弱，所以需要通过其他的测量方法来评价地层的孔隙结构。
4.针对目前技术存在的普通电阻率无法对杂储层进行准确评价的问题，需要一种基于毛管压力的岩石电阻率频谱测试方法、联测装置。


技术实现要素：

5.为解决上述现有技术的问题，本文实施例提供了一种基于毛管压力的岩石电阻率频谱测试方法、联测装置。
6.本文实施例提供了一种基于毛管压力的岩石电阻率频谱测试方法，所述方法包括：获取岩心的饱含水电阻率频谱数据；将所述饱含水电阻率频谱数据输入至预先构建的毛管压力类型预测模型中，确定所述饱含水电阻率频谱数据对应的毛管压力曲线类型，所述毛管压力曲线类型由毛管压力曲线的特征点数量及所述特征点对应的压差值表示；根据所述特征点确定压差测量序列，包括最小压差值、最大压差值、压差门槛值；测量压差测量序列中各点的电阻率频谱数据。
7.根据本文实施例的一个方面，所述毛管压力类型预测模型的训练过程如下：获取不同岩心的历史饱含水电阻率频谱数据及所述历史饱含水电阻率频谱样本数据对应的历史毛管压力曲线数据；从所述历史毛管压力曲线数据中选择随机数量的特征点作为初始特征点，所述初始特征点不包括历史毛管压力曲线的首尾端点，所述初始特征点包括其对应的含水饱和度及压差值；将所述初始特征点、所述历史饱含水电阻率频谱数据确定为训练样本集；根据所述训练样本集，构建初始毛管压力曲线；计算所述初始毛管压力曲线与历史毛管压力曲线数据的误差；根据所述误差及最大迭代次数，迭代调整所述初始毛管压力曲线的初始特征点的数量，当所述误差满足预设误差值，确定毛管压力类型预测模型。
8.根据本文实施例的一个方面，所述根据误差及所述最大迭代次数，迭代调整初始毛管压力曲线的特征点数量，确定毛管压力类型预测模型包括：判断当前迭代次数是否超过最大迭代次数；若否，进一步判断所述误差是否大于预设误差值；若所述误差小于等于预
设误差值，确定初始毛管压力曲线的最小特征点数，确定所述最小特征点数对应的毛管压力类型预测模型；若所述误差大于所述预设误差值，调整所述初始特征点的数量，重新执行迭代所述初始毛管压力曲线及判断误差与预设误差值的步骤，直到当前迭代次数超过最大迭代次数；若是，调整所述初始特征点的数量，重新执行迭代所述初始毛管压力曲线及判断误差与预设误差值的步骤。
9.根据本文实施例的一个方面，所述调整所述初始特征点的数量，重新执行判断误差与预设误差值的步骤包括：依次增加初始特征点的数量，计算所述初始毛管压力曲线与历史毛管压力曲线数据的误差；或依次减少初始特征点的数量，计算所述初始毛管压力曲线与历史毛管压力曲线数据的误差。
10.根据本文实施例的一个方面，根据所述特征点及所述压差值，确定压差测量序列；当压差测量序列中的各压差测量点满足压差平衡条件，测量所述压差平衡条件下的电阻率频谱，所述压差平衡条件包括压差平衡时间、平衡出液量。
11.根据本文实施例的一个方面，所述方法进一步包括：对历史饱含水电阻率频谱数据进行预处理，包括：滤波降噪、曲线平滑、标准化处理。
12.本文实施例提供了一种岩石电阻率频谱与毛管压力联测装置，所述装置包括岩心夹持器、电频谱测量系统、液压系统，所述岩心夹持器包括上下夹持活塞、夹持岩心室、上下出液管、上下泄压阀、半渗透隔板，上下活塞连接两供电电极，岩心室侧端连接两测量电极，用于夹持岩心；所述电频谱测量系统包括阻抗分析仪、供电电极和测量电极，用于测量电阻率频谱数据；所述液压系统包括上游压力传感器、下游压力传感器、压差传感器、计量泵、驱替泵、围压泵及泵控制器，用于对岩心施加不同大小的压力差。
13.根据本文实施例的一个方面，电频谱测量系统进一步用于，根据毛管压力类型预测模型预测输出的特征点数量、所述特征点对应的压差值，测量各特征点的电阻率频谱数据。
14.本文实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述基于毛管压力的岩石电阻率频谱测试方法。
15.本文实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现所述基于毛管压力的岩石电阻率频谱测试方法。
16.本方案通过预测毛管压力曲线类型，减少测量压差点数，提升测量效率；本文通过电阻率与毛管压力联测，实现在不同的压力测量点测量不同的电阻率频谱曲线，比单一频谱的电阻率提供更多地层信息，有利于未来地质研究的发展。
附图说明
17.为了更清楚地说明本文实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本文的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
18.图1所示为本文实施例一种基于毛管压力的岩石电阻率频谱测试方法的流程图；
19.图2所示为本文实施例一种确定毛管压力类型预测模型的方法流程图；
20.图3所示为本文实施例一种确定毛管压力类型预测模型的方法流程图；
21.图4所示为本文实施例一种测量压差测量序列中电阻率频谱数据的方法流程图；
22.图5所示为本文实施例一种岩石电阻率频谱与毛管压力联测装置的结构示意图；
23.图6所示为本文实施例一种岩心夹持器的剖面示意图；
24.图7a所示为本文实施例一种孔隙结构的毛管压力曲线示意图；
25.图7b所示为本文实施例又一种孔隙结构的毛管压力曲线示意图；
26.图8所示为本文实施例一种计算机设备的结构示意图。附图符号说明：
27.501、驱替泵；
28.502、围压泵；
29.503、计量泵；
30.504、驱替泵吸液瓶；
31.505、围压泵吸液瓶；
32.506、计量泵吸液瓶；
33.507、岩心室上端压力传感器；
34.508、岩心室下端压力传感器；
35.509、压差传感器；
36.510、管道控制开关；
37.511、温度传感器；
38.512、岩心室下端泄压阀；
39.513、驱替泵储油罐；
40.514、驱替泵储水罐；
41.515、数据采集板；
42.516、计算机显示屏；
43.517、计算机；
44.518、阻抗分析仪；
45.519、恒温箱；
46.520、压力泵手动控制仪；
47.521、岩心室上端泄压阀；
48.522、岩心室上端温度传感器；
49.523、岩心夹持器；
50.601、岩心室上端夹持旋钮；
51.602、上端供电电极；
52.603、驱替端口；
53.604、上端排气口；
54.605、岩心样品；
55.606、围压泄压阀；
56.607、半渗透隔板；
57.608、下端排气口；
58.609、计量端口；
59.610、下端供电电极；
60.611、岩心室下端夹持旋钮；
61.612、橡胶套；
62.613、第一测量电极；
63.614、第二测量电极；
64.615、围压输入端口；
65.616、夹持器金属外壳；
66.802、计算机设备；
67.804、处理器；
68.806、存储器；
69.808、驱动机构；
70.810、输入/输出模块；
71.812、输入设备；
72.814、输出设备；
73.816、呈现设备；
74.818、图形用户接口；
75.820、网络接口；
76.822、通信链路；
77.824、通信总线。
具体实施方式
78.为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本文实施例中的附图，对本文实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本文一部分实施例，而不是全部的实施例。基于本文中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本文保护的范围。
79.需要说明的是，本文的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本文的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、装置、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
80.本说明书提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的系统或装置产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行。
81.需要说明的是，本文的方法可用于交通运输领域，也可用于除交通运输领域之外的任意领域，本文方法及装置的应用领域不做限定。
82.图1所示为本文实施例一种基于毛管压力的岩石电阻率频谱测试的方法流程图，具体包括如下步骤：
83.步骤101，获取岩心的饱含水电阻率频谱数据。在本说明书的一些实施例中，饱含水电阻率频谱数据为含水饱和度为100％的电阻率频谱数据。当岩心所受到的外界施加的压力为0，岩心内部孔隙没有驱替出液体时，测量岩心当前状态下的电阻率频谱数据。饱含水电阻率频谱数据可以准确表示毛管压力曲线的曲线走向。进一步的，具有不同孔隙结构的岩心或岩样的饱含水电阻率频谱数据是不同的。本说明书中，饱含水电阻率频谱数据可以以连续曲线的形式表示，电阻率值随着频率的变化而发生变化，电阻率的变化从微观上反映地层孔隙结构。具体的，随着频率增加，饱含水电阻率频谱的电阻率逐渐下降。饱含水电阻率频谱的变化信息与岩心的孔隙结构具有较大的相关性，可以从微观反映孔隙结构的具体形态。因此，饱含水电阻率频谱数据一定程度上可以反映岩心/岩样的孔隙结构。
84.步骤102，将所述饱含水电阻率频谱数据输入至预先构建的毛管压力类型预测模型中，确定所述饱含水电阻率频谱数据对应的毛管压力曲线类型，所述毛管压力曲线类型由毛管压力曲线的特征点数量及所述特征点对应的压差值表示。例如，如图7a所示，某毛管压力曲线的特征点数量为7点，这7个特征点对应的含水饱和度分别为：0.45、0.5、0.6、0.8、0.92、0.99、1；这7个特征点对应的压差值分别为120pc、90pc、20pc、10pc、9pc、7pc、0pc，该毛管压力曲线表示某一孔隙结构的岩样；又例如，如图7b所示，某毛管压力曲线的特征点数量为5点，这5个特征点对应的含水饱和度分别为：0.85、0.91、0.94、0.99、1；这5个特征点对应的压差值分别为250pc、90pc、60pc、35pc、0pc，该毛管压力曲线表示具有另一孔隙结构的岩样。
85.在本步骤中，预先构建的毛管压力类型预测模型是由历史饱含水电阻率频谱数据、与历史饱含水电阻率频谱数据对应的历史毛管压力曲线数据训练得到的。关于毛管压力类型预测模型的训练过程详见图2的具体描述。
86.在本说明书中，将步骤101中获取的饱含水电阻率频谱数据作为毛管压力类型预测模型的输入，毛管压力类型预测模型输出与饱含水电阻率频谱数据对应的岩心所对应的毛管压力曲线类型。其中，毛管压力曲线类型通过毛管压力曲线的特征点数量及这些特征点的压差值表征。例如，某岩样的孔隙结构较为疏松，通过该岩样的饱含水电阻率频谱数据预测得到的毛管压力类型可以仅由5个特征点表示；而某岩样的孔隙结构较为紧密，通过该岩样的饱含水电阻率频谱数据预测得到的毛管压力类型需要由7个特征点表示。进一步的，预测得到的毛管压力曲线中的每个特征点具有对应的压差值及含水包含度。
87.步骤103，根据所述特征点确定压差测量序列，包括最小压差值、最大压差值、压差门槛值。本步骤中，压差测量序列包括多个压差测量点，每个压差测量点为步骤102中确定的特征点，每一特征点对应有压差值信息。压差测量序列中，至少包括最大压差值、最小压差值、压差门槛值。进一步的，压差测量序列还可以包括其他压差测量值。
88.步骤104，测量压差测量序列中各点的电阻率频谱数据。本步骤结合图5所述的岩石电阻率频谱与毛管压力联测装置进行电阻率频谱数据的测量。
89.图2所示为本文实施例一种确定毛管压力类型预测模型的方法流程图。本说明书通过训练得到可以预测毛管压力曲线类型的模型，利用最少的特征点表示毛管压力曲线，进而减少电阻率频谱的测量时间。具体包括如下步骤：
90.步骤201，获取不同岩心的历史饱含水电阻率频谱数据及所述历史饱含水电阻率频谱样本数据对应的历史毛管压力曲线数据。本步骤中，不同岩心的历史饱含水电阻率频谱数据及历史饱含水电阻率频谱数据可以是历史采集或测量得到的数据。
91.步骤202，从所述历史毛管压力曲线数据中选择随机数量的特征点作为初始特征点，所述初始特征点不包括历史毛管压力曲线的首尾端点，所述初始特征点包括其对应的含水饱和度及压差值。
92.在本说明书的一些实施例中，从历史毛管压力曲线数据中选择随机数量的特征点，作为初始特征点，同时也作为毛管压力类型预测模型的训练样本的标签。其中，可以从历史毛管压力曲线数据中随机选择3个、4个、5个、6个、7个、8个甚至更多数量的特征点作为初始特征点。其中，随机选择的初始特征点不为历史毛管压力曲线的曲线起始点。在本步骤中，随机数量的特征点(即，初始特征点)为少数特征点。
93.步骤203，将所述初始特征点、所述历史饱含水电阻率频谱数据确定为训练样本集。其中，历史饱含水电阻率频谱数据为训练样本集中的特征数据，初始特征点为训练样本数据集中的标签。其中，训练样本集可以与一种孔隙度结构的岩样对应。例如，具有紧密孔隙的岩样a1对应的历史饱含水电阻率频谱数据、历史毛管压力数据中的5个初始特征点，作为岩样a1对应的训练样本集；具有中等紧密程度的岩样a2的历史饱含水电阻率频谱数据、历史毛管压力数据中的7个初始特征点，作为岩样a2的训练样本集；具有疏松孔隙的岩样a3对应的历史饱含水电阻率频谱数据、历史毛管压力数据中的9个初始特征点，作为岩样a3的训练样本集。本说明书中，对具有不同孔隙结构的岩样，利用该岩样的初始特征点作为训练样本集，训练模型。
94.步骤204，根据所述训练样本集，构建初始毛管压力曲线。利用历史包饱含水电阻率频谱数据及随机数量的训练样本集拟合初始毛管压力曲线。在本步骤中，利用二项式拟合以构建毛管压力曲线。
95.步骤205，计算所述初始毛管压力曲线与历史毛管压力曲线数据的误差。在本步骤中，二次多项式拟合计算拟合的初始毛管压力曲线与原本的历史毛管压力曲线的误差。其中，误差可以是平均相对误差，为构建的初始毛管压力曲线中各个初始特征点与历史毛管压力曲线中对应的各点之间的平均相对误差累加和，即为，累积平均相对误差。
96.步骤206，根据所述误差及最大迭代次数，迭代调整所述初始毛管压力曲线的初始特征点的数量，当误差满足预设误差值，确定毛管压力类型预测模型。在本步骤中，最大迭代次数为预先设定的初始毛管压力曲线迭代次数上限，例如，设定最大迭代次数为3次、5次、6次、10次等。当拟合的初始毛管压力曲线中各初始特征点与历史毛管压力曲线的累积平均相对误差满足误差标准，即，各初始特征点的累积平均相对误差小于等于误差标准，则可以确定当前拟合得到的初始毛管压力曲线较为接近历史毛管压力曲线。
97.图3所示为本文实施例一种确定毛管压力类型预测模型的方法流程图，具体包括如下步骤：
98.步骤301，判断当前迭代次数是否超过最大迭代次数。如步骤206所述，最大迭代次数为预先设定的初始毛管压力曲线迭代次数上限，用于控制模型训练的程度。
99.步骤302，若否，进一步判断所述误差是否大于预设误差值。
100.步骤303，若所述误差小于等于预设误差值，确定初始毛管压力曲线的最小特征点
数，确定所述最小特征点数对应的毛管压力类型预测模型。在本步骤中，对同一孔隙类型的岩样，可以使用不同数量的初始特征点构建初始毛管压力曲线。例如，对同一孔隙类型的岩样，使用3个初始特征点和使用5个、7个初始特征点构建的初始毛管压力曲线存在不同。则在利用不同数量的初始特征点计算与历史毛管压力曲线的误差的时候，累加误差值也不同。例如，同一孔隙类型的岩样，利用3个初始特征点计算的误差为3％；利用5个初始特征点计算的与历史毛管压力曲线的误差为4％；利用7个初始特征点计算的与历史毛管压力曲线的误差为5％。在本步骤中，当同一孔隙类型的岩样的不同数量的初始特征点与历史毛管压力曲线的误差均满足小于等于预设误差值的条件时，选择最小数量的初始特征点数，将最小数量的初始特征点数对应的毛管压力曲线确定为毛管压力类型预测模型。例如，预设误差值为5％，而同一孔隙类型的岩样，取3个、5个、7个初始特征点构建的初始毛管压力曲线与历史毛管压力曲线的累积相对评价误差均小于或等于5％，则选择3个初始特征点对应的初始毛管压力曲线作为毛管压力类型预测模型。
101.步骤304，若所述误差大于所述预设误差值，调整所述初始特征点的数量，重新执行迭代所述初始毛管压力曲线及判断误差与预设误差值的步骤，直到当前迭代次数超过最大迭代次数。当误差大于所述预设误差值，认为该孔隙类型的岩样所取的初始特征点的数量不符合曲线特征，需要调整初始特征点的数量，重新构建初始毛管压力曲线。具体的，可以通过依次增加初始特征点的数量或依次减少初始特征点的数量以构建初始毛管压力曲线。直到构建出的新的初始毛管压力曲线与原始的历史毛管压力曲线之间满足误差条件。则认为构建好的毛管压力曲线的特征点数量为所取点的数量。
102.步骤305，若是，调整所述初始特征点的数量，重新执行迭代所述初始毛管压力曲线及判断误差与预设误差值的步骤。若当前迭代次数超过最大迭代次数，说明当前迭代初始毛管压力曲线的过程已不符合训练标准或训练要求，需要通过调整初始特征点的数量，重新迭代毛管压力曲线。
103.在本说明书的一些实施例中，调整所述初始特征点的数量，重新执行判断误差与预设误差值的步骤包括：依次增加初始特征点的数量，计算所述初始毛管压力曲线与历史毛管压力曲线数据的误差；或依次减少初始特征点的数量，计算所述初始毛管压力曲线与历史毛管压力曲线数据的误差。
104.在本说明书的一些实施例中，还包括对历史饱含水电阻率频谱数据进行预处理，包括：滤波降噪、曲线平滑、标准化处理。
105.图4所示为本文实施例一种测量压差测量序列中电阻率频谱数据的方法流程图，具体包括如下步骤：
106.步骤401，根据所述特征点及所述压差值，确定压差测量序列。根据特征点对应的含水饱和度、压差值，可以确定压差测量序列，压差测量序列包括两个或两个以上的特征点。
107.步骤402，当压差测量序列中的各压差测量点满足压差平衡条件，测量所述压差平衡条件下的电阻率频谱，所述压差平衡条件包括压差平衡时间、平衡出液量。
108.在本说明书的一些实施例中，在测量电阻率频谱时，需要确保每一压差测量点达到的平衡条件。其中，平衡条件要求平衡检测时间和监测时间的出液量之差小于平衡出液量。具体的，预先通过计算机软件设定压差平衡监测时间、平衡出液量以及不同的压差测量
点。当每一压差条件下满足平衡监测时间以及监测时间的出液量之差小于平衡出液量，则认为已经达到该压差的稳定点。此时记录总出液量及压差值，并测量当前条件下的岩样的复电阻率频谱，则测量结束。接着测量下一个设定的压差测量点并重复上述操作，直到所有设定压差测量点均完成实验测量。保存实验数据并生成关于毛管压力与复电阻率频谱对应关系的实验测量报告。
109.图5所示为本文实施例一种岩石电阻率频谱与毛管压力联测装置的结构示意图，所述装置包括岩心夹持器、电频谱测量系统、液压系统。
110.具体的，该装置中包括为驱替泵501，围压泵502，计量泵503，驱替泵吸液瓶504，围压泵吸液瓶505，计量泵吸液瓶506，岩心室上端压力传感器507，岩心室下端压力传感器508，压差传感器509，管道控制开关510，温度传感器511，岩心室下端泄压阀512，驱替泵储油罐513，驱替泵储水罐514，数据采集板515，计算机显示屏516，计算机517，阻抗分析仪518，恒温箱519，压力泵手动控制仪520，岩心室上端泄压阀521，岩心室上端温度传感器522，岩心夹持器523。其中，阻抗分析仪518用于测量岩石样品的复电阻率频谱。测量频率范围为40hz-110mhz，支持扫频测量。压力泵手动控制仪520用于控制计量泵、围压泵以及驱替泵产生压力，压力的数值范围为0-10000psi。
111.其中，岩心夹持器523包括上下夹持活塞、夹持岩心室、上下出液管、上下泄压阀、半渗透隔板，上下活塞连接两供电电极，夹持岩心室侧端连接两测量电极，岩心夹持器523用于夹持岩心。在本说明书的一些实施例中，通过在岩心夹持器内加围压及产生上下端压力差的方式，使半渗透隔板驱水并计算出驱替出的水量，采用四极法测量在该压差条件下的电阻率频谱。
112.在本说明书的一些实施例中，岩石电阻率频谱与毛管压力联测装置还包括恒温控制系统，恒温控制系统包括恒温箱519、温度传感器、温度采集板、测温电路以及电加热器。恒温系统测量温度数据通过导线连接至计算机517。
113.在本说明书的一些实施例中，电频谱测量系统包括阻抗分析仪518、供电电极和测量电极，用于测量电阻频谱数据。电极通过导线与岩心夹持器连接，阻抗分析仪518通过导线与计算机517连接。
114.在本说明书的一些实施例中，液压系统包括上游压力传感器、下游压力传感器、压差传感器、计量泵、驱替泵、围压泵以及泵控制器，液压系统用于对岩心施加不同大小的压力差。计量泵、驱替泵及围压泵均可由泵控制器手动操作，泵控制器通过导线与计算机517连接。其中，上游压力传感器、下游压力传感器、压差传感器、温度测量电路通过导线与多路采集控制板连接，多路采集控制板安装在机柜内并通过导线与计算机517连接，各个系统均可采用计算机517控制，所测量得到的数据包括阻抗分析仪的复电阻率频谱数据，毛管压力测量部分的计量泵、驱替泵、围压泵溶液容积记录数据，岩心上下端压力数据以及温度控制系统的温度记录数据，这些数据进一步通过计算机517进行进一步处理。
115.图6所示为本文实施例一种岩心夹持器的剖面示意图。其中601为岩心室上端夹持旋钮，602为上端供电电极，603为驱替端口，604为上端排气口，605为岩心样品，606为围压泄压阀，607为半渗透隔板，608为下端排气口，609为计量端口，610为下端供电电极，611为岩心室下端夹持旋钮，612为橡胶套，613为第一测量电极，614为第二测量电极，615为围压输入端口，616为夹持器金属外壳。
116.在本说明书的一些实施例中，在进行电阻率频谱数据测量实验之前，需要对计量泵、驱替泵、围压泵进行清洗，尤其是计量泵，以将计量泵中的液体排出。还需要对管线进行清洗，并将管线中的存在的空气排出。并用蒸馏水将泵的液缸清洗两遍，用实验所用盐溶液再清洗两遍，然后再吸入实验用盐溶液至泵的液缸内。另外，驱替泵和储液罐内也要确保排空空气。
117.在本说明书的一些实施例中，在装载半渗透隔板与岩心样品的步骤中，半渗透隔板在使用前需与待测岩心样品一起经过加压饱和处理。在装载过程中，先装载半渗透隔板到岩心夹持器内，确保半渗透隔板与夹持器两端电极接触良好，再放置岩心确保样品端面与半渗透隔板和电极紧密接触，拧紧夹持器两端电极。
118.在本说明书的一些实施例中，管线中的气体已在前述步骤排出，但夹持器两端仍存在少量气体，在实验前需排出。具体的，需先加围压，岩心样品与夹持器内壁间可能存在间隙，如果不加围压直接进行夹持器两端端口的排气操作可能会使两端液体沿夹持器内壁流过。待围压使得橡胶套变形后，进行排气操作。排气时打开排气口阀门，设定一定的压力，排气口会现有气泡冒出，然后有液体排出，待液体稳定排出没有气泡时关闭排气口阀门。
119.图7a所示为本文实施例一种孔隙结构的毛管压力曲线示意图，图中毛管压力曲线的横坐标为含水饱和度sw，纵坐标为压力差值(也为，压差值)pc/psi。图中毛管压力曲线中的特征点数量为7，每个特征点都能反映该点的含水饱和度及压差值；图7b所示为本文实施例又一种孔隙结构的毛管压力曲线示意图，图中毛管压力曲线的横坐标为sw，纵坐标为压力差值(也为，压差值)。图中毛管压力曲线中的特征点数量为5，都能反映该点的含水饱和度及压差值。
120.如图8所示，为本文实施例提供的一种计算机设备，所述计算机设备802可以包括一个或多个处理器804，诸如一个或多个中央处理单元(cpu)，每个处理单元可以实现一个或多个硬件线程。计算机设备802还可以包括任何存储器806，其用于存储诸如代码、设置、数据等之类的任何种类的信息。非限制性的，比如，存储器806可以包括以下任一项或多种组合：任何类型的ram，任何类型的rom，闪存设备，硬盘，光盘等。更一般地，任何存储器都可以使用任何技术来存储信息。进一步地，任何存储器可以提供信息的易失性或非易失性保留。进一步地，任何存储器可以表示计算机设备802的固定或可移除部件。在一种情况下，当处理器804执行被存储在任何存储器或存储器的组合中的相关联的指令时，计算机设备802可以执行相关联指令的任一操作。计算机设备802还包括用于与任何存储器交互的一个或多个驱动机构808，诸如硬盘驱动机构、光盘驱动机构等。
121.计算机设备802还可以包括输入/输出模块810(i/o)，其用于接收各种输入(经由输入设备812)和用于提供各种输出(经由输出设备814)。一个具体输出机构可以包括呈现设备816和相关联的图形用户接口(gui)818。在其他实施例中，还可以不包括输入/输出模块810(i/o)、输入设备812以及输出设备814，仅作为网络中的一台计算机设备。计算机设备802还可以包括一个或多个网络接口820，其用于经由一个或多个通信链路822与其他设备交换数据。一个或多个通信总线824将上文所描述的部件耦合在一起。
122.通信链路822可以以任何方式实现，例如，通过局域网、广域网(例如，因特网)、点对点连接等、或其任何组合。通信链路822可以包括由任何协议或协议组合支配的硬连线链路、无线链路、路由器、网关功能、名称服务器等的任何组合。
123.对应于图1至图4中的方法，本文实施例还提供了一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时执行上述方法的步骤。
124.本文实施例还提供一种计算机可读指令，其中当处理器执行所述指令时，其中的程序使得处理器执行如图1至图4所示的方法。
125.应理解，在本文的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本文实施例的实施过程构成任何限定。
126.还应理解，在本文实施例中，术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系。例如，a和/或b，可以表示：单独存在a，同时存在a和b，单独存在b这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。
127.本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本文的范围。
128.所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
129.在本文所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口、装置或单元的间接耦合或通信连接，也可以是电的，机械的或其它的形式连接。
130.所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本文实施例方案的目的。
131.另外，在本文各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
132.所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本文的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本文各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的
介质。
133.本文中应用了具体实施例对本文的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本文的方法及其核心思想；同时，对于本领域的一般技术人员，依据本文的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本文的限制。