《数字岩芯CT图像的二值化处理方法》.pdf

ICS

CCS点击此处添加CCS号

GRM

中关村绿色矿山产业联盟团体标准

T/GRM001—2024

数字岩芯CT图像的二值化处理方法

BinarizationstandardfordigitalcoreCTimages

（征求意见稿）

在提交反馈意见时，请将您知道的相关专利连同支持性文件一并附上。

2024-08-01发布2024-12-31实施

中关村绿色矿山产业联盟  发布

T/GRM001—2024

数字岩芯CT图像的二值化处理方法

1范围

本文件规定了数字岩芯CT图像的二值化处理方法的原理、处理流程、计算方法、准确度和试验报告

等。

本文件适用于数字岩芯CT图像的二值化处理方法。

2规范性引用文件

下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，

仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本（包括所有的修改单）适用于本

文件。

GB/T4756石油液体手工取样法

GB/T6682分析实验室用水规格和实验方法

3术语和定义

下列术语和定义适用于本文件。

3.1绝对孔隙度

岩石的绝对孔隙体积与岩石的表面体积的比值。

3.2阈值灰度

将图像中的所有亮度值根据指定的亮度值分成高于阈值和低于阈值的两类。

4处理流程

在进行了高分辨率的全岩心建模后，就可以应用岩心实测氦气孔隙度与核磁T2谱进行标定，反向确

定阈值，达到计算孔隙度与实测孔隙度的统一，但由于切片上往往存在大量像素点的灰度等于阈值灰度，

为尽可能合理的完成计算孔隙度与实际孔隙含量的完全统一，本标准建立了如下分割方案：

a)通过全岩心CT扫描在一定分辨率下获得覆盖整块岩心视域的三维CT数据；

b)通过实测氦气孔隙度与核磁T2谱反映的孔径分布，可以获得在此扫描分辨率下的实际孔隙

含量；

c)通过确定三维CT数据体的阈值灰度，将灰度小于阈值灰度的所有像素点均计为孔隙点；

d)在全岩心模型内寻找所有灰度等于阈值灰度的像素点，并统计这些像素点相邻的26个体素单

元内（正方体区域内，如图1）的孔隙点个数；

e)按照统计的孔隙点个数（0到26个）将灰度等于阈值灰度的像素点分为27类，周围孔隙点个

数越多的像素将其优先计为孔隙点，依此确定阈值类别N，高于这一类别N的均计为孔隙（如

图2）；

f)由于阈值类别N内依旧存在多个像素点待分割（如图3红框内），以实现计算孔隙度与实际

孔隙含量完全吻合，但该类内的像素点周围孔隙点数已相等，这里仅将类内像素点做随机选

取，直至计算孔隙度与实际孔隙含量完全相等，并以此为二值化处理效果的判别标准，最终

输出二值化图像（如图3）。

1

T/GRM001—2024

图1领域搜索范围示意图

图2统计结果选取示意图

图3某砂岩样品的CT数据的二值化处理结果，基质（左）与孔隙结构（右）

5计算方法

用式（1）计算当前分辨率下CT数据的实际孔隙含量：