资源描述:
2 0 1 6年 1 1月 N o v e m b e r 2 0 1 6 岩 矿 测 试 R O C KA N DM I N E R A LA N A L Y S I S V o l . 3 5 ,N o . 6 5 9 5~ 6 0 2 收稿日期 2 0 1 6- 0 5- 2 3 ;修回日期 2 0 1 6- 0 8- 1 9 ;接受日期 2 0 1 6- 0 8- 3 0 基金项目中国科学院战略性先导科技专项( B类) “ 页岩三维成像实验技术和数据获取技术” ( X D B 1 0 0 2 0 1 0 2 ) ; 国家杰出 青年科学基金资助项目( 4 1 3 2 5 0 1 6 ) 作者简介王羽, 助理研究员, 从事同步辐射技术在地质地震领域的应用研究。E - m a i l y u w a n g @s i n a p . a c . c n 。 通讯作者汪丽华, 副研究员, 从事同步辐射技术在地质考古领域的应用研究。E - m a i l l h w a n g @s i n a p . a c . c n 。 王羽,金婵,汪丽华, 等. 基于 S E M图像灰度水平的页岩孔隙分割方法研究[ J ] . 岩矿测试, 2 0 1 6 , 3 5 ( 6 ) 5 9 5- 6 0 2 . WA N GY u ,J I NC h a n ,WA N GL i - h u a ,e t a l . P o r e S e g m e n t a t i o nM e t h o d s B a s e do nG r a y S c a l e o f S c a n n i n g E l e c t r o nM i c r o s c o p y I m a g e s [ J ] . R o c ka n dM i n e r a l A n a l y s i s , 2 0 1 6 , 3 5 ( 6 ) 5 9 5- 6 0 2 . 【 D O I 1 0 . 1 5 8 9 8 / j . c n k i . 1 1- 2 1 3 1 / t d . 2 0 1 6 . 0 6 . 0 0 5 】 ★ 基于 S E M图像灰度水平的页岩孔隙分割方法研究 王 羽1 , 2,金 婵1 , 2,汪丽华1 , 2 *,王建强1 , 2,姜 政1 , 2,王彦飞3 ( 1 . 中国科学院微观界面物理与探测重点实验室,上海 2 0 1 8 0 0 ; 2 . 中国科学院上海应用物理研究所上海光源,上海 2 0 1 2 0 4 ; 3 . 中国科学院地质与地球物理研究所,北京 1 0 0 0 2 9 ) 摘要微观孔隙结构是研究页岩气吸附运移机制和建立地质模型的基础, 氩离子抛光 -扫描电子显微镜 ( S E M) 技术是开展此项研究的主要实验方法, 但已有的研究大多是关注页岩孔隙分类, 较少从定量角度表 征其特征。为开展页岩微观孔隙结构定量研究, 提高孔隙分割质量, 本研究分别利用边缘检测分割法、 流域 分割法、 手动和自动阈值分割法对页岩无机孔和有机孔二次电子图像进行分割实验, 对比不同方法的分割效 果。结果表明, 通过选取合适的分割阈值, 基于 S E M图像的手动阈值分割法能够表征 1n m以上的孔隙, 准确地识别有机质与脆性矿物边缘、 孔隙与有机质边缘, 使得页岩孔隙提取结果趋近于真实, 能更有效地对 页岩孔隙结构进行定量分析。 关键词富有机质页岩;S E M图像;孔隙分割;阈值法 中图分类号P 6 1 9 . 2 2 7 ;P 5 7 5 . 2文献标识码A 微观孔隙结构是研究页岩气吸附运移机制和建 立地质模型的基础[ 1 - 2 ]。氩离子抛光 - 扫描电子显 微镜( A r I o nM i l l i n g -S E M) 技术, 聚焦离子束扫描 电子显微镜( F I B-S E M) 技术和 X射线断层扫描 ( X- C T ) 技术是当前页岩微观孔隙结构研究的主要 测试技术, 不仅能够直观地定性反映页岩的孔隙类 型、 大小和形状, 同时结合相关软件能够定量提取孔 隙结构参数( 孔径分布、 孔隙度、 配位数、 连通率 等) , 实现页岩微纳米级孔隙结构的定性描述与定 量分析相结合[ 3 - 8 ]。采用上述方法进行页岩孔隙结 构定量分析的前提均涉及合理的图像分割, 如何提 高孔隙识别效率是进行页岩孔隙结构图像分析和岩 石物理参数提取的关键。目前边缘检测分割法、 流 域分割法和阈值分割法是岩石孔隙提取最常用的方 法。刘娜等[ 9 ]提出基于形态重构的流域分割法, 能 够对岩石颗粒图像进行较好分割; H o u b e n等[ 1 0 ]分 别利用不同图像分割方法, 对 O p a l i n u s 黏土岩进行 孔隙分割实验, 得出阈值法与边缘检测法相结合是 进行孔隙分割的有效手段; J i 等[ 1 1 ]采用 O t s u阈值 算法对 M a j e l l a 灰岩 C T图像进行分割, 实现了该灰 岩宏孔隙与微孔隙的定量描述。与上述分割对象相 比较, 页岩孔隙类型众多, 形态复杂, 且其边缘有一 定范围, 不同图像分割方法的提取结果可能存在较 大误差, 有必要针对不同分割方法提取页岩孔隙的 有效性展开研究。 高分辨 S E M二次电子图像可以清楚显示页岩 孔隙、 无机矿物基质和有机质的灰度差异[ 1 0 ]。利用 图像处理技术可将灰度值相对低的孔隙从 S E M图 像中分割出来, 即能够精确识别页岩有效孔隙并提 取其特征参数。因此, 本文分别采用边缘检测分割 法、 流域分割法、 手动和自动阈值分割法对页岩 595 ChaoXing S E M二次电子图像中的无机孔和有机孔开展孔隙 识别实验, 对比不同方法的分割效果, 探究适用于页 岩孔隙识别的最佳方法, 以期为页岩二维与三维孔 喉结构重建提供数据参考。 1 实验部分 1 . 1 样品采集与处理 页岩样品采自四川盆地中部石柱县打风坳剖面 志留系龙马溪组, X射线衍射实验表明其主要矿物包 含石英、 长石、 方解石、 白云石、 黄铁矿和黏土矿物。 取页岩薄片样品进行初步机械抛光与进一步氩离子 抛光获得高质量表面后, 应用 S E M技术发现无机孔 和有机孔两种最主要的孔隙类型。无机孔发育在矿 物基质中, 孔径从数百纳米至数微米, 形状复杂多样; 有机孔主要发育在有机质内部, 呈片麻状或蜂窝状。 1 . 2 孔隙识别原理 二次电子图像不仅与样品形貌相关, 还与样品 组成元素的原子序数关系密切, 能够清楚地显示孔 隙与矿物基质在灰度值上的差异( 图 1 ) 。灰度值分 布图表明黑色凹陷区域为孔隙, 由于电子聚集效应, 其边缘通常较为明亮。 ①电子聚集造成的明亮边缘; ②形貌起伏所造成的假孔隙。 图 1 基于 S E M图像的孔隙识别 F i g . 1 P o r ei d e n t i f i c a t i o nu s i n gs e c o n d a r yS E Mi m a g e 1 . 3 孔隙分割方法 原始二次电子图像经过裁剪后, 首先进行背景 灰度扣除和平滑, 以消除图像亮度不均匀性和噪声; 不同的图像分割方法在 V C++开发环境下使用面 向对象的方法进行开发。 1 . 3 . 1 边缘检测分割方法 边缘检测分割方法是利用一阶和二阶偏导数识 别图像中的不连续点( 梯度方法) [ 1 2 ], 其定义如下 G ( x , y )= ∂ f ∂ x , ∂ f ∂ [] y T , 其中 f ( x , y ) 为灰度图像。本方 法主要用以标识矿物颗粒或孔隙的边缘。本文采用 s o b e l 算子能够加强中心像素上下左右 4个方向像 素的权重, 图像边缘亮度较大, 并能抑制孤立噪声。 1 . 3 . 2 流域分割方法 基于拓扑理论的流域分割方法( 分水岭分割方 法) 对灰度敏感, 分割得到的孔隙边界定位准确且 连续闭合。其中心思想是将图像视为地质学中的拓 扑地貌, 图像中每点像素的灰度值表示该点的海拔 高度, 每一个局部极小值及其影响区域称为集水盆, 而集水盆的边界视为分水岭。流域分割方法的迭代 标注计算过程分为排序和淹没两步, 最终得到输入 图像的集水盆图像, 集水盆之间的边界点即为分水 岭[ 1 3 ]。具体到页岩二次电子 S E M图像中, 可将孔 隙边界看作分水岭。进行孔隙分割时, 先利用 s o b e l 算子将原始图像转化成梯度图像, 把梯度图像作为 输入图像, 以获取孔隙的边缘信息。 1 . 3 . 3 阈值分割方法 阈值分割方法是基于图像中要提取的目标与背 景的灰度差异性, 依据像素的灰度值进行孔隙的提 取[ 1 4 ], 该算法一般描述为 T=T [ x , y , f ( x , y ) , g ( x , y ) ] , 其中 f ( x , y ) 代表像素点( x , y ) 处的灰度值, g ( x , y ) 代表该点邻域的某种局部性质。在页岩孔 隙分割中, 假若页岩 S E M图像像素的灰度 T分布范 围 0 ( 黑色)~ 2 5 5 ( 白色) , 根据视觉判断确定一个 临界灰度值 T p, 当像素的灰度值小于 Tp时, 确定为 孔隙, 反之则确定为矿物基质。 2 结果与讨论 2 . 1 不同方法孔隙分割效果 分别采用边缘检测分割法、 流域分割法、 自动阈 值分割法和手动阈值分割法对页岩气储集空间无机 孔和有机孔二次电子扫描图像进行孔隙分割实验, 实验对比结果如图 2和图 3所示。 图 2 a 为页岩无机粒间孔的 S E M二次电子图 像, 边缘检测法能够识别无机孔隙, 但是受形貌起伏 695 第 6期 岩 矿 测 试 h t t p ∥w w w . y k c s . a c . c n 2 0 1 6年 ChaoXing 图 2 无机孔的不同孔隙分割方法效果对比图 F i g . 2 C o m p a r i s o n s o f m i n e r a l m a t r i xp o r es e g m e n t a t i o ne f f e c t s w i t hd i f f e r e n t s e g m e n t a t i o nm e t h o d s 影响, 误差较大( 图 2 b ) 。流域分割法能够有效地识 别大孔, 但是忽视了狭缝和孤立的小孔( 图 2 c ) 。阈 值法能够较为精确地识别大多数无机孔隙, 但是相 较而言, 自动阈值法边缘信息不明确( 图 2 d , e ) 。 图 3 a 为页岩有机孔的 S E M二次电子图像, 边 缘检测分割法( 图 3 b ) 和流域分割法( 图 3 c ) 均能够 识别有机质与孔隙, 然而边缘检测分割法受孔隙形 貌起伏影响, 造成较大误差。流域分割法虽然定位 相对准确, 但是边缘信息不明确。相较而言, 阈值分 割法对灰度十分敏感, 能在一定程度上改善边缘检 测分割法和流域分割法中丢失微弱边缘的缺点( 图 3 d , e ) , 适用于进行页岩孔隙提取。自动阈值分割 法操作便捷, 但是其结果精度不足以满足孔隙定量 分析的需求。手动阈值分割法通过自主选择合理的 分割阈值, 能够使结果趋近于真实的孔隙度。 2 . 2 阈值分割法在页岩孔隙分割中的应用示例 作为一种典型的多相致密储层, 页岩中发育有 多种不同类型的孔隙, 如黏土矿物层间孔、 粒间孔、 795 第 6期王羽, 等 基于 S E M图像灰度水平的页岩孔隙分割方法研究第 3 5卷 ChaoXing 图 3 有机孔的不同孔隙分割方法效果对比图 F i g . 3 C o m p a r i s o n s o f o r g a n i cm a t t e r p o r es e g m e n t a t i o ne f f e c t s w i t hd i f f e r e n t s e g m e n t a t i o nm e t h o d s 895 第 6期 岩 矿 测 试 h t t p ∥w w w . y k c s . a c . c n 2 0 1 6年 ChaoXing 溶蚀孔、 晶间孔和有机孔等[ 1 5 - 1 6 ]。页岩气储层性能 的优劣在很大程度上取决于微纳孔隙的大小、 分布、 形貌和连通性等。利用手动阈值分割法对页岩中普 遍发育的孔隙结构进行孔隙提取, 分割阈值根据统 计学原理, 利用高分辨扫描电子图像获取。 图 4 a 是页岩微观结构的高分辨二次电子图像, 根据灰度差异可以清楚地区分黄铁矿、 脆性矿物、 有 机质和孔隙四大类组分。将图像放大至 2 5 0 0倍, 可 清晰地观测到边界阈值出现梯度变化, 利用 M a t l a b 中的 i m t o o l 工具依次读取两种不同组分的边界阈值 3 0 0组, 该组数据的平均值视为分割阈值。如图 4 b 和 c 所示, 分别展示了有机质与孔隙、 有机质与脆性 矿物的边界。阈值读取结果表明有机孔的阈值范围 为 2 2~ 2 7 , 有机质的阈值范围为 4 8~ 5 3 , 从有机孔 过渡到有机质, 阈值呈现出梯度上升, 其范围基本涵 盖 2 8~ 4 9 , 依据统计学原理计算将 3 2视为有机孔 与有机质的分割阈值( 图 4 b ) 。同理, 图 4 c 中根据 阈值统计计算结果, 可将 7 1视为有机质与无机矿物 基质的分割阈值。 图 5为应用手动阈值分割法对页岩不同类型孔 隙进行孔隙提取的结果, 其中黑色区域部分代表孔 隙。对粒间孔、 层间孔、 解理孔、 溶蚀孔和微裂缝等 无机孔而言, 当分割阈值在合理范围微小变动时, 分 割结果相近, 对孔隙定量分析结果影响较小。以 图 5 e 所示的微裂缝为例, 当分割阈值分别设为 2 5 、 3 0与 3 5时, 面孔率计算结果依次为 7 . 5 %、 8 . 2 %与 8 . 5 %, 差异较小。有机孔与有机质灰度相似, 分割 阈值的微弱变动对分割结果影响相对大, 根据不同 分割阈值进行孔隙定量分析时会产生较大误差。 图 5 g 是有机孔孔隙分割结果, 结合形貌学特征, 该 分割阈值设为 1 8 , 面孔率计算结果为 2 5 %。但是当 分割阈值分别设为 2 0和 2 5时, 面孔率计算结算分 别为 2 7 . 6 %与 2 9 . 3 %, 波动较大。上述结果表明, 应用手动阈值分割法能够对页岩不同类型孔隙结构 进行有效的孔隙提取, 其中合理选择分割阈值是孔 隙提取的关键。 2 . 3 手动阈值法对 C T断层扫描定量分析页岩 孔隙的启示 手动阈值法能够准确地识别有机质与脆性矿物 边缘、 孔隙与有机质边缘, 更加适合页岩孔隙提取。 基于 S E M图像的手动阈值分割法不仅能够有效地 对页岩孔隙结构进行定量分析, 亦可为 C T图像的 ( a ) 页岩微观结构的二次电子图像; ( b ) 有机质与孔隙的阈值边界; ( c ) 有机质与脆性矿物的阈值边界。( b ) 和( c ) 中红色十字和白色 数字分别表示像素点位置及该像素点的阈值。 图 4 页岩孔隙分割中分割阈值的获取方法示例 F i g . 4 I m a g e ss h o w i n gt h r e s h o l dv a l u ea c q u i s i t i o nm e t h o di n s h a l ep o r es e g m e n t a t i o n s 孔隙分割提供参考。 C T无损扫描技术被广泛应用于岩石孔隙结构 研究。随着测试技术的改进, C T图像质量虽然有了 一定提高, 但是生成高质量的纳米级 C T图像依然 存在难点, 使得微纳米级孔隙的精确提取十分困难。 作为典型的非常规致密多相储层, 页岩孔隙三维结 构的定量化分析一直是研究的重点与难点。国内外 研究者对页岩三维孔隙结构的研究多采用 F I B- S E M和 C T无损扫描技术, 前者主要基于 S E M图像 进行孔隙和不同矿物的分割, 后者则基于 X射线切 片进行图像分割, 两者的灰度深浅均与样品组成元 素的原子序数线性相关。因此, 可通过 S E M图像确 定不同组分的分割阈值, 并计算分割阈值的比例关 系, 进而以该比例关系为参考, 对页岩 C T扫描图像 进行孔隙提取, 以获取更为准确的定量结果, 为页岩 气储量评估和渗流模拟提供基础数据支持。 995 第 6期王羽, 等 基于 S E M图像灰度水平的页岩孔隙分割方法研究第 3 5卷 ChaoXing 图 5 手动阈值分割法对页岩不同类型孔隙进行孔隙提取的结果 F i g . 5 C o m p a r i s o n s o f o r g a n i cm a t t e r p o r es e g m e n t a t i o ne f f e c t s d e v e l o p i n gi ns h a l ew i t hd i f f e r e n t s e g m e n t a t i o nm e t h o d s 3 结论 本研究分别利用边缘检测分割法、 流域分割法 和手动或自动阈值分割法进行页岩孔隙结构提取实 验, 结果表明手动阈值分割法通过自主选择合理的 分割阈值, 能够有效区分无机孔、 矿物、 有机质和有 机孔, 使得页岩孔隙提取结果趋近于真实, 为页岩微 观孔隙结构的定量化研究提供了新的思路与参考。 同时, 基于扫描电子图像获取的分割阈值比例亦可 为 C T图像的孔隙分割提供参考。依据统计学原理 确定分割阈值的方法操作简单, 但是分割阈值的微 弱变动对孔隙分割影响大, 尤其是有机孔。因而采 取阈值法进行孔隙分割时还须充分考虑孔隙的真实 空间形态, 合理选择分割阈值。 4 参考文献 [ 1 ] 邹才能, 朱如凯, 白斌, 等. 中国油气储层中纳米孔首 次发现及其科学价值[ J ] . 岩石学报, 2 0 1 1 , 2 7 ( 6 ) 1 8 5 7- 1 8 6 4 . Z o uCN , Z h uRK , B a i B , e t a l . F i r s t D i s c o v e r y o f N a n o - p o r eT h r o a t i nO i l a n dG a sR e s e r v o i ri nC h i n aa n dI t s S c i e n t i f i cV a l u e [ J ] . A c t aP e t r o l o g i c aS i n i c a , 2 0 1 1 , 2 7 ( 6 ) 1 8 5 7- 1 8 6 4 . [ 2 ] C l a r k s o nCR , J e n s e nJL , B l a s i n g a m eTA . R e s e r v o i r E n g i n e e r i n gf o r U n c o n v e n t i o n a l G a s R e s e r v o i r s Wh a t D o We H a v e t o C o n s i d e r [C] / / S P E P a p e r 1 4 5 0 8 0 P r e s e n t e da tt h eS o c i e t yo fP e t r o l e u m E n g i n e e r sN o r t h A m e r i c a nU n c o n v e n t i o n a l G a sC o n f e r e n c ea n d E x h i b i t i o n . Wo o d l a n d s , T e x a s , 2 0 1 1 . [ 3 ] L o u c k sRG , R e e dRM , R u p p e l SC , e t a l . M o r p h o l o g y , G e n e s i s , a n d D i s t r i b u t i o n o fN a n o m e t e r - s c a l eP o r e si n S i l i c e o u s M u d s t o n e s o f t h e M i s s i s s i p p i a nB a r n e t t S h a l e [ J ] . J o u r n a l o f S e d i m e n t a r y R e s e a r c h , 2 0 0 9 , 7 9 ( 1 2 ) 8 4 8 - 8 6 1 . [ 4 ] K l a v e r J , D e s b o i s G , L i t t k e R , e t a l . B I B - S E MC h a r a c t e r i - z a t i o no fP o r eS p a c eM o r p h o l o g ya n dD i s t r i b u t i o ni n P o s t m a t u r et oO v e r m a t u r eS a m p l e sf r o mt h eH a y n e s v i l l e a n dB o s s i e r S h a l e s [ J ] . M a r i n ea n dP e t r o l e u mG e o l o g y , 2 0 1 5 , 5 9 4 5 1- 4 6 6 . [ 5 ] 王羽, 金婵, 汪丽华, 等. 应用氩离子抛光 -扫描电镜 方法研究四川九老洞组页岩微观孔隙特征[ J ] . 岩矿 测试, 2 0 1 5 , 3 4 ( 3 ) 2 7 8- 2 8 5 . Wa n gY , J i nC , Wa n gLH , e t a l . C h a r a c t e r i z a t i o no f P o r e S t r u c t u r e s o f J i u l a o d o n gF o r m a t i o nS h a l ei nt h eS i c h u a n B a s i nb y S E Mw i t hA r - i o nM i l l i n g [ J ] . R o c ka n dM i n e r a l A n a l y s i s , 2 0 1 5 , 3 4 ( 3 ) 2 7 8- 2 8 5 . [ 6 ] C u r t i sM E , S o n d e r g e l dC H , A m b r o s eR J , e ta l . M i c r o s t r u c t u r a l I n v e s t i g a t i o no f G a sS h a l e si nT w oa n d 006 第 6期 岩 矿 测 试 h t t p ∥w w w . y k c s . a c . c n 2 0 1 6年 ChaoXing T h r e e D i m e n s i o n sU s i n g N a n o m e t e r - s c a l e R e s o l u t i o n I m a g i n g [ J ] . A A P GB u l l e t i n , 2 0 1 2 , 9 6 ( 4 ) 6 6 5- 6 7 7 . [ 7 ] D e w e r sTA , H e a t hJ , E w yR , e t a l . T h r e e - d i m e n s i o n a l P o r eN e t w o r k sa n dT r a n s p o r t P r o p e r t i e so f aS h a l eG a s F o r m a t i o nD e t e r m i n e df r o m F o c u s e dI o nB e a m S e r i a l I m a g i n g [ J ] . I n t e r n a t i o n a l J o u r n a l o f O i l G a sa n dC o a l T e c h n o l o g y , 2 0 1 2 , 5 ( 2- 3 ) 2 2 9- 2 4 8 . [ 8 ] Wa n gY , P uJ , Wa n gLH , e ta l . C h a r a c t e r i z a t i o no f T y p i c a l 3 DP o r eN e t w o r k s o f J i u l a o d o n gF o r m a t i o nS h a l e U s i n gN a n o - t r a n s m i s s i o nX - r a yM i c r o s c o p y [ J ] . F u e l , 2 0 1 6 , 1 7 0 8 4- 9 1 . [ 9 ] 刘娜, 郭连军, 赵楠楠. 基于形态重构的分水岭岩石图 像分割方法[ J ] . 辽宁科技大学学报, 2 0 1 0 , 3 3 ( 5 ) 4 9 5- 4 9 8 . L i uN , G u oLJ , Z h a oNN . S e g m e n t a t i o nA l g o r i t h m o f R o c kI m a g ew i t hM o r p h o l o g i c a lR e c o n s t r u c t i o n[ J ] . J o u r n a lo f U n i v e r s i t y o f S c i e n c e a n d T e c h n o l o g y L i a o n i n g , 2 0 1 0 , 3 3 ( 5 ) 4 9 5- 4 9 8 . [ 1 0 ] H o u b e nME , D e s b o i s G , U r a i J L . P o r eM o r p h o l o g ya n d D i s t r i b u t i o ni nt h eS h a l yF a c i e s o f O p a l i n u s C l a y( M o n t T e r r i , S w i t z e r l a n d ) I n s i g h t s f r o mR e p r e s e n t a t i v e 2 DB I B - S E MI n v e s t i g a t i o n s o nm mt o n mS c a l e [ J ] . A p p l i e dC l a y S c i e n c e , 2 0 1 3 , 7 1 8 2- 9 7 . [ 1 1 ] J i YT , H a l l SA , B a u dP , e t a l . C h a r a c t e r i z a t i o no f P o r e S t r u c t u r ea n dS t r a i nL o c a l i z a t i o ni nM a j e l l a L i m e s t o n e b y X - r a y C o m p u t e d T o m o g r a p h y a n d D i g i t a l I m a g e C o r r e l a t i o n [ J ] . G e o p h y s i c a l J o u r n a l I n t e r n a t i o n a l , 2 0 1 5 , 2 0 0 ( 2 ) 6 9 9- 7 1 7 . [ 1 2 ] 何斌, 马天予, 王运坚等编著. V i s u a 1 C++数字图像 处理( 第二版) [ M] . 北京 人民邮电出版社, 2 0 0 2 3 9 4- 3 9 8 . H eB , M aTY , Wa n gYJ , e t a l . V i s u a l CD i g i t a l I m a g e P r o c e s s i n g( S e c o n dE d i t i o n )[ M] .B e i j i n g P o s t s& T e l e c o mP r e s s , 2 0 0 2 3 9 4- 3 9 8 . [ 1 3 ] G o n z a l e zRC , Wo o d sRE . D i g i t a lI m a g eP r o c e s s i n g ( S e c o n dE d i t i o n )[ M] . B e i j i n g P u b l i s h i n gH o u s eo f E l e c t r o n i c s I n d u s t r y , 2 0 0 5 . [ 1 4 ] N i k h i l RP , S a n k a r KP . AR e v i e wo nI m a g eS e g m e n t - a t i o n T e c h n i q u e s[J ] .T h e J o u r n a lo ft h e P a t t e r n R e c o g n i t i o nS o c i e t y , 1 9 9 3 , 2 6 ( 9 ) 1 2 7 7- 1 2 9 4 . [ 1 5 ] L o u c k s RG , R e e dRM, R u p p e l SC , e t a l . S p e c t r u mo f P o r eT y p e sf o rM a t r i x - r e l a t e dM u dP o r e s [ J ] . A A P G B u l l e t i n , 2 0 1 2 , 9 6 ( 6 ) 1 0 7 1- 1 0 9 8 . [ 1 6 ] 于炳松. 页岩气储层孔隙分类与表征[ J ] . 地学前缘, 2 0 1 3 , 2 0 ( 4 ) 2 1 1- 2 2 0 . Y uBS . C l a s s i f i c a t i o na n dC h a r a c t e r i z a t i o no ft h eG a s S h a l eP o r eS y s t e m [ J ] . E a r t hS c i e n c eF r o n t i e r s , 2 0 1 3 , 2 0 ( 4 ) 2 1 1- 2 2 0 . 106 第 6期王羽, 等 基于 S E M图像灰度水平的页岩孔隙分割方法研究第 3 5卷 ChaoXing P o r eS e g me n t a t i o nMe t h o d sB a s e do nG r a yS c a l eo fS c a n n i n gE l e c t r o n Mi c r o s c o p yI ma g e s W A N GY u 1 , 2 , J I NC h a n 1 , 2 , W A N GL i - h u a 1 , 2 *, W A N GJ i a n - q i a n g 1 , 2 , J I A N GZ h e n g 1 , 2 , W A N GY a n - f e i 3 ( 1 . K e yL a b o r a t o r y o f I n t e r f a c i a l P h y s i c s a n dT e c h n o l o g y ,S h a n g h a i I n s t i t u t e o f A p p l i e dP h y s i c s ,C h i n e s e A c a d e m y o f S c i e n c e s ,S h a n g h a i 2 0 1 8 0 0 ,C h i n a ; 2 . S h a n g h a i S y n c h r o t r o nR a d i a t i o nF a c i l i t y ,S h a n g h a i I n s t i t u t e o f A p p l i e dP h y s i c s ,C h i n e s e A c a d e m y o f S c i e n c e s , S h a n g h a i 2 0 1 2 0 4 ,C h i n a ; 3 . I n s t i t u t eo f G e o l o g ya n dG e o p h y s i c s ,C h i n e s eA c a d e m yo f S c i e n c e s ,B e i j i n g 1 0 0 0 2 9 ,C h i n a ) H i g h l i g h t s P o r es e g m e n t a t i o nw a s p e r f o r m e db a s e do ng r a ys c a l eo f S E Mi m a g e s . S e g m e n t a t i o n r e s u l t sw e r ec o m p a r a t i v e l ya n a l y z e d u s i n ge d g ed e t e c t i o n , w a t e r s h e d , a u t oa n d m a n u a l t h r e s h o l d i n gm e t h o d s . T h em a n u a l t h r e s h o l d i n gm e t h o dw a s m o r es u i t a b l ef o r s h a l ep o r es e g m e n t a t i o n . A b s t r a c t M i c r o s c o p i cp o r es t r u c t u r e so f s h a l ea r et h eb a s i sf o ri n v e s t i g a t i n gt h e a d s o r p t i o n a n d m i g r a t i o n m e c h a n i s m o f s h a l eg a s a n db u i l d i n gag e o l o g i c a l m o d e l . A ri o n m i l l i n gc o m b i n e d w i t h S c a n n i n g E l e c t r o nM i c r o s c o p y( S E M)i st h em a i n t e c h n i q u ef o ra n a l y z i n gm i c r o s c o p i cp o r e s t r u c t u r e so fs h a l e . H o w e
展开阅读全文