一种基于数字岩石图像分析的复杂岩性识别方法
郝洋①,, 杜焕福①,, 万亚旗①,, 曾翔①,, 王鑫①,
①中石化经纬有限公司地质测控技术研究院
②中国石化测录井重点实验室

作者简介:郝洋 研究实习员(自然科学),1999年生,2024年毕业于中国石油大学(华东)控制科学与工程专业,现在中石化经纬有限公司地质测控技术研究院主要从事录井设备数智化研发工作。通信地址:266001 山东省青岛市市南区台湾路4号。E-mail:haoyanghy@yeah.net

摘要

随着录井分析技术的长足进展,以RoqSCAN为代表的岩屑矿物识别技术已成为录井分析的重要工具。然而,该技术目前仅能测试岩屑颗粒的矿物组成,无法实现岩性自动识别,对复杂岩性更难以鉴定。为此,提出了一种基于数字岩石图像分析的复杂岩性识别方法:利用边缘检测与图像分割技术,实现岩屑颗粒的自动提取;结合颗粒形态和矿物信息,开展颗粒级别的岩性识别;通过对识别结果进行统计与可视化分析,为复杂岩性的识别与分类提供科学、可靠的数据支持。目前,该识别方法已在鄂尔多斯盆地D 2井成功应用,显著提升了复杂地层的识别准确度和效率,为录井现场岩屑高效识别提供了新的技术方法。

关键词: 复杂岩性; RoqSCAN; 数字岩石图像; 图像分割; 颗粒形态
中图分类号:TE132.1 文献标志码:A
A method for complex lithology identification based on digital rock image analysis
HAO Yang①,, DU Huanfu①,, WAN Yaqi①,, ZENG Xiang①,, WANG Xin①,
①Geological Measurement & Control Technology Research Institute, Sinopec Matrix Corporation, Qingdao, Shandong 266001, China
②Sinopec Key Laboratory, Qingdao, Shandong 266500, China
Abstract

With the significant progress in mud logging analysis technology,cuttings mineralogical identification technology represented by RoqSCAN has become an important tool in mud logging analysis. However,this technology currently can only test the mineral composition within cuttings particles and cannot achieve automatic lithology identification,making it even more difficult to identify complex lithologies. Therefore, a complex lithology identification method based on digital rock image analysis was proposed, i.e., using edge detection and image segmentation techniques to achieve the automatic extraction of rock particles, combining particle morphology and mineral information to conduct particle-level lithology identification, and providing scientific and reliable data support for the identification and classification of complex lithologies through statistical and visual analysis of the identification results. At present, this recognition method has been successfully applied to well D 2 in Ordos Basin, significantly improving the accuracy and efficiency of identifying complex strata, and providing a new technical method for efficient identification of the cuttings in mud logging sites.

Keyword: complex lithology; RoqSCAN; digital rock image; image segmentation; particle morphology
0 引言

岩屑的准确识别是实现地层划分和储层评价的关键环节[1, 2, 3]。然而随着勘探开发向复杂储层拓展, 岩性愈加复杂, 给现场录井岩性识别带来巨大挑战。RoqSCAN多功能矿物分析系统[4]作为录井行业新兴分析技术, 利用扫描电镜和能谱探头扫描岩屑颗粒, 获取颗粒表面的结构特征与矿物成分信息, 可实现岩屑颗粒的矿物定量识别[5]。尽管该技术在岩屑的矿物分析中表现出强大能力[6], 但输出数据却涵盖图像、矿物含量等多个维度, 结果解释仍高度依赖人工, 对操作人员的专业性和经验要求较高, 难以满足现场高频次、快速决策的应用需求。

图像分割技术[7]能够精准高效地从大量岩屑样本中提取单个颗粒的轮廓与边界, 实现颗粒尺度级别观察和分析[8]。利用图像识别技术, 结合纹理分析与矿物成分信息, 可对海量岩屑颗粒岩性进行自动分类与识别, 极大地提升了识别的速度和准确性, 减少了对人工操作的依赖[9], 实现全过程的数字化管理与处理。同时, 基于数据可视化和交互式分析平台的构建[10], 能够将复杂的多个维度的数据以直观、生动的形式展示, 帮助地质人员快速理解和判断, 提升决策的科学性和实时性[11, 12]

本文针对 RoqSCAN 在实际操作中存在的问题, 结合当前蓬勃发展的图像分析技术与工作实践, 提出一种基于数字岩石图像分析的复杂岩性识别新方法。该方法通过边缘检测与图像分割实现岩屑颗粒的自动提取, 结合颗粒形态及其矿物组成信息开展岩性智能识别, 构建可视化展示与统计分析模块, 形成一套面向复杂地层的智能识别体系, 不仅提升了岩性识别的自动化程度与精度, 也有效降低了解析成本, 为深层与复杂区油气目标的快速识别与评价提供了技术支撑。

1 技术思路

依托 RoqSCAN多功能矿物分析系统获取高分辨率岩屑样本图像及矿物成分图像, 开展图像处理与联合识别工作。为实现岩屑样本岩性的准确分析与复杂储层的精准评价, 本研究采用基于颗粒尺度的分析思路。

首先, 对岩屑样本图像与矿物成分图像进行图像增强与配准处理, 以提升图像质量并确保空间对齐; 其次, 应用图像分割算法提取颗粒目标区域, 并结合颗粒的形貌特征与其对应的矿物组成, 实现单颗粒岩性的识别与分类; 最后, 对颗粒级岩屑分类结果进行统计分析与可视化, 定量给出样本中各类岩性颗粒的数量及其占比, 为样本尺度的岩性分析与描述提供可靠依据。完整的技术路线如图1所示。

2 数字岩石图像分析方法
2.1 RoqSCAN岩屑样本图像采集

RoqSCAN多功能矿物分析系统(图2)是一种能够对岩屑颗粒进行矿物定量分析的先进系统, 广泛应用于地质录井和储层评价领域[13]。该系统利用扫描电镜和能谱分析, 快速获取岩屑的矿物组成、颗粒结构、孔隙度和元素分布等多维信息。相比传统人工识别方法, RoqSCAN不仅能提高识别的准确性和一致性, 还可以实现岩性分析的数字化和标准化。特别是在复杂地层条件下, RoqSCAN可作为图像识别的重要补充, 为构建岩性智能识别体系提供了关键的矿物成分数据支撑。

图2 RoqSCAN多功能矿物分析系统

RoqSCAN多功能矿物分析系统通过扫描电镜获取高分辨率的岩屑颗粒图像, 可清晰呈现颗粒的形貌、边界及表面特征。在此基础上, 利用能谱仪对岩屑颗粒表面的矿物类型进行识别, 从而生成反映样本矿物成分分布图。最终可以得到岩屑样本图像(图3a)和岩屑矿物成分图像(图3b)。

图3 RoqSCAN岩屑分析图像

2.2 图像预处理

岩屑样本图像中黑色背景为基质, 灰、白色等亮色区域为岩屑颗粒。本研究采用颗粒尺度的岩性分析方法, 首先对图像进行预处理, 以便去除背景基质对后续颗粒分割的干扰[14]; 随后基于颗粒轮廓特征应用图像分割技术, 对不同尺寸的岩屑颗粒进行精确分离。为提升分析准确性, 分割前可设定最小颗粒面积阈值, 自动筛除尺寸较小、无分析价值的岩屑颗粒。具体预处理过程如下:

二值化处理:将岩屑样本图像转换为灰度图, 并进行二值化处理, 设定固定阈值(如50), 将灰度图进行阈值分割, 生成黑白二值图像, 从而有效凸显岩屑颗粒区域[15]

颗粒筛选:手动输入筛选最小颗粒直径与比例尺, 按比例尺计算最小颗粒面积阈值, 筛选出面积大于该阈值的颗粒, 保留颗粒主体并将背景部分设置为透明。

经过图像预处理后得到的岩屑样本图像如图4所示。矿物成分图像仅进行颗粒筛选预处理, 其核心目标在于实现图像对齐, 确保在岩屑样本图像中提取颗粒后, 能够在对应位置准确获取其矿物组成信息。预处理过程不仅有效去除了背景基质的干扰, 还对直径较小的颗粒进行了筛选与剔除, 从而保留了更具代表性的颗粒区域, 为后续颗粒级岩性分析提供可靠基础。

图4 岩屑样本图像预处理前后对比

2.3 图像分割

利用边缘检测算法对岩屑样本图像进行处理, 可有效识别图像中的外部轮廓, 初步定位各个岩屑颗粒的位置[16]。在此基础上, 生成对应的掩膜图像, 进一步从原始图像中提取出各颗粒所对应的区域, 形成清晰、准确的颗粒分割结果图像。然后, 根据每一个颗粒的空间位置信息, 在矿物成分图像中提取对应区域, 从而获取该颗粒的详细矿物组成。这一过程实现了从“颗粒形态”到“矿物成分”的一一对应, 具备较高的空间精度和分析效率。

通过该方法, 可以实现对岩屑颗粒的精准分割, 并且利用空间位置信息实现岩屑颗粒与其对应的矿物信息进行匹配, 分割后的岩屑样本图像及其对应的矿物成分图像如图5所示。

图5 分割后的岩屑样本图像及其对应的矿物成分信息展示

2.4 岩屑颗粒岩性识别与分类

在完成岩屑颗粒的图像分割后, 本研究基于岩屑样本图像与矿物成分图像的联合分析, 实现颗粒级岩性的识别[17]。岩屑样本图像提供了岩屑的岩石学特征(如粒度、结构), 而矿物成分图像则揭示了岩屑的矿物种类。

首先, 依据RoqSCAN获得的矿物组分数据, 将岩屑颗粒按主导矿物成分划分为3大类岩性:以石英矿物为主的石英类岩屑、以碳酸盐矿物为主的碳酸盐类岩屑, 以及以泥质矿物为主的泥岩类岩屑。3类主导矿物采用不同颜色加以区分:石英矿物标示为黄色, 碳酸盐矿物标示为蓝色, 泥质矿物标示为绿色。为提升颜色判别的稳定性, 本研究将原始RGB图像转换至HSV颜色空间。HSV以色调、饱和度和明度表征颜色, 更契合人眼感知特征, 能够有效增强颜色阈值分割的鲁棒性[18]。依据预设的色调范围提取对应矿物类型, 统计各类颜色像素在单个颗粒中的占比, 并将占比最大的颜色像素所对应的矿物类型确定为该颗粒的主导矿物, 进而完成颗粒岩性的判定; 当与主导矿物对应的颜色像素占比低于10%时, 即认为该颗粒内部矿物特征分散、主导性不足, 将其统一归类为“其他”, 以降低噪声干扰。

虽完成了岩性大类的识别, 但在石英类岩屑颗粒中, 因硅质岩和砂岩的矿物组成高度相似, 难以仅凭成分实现准确区分[19]。为提高岩性划分的准确性与实用性, 还需进一步结合岩屑样本图像的结构特征(如颗粒形态、纹理), 对石英类岩屑颗粒进行更精细的分类与识别。

通过研究大量的岩屑样本发现, 硅质岩、砂岩和泥质砂岩3类岩屑颗粒的矿物成分大部分为石英, 但在岩屑样本图像的灰度分布特征上具有明显差异, 这为岩性细分类提供了可量化的依据。具体而言:硅质岩颗粒表面相对平整、质地均匀, 其灰度直方图通常仅出现一个明显主峰, 反映出图像灰度分布的集中性[20]; 砂岩由于表面粗糙且夹杂多种矿物杂质, 灰度变化较为复杂, 通常出现两个主峰, 且两个峰值的高度差距较小, 峰值比小于10[21]; 泥质砂岩同样具有双峰特征, 但因其杂质含量相对较低, 灰度分布更具倾向性, 峰值比通常大于或等于10[22]。因此, 可将灰度直方图的主峰数量与峰值比, 作为石英类岩屑细分类的重要依据(图6)。

图6 石英类岩屑岩性细分示意图

基于上述流程, 岩屑颗粒首先依据矿物成分划分为石英类、碳酸盐类与泥岩类3大类, 随后对石英类岩屑进一步按灰度特征细分为硅质岩、砂岩和泥质砂岩, 从而实现颗粒级的精细岩性识别。

2.5 统计分析与可视化

通过对岩屑颗粒的识别与分类能够获取每个颗粒的岩性信息, 但颗粒层面的数据无法完整反映该深度区间内地层的整体特征, 因此需要对该深度区间内岩屑样本中所有颗粒的岩性信息进行统计与分析。基于此, 本研究对不同深度的岩屑样本进行系统统计与可视化处理, 提取颗粒数量、颗粒岩性类型及其占比等定量结构信息。通过生成岩性分类图、颗粒占比图和岩性统计Excel表格, 可直观展现岩性随深度的变化趋势, 为后续的地质分析与复杂储层评价提供有效的数据支持。

对每个岩屑样本, 通过图像分析方法识别其中每个颗粒的岩性, 获取颗粒的岩性信息和面积大小, 为后续统计分析提供基础。基于识别结果, 生成该样本的颗粒岩性分类图(图7)。在颗粒岩性分类图中, 每行对应同一岩性类别的颗粒, 颗粒按面积大小排序, 可直观展示不同岩性颗粒的数量和分布情况, 不同岩性的矿物成分图则用于辅助验证分类准确性。通过颗粒岩性分类图, 可以快速了解样本中各岩性颗粒的占比和分布, 为深入分析岩性特征及储层评价提供数据支撑。

图7 颗粒岩性分类图

岩屑颗粒分类图只能对不同岩性颗粒数量及分布提供直观感受, 难以实现准确的定量分析。为此, 在完成颗粒级岩性分类的基础上, 对样本中不同岩性颗粒的数量进行统计汇总, 计算各类岩性颗粒在样本中的占比, 如图8所示。该统计结果以扇形展示, 能够清晰、准确地反映样本中各类岩性的颗粒数量比例, 为后续地层岩性评价提供可靠的定量依据。

图8 各类岩性颗粒占比图

目前, 已对单个岩屑样本获取了颗粒岩性信息, 并统计各类岩性颗粒的数量及占比。为了分析整段地层的岩性特征, 需要对多个连续样本进行系统汇总。具体做法是, 将每个样本中各类岩性颗粒的数量及占比整理输入至Excel表格中, 通过该表可以直观了解各类岩性在不同深度样本中的分布情况及相对占比。这些统计结果反映每个样本中不同岩性的颗粒占比, 是客观的原始数据。通过对Excel表中的占比变化进行分析, 可以识别地层中岩性的主导类型, 观察各岩性随深度的变化趋势, 并发现岩性演化的规律和异常点, 为后续的地质分析和复杂储层评价提供可靠的数据支持。

3 实例应用
3.1 工程概况

本文将数字岩石图像分析方法应用于鄂尔多斯盆地的D 2井。该井在井段5 500~6 427 m钻遇寒武系至长城系地层, 由于地层连续沉积, 岩性变化规律不明显, 导致现场岩屑识别难以准确卡取层位。本研究尝试利用RoqSCAN分析和图像识别来解决此问题。

在上述重点钻探井段内, 采集了共计61个代表性岩屑, 利用RoqSCAN 多功能矿物分析系统对样本进行高精度扫描和分析, 得到了61个样本的高分辨率岩屑样本图像和矿物成分图像。随后, 通过自主研发的软件系统, 结合岩屑样本图像和矿物成分信息, 通过图像处理技术, 对岩屑样本进行颗粒尺度下的分割与岩性识别。通过应用软件最终生成了各深度点具体的岩性组分占比图及全井段综合统计结果。

3.2 岩性分析

基于图像处理技术, 可以对所有岩屑样本中的颗粒进行岩性分类, 并形成可统计的数据结果, 从而直观展示各类岩性颗粒在样本中的占比, 进而得到该深度下较为客观的岩性判别结果。为了进一步说明方法的应用效果, 以下将分别以井深5 500、6 420 m处的岩性分析为例进行展示。

(1)采用数字图像分析方法对5 500 m处样本开展岩性分析。颗粒岩性分类图(图9)表明, 样本以碳酸盐岩颗粒为主, 并含少量泥岩、砂岩等碎屑成分。基于颗粒级统计结果得到的各类岩性颗粒占比图(图10)显示, 该深度样本中碳酸盐岩颗粒占比超过85%, 占据绝对优势。基于上述定量统计结果, 并结合样本原始图像中颗粒形态与结构特征及RoqSCAN矿物组分分析结果综合判识, 碳酸盐岩以白云石为主、该深度地层属于徐庄组。

图9 5 500 m处颗粒岩性分类图

图10 5 500 m处岩性颗粒占比图

(2)对井深6 420 m处的岩性进行判定。岩屑颗粒分类图(图11)显示, 该深度以石英类岩屑颗粒为主, 且整体粒径较大。由分类统计结果得到各类岩性颗粒占比图(图12)表明, 硅质岩、砂岩及泥质砂岩合计占比约70%, 泥岩占比约25%, 碳酸盐岩含量仅为5%。通过岩性含量对比可知, 此样品碳酸盐岩含量大幅降低, 地层以砂岩及硅质岩为主, 其岩性特征与中元古界长城系地层的岩性吻合度高, 因此判定地层进入长城系。

图11 6 420 m处颗粒岩性分类图

图12 6 420 m处岩性颗粒占比图

4 总结与展望

(1)提出了一种基于图像处理的颗粒自动提取方法。采用边缘检测与图像分割技术, 实现了高分辨率岩屑图像中颗粒区域的精准提取与轮廓定位, 为后续岩性识别奠定了数据基础。

(2)构建了岩屑颗粒智能识别模型。融合矿物成分、纹理结构等图像特征, 结合颜色空间分析与分类算法, 实现了颗粒尺度下岩屑典型岩性的识别与分类。

(3)开发了岩性统计与可视化分析模块:基于识别结果, 构建了颗粒岩性组成比例统计模型, 并实现了图形化展示, 支持对不同深度样本的岩性占比纵向对比分析, 有效揭示岩性分布规律与储层层序特征。

(4)当前研究主要针对样本中岩屑颗粒级别的岩性分类, 在此基础上, 后续将进一步开展样本层面的综合分析, 并结合人工解释结果对分类结果进行对比验证, 以逐步完善岩性分类模型并评估识别效果, 从而提升该方法在地质分析中的应用价值。

编辑 郑春生

参考文献
[1] HUO F C, LI A, ZHAO X Q, et al. Novel lithology identification method for drilling cuttings under PDC bit condition[J]. Journal of Petroleum Science and Engineering, 2021, 205: 108898. [本文引用:1]
[2] TAMAAZOUSTI Y, FRANÇOIS M, KHERROUBI J. Automated identification and quantification of rock types from drill cuttings[C]//SEG Technical Program Expand ed Abstracts 2020. Tulsa: Society of Exploration Geophysicists, 2020: 1591-1595. [本文引用:1]
[3] 中国石油化工集团有限公司, 中石化石油工程技术服务有限公司, 中石化经纬有限公司, . 基于激光元素信息的砂岩-泥岩识别图版识别岩性方法: 201911418647. 5[P]. 2021-07-16.
China Petrochemical Corporation, Sinopec Oilfield Service Corporation, Sinopec Matrix Corporation, et al. Identifying lithology method based on laser element information for sand stone-mudstone identification chart: 201911418647. 5[P]. 2021-07-16. [本文引用:1]
[4] WANG X, HE J Y, QIU X Q, et al. An improved acid fracturing design through RoqSCAN technology: A case study from Daniudi Gas Field in Ordos Basin, China[J]. Processes, 2023, 11(5): 1475. [本文引用:1]
[5] ASHTON T, LY C V, SPENCE G. Application of real-time wellsite tool for enhanced geosteering and reservoir characterization[C]//Gulf Coast Association of Geological Societies Transactions. Houston: Gulf Coast Association of Geological Societies, 2013: 13-22. [本文引用:1]
[6] 王志战, 朱祖扬, 李丰波, . 便携式岩屑声波录井系统研制与测试[J]. 石油钻探技术, 2020, 48(6): 109-115.
WANG Zhizhan, ZHU Zuyang, LI Fengbo, et al. Development and testing of a portable acoustic logging system on cuttings[J]. Petroleum Drilling Techniques, 2020, 48(6): 109-115. [本文引用:1]
[7] CHENG Y, LI B J. Image segmentation technology and its application in digital image processing[C]//2021 IEEE Asia-Pacific Conference on Image Processing, Electronics and Computers (IPEC). Dalian, China: IEEE, 2021: 1174-1177. [本文引用:1]
[8] 潘显珊, 王正勇, 罗彬彬, . 基于DAF-STDC实时模型的岩屑图像语义分割识别[J]. 计算机系统应用, 2024, 33(12): 222-230.
PAN Xianshan, WANG Zhengyong, LUO Binbin, et al. Semantic segmentation and recognition of rock debris image based on real-time model DAF-STDC[J]. Computer Systems & Applications, 2024, 33(12): 222-230. [本文引用:1]
[9] 杜焕福, 孙鑫, 王春伟, . 东营凹陷页岩油双甜点录井解释评价方法研究[J]. 矿产勘查, 2023, 14(3): 480-490.
DU Huanfu, SUN Xin, WANG Chunwei, et al. Study on mud logging interpretation and evaluation method of geological and engineering sweet spots for shale oil in Dongying Sag[J]. Mineral Exploration, 2023, 14(3): 480-490. [本文引用:1]
[10] 罗崇兴, 师明元, 王正勇, . 一种改进DeepLabV3+的岩屑图像语义分割算法[J]. 现代计算机, 2023, 29(2): 23-31.
LUO Chongxing, SHI Mingyuan, WANG Zhengyong, et al. An improved cuttings image semantic segmentation algorithm based on DeepLabV3+[J]. Modern Computer, 2023, 29(2): 23-31. [本文引用:1]
[11] SINGH H, LI C X, CHENG P, et al. Automated real-time formation evaluation from cuttings and drilling data analysis: State of the art[J]. Advances in Geo-Energy Research, 2023, 8(1): 19-36. [本文引用:1]
[12] DONG Z K, YU H W, TIAN G S, et al. Data visualization strategy for predicting rock mass classification with TBM key rock-fragmentation parameters[J]. Transportation Geotechnics, 2025, 51: 101489. [本文引用:1]
[13] LIU C C, WANG Y J, PING X L, et al. Optimization of fracturing stages and perforation clusters based on RoqSCAN cutting mineralogical analysis[J]. Journal of Physics: Conference series, 2025, 3048: 012116. [本文引用:1]
[14] ZHAO Z M, GAO W Y, CHANG J B, et al. An analysis of Deep Neural Network model in recognition of mud cuttings image for practical applications[J]. Journal of Petroleum Science and Technology, 2022, 12(4): 53-61. [本文引用:1]
[15] ELEN A, DÖNMEZ E. Histogram-based global thresholding method for image binarization[J]. Optik, 2024, 306: 171814. [本文引用:1]
[16] ZULUNOV R, SOLIYEV B, KAYUMOV A, et al. Detecting mobile objects with ai using edge detection and background subtraction techniques[C]//E3S Web of Conferences, International Conference on Green Energy: Intelligent Transport Systems-Clean Energy Transitions (GreenEnergy 2023). Les Ulis: EDP Sciences, 2024: 03004. [本文引用:1]
[17] ZHANG Y H, LI Z B, HAN X. Rock cutting image recognition based on color and texture feature fusion[J]. Journal of Physics: Conference Series, 2024, 2901: 012026. [本文引用:1]
[18] ZHAI H Y, CHEN M, YANG X X, et al. Multi-scale HSV color feature embedding for high-fidelity NIR-to-RGB spectrum translation[C]//Proceedings of the 2024 IEEE International Conference on Systems, Man, and Cybernetics (SMC). Kuching, Malaysia: IEEE, 2024: 3287-3292. [本文引用:1]
[19] 雄贝(上海)科技有限公司. 一种基于矿物成分的岩屑岩性激光识别方法: 201911327098. 0[P]. 2019-12-20.
Xiongbei Tech (Shanghai) Co. , Ltd. A mineral composition-based laser identification method for cuttings lithology: 201911327098. 0[P]. 2019-12-20. [本文引用:1]
[20] YANG Z Q, ZHAO Q H, GAN J, et al. Damage evolution characteristics of siliceous slate with varying initial water content during freeze-thaw cycles[J]. Science of The Total Environment, 2024, 950: 175200. [本文引用:1]
[21] JIANG B R, YE W J, ZHAO Y Y. Nanopore structure and fractal characteristics of sand stone after high-temperature action[J]. Environmental Earth Sciences, 2024, 83(10): 300. [本文引用:1]
[22] JUN J, LIANG W. Investigation of the pore structure characteristics and fluid components of Quaternary mudstone biogas reservoirs: A case study of the Qaidam Basin in China[J]. Scientific Reports, 2024, 14(1): 26512. [本文引用:1]