两种图像分割算法在FMI成像资料中的应用.docx

两种图像分割算法在 FM I 成像资料中的应用张丽莉, 刘瑞林(江汉石油学院地球物理勘探系, 湖北 荆州 434102)[ 摘要 ] 为了从地层微电阻率扫描成像测井 (FM I) 图像中分割出主要反映裂缝、孔洞的子图像 两种图像分割算法在 FM I 成像资料中的应用 张丽莉, 刘瑞林 (江汉石油学院地球物理勘探系, 湖北 荆州 434102) [ 摘要 ] 为了从地层微电阻率扫描成像测井 (FM I) 图像中分割出主要反映裂缝、孔洞的子图像, 介绍了两 种图像分割算法: 一种是基于图像间模糊散度的阈值化算法, 另一种是基于过渡区的分割算法。并分别用 这两种分割算法对塔里木盆地中央隆起 T Z45 井 FM I 图像进行分割。研究结果表明, 这两种算法对 FM I 成 像测井资料获得了较好的分割效果。 [ 主题词 ] 地层微扫描测井; 图像处理; 模糊数学; (模糊散度) ; (过渡区) [ 中图分类号 ] P 631184 [ 文献标识码 ] A [ 文章编号 1000 9752 (1999) 04 0088 03 目前, 对地层微电阻率扫描成像测井 (FM I) 资料的应用主要局限于孔洞、裂缝的定性识别和地层倾 角解释等方面。对于 FM I 资料的图像处理和定量参数提取方面, 国内外都处于研究、试验阶段。为了从 FM I 资料中提取定量信息, 如面孔率、裂缝长度等参数, 一个重要的步骤就是要对 FM I 图像数据进行分 割, 以期从实际 FM I 资料中分离出反映孔洞、裂缝的子图像。FM I 图像分割的结果直接影响参数计算的 准确性 。然而, 有关图像分割的技术与方法主要针对光学扫描图像进行的1 。FM I 成像测井资料与光学 扫描图像不同。FM I 图像值是所测量的地质体的电导率值。在不同地质体 (如基岩与孔洞) 的边界处, 测 量值不是陡变的, 而是存在明显的过渡区。另外, 含孔洞、裂缝的地层的 FM I 成像测井图像可以认为由 目标 (孔洞等) 和背景 (基岩) 两类组成, 目标与背景的灰度值主要集中在相应的均值附近。因此, FM I 测井图像可以通过灰度特征进行分割。笔者根据 FM I 图像的这一特征, 应用两种图像分割算法 (一种是 基于图像间模糊散度的阈值化算法, 另一种是基于过渡区的分割算法) , 对 FM I 图像进行了试验。 1 两种图像分割算法的基本原理 111 图像间模糊散度的阈值化分割算法 图像间模糊散度的阈值化分割算法是一种采用模糊集合分别表达分割前后的图像, 通过最小模糊散 度准则2 实现图像分割中最优阈值提取的算法。 最小模糊散度准则: 选择最优阈值使其满足使原始图和分割图的信息距离量最小的条件, 即要求最 优阈值满足使原始图与分割图之间的模糊散度最小。 用 P 与 Q 分别表征分割前后的图像。若将 P 与 Q 定义成论域 S 上的两个模糊集合, 则 {u p (s) f s ∈ S } {u q (s) f s ∈ S } P = Q = 式中, u p (s) , u q (s) 分别是 P , Q 的隶属度函数, 满足 0 ≤ u p (s) , u q (s) ≤ 1。在图像分割的应用中, 论域S 取灰 度轴, 元素 s 取像素点的灰度。在一幅图像内, 目标与背景的灰度主要集中在目标类与背景类的均值附近, 即均值的隶属度最高 (为 1) , 类内像素的隶属度和该像素到均值的距离成反比, 而且类内 (外) 所有像素的 隶属度应该大于 (小于) 015。对未分割的原始图 G 而言, 满足上述思想的目标的隶属度函数 u 1 (x ; t, G ) 和 背景的隶属度函数 u 2 (x ; t, G ) 分别为 x - m 1 ( t) a 1 0 ≤ x ≤ t 1 - 2 C u 1 (x ; t, G ) = u 2 (x ; t, G ) = 1 - u 1 (x ; t, G ) (1) a 1 x - m 2 ( t) 0 x ≤L 2 C [ 作者简介] 张丽莉 ( 1975 ) , 女, 1997 江汉石油学院毕业, 在读硕士生, 现从事成像测井应用的研究工作。 刘瑞林, 钟广法, 周云才 1FM I、A R I 成像测井资料定量计算储层参数方法研究及应用 1 塔里木油田公司与江汉石油学院研究报告 119991 第 21 卷第 4 期张丽莉等: 两种图像分割算法在 FM I 成像资料中的应用·89·式中, t 为待选阈值; L 为灰度上界; a 为类内元素对均值的紧致度参数, 可描述诸如大约、接近、相似等