[发明专利]岩心孔隙度参数场的获取方法

说明书

本发明公开了一种岩心孔隙度参数场的获取方法，属于石油地球物理勘探领域。该方法包括：在k张岩心二值化图像上切割边长为L的方形截面，并将其对应的多个二维矩阵按三维顺序叠加，形成三维矩阵并封装至单元数组的各个元胞数组中，获取每个元胞数组对应的孔隙度，进而获得单元数组的孔隙度参数场。根据单元数组的孔隙度参数场计算三维切割样品的平均孔隙度模拟值，测定岩心样品的平均孔隙度实测值，当平均孔隙度模拟值和平均孔隙度实测值的误差小于5％时，确定L和k的取值为切割参数，并且所获取的单元数组的孔隙度参数场作为岩心样品的岩心孔隙度参数场。该方法能保证获得的岩心孔隙度参数场真实反映岩心样品孔隙度，且减少数据量。

技术领域

本发明涉及石油地球物理勘探领域，特别涉及一种岩心孔隙度参数场的获取方法。

背景技术

数字岩心是基于二维扫描电镜图像或三维CT扫描图像，运用计算机图像处理技术，通过一定的算法完成数字岩心重构，其能够最大程度地反映出地层信息，对于岩石物理的研究具有重要的意义。通过数字岩心所建立的孔隙介质模型的孔隙结构最大程度地与实际岩石相同，这不仅可用来精确地预测岩石的物理特性，还可以用来了解不同的物理响应之间的内在关系。所以，提供一种通过构建三维数字岩心来获取岩心孔隙度的方法十分必要。

现有技术(201510526748.x)提供了一种岩心的孔隙度矩阵反演方法，该方法包括以下步骤：利用X射线对岩心样品进行CT扫描，得到岩心CT图像；对所述岩心CT图像依次进行灰度化处理和对比度增强处理，得到岩心灰度图；对所述岩心灰度图进行二值化处理，得到多个二值化矩阵，并将所述多个二值化矩阵按三维顺序叠加成三维矩阵；将所述三维矩阵封装至单元数组的各个元胞数组中，通过对每个元胞数组的孔隙度进行局部反演，随后对所有元胞数组的孔隙度进行反演，实现岩心孔隙度矩阵反演，得到岩心孔隙度参数场。其中，所述的反演即为将图像转化为孔隙度数值矩阵的过程。该方法通过获取能精确表征岩心样品孔隙结构信息的三维矩阵，并将所得到的三维矩阵构建成多个元胞数值矩阵，并对其中的孔隙度进行反演，从而获取岩心孔隙度参数场。

发明人发现现有技术至少存在以下问题：

现有技术提供的方法得到的岩心孔隙度参数场数据量大，在利用其模拟计算岩心样品的平均孔隙度时其计算量大，同时将其作为岩心数值试验研究的基础数据进行相应计算时，其计算量也很大。

发明内容

本发明实施例所要解决的技术问题在于，提供了一种既能保证获得的岩心孔隙度参数场能真实反映岩心样品孔隙度特征，又能减少其平均孔隙度模拟计算量的岩心孔隙度参数场的获取方法。具体技术方案如下：

一种岩心孔隙度参数场的获取方法，包括：利用X射线沿岩心样品的轴向间隔预定距离对所述岩心样品进行CT扫描，得到k张岩心CT图像；对所述k张岩心CT图像依次进行灰度化处理、对比度增强处理、二值化处理，得到k张岩心二值化图像。进一步地，所述方法还包括：

步骤a、在k张所述岩心二值化图像上分别切割边长为L的方形截面，并将k张所述方形截面所对应的多个二维矩阵按三维顺序叠加，形成三维矩阵；

步骤b、将所述三维矩阵封装至单元数组的各个元胞数组中，通过对每个所述元胞数组中的孔隙度参数进行局部反演，获得每个所述元胞数组对应的孔隙度，进而获得所述单元数组的孔隙度参数场；

步骤c、根据所述单元数组的孔隙度参数场计算三维切割样品的平均孔隙度模拟值，测定所述岩心样品的平均孔隙度实测值，当所述平均孔隙度模拟值和所述平均孔隙度实测值之间的误差小于5％时，确定所述L和所述k的取值为切割参数，并且将基于所述切割参数所获取的单元数组的孔隙度参数场作为所述岩心样品的岩心孔隙度参数场。

具体地，所述步骤c还包括：确定所述平均孔隙度模拟值与所述平均孔隙度实测值之间的误差是否大于5％；

如果否，则确定所述L和所述k的原值为所述切割参数；