CN115506758B - 一种致密储层排驱压力确定方法、装置、设备和存储介质 - Google Patents

Abstract

本文提供一种致密储层排驱压力确定方法、装置、设备和存储介质，其中方法包括：获取岩心样本的压汞数据，所述压汞数据包括多个进汞压力和各个进汞压力对应的进汞饱和度；根据所述进汞压力和所述进汞饱和度，计算相邻两个进汞压力对应的单位增压饱和度增量；选择满足筛选条件的单位增压饱和度增量对应的进汞压力作为排驱压力值。本文提供的致密储层排驱压力确定方法从排驱压力的定义出发，根据单位增压饱和度增量来确定排驱压力值，结果更加准确，能够为储层研究和开发提供可靠的理论支撑；且求取过程简单便捷，提高了对致密储层排驱压力确定的效率。

Description

一种致密储层排驱压力确定方法、装置、设备和存储介质

技术领域

本发明涉及石油勘探开发技术领域，尤其是一种致密储层排驱压力确定方法、装置、设备和存储介质。

背景技术

排驱压力，也称为阈压或门槛压力，它是非润湿相开始连续进入岩样最大喉道时所对应的毛管压力。排驱压力是研究储层储集性能和渗流性能的一个重要指标，排驱压力的确定对储层评价具有重要的意义。

排驱压力的确定方法可分为毛管压力曲线法、实验室直接测量法、地球物理评价法和影响因素综合分析法这四类。其中，毛管压力曲线切线法是最为常用的、且能快速确定排驱压力的方法：在半对数坐标中沿着毛管压力曲线平坦部分的第一个拐点做切线，切线延长线与纵坐标轴相交的压力点即为排驱压力。

然而，在致密油储层岩心压汞资料中，毛管压力曲线在开始进汞阶段大多没有明显的平坦部分，以曲线开始稳定时对应的毛细管压力来确定排驱压力的方法难以适用，很难精确地做出它的切线。因此，对毛管压力曲线做切线的方法具有较大的局限性。当拐点难以确定时，还可以以相邻两个测点进汞饱和度差大于等于2％回退一个测压点的压力作为排驱压力。但这种方法获得的排驱压力确定值通常较高于实际值，具有较大误差；加之压汞发具有的滞后效应，使得获得的数据不能准确反映实际排驱压力特征。

有鉴于此，本文旨在提供一种致密储层排驱压力确定方法、装置、设备和存储介质，以解决现有技术中对致密储层排驱压力获取不准确的问题。

发明内容

针对现有技术的上述问题，本文的目的在于，提供一种致密储层排驱压力确定方法、装置、设备和存储介质，以解决现有技术中获取排驱压力不便捷和不准确的问题。

为了解决上述技术问题，本文的具体技术方案如下：

一方面，本文提供一种致密储层排驱压力确定方法，包括：

获取岩心样本的压汞数据，所述压汞数据包括多个进汞压力和各个进汞压力对应的进汞饱和度；

根据所述进汞压力和所述进汞饱和度，计算相邻两个进汞压力对应的单位增压饱和度增量；

选择满足筛选条件的单位增压饱和度增量对应的进汞压力作为排驱压力值。

具体地，所述选择满足筛选条件的单位增压饱和度增量对应的进汞压力作为排驱压力值，进一步为：

获取岩心样本的参数数据，所述参数数据包括所述岩心样本的孔喉分布特征和渗透率；

根据所述参数数据，确定排驱压力范围；

选择满足所述排驱压力范围的单位增压饱和度增量对应的进汞压力作为排驱压力值。

进一步地，所述选择满足筛选条件的单位增压饱和度增量对应的进汞压力作为排驱压力值，进一步为：

根据所述进汞压力值和所述进汞饱和度，绘制毛管压力曲线；

根据所述毛管压力曲线的形态和值域，将所述毛管压力曲线划出分麻皮效应区和喉道初始进汞区；

基于划分出的所述麻皮效应区和所述喉道初始进汞区，选择满足所述排驱压力范围的单位增压饱和度增量对应的进汞压力作为排驱压力值。

进一步地，计算相邻两个进汞压力对应的单位增压饱和度增量的公式为：

SHg′|Pc＝△SHg/△Pc

其中，SHg′|Pc为单位增压饱和度增量，△SHg为饱和度增量，△Pc为进汞压力增量。

进一步地，所述基于划分出的所述麻皮效应区和所述喉道初始进汞区，选择满足所述排驱压力范围的单位增压饱和度增量对应的进汞压力作为排驱压力值为：

基于划分出的所述麻皮效应区和所述喉道初始进汞区，选择满足所述排驱压力范围且单位增压饱和度增量大于等于1％/MPa的进汞压力作为排驱压力值。

优选地，所述方法还包括：

对所述排驱压力值进行修正。

优选地，所述根据所述进汞压力和所述进汞饱和度，计算相邻两个进汞压力对应的单位增压饱和度增量，还包括：

增加进汞压力和进汞饱和度的数量。

另一方面，本说明书实施例还提供一种致密储层排驱压力确定装置，包括：

获取模块，用于获取岩心样本的压汞数据，所述压汞数据包括多个进汞压力和各个进汞压力对应的进汞饱和度；

计算模块，根据所述进汞压力和所述进汞饱和度，计算相邻两个进汞压力对应的单位增压饱和度增量；

选择模块，用于选择满足筛选条件的单位增压饱和度增量对应的进汞压力作为排驱压力值。

另一方面，本说明书实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述技术方案提供的方法步骤。

另一方面，本说明书实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行计算机程序，所述计算机程序被处理器执行时实现如上述技术方案提供的方法步骤。

采用上述技术方案，本文提供的一种致密储层排驱压力确定方法、装置、设备和存储介质从排驱压力的定义出发，根据单位增压饱和度增量来确定致密储层的排驱压力值，结果更加准确，能够为储层研究和开发提供可靠的理论支撑；且求取过程简单便捷，提高了对致密储层排驱压力确定的效率。

为让本文的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。

附图说明

为了更清楚地说明本文实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本文的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1示出了本文实施例提供的一种致密储层排驱压力确定方法的流程示意图；

图2示出了一种确定排驱压力的方法流程示意图；

图3(a)至图3(c)示出了岩心样本的CT图像；

图4(a)和图4(b)示出了不同岩心样本的孔喉分布特征；

图5示出了渗透率与排驱压力的关系曲线图；

图6示出了岩心样本的毛管压力曲线示意图；

图7示出了单位增压饱和度增量与进汞压力的关系曲线；

图8示出了本实施例提供的一种致密储层排驱压力确定装置的结构示意图；

图9示出了本实施例提供的一种计算机设备的结构示意图。

附图符号说明：

10、获取模块；

20、计算模块；

30、选择模块；

902、计算机设备；

904、处理器；

906、存储器；

908、驱动机构；

910、输入/输出模块；

912、输入设备；

914、输出设备；

916、呈现设备；

918、图形用户接口；

920、网络接口；

922、通信链路；

924、通信总线。

具体实施方式

下面将结合本文实施例中的附图，对本文实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本文一部分实施例，而不是全部的实施例。基于本文中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本文保护的范围。

需要说明的是，本文的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本文的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、装置、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。

自1921年Washburn提出用非润湿相流体测量孔隙结构以来，压汞法测量多孔介质孔隙结构特征的压汞仪器已被广泛应用，以快速、可靠、测量范围宽阔等优势，满足不同行业工业化研究和应用。当前，最先进的压汞仪器可实现全自动在线测试，压力测量范围在1psi至60000psi，压力分辨率可控制在1psi以内；测量孔直径范围在950μm至3.6nm，测量范围与汞的特征、样品及孔隙结构特征等因素相关。压汞法测量孔径时的基本假设为：(1)孔隙是圆柱形的，运用毛细管束理论描述；(2)所有的孔都向外延伸至试样的外表面。

排驱压力(Displacement Pressure，简写为Pd或Pcd)，也称为阈压或门槛压力，它是指某一岩样中的润湿相流体被非润湿相流体开始排替所需的最低压力。排驱压力的大小，与孔隙喉道大小有直接关系，是研究储层储集性能和渗流性能的一个重要指标，排驱压力的确定对储层评价具有重要的意义。

排驱压力的确定方法大致可分为毛管压力曲线法、实验室直接测量法、地球物理评价法和影响因素综合分析法这四类。其中，毛管压力曲线切线法是最为常用的、且能快速确定排驱压力的方法：在半对数坐标中沿着毛管压力曲线平坦部分的第一个拐点做切线，切线延长线与纵坐标轴相交的压力点即为排驱压力。

然而，在致密油储层岩心压汞资料中，毛管压力曲线在开始进汞阶段大多没有明显的平坦部分，以曲线开始稳定时对应的毛细管压力来确定排驱压力的方法难以适用，很难精确地做出它的切线。因此，对毛管压力曲线做切线的方法具有较大的局限性。当拐点难以确定时，还可以以相邻两个测点进汞饱和度差大于等于3％回退一个测压点的压力作为排驱压力。但这种方法获得的排驱压力确定值通常高于实际值，具有较大误差；加之压汞法具有的滞后效应，使得获得的数据不能准确反映实际排驱压力特征。

为了解决上述问题，本文实施例提供一种致密储层排驱压力确定方法，能够更加方便、可靠的获得排驱压力。图1是本文实施例提供的一种致密储层排驱压力确定方法的步骤示意图，本说明书提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的系统或装置产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行。具体的如图1所示，所述方法可以包括：

S110：获取岩心样本的压汞数据，所述压汞数据包括多个进汞压力和各个进汞压力对应的进汞饱和度；

汞是常见流体中表面张力最大的一种，常利用这一特点进行压汞试验以获取岩样的孔喉特征和排驱压力。压汞试验时，获取多个测试点，每个测试点的横坐标和纵坐标分别为进汞饱和度和进汞压力，即进汞压力和进汞饱和度是一一对应的。

S120：根据所述进汞压力和所述进汞饱和度，计算相邻两个进汞压力对应的单位增压饱和度增量；压汞仪器可显示有各测试点对应的进汞压力值，进汞饱和度的计算可以采用本领域常规手段，例如可以在岩样预处理结束后，高压压汞实验开始前，根据天然气行业标准分别测得岩样的气测渗透率和气测孔隙度，借助于这两个参数可以确定进汞饱和度。

样本1采集的各测试点对应的进汞压力值和进汞饱和度见表1。

表1

从表1可以看出。1号样本的压汞数据测试点有30个，其中进汞饱和度为非零值的点有29个(即有效测试点)，测试点数和有效测试点数均比较充足。其中，各测试点进汞压力从大到小排列，分段进汞饱和度(％)是两个测试点之间进汞饱和度之差，则根据下述公式计算相邻两个测试点对应的单位增压饱和度增量：

SHg′|Pc＝△SHg/△Pc

其中，SHg′|Pc为单位增压饱和度增量，△SHg为饱和度增量，△Pc为进汞压力增量。单位增压饱和度增量的量纲为％/MPa。

例如，第2个测试点和第3个测试点的单位增压饱和度增量为29.910％/MPa，由计算得到。

S130：筛选满足筛选条件的单位增压饱和度增量对应的进汞压力值作为排驱压力值，可记为Pcd0。

本说明书实施例提供了一种新的确定致密储层排驱压力的方法：根据相邻两个测试点之间单位增压饱和度增量来确定排驱压力，可避免现有的排驱压力确定方法带来的排驱压力偏高、排驱压力与渗透率关系出现严重矛盾的情况，且本说明书实施例提供的致密储层排驱压力确定方法能适用于处理常规压汞和高压压汞的排驱压力确定，适用范围广。

在一些可行的实施例中，如图2所示，步骤S130：筛选满足筛选条件的单位增压饱和度增量对应的进汞压力值作为排驱压力值，进一步为：

S210：获取岩心样本的参数数据，所述参数数据包括孔喉分布特征和渗透率；

孔喉是孔隙和喉道的总称，孔隙是流体赋存于岩石中的基本储集空间，而喉道则是流体在岩石中渗流的重要通道；孔喉的大小和分布对渗透率有很大影响，孔喉越大流体越容易在储层中流动，而排驱压力更是与孔隙喉道大小有直接关系。

除此之外，所述参数数据还包括根据岩心CT图像、铸体薄片观测得到的数据，以及恒速压汞和离心核磁测得的数据。

如图3(a)至图3(c)所示是1号岩心样本的CT图像，图3(a)是岩心样本的在俯视视角下的CT图像，图3(b)是正视视角下岩心样本的CT图像，图3(c)是岩心样本的三维图像，CT图像可用于解释样本的孔喉特征。

如图4(a)和图4(b)所示是不同岩心样本的孔喉分布特征，图4(a)中，喉道的尺寸分布与孔隙的尺寸分布没有交集，说明孔隙的尺寸与喉道的尺寸差异较大；而图4(b)中，喉道的尺寸分布与孔隙的尺寸分布存在交集，即有些孔隙的尺寸与喉道的尺寸大小相当。从而在压汞试验时，压汞初始阶段汞进入尺寸较大的孔隙到逐渐占据尺寸较小的孔隙的过程中，随压力增加进汞速度会逐渐降低甚至降低至0(孔隙和喉道转换区)，最大喉道突破后，随进汞压力升高，进汞速度又快速增加。当孔隙的尺寸与喉道的尺寸具有较大差距时，进汞速度上的变化更加明显；总而言之，孔喉分布特征对压汞测试和排驱压力的确定具有指导意义。

S220：根据所述参数数据，确定排驱压力范围；

如图5所示，是一种典型的岩心样本渗透率与排驱压力的关系曲线图。渗透率受孔隙喉道的制约，较大尺寸的喉道控制渗透率大小，排驱压力对应的最大喉道更是占据重要的地位，通常情况下，排驱压力与渗透率具有较好的相关关系。渗透率与排驱压力的关系曲线可作为约束和准确确定致密储层排驱压力的参考依据。

测得1号样本的孔隙度和渗透率分别为：8.477％和0.0848mD。

根据上述参数数据，采用定量描述和统计方法获得致密储层岩心样本排驱压力的基本分布特征，以便确定最大喉道分布的上限，为准确确定排驱压力提供基础性约束。根据所述参数数据，可以计算得到1号岩心样本的排驱压力在0.5MPa至2.94MPa范围内。

S230：选择满足所述排驱压力范围的单位增压饱和度增量对应的进汞压力作为排驱压力值。

进一步地，步骤230：选择满足所述排驱压力范围的单位增压饱和度增量对应的进汞压力作为排驱压力值，进一步为：

S310：根据所述进汞压力值和所述进汞饱和度，在双对数坐标系下，绘制毛管压力曲线；如图6所示，是双对数坐标系下1号岩心样本的毛管压力曲线示意图。

S320：根据所述毛管压力曲线的形态和值域，将所述毛管压力曲线划出分麻皮效应区和喉道初始进汞区，还可以划分出快速进汞区和高压减速进汞区。

在进行压汞试验时，汞在岩心样本表面的粗糙坑凹处、拐角处和棱线上形成一些不贴合的汞蒸气空腔，这种现象就是“麻皮效应”。需要说明的是，麻皮效应区主要反映岩心的孔隙大小特征，缺少孔隙发育的岩心，可能看不到麻皮效应。麻皮效应越显著，代表孔隙越发育，储集能力越好。

随着压力逐渐增大，空腔会逐渐缩小最后被汞填满而消失。此时，汞并未真正进入孔喉系统。当压力进一步增大开始促使汞进入最大的孔喉系统中时，此时的压力被称为排驱压力，麻皮效应结束在排驱压力点。因此，可根据麻皮效应区的值域和/或喉道初始进汞缺的值域确定排驱压力修正值。

快速进汞区，主要反映连通喉道最发育的部分，一般接近喉道半径中值(非正态分布)或平均值(正态分布)；高压减速进汞区，随着进汞压力增高，喉道变得更为细微并且减少，因此，进汞速度降低。

根据表1中的压汞数据可以看出，1号样本的麻皮效应显著，麻皮效应区对应的进汞压力在0.01369MPa至0.4778MPa范围内，麻皮效应反映的主要孔隙半径应在1.5μm至53.7μm；喉道初始进汞区对应的进汞压力在0.11MPa至2.746MPa范围内。需要说明的是，此处麻皮效应区和喉道初始进汞区有部分重叠，初始进汞区与麻皮效应区本就是一个连续的压力区间，并且此时是在初定区间阶段，可以存在重叠区，待最终确定排驱压力后再清晰划分开即可。

S330：根据划分出的所述麻皮效应区和所述喉道初始进汞区，选择满足所述排驱压力范围的单位增压饱和度增量对应的进汞压力作为排驱压力值。

优选地，步骤S330进一步为：基于划分出的区域，选择位于排驱压力范围内且所述单位增压饱和度增量大于等于1％/MPa的进汞压力作为排驱压力值。

则根据表1中的压汞数据可知，有四个测试点(即测压点10至13)位于排驱压力范围(0.5MPa至2.94MPa)内；且喉道初始进汞区对应的进汞压力在0.11MPa至2.746MPa；测试点13的单位增压饱和度增量为1.300％满足大于等于1％/MPa要求，由于1.300％是测试点13和测试点12间的单位增压饱和度增量，因此确定测试点12对应的进汞压力为排驱压力，即1号样本的排驱压力为2.048MPa。

需要说明的是，本说明书实施例中，上述选择单位增压饱和度增量大于等于1％/MPa的进汞压力作为排驱压力值是指位于排驱压力范围内，符合喉道初始进汞区进汞压力区间，第一个突破1％/MPa的测试点的进汞压力值。

优选地，在一些实施例中，所述方法还包括：

对所述排驱压力值进行修正。

如表1中压汞数据可知，测试点14的进汞压力和测试点13的进汞压力之间差距较大，使得测试点13和测试点14之间的单位增压饱和度增量(5.357％/MPa)与测试点13与测试点12之间的单位增压饱和度增量(1.300％/MPa)相比差距较大，因此，优选地，对根据上述步骤获得的排驱压力进行修正。

在一些可行的实施例中，可以沿毛管压力曲线的形态和趋势增设新的测试点数据(插值法)，例如可确定为(3.25MPa，6.36％)；修正后排驱压力值仍为2.048MPa，即根据本说明书实施例提供的致密储层排驱压力确定方法得到的排驱压力值准确性较高。

当然，还可以通过一些其他的方法对获取到的排驱压力值进行修正。例如，切线法：根据毛管压力曲线中第一个拐点做切线，获取切线延长线与纵坐标的交点对应的压力值为排驱压力修正值。还可以通过饱和度增量3％回退法(或2％回退法)获得排驱压力修正值。进而根据排驱压力修正值修正所述排驱压力值。

根据表1中的压汞数据，根据现有的饱和度增量2％回退法或3％回退法，得到的排驱压力均为2.746MPa(即测试点14对应的进汞饱和度为7.364％，大于2％和3％)。然而7.364％远大于2％或3％；且2.746MPa作为排驱压力时，此时对应的最大孔喉半径为为0.267μm，低于有效致密油甜点标准，与孔隙度反映的甜点不一致，因此，2.746MPa偏高不宜作为排驱压力值。也就是说，根据本说明书实施例提供的致密储层排驱压力确定方法，与现有技术相比，能够提高排驱压力值确定的准确性。

如图7所示，还可以根据步骤S120计算得到的相邻两个测试点对应的单位增压饱和度增量，绘制单位增压饱和度增量与进汞压力的关系曲线。图中PCd0指示的点对应的进汞压力即为排驱压力。该曲线可以利用喉道突破后进汞饱和度的快速增加特征，发现和确定突破喉道所对应的压力，也就是排驱压力。

综上所述，本说明书实施例提供的一种致密储层排驱压力确定方法，是从“压汞突破最大连通孔隙的喉道开始连续进入岩石孔隙”这一排驱压力定义出发，根据单位增压饱和度增量，综合了孔喉分布特征、排驱压力与渗透率基本关系，以及麻皮效应、喉道初始进汞区等致密储层特征来确定排驱压力，结果更加准确，能够为储层研究和开发提供可靠的理论支撑；且求取过程简单便捷，有利于提高对致密储层排驱压力确定的效率。

在此基础上，还参考饱和度增量3％回退法等、切线法等现有方法确定的排驱压力，可对排驱压力值进行修正优化。

如表2所示，为2号岩心样本的压汞数据。并测得2号岩心样本的基础数据包括孔隙度10.17％，渗透率为0.137mD。

表2

从表2中的压汞数据可知，测压点总数为24个，汞饱和度非0值点仅11个，属于有效测试点数量偏少的情况。

通过孔隙度、渗透率等参数数据初步预估排驱压力范围应在0.5MPa至2.0MPa之间。观察该范围区间的测试点情况，仅有四个测试点在预期的排驱压力范围附近。计算饱和度随压力增加而增大的速度，有一部分麻皮效应特征显示，反映出相应孔隙半径对应压力0.523MPa，主要孔隙半径应在1.4μm以上。

计算相邻两个测试点件的单位增压饱和度增量，并结合麻皮效应区、喉道初始进汞区的值域，以及排驱压力范围，得到排驱压力值为1.4928MPa(由测试点16与测试点17之间的单位增压饱和度增量为6.311％/MPa确定得到，测试点13与测试点14之间的单位增压饱和度增量为2.318％/MPa，也大于1％/MPa，但该测试点位于麻皮效应区，因此该测试点的进汞压力不应作为排驱压力)。由于排驱压力附近测压点较稀疏，单位增压饱和度增量从0.563％/MPa到6.311％/MPa突变较大，因此，确定的排驱压力值可能仍存在不够精确的问题；实际的排驱压力应当在1.4928MPa之上略高一点位置，通过插值法对排驱压力进行修正，最终确定排驱压力的值为1.613MPa。

当样本的压汞数据在预估的排驱压力范围区间的分布不够均匀时，可以重新考虑布置测压点，增加进汞压力和进汞饱和度的数量，特别是使得排驱压力附近的测点密集，以便精确确定排驱压力值。

表3是3号岩心样本的压汞数据，并测得3号岩心样本的基础数据包括孔隙度10.76％，渗透率为0.421mD。

表3

从表3中的压汞数据可知，测试点总数为14个，通过孔隙度、渗透率等参数数据初步预估排驱压力范围应在0.25MPa至0.735MPa之间，则仅有三个测压点在该范围内，测压点明显偏少；计算饱和度随压力增加而增大的速度，没有麻皮效应特征显示。

根据本说明书实施例中，选择单位增压饱和度增量大于等于1％/MPa的进汞压力作为排驱压力值确定的排驱压力为0.294MPa。上述表格中有多个测试点既满足排驱压力范围，又满足单位增压饱和度增量满足大于等于1％/MPa(测试点1至测试点3)，但选择第一个满足该条件的测试点(即测试点1)对应的进汞压力值作为排驱压力值。

但由于压汞数据测压点较粗，也不精确，因此确定的排驱压力可能存在不精确的问题。按照现有的饱和度3％回退法确定的排驱压力为0.459Mpa，但该方法确定的排驱压力有些偏高，原因是进汞饱和度由2.577％升到9.74％跨度远远大于3％，进入到快速进汞区，这点从饱和度压力增量参数可以得到证实。

按照喉道半径分布特征，基本上满足对数正态分布，采用插值法对根据本说明书实施例得到的排驱压力(0.294MPa)进行修正，得到修正后的排驱压力为0.3679MPa，对应的最大喉道半径为1.99微米，综合地区规律和实验数据确定排驱压力为0.3679MPa。

因此，当压汞数据测试点数量不足时，相邻两个测试点间的进汞压力差、饱和度差均较大，则会影响根据进汞饱和度差回退法确定的排驱压力的精度，步骤S120：所述根据所述进汞压力和所述进汞饱和度，计算相邻两个进汞压力对应的单位增压饱和度增量，还包括：

增加进汞压力值和进汞饱和度(即测试点)的数量。

在根据本说明书实施例中提供的根据单位增压饱和度增量确定排驱压力值的基础上，采用插值法或其他方法对排驱压力值加以修正，能够适用于测试点数及有效测试点数较少情况下致密储层排驱压力的确定，满足在有效测试点数据匮乏情况下对排驱压力的确定需求。

综上所述，本说明书实施例提供的致密储层排驱压力确定方法利用了排驱压力的定义，孔隙喉道分布规律，致密油甜点特征，渗透率与排驱压力的地区规律等参数确定的排驱压力更为准确，能够适用于各类致密储层的排驱压力确定，且更加方便、可靠。

如图8所示，本说明书实施例还提供一种致密储层排驱压力确定装置，包括：

获取模块10，用于获取岩心样本的压汞数据，所述压汞数据包括多个进汞压力和各个进汞压力对应的进汞饱和度；

计算模块20，根据所述进汞压力和所述进汞饱和度，计算各进汞压力值对应的单位增压饱和度增量；

选择模块30，用于选择满足筛选条件的单位增压饱和度增量对应的进汞压力作为排驱压力值。

通过本说明书实施例提供的装置所取得的有益效果和上述方法所取得的有益效果相一致，此处不再赘述。

如图9所示，为本文实施例提供的一种计算机设备，所述计算机设备902可以包括一个或多个处理器904，诸如一个或多个中央处理单元(CPU)，每个处理单元可以实现一个或多个硬件线程。计算机设备902还可以包括任何存储器906，其用于存储诸如代码、设置、数据等之类的任何种类的信息。非限制性的，比如，存储器906可以包括以下任一项或多种组合：任何类型的RAM，任何类型的ROM，闪存设备，硬盘，光盘等。更一般地，任何存储器都可以使用任何技术来存储信息。进一步地，任何存储器可以提供信息的易失性或非易失性保留。进一步地，任何存储器可以表示计算机设备902的固定或可移除部件。在一种情况下，当处理器904执行被存储在任何存储器或存储器的组合中的相关联的指令时，计算机设备902可以执行相关联指令的任一操作。计算机设备902还包括用于与任何存储器交互的一个或多个驱动机构908，诸如硬盘驱动机构、光盘驱动机构等。

计算机设备902还可以包括输入/输出模块910(I/O)，其用于接收各种输入(经由输入设备912)和用于提供各种输出(经由输出设备914)。一个具体输出机构可以包括呈现设备916和相关联的图形用户接口918(GUI)。在其他实施例中，还可以不包括输入/输出模块910(I/O)、输入设备912以及输出设备914，仅作为网络中的一台计算机设备。计算机设备902还可以包括一个或多个网络接口920，其用于经由一个或多个通信链路922与其他设备交换数据。一个或多个通信总线924将上文所描述的部件耦合在一起。

通信链路922可以以任何方式实现，例如，通过局域网、广域网(例如，因特网)、点对点连接等、或其任何组合。通信链路922可以包括由任何协议或协议组合支配的硬连线链路、无线链路、路由器、网关功能、名称服务器等的任何组合。

对应于图1至图2中的方法，本文实施例还提供了一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时执行上述方法的步骤。

本文实施例还提供一种计算机可读指令，其中当处理器执行所述指令时，其中的程序使得处理器执行如图1至图2所示的方法。

应理解，在本文的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本文实施例的实施过程构成任何限定。

还应理解，在本文实施例中，术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系。例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。

本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本文的范围。

所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

在本文所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口、装置或单元的间接耦合或通信连接，也可以是电的，机械的或其它的形式连接。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本文实施例方案的目的。

另外，在本文各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。

所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本文的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本文各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

本文中应用了具体实施例对本文的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本文的方法及其核心思想；同时，对于本领域的一般技术人员，依据本文的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本文的限制。

Claims (7)

1.一种致密储层排驱压力确定方法，其特征在于，包括：

获取岩心样本的压汞数据，所述压汞数据包括多个进汞压力和各个进汞压力对应的进汞饱和度；

根据所述进汞压力和所述进汞饱和度，计算相邻两个进汞压力对应的单位增压饱和度增量；

获取岩心样本的参数数据，所述参数数据包括所述岩心样本的孔喉分布特征和渗透率；

根据所述参数数据，采用定量描述方法和统计方法获得致密储层岩心样本排驱压力的基本分布特征，以确定最大喉道分布的上限进而确定排驱压力范围；

根据所述进汞压力值和所述进汞饱和度，在双对数坐标系下，绘制毛管压力曲线；

根据所述毛管压力曲线的形态和值域，将所述毛管压力曲线划出分麻皮效应区和喉道初始进汞区；

基于划分出的麻皮效应区和喉道初始进汞区，选择位于排驱压力范围内、符合喉道初始进汞区进汞压力区间且所述单位增压饱和度增量第一个突破1%/MPa的测试点的进汞压力值作为排驱压力值。

2.根据权利要求1所述的方法，其特征在于，计算相邻两个进汞压力对应的单位增压饱和度增量的公式为：

SHg′|Pc=△SHg/△Pc

其中，SHg′|Pc为单位增压饱和度增量，△SHg为饱和度增量，△Pc为进汞压力增量。

3.根据权利要求1所述的方法，其特征在于，所述方法还包括：

对所述排驱压力值进行修正。

4.根据权利要求1所述的方法，其特征在于，所述根据所述进汞压力和所述进汞饱和度，计算相邻两个进汞压力对应的单位增压饱和度增量，还包括：

增加进汞压力和进汞饱和度的数量。

5.一种致密储层排驱压力确定装置，其特征在于，所述装置实施如权利要求1至4任意一项所述的致密储层排驱压力确定方法，所述装置包括：

获取模块，用于获取岩心样本的压汞数据，所述压汞数据包括多个进汞压力和各个进汞压力对应的进汞饱和度；

计算模块，根据所述进汞压力和所述进汞饱和度，计算相邻两个进汞压力对应的单位增压饱和度增量；

选择模块，用于获取岩心样本的参数数据，所述参数数据包括所述岩心样本的孔喉分布特征和渗透率；根据所述参数数据，采用定量描述方法和统计方法获得致密储层岩心样本排驱压力的基本分布特征，以确定最大喉道分布的上限进而确定排驱压力范围；根据所述进汞压力值和所述进汞饱和度，在双对数坐标系下，绘制毛管压力曲线；根据所述毛管压力曲线的形态和值域，将所述毛管压力曲线划出分麻皮效应区和喉道初始进汞区；基于划分出的麻皮效应区和喉道初始进汞区，选择位于排驱压力范围内、符合喉道初始进汞区进汞压力区间且所述单位增压饱和度增量第一个突破1%/MPa的测试点的进汞压力值作为排驱压力值。

6.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至4任意一项所述的方法步骤。

7.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行计算机程序，所述计算机程序被处理器执行时实现如权利要求1至4任意一项所述的方法步骤。