CN107329182A - 一种确定储层渗透率的方法及装置 - Google Patents

Abstract

本申请实施例公开了一种确定储层渗透率的方法及装置。所述方法包括：基于毛管压力曲线信息，确定孔喉类型划分界限；基于孔喉类型划分界限、岩心核磁T2谱和压汞孔喉半径分布曲线，确定与岩心核磁T2谱对应的孔隙分量划分界限；针对多个岩心中指定岩心，确定指定岩心的与孔隙分量划分界限对应的多种孔喉类型分别在指定岩心的核磁T2谱上的孔隙分量；建立指定岩心的渗透率与指定岩心的多个孔隙分量的关联关系；对采样点位置处的核磁测井T2谱进行油气校正，得到油气校正后的核磁测井T2谱；基于油气校正后的核磁测井T2谱和关联关系，确定采样点位置处的渗透率。本申请实施例提供的技术方案，可以提高所确定的储层渗透率的精度。

Description

一种确定储层渗透率的方法及装置

技术领域

本申请涉及储层评价技术领域，特别涉及一种确定储层渗透率的方法及装置。

背景技术

储层渗透率是衡量流体在压力差下通过多孔岩石有效孔隙能力的一种量值。准确求取储层渗透率是储层评价的基础。

Tim-Coates模型是确定储层渗透率最常用的方法之一。通过Tim-Coates模型可以建立储层渗透率K与可动流体体积、束缚水体积、孔隙度的关联关系。

然而，Tim-Coates模型仅将孔隙体积划分为可动流体体积和束缚水体积，难以刻画孔隙结构的复杂程度，从而可能影响储层渗透率的计算精度。

发明内容

本申请实施例的目的是提供一种确定储层渗透率的方法及装置，以提高所确定的储层渗透率的精度。

为解决上述技术问题，本申请实施例提供一种确定储层渗透率的方法及装置是这样实现的：

一种确定储层渗透率的方法，提供有目的储层在目标井位置处的核磁测井数据；其中，所述目的储层在目标井位置处包括多个采样点，所述核磁测井数据包括多个采样点位置处分别对应的核磁测井T2谱；还提供有所述目的储层中的多个岩心在饱和水状态下分别对应的岩心核磁T2谱、毛管压力曲线信息和压汞孔喉半径分布曲线；所述方法包括：

基于所述毛管压力曲线信息，确定孔喉类型划分界限；

基于所述孔喉类型划分界限、所述岩心核磁T2谱和所述压汞孔喉半径分布曲线，确定与所述岩心核磁T2谱对应的孔隙分量划分界限；

针对所述多个岩心中指定岩心，确定与所述孔隙分量划分界限对应的多种孔喉类型分别在所述指定岩心的岩心核磁T2谱上的孔隙分量；

建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系；

对所述采样点位置处的核磁测井T2谱进行油气校正，得到油气校正后的核磁测井T2谱；

基于所述油气校正后的核磁测井T2谱和所述关联关系，确定所述采样点位置处的渗透率。

优选方案中，所述建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系，包括：

获取所述指定岩心的渗透率和孔隙度，根据所述指定岩心的渗透率与所述指定岩心的多个孔隙分量，确定多条拟合直线以及与多条拟合直线分别对应的关联参数；其中，所述多条分段拟合直线中一条分段拟合直线与所述指定岩心的渗透率、所述指定岩心的一个孔隙分量相对应；

基于预设判别条件和所述关联参数，确定所述指定岩心的渗透率分别与所述指定岩心的多个孔隙分量之间的相关性；

根据所述相关性、所述指定岩心的多个孔隙分量和所述指定岩心的渗透率，建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系。

优选方案中，所述根据所述相关性、所述指定岩心的多个孔隙分量和所述指定岩心的渗透率，建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系，包括：

根据所述相关性，确定与所述指定岩心的渗透率、所述指定岩心的多个孔隙分量相关联的拟合模型，所述拟合模型中包括多个拟合参数；

基于所述指定岩心的多个孔隙分量和所述指定岩心的渗透率，利用最小二乘法对所述拟合模型进行拟合，以确定所述拟合模型中所述拟合参数的数值。

优选方案中，所述根据所述指定岩心的渗透率与所述指定岩心的多个孔隙分量，确定多条拟合直线以及与多条拟合直线分别对应的关联参数，包括：

对所述指定岩心的渗透率和所述指定岩心的多个孔隙分量进行线性拟合处理，确定所述多条分段拟合直线以及与所述多条分段拟合直线分别对应的关联参数；所述关联参数包括：斜率和线性拟合系数。

优选方案中，所述预设判别条件包括：

当一条分段拟合直线对应的关联参数中线性相关系数的绝对值大于或等于0.6且斜率大于0时，所述指定岩心的渗透率与该分段拟合直线对应的孔隙分量之间的相关性为正相关；或者，

当一条分段拟合直线对应的关联参数中线性相关系数的绝对值大于或等于0.6且斜率的小于0时，所述指定岩心的渗透率与该分段拟合直线对应的孔隙分量之间的相关性为负相关；或者，

当一条分段拟合直线对应的关联参数中线性相关系数的绝对值小于0.6时，所述指定岩心的渗透率与该分段拟合直线对应的孔隙分量之间的相关性为不相关。

优选方案中，所述对所述采样点位置处的核磁测井T2谱进行油气校正，包括：

基于预设T2截止值，确定所述核磁测井数据中指定采样点位置处的核磁测井T2谱对应的采样点束缚水饱和度，以及将所述采样点位置处的核磁测井T2谱划分为原始束缚水谱和原始可动流体谱；

基于所述预设T2截止值和所述岩心核磁T2谱，确定与所述岩心核磁T2谱对应的岩心束缚水饱和度和岩心T2几何平均值，并建立与所述岩心核磁T2谱对应的岩心束缚水饱和度和岩心T2几何平均值的关联关系；

根据所述采样点束缚水饱和度，以及所述岩心束缚水饱和度和岩心T2几何平均值的关联关系，确定与所述采样点束缚水饱和度对应的采样点T2几何平均值；

获取所述采样点位置处的孔隙度，构造所述采样点位置处的目标可动流体谱，以使得所述目标可动流体谱的振幅与所述原始束缚水谱的振幅之和等于所述采样点位置处的孔隙度，且所述采样点位置处的油气校正后的核磁测井T2谱对应的T2几何平均值等于所述采样点T2几何平均值；其中，所述油气校正后的核磁测井T2谱包括所述目标可动流体谱以及所述原始束缚水谱。

优选方案中，所述基于所述油气校正后的核磁测井T2谱和所述关联关系，确定所述采样点位置处的渗透率，包括：

针对所述多个采样点中指定采样点，确定所述指定采样点位置处的多个目标孔隙分量；其中，所述目标孔隙分量包括与所述孔隙分量划分界限对应的一种孔喉类型在所述油气校正后的核磁测井T2谱上的孔隙分量；其中，所述孔喉类型与所述目标孔隙分量一一对应；

基于所述目标孔隙分量和所述关联关系，确定所述指定采样点位置处的渗透率。

优选方案中，所述基于所述孔喉类型划分界限、所述岩心核磁T2谱和所述压汞孔喉半径分布曲线，确定与所述岩心核磁T2谱对应的孔隙分量划分界限，包括：

基于所述岩心核磁T2谱的主峰和所述压汞孔喉半径分布曲线的主峰，建立孔喉半径与T2弛豫时间之间的幂函数关系；

根据所述幂函数关系，将所述孔喉类型划分界限转换为与所述岩心核磁T2谱对应的孔隙分量划分界限。

优选方案中，所述基于所述毛管压力曲线信息，确定孔喉类型划分界限，包括：

基于所述毛管压力曲线信息，对所述岩心的孔隙结构进行分类，得到多种岩心孔隙结构类型；其中，所述毛管压力曲线信息包括：毛管压力曲线特征和进汞压力；

基于预设的进汞压力与孔喉半径的关联关系，以及与所述多种孔隙结构类型分别对应的毛管压力曲线中起始下凹处的进汞压力，确定多个孔喉半径分界值，并将所述多个孔喉半径分界值作为孔喉类型划分界限。

一种确定储层渗透率的装置，所述装置提供目的储层在目标井位置处的核磁测井数据；其中，所述目的储层在目标井位置处包括多个采样点，所述核磁测井数据包括多个采样点位置处分别对应的核磁测井T2谱；所述装置还提供有所述目的储层中的多个岩心在饱和水状态下分别对应的岩心核磁T2谱、毛管压力曲线和压汞孔喉半径分布曲线；所述装置包括：孔喉界限确定模块、孔隙分量界限确定模块、孔隙分量确定模块、关联关系建立模块、油气校正模块和渗透率确定模块；其中，

所述孔喉界限确定模块，用于基于所述毛管压力曲线，确定孔喉类型划分界限；

所述孔隙分量界限确定模块，用于基于所述孔喉类型划分界限、所述岩心核磁T2谱和所述压汞孔喉半径分布曲线，确定与所述岩心核磁T2谱对应的孔隙分量划分界限；

所述孔隙分量确定模块，用于针对所述多个岩心中指定岩心，确定与所述孔隙分量划分界限对应的多种孔喉类型分别在所述指定岩心的岩心核磁T2谱上的孔隙分量；

所述关联关系建立模块，用于建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系；

所述油气校正模块，用于对所述采样点位置处的核磁测井T2谱进行油气校正，得到油气校正后的核磁测井T2谱；

所述渗透率确定模块，用于基于所述油气校正后的核磁测井T2谱和所述关联关系，确定所述采样点位置处的渗透率。

本申请实施例提供了一种确定储层渗透率的方法及装置，基于毛管压力曲线信息，可以确定孔喉类型划分界限；基于孔喉类型划分界限、岩心核磁T2谱和压汞孔喉半径分布曲线，可以确定与岩心核磁T2谱对应的孔隙分量划分界限；针对多个岩心中指定岩心，可以确定指定岩心的与孔隙分量划分界限对应的多种孔喉类型分别在指定岩心的核磁T2谱上的孔隙分量；可以建立指定岩心的渗透率与指定岩心的多个孔隙分量的关联关系；可以对采样点位置处的核磁测井T2谱进行油气校正，得到油气校正后的核磁测井T2谱；基于油气校正后的核磁测井T2谱和关联关系，可以确定采样点位置处的渗透率。考虑了油气校正后的核磁测井T2谱上不同孔隙分量对渗透率的贡献，从而可以提高所确定的储层渗透率的精度。

附图说明

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

图1是本申请一种确定储层渗透率的方法实施例的流程图；

图2是本申请实施例中岩心孔隙结构类型对应的毛管压力曲线的示意图；

图3是本申请实施例中岩心铸体薄片图像；

图4是本申请实施例中孔喉半径与T2弛豫时间之间的交会示意图；

图5是本申请实施例中孔隙分量与渗透率的交会示意图；

图6是本申请实施例中分别采用不同方法所确定的渗透率结果的对比示意图；

图7是本申请确定储层渗透率的装置实施例的组成结构图。

具体实施方式

本申请实施例提供一种确定储层渗透率的方法及装置。

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

本申请实施例提供了一种确定储层渗透率的方法。所述方法提供有目的储层在目标井位置处的核磁测井数据；其中，所述目的储层在目标井位置处包括多个采样点，所述核磁测井数据包括多个采样点位置处分别对应的核磁测井T2谱。所述方法还提供有所述目的储层中的多个岩心在饱和水状态下分别对应的岩心核磁T2谱、毛管压力曲线信息和压汞孔喉半径分布曲线。

在本实施方式中，所述目的储层可以是储层渗透率尚未确定的储层。所述目的储层中的目标井位置可以指通过在所述目的储层中已经开设的钻井，能够检测到储层的核磁测井数据的位置。所述目的储层可以是油气勘探过程中需要勘探开发的储层。所述目的储层可以是碳酸盐岩储层。

在本实施方式中，可以通过核磁共振仪器在所述目的储层的目标井位置处进行核磁测井处理，得到所述目的储层在目标井位置处的核磁测井数据。

在本实施方式中，可以通过钻井的方式，获取所述目的储层中的多个岩心。还可以对所述多个岩心中每一个岩心在百分之白的饱和水状态下分别依次进行压汞实验、核磁共振实验、孔隙度测试和渗透率测试，从而可以得到所述目的储层中的多个岩心在饱和水状态下分别对应的岩心核磁T2谱、孔隙度、渗透率、毛管压力曲线信息和压汞孔喉半径分布曲线。

在本实施方式中，所述目的储层在目标井位置处可以包括多个采样点。所述核磁测井数据可以包括多个采样点位置处分别对应的核磁测井T2谱。

图1是本申请一种确定储层渗透率的方法实施例的流程图。如图1所示，所述方法，包括以下步骤。

步骤S101：基于所述毛管压力曲线信息，确定孔喉类型划分界限。

在本实施方式中，基于所述毛管压力曲线信息，确定孔喉类型划分界限，具体可以包括，可以基于所述毛管压力曲线信息，对所述岩心的孔隙结构进行分类，可以得到多种岩心孔隙结构类型。其中，所述毛管压力曲线信息可以包括：毛管压力曲线特征和进汞压力。可以基于预设的进汞压力与孔喉半径的关联关系，以及与所述多种孔隙结构类型分别对应的毛管压力曲线中起始下凹处的进汞压力，可以确定多个孔喉半径分界值，并将所述多个孔喉半径分界值作为孔喉类型划分界限。

例如，可以根据所述多个岩心中每一个岩心制作对应的铸体薄片样品，并获取所述铸体薄片样品的铸体薄片图像。基于所述岩心的铸体薄片样品的铸体薄片图像，以及所述毛管压力曲线信息中毛管压力曲线特征和进汞压力，可以得到四种岩心孔隙结构类型，分别为Ⅰ类孔隙结构、Ⅱ类孔隙结构、Ⅲ类孔隙结构和Ⅳ类孔隙结构。

图2中(a)、(b)、(c)和(d)分别是所述四种岩心孔隙结构类型对应的毛管压力曲线的示意图。图2中横坐标和纵坐标分别表示进汞压力Pc和非润湿相流体饱和度S的log函数值。其中，所述进汞压力可以反映孔隙结构中孔喉半径的大小，进汞压力越大，孔喉半径越小。所述非润湿相流体饱和度可以反映各个所述孔喉之间的连通性。图3中(a)、(b)、(c)和(d)分别是所述四种岩心孔隙结构类型对应的岩心铸体薄片图像。如图2和图3中(a)、(b)、(c)和(d)的方框区域所示，所述四种孔隙结构的毛管压力曲线特征，以及铸体薄片图像中铸体薄片视域内面孔率和孔隙类型分别为：

(1)所述Ⅰ类孔隙结构对应的毛管压力曲线特征为毛管压力曲线开口向下，平直段在所述四种岩心孔隙结构类型分别对应的毛管压力曲线特征的平直段中最短，无明显下凹曲线段；所述Ⅰ类孔隙结构对应的孔喉半径在所述四种岩心孔隙结构对应的孔喉半径中最大；所述Ⅰ类孔隙结构对应的铸体薄片视域内面孔率较大，孔隙以铸模孔、粒内孔为主；

(2)所述Ⅱ类孔隙结构对应的毛管压力曲线特征为毛管压力曲线开口向上，平直段较短，有明显下凹曲线段；所述Ⅱ类孔隙结构对应的孔喉半径较大；所述Ⅱ类孔隙结构对应的铸体薄片视域内面孔率较大，孔隙以铸模孔、粒内孔为主，但有部分微孔；

(3)所述Ⅲ类孔隙结构对应的毛管压力曲线特征为毛管压力曲线开口向上，平直段较长，有明显下凹曲线段，但下凹曲度比所述Ⅱ类孔隙结构对应的下凹曲线段的下凹曲度小；所述Ⅲ类孔隙结构对应的孔喉半径较大；所述Ⅲ类孔隙结构对应的铸体薄片视域内面孔率较小，孔隙以铸模孔、粒内孔和微孔为主；

(4)所述Ⅳ类孔隙结构对应的毛管压力曲线特征为毛管压力曲线开口向下，平直段最长，有下凹曲线段；所述Ⅳ类孔隙结构对应的孔喉半径较小；所述Ⅳ类孔隙结构对应的铸体薄片视域内孔隙以微孔为主。

在本实施方式中，可以采用下述公式表征所述预设的进汞压力与孔喉半径的关联关系：

其中，Pc表示进汞压力，r表示孔喉半径，σ表示汞与空气之间的界面张力，通常取值为480达因/厘米(dyn/cm)，θ表示汞与储层中岩石的润湿角，通常取值为140度。可以基于预设的进汞压力与孔喉半径的关联关系，以及与所述四种孔隙结构类型分别对应的毛管压力曲线中起始下凹处的进汞压力，可以确定三个孔喉半径分界值r₁、r₂和r₃，分别为0.15微米(μm)、1μm和5μm，即微喉、细喉、中喉和粗喉，四种孔喉类型的孔喉半径的取值范围分别为：r＜r₁、r₁≤r≤r₂、r₂＜r≤r₃和r＜r₃。并将所述多个孔喉半径分界值作为孔喉类型划分界限。

步骤S102：基于所述孔喉类型划分界限、所述岩心核磁T2谱和所述压汞孔喉半径分布曲线，确定与所述岩心核磁T2谱对应的孔隙分量划分界限。

在本实施方式中，确定与所述岩心核磁T2谱对应的孔隙分量划分界限，具体可以包括，基于所述岩心核磁T2谱的主峰和所述压汞孔喉半径分布曲线的主峰，可以建立孔喉半径与T2弛豫时间之间的幂函数关系。根据所述幂函数关系，可以将所述孔喉类型划分界限转换为与所述岩心核磁T2谱对应的孔隙分量划分界限。

例如，图4是本申请实施例中孔喉半径与T2弛豫时间之间的交会示意图。可以根据图4所示的多个数据点，每个数据点均对应一个岩心的孔喉半径与T2弛豫时间，可以采用最小二乘拟合方法，建立孔喉半径与T2弛豫时间之间的幂函数关系。根据所述幂函数关系，可以将所述孔喉类型划分界限转换为与所述岩心核磁T2谱对应的孔隙分量划分界限。所述孔隙分量划分界限可以为由弛豫时间分界值t₁、t₂和t₃构成，分别为30毫秒(ms)、90ms、200ms。

步骤S103：针对所述多个岩心中指定岩心，确定所述指定岩心的与所述孔隙分量划分界限对应的多种孔喉类型分别在所述指定岩心的核磁T2谱上的孔隙分量。

在本实施方式中，可以采用下述公式分别计算所述指定岩心的与所述孔隙分量划分界限对应的多种孔喉类型分别在所述指定岩心的核磁T2谱上的孔隙分量为：

其中，S1、S2、S3和S4分别表示微喉、细喉、中喉和粗喉在所述指定岩心的核磁T2谱上的孔隙分量；t₁、t₂和t₃分别表示孔隙分量S1和S2、S2和S3、S3和S4之间的驰豫时间分界值，t₀和t_e分别表示所述指定岩心的核磁T2谱的驰豫时间起始值和驰豫时间终止值；A表示所述指定岩心的核磁T2谱的幅度值。所述指定岩心可以为所述多个岩心中的任一岩心。

步骤S104：建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系。

在本实施方式中，建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系，具体可以包括，可以获取所述指定岩心的渗透率和孔隙度。根据所述指定岩心的渗透率与所述指定岩心的多个孔隙分量，可以确定多条拟合直线以及与多条拟合直线分别对应的关联参数。其中，所述多条分段拟合直线中一条分段拟合直线与所述指定岩心的渗透率、所述指定岩心的一个孔隙分量相对应。基于预设判别条件和所述关联参数，可以确定所述指定岩心的渗透率分别与所述指定岩心的多个孔隙分量之间的相关性。根据所述相关性、所述指定岩心的多个孔隙分量和所述指定岩心的渗透率，可以建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系。

在本实施方式中，所述预设判别条件可以包括：当一条分段拟合直线对应的关联参数中线性相关系数的绝对值大于或等于0.6且斜率大于0时，所述指定岩心的渗透率与该分段拟合直线对应的孔隙分量之间的相关性为正相关；或者，当一条分段拟合直线对应的关联参数中线性相关系数的绝对值大于或等于0.6且斜率的小于0时，所述指定岩心的渗透率与该分段拟合直线对应的孔隙分量之间的相关性为负相关；或者，当一条分段拟合直线对应的关联参数中线性相关系数的绝对值小于0.6时，所述指定岩心的渗透率与该分段拟合直线对应的孔隙分量之间的相关性为不相关。

在本实施方式中，根据所述相关性、所述指定岩心的多个孔隙分量和所述指定岩心的渗透率，建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系，具体可以包括，根据所述相关性，可以确定与所述指定岩心的渗透率、所述指定岩心的多个孔隙分量相关联的拟合模型，所述拟合模型包括多个拟合参数。基于所述指定岩心的多个孔隙分量和所述指定岩心的渗透率，利用最小二乘法对所述拟合模型进行拟合，以确定所述拟合模型中所述拟合参数的数值。

在本实施方式中，根据所述指定岩心的渗透率与所述指定岩心的多个孔隙分量，确定多条拟合直线以及与多条拟合直线分别对应的关联参数，具体可以包括，可以对所述指定岩心的渗透率和所述指定岩心的多个孔隙分量进行线性拟合处理，确定所述多条分段拟合直线以及与所述多条分段拟合直线分别对应的关联参数。其中，所述关联参数可以包括：斜率和线性拟合系数。

例如，图5中(a)、(b)、(c)和(d)分别是所述指定岩心的四个孔隙分量与所述指定岩心的渗透率的交会示意图。图5中(a)、(b)、(c)和(d)的横坐标分别表示微喉、细喉、中喉和粗喉对应的孔隙分量，即微孔占比、小孔占比、中孔占比和大孔占比，纵坐标表示岩心渗透率。可以对所述指定岩心的渗透率和所述指定岩心的四个孔隙分量进行线性拟合处理，可以确定这四条分段拟合直线以及与所述多条分段拟合直线分别对应的关联参数，即孔隙分量S1和S2对应的斜率小于0，孔隙分量S3和S4对应的斜率大于0，孔隙分量S1、S2、S3和S4的线性拟合系数均大于0.6。如此，可以采用下述公式表征所述拟合模型：

其中，K表示所述指定岩心的渗透率，表示所述指定岩心的孔隙度，a、b、c、d和e均为拟合参数。所述拟合参数可以通过最小二乘拟合方法来确定。

步骤S105：对所述采样点位置处的核磁测井T2谱进行油气校正，得到油气校正后的核磁测井T2谱。

在本实施方式中，对所述采样点位置处的核磁测井T2谱进行油气校正，具体可以包括，基于预设T2截止值，可以确定所述核磁测井数据中指定采样点位置处的核磁测井T2谱对应的采样点束缚水饱和度，以及将所述采样点位置处的核磁测井T2谱划分为原始束缚水谱和原始可动流体谱。基于所述预设T2截止值和所述岩心核磁T2谱，可以确定与所述岩心核磁T2谱对应的岩心束缚水饱和度和岩心T2几何平均值，并建立与所述岩心核磁T2谱对应的岩心束缚水饱和度和岩心T2几何平均值的关联关系。根据所述采样点束缚水饱和度，以及所述岩心束缚水饱和度和岩心T2几何平均值的关联关系，可以确定与所述采样点束缚水饱和度对应的采样点T2几何平均值。可以获取所述采样点位置处的孔隙度。可以采用波动函数构造所述采样点位置处的目标可动流体谱，以使得所述目标可动流体谱的振幅与所述原始束缚水谱的振幅之和等于所述采样点位置处的孔隙度，且所述采样点位置处的油气校正后的核磁测井T2谱对应的T2几何平均值等于所述采样点T2几何平均值。其中，所述油气校正后的核磁测井T2谱包括所述目标可动流体谱以及所述原始束缚水谱。例如，可以采用下述公式表征所述波动函数：

其中，y表示所述波动函数的函数值，即振幅，t表示弛豫时间，ω表示角频率，和表示相位，A₁和A₂表示振幅系数。

步骤S106：基于所述油气校正后的核磁测井T2谱和所述关联关系，确定所述采样点位置处的渗透率。

在本实施方式中，基于所述油气校正后的核磁测井T2谱和所述关联关系，确定所述采样点位置处的渗透率，具体可以包括，针对所述多个采样点中指定采样点，确定所述指定采样点位置处的多个目标孔隙分量。其中，所述目标孔隙分量为与所述孔隙分量划分界限对应的一种孔喉类型在所述油气校正后的核磁测井T2谱上的孔隙分量其中，所述孔喉类型与所述目标孔隙分量一一对应。基于所述目标孔隙分量和所述关联关系，可以确定所述指定采样点位置处的渗透率。

例如，图6是本申请实施例中分别采用不同方法所确定的渗透率结果的对比示意图。图6中的渗透率曲线对应的横坐标分别为新模型渗透率、KTIM和KSDR，纵坐标表示地层深度。其中，KTIM和KSDR分别表示采用Tim-Coates模型和SDR模型确定的渗透率。图6中黑色实线、灰色实线和灰色虚线分别表示采用新模型(即本申请方法)、SDR模型和Tim-Coates模型确定的目的储层在目标井位置处的各个采样点位置处的渗透率结果。图6中多条灰色直线表示采用实验室岩心渗透率测试方法确定的渗透率结果。不仅如此，图6中还包括所述油气校正后的核磁测井T2谱，以及与油气校正后的核磁测井T2谱对应不同孔喉类型的孔隙分量，即微孔占比S1、小孔占比S2、中孔占比S3和大孔占比S4。可以分别计算采用本申请方法所确定的渗透率结果、采用SDR模型所确定的渗透率结果和采用Tim-Coates模型所确定的渗透率结果分别与采用实验室岩心渗透率测试方法所确定的渗透率结果之间的皮尔森(Pearson)相关系数，计算结果分别为0.69、0.53和0.56。其中，Pearson相关系数越接近1，表示相关程度越高。同时，如图5所示的渗透率对比结果，可以看出，采用本申请方法所确定的渗透率结果的精度较高。

所述确定储层渗透率的方法实施例，基于毛管压力曲线信息，可以确定孔喉类型划分界限；基于孔喉类型划分界限、岩心核磁T2谱和压汞孔喉半径分布曲线，可以确定与岩心核磁T2谱对应的孔隙分量划分界限；针对多个岩心中指定岩心，可以确定指定岩心的与孔隙分量划分界限对应的多种孔喉类型分别在指定岩心的核磁T2谱上的孔隙分量；可以建立指定岩心的渗透率与指定岩心的多个孔隙分量的关联关系；可以对采样点位置处的核磁测井T2谱进行油气校正，得到油气校正后的核磁测井T2谱；基于油气校正后的核磁测井T2谱和关联关系，可以确定采样点位置处的渗透率。考虑了油气校正后的核磁测井T2谱上不同孔隙分量对渗透率的贡献，从而可以提高所确定的储层渗透率的精度。

图7是本申请确定储层渗透率的装置实施例的组成结构图。所述装置提供目的储层在目标井位置处的核磁测井数据；其中，所述目的储层在目标井位置处包括多个采样点，所述核磁测井数据包括多个采样点位置处分别对应的核磁测井T2谱；所述装置还提供有所述目的储层中的多个岩心在饱和水状态下分别对应的岩心核磁T2谱、毛管压力曲线和压汞孔喉半径分布曲线。如图7所示，所述装置可以包括：孔喉界限确定模块100、孔隙分量界限确定模块200、孔隙分量确定模块300、关联关系建立模块400、油气校正模块500和渗透率确定模块600。

所述孔喉界限确定模块100，可以用于基于所述毛管压力曲线，确定孔喉类型划分界限。

所述孔隙分量界限确定模块200，用于基于所述孔喉类型划分界限、所述岩心核磁T2谱和所述压汞孔喉半径分布曲线，确定与所述岩心核磁T2谱对应的孔隙分量划分界限。

所述孔隙分量确定模块300，可以用于针对所述多个岩心中指定岩心，确定与所述孔隙分量划分界限对应的多种孔喉类型分别在所述指定岩心的岩心核磁T2谱上的孔隙分量。

所述关联关系建立模块400，可以用于建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系。

所述油气校正模块500，可以用于对所述采样点位置处的核磁测井T2谱进行油气校正，得到油气校正后的核磁测井T2谱。

所述渗透率确定模块600，可以用于基于所述油气校正后的核磁测井T2谱和所述关联关系，确定所述采样点位置处的渗透率。

所述确定储层渗透率的装置实施例与所述确定储层渗透率的方法实施例相对应，可以实现所述确定储层渗透率的方法实施例，并取得方法实施例的技术效果。

在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。

本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

上述实施例阐明的装置、模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。

为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。

通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中，内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括短暂电脑可读媒体(transitory media)，如调制的数据信号和载波。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (10)

1.一种确定储层渗透率的方法，其特征在于，提供有目的储层在目标井位置处的核磁测井数据；其中，所述目的储层在目标井位置处包括多个采样点，所述核磁测井数据包括多个采样点位置处分别对应的核磁测井T2谱；还提供有所述目的储层中的多个岩心在饱和水状态下分别对应的岩心核磁T2谱、毛管压力曲线信息和压汞孔喉半径分布曲线；所述方法包括：

基于所述毛管压力曲线信息，确定孔喉类型划分界限；

基于所述孔喉类型划分界限、所述岩心核磁T2谱和所述压汞孔喉半径分布曲线，确定与所述岩心核磁T2谱对应的孔隙分量划分界限；

针对所述多个岩心中指定岩心，确定与所述孔隙分量划分界限对应的多种孔喉类型分别在所述指定岩心的岩心核磁T2谱上的孔隙分量；

建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系；

对所述采样点位置处的核磁测井T2谱进行油气校正，得到油气校正后的核磁测井T2谱；

基于所述油气校正后的核磁测井T2谱和所述关联关系，确定所述采样点位置处的渗透率。

2.根据权利要求1所述的一种确定储层渗透率的方法，其特征在于，所述建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系，包括：

获取所述指定岩心的渗透率和孔隙度，根据所述指定岩心的渗透率与所述指定岩心的多个孔隙分量，确定多条拟合直线以及与多条拟合直线分别对应的关联参数；其中，所述多条分段拟合直线中一条分段拟合直线与所述指定岩心的渗透率、所述指定岩心的一个孔隙分量相对应；

基于预设判别条件和所述关联参数，确定所述指定岩心的渗透率分别与所述指定岩心的多个孔隙分量之间的相关性；

根据所述相关性、所述指定岩心的多个孔隙分量和所述指定岩心的渗透率，建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系。

3.根据权利要求2所述的一种确定储层渗透率的方法，其特征在于，所述根据所述相关性、所述指定岩心的多个孔隙分量和所述指定岩心的渗透率，建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系，包括：

根据所述相关性，确定与所述指定岩心的渗透率、所述指定岩心的多个孔隙分量相关联的拟合模型，所述拟合模型中包括多个拟合参数；

基于所述指定岩心的多个孔隙分量和所述指定岩心的渗透率，利用最小二乘法对所述拟合模型进行拟合，以确定所述拟合模型中所述拟合参数的数值。

4.根据权利要求2所述的一种确定储层渗透率的方法，其特征在于，所述根据所述指定岩心的渗透率与所述指定岩心的多个孔隙分量，确定多条拟合直线以及与多条拟合直线分别对应的关联参数，包括：

对所述指定岩心的渗透率和所述指定岩心的多个孔隙分量进行线性拟合处理，确定所述多条分段拟合直线以及与所述多条分段拟合直线分别对应的关联参数；所述关联参数包括：斜率和线性拟合系数。

5.根据权利要求4所述的一种确定储层渗透率的方法，其特征在于，所述预设判别条件包括：

当一条分段拟合直线对应的关联参数中线性相关系数的绝对值大于或等于0.6且斜率大于0时，所述指定岩心的渗透率与该分段拟合直线对应的孔隙分量之间的相关性为正相关；或者，

当一条分段拟合直线对应的关联参数中线性相关系数的绝对值大于或等于0.6且斜率的小于0时，所述指定岩心的渗透率与该分段拟合直线对应的孔隙分量之间的相关性为负相关；或者，

当一条分段拟合直线对应的关联参数中线性相关系数的绝对值小于0.6时，所述指定岩心的渗透率与该分段拟合直线对应的孔隙分量之间的相关性为不相关。

6.根据权利要求1所述的一种确定储层渗透率的方法，其特征在于，所述对所述采样点位置处的核磁测井T2谱进行油气校正，包括：

基于预设T2截止值，确定所述核磁测井数据中指定采样点位置处的核磁测井T2谱对应的采样点束缚水饱和度，以及将所述采样点位置处的核磁测井T2谱划分为原始束缚水谱和原始可动流体谱；

基于所述预设T2截止值和所述岩心核磁T2谱，确定与所述岩心核磁T2谱对应的岩心束缚水饱和度和岩心T2几何平均值，并建立与所述岩心核磁T2谱对应的岩心束缚水饱和度和岩心T2几何平均值的关联关系；

根据所述采样点束缚水饱和度，以及所述岩心束缚水饱和度和岩心T2几何平均值的关联关系，确定与所述采样点束缚水饱和度对应的采样点T2几何平均值；

获取所述采样点位置处的孔隙度，构造所述采样点位置处的目标可动流体谱，以使得所述目标可动流体谱的振幅与所述原始束缚水谱的振幅之和等于所述采样点位置处的孔隙度，且所述采样点位置处的油气校正后的核磁测井T2谱对应的T2几何平均值等于所述采样点T2几何平均值；其中，所述油气校正后的核磁测井T2谱包括所述目标可动流体谱以及所述原始束缚水谱。

7.根据权利要求1所述的一种确定储层渗透率的方法，其特征在于，所述基于所述油气校正后的核磁测井T2谱和所述关联关系，确定所述采样点位置处的渗透率，包括：

针对所述多个采样点中指定采样点，确定所述指定采样点位置处的多个目标孔隙分量；其中，所述目标孔隙分量包括与所述孔隙分量划分界限对应的一种孔喉类型在所述油气校正后的核磁测井T2谱上的孔隙分量；其中，所述孔喉类型与所述目标孔隙分量一一对应；

基于所述目标孔隙分量和所述关联关系，确定所述指定采样点位置处的渗透率。

8.根据权利要求1所述的一种确定储层渗透率的方法，其特征在于，所述基于所述孔喉类型划分界限、所述岩心核磁T2谱和所述压汞孔喉半径分布曲线，确定与所述岩心核磁T2谱对应的孔隙分量划分界限，包括：

基于所述岩心核磁T2谱的主峰和所述压汞孔喉半径分布曲线的主峰，建立孔喉半径与T2弛豫时间之间的幂函数关系；

根据所述幂函数关系，将所述孔喉类型划分界限转换为与所述岩心核磁T2谱对应的孔隙分量划分界限。

9.根据权利要求1所述的一种确定储层渗透率的方法，其特征在于，所述基于所述毛管压力曲线信息，确定孔喉类型划分界限，包括：

基于所述毛管压力曲线信息，对所述岩心的孔隙结构进行分类，得到多种岩心孔隙结构类型；其中，所述毛管压力曲线信息包括：毛管压力曲线特征和进汞压力；

基于预设的进汞压力与孔喉半径的关联关系，以及与所述多种孔隙结构类型分别对应的毛管压力曲线中起始下凹处的进汞压力，确定多个孔喉半径分界值，并将所述多个孔喉半径分界值作为孔喉类型划分界限。

10.一种确定储层渗透率的装置，其特征在于，所述装置提供目的储层在目标井位置处的核磁测井数据；其中，所述目的储层在目标井位置处包括多个采样点，所述核磁测井数据包括多个采样点位置处分别对应的核磁测井T2谱；所述装置还提供有所述目的储层中的多个岩心在饱和水状态下分别对应的岩心核磁T2谱、毛管压力曲线和压汞孔喉半径分布曲线；所述装置包括：孔喉界限确定模块、孔隙分量界限确定模块、孔隙分量确定模块、关联关系建立模块、油气校正模块和渗透率确定模块；其中，

所述孔喉界限确定模块，用于基于所述毛管压力曲线，确定孔喉类型划分界限；

所述孔隙分量界限确定模块，用于基于所述孔喉类型划分界限、所述岩心核磁T2谱和所述压汞孔喉半径分布曲线，确定与所述岩心核磁T2谱对应的孔隙分量划分界限；

所述孔隙分量确定模块，用于针对所述多个岩心中指定岩心，确定与所述孔隙分量划分界限对应的多种孔喉类型分别在所述指定岩心的岩心核磁T2谱上的孔隙分量；

所述关联关系建立模块，用于建立所述指定岩心的渗透率与所述指定岩心的多个孔隙分量的关联关系；

所述油气校正模块，用于对所述采样点位置处的核磁测井T2谱进行油气校正，得到油气校正后的核磁测井T2谱；

所述渗透率确定模块，用于基于所述油气校正后的核磁测井T2谱和所述关联关系，确定所述采样点位置处的渗透率。