基于模拟压汞法获取岩石的孔径分布的方法和装置制造方法及图纸

本申请实施例提供一种基于模拟压汞法获取岩石的孔径分布的方法和装置，该方法包括：获取待测定岩石的岩心的至少部分的三维数字图像，根据该三维数字图像，获取岩心的至少部分包括的孔隙的骨架图像，获取该骨架图像包括的所有点各自对应的孔径值，点对应的孔径值为点距离孔隙边界的最小距离的两倍，根据骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布。本申请的方法，可测定的孔隙尺寸的范围不受限制，可获取各种待测定岩石的孔隙的孔径分布。

Method and device for obtaining pore size distribution of rock based on simulated mercury injection method

The embodiment of the present application provides a method and apparatus for acquiring pore size distribution of rocks based on an analog mercury intrusion method, which includes acquiring a three-dimensional digital image of at least a part of the core to be measured, acquiring a skeleton image of at least a part of the pores included in the core from the three-dimensional digital image, and acquiring the skeleton image. The pore size of each point is twice the minimum distance between the point and the pore boundary. According to the pore size of each point in the skeleton image, the pore size distribution of the rock to be measured is obtained by simulating mercury intrusion method. The method of the present application can be used to obtain the pore size distribution of various rocks to be determined without restriction on the range of pore size to be measured.

【技术实现步骤摘要】
基于模拟压汞法获取岩石的孔径分布的方法和装置
本申请实施例涉及石油天然气地质勘探技术，尤其涉及一种基于模拟压汞法获取岩石的孔径分布的方法和装置。
技术介绍
储集层岩石的微观孔隙结构特征(孔径的大小、分布、孔隙的几何形状及相互联通情况)直接影响储集层的油气资源的聚集、流动和产出，是评价储集层资源潜力的关键参数之一。储集层岩石的孔径分布对油气的勘探开发至关重要，直接关系到储集层岩石的储油和储气能力、油和气的赋存状态、运移机理和能力以及储集层岩石的沉积演化过程，最终影响油气开采的难易程度和勘探开发价值。因此精确厘定储集层岩石的孔径分布特征意义重大。只有准确了解了储集层岩石的孔径分布特征才能正确反应储集层岩石的渗流能力，并制定出正确的勘探开发方案。目前用来表征储集层孔隙特征的方法主要有：吸附法、压汞法、扫描电镜、光学显微镜、CT扫描、薄片分析等方法。不同的表征方法可监测的孔隙范围大不相同，其中，氮气吸附法可以测定小于100nm孔径的孔隙，压汞法可以测定的孔径范围为3nm-100um，上述两种方法对可测定的孔隙的孔径具有一定的限制。扫描电镜和CT扫描等成像方法虽然可以直接观察孔隙结构，但无法获知孔隙的尺寸特征。
技术实现思路
本申请实施例提供一种基于模拟压汞法获取岩石的孔径分布的方法和装置，可测定的孔隙尺寸的范围不受限制，可获取各种待测定岩石的孔隙的孔径分布。第一方面，本申请提供一种基于模拟压汞法获取岩石的孔径分布的方法，其特征在于，包括：获取待测定岩石的岩心的至少部分的三维数字图像；根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像；获取所述骨架图像包括的各点对应的孔径值；其中，点对应的孔径值为所述点距离孔隙边界的最小距离的两倍；根据所述骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布。在一种可能的设计中，所根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像，包括：对所述三维数字图像进行图像分割，得到所述岩心的至少部分包括的孔隙的图像；根据所述孔隙的图像，获取所述骨架图像。在一种可能的设计中，所述根据所述孔隙的图像，获取所述骨架图像，包括：对所述孔隙的图像进行距离变换，得到第一图像；对所述第一图像进行三维细化处理，得到所述骨架图像。在一种可能的设计中，根据所述骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布，包括：确定注汞端面上所有的第一点对应的孔径值中，最大的第一孔径值，以及所述各点对应的孔径值中，最小的第二孔径值；所述注汞端面为所述骨架图像的第一底面；从所述第一孔径值和所述第二孔径值之间的孔径值中确定出M个目标孔径值，所述M个目标孔径值包括所述第一孔径值和所述第二孔径值；对于每个目标孔径值，通过模拟注汞法，获取所述孔径值对应的进汞体积和，所述进汞体积和为相应孔径值对应的孔隙体积和；根据所述M个目标孔径值和每个目标孔径值对应的孔隙体积和，得到所述待测定岩石的孔隙的孔径分布。在一种可能的设计中，所述通过模拟注汞法，获取所述目标孔径值对应的进汞体积和，包括：对于所述注汞端面的每个第一点，获取包括所述第一点的至少一条第一骨架线；对于每条第一骨架线，确定每条第一骨架线包括的至少一个目标点，所述目标点为第一点和第二点之间的点，所述第二点为在第一方向上，所述第一骨架线上第一个对应的孔径值小于所述目标孔径值的点；所述第一方向为从所述注汞端面指向所述骨架图像的第二底面的方向，至少一个目标点包括第一点和所述第二点；根据所有第一骨架线包括的所有目标点，各自对应的孔径值，获取所述第一点对应的进汞体积和；根据所有第一点各自对应的孔隙体积和，得到所述目标孔径值对应的进汞体积和。在一种可能的设计中，在所述根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像之前，还包括：对所述三维数字图像进行滤波，得到滤波后的三维数字图像；则所述根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像，包括：根据所述滤波后的三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像。第二方面，本申请提供一种基于模拟压汞法获取岩石的孔径分布的装置，包括：三维数字图像获取模块，用于获取待测定岩石的岩心的至少部分的三维数字图像；骨架图像获取模块，用于根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像；骨架图像获取模块，还用于获取所述骨架图像包括的各点对应的孔径值；点对应孔径值为所述点距离孔隙边界的最小距离的两倍孔径分布获取模块，用于根据所述骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布。在一种可能的设计中，所述孔径分布获取模块，具体用于：确定注汞端面上所有的第一点对应的孔径值中，最大的第一孔径值，以及所述各点对应的孔径值中，最小的第二孔径值；所述注汞端面为所述骨架图像的第一底面；从所述第一孔径值和所述第二孔径值之间的孔径值中确定出M个目标孔径值，所述M个目标孔径值包括所述第一孔径值和所述第二孔径值；对于每个目标孔径值，通过模拟注汞法，获取所述孔径值对应的进汞体积和，所述进汞体积和为相应孔径值对应的孔隙体积和；根据所述M个目标孔径值和每个目标孔径值对应的孔隙体积和，得到所述待测定岩石的孔隙的孔径分布。在一种可能的设计中，所述孔径分布获取模块，具体用于：对于所述注汞端面的每个第一点，获取包括所述第一点的至少一条第一骨架线；对于每条第一骨架线，确定每条第一骨架线包括的至少一个目标点，所述目标点为第一点和第二点之间的点，所述第二点为在第一方向上，所述第一骨架线上第一个对应的孔径值小于所述目标孔径值的点；所述第一方向为从所述注汞端面指向所述骨架图像的第二底面的方向，至少一个目标点包括第一点和所述第二点；根据所有第一骨架线包括的所有目标点，各自对应的孔径值，获取所述第一点对应的进汞体积和；根据所有第一点各自对应的孔隙体积和，得到所述目标孔径值对应的进汞体积和。在一种可能的设计中，所述装置还包括滤波模块，所述滤波模块用于在根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像之前，对所述三维数字图像进行滤波，得到滤波后的三维数字图像；则骨架图像获取模块，用于根据所述滤波后的三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像。本申请的基于模拟压汞法获取岩石的孔径分布的方法，通过获取待测定岩石的三维数字图像，对三维数字图像进行处理，得到待测定岩心对应的孔隙的骨架图像，根据骨架图像包括的各点各自对应的孔径值，得到待测定岩石的孔隙的孔径分布，可获取待测定岩石包括的所有孔径值对应的孔隙体积。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本申请提供的基于模拟压汞法获取岩石的孔径分布的方法的流程图一；图2为本申请实施例提供的岩心的至少部分的CT扫描图像的横截面图像；图3为本申请实施例提供的岩心的至少部分包括孔隙的图像的效果图；图4为本申请实施例提供的岩心的至少部分包括的孔隙的骨架线图像本文档来自技高网...

【技术保护点】
1.一种基于模拟压汞法获取岩石的孔径分布的方法，其特征在于，包括：获取待测定岩石的岩心的至少部分的三维数字图像；根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像；获取所述骨架图像包括的各点对应的孔径值；其中，点对应的孔径值为所述点距离孔隙边界的最小距离的两倍；根据所述骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布。

【技术特征摘要】
1.一种基于模拟压汞法获取岩石的孔径分布的方法，其特征在于，包括：获取待测定岩石的岩心的至少部分的三维数字图像；根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像；获取所述骨架图像包括的各点对应的孔径值；其中，点对应的孔径值为所述点距离孔隙边界的最小距离的两倍；根据所述骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布。2.根据权利要求1所述的方法，其特征在于，所根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像，包括：对所述三维数字图像进行图像分割，得到所述岩心的至少部分包括的孔隙的图像；根据所述孔隙的图像，获取所述骨架图像。3.根据权利要求2所述的方法，其特征在于，所述根据所述孔隙的图像，获取所述骨架图像，包括：对所述孔隙的图像进行距离变换，得到第一图像；对所述第一图像进行三维细化处理，得到所述骨架图像。4.根据权利要求1所述的方法，其特征在于，根据所述骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布，包括：确定注汞端面上所有的第一点对应的孔径值中，最大的第一孔径值，以及所述各点对应的孔径值中，最小的第二孔径值；所述注汞端面为所述骨架图像的第一底面；从所述第一孔径值和所述第二孔径值之间的孔径值中确定出M个目标孔径值，所述M个目标孔径值包括所述第一孔径值和所述第二孔径值；对于每个目标孔径值，通过模拟注汞法，获取所述孔径值对应的进汞体积和，所述进汞体积和为相应孔径值对应的孔隙体积和；根据所述M个目标孔径值和每个目标孔径值对应的孔隙体积和，得到所述待测定岩石的孔隙的孔径分布。5.根据权利要求4所述的方法，其特征在于，所述通过模拟注汞法，获取所述孔径值对应的进汞体积和，包括：对于所述注汞端面的每个第一点，获取包括所述第一点的至少一条第一骨架线；对于每条第一骨架线，确定每条第一骨架线包括的至少一个目标点，所述目标点为第一点和第二点之间的点，所述第二点为在第一方向上，所述第一骨架线上第一个对应的孔径值小于所述目标孔径值的点；所述第一方向为从所述注汞端面指向所述骨架图像的第二底面的方向，至少一个目标点包括第一点和所述第二点；根据所有第一骨架线包括的所有目标点，各自对应的孔径值，获取所述第一点对应的进汞体积和；根据所有第一点各自对应的孔隙体积和，得到所述目标孔径值对应的进汞体积和。6.根据权利要求1～5任一项所述的方法，其特征在于，在所述根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像之前...

【专利技术属性】
技术研发人员：李国梁，杨继进，张玉星，
申请(专利权)人：中国科学院地质与地球物理研究所，
类型：发明
国别省市：北京,11