周鳳鳴, 趙岳, 莊東志, 謝偉彪, 馬明明

(1.冀東油田勘探開發(fā)研究院, 河北 唐山 063004; 2.中國石油大學(xué)(華東), 山東 青島 266000)

0 引 言

橫波速度是研究地層性質(zhì)的重要參數(shù),由于其測量比較困難,因此進行地層橫波速度預(yù)測具有重要意義[1-3]。利用經(jīng)驗公式進行飽含水的碎屑硅酸鹽、石灰?guī)r、白云巖、砂巖等巖性的橫波預(yù)測較為常見,并應(yīng)用于AVO技術(shù)的理論研究及應(yīng)用[4-7],經(jīng)驗公式的方法主要是基于巖石孔隙度這個橋梁,沒有考慮巖石孔隙形態(tài)、尺度等孔隙結(jié)構(gòu)對關(guān)系模型建立的影響。利用理論模型進行橫波速度預(yù)測法也是常用方法之一,1995年Xu和White 聯(lián)合K-T模型和微分等效介質(zhì)理論和Gassmann理論,提出了砂、泥巖混合巖石的理論模型。為了減小Xu-White模型參數(shù)提取誤差對預(yù)測精度的影響,喬悅東等[8]通過引進虛擬孔隙度變量和非線性優(yōu)化算法。孫福利等[9]采用Biot-Gassmann基于基質(zhì)彈性模量與骨架彈性模量關(guān)系,通過迭代計算方法得到橫波速度。這些理論計算模型,考慮了孔隙形狀等參數(shù)對橫波速度的影響,但是需要給定的參數(shù)較多,計算時容易出現(xiàn)多解性。

大量研究證實,地層橫波速度不僅與地層巖石孔隙度有關(guān),也與孔隙結(jié)構(gòu)關(guān)系緊密,特別受孔隙形狀影響。核磁共振橫向弛豫時間T2幾何平均值是反映孔隙結(jié)構(gòu)特征的參數(shù),與巖石孔隙比表面、形狀、孔喉半徑等有關(guān)。因此,兩者之間有著必然的相關(guān)性。本文以理論分析和巖石物理實驗為基礎(chǔ),確定核磁共振T2幾何平均值與巖石橫波波速的關(guān)系。通過核磁共振和橫波速度測量實驗,建立核磁共振T2幾何平均值與橫波波速之間關(guān)系。實驗結(jié)果表明,核磁共振T2幾何平均值與橫波波速之間存在冪函數(shù)關(guān)系,實例驗證表明該方法模型簡單實用,效果較好。

1 核磁共振T2譜與橫波速度預(yù)測基本原理

巖石孔隙度倒數(shù)和橫波速度之間存在線性關(guān)系,即橫波速度可表示為

(1)

式中,vs為橫波速度;φ為巖石孔隙度;a和b為系數(shù)。

根據(jù)核磁共振橫向弛豫時間T2與孔隙半徑r之間的關(guān)系[10]

(2)

ri=C1T2i

式中,Fs為孔隙形狀因子,球形孔隙取Fs=3,柱狀和板狀喉道取Fs=2;ri為孔隙半徑;T2i為第i種孔隙的橫向弛豫時間;ρ為表面弛豫強度;C1為系數(shù)。

實驗表明,平均孔喉半徑r與T2幾何平均值成線性正比關(guān)系[10]

r=C2T2gm

(3)

根據(jù)巖石物理體積模型,孔隙平均孔喉半徑近似表示為[11-12]

r=2φ/S

(4)

式中,r為孔隙平均孔喉半徑;S為巖石比面;系數(shù)2是一個理論值,實際應(yīng)用時是一個隨地區(qū)不同而變化的數(shù)值,應(yīng)按具體資料統(tǒng)計確定。

結(jié)合式(1)、式(3)、式(4)得到理想模型

(5)

將式(5)化簡可以得到

(6)

真實地層中關(guān)系式為

(7)

理論證明橫波波速與T2幾何平均值之間存在冪函數(shù)關(guān)系。式(6)是在理想的巖石物理體積模型下得出的,A為化簡后的系數(shù)。式(7)是真實地層中T2幾何平均值與橫波波速關(guān)系式,系數(shù)A以及T2gm的指數(shù)大小a與孔隙形狀因子、比表面、巖性等有關(guān),需要根據(jù)不同地區(qū)實驗數(shù)據(jù)擬合得出。在真實地層中,孔隙結(jié)構(gòu)通常較復(fù)雜,故T2的指數(shù)的絕對值一般小于1。

2 實驗結(jié)果及分析

對工區(qū)的巖石樣品使用主頻為2 MHz的IPC-801核磁共振信號巖心分析儀,回波間隔為0.37 ms,等待時間為6 000 ms,累加次數(shù)為128次,回波個數(shù)為4 096,在室溫25 ℃條件下進行核磁共振信號采集;采用SIRT反演方法,迭代次數(shù)5 000進行T2譜反演。通過得到的核磁共振T2譜、T2分布點和對應(yīng)孔隙度分量,進而計算得到T2幾何平均值。對工區(qū)的巖石樣品使用計算機脈沖信號發(fā)生接收器產(chǎn)生和接受脈沖信號,應(yīng)用Olympus橫波探頭測量橫波信號;使用混合信號示波器采集橫波波形,根據(jù)橫波波形判斷橫波到時,并由此計算橫波速度。

2.1 飽含水巖樣核磁共振T2與橫波速度關(guān)系

圖1 砂巖vs—T2gm關(guān)系曲線

通過對22塊砂巖巖樣進行核磁共振與橫波速度測量,22塊巖樣為洗油后100%水飽和巖樣。表1顯示了T2幾何平均值和橫波速度存在良好的相關(guān)關(guān)系,圖1顯示T2幾何平均值增大,橫波波速降低。分析表明,核磁共振T2幾何平均值與橫波波速之間存在冪函數(shù)關(guān)系

(8)

式中,vs為橫波波速;T2gm為幾何平均值。

同樣,對11塊含礫不等粒砂巖進行核磁共振與橫波速度測量,二者也具有良好相關(guān)系(見圖2),實驗結(jié)果見表2。建立核磁共振T2幾何平均值與橫波波速之間的冪函數(shù)關(guān)系

(9)

表1 飽含水砂巖核磁共振與橫波速度實驗數(shù)據(jù)表

*非法定計量單位,1 mD=9.87×10-4μm2,下同

圖2 含礫不等粒砂巖vs—T2gm關(guān)系曲線

表2飽含水含礫不等粒砂巖核磁共振與橫波速度實驗數(shù)據(jù)表

巖樣編號T2幾何平均值/ms橫波速度/(m·s-1)孔隙度/%滲透率/mD111 582831948 0138812 434 2706215 310181849 3678311 018 92936 936422039 667127 390 296945 203512345 050095 750 0566521 739411740 6627410 7214 0515618 704731701 5125512 3827 8336720 038111639 0696213 6544 279283 66272481 982774 5720 0034910 172971851 6638212 745 8492107 698041843 8743313 642 60921115 217511743 6712911 3115 7831

2.2 束縛水條件巖樣核磁共振T2與橫波速度關(guān)系

橫波速度經(jīng)常用來進行儲層流體識別,縱橫波速比、泊松比、孔隙壓縮系數(shù)等[13-14]是進行天然氣儲層識別常用手段,其中均離不開橫波速度這個參數(shù)。在典型天然氣儲層中,孔隙流體除天然氣外,主要是不可動的束縛水,同時考慮到氣體性質(zhì)對橫波速度的影響很小,而且天然氣的含氫指數(shù)一般較小,特別是在中深層儲層中,束縛水含量較高,導(dǎo)致氣體性質(zhì)變化對儲層的T2幾何平均值影響較小。因此,在實際利用儲層的核磁共振T2幾何均值預(yù)測橫波和巖樣核磁共振實驗的時候,一般可用空氣替代天然氣,其實驗結(jié)果有一定的可替代性;在一定條件下研究對核磁共振T2譜的含氣和含烴校正,也能提高利用核磁共振測井預(yù)測橫波的精度。

圖3 束縛水條件vs—T2gm關(guān)系曲線

采用GL-21M高速冷凍離心機,使用高速離心方法,對巖樣進行減水飽和度處理,直至巖樣含水達到束縛水飽和度條件。束縛水巖樣核磁共振T2譜和巖石橫波速度測量實驗數(shù)據(jù)見表3,顯示出二者具有良好相關(guān)性(見圖3)。核磁共振T2幾何平均值與橫波波速之間的關(guān)系為

(10)

圖4 某井砂巖段橫波波速預(yù)測成果圖*非法定計量單位,1 ft=12 in=0.304 8 m,下同

表3 砂巖(束縛水)核磁共振與橫波速度測量實驗數(shù)據(jù)表

3 應(yīng)用實例

基于目標區(qū)巖石物理實驗數(shù)據(jù),建立T2幾何平均值與橫波波速關(guān)系模型,進行橫波速度預(yù)測。如圖4所示,某泥質(zhì)砂巖井段,儲層孔隙度約3%～15%,滲透率小于1×10-3μm2,電阻率3～40 Ω·m,對4個層段砂巖進行橫波速度預(yù)測,并將其同偶極橫波測井實測橫波速度進行對比,顯示二者非常吻合。

不同區(qū)塊、不同巖性得到的T2幾何平均值與橫波波速關(guān)系模型不同,橫波預(yù)測需要根據(jù)不同的巖性特征和流體性質(zhì)分別建立,以提高橫波波速預(yù)測精度。

4 結(jié) 論

(1) 從理論上分析核磁共振T2幾何平均值與橫波波速之間存在冪函數(shù)關(guān)系,T2幾何平均值增大,橫波波速降低;基于巖石物理實驗數(shù)據(jù)驗證了該關(guān)系的正確性,進而證實了利用核磁共振T2幾何均值和橫波波速之間的關(guān)系建立區(qū)域模型預(yù)測橫波波速的可行性。

(2) 核磁共振T2幾何平均值不僅包含孔隙度的影響,而且還考慮了孔隙尺度、比表面等孔隙結(jié)構(gòu)的影響。利用T2幾何平均值預(yù)測橫波波速的方法巖石物理基礎(chǔ)更完善,物理意義更明確,而且模型建立方法簡單、實用,豐富了橫波預(yù)測的方法。

(3) 橫波波速受巖性影響較大,對工區(qū)建立模型時,應(yīng)劃分不同的巖性結(jié)合不同的流體性質(zhì)進行建立,以提高橫波波速預(yù)測精度。

參考文獻:

[1] 郭棟, 王興謀, 印興耀, 等. 橫波速度計算方法與應(yīng)用 [J]. 油氣地球物理, 2007, 5(3): 18-22.

[2] 孟憲軍, 姜秀娣, 黃捍東, 等. 疊前AVA廣義非線性縱、橫波速度反演 [J]. 石油地球物理勘探, 2004, 39(6): 645-650.

[3] 熊曉軍, 林凱, 賀振華. 基于等效彈性模量反演的橫波速度預(yù)測方法 [J]. 石油地球物理勘探, 2012, 47(5): 723-728.

[4] Castagna J P, Batzle M L, Eastwood R L. Relationships Between Compressional-wave and Shear-wave Velocities in Clastic Silicate Rocks [J]. Geophysics, 1985, 50(4): 571-581.

[5] Han D, Nur A, Morgan D. Effects of Porosity and Clay Content on Wave Velocities in Sandstones [J]. Geophysics, 1986, 51(11): 2093-2107.

[6] Castagna J P, Backus M M. Offset-dependent Reflectivity——Theory and Practice of AVO Analysis [J].

Society of Exploration Geophysicists, 1993(8): 135-181.

[7] 侯波, 陳小宏, 張孝珍, 等. 基于近似流體替代的快速無迭代橫波速度預(yù)測方法 [J]. 科技導(dǎo)報, 2011, 29(31): 34-37.

[8] 喬悅東, 高云峰, 安鴻偉. 基于Xu-White模型的優(yōu)化測井橫波速度預(yù)測技術(shù)研究與應(yīng)用 [J]. 石油天然氣學(xué)報, 2007, 29(5): 100-102.

[9] 孫福利, 楊長春, 麻三懷, 等. 橫波速度預(yù)測方法 [J]. 地球物理學(xué)進展, 2008, 23(2): 470-474.

[10] 李天降, 李子豐, 趙彥超, 等. 核磁共振與壓汞法的孔隙結(jié)構(gòu)一致性研究 [J]. 天然氣工業(yè), 2006, 26(10): 57-59.

[11] 原海涵. 用測井方法研究巖石孔隙的平均孔道半徑 [J]. 測井技術(shù), 1983, 7(2): 21-26.

[12] 岳文正, 趙克超. 用核磁共振及巖石物理實驗求地層束縛水飽和度及平均孔隙半徑 [J]. 測井技術(shù), 2002, 26(1): 22-25.

[13] 楊孛, 李瑞, 康坤. 縱、橫波速度和密度重構(gòu)算法在蘇里格地區(qū)的應(yīng)用 [J]. 西安石油大學(xué)學(xué)報, 2012, 27(4): 1-7.

[14] 李宏兵, 張佳佳, 姚逢昌. 巖石的等效孔隙縱橫比反演及其應(yīng)用 [J]. 地球物理學(xué)報, 2013, 56(2): 608-615.