CN108492331B - 基于模拟压汞法获取岩石的孔径分布的方法和装置 - Google Patents
基于模拟压汞法获取岩石的孔径分布的方法和装置 Download PDFInfo
- Publication number
- CN108492331B CN108492331B CN201810227776.5A CN201810227776A CN108492331B CN 108492331 B CN108492331 B CN 108492331B CN 201810227776 A CN201810227776 A CN 201810227776A CN 108492331 B CN108492331 B CN 108492331B
- Authority
- CN
- China
- Prior art keywords
- point
- image
- skeleton
- aperture value
- target
- 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
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 title claims abstract description 93
- 229910052753 mercury Inorganic materials 0.000 title claims abstract description 93
- 239000011435 rock Substances 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 79
- 238000009826 distribution Methods 0.000 title claims abstract description 62
- 238000002347 injection Methods 0.000 title claims abstract description 43
- 239000007924 injection Substances 0.000 title claims abstract description 43
- 238000004088 simulation Methods 0.000 title abstract 2
- 239000011148 porous material Substances 0.000 claims abstract description 228
- 238000001914 filtration Methods 0.000 claims description 19
- 230000009466 transformation Effects 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 6
- 238000003709 image segmentation Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 abstract description 5
- 238000013461 design Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 11
- 239000000243 solution Substances 0.000 description 10
- 238000002591 computed tomography Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 9
- 238000002459 porosimetry Methods 0.000 description 7
- 239000007789 gas Substances 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000003825 pressing Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000001179 sorption measurement Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Generation (AREA)
- Image Processing (AREA)
- Processing Or Creating Images (AREA)
Abstract
本申请实施例提供一种基于模拟压汞法获取岩石的孔径分布的方法和装置，该方法包括：获取待测定岩石的岩心的至少部分的三维数字图像，根据该三维数字图像，获取岩心的至少部分包括的孔隙的骨架图像，获取该骨架图像包括的所有点各自对应的孔径值，点对应的孔径值为点距离孔隙边界的最小距离的两倍，根据骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布。本申请的方法，可测定的孔隙尺寸的范围不受限制，可获取各种待测定岩石的孔隙的孔径分布。
Description
技术领域
本申请实施例涉及石油天然气地质勘探技术，尤其涉及一种基于模拟压汞法获取岩石的孔径分布的方法和装置。
背景技术
储集层岩石的微观孔隙结构特征(孔径的大小、分布、孔隙的几何形状及相互联通情况)直接影响储集层的油气资源的聚集、流动和产出，是评价储集层资源潜力的关键参数之一。储集层岩石的孔径分布对油气的勘探开发至关重要，直接关系到储集层岩石的储油和储气能力、油和气的赋存状态、运移机理和能力以及储集层岩石的沉积演化过程，最终影响油气开采的难易程度和勘探开发价值。因此精确厘定储集层岩石的孔径分布特征意义重大。只有准确了解了储集层岩石的孔径分布特征才能正确反应储集层岩石的渗流能力，并制定出正确的勘探开发方案。
目前用来表征储集层孔隙特征的方法主要有：吸附法、压汞法、扫描电镜、光学显微镜、CT扫描、薄片分析等方法。不同的表征方法可监测的孔隙范围大不相同，其中，氮气吸附法可以测定小于100nm孔径的孔隙，压汞法可以测定的孔径范围为3nm-100um，上述两种方法对可测定的孔隙的孔径具有一定的限制。扫描电镜和CT扫描等成像方法虽然可以直接观察孔隙结构，但无法获知孔隙的尺寸特征。
发明内容
本申请实施例提供一种基于模拟压汞法获取岩石的孔径分布的方法和装置，可测定的孔隙尺寸的范围不受限制，可获取各种待测定岩石的孔隙的孔径分布。
第一方面，本申请提供一种基于模拟压汞法获取岩石的孔径分布的方法，其特征在于，包括：
获取待测定岩石的岩心的至少部分的三维数字图像；
根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像；
获取所述骨架图像包括的各点对应的孔径值；其中，点对应的孔径值为所述点距离孔隙边界的最小距离的两倍；
根据所述骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布。
在一种可能的设计中，所根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像，包括：
对所述三维数字图像进行图像分割，得到所述岩心的至少部分包括的孔隙的图像；
根据所述孔隙的图像，获取所述骨架图像。
在一种可能的设计中，所述根据所述孔隙的图像，获取所述骨架图像，包括：
对所述孔隙的图像进行距离变换，得到第一图像；
对所述第一图像进行三维细化处理，得到所述骨架图像。
在一种可能的设计中，根据所述骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布，包括：
确定注汞端面上所有的第一点对应的孔径值中，最大的第一孔径值，以及所述各点对应的孔径值中，最小的第二孔径值；所述注汞端面为所述骨架图像的第一底面；
从所述第一孔径值和所述第二孔径值之间的孔径值中确定出M个目标孔径值，所述M个目标孔径值包括所述第一孔径值和所述第二孔径值；
对于每个目标孔径值，通过模拟注汞法，获取所述孔径值对应的进汞体积和，所述进汞体积和为相应孔径值对应的孔隙体积和；
根据所述M个目标孔径值和每个目标孔径值对应的孔隙体积和，得到所述待测定岩石的孔隙的孔径分布。
在一种可能的设计中，所述通过模拟注汞法，获取所述目标孔径值对应的进汞体积和，包括：
对于所述注汞端面的每个第一点，获取包括所述第一点的至少一条第一骨架线；
对于每条第一骨架线，确定每条第一骨架线包括的至少一个目标点，所述目标点为第一点和第二点之间的点，所述第二点为在第一方向上，所述第一骨架线上第一个对应的孔径值小于所述目标孔径值的点；所述第一方向为从所述注汞端面指向所述骨架图像的第二底面的方向，至少一个目标点包括第一点和所述第二点；
根据所有第一骨架线包括的所有目标点，各自对应的孔径值，获取所述第一点对应的进汞体积和；
根据所有第一点各自对应的孔隙体积和，得到所述目标孔径值对应的进汞体积和。
在一种可能的设计中，在所述根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像之前，还包括：
对所述三维数字图像进行滤波，得到滤波后的三维数字图像；
则所述根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像，包括：
根据所述滤波后的三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像。
第二方面，本申请提供一种基于模拟压汞法获取岩石的孔径分布的装置，包括：
三维数字图像获取模块，用于获取待测定岩石的岩心的至少部分的三维数字图像；
骨架图像获取模块，用于根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像；
骨架图像获取模块，还用于获取所述骨架图像包括的各点对应的孔径值；点对应孔径值为所述点距离孔隙边界的最小距离的两倍
孔径分布获取模块，用于根据所述骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布。
在一种可能的设计中，所述孔径分布获取模块，具体用于：
确定注汞端面上所有的第一点对应的孔径值中，最大的第一孔径值，以及所述各点对应的孔径值中，最小的第二孔径值；所述注汞端面为所述骨架图像的第一底面；
从所述第一孔径值和所述第二孔径值之间的孔径值中确定出M个目标孔径值，所述M个目标孔径值包括所述第一孔径值和所述第二孔径值；
对于每个目标孔径值，通过模拟注汞法，获取所述孔径值对应的进汞体积和，所述进汞体积和为相应孔径值对应的孔隙体积和；
根据所述M个目标孔径值和每个目标孔径值对应的孔隙体积和，得到所述待测定岩石的孔隙的孔径分布。
在一种可能的设计中，所述孔径分布获取模块，具体用于：
对于所述注汞端面的每个第一点，获取包括所述第一点的至少一条第一骨架线；
对于每条第一骨架线，确定每条第一骨架线包括的至少一个目标点，所述目标点为第一点和第二点之间的点，所述第二点为在第一方向上，所述第一骨架线上第一个对应的孔径值小于所述目标孔径值的点；所述第一方向为从所述注汞端面指向所述骨架图像的第二底面的方向，至少一个目标点包括第一点和所述第二点；
根据所有第一骨架线包括的所有目标点，各自对应的孔径值，获取所述第一点对应的进汞体积和；
根据所有第一点各自对应的孔隙体积和，得到所述目标孔径值对应的进汞体积和。
在一种可能的设计中，所述装置还包括滤波模块，所述滤波模块用于在根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像之前，对所述三维数字图像进行滤波，得到滤波后的三维数字图像；
则骨架图像获取模块，用于根据所述滤波后的三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像。
本申请的基于模拟压汞法获取岩石的孔径分布的方法，通过获取待测定岩石的三维数字图像，对三维数字图像进行处理，得到待测定岩心对应的孔隙的骨架图像，根据骨架图像包括的各点各自对应的孔径值，得到待测定岩石的孔隙的孔径分布，可获取待测定岩石包括的所有孔径值对应的孔隙体积。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本申请提供的基于模拟压汞法获取岩石的孔径分布的方法的流程图一；
图2为本申请实施例提供的岩心的至少部分的CT扫描图像的横截面图像；
图3为本申请实施例提供的岩心的至少部分包括孔隙的图像的效果图；
图4为本申请实施例提供的岩心的至少部分包括的孔隙的骨架线图像的效果图；
图5为本申请实施例提供的骨架图像中包括的各点的示意图；
图6为本申请实施例提供的包括第一点的两条第一骨架线的示意图一；
图7为本申请实施例提供的包括第一点的两条第一骨架线的示意图二；
图8为本申请实施例提供的孔径-孔隙体积和曲线；
图9为本申请实施例提供的孔径-孔隙体积曲线；
图10为本申请实施例提供的孔径分布柱状图；
图11为本申请实施例提供的基于模拟压汞法获取岩石的孔径分布的方法的流程图二；
图12为对图2中的图像滤波后的图像；
图13为岩心的三维数字图像滤波后的立体图像效果图；
图14为本申请实施例提供的基于模拟压汞法获取岩石的孔径分布的装置的结构示意图一；
图15为本申请实施例提供的基于模拟压汞法获取岩石的孔径分布的装置的结构示意图二。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
图1为本申请提供的基于模拟压汞法获取岩石的孔径分布的方法的流程图一，如图1所示，本实施例的方法可以包括：
步骤S101、获取待测定岩石的岩心的至少部分的三维数字图像；
步骤S102、根据该三维数字图像，获取岩心的至少部分包括的孔隙的骨架图像；
步骤S103、获取该骨架图像包括的各点各自对应的孔径值；其中，点对应的孔径值为所述点距离孔隙边界的最小距离的两倍；
步骤S104、根据骨架图像包括的各点对应的孔径值，通过模拟注汞法获取待测定岩石的孔隙的孔径分布。
具体地，本实施例中的岩心可为根据地质勘查工作或工程的需要，使用环状岩心钻头及其他取心工具，从孔内取出的圆柱状岩石样品。也就是说，岩心是待测定岩石的一部分，可为圆柱体。其中，待测定岩石可为储存石油或者天然气等资源的储集层岩石。
对于步骤S101、采用扫描设备对待测定岩石的岩心的至少部分进行扫描，得到该岩心的至少部分三维数字图像，也就是说岩心的至少部分为成像区域，成像区域所成的像为本实施例中的三维数字图像。其中，岩心的至少部分可为整个岩心，也可为岩心的一部分，岩心的至少部分也为圆柱体。
扫描设备可为电子计算机断层扫描(Computed Tomography，简称CT)设备，还可为聚焦粒子束-扫描电子显微镜(scanning electron microscope，简称SEM)。其中，得到的图像大小不限制，每个像素大小也不限制。
图2为本申请实施例提供的岩心的至少部分的CT扫描图像的横截面图像。
若扫描设备可为CT设备，则获取的岩心的至少部分的三维数字图像可为CT图像，获取的图像大小可为：550×550×800，若每个像素代表3um，此时对应的扫描区域的长为1.65mm，宽为1.65mm，高为2.4mm。岩心的至少部分的CT图像如图2所示，其中，图2中的21为孔隙，22为矿物，也就是岩石的本体。
对于步骤S102和步骤S103、根据岩心的至少部分的三维数字图像，获取岩心的至少部分包括的孔隙的骨架图像；获取该骨架图像包括的各点各自对应的孔径值；其中，点对应的孔径值为点距离孔隙边界的最小距离的两倍，孔隙边界为相应点所在的孔隙的边界。
具体地，图3为本申请实施例提供的岩心的至少部分包括孔隙的图像的效果图；图4为本申请实施例提供的岩心的至少部分包括的孔隙的骨架线图像的效果图；图5为本申请实施例提供的骨架图像中包括的各点的示意图。
步骤“根据岩心的至少部分的三维数字图像，获取岩心的至少部分包括的孔隙的骨架图像”，包括：对三维数字图像进行图像分割，得到岩心的至少部分包括的孔隙的图像；根据岩心的至少部分包括的孔隙的图像，获取岩心至少部分包括的孔隙的骨架图像。
其中，对上述三维数字图像进行图像分割的方法，可为阈值分割的方法，即采用阈值分割的方法，将三维数字图像中的孔隙和矿物区分开来，得到岩心的至少部分包括的孔隙的图像，或者说岩心的成像区域包括的孔隙的图像，如图3所示。
进一步地，“根据岩心的至少部分包括的孔隙的图像，获取岩心至少部分包括的孔隙的骨架图像”包括：“对该孔隙的图像进行距离变换，得到第一图像；对该第一图像进行三维细化处理，得到岩心至少部分包括的孔隙的骨架图像”。
其中，岩心至少部分包括的孔隙的骨架图像可如图4所示，骨架图像中包括的各点如图5所示，参见图4和图5，该骨架图像由多条骨架线组成，组成孔隙的骨架图像的每个点的尺寸为一个体素大小，也就是骨架线上的每个点的尺寸为一个体素大小。其中，骨架图像的端点和组成所述骨架图像的骨架线之间的交点可称为骨架图像包括的节点，“骨架图像的端点”是指骨架线的端点，包括骨架图像的两个底面上的点。
距离变换可为欧拉变换，具体变换算法可采用现有的算法，三维细化算法可采用现有的算法，本实施例中不再赘述。
在“对该孔隙的图像进行距离变换，得到第一图像；对该第一图像进行三维细化处理，得到岩心至少部分包括的孔隙的骨架图像”的同时，便可获取到该骨架图像包括的所有点各自对应的孔径值，获取方法为现有技术中的方法，本实施例中不再赘述。
由于孔隙的骨架图像是由多条骨架线组成，而骨架线是孔隙所覆盖区域的中轴线，孔径值为所述骨架线上点距离孔隙边界的最小距离的两倍。
对于步骤S104、根据骨架图像包括的各点各自对应的孔径值，获取孔径分布，包括：
A1、确定注汞端面上所有的第一点对应的孔径值中，最大的第一孔径值，以及骨架图像包括的各点对应的孔径值中，最小的第二孔径值；注汞端面为骨架图像的第一底面；
A2、从第一孔径值和第二孔径值之间的孔径值中确定出M个目标孔径值，M个目标孔径值包括第一孔径值和第二孔径值；
A3、对于每个目标孔径值，通过模拟注汞过程，获取目标孔径值对应的进汞体积和，进汞体积和为相应目标孔径值对应的孔隙体积和；
A4、根据M个目标孔径值和每个目标孔径值对应的孔隙体积和，得到孔径分布。
具体地，对于A1，将骨架图像的一个底面作为模拟注汞过程的注汞端面，比如第一底面。其中，岩心为圆柱体，岩心的至少部分的三维数字图像也为圆柱体，因此，骨架图像也为圆柱体。
注汞端面上所有的第一点为骨架图像包括的端点中位于该注汞端面上的点，根据各第一点各自对应的孔径值，确定出最大的孔径值，此处称为第一孔径值；根据骨架图像包括的所有点对应的孔径值，确定出最小的孔径值，此处称为第二孔径值。
对于A2，从第一孔径值和第二孔径值之间的孔径值中确定出M个目标孔径值，M个目标孔径值包括第一孔径值和第二孔径值；
比如：第一孔径值为43.3um，第二孔径值为2.7um，从第二孔径值开始，每间隔一定数值，确定一个孔径值，得到24个孔径值：2.7、2.9、3.1、3.2、3.4、3.6、4.5、4.7、4.9、5.1、5.2、5.4、6.3、7.2、9.0、10.8、11.7、16.2、20.8、25.3、29.8、34.3、38.8、43.3。
本领域技术人员应当明白，M个孔径值在第一孔径值和第二孔径值之间应当尽量均匀，既不能大部分都靠近第一孔径值，也不能大部分都靠近第二孔径值。
对于A3，对于每个目标孔径值，通过模拟注汞法，获取目标孔径值对应的进汞体积和，进汞体积和为相应目标孔径值对应的孔隙体积和；
其中，获取一个目标孔径值对应的孔隙体积和的具体过程如下：
a1、对于注汞端面的每个第一点：(1)获取包括第一点的至少一条第一骨架线；
骨架图像中的骨架线是纵横交错的，第一点可能包括在多条骨架线中，本实施例中的第一骨架线是指包括第一点的一条没有分支的曲线，也就是说每个第一点对应的至少一条第一骨架线至少有一个交点—相应的第一点。
图6为本申请实施例提供的包括第一点的两条第一骨架线的示意图一，图7为本申请实施例提供的包括第一点的两条第一骨架线的示意图二。
参见图6，c-c1为一条第一骨架线，c-c2为一条第一骨架线，c-c1和c-c2的交点为第一点c。
参见图7，c-c3为一条第一骨架线，c-c4为一条第一骨架线，c-c3和c-c4之间具有重叠部分c-c5，其中，c为第一点。
(2)对于每条第一骨架线，确定每条第一骨架线包括的至少一个目标点，目标点为第一点和第二点之间的点，所述第二点为在第一方向上，第一骨架线上第一个对应的孔径值小于目标孔径值的点；第一方向为从注汞端面指向骨架图像的第二底面的方向，至少一个目标点包括第一点和所述第二点；
(3)根据所有第一骨架线包括的所有目标点，各自对应的孔径值，获取第一点对应的进汞体积和。
具体地，若目标孔径值为5.2um，包括当前第一点A的各第一骨架线为：第一骨架线1、第一骨架线2；
对于第一骨架线1，从当前第一点开始，沿第一方向，依次比较第一骨架线1上的点对应的孔径值与目标孔径值5.2um的大小，直至出现第一骨架线1上的第二点对应的孔径值小于目标孔径值5.2um时，停止比较。
比如，第一骨架线1在第一方向上依次具有包括点A1、A2、A3、A4、A5、A6、……、k、…、K，其中，点A1为当前的第一点A，也是第一骨架线1上的第一个点，点A2是第一骨架线1上的第二个点、以此类推，点k是第一骨架线1上的第k个点，点K是第一骨架线1上的第K个点
点A1、A2、各自对应的孔径值均大于或等于目标孔径值5.2um，点A3对应的孔径值小于目标孔径值5.2um，则点A3便为第二点。点A1、A2、A3为第一骨架线1对应的目标点。
第一骨架线2在第一方向上依次具有包括点B1、B2、B3、B4、B5、B6、……、l、…、L，其中，点B1为当前的第一点A，与点A1相同，也是第一骨架线2上的第一个点，点B2是第一骨架线2上的第二个点、以此类推，点l是第一骨架线1上的第k个点，点L是第一骨架线1上的第L个点。
点B1、B2、B3、各自对应的孔径值均大于或等于目标孔径值5.2um，点B4对应的孔径值小于目标孔径值5.2um，则点B4便为第一骨架线2上的第二点。点B1、B2、B3、B4为第一骨架线2对应的目标点。
若A1、A2、A3、B1、B2、B3、B4中，除了A1与B1为同一点外、A2与B2为同一点，其余的点不相同，则与第一点A对应的所有第一骨架线的所有目标点有A1、A2、A3、B3、B4。
第一点A对应的进汞体积和为：
其中，为第一骨架线1上第一个点(点A1、B1、A)对应的孔径值，为点A2对应的孔径值，为点B3对应的孔径值，为点B4对应的孔径值，为点A1和点A2之间的直线距离，为点A2和点A3之间的直线距离，点A2和点B3之间的直线距离，点B3和点B4之间的直线距离。
按照获取第一点A对应的孔隙体积的方法(步骤(1)、(2)、(3)的方法)，依次获取注汞端面上其它的第一点各自对应的孔隙体积和。
a2、根据所有第一点各自对应的孔隙体积和，得到该目标孔径值对应的孔隙体积和。
将注汞端面上所有第一点各自对应的孔隙体积和相加，得到当前目标孔径值对应的孔隙体积和。
按照上述a1～a2的方法，得到各目标孔径值各自对应的孔隙体积和。
对于A4、根据M个目标孔径值和每个目标孔径值对应的孔隙体积和，得到孔径分布。
具体地，图8为本申请实施例提供的孔径-孔隙体积和曲线，图9为本申请实施例提供的孔径-孔隙体积曲线，图10为本申请实施例提供的孔径分布柱状图。
比如，根据M个目标孔径值和每个目标孔径值对应的孔隙体积和，做孔径值-孔隙体积和曲线，如图8所示；图8的获取方法以孔径值为横坐标，孔径值对应的孔隙体积和作为纵坐标，将A4中得到的每个目标孔径值以及对应的孔隙体积确定的点连接成线得到的。
参见图8，根据图8中的曲线可以得到大于某一孔径值的所有的孔隙的体积之和与等于该孔径值对应的孔隙的体积的和。
对图8所示的曲线求导，得到图9，参见图9，图9中的曲线显示的信息为孔径从小到大变化时，孔径对应的孔隙的体积的变化。
参见图10，图10中显示的信息为各孔径范围内的孔隙的总体积，也就是孔径分布。比如根据图8中的曲线得到孔径10对应的孔隙体积和减去根据图8中的曲线得到孔径20对应的孔隙体积和，得到的差值为即为在10～20的孔径范围内的孔隙的总体积。
本实施例的方法，通过获取待测定岩石的三维数字图像，对三维数字图像进行处理，得到待测定岩心对应的孔隙的骨架图像，根据骨架图像包括的各点各自对应的孔径值，得到待测定岩石的孔隙的孔径分布，可测定的孔隙尺寸的范围不受限制，可获取各种待测定岩石的孔隙的孔径分布。
为了增加获得的孔径分布的准确性，本申请实施例在上一实施例的基础上做了进一步地改进。图11为本申请实施例提供的基于模拟压汞法获取岩石的孔径分布的方法的流程图二，如图11所示，本实施例的方法可以包括：
步骤S201、获取待测定岩石的岩心的至少部分的三维数字图像；
步骤S202、对岩心的至少部分的三维数字图像进行滤波，得到滤波后的三维数字图像；
步骤S203、根据该滤波后的三维数字图像，获取岩心的至少部分包括的孔隙的骨架图像；
步骤S204、获取该骨架图像包括的各点各自对应的孔径值；其中，点对应的孔径值为所述点距离孔隙边界的最小距离的两倍；
步骤S205、根据骨架图像包括的各点各自对应的孔径值，通过模拟注汞法获取待测定岩石的孔隙的孔径分布。
具体地，步骤S201参照上一实施例中的步骤S101。
对于步骤S202、图12为对图2中的图像滤波后的图像，图13为岩心的三维数字图像滤波后的立体图像效果图。
滤波方法可为均值滤波、非局部均值滤波、中值滤波、维纳滤波、高斯滤波等。
比如采用非局部均值滤波方法对图2中的三维数字图像进行滤波，得到滤波后的三维数字图像，如图12和图13所示。
对于步骤S203～步骤S205，参照上一实施例中的步骤S102～步骤S104，与步骤S102～步骤S104不相同的是，获取岩心的孔隙的骨架图像基于的图像为滤波后的三维数字图像，这样获取到的岩心的孔隙的骨架图像更准确，最终得到的孔径分布也更准确。
本实施例中的基于模拟压汞法获取岩石的孔径分布的方法，不受孔隙大小的限制，可以获取任何岩心的孔径分布，且得到的孔径分布的准确度更高。
图14为本申请实施例提供的取岩石的孔径分布的装置的结构示意图一，如图14所示，本实施例的装置可以包括：三维数字图像获取模块41、骨架图像获取模块42和孔径分布获取模块43。
其中，所述三维数字图像获取模块41，用于获取待测定岩石的岩心的至少部分的三维数字图像；
所述骨架图像获取模块42，用于根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像；
所述骨架图像获取模块42，还用于获取所述骨架图像包括的各点各自对应的孔径值；所述点的孔径值为所述点距离孔隙边界的最小距离的两倍；
所述孔径分布获取模块43，用于根据所述骨架图像包括的各点各自对应的孔径值，获取待测定岩石的孔隙的孔径分布。
本实施例的装置，可以用于执行上述方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
在一种可能的设计中，所述骨架图像获取模块42，具体用于：
对所述三维数字图像进行图像分割，得到所述岩心的至少部分包括的孔隙的图像；
根据所述孔隙的图像，获取所述骨架图像。
在一种可能的设计中，所述骨架图像获取模块42，具体用于：
对所述孔隙的图像进行距离变换，得到第一图像；
对所述第一图像进行三维细化处理，得到所述骨架图像。
在一种可能的设计中，所述孔径分布获取模块43，具体用于：
确定注汞端面上所有的第一点对应的孔径值中，最大的第一孔径值，以及所述各点对应的孔径值中，最小的第二孔径值；所述注汞端面为所述骨架图像的第一底面；
从所述第一孔径值和所述第二孔径值之间的孔径值中确定出M个目标孔径值，所述M个目标孔径值包括所述第一孔径值和所述第二孔径值；
对于每个目标孔径值，通过模拟注汞法，获取所述孔径值对应的进汞体积和，所述进汞体积和为相应孔径值对应的孔隙体积和；
根据所述M个目标孔径值和每个目标孔径值对应的进汞体积和，得到所述待测定岩石的孔隙的孔径分布。
在一种可能的设计中，所述孔径分布获取模块43，具体用于：
对于所述注汞端面的每个第一点，获取包括所述第一点的至少一条第一骨架线；
对于每条第一骨架线，确定每条第一骨架线包括的至少一个目标点，所述目标点为第一点和第二点之间的点，所述第二点为在第一方向上，所述第一骨架线上第一个对应的孔径值小于所述目标孔径值的点；所述第一方向为从所述注汞端面指向所述骨架图像的第二底面的方向，至少一个目标点包括第一点和所述第二点；
根据所有第一骨架线包括的所有目标点，各自对应的孔径值，获取所述第一点对应的进汞体积和；
根据所有第一点各自对应的孔隙体积和，得到所述目标孔径值对应的进汞体积和。
本实施例的装置，可以用于执行上述方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
图15为本申请实施例提供的取岩石的孔径分布的装置的结构示意图二，如图15所示，本实施例的装置在图14所示的装置的基础上还可以包括：滤波模块44；
滤波模块44用于在根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像之前，对所述三维数字图像进行滤波，得到滤波后的三维数字图像；
则骨架图像获取模块42，用于根据所述滤波后的三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像。
本实施例的装置，可以用于执行上述方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是：以上各实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述各实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的范围。
Claims (8)
1.一种基于模拟压汞法获取岩石的孔径分布的方法，其特征在于，包括：
获取待测定岩石的岩心的至少部分的三维数字图像；
根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像；
获取所述骨架图像包括的各点对应的孔径值；其中，点对应的孔径值为所述点距离孔隙边界的最小距离的两倍；
根据所述骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布，具体包括：
确定注汞端面上所有的第一点对应的孔径值中，最大的第一孔径值，以及所述各点对应的孔径值中，最小的第二孔径值；所述注汞端面为所述骨架图像的第一底面；
从所述第一孔径值和所述第二孔径值之间的孔径值中确定出M个目标孔径值，所述M个目标孔径值包括所述第一孔径值和所述第二孔径值；
对于每个目标孔径值，通过模拟压汞法，获取所述孔径值对应的进汞体积和，所述进汞体积和为相应孔径值对应的孔隙体积和；
根据所述M个目标孔径值和每个目标孔径值对应的孔隙体积和，得到所述待测定岩石的孔隙的孔径分布。
2.根据权利要求1所述的方法，其特征在于，所根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像，包括：
对所述三维数字图像进行图像分割，得到所述岩心的至少部分包括的孔隙的图像；
根据所述孔隙的图像，获取所述骨架图像。
3.根据权利要求2所述的方法，其特征在于，所述根据所述孔隙的图像，获取所述骨架图像，包括：
对所述孔隙的图像进行距离变换，得到第一图像；
对所述第一图像进行三维细化处理，得到所述骨架图像。
4.根据权利要求1所述的方法，其特征在于，所述通过模拟注汞法，获取所述孔径值对应的进汞体积和，包括：
对于所述注汞端面的每个第一点，获取包括所述第一点的至少一条第一骨架线；
对于每条第一骨架线，确定每条第一骨架线包括的至少一个目标点，所述目标点为第一点和第二点之间的点，所述第二点为在第一方向上，所述第一骨架线上第一个对应的孔径值小于所述目标孔径值的点；所述第一方向为从所述注汞端面指向所述骨架图像的第二底面的方向，至少一个目标点包括第一点和所述第二点；
根据所有第一骨架线包括的所有目标点，各自对应的孔径值，获取所述第一点对应的进汞体积和；
根据所有第一点各自对应的孔隙体积和，得到所述目标孔径值对应的进汞体积和。
5.根据权利要求1～4任一项所述的方法，其特征在于，在所述根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像之前，还包括：
对所述三维数字图像进行滤波，得到滤波后的三维数字图像；
则所述根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像，包括：
根据所述滤波后的三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像。
6.一种基于模拟压汞法获取岩石的孔径分布的装置，其特征在于，包括：
三维数字图像获取模块，用于获取待测定岩石的岩心的至少部分的三维数字图像；
骨架图像获取模块，用于根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像；
骨架图像获取模块，还用于获取所述骨架图像包括的各点各自对应的孔径值；点对应孔径值为所述点距离孔隙边界的最小距离的两倍；
孔径分布获取模块，用于根据所述骨架图像包括的各点各自对应的孔径值，通过模拟压汞法获取待测定岩石的孔隙的孔径分布；
所述孔径分布获取模块，具体用于：
确定注汞端面上所有的第一点对应的孔径值中，最大的第一孔径值，以及所述各点对应的孔径值中，最小的第二孔径值；所述注汞端面为所述骨架图像的第一底面；
从所述第一孔径值和所述第二孔径值之间的孔径值中确定出M个目标孔径值，所述M个目标孔径值包括所述第一孔径值和所述第二孔径值；
对于每个目标孔径值，通过模拟压汞法，获取所述孔径值对应的进汞体积和，所述进汞体积和为相应孔径值对应的孔隙体积和；
根据所述M个目标孔径值和每个目标孔径值对应的孔隙体积和，得到所述待测定岩石的孔隙的孔径分布。
7.根据权利要求6所述的装置，其特征在于，所述孔径分布获取模块，具体用于：
对于所述注汞端面的每个第一点，获取包括所述第一点的至少一条第一骨架线；
对于每条第一骨架线，确定每条第一骨架线包括的至少一个目标点，所述目标点为第一点和第二点之间的点，所述第二点为在第一方向上，所述第一骨架线上第一个对应的孔径值小于所述目标孔径值的点；所述第一方向为从所述注汞端面指向所述骨架图像的第二底面的方向，至少一个目标点包括第一点和所述第二点；
根据所有第一骨架线包括的所有目标点，各自对应的孔径值，获取所述第一点对应的进汞体积和；
根据所有第一点各自对应的孔隙体积和，得到所述目标孔径值对应的进汞体积和。
8.根据权利要求6～7任一所述的装置，其特征在于，所述装置还包括滤波模块，所述滤波模块用于在根据所述三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像之前，对所述三维数字图像进行滤波，得到滤波后的三维数字图像；
则所述骨架图像获取模块，用于根据所述滤波后的三维数字图像，获取所述岩心的至少部分包括的孔隙的骨架图像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810227776.5A CN108492331B (zh) | 2018-03-20 | 2018-03-20 | 基于模拟压汞法获取岩石的孔径分布的方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810227776.5A CN108492331B (zh) | 2018-03-20 | 2018-03-20 | 基于模拟压汞法获取岩石的孔径分布的方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108492331A CN108492331A (zh) | 2018-09-04 |
| CN108492331B true CN108492331B (zh) | 2019-04-23 |
Family
ID=63318497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810227776.5A Active CN108492331B (zh) | 2018-03-20 | 2018-03-20 | 基于模拟压汞法获取岩石的孔径分布的方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108492331B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114428043B (zh) * | 2020-09-18 | 2024-06-21 | 中国石油化工股份有限公司 | 多孔介质孔径分布表征方法及电子设备 |
| CN113176187A (zh) * | 2021-03-23 | 2021-07-27 | 中国地质大学（武汉） | 一种表征岩石孔隙网络连通性的方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101802649A (zh) * | 2008-04-10 | 2010-08-11 | 普拉德研究及开发股份有限公司 | 利用井眼图像、数字岩石样品以及多点统计算法生成数值假岩心的方法 |
| CN103822865A (zh) * | 2014-03-20 | 2014-05-28 | 中国石油大学(华东) | 一种高分辨率三维数字岩心建模方法 |
| CN104634718A (zh) * | 2015-03-05 | 2015-05-20 | 中国石油大学(华东) | 应用核磁共振表征致密砂岩孔径分布的标定方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105372168B (zh) * | 2015-12-08 | 2018-03-23 | 清华大学 | 一种混凝土孔隙三维结构的构建方法 |
| CN106383078B (zh) * | 2016-09-20 | 2019-05-07 | 中国石油天然气股份有限公司 | 岩石水驱效率的确定方法及装置 |
-
2018
- 2018-03-20 CN CN201810227776.5A patent/CN108492331B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101802649A (zh) * | 2008-04-10 | 2010-08-11 | 普拉德研究及开发股份有限公司 | 利用井眼图像、数字岩石样品以及多点统计算法生成数值假岩心的方法 |
| CN103822865A (zh) * | 2014-03-20 | 2014-05-28 | 中国石油大学(华东) | 一种高分辨率三维数字岩心建模方法 |
| CN104634718A (zh) * | 2015-03-05 | 2015-05-20 | 中国石油大学(华东) | 应用核磁共振表征致密砂岩孔径分布的标定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108492331A (zh) | 2018-09-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10540760B2 (en) | Method of separating, identifying and characterizing cracks in 3D space | |
| CN106780744B (zh) | 利用不同分辨率ct图像构建多尺度三维数字岩心的方法 | |
| CN108802073B (zh) | 一种基于数字岩心的岩石电性参数获取方法及装置 | |
| CN103698803B (zh) | 一种岩石孔隙结构表征方法及装置 | |
| Bauer et al. | From computed microtomography images to resistivity index calculations of heterogeneous carbonates using a dual-porosity pore-network approach: Influence of percolation on the electrical transport properties | |
| CN105074456B (zh) | 用于改进由岩石样品的材料特性的直接数值模拟以及确定材料特性的不确定性的系统和方法 | |
| CN107709699B (zh) | 生成多孔岩石样本的三维微型模型 | |
| CN109242985B (zh) | 一种从三维图像确定孔隙结构关键参数的方法 | |
| WO2017041281A1 (en) | Porous media anaylysis system and method | |
| CN115398216A (zh) | 用于估计岩石的流体饱和度的方法 | |
| EP2909658A2 (en) | Method for modeling a reservoir using 3d multiple-point simulations with 2d training images | |
| CN108492331B (zh) | 基于模拟压汞法获取岩石的孔径分布的方法和装置 | |
| CN109448104B (zh) | 一种基于岩心图像的拓展重构方法及装置 | |
| CN111833432A (zh) | 一种基于岩心二维灰度图像的三维重建方法 | |
| US20240070811A1 (en) | Methods and devices for enhancement of porous media images | |
| US20240070356A1 (en) | Methods and devices for stencil-based pore network construction | |
| Dong et al. | Verification of the carbonate double-porosity conductivity model based on digital cores | |
| CN114705606B (zh) | 一种基于网络化分析的岩石内部关键渗流节点的封堵方法 | |
| CN113592829B (zh) | 基于分割重组的深度学习泥沙颗粒识别方法及装置 | |
| CN113781443A (zh) | 一种基于数字图像颜色梯度的岩层裂隙分布规律获取方法、系统、终端及可读存储介质 | |
| RU2774959C1 (ru) | Способ определения фильтрационных свойств неоднородных пористых образцов | |
| CN111986078A (zh) | 一种基于引导数据的多尺度岩心ct图像融合重建的方法 | |
| CN114022631B (zh) | 基于深度卷积生成对抗网络的三维地质模型自动重构方法 | |
| Schembre-McCabe et al. | Two examples of adding value through Digital Rock Technology | |
| Zhang | Comparative Study of Pore Structure Characterization for Different Rock Samples Using Digital Images |
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 |