馬霄一 李呈呈 白 俊 馬中高

(中國石化石油物探技術(shù)研究院，江蘇南京 211103)

0 引言

大部分沉積巖表現(xiàn)出一定程度的各向異性特征，而具有層狀結(jié)構(gòu)的頁巖各向異性特征更明顯。頁巖的強各向異性特征導(dǎo)致頁巖脆性和有機質(zhì)含量預(yù)測精度難以保證。因此，建立有機質(zhì)含量、脆性等儲層物性參數(shù)與彈性參數(shù)的巖石物理關(guān)系，是應(yīng)用頁巖甜點地震預(yù)測技術(shù)的關(guān)鍵環(huán)節(jié)。

美國和加拿大等國的頁巖氣已進入商業(yè)化開發(fā)階段，有關(guān)頁巖的研究，主要集中在頁巖內(nèi)部成分、應(yīng)力等對頁巖各向異性的影響。中國的頁巖氣開發(fā)還處于初級階段，對于頁巖各向異性的研究相對較少。雖然頁巖一般表現(xiàn)為強各向異性，但業(yè)界通常還是采用Thomsen參數(shù)表征頁巖的各向異性特征[1]。相關(guān)研究表明，沉積巖的縱橫波各向異性參數(shù)間存在一定程度的線性相關(guān)性，通過線性回歸得到的ε和γ間的經(jīng)驗公式具有實際意義[1]?？傆袡C碳含量(TOC)是表征頁巖儲層特征的重要指標，它能夠有效地評價頁巖有機質(zhì)豐度。隨著頁巖TOC的增加，不僅能增大頁巖的儲集空間，而且在一定程度上也有利于頁巖游離氣的吸附。因此，TOC預(yù)測和評價至關(guān)重要。

Wang[1]將不同沉積巖(包括頁巖、砂巖和碳酸鹽巖)的縱橫波各向異性參數(shù)ε和γ無差別地交會于同一坐標系下，獲得了最佳擬合曲線表達式，認為沉積巖的縱橫波各向異性參數(shù)可以通過該表達式相互預(yù)測，且這種預(yù)測獨立于圍壓、孔隙流體以及巖性等因素。Vernik等[2-3]研究頁巖的聲波速度各向異性，分析和探討在干燥、鹽水飽和以及不同圍壓條件下ε和γ間的相關(guān)性；Sierra等[4]通過實驗研究頁巖的彈性波速度各向異性，發(fā)現(xiàn)頁巖定向排列的黏土礦物成分與各向異性具有較好的相關(guān)性；Prasad等[5]認為有機質(zhì)含量減少會增強頁巖骨架堅硬程度；金吉能等[6]、尹俊等[7]利用地震多屬性反演技術(shù)預(yù)測TOC的平面分布，預(yù)測值與實測值吻合程度較高，但多解性較強。

伍宇明等[8]以龍馬溪頁巖為例，建立了頁巖彈性各向異性參數(shù)與礦物分布之間關(guān)系，通過引入變異系數(shù)描述礦物或孔裂隙在不同方向的差異性，并通過兩個正交方向的變異系數(shù)建立微觀異質(zhì)性指標，用于表征微觀尺度上礦物或孔裂隙在不同方向的差異程度。宏觀波速與礦物的微觀變異系數(shù)明顯相關(guān)，暗示巖石彈性各向異性與礦物分布直接相關(guān)。王鑫朋等[9]研究了川東龍馬溪組頁巖滲透率各向異性和應(yīng)力敏感性，發(fā)現(xiàn)巖心滲透率隨著與層理面夾角增大明顯降低，不同巖樣的滲透率隨著有效應(yīng)力增大呈指數(shù)遞減，平行層理的頁巖滲透率應(yīng)力敏感性相對更強。袁和義等[10]基于橫觀各向同性介質(zhì)模型，通過測量巖樣不同方向的縱波波速和快、慢橫波速度，計算了長寧區(qū)塊龍馬溪組巖樣的剛度系數(shù)和彈性系數(shù)，發(fā)現(xiàn)頁巖的縱波波速和快、慢橫波速度均表現(xiàn)出各向異性，并且隨圍壓升高，波速升高速率變緩。

中國東北松遼盆地南部嫩江組是重要的烴源巖層系，近年鉆遇嫩江組的探井部分具有較好的油氣顯示。因此，研究該層系的頁巖巖石物理特征，對認識研究區(qū)頁巖油分布規(guī)律和指導(dǎo)未來頁巖油的勘探和開發(fā)具有重要的意義。在前人研究的基礎(chǔ)上，本文選取通化盆地幾組嫩江組巖心，測試不同圍壓條件下的縱、橫波速度，分析ε和γ間的相關(guān)性，探討TOC與頁巖縱橫波速度比的相關(guān)性。

1 實驗方案

實驗巖樣(圖1)取自中國東北通化盆地，巖性主要為黑色泥頁巖、灰黑色泥巖。巖心層理發(fā)育，天然裂縫較少，以鉆井誘導(dǎo)裂縫為主。根據(jù)已知資料(表1)，巖樣6/90和8/91有機碳含量分別為3.18%、3.62%，黏土礦物含量分別為35%、37%，脆性礦物含量均在50%以上，為有利儲層。

為了分析頁巖彈性參數(shù)在不同方向的變化規(guī)律，實驗巖樣從兩個方向鉆?。孩俅怪庇趯永矸较?，層理與軸向應(yīng)力的加載方向呈90°；②平行于層理方向，層理與軸向應(yīng)力的加載方向呈0°。實驗巖樣的直徑為25mm，高度為50mm。為了使實驗數(shù)據(jù)更可靠，分別取得8組、共16塊不同井位的巖心巖樣，每塊巖樣編號前面數(shù)字表示不同井段，后面的數(shù)字“0”或“90”，表示垂直于層理或平行于層理(圖1a和表1)。

表1 巖樣巖石物理參數(shù)

圖1 巖心樣本(a)及其層理(b)

本文實驗使用美特斯工業(yè)系統(tǒng)公司生產(chǎn)的MTS-815三軸應(yīng)力測試儀器，巖樣受力方式如圖2所示。實驗全過程巖樣處于干燥狀態(tài)，軸向應(yīng)力通過位移方式加載，圍壓通過液壓控制方式加載?？v、橫波換能器的頻率為500kHz，發(fā)射和接收探頭與試件直接接觸，采用垂直透射法，測試精度控制在1%左右。

圖2 實驗加載示意圖

2 頁巖各向異性參數(shù)表征

頁巖是一種強各向異性介質(zhì)，通?？梢悦枋鰹榫哂行D(zhuǎn)對稱軸的VTI介質(zhì)[11]。Thomsen[12]給出了表征VTI介質(zhì)的彈性參數(shù)，可以直觀表示巖石縱、橫波各向異性，即

(1)

(2)

式中：vP為縱波速度；vS為橫波速度；0°表示波沿對稱軸方向傳播，90°表示波垂直于對稱軸傳播。Thomsen理論雖然是建立在弱各向異性的假設(shè)基礎(chǔ)上，而頁巖通常表現(xiàn)為強各向異性特征，但Thomsen參數(shù)對頁巖各向異性的表征仍然具有較強的適用性。

3 不同壓力縱橫波速度測試

在干燥環(huán)境下，分別在0、5、15、25、35、45、55MPa壓力下對16塊巖樣進行超聲縱、橫波速度測試(測試過程中保持圍壓與軸壓一致)，結(jié)果分別如圖3和圖4所示。由圖可以看出，縱、橫波速度隨著壓力增加均增大，且在低圍壓下速度增大更為突出；縱、橫波速度在平行于層理方向大于垂直于層理方向。起始壓力(0MPa)條件下，在垂直于層理方向縱波速度為2694.84～4791.95m/s，橫波速度為1750.51～2817.18m/s；在平行于層理方向，縱波速度為4002.59～5486.42m/s，橫波速度為2427.51～2990.74m/s。在其他壓力條件下，速度測試也有類似結(jié)果。波的傳播速度與傳播方向有直接關(guān)系，垂直于層理時，遇到的層理面最多，測試結(jié)果反映所有層理對傳播的影響，速度較小；平行于層理傳播時，波傳播只會從其中某層通過(最快路徑)，速度較大。根據(jù)圖3、圖4的測試結(jié)果，計算8組巖樣不同壓力環(huán)境下兩個方向縱、橫波速度差(圖5)，發(fā)現(xiàn)有利儲層巖樣6/90、8/91的速度差異最大，這為判別有利儲層提供了新的依據(jù)。

圖3 16塊巖樣縱波速度隨壓力的變化曲線

圖4 16塊巖樣橫波速度隨壓力的變化曲線

圖5 8組巖樣不同壓力環(huán)境下兩個方向縱(a)、橫(b)波速度差統(tǒng)計直方圖

4 實驗結(jié)果分析

4.1 各向異性參數(shù)關(guān)系

根據(jù)圖3、圖4的實驗測量結(jié)果，應(yīng)用式(1)、式(2)計算不同壓力環(huán)境下8組巖樣的ε和γ，其交會分析結(jié)果如圖6所示。可以發(fā)現(xiàn)，所有的頁巖巖樣一般ε略大于γ；同一巖樣，隨壓力增大，ε和γ參數(shù)都有減小的趨勢；有利儲層的ε和γ要遠遠大于非有利儲層，這是否是該地區(qū)頁巖的必然規(guī)律，還需要大量實驗驗證。根據(jù)周楓等[13]的研究結(jié)果，頁巖各向異性參數(shù)與黏土礦物含量、有機質(zhì)含量呈正相關(guān)。本文選取的巖樣中有利儲層具有明顯的高TOC特征，所以可以認為TOC是影響巖石各向異性強度的重要因素之一。

圖6 不同壓力環(huán)境下8組頁巖巖樣的ε -γ交會圖

4.2 壓力對ε與γ間關(guān)系的影響

Vernik等[3]對Bazhenov、Bakken、Monterey、Niobrara和North Sea五個地區(qū)頁巖巖樣建立了實驗室測量數(shù)據(jù)庫，選取干燥條件下測試的所有縱、橫波參數(shù)ε與γ進行交會分析。圖7是圍壓為70MPa時頁巖的ε與γ交會分析結(jié)果，可見二者存在較好的線性相關(guān)性，且通過線性回歸可以得到該條件下的線性回歸方程。該線性擬合為多地區(qū)擬合，故相關(guān)系數(shù)(R)不高。Johnston等[14]對Devonian-Mississippian黑色頁巖，測量了三個不同圍壓(10、50、100MPa)環(huán)境下的縱、橫波速度，計算各向異性參數(shù)，其ε與γ交會分析結(jié)果如圖8所示。

圖7 五個地區(qū)頁巖巖樣ε與γ的交會分析結(jié)果[3]

應(yīng)用本文實驗7個壓力環(huán)境下得到的全部縱橫波各向異性參數(shù)數(shù)據(jù)進行ε與γ線性回歸分析，結(jié)果如圖9所示。與圖8對比可知，本文實驗的回歸方程與文獻[14]的結(jié)果較一致。從圖9可以看出，在單個壓力環(huán)境下，ε和γ在一定程度上具有線性相關(guān)性。縱波各向異性參數(shù)始終大于橫波各向異性參數(shù)，且在壓力較小時，這一現(xiàn)象更明顯。隨著壓力增大，縱、橫波各向異性參數(shù)間的差異逐漸減小。綜合分析不同壓力下線性回歸表達式，可以發(fā)現(xiàn)，隨著壓力的增大，斜率增大，從零壓力下的0.8638上升到55MPa壓力下的1.135；截距大致相同，約-0.041～-0.022，反映了這些巖樣來自于同一地區(qū)，具有相似的本征特征。隨壓力增大，縱橫波各向異性值均減小，但ε減小的幅度大于γ，這一規(guī)律已被Sondergeld等[15]的實驗數(shù)據(jù)分析證實。

圖8 Devonian-Mississippian黑色頁巖不同壓力環(huán)境下的ε與γ交會分析結(jié)果[14]

圖9 本文頁巖巖樣不同壓力環(huán)境下ε與γ的線性回歸分析

綜上所述，頁巖的ε與γ具有線性相關(guān)性。實際測量中，可利用這些關(guān)系式，在已知一個參量的情況下預(yù)測另一個參量，或者為另一個參量的測量提供參考。

4.3 壓力對縱橫波速度比的影響

目前，國內(nèi)各大油氣田在地震屬性預(yù)測儲層方面都有較廣泛的應(yīng)用[16-17]。建立縱橫波速度比與非儲層、氣層、水層、氣水同層之間的關(guān)系，采用交會分析技術(shù)，建立儲層物性與縱橫波速度比的關(guān)系，提取有利儲層的平面分布，刻畫優(yōu)質(zhì)儲層的空間展布特征，而縱橫波速度比經(jīng)常作為儲層敏感參數(shù)的重要指示。根據(jù)上述測試及分析結(jié)果，選取1組有利儲層(6/90)和3組非有利儲層(6/45、5/62、5/68)進行不同壓力下的縱橫波速度比分析。

圖10為本文選取的頁巖巖樣在不同壓力環(huán)境下的縱橫波速度比變化，可以看出，所有縱橫波速度比值隨著壓力增加而變大。由圖3和圖4可知，壓力增大，縱、橫波速度均增大，相同條件下，縱波速度的增加幅度明顯大于橫波速度，因此縱橫波速度比隨著壓力增加而增大，范圍為1.50～1.85；平行于層理方向的縱橫波速度比均大于垂直于層理方向，體現(xiàn)了縱橫波速度比也有明顯的各向異性。

圖10 不同壓力環(huán)境下4組頁巖巖樣兩個方向的縱橫波速度比變化

仔細觀察圖10可以發(fā)現(xiàn)，右側(cè)紅色標記的有利儲層6/90巖樣，縱橫波速度比隨壓力變化最大。其余巖樣壓力變化導(dǎo)致的縱橫波速度比變化基本都在2.7%以下，而6/90巖樣縱橫波速度比變化在平行于層理方向為3.3%，在垂直于層理方向可達8.7%，且有利儲層6/90巖樣具有較低的縱橫波速度比。

有利儲層具備較低縱橫波速度比，且對壓力敏感(隨壓力變化較大)，這可能與其TOC較高有關(guān)?？v橫波速度比表現(xiàn)出的各向異性特征對判斷TOC高低、確定是否為有利儲層有很大幫助。

4.4 TOC與vP/vS的關(guān)系

Vernik等[11]指出，TOC和密度間存在較好的相關(guān)性。相對黏土和方解石(密度高達2.7g/cm3)，有機質(zhì)具有相對低的密度。此外，幾乎所有頁巖巖樣孔隙度均在1%～5%范圍內(nèi)，TOC會在很大程度上影響頁巖的密度。一般來說，密度越小，TOC越高，因此可以尋找一個動態(tài)參量，通過密度建立與頁巖TOC的對應(yīng)關(guān)系[18]。測井數(shù)據(jù)已經(jīng)證實可以利用縱橫波速度比區(qū)分并識別巖性。由本文實驗分析可知，vP/vS對頁巖TOC具有一定指示作用。

根據(jù)實際儲層環(huán)境選擇45MPa作為儲層壓力，將該壓力下實驗獲得的vP/vS數(shù)據(jù)與零壓下的密度做交會分析，結(jié)果如圖11所示?？梢妚P/vS隨密度的增大而增大，因此vP/vS與TOC呈負相關(guān)。有利儲層6/90、8/91巖樣屬高有機質(zhì)含量頁巖(表1)，其兩個方向的vP/vS都小于1.7，這與Vernik等[11]和Yan等[19]把vP/vS<1.6作為高TOC指標大致相符。

圖11 縱橫波速度比與密度的交會分析結(jié)果

4.5 彈性模量的影響因素

根據(jù)實驗測量的兩個方向的縱、橫波速度和孔隙度等數(shù)據(jù)，計算頁巖巖樣的楊氏模量和泊松比，如圖12所示。其中平行于層理方向巖樣位于右上方(高彈性模量)紅色橢圓內(nèi)；垂直于層理方向巖樣位于左下方(低彈性模量)藍色橢圓內(nèi)。平行于層理方向巖樣孔隙度隨彈性模量(楊氏模量、泊松比)增大而增高(紅色箭頭所示)，垂直于層理方向巖樣孔隙度隨彈性模量(楊氏模量、泊松比)增大而降低(藍色箭頭所示)。

圖12 不同方向巖樣的彈性模量與孔隙度交會圖

總的來說，頁巖彈性模量與測量方向、孔隙度都有關(guān)系，且表現(xiàn)出一定的規(guī)律，可為后續(xù)的頁巖脆性和甜點預(yù)測提供地震定量解釋依據(jù)[20-21]。

5 結(jié)論

根據(jù)三軸壓縮實驗和超聲波測試，獲得了頁巖巖樣垂直于層理方向和平行于層理方向的縱、橫波速度，分析結(jié)果表明：

(1)有利儲層各向異性參數(shù)值更大，證明了有機質(zhì)含量顯著影響頁巖的各向異性程度；

(2)應(yīng)用線性回歸方法獲得了ε和γ的經(jīng)驗表達式，從而為頁巖縱、橫波各向異性參數(shù)的相互預(yù)測提供一種簡捷而實用的方法；

(3)有機質(zhì)含量與vP/vS存在負相關(guān)性，有利儲層vP/vS<1.7，可為有利儲層預(yù)測提供依據(jù)；

(4)不同方向的彈性模量隨孔隙度變化呈不同的變化趨勢，平行于層理方向的彈性模量隨孔隙度增高變大，垂直于層理方向的彈性模量隨孔隙度增高變小，可為儲層脆性評價和壓裂等提供指導(dǎo)。