戴 宗 羅東紅 謝明英 唐 放 李海龍 涂志勇 王晨晨 陳 強

(1. 中海石油(中國)有限公司深圳分公司研究院 廣東深圳 518000； 2. 數(shù)巖科技(廈門)股份有限公司 福建廈門 361000)

在剩余油潛力較大的老油田再開采過程中，由于部分室內(nèi)實驗基礎(chǔ)數(shù)據(jù)時間久遠(yuǎn)參考意義不大；而巖心庫中僅存的儲層代表性巖心往往由于老化嚴(yán)重、剩余量少，不具備重新進(jìn)行物性參數(shù)、粒徑分析等室內(nèi)實驗測試的條件，進(jìn)而無法準(zhǔn)確體現(xiàn)儲層巖石物理特征，這對老油田挖潛帶來了較大的困難和挑戰(zhàn)。隨著計算機圖像處理技術(shù)、計算方法的發(fā)展，以及CT掃描等技術(shù)在巖石物理實驗中的應(yīng)用，可以構(gòu)建能夠反映儲層巖石真實孔隙空間結(jié)構(gòu)特征的數(shù)字巖心，為巖石微觀結(jié)構(gòu)特征和滲流模擬提供研究平臺[1-4]。數(shù)字巖心建模方法主要包括以CT掃描為代表的物理建模法和數(shù)值重建法[5-8]。常見的數(shù)值重建算法主要有隨機法和地質(zhì)過程模擬法兩大類。隨機重建法是一種借助少量的巖心平面圖像資料，通過圖像分析技術(shù)提取建模信息，然后應(yīng)用某種數(shù)值計算方法建立數(shù)字巖心的方法，常用的隨機算法主要有：完全隨機法、高斯場法、順序指示模擬法、模擬退火算法、多點統(tǒng)計法和馬爾科夫鏈-蒙特卡洛算法等。地質(zhì)過程模擬法通過對顆粒的沉積過程及后續(xù)的壓實作用、成巖作用進(jìn)行模擬建立數(shù)字巖心[9-10]。

對于不具備完整的三維空間結(jié)構(gòu)特征和二維孔隙空間特征的巖屑樣品，不能直接通過CT掃描或者隨機重建方法來構(gòu)建三維數(shù)字巖心，因此需要通過地質(zhì)過程模擬法來構(gòu)建代表性的三維數(shù)字巖心。開展地質(zhì)過程模擬前，需要首先測量巖心的粒度組成，在此基礎(chǔ)上模擬顆粒的沉積、壓實和成巖過程，構(gòu)建三維數(shù)字巖心模型。本文針對儲層代表性的巖屑樣品，通過掃描電鏡成像可生成MAPS，在此基礎(chǔ)上提出一種粒徑分析新方法獲取巖屑粒徑分布，并通過地質(zhì)過程模擬法構(gòu)建代表性的三維數(shù)字巖心，分析儲層微觀孔隙結(jié)構(gòu)特征，為老油田剩余油挖潛提供基礎(chǔ)數(shù)據(jù)。

1 基于3D距離圖的分水嶺算法

MAPS圖像是在選定區(qū)域內(nèi)排布掃描出一系列連續(xù)且邊緣重疊的大量高分辨率的小圖像，掃描完成后將這些小圖像進(jìn)行拼接進(jìn)而得到一張超高分辨率、超大面積的二維背散射電子圖像?；贛APS圖像通過分水嶺算法進(jìn)行分割，可以研究不同尺度下的巖石顆粒及孔隙之間的關(guān)聯(lián)關(guān)系。分水嶺算法是一種基于拓?fù)淅碚摰臄?shù)學(xué)形態(tài)學(xué)的分割方法，其基本思想是把圖像看作是測地學(xué)上的拓?fù)涞孛玻瑘D像中每一點像素的灰度值表示該點的海拔高度，每一個局部極小值及其影響區(qū)域稱為集水盆，而集水盆的邊界則形成分水嶺。

然而對于巖屑樣品，基于傳統(tǒng)的分水嶺算法對MAPS圖像分割不能有效地識別出不同的顆粒邊緣。本文通過計算圖像的3D距離圖，提出一種基于3D距離圖的分水嶺算法，利用該算法可以對巖屑MAPS圖像粒徑進(jìn)行精確分割和計算，具體步驟如下：根據(jù)照片中的特征灰度值或形態(tài)，對原始照片進(jìn)行特征分割，劃分出顆粒相、基質(zhì)相和孔隙相；對要統(tǒng)計的特征進(jìn)行標(biāo)記，如統(tǒng)計礦物顆粒時，需要對顆粒相進(jìn)行標(biāo)記；對特征屬性進(jìn)行計算，如顆粒大小(面積、體積、周長等參數(shù))；對特征參數(shù)進(jìn)行分組統(tǒng)計，如顆粒的等效直徑等[11-12]。下面分別對以上過程進(jìn)行闡述。

1.1 特征分割

圖1 基于改進(jìn)分水嶺算法的粒徑分割工作流程Fig .1 Particle segmentation workflow with modified watershed algorithm

1.2 特征標(biāo)記

特征分割后，需要對圖中的各個特征(每個連通區(qū)域為一單獨的特征)進(jìn)行編號，然后再分別計算各個特征的屬性。Haralick給出了特征標(biāo)記算法[11]。圖2為特征標(biāo)記示意圖，圖中分別用1、2、3標(biāo)記了三個特征區(qū)域;此外,0為背景區(qū)域，作為特征之間的分隔線。

1.3 特征屬性計算

二維圖像特征通常有以下屬性，如面積、邊界周長、形狀因子、重心位置等。在分析顆粒相時，常用等效圓半徑R作為顆粒大小的一個度量。其定義為

(1)

式(1)中：A為顆粒面積。

圖2 特征標(biāo)記示意圖Fig .2 Characteristic mark schematic diagram

1.4 特征參數(shù)分組統(tǒng)計

對圖像中所有顆粒特征的等效圓半徑進(jìn)行計算，得到等效圓半徑列表。然后對顆粒的等效圓半徑按照大小分成若干組，利用直方圖工具開展統(tǒng)計，獲得顆粒相尺寸的頻率分布和累積分布圖。

2 實例驗證與分析

圖3(左)為南海東部海域某老油田A油田目標(biāo)區(qū)塊巖屑樣品的背散射MAPS圖像，圖像中顏色較深的部分為孔隙，灰色的部分為礦物顆粒。背散射圖像的灰度與礦物的平均原子序數(shù)相關(guān)，平均原子序數(shù)較小的其灰度值也較小，可以據(jù)此來觀察/粗估礦物/孔隙在巖心表面的分布情況。圖3(右上)是選取的一小塊區(qū)域進(jìn)行觀察，并進(jìn)行方法及參數(shù)測試。測試可行后應(yīng)用到整幅圖像中，獲得整張MAPS圖像的統(tǒng)計分析結(jié)果。對這一小塊區(qū)域進(jìn)行特征分割，獲得了顆粒相(圖3右下白色區(qū)域)和孔隙相(黑色區(qū)域)。需要注意的是，由于部分顆粒之間有接觸，在劃分顆粒相時采用常規(guī)的閾值分割無法分離這些有接觸的礦物碎屑顆粒，這時如果不進(jìn)行處理，勢必會高估礦物顆粒的大小，并低估顆粒的數(shù)量，從而影響計算結(jié)果。因此，在進(jìn)行顆粒相的統(tǒng)計時，還需要進(jìn)行分離接觸的顆粒相這一操作。

圖3 南海東部A油田巖心樣品的背散射MAPS圖像及其孔隙、顆粒分割圖(MAPS圖像大小16 mm×15 mm)Fig .3 Back-scattered MAPS and pore-grain image segmentation of A oilfield in eastern South China Sea (MAPS image size 16 mm×15 mm)

圖4為示例背散射圖像中礦物顆粒的分離效果圖，左圖為原始背散射圖像，中圖為顆粒和孔隙相分割圖(其中白色為顆粒相、黑色為孔隙相)，右圖為采用基于3D距離圖的分水嶺方法得到的顆粒分離圖，互相接觸的顆粒之間已經(jīng)用分水嶺線分隔開。

圖4 南海東部A油田MAPS圖像中礦物顆粒分離Fig .4 Grain separation of MAPS image of A oilfield in eastern South China Sea

基于分離的顆粒圖像，進(jìn)行特征分割、特征標(biāo)記、特征屬性計算、特征參數(shù)分組統(tǒng)計，最后應(yīng)用到整張MAPS圖像中可以獲取巖屑樣品的顆粒粒徑頻率分布圖(圖5)。同時，利用傳統(tǒng)室內(nèi)實驗粒徑分析結(jié)果對其進(jìn)行了驗證，兩者較為吻合(圖5)。

圖5 南海東部A油田顆粒粒徑的頻率分布圖Fig .5 Frequency distribution of particle size of A oilfield in eastern South China Sea

基于巖屑樣品的顆粒粒徑分布，通過地質(zhì)過程模擬法模擬顆粒的沉積過程、壓實作用和成巖作用，通過已有的地層孔隙度數(shù)據(jù)建立儲層代表性的三維數(shù)字巖心(圖6a)，并提取相應(yīng)的孔隙網(wǎng)絡(luò)模型[13](圖6b)。

圖6 基于地質(zhì)過程模擬法構(gòu)建的南海東部A油田三維 數(shù)字巖心(模型大?。?0 mm×10 mm×10 mm)Fig .6 Reconstruction of 3D digital rock of A oilfield in eastern South China Sea with geological process-based method(model size:10 mm×10 mm×10 mm)

以數(shù)字巖心和孔隙網(wǎng)絡(luò)模型為研究平臺，可以計算出儲層的孔隙度為24.36%，滲透率為8 105 mD，并分析儲層微觀孔隙結(jié)構(gòu)特征。儲層孔喉半徑分布11～200 μm(圖7a)，主平均值為40 μm；形狀因子分布在0.022～0.079(圖7b)，平均值為0.055；配位數(shù)分布在1～12(圖7c)，平均值為4.48。其中，形狀因子是為了定量表征孔隙空間中孔隙、喉道單元體的形狀特征，其定義如下：

G=Ap/P2

(2)

式(2)中：Ap為孔隙、喉道的橫截面面積，m2；P為孔隙、喉道單元體橫截面的周長，m。

圖7 南海東部A油田孔隙結(jié)構(gòu)特征分布Fig .7 Pore structure characteristic distribution of A oilfield in eastern South China Sea

3 結(jié)論

通過改進(jìn)分水嶺算法對顆粒圖像進(jìn)行計算，包括特征分割、特征標(biāo)記、特征屬性計算、特征參數(shù)分組統(tǒng)計等，最后應(yīng)用到整張MAPS圖像中可以獲取巖屑樣品的顆粒粒徑頻率分布圖。該方法有效解決了巖屑粒徑測量困難的現(xiàn)狀，并通過地質(zhì)過程模擬法建立儲層代表性的三維數(shù)字巖心，提取相應(yīng)的孔隙網(wǎng)絡(luò)模型，計算孔滲物性參數(shù)和孔隙-喉道結(jié)構(gòu)特征，可為缺少巖心數(shù)據(jù)的老油田剩余油挖潛提供基礎(chǔ)數(shù)據(jù)。