CN110018108B - 岩石孔隙孔径的测定方法及设备 - Google Patents
岩石孔隙孔径的测定方法及设备 Download PDFInfo
- Publication number
- CN110018108B CN110018108B CN201910397472.8A CN201910397472A CN110018108B CN 110018108 B CN110018108 B CN 110018108B CN 201910397472 A CN201910397472 A CN 201910397472A CN 110018108 B CN110018108 B CN 110018108B
- Authority
- CN
- China
- Prior art keywords
- pore
- size
- rock
- pores
- volume
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 466
- 239000011435 rock Substances 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000004364 calculation method Methods 0.000 claims abstract description 33
- 238000000605 extraction Methods 0.000 claims description 10
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 9
- 239000011707 mineral Substances 0.000 claims description 9
- 230000011218 segmentation Effects 0.000 claims description 8
- 238000003709 image segmentation Methods 0.000 claims description 7
- 238000012216 screening Methods 0.000 claims description 6
- 238000013461 design Methods 0.000 description 21
- 238000010586 diagram Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 6
- 239000011800 void material Substances 0.000 description 6
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 239000011343 solid material Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供一种岩石孔隙孔径的测定方法及设备，该方法包括：获取待测定岩石对应的三维岩心图像，并提取三维岩心图像中的所有孔隙，获取预设的几何体的尺寸，分别根据几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目，分别根据各孔隙对应的最多填充数目和几何体的尺寸计算各孔隙对应的孔隙体积，分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值。本发明提供的岩石孔隙孔径的测定方法可以提高孔径计算的准确度。
Description
技术领域
本发明实施例涉及石油天然气地质勘探技术领域，尤其涉及一种岩石孔隙孔径的测定方法及设备。
背景技术
储集层岩石中未被矿物颗粒、胶结物或者其它固体物质填集的空间称为孔隙。由于孔隙的几何形态千差万别，导致储层储油和储气能力不同，油和气的赋存状态、运移机理和能力也不相同。因此，准确测定不同几何形态的孔隙的孔径大小至关重要。
目前，现有技术在测定不同集合形态的孔隙的孔径时，一般采用等效体积法进行测定，即将数字岩心中的孔隙直接按体积等效为球体，该球体的直径即代表该孔隙的孔径值。
然而，发明人发现现有技术中至少存在如下问题：对于几何形态不同孔隙，采用等效体积法计算出的孔隙体积会包含孔隙的狭长部分的体积，会导致计算出的孔径值偏大。
发明内容
本发明实施例提供一种岩石孔隙孔径的测定方法及设备，以解决现有技术中孔隙孔径值计算结果偏大的问题。
第一方面，本发明实施例提供一种岩石孔隙孔径的测定方法，包括：
获取待测定岩石对应的三维岩心图像，并提取所述三维岩心图像中的所有孔隙；
获取预设的几何体的尺寸；
分别根据所述几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目；
分别根据各孔隙对应的最多填充数目和所述几何体的尺寸计算各孔隙对应的孔隙体积；
分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值。
在一种可能的设计中，在所述提取所述三维岩心图像中的所有孔隙之前，还包括：
对所述三维岩心图像进行图像阈值分割得到所述三维岩心图像中的所有孔隙区域和矿物区域；
相应地，所述提取所述三维岩心图像中的所有孔隙，包括：
分别根据预设的分水岭算法对各孔隙区域进行分割，得到各孔隙区域中的所有孔隙。
在一种可能的设计中，在所述分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值之后，还包括：
根据各孔隙的孔径值得到所述待测定岩石中的孔隙的孔径分布。
在一种可能的设计中，所述根据各孔隙的孔径值得到所述待测定岩石中的孔隙的孔径分布，包括：
获取多个指定孔径值；
针对各个指定孔径值，统计孔径值为该指定孔径值的孔隙的个数；
在预设孔径数目坐标系中，对各个指定孔径值和各个指定孔径值对应的个数进行曲线拟合，生成孔径数目与孔径的关系曲线。
在一种可能的设计中，所述分别根据所述几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目，包括：
分别根据各孔隙包含的体素的坐标信息确定各孔隙的尺寸；
分别根据各孔隙的尺寸和所述几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目。
在一种可能的设计中，所述分别根据预设的分水岭算法对各孔隙区域进行分割，得到各孔隙区域中的所有孔隙之前，还包括：
获取各孔隙区域包含的体素数目；
筛选出体素数目大于预设体素数目的孔隙区域。
在一种可能的设计中，所述分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值，包括：
分别根据各孔隙对应的孔隙体积和预设的球体体积公式计算各孔隙对应的球体直径；
将各孔隙对应的球体直径作为对应的孔隙的孔径值。
第二方面，本发明实施例提供一种岩石孔隙孔径的测定设备，包括：
孔隙提取模块，用于获取待测定岩石对应的三维岩心图像，并提取所述三维岩心图像中的所有孔隙；
尺寸获取模块，用于获取预设的几何体的尺寸；
填充数目计算模块，用于分别根据所述几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目；
孔隙体积计算模块，用于分别根据各孔隙对应的最多填充数目和所述几何体的尺寸计算各孔隙对应的孔隙体积；
孔径值计算模块，用于分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值。
在一种可能的设计中，所述岩石孔隙孔径的测定设备还包括：图像分割模块；
所述第一图像分割模块，用于对所述三维岩心图像进行图像阈值分割得到所述三维岩心图像中的所有孔隙区域和矿物区域；
在一种可能的设计中，所述孔隙提取模块具体用于分别根据预设的分水岭算法对各孔隙区域进行分割，得到各孔隙区域中的所有孔隙。
在一种可能的设计中，所述岩石孔隙孔径的测定设备还包括：孔径分布生成模块；
所述孔径分布生成模块，用于根据各孔隙的孔径值得到所述待测定岩石中的孔隙的孔径分布。
在一种可能的设计中，所述孔径分布生成模块具体用于：
获取多个指定孔径值；
针对各个指定孔径值，统计孔径值为该指定孔径值的孔隙的个数；
在预设孔径数目坐标系中，对各个指定孔径值和各个指定孔径值对应的个数进行曲线拟合，生成孔径数目与孔径的关系曲线。
在一种可能的设计中，所述填充数目计算模块具体用于：分别根据各孔隙包含的体素的坐标信息确定各孔隙的尺寸；分别根据各孔隙的尺寸和所述几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目。
在一种可能的设计中，所述孔隙提取模块还用于：在分别根据预设的分水岭算法对各孔隙区域进行分割，得到各孔隙区域中的所有孔隙之前，获取各孔隙区域包含的体素数目；筛选出体素数目大于预设体素数目的孔隙区域。
在一种可能的设计中，所述孔径值计算模块具体用于：分别根据各孔隙对应的孔隙体积和预设的球体体积公式计算各孔隙对应的球体直径；将各孔隙对应的球体直径作为对应的孔隙的孔径值。
第三方面，本发明实施例提供一种岩石孔隙孔径的测定设备，包括：至少一个处理器和存储器；
所述存储器存储计算机执行指令；
所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述至少一个处理器执行如上第一方面以及第一方面各种可能的设计所述的岩石孔隙孔径的测定方法。
第四方面，本发明实施例提供一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机执行指令，当处理器执行所述计算机执行指令时，实现如上第一方面以及第一方面各种可能的岩石孔隙孔径的测定方法。
本实施例提供的岩石孔隙孔径的测定方法及设备，该方法通过提取待测定岩石对应的三维岩心图像中的所有孔隙，计算各孔隙最多可填充的预设尺寸的几何体的数目，得到各孔隙各自对应的最多填充数目，基于各孔隙各自对应的最多填充数目和该几何体的尺寸，得到各孔隙各自对应的孔隙体积，对各孔隙各自对应的孔隙体积进行等体积计算，得到各孔隙各自的孔径值，当孔隙中的部分区域较狭长，填充不了几何体时，计算出的孔隙体积不包含该狭长区域，因此利用最多填充数目计算出的孔径值更加准确，更能反映出待测定岩石的储存能力以及运移能力，以帮助相关研究人员制定出正确的勘探方案。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的孔隙孔径的计算方法的示意图；
图2为本发明实施例提供的岩石孔隙孔径的测定方法的流程示意图一；
图3为本发明实施例提供的滤波后的三维岩心图像的示意图；
图4为本发明实施例提供的岩石孔隙孔径的测定方法的流程示意图二；
图5为本发明实施例提供的孔径数目与孔径的关系曲线的示意图；
图6为本发明实施例提供的孔隙体积占比与孔径的关系曲线的示意图；
图7为本发明实施例提供的岩石孔隙孔径的测定设备的结构示意图一；
图8为本发明实施例提供的岩石孔隙孔径的测定设备的结构示意图二；
图9为本发明实施例提供的岩石孔隙孔径的测定设备的硬件结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明实施例提供的孔隙孔径的计算方法的示意图。如图1中图(a)所示，孔隙的几何形态不是规则的，当使用等效体积法测定数字岩心中的孔隙孔径时，会扩大狭长孔隙区域在孔径计算中的作用，致使计算出的孔径值偏大(参见图1中的图(b))，本实施例提供一种岩石孔隙孔径的测定方法，其是基于最大格子法(参见图1中的图(c))，即计算数字岩心中的孔隙空间最多可填充预设固定大小的格子的数目，将填充的所有格子的体积转化为等体积的球体，该球体直径代表该孔隙的孔径值，当使用最大格子法计算孔隙孔径时，由于狭长孔隙区域可能无法填充格子，使得计算出孔径值更加准确，提高孔隙孔径的计算的准确度。下面采用详细的实施例进行详细说明。
图2为本发明实施例提供的岩石孔隙孔径的测定方法的流程示意图一，如图2所示，该方法包括：
S201、获取待测定岩石对应的三维岩心图像，并提取三维岩心图像中的所有孔隙。
在本实施例中，待测定岩石可为储存石油或者天然气等资源的储集层岩石。岩心可为根据地质勘查工作或工程的需要，使用环状岩心钻头及其他取芯工具，从孔内取出的圆柱状岩石样品，即岩心是待测定岩石的一部分，其可为圆柱体。当需要测定岩心中的孔隙孔径时，利用岩心扫描设备对岩心进行扫描，得到三维岩心图像。
为了测定岩心中的孔隙，需要先从三维岩心图像中提取其包含的孔隙，在提取孔隙之前，需要先对三维岩心图像进行分割，确定该三维岩心图像中的孔隙区域和矿物区域，具体过程为：对三维岩心图像进行图像阈值分割得到三维岩心图像中的所有孔隙区域和矿物区域。
其中，阈值分割是一种常用的图像分割方法，采用阈值分割对图像进行分割时实现简单、计算量较小、性能较稳定，其实现原理是：通过设定不同的特征阈值，将图像的象素值分成若干类，得到多个像素集合，得到的每个集合形成一个与现实景物相对应的区域，各个区域内部具有一致的属性，而相邻区域不具有这种一致属性，从而确定图像中的目标区域和背景区域。
为了减少孔隙孔径的计算量，可以在从三维岩心图像中的孔隙区域提取孔隙前，将异常孔隙区域剔除，该异常孔隙区域表示体积较小的孔隙区域，具体过程可以为：获取各孔隙区域包含的体素数目。筛选出体素数目大于预设体素数目的孔隙区域。由于孔隙区域是由至少一个体素组成，当孔隙区域的体积较小时，其包含的体素也相应较少，因此，可以通过孔隙区域包含的体素的数目剔除体积较小的孔隙区域。
当确定出三维岩心图像中的孔隙区域和矿物区域后，可以根据分水岭算法提取孔隙区域中的所有孔隙，其具体过程为：分别根据预设的分水岭算法对各孔隙区域进行分割，得到各孔隙区域中的所有孔隙。
为了提高后续孔隙孔径计算的准确性，还可以在获得待测定岩石对应的三维岩心图像后，对该三维岩心图像进行图像滤波处理，得到滤波后的三维岩心图像(如图3所示)，在后续计算孔隙孔径过程中，利用滤波后的三维岩心图像。通过对图像进行滤波处理可以去除该三维岩心图像中的噪声点且可以保留该三维图像的细节特征。
S202、获取预设的几何体的尺寸。
在本实施例中，预设的几何体可以是长方体、也可以为正方体，或其它几何体。相关人员可以根据实际需要将该几何体的尺寸设定为该三维岩心图像中的体素的若干倍，例如，本实施例中将几何体的尺寸设定为体素的两倍，例如，体素的大小为3um×3um×3um，则几何体的尺寸为6um×6um×6um。
S203、分别根据几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目。
本实施例中基于最大格子法测定孔隙孔径，因此需要计算各个孔隙最多可填充几何体的数目，其计算过程可以为：分别根据各孔隙包含的体素的坐标信息确定各孔隙的尺寸。分别根据各孔隙的尺寸和几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目。
可以根据孔隙包含的体素的坐标信息以及几何体的尺寸计算该孔隙的孔隙体积，该坐标信息为体素在x方向上的第一坐标值、在y方向上的第二坐标值以及在z方向上的第三坐标值，以任一坐标轴为基准，例如，以z轴为基准，将孔隙划分为多个孔隙区域，每个孔隙区域的高度为该几何体的高度，即该几何体在z方向上的大小。针对各孔隙区域，在该孔隙区域中查找距离z轴最近的体素以及距离z轴最远的体素，获取最近的体素的坐标信息以及最远的体素的坐标信息，根据最远的体素的坐标信息和最近的坐标信息，得到该孔隙区域的尺寸，将该孔隙区域的尺寸除以该几何体的尺寸得到该孔隙区域对应的最多填充数目。将一个孔隙对应的所有孔隙区域所对应的最多填充数目进行相加，得到该孔隙对应的最多填充数目。
其中，由于该孔隙区域的高度已确定，相应地，根据最远的体素的坐标信息和最近的坐标信息，得到该孔隙区域的尺寸的过程可以为：将最远的坐标信息中的第一坐标值减去最近的坐标信息中的第一坐标值，得到该孔隙区域的长度，将最远的坐标信息中的第二坐标值减去最近的坐标信息中的第二坐标值，得到该孔隙区域的宽度，将该孔隙区域的长度×该孔隙区域的宽度×该孔隙区域的高度，得到该孔隙区域的尺寸。
为了减少计算量，在计算孔隙的最多可填充的几何体的数目时，也可以直接统计该孔隙包含的体素的个数，由于几何体的尺寸是体素尺寸的倍数，因此可以直接将体素的个数除以倍数得到该孔隙对应的最多可填充的几何体的数目。
S204、分别根据各孔隙对应的最多填充数目和几何体的尺寸计算各孔隙对应的孔隙体积。
在本实施例中，计算第一孔隙对应的最多填充数目与几何体的尺寸的乘积，得到第一孔隙对应的孔隙体积，其中该第一孔隙为多个孔隙中的任意一个孔隙。
S205、分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值。
可选地，步骤S205包括：分别根据各孔隙对应的孔隙体积和预设的球体体积公式计算各孔隙对应的球体直径。将各孔隙对应的球体直径作为对应的孔隙的孔径值。
在本实施例中，球体体积公式为V＝3/4*π*R3其中，V为球体体积，R为球体的半径，将第一孔隙的孔隙体积作为球体体积代入球体体积公式中，计算得到对应的球体半径，然后将球体半径乘以2得到球体直径，将该球体直径作为该第一孔隙的孔径值。
从上述描述可知，通过提取待测定岩石对应的三维岩心图像中的所有孔隙，计算各孔隙最多可填充的预设尺寸的几何体的数目，得到各孔隙各自对应的最多填充数目，基于各孔隙各自对应的最多填充数目和该几何体的尺寸，得到各孔隙各自对应的孔隙体积，对各孔隙各自对应的孔隙体积进行等体积计算，得到各孔隙各自的孔径值，当孔隙中的部分区域较狭长，填充不了几何体时，计算出的孔隙体积不包含该狭长区域，因此利用最多填充数目计算出的孔径值更加准确，更能反映出待测定岩石的储存能力以及运移能力，以帮助相关研究人员制定出正确的勘探方案。
图4为本发明实施例提供的岩石孔隙孔径的测定方法的流程示意图二，本实施例在图2实施例的基础上，对本实施例的具体实现过程进行了详细说明。如图4所示，该方法包括：
S401、获取待测定岩石对应的三维岩心图像，并提取三维岩心图像中的所有孔隙。
S402、获取预设的几何体的尺寸。
S403、分别根据几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目。
S404、分别根据各孔隙对应的最多填充数目和几何体的尺寸计算各孔隙对应的孔隙体积。
S405、分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值。
在本实施例中，S401至S405与图2实施例中的S201至S205类似，本实施例此处不再进行赘述。
S406、根据各孔隙的孔径值得到待测定岩石中的孔隙的孔径分布。
在本实施例中，可以基于各孔隙的孔径值得到孔隙的孔径分布规律，该孔径分布规律可以为孔径数目与孔径值的分布规律，其过程可以为：获取多个指定孔径值，针对各个指定孔径值，统计孔径值为该指定孔径值的孔隙的个数。在预设孔径数目坐标系中，对各个指定孔径值和各个指定孔径值对应的个数进行曲线拟合，生成孔径数目与孔径的关系曲线。
在本实施例中，相关人员可以按照实际需求确定多个指定孔径值，针对各个指定孔径值，统计孔径值为该指定孔径值的孔隙的个数。分别将各个指定孔径值以及各个指定孔径值各自对应的个数作为一个坐标点，并标注在横轴为孔径值，纵轴为孔径数目的预设孔径数目坐标系中，对标注后的预设孔径数目坐标系中的各个坐标点进行曲线拟合，得到孔径数目与孔径的关系曲线。
为了提高数据标注的准确度以及更好地表示孔径与数目的分布规律，可以在得到各个指定孔径值各自对应的孔隙的个数后，对个数取对数，得到各个指定孔径值各自对应的个数的对数值，将各个指定孔径值以及各个指定孔径值各自对应的对数值作为一个坐标点，并标注在横轴为孔径值，纵轴为孔径对数目的预设孔径数目坐标系中，对标注后的预设孔径数目坐标系中的各个坐标点进行曲线拟合，得到孔径数目与孔径的关系曲线(如图5所示)。
如图5所示，孔隙数目随孔径的增大呈现递减趋势，大多数孔隙的孔径分布在40um以下。
在本实施例中，根据各孔隙的孔径值得到待测定岩石中的孔隙的孔径分布还可以包括：统计提取出的所有孔隙的孔隙总数目，针对各个指定孔径值，获取该指定孔径值对应的孔隙体积，并统计体积为该孔隙体积的第一孔隙数目，计算该第一孔隙数目与孔隙总数目的比值，得到该指定孔径值对应的体积占比，在预设孔径体积坐标系中，对各个指定孔径值和各个指定孔径值对应的体积占比进行曲线拟合，生成孔隙体积占比与孔径的关系曲线。
在本实施例中，针对各个指定孔径值，获取所有孔隙的孔隙总数目，针对各个指定孔径值，获取该指定孔径值对应的孔隙体积，统计体积为该孔隙体积的第一孔隙数目，计算该第一孔隙数目与孔隙总数目的比值，得到该指定孔径值对应的体积占比。分别将各个指定孔径值以及各个指定孔径值各自对应的体积占比作为一个坐标点，并标注在横轴为孔径值，纵轴为体积的预设孔径体积坐标系中，对标注后的预设孔径体积坐标系中的各个坐标点进行曲线拟合，得到孔隙体积占比与孔径的关系曲线(如图6所示)。
其中，体积占比为体积所占比例。
如图6所示，不同孔径对应的孔隙体积占比呈现两头小中间大的趋势，大多数的孔隙体积位于孔径90-160um范围内。
在本实施例中，根据计算出的孔隙体积以及孔隙的孔径还可以发现，最大的孔隙体积以及该最大的孔隙体积对应的最大孔径值，最小的孔隙体积以及该最小的孔隙体积对应的最小孔径值。
在本实施例中，通过根据计算出的孔隙的孔径值得到待测定岩石的孔径分布特征，研究人员可以通过研究该待测定岩石的孔径分布特征得到待测定岩石的渗流能力，制定出正确的勘探开发方案。
图7为本发明实施例提供的岩石孔隙孔径的测定设备的结构示意图一。如图7所示，该岩石孔隙孔径的测定设备70包括：孔隙提取模块701、尺寸获取模块702、填充数目计算模块703、孔隙体积计算模块704以及孔径值计算模块705。
孔隙提取模块701，用于获取待测定岩石对应的三维岩心图像，并提取三维岩心图像中的所有孔隙。
尺寸获取模块702，用于获取预设的几何体的尺寸。
填充数目计算模块703，用于分别根据几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目。
孔隙体积计算模块704，用于分别根据各孔隙对应的最多填充数目和几何体的尺寸计算各孔隙对应的孔隙体积。
孔径值计算模块705，用于分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值。
在一种可能的设计中，岩石孔隙孔径的测定设备还包括：图像分割模块。
第一图像分割模块，用于对三维岩心图像进行图像阈值分割得到三维岩心图像中的所有孔隙区域和矿物区域。
在一种可能的设计中，孔隙提取模块701具体用于分别根据预设的分水岭算法对各孔隙区域进行分割，得到各孔隙区域中的所有孔隙。
在一种可能的设计中，填充数目计算模块703具体用于：分别根据各孔隙包含的体素的坐标信息确定各孔隙的尺寸。分别根据各孔隙的尺寸和几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目。
在一种可能的设计中，孔隙提取模块704还用于：在分别根据预设的分水岭算法对各孔隙区域进行分割，得到各孔隙区域中的所有孔隙之前，获取各孔隙区域包含的体素数目。筛选出体素数目大于预设体素数目的孔隙区域。
在一种可能的设计中，孔径值计算模块705具体用于：分别根据各孔隙对应的孔隙体积和预设的球体体积公式计算各孔隙对应的球体直径。将各孔隙对应的球体直径作为对应的孔隙的孔径值。
本实施例提供的设备，可用于执行上述方法实施例的技术方案，其实现原理和技术效果类似，本实施例此处不再赘述。
图8为本发明实施例提供的岩石孔隙孔径的测定设备的结构示意图二。如图8所示，本实施例在图7实施例的基础上，还包括：孔径分布生成模块706。
在一种可能的设计中，孔径分布生成模块706，用于根据各孔隙的孔径值得到待测定岩石中的孔隙的孔径分布。
在一种可能的设计中，孔径分布生成模块706具体用于：获取多个指定孔径值。针对各个指定孔径值，统计孔径值为该指定孔径值的孔隙的个数。在预设孔径数目坐标系中，对各个指定孔径值和各个指定孔径值对应的个数进行曲线拟合，生成孔径数目与孔径的关系曲线。
本实施例提供的设备，可用于执行上述方法实施例的技术方案，其实现原理和技术效果类似，本实施例此处不再赘述。
图9为本发明实施例提供的岩石孔隙孔径的测定设备的硬件结构示意图。如图9所示，本实施例的岩石孔隙孔径的测定设备90包括：处理器901以及存储器902；其中存储器902，用于存储计算机执行指令；
处理器901，用于执行存储器存储的计算机执行指令，以实现上述实施例中接收设备所执行的各个步骤。具体可以参见前述方法实施例中的相关描述。
可选地，存储器902既可以是独立的，也可以跟处理器901集成在一起。
当存储器902独立设置时，该语音交互设备还包括总线903，用于连接所述存储器902和处理器901。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，当处理器执行所述计算机执行指令时，实现如上所述的岩石孔隙孔径的测定方法。
在本发明所提供的几个实施例中，应该理解到，所揭露的设备和方法，可以通过其它的方式实现。例如，以上所描述的设备实施例仅仅是示意性的，例如，所述模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或模块的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的模块可以是或者也可以不是物理上分开的，作为模块显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能模块可以集成在一个处理单元中，也可以是各个模块单独物理存在，也可以两个或两个以上模块集成在一个单元中。上述模块成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
上述以软件功能模块的形式实现的集成的模块，可以存储在一个计算机可读取存储介质中。上述软件功能模块存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)或处理器(英文：processor)执行本申请各个实施例所述方法的部分步骤。
应理解，上述处理器可以是中央处理单元(英文：Central Processing Unit，简称：CPU)，还可以是其他通用处理器、数字信号处理器(英文：Digital Signal Processor，简称：DSP)、专用集成电路(英文：Application Specific Integrated Circuit，简称：ASIC)等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合发明所公开的方法的步骤可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。
存储器可能包含高速RAM存储器，也可能还包括非易失性存储NVM，例如至少一个磁盘存储器，还可以为U盘、移动硬盘、只读存储器、磁盘或光盘等。
总线可以是工业标准体系结构(Industry Standard Architecture，ISA)总线、外部设备互连(Peripheral Component，PCI)总线或扩展工业标准体系结构(ExtendedIndustry Standard Architecture，EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，本申请附图中的总线并不限定仅有一根总线或一种类型的总线。
上述存储介质可以是由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(SRAM)，电可擦除可编程只读存储器(EEPROM)，可擦除可编程只读存储器(EPROM)，可编程只读存储器(PROM)，只读存储器(ROM)，磁存储器，快闪存储器，磁盘或光盘。存储介质可以是通用或专用计算机能够存取的任何可用介质。
一种示例性的存储介质耦合至处理器，从而使处理器能够从该存储介质读取信息，且可向该存储介质写入信息。当然，存储介质也可以是处理器的组成部分。处理器和存储介质可以位于专用集成电路(Application Specific Integrated Circuits，简称：ASIC)中。当然，处理器和存储介质也可以作为分立组件存在于电子设备或主控设备中。
本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (8)
1.一种岩石孔隙孔径的测定方法，其特征在于，包括：
获取待测定岩石对应的三维岩心图像，并提取所述三维岩心图像中的所有孔隙；
获取预设的几何体的尺寸；
分别根据所述几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目；
分别根据各孔隙对应的最多填充数目和所述几何体的尺寸计算各孔隙对应的孔隙体积；
分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值；
其中，在所述提取所述三维岩心图像中的所有孔隙之前，还包括：
对所述三维岩心图像进行图像阈值分割得到所述三维岩心图像中的所有孔隙区域和矿物区域；其中所述孔隙区域是由至少一个体素组成；
其中，所述提取所述三维岩心图像中的所有孔隙，包括：
分别根据预设的分水岭算法对各孔隙区域进行分割，得到各孔隙区域中的所有孔隙；
其中所述几何体的尺寸是所述体素的尺寸的至少两倍；
其中，所述分别根据预设的分水岭算法对各孔隙区域进行分割，得到各孔隙区域中的所有孔隙之前，还包括：
获取各孔隙区域包含的体素数目；
筛选出体素数目大于预设体素数目的孔隙区域。
2.根据权利要求1所述的方法，其特征在于，在所述分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值之后，还包括：
根据各孔隙的孔径值得到所述待测定岩石中的孔隙的孔径分布。
3.根据权利要求2所述的方法，其特征在于，所述根据各孔隙的孔径值得到所述待测定岩石中的孔隙的孔径分布，包括：
获取多个指定孔径值；
针对各个指定孔径值，统计孔径值为该指定孔径值的孔隙的个数；
在预设孔径数目坐标系中，对各个指定孔径值和各个指定孔径值对应的个数进行曲线拟合，生成孔径数目与孔径的关系曲线。
4.根据权利要求1所述的方法，其特征在于，所述分别根据所述几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目，包括：
分别根据各孔隙包含的体素的坐标信息确定各孔隙的尺寸；
分别根据各孔隙的尺寸和所述几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目。
5.根据权利要求1-4任一项所述的方法，其特征在于，所述分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值，包括：
分别根据各孔隙对应的孔隙体积和预设的球体体积公式计算各孔隙对应的球体直径；
将各孔隙对应的球体直径作为对应的孔隙的孔径值。
6.一种岩石孔隙孔径的测定设备，其特征在于，包括：
孔隙提取模块，用于获取待测定岩石对应的三维岩心图像，并提取所述三维岩心图像中的所有孔隙；
尺寸获取模块，用于获取预设的几何体的尺寸；
填充数目计算模块，用于分别根据所述几何体的尺寸计算各孔隙最多可填充的几何体的数目，得到各孔隙对应的最多填充数目；
孔隙体积计算模块，用于分别根据各孔隙对应的最多填充数目和所述几何体的尺寸计算各孔隙对应的孔隙体积；
孔径值计算模块，用于分别对各孔隙对应的孔隙体积进行等体积计算，得到各孔隙的孔径值；
其中，所述岩石孔隙孔径的测定设备还包括：图像分割模块；
所述图像分割模块，用于对所述三维岩心图像进行图像阈值分割得到所述三维岩心图像中的所有孔隙区域和矿物区域；其中所述孔隙区域是由至少一个体素组成；
其中，所述孔隙提取模块具体用于分别根据预设的分水岭算法对各孔隙区域进行分割，得到各孔隙区域中的所有孔隙；
其中所述几何体的尺寸是所述体素的尺寸的至少两倍；
其中，所述孔隙提取模块还用于：在分别根据预设的分水岭算法对各孔隙区域进行分割，得到各孔隙区域中的所有孔隙之前，获取各孔隙区域包含的体素数目；筛选出体素数目大于预设体素数目的孔隙区域。
7.一种岩石孔隙孔径的测定设备，其特征在于，包括：至少一个处理器和存储器；
所述存储器存储计算机执行指令；
所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述至少一个处理器执行如权利要求1至5任一项所述的岩石孔隙孔径的测定方法。
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机执行指令，当处理器执行所述计算机执行指令时，实现如权利要求1至5任一项所述的岩石孔隙孔径的测定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910397472.8A CN110018108B (zh) | 2019-05-14 | 2019-05-14 | 岩石孔隙孔径的测定方法及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910397472.8A CN110018108B (zh) | 2019-05-14 | 2019-05-14 | 岩石孔隙孔径的测定方法及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110018108A CN110018108A (zh) | 2019-07-16 |
| CN110018108B true CN110018108B (zh) | 2020-04-28 |
Family
ID=67193635
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910397472.8A Active CN110018108B (zh) | 2019-05-14 | 2019-05-14 | 岩石孔隙孔径的测定方法及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110018108B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112233166A (zh) * | 2020-09-11 | 2021-01-15 | 安徽理工大学 | 一种基于多孔介质三维孔隙空间图像的孔径分布评估方法 |
| CN112146957B (zh) * | 2020-09-25 | 2022-04-15 | 东北石油大学 | 一种基于数字岩心实现人造岩心定量制作方法 |
| CN113177949B (zh) * | 2021-04-16 | 2023-09-01 | 中南大学 | 一种大尺寸岩石颗粒特征识别方法及装置 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102222359B (zh) * | 2011-05-24 | 2012-09-26 | 中国石油天然气股份有限公司 | 一种岩心三维孔隙结构重构方法 |
| BR112016002538B8 (pt) * | 2013-08-08 | 2022-08-02 | Ingrain Inc | Método para gerar propriedades de formação de um ambiente de fundo de poço pressurizado e sistema para gerar propriedades de formação de um ambiente de fundo de poço pressurizado |
| CN104794709B (zh) * | 2015-04-10 | 2018-03-09 | 四川大学 | 一种三维岩心图像孔隙和喉道的分割方法 |
| CN105352873B (zh) * | 2015-11-26 | 2018-12-21 | 中国石油大学(北京) | 页岩孔隙结构的表征方法 |
| CN106324002A (zh) * | 2016-10-18 | 2017-01-11 | 中国石油大学(华东) | 一种基于岩石分类和多尺度数字岩心的碳酸盐岩孔隙结构表征方法 |
| CN106803253B (zh) * | 2017-01-17 | 2019-06-04 | 四川大学 | 一种三维岩石图像裂缝识别方法 |
| CN107067379B (zh) * | 2017-03-16 | 2018-07-20 | 中国科学院地质与地球物理研究所 | 基于三维fib-sem图像的页岩孔隙定量表征方法 |
| CN108491677A (zh) * | 2018-07-04 | 2018-09-04 | 河海大学 | 基于改进最大球法的微观孔隙模型的孔隙特征统计方法 |
| CN109142875B (zh) * | 2018-09-30 | 2021-08-10 | 西南石油大学 | 一种利用数字岩心获取致密砂岩岩石电学特性的方法 |
| CN109242985B (zh) * | 2018-10-29 | 2020-06-05 | 中国科学院力学研究所 | 一种从三维图像确定孔隙结构关键参数的方法 |
-
2019
- 2019-05-14 CN CN201910397472.8A patent/CN110018108B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110018108A (zh) | 2019-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110018108B (zh) | 岩石孔隙孔径的测定方法及设备 | |
| Akinci et al. | Parallel surface reconstruction for particle‐based fluids | |
| US11590708B2 (en) | Three-dimensional fluid micromodels | |
| CN104619952A (zh) | 具可靠多相渗透性决定数字岩石分析系统及方法 | |
| CN113190515A (zh) | 基于异构并行计算的城市级海量点云坐标转换方法 | |
| CN108663287B (zh) | 一种利用ct图像精确计算煤岩密度的方法 | |
| CN104573339B (zh) | 页岩气储层的地质参数确定方法和装置 | |
| CN116630565B (zh) | 基于多段线缓冲的地质钻孔三维模型生成方法及装置 | |
| CN115880436B (zh) | 一种cad模型可见性确定方法、系统及电子设备 | |
| CN111707695A (zh) | 基于核磁共振数据确定储层渗透率的方法、装置及设备 | |
| Schütz et al. | GPU‐Accelerated LOD Generation for Point Clouds | |
| CN110135382A (zh) | 一种人体检测方法和装置 | |
| CN115237355A (zh) | 一种基于闪存数据类识别的定向交换方法及系统 | |
| WO2021231170A1 (en) | Systems and methods for generating vertical and lateral heterogeneity indices of reservoirs | |
| CN115081210A (zh) | 一种受限制游走模拟以计算多孔介质内流动迂曲度的方法 | |
| CN112184900B (zh) | 高程数据的确定方法、装置及存储介质 | |
| Silvennoinen et al. | Occluder simplification using planar sections | |
| CN109993786B (zh) | 迂曲度获取方法、装置、设备以及存储介质 | |
| CN114187423A (zh) | 用于三维模拟试验中围岩裂隙重构方法、电子设备及存储介质 | |
| CN115410049B (zh) | 一种岩体溶蚀程度的分类评估方法及装置 | |
| JP7333956B2 (ja) | 三次元の片刃分布の推定プログラム、方法及び装置 | |
| CN112885402B (zh) | 固态硬盘坏块查询方法、装置及计算机可读存储介质 | |
| CN117807154B (zh) | 一种用于展示系统的时序数据可视化方法、设备和介质 | |
| CN114460642A (zh) | 走滑断层分段识别方法及设备 | |
| CN111487685B (zh) | 岩石储集层的电阻率计算方法、装置及存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |