數字巖心EDS彩色圖像分割

謝 冰

（長江大學地球物理與石油資源學院，湖北 武漢430100）

0 引言

隨著我國對油氣資源的需求量不斷增加，非常規油氣的勘探與開發越來越備受關注。巖石物理研究在油氣藏評價中占據著重要地位。針對頁巖、致密砂巖、碳酸鹽巖等非常規油氣復雜儲層，傳統的巖石物理實驗遇到了諸多困難[1]。數字巖石物理的出現很好地解決了這一難題。數字巖石物理通過建立三維數字巖心并借助數值算法對巖石的聲學特性、電性、核磁共振特性及滲流特性等進行數值模擬實驗。數字巖石物理實驗與傳統的巖石物理實驗相比有很多優點[1]。三維數字巖心重建方法包括物理方法和數值重建方法兩大類。物理方法借助CT、SEM、FIB、EDS等儀器獲取巖心的多幅二維圖像直接重建三維數字巖心。數值重建方法則借助少量的巖心二維圖像，通過圖像分析提取建模信息，采用各種重建算法[2]建立三維數字巖心。頁巖的礦物成分及孔隙分布相當復雜，通過EDS圖像可以有效地分析頁巖的各種礦物及孔隙分布。本文通過對Eagle Ford頁巖EDS彩色圖像進行分割，得到了各種礦物及孔隙的二值圖像，為后期的三維數字巖心重建及巖石物理數值模擬奠定了數據基礎。

1 X射線能譜儀（EDS）

X射線能譜儀（EDS）是依據電子與物質的相互作用。當一束高能的入射電子轟擊物質表面時，被激發的區域將產生二次電子、俄歇電子、特征x射線和連續譜X射線、背散射電子、透射電子，以及在可見、紫外、紅外光區域產生的電磁輻射。原則上講，利用電子和物質的相互作用，可以獲取被測樣品本身的各種物理、化學性質的信息，如形貌、組成、晶體結構、電子結構和內部電場或磁場等等。EDS根據上述不同信息產生的機理，對二次電子、背散射電子的采集，可得到有關物質微觀形貌的信息，對x射線的采集，可得到物質化學成分的信息。

1.1 基本原理[3]

1.2 分析技術

EDS分析包括X射線的測量、空間分辨率、峰/背比(P/B)、定性分析、定量分析、元素的面分布分析方法。電子束只打到試樣上一點，得到這一點的X射線譜的分析方法是點分析方法。與此不同的是，用掃描像觀察裝置，使電子束在試樣上做二維掃描，測量特征X射線的強度，使與這個強度對應的亮度變化與掃描信號同步在陰極射線管CRT上顯示出來，就得到特征X射線強度的二維分布的像。這種觀察方法稱為元素的面分布分析方法，它是一種測量元素二維分布的非常方便的方法。利用特定的方法可以將元素二維分布圖像轉化為礦物及孔隙分布圖像。

2 EDS圖像分割

2.1 分割方法

圖1是Eagle Ford頁巖氣儲層某塊頁巖巖心的EDS彩色圖像[4]，圖中各種顏色對應的礦物名稱以及相應的RGB值如表1所示。本文是基于RGB色彩空間對EDS圖像進行分割。EDS彩色圖像中每個像素對應一個RGB值。首先利用Matlab中的imtool函數得到每種礦物對應顏色的RGB值；然后利用三維歐式距離空間對每個像素的RGB值進行聚類，得到聚類后的EDS彩色圖像（圖2）；最后根據聚類后的EDS彩色圖像進行二值分割。

圖1 Eagle Ford頁巖EDS圖像

圖2 聚類后的EDS圖像

表1 圖1中各種顏色對應的礦物名稱和GRB值

2.2 分割結果

圖3

圖3（a）-（g）分別是對不同礦物以及孔隙進行二值分割后得到的二值圖像。由圖可知，基于RGB色彩空間的歐式距離空間聚類并進行二值分割的方法，除了有機質，其它礦物及孔隙的二值分割效果很好。有機質的二值分割主要是由于不同顏色礦物接觸邊界的影響。

2.3 源程序

參考http：//user.qzone.qq.com/1057206671/2。

3 結論

本文提出的基于RGB色彩空間的歐式距離空間聚類并進行二值分割的方法能夠有效地對數字巖心EDS彩色圖像進行分割，從而得到各種礦物及孔隙的二值分割圖像，為三維數字巖心重建提供數據基礎。但是該方法存在的缺陷在于不能對不同顏色礦物接觸邊界處存在的過渡色進行準確的聚類，從而會影響最終的圖像分割結果。建議可以從HSV、L*a*b*等其它色彩空間對數字巖心EDS彩色圖像進行二值分割。

［1］孫建孟,姜黎明,等.數字巖心技術測井應用與展望[J].測井技術,2012,36(1)：1-7.

［2］劉學鋒,張偉偉,等.三維數字巖心建模方法綜述[J].地球物理學進展,2013,28(6)：3066-3067.

［3］http：//www.mttlab.com/g/testing/morphology/sem_eds/index.html[OL].

［4］Joel D.Walls,Steven W.Sinclair.Eagle Ford Shale Reservoir Properties From Digital Rock Physics[J].EAGE,2011,29(6)：97-101.