成鎖，簡世凱，3，譚楊，竇齊豐，劉云祥

（1.中國石油塔里木油田分公司勘探開發(fā)研究院，新疆 庫爾勒 841000；2.中國石油超深層復雜油氣藏勘探開發(fā)技術研發(fā)中心，新疆 庫爾勒 841000；3.深層油氣全國重點實驗室（中國石油大學（華東）），山東 青島 266580；4.長江大學，湖北 武漢 430000）

0 引言

塔里木盆地近40%的石油和天然氣儲量發(fā)現(xiàn)于碳酸鹽巖儲層， 其中塔北-塔中地區(qū)是重要的勘探地區(qū)。奧陶系良里塔格組、一間房組和鷹山組為重要的勘探層系。在多期次的構(gòu)造和巖溶作用的疊加下，形成了具有基質(zhì)孔隙、 溶蝕孔隙及裂縫等不同孔隙類型的碳酸鹽巖儲層［1-4］。 這些地質(zhì)特征，尤其孔隙結(jié)構(gòu)的變化引起復雜的孔隙度與地震波速度的關系， 給地震表征儲層孔隙類型造成很大挑戰(zhàn)。 表征儲層孔隙類型變化的巖石物理模型較多， 如最早的Wood 公式， 其對淺層、弱成巖高孔隙度碎屑巖儲層的孔隙度-地震波速度關系刻畫效果較好，可反映單一粒間孔的孔隙類型［5］。目前最常用的是威利公式刻畫速度-孔隙度關系，進而刻畫孔隙類型，且該模型對淺—中埋深地層、巖性較純的砂巖儲層晶間孔的刻畫效果比較好［5］。 為了刻畫復雜孔隙結(jié)構(gòu)儲層，Robert 等［6］對威利公式進行了修改，修改的威利公式對粒內(nèi)孔、 晶間孔及裂縫等3 種孔隙類型進行刻畫， 在碎屑巖儲層孔隙類型刻畫中取得了一定的效果［7］。

上述3 種方法都是針對比較單一、 巖性較純的碎屑巖孔隙類型的刻畫，沒有針對碳酸鹽巖儲層。通過對Biot 理論模型進行修改和延展，Saleh 等［7-8］提出了一個新的巖石物理模型， 在該模型中孔構(gòu)參數(shù)被引入來表征孔隙類型的變化，該參數(shù)是縱波速度、橫波速度及密度等地震彈性參數(shù)的函數(shù)， 該模型為碳酸鹽巖儲層孔隙類型表征的巖石物理模型， 并在巖心尺度得到應用和驗證［5］。 在Sun 等研究的基礎上，Dou 等［9-11］在測井尺度實現(xiàn)了孔構(gòu)參數(shù)的反演， 并在美國二疊系盆地碳酸鹽巖儲層孔隙結(jié)構(gòu)變化刻畫中得到很好的應用，將該模型從數(shù)學模型建立到實際應用推進了一步。

塔里木盆地縫洞型碳酸鹽巖儲層礦物成分、 孔隙類型復雜多樣，導致其體積-剪切模量與孔隙度交會圖十分發(fā)散，即其體積-剪切模量在孔隙度恒定時，體積-剪切模量存在較大差異， 傳統(tǒng)的巖石物理模型不能準確刻畫二者之間的關系［12-21］。 因此，基于研究區(qū)大量的常規(guī)測井、偶極聲波測井、成像測井及巖心分析數(shù)據(jù)，本次研究利用Sun 巖石物理模型對研究區(qū)碳酸鹽巖儲層的孔隙類型、溶洞直徑、溶洞充填程度、溶洞面孔率及裂縫進行表征， 為研究區(qū)復雜縫洞體優(yōu)質(zhì)儲層的預測奠定基礎。

1 方法原理

地震波在地層中傳播的速度是油藏地球物理中重點研究的一個參數(shù)， 它是連接油藏參數(shù)和巖石聲學特征的橋梁。 Sun 巖石物理模型是目前進行儲層孔隙類型變化表征的定量模型之一， 相對于常規(guī)的巖石物理模型能較好地定量表征碳酸鹽巖儲層孔隙類型變化，該模型中孔構(gòu)參數(shù)相對獨立于孔隙度參數(shù)， 同時表征在孔隙度不變的條件下縱波速度、橫波速度的變化，進而表征孔喉半徑及孔隙的形態(tài)變化。 該模型可簡化為

式中：vp，vs分別為縱、橫波波速；ρ 為體積密度；K，μ 分別為體積模量和剪切模量；ρs為固體基質(zhì)的密度；Ks為體積模量；ρf為流體密度；Kf為流體體積模量；? 為孔隙度；?k為孔隙彈性參數(shù)；f，γ，fu，γu分別為骨架剛度和柔度的不同參數(shù)。

式（6）中γ 不僅獨立于孔隙度大小，而且跟孔隙的形態(tài)有關。 理論上，γ 是孔隙長短軸的比值，即圓的孔隙γ 低，扁的孔隙γ 高，線狀的裂縫γ 值最高。 因此，這個孔隙彈性理論模型在巖心和測井尺度可以很好表征儲層內(nèi)部不同孔隙結(jié)構(gòu)類型的變化及孔隙間的連通性。 因為它是縱波速度、橫波速度及密度的函數(shù)，可以將儲層孔隙類型的變化和地球物理彈性參數(shù)聯(lián)系起來，達到地震反演的目的。測井尺度的孔構(gòu)參數(shù)γ 的計算式為

2 實例分析

2.1 成像測井二次處理及解釋

對于溶蝕及裂縫油氣藏， 成像測井是計算與評價其定量評價參數(shù)的有效手段之一。 通過對成像測井數(shù)據(jù)的處理，首先得到了孔洞儲層的孔洞數(shù)量、直徑、面孔率，裂縫儲層長度、寬度、孔隙度等定量參數(shù)，這些是利用Sun 巖石物理模型進行孔隙類型表征的重要數(shù)據(jù)。 本次研究利用電阻率成像測井， 測井儀器記錄的150 條微電阻率曲線為陣列數(shù)據(jù)。 首先，這些陣列數(shù)據(jù)經(jīng)過處理（具體包括深度校正、加速度校正、均衡處理等）， 就把全部微電阻率值進行了灰度或彩色刻度，將每個電極的每個采樣點電阻率變成一個色元；然后，采用色板（通常的色板為黑—棕—黃—白，暗色代表低電阻率，淺色或亮色代表高電阻率）進行配色，形成彩色圖像，圖像顏色代表所測地層電阻率的高低。利用圖像色彩變化將電阻率的高低轉(zhuǎn)化為地層特征變化的地質(zhì)信息，即根據(jù)地層電阻率的變化，來進行裂縫、孔洞發(fā)育程度及孔洞充填等的定量和定性解釋。 圖1 為研究區(qū)內(nèi)A 井成像測井二次處理解釋結(jié)果。

圖1 研究區(qū)A 井裂縫段成像測井二次處理解釋成果Fig.1 Interpretation result of imaging logging secondary processing of fracture sections of Well A in study area

由圖1 可以看出：在6 177 m 時裂縫最為發(fā)育。 裂縫長度為17 m，裂縫密度為2 條/m，裂縫寬度為0.017 mm，裂縫孔隙度為0.025，在該深度都達到最大。

2.2 Sun 巖石物理模型儲層評價

由于孔隙類型及巖相等地質(zhì)因素影響， 地震波速度與地層孔隙度表現(xiàn)出十分復雜的關系。 由研究區(qū)A井的孔隙度和縱波速度的相關研究可知， 隨孔隙度的增大，縱波速度在降低。這與威利公式里刻畫的孔隙度與縱波速度的關系一致。但在孔隙度不變的情況下，縱波速度表現(xiàn)出很大的變化。 如孔隙度為8%左右的樣品，縱波速度變化范圍為3 900～5 900 m/s。這種恒定孔隙度下速度的差異是受不同地區(qū)地質(zhì)因素的影響造成的。 地質(zhì)因素主要為礦物組分、巖性及孔隙類型等。 本次研究綜合巖心分析數(shù)據(jù)及成像測井二次定量解釋結(jié)果， 利用Sun 巖石物理模型反演孔構(gòu)參數(shù)分析孔隙類型、溶洞特征及裂縫的變化。

成像結(jié)果分析表明， 孔構(gòu)參數(shù)大于7 的層段基本是泥質(zhì)充填的孔洞（見圖2）。 圖2a 是縱波速度和孔隙度的交會圖， 粉色框內(nèi)的點是孔構(gòu)參數(shù)在7 以上的層段。 分析這些數(shù)據(jù)點在FMI 上的特征可以得出，該段為溶蝕孔洞段，孔洞內(nèi)有片狀的泥質(zhì)條帶分布，溶蝕孔隙雖然規(guī)模很大，但都被泥質(zhì)充填，為無效孔隙。 這說明溶蝕孔隙儲層內(nèi)，模型中高孔構(gòu)參數(shù)反映的是溶蝕、充填的孔隙類型。

圖2 測井孔隙度-縱波速度交會圖及成像特征Fig.2 Well logging porosity-P-wave velocity crossplot and image characteristics

孔構(gòu)參數(shù)小于5 的樣品， 推測溶洞直徑在70 mm以上，溶洞面孔率在0.15 以上。

孔構(gòu)參數(shù)在5～6 的樣品，孔隙度在10%～15%。 根據(jù)成像測井解釋，確定其為半充填、弱充填溶洞，孔洞數(shù)量平均為15 個， 最大可達25 個， 最大溶洞直徑約160 mm，平均在80 mm，成像面孔率最大達0.6，平均在0.15。 這說明該段溶洞比較發(fā)育，具體成像見圖3。

圖3 孔隙度-縱波速度交會圖及孔構(gòu)參數(shù)5～6 段的成像處理特征Fig.3 Porosity-P-wave velocity crossplot and image processing characteristics of 5-6 structure parameter sections

孔構(gòu)參數(shù)在6～7 的樣品， 孔隙度在4.5%～10.0%。根據(jù)成像測井，確定其為半充填溶洞，且其規(guī)模發(fā)育中等偏小，孔洞數(shù)量平均為4 個，最大可達15 個，最大溶洞直徑約為70 mm，均值在15 mm，最大成像面孔率約為0.05，平均在0.025。相對于孔構(gòu)參數(shù)大于7 的層段，孔構(gòu)參數(shù)在6～7 的層段是溶蝕孔洞相對發(fā)育較弱的溶洞段（見圖4）。

圖4 孔構(gòu)參數(shù)6～7 段的成像處理特征Fig.4 Image processing characteristics of 6-7 structure parameter sections

對于具有高孔隙度的溶洞型儲層， 孔構(gòu)參數(shù)可以 刻畫溶洞發(fā)育的規(guī)模。 孔隙度2.0%～4.5%的基質(zhì)孔隙段，孔構(gòu)參數(shù)基本在7 以上，隨著其增大，孔隙類型從晶間孔變化到微晶孔， 儲層變得更致密； 孔隙度小于2%、孔構(gòu)參數(shù)大于12 的層段，表現(xiàn)為裂縫特征。 最終解釋圖版見圖5。

圖5 孔隙類型及溶洞大小孔構(gòu)參數(shù)解釋圖版Fig.5 Interpretation chart of pore type,vug size and pore structure parameter

3 結(jié)論

1） 針對塔里木奧陶系具有復雜的孔隙類型的儲層，本文提出的基于Sun 巖石物理模型表征方法，基本達到在溶蝕孔隙、 基質(zhì)孔隙及裂縫段的溶蝕孔徑尺度及充填程度、孔隙類型表征的目的。

2）對于孔隙度大于4.5%的溶蝕孔隙，在孔構(gòu)參數(shù)大于7 的樣品，基本為泥質(zhì)充填的孔隙。 對于弱充填、半充填或未充填的溶蝕孔， 孔構(gòu)參數(shù)能達到表征溶蝕孔隙半徑、面孔率的目的，即隨著孔構(gòu)參數(shù)的降低，面孔率及溶蝕孔隙半徑增大，溶蝕儲層的質(zhì)量變好。對于孔隙度2.0%～4.5%基質(zhì)孔隙， 孔構(gòu)參數(shù)也能表征孔隙類型的變化，即孔構(gòu)參數(shù)增大，孔隙類型具有從晶間孔到晶間微孔及微裂縫的發(fā)育特征。 裂縫基本發(fā)育于孔隙度小于2% 的非儲層段， 表現(xiàn)為低速及高孔構(gòu)參數(shù)的特征。

3）本文的研究成果對類似區(qū)儲層孔隙類型的刻畫有一定的指導意義。 雖然該模型一定程度實現(xiàn)了孔隙類型的刻畫，但沒有考慮各向異性，對裂縫的刻畫需要進行深入的研究。