田曉冬，劉得芳，殷秋麗，莊東志

中國(guó)石油冀東油田分公司 勘探開(kāi)發(fā)研究院，河北 唐山 063000

0 引言

與沉積巖儲(chǔ)層相比，火山巖儲(chǔ)層是一類 “特殊”的儲(chǔ)層類型。自發(fā)現(xiàn)以來(lái)，人們對(duì)火山巖儲(chǔ)層的形成、分布和測(cè)井評(píng)價(jià)十分關(guān)注[1]。幾十年來(lái)，對(duì)火山巖儲(chǔ)層控制因素的研究不斷深入，但總體都是圍繞酸性火山巖展開(kāi)，隨著國(guó)內(nèi)眾多中基性火山巖具有優(yōu)良的油氣藏儲(chǔ)集性的發(fā)現(xiàn)，中基性火山巖儲(chǔ)層測(cè)井評(píng)價(jià)受到越來(lái)越多的國(guó)內(nèi)研究人員的關(guān)注[2-3]。火山巖儲(chǔ)層巖石類型、孔隙結(jié)構(gòu)復(fù)雜多樣，學(xué)者們對(duì)火山巖儲(chǔ)層做了大量的研究工作，包括巖性識(shí)別、物性評(píng)價(jià)及流體識(shí)別等，但研究對(duì)象多為未蝕變的火山巖，而對(duì)于蝕變火山巖的測(cè)井解釋研究相對(duì)較少。王春燕等[4]利用巖性指數(shù)和蝕變指數(shù)建立火山巖蝕變判別方法，校正了蝕變后的電阻率曲線，根據(jù)校正后的電阻率與孔隙度建立了一種蝕變火山巖氣水層識(shí)別的方法；朱學(xué)娟等[5]為提高儲(chǔ)層儲(chǔ)集性能評(píng)價(jià)的實(shí)際性，對(duì)蝕變玄武巖儲(chǔ)層的有效性進(jìn)行了分析；申波等[6]利用主成分分析(PCA)技術(shù)，通過(guò)連續(xù)提取巖石的變骨架參數(shù)計(jì)算了蝕變地層的孔隙度；張麗華等[7]建立了包含基質(zhì)、裂縫和非連通孔洞的三重孔隙模型，構(gòu)建了包含黏土影響火山巖儲(chǔ)層的飽和度方程，但未給出黏土含量的計(jì)算方法；楊雪[8]研究了各類原巖能夠發(fā)生的蝕變種類，構(gòu)造了以單巖性骨架礦物、孔隙和黏土礦物三組分體積模型，結(jié)合變骨架中子-密度交會(huì)圖，研究了不同巖性火山熔巖的孔隙度和蝕變黏土含量；孫茹雪[9]建立了考慮蝕變的中基性火山巖儲(chǔ)層模型，利用五組分人工蜂群算法(ABC)對(duì)火山巖儲(chǔ)層進(jìn)行了評(píng)價(jià)計(jì)算；高衍武等[10]利用黏土蝕變程度指數(shù)對(duì)黏土化蝕變火山巖的黏土含量進(jìn)行定量計(jì)算，利用體積模型對(duì)黏土化蝕變的火山巖孔隙度進(jìn)行了計(jì)算。

前人關(guān)于蝕變火山巖孔隙度計(jì)算的研究較多[6-10]，而針對(duì)蝕變火山巖飽和度的研究則相對(duì)較少。南堡地區(qū)沙河街組中基性火山巖發(fā)育，由薄片分析資料發(fā)現(xiàn)，研究區(qū)A、B、C 3口井存在不同程度的蝕變現(xiàn)象?；鹕綆r蝕變會(huì)直接影響到儲(chǔ)層的孔隙結(jié)構(gòu)，產(chǎn)生的黏土礦物對(duì)測(cè)井響應(yīng)影響較大，如果忽視蝕變對(duì)儲(chǔ)層的影響，將嚴(yán)重影響到儲(chǔ)層測(cè)井解釋的準(zhǔn)確性，進(jìn)而導(dǎo)致儲(chǔ)層劃分不準(zhǔn)而造成一定的經(jīng)濟(jì)損失。因此，筆者通過(guò)研究蝕變作用對(duì)儲(chǔ)層物性、測(cè)井響應(yīng)特征等帶來(lái)的影響，構(gòu)建能夠定性判別蝕變程度的指標(biāo)因子，并根據(jù)自然伽馬、深側(cè)向電阻率、聲波、密度和中子5條測(cè)井曲線擬合預(yù)測(cè)黏土含量，定量評(píng)價(jià)儲(chǔ)層的蝕變程度。通過(guò)巖電實(shí)驗(yàn)發(fā)現(xiàn)火山巖蝕變前后阿爾奇公式中的飽和度指數(shù)發(fā)生變化，并對(duì)儲(chǔ)層飽和度的評(píng)價(jià)產(chǎn)生了影響。該研究對(duì)增強(qiáng)火山巖儲(chǔ)層的認(rèn)知水平具有實(shí)際意義，對(duì)火山巖儲(chǔ)層油氣藏精細(xì)評(píng)價(jià)及后期勘探開(kāi)發(fā)提供理論參考依據(jù)。

1 火山巖蝕變及產(chǎn)物

1.1 火山巖蝕變機(jī)理

巖石的蝕變作用，是指巖石在一定的環(huán)境下，經(jīng)過(guò)一系列的化學(xué)和物理變化導(dǎo)致微觀組成和結(jié)構(gòu)構(gòu)造發(fā)生改變的作用?；鹕綆r發(fā)生蝕變可能由2種作用引起：①風(fēng)化淋濾作用，②熱液交代蝕變作用。

風(fēng)化淋濾作用包括物理和化學(xué)風(fēng)化作用，物理風(fēng)化作用是指巖石只發(fā)生機(jī)械破碎而化學(xué)成分未改變的風(fēng)化作用。物理風(fēng)化的結(jié)果是使巖石破碎，形成各種碎屑物質(zhì)?；瘜W(xué)風(fēng)化作用是巖石與礦物在一定的環(huán)境條件下(水、氧氣和二氧化碳)，經(jīng)過(guò)一系列的化學(xué)反應(yīng)，發(fā)生分解或溶蝕，改變巖石的化學(xué)成分，形成次生孔隙和溶蝕孔洞等[11]。

熱液交代蝕變作用是熱液在開(kāi)放或半封閉系統(tǒng)內(nèi)運(yùn)動(dòng)過(guò)程中產(chǎn)生的，熱液物質(zhì)會(huì)以離子、分子(團(tuán))等形式運(yùn)動(dòng)，沿各種大小裂隙，包括粒間孔隙、毛細(xì)裂隙和礦物的節(jié)理等對(duì)原巖進(jìn)行交代或物質(zhì)置換。熱液活動(dòng)的存在，使巖石內(nèi)部結(jié)構(gòu)、構(gòu)造和礦物共生關(guān)系變得非常復(fù)雜。在熱液交代過(guò)程中，溶解和沉淀幾乎是同時(shí)進(jìn)行的，通常遵循等體積定律，即巖石體積一般不會(huì)發(fā)生改變[12]。在南堡研究區(qū)，中基性火山巖多發(fā)生熱液交代蝕變作用。

1.2 火山巖蝕變類型及產(chǎn)物

火山巖熱液交代蝕變主要分為黏土化和非黏土化蝕變兩種(表1)。在水和氣等熱液作用下，火山巖中的橄欖石常蝕變?yōu)橐晾头浇馐?，斜長(zhǎng)石蝕變?yōu)楦邘X石，輝石和角閃石蝕變?yōu)榫G泥石，而綠泥石及方解石還可進(jìn)一步蝕變?yōu)檠趸F等[4]。綜合來(lái)看，火山巖的蝕變具有兩面性：一方面可能有利于成儲(chǔ)，另一方面可能會(huì)破壞儲(chǔ)層孔隙的有效性。因此，對(duì)于蝕變火山巖儲(chǔ)層分析時(shí)，要充分考慮這兩種蝕變對(duì)儲(chǔ)層造成的影響。

表1 火山巖蝕變類型

原巖發(fā)生蝕變作用的強(qiáng)度與范圍主要取決于原巖和流體的物理化學(xué)性質(zhì)。二者差異越大，蝕變作用越強(qiáng)烈，蝕變作用越廣泛。常見(jiàn)的蝕變類型很多，根據(jù)蝕變作用產(chǎn)生的主要礦物進(jìn)行分類，如綠泥石化、黑云母化、綠簾石化和絹云母化等；根據(jù)蝕變后的巖石分類，如云英巖化、矽卡巖化和青盤巖化等；以特征性的交代元素、化學(xué)組分或化合物分類，如鉀化、鈉化、硅化和碳酸鹽化等[8]。研究區(qū)火山巖發(fā)生的大都是綠泥石化，蝕變的黏土產(chǎn)物多為綠泥石。綠泥石主要是由鐵、鎂硅酸鹽礦物直接分解形成或由中、低溫?zé)嵋簬腓F、鎂組分發(fā)生交代蝕變而成。

2 火山巖蝕變層測(cè)井響應(yīng)及評(píng)價(jià)

2.1 火山巖蝕變及儲(chǔ)層空間特征

圖1是研究區(qū)內(nèi)3口井(井A、B、C)20塊火山巖巖樣(A井8塊、B井7塊、C井5塊)的黏土含量直方圖，圖2是這些巖樣的主要礦物含量直方圖。從圖1可以看出，A井巖樣的黏土含量為3.4%～36.5%，平均為18.7%；B井為12%～36.4%，平均為20.3%；C井為14.4%～26%，平均為19.1%。由此可知這3口井火山巖的黏土化程度較高。由圖2可知碳酸鹽巖(方解石)含量，A井平均為4.1%，B井為9.1%～18.5%,平均為14.1%，C井13.1%～28.1%，平均為14.8%。可見(jiàn)研究區(qū)內(nèi)火山巖黏土化明顯，碳酸巖化并不是特別普遍。研究區(qū)內(nèi)火山巖主要的礦物為斜長(zhǎng)石、石英和輝石，蝕變產(chǎn)物以黏土礦物為主，并有一定量的方解石。

圖1 研究區(qū)A、B、C井黏土含量直方圖Fig.1 Histogram of clay content in wells A, B and C in study area

圖2 礦物含量直方圖Fig.2 Histogram of mineral content

表2為X-衍射黏土礦物相對(duì)含量分析數(shù)據(jù)表，可見(jiàn)除了B井含蒙脫石外，A井不含蒙脫石，C井含伊蒙混層。伊利石普遍較少，基本<15%。各井的綠泥石相對(duì)含量較多，尤其是C井綠泥石相對(duì)含量較高，達(dá)到28%～74%。

表2 X--衍射黏土礦物分析

表3為研究區(qū)8塊巖樣掃描電鏡綜合特征以及儲(chǔ)集性能描述。由表3可見(jiàn)，A井的巖樣較致密，且多見(jiàn)蝕變的綠泥石。B井孔隙發(fā)育，綠泥石和伊利石多見(jiàn)。C井的2塊巖樣，黏土礦物主要為綠泥石、蒙脫石和孔隙，均較為發(fā)育。3口井中均發(fā)育黏土礦物，多見(jiàn)綠泥石，綠泥石是黏土化蝕變常見(jiàn)的產(chǎn)物之一，由此也可驗(yàn)證研究區(qū)蝕變作用多以黏土化為主。

2.2 蝕變火山巖儲(chǔ)層測(cè)井響應(yīng)特征

在前人研究的基礎(chǔ)上[13-17]，結(jié)合南堡地區(qū)薄片及全巖分析資料，對(duì)研究區(qū)蝕變火山巖儲(chǔ)層測(cè)井響應(yīng)特征進(jìn)行分析。測(cè)井響應(yīng)能夠反映出儲(chǔ)層多方面性質(zhì)，如儲(chǔ)層巖性、孔隙結(jié)構(gòu)和孔隙流體等。研究發(fā)現(xiàn)火山巖發(fā)生蝕變時(shí)，測(cè)井曲線特征發(fā)生變化。

電阻率測(cè)井特征一般來(lái)說(shuō)，火山巖骨架為非導(dǎo)體，電阻率測(cè)井曲線能夠指示儲(chǔ)層孔隙結(jié)構(gòu)及孔隙中流體類型。當(dāng)火山巖儲(chǔ)層發(fā)生黏土化蝕變時(shí)，某些骨架礦物蝕變生成一定量的黏土礦物，這些黏土礦物含水較多、導(dǎo)電性較強(qiáng)，對(duì)電阻率測(cè)井有附加影響[13]。此時(shí)電阻率測(cè)井值不僅反映了孔隙流體性質(zhì)，還包括了蝕變黏土的導(dǎo)電性。由于影響儲(chǔ)層導(dǎo)電性的因素復(fù)雜多樣，僅靠電阻率測(cè)井很難直接反映蝕變的影響，就本研究區(qū)而言，蝕變層較未蝕變層電阻率測(cè)井值略有減小。

表3 掃描電鏡綜合特征描述

中子測(cè)井特征中子測(cè)井是地層含氫量的直接反映。常規(guī)儲(chǔ)層中，含有氫核的油或水充填在儲(chǔ)集空間中，因此儲(chǔ)集層的含氫量多少反映了巖石孔隙度的大小?；鹕綆r蝕變層段通常會(huì)出現(xiàn)中子測(cè)井值異常高的現(xiàn)象[14]，研究發(fā)現(xiàn)火山巖發(fā)生蝕變后產(chǎn)生的綠泥石、濁沸石和絹云母等礦物中存在大量的結(jié)晶水和結(jié)構(gòu)水，這些水分的存在能夠直接導(dǎo)致中子測(cè)井響應(yīng)的高異常[15]。尤其是在發(fā)生特別嚴(yán)重的蝕變時(shí)，綠泥石化后引起中子測(cè)井值異常，使得中子測(cè)井值明顯高于同類未發(fā)生蝕變的火山巖。本研究區(qū)火山巖層段，未蝕變層的中子測(cè)井平均值≈8%，蝕變層段的均值>16%。

密度、聲波測(cè)井特征未蝕變的火山巖基質(zhì)孔不發(fā)育，骨架密度大。發(fā)生蝕變以后，基質(zhì)孔不斷發(fā)育，長(zhǎng)石及暗色礦物蝕變形成黏土，具有泥質(zhì)巖類性質(zhì)，與未蝕變的中基性火山巖比較，巖石密度降低、聲波時(shí)差增大。隨著蝕變程度的增大，儲(chǔ)層黏土含量增加，儲(chǔ)層物性、含氣性均變差[16-17]。研究區(qū)未蝕變層的密度、聲波測(cè)井均值分別為2.55 g/cm3、60 μs/ft，蝕變層為2.50 g/cm3、68 μs/ft。

ECS測(cè)井特征中基性火山巖發(fā)生蝕變會(huì)降低國(guó)際上通用的TAS模版識(shí)別火山巖巖性的可行性。根據(jù)經(jīng)薄片證實(shí)的火山巖巖性所對(duì)應(yīng)的TAS圖，表明基性玄武巖發(fā)生綠泥石化時(shí)SiO2含量增大，導(dǎo)致在TAS圖上數(shù)據(jù)點(diǎn)落到中性火山巖區(qū)域，混淆了火山巖的巖性，若不結(jié)合其他資料，則會(huì)把基性火山巖錯(cuò)誤地判斷為中性火山巖。研究區(qū)內(nèi)的C井玄武巖在TAS圖上表現(xiàn)為中性而不是基性，就是由于玄武巖蝕變導(dǎo)致SiO2含量增加引起的。

2.3 火山巖儲(chǔ)層蝕變程度評(píng)價(jià)

火山巖蝕變過(guò)程復(fù)雜多樣，蝕變后給儲(chǔ)層評(píng)價(jià)帶來(lái)的困難較大。筆者針對(duì)蝕變嚴(yán)重的A井，選取了蝕變段與未蝕變段進(jìn)行對(duì)比評(píng)價(jià)。在前人研究[4]的基礎(chǔ)上，提出了一個(gè)蝕變指數(shù)IA來(lái)反映蝕變作用的強(qiáng)弱程度：

(1)

式中：CNL為中子孔隙度(%)；RLLD是深側(cè)向電阻率(Ω·m)。

蝕變作用發(fā)生時(shí)，產(chǎn)生的黏土礦物含有大量結(jié)晶水或結(jié)構(gòu)水，使中子孔隙度增大，電阻率測(cè)井值減小，因此式(1)可以用來(lái)描述蝕變程度，蝕變指數(shù)IA越大，表明巖石蝕變程度越高；IA越小，蝕變程度越低。

由式(1)計(jì)算出A井全井段的IA值，圖3為蝕變層段與未蝕變層巖樣對(duì)IA大小的柱狀圖。通過(guò)蝕變指數(shù)的大小，可初步估計(jì)火山巖儲(chǔ)層的蝕變情況。

圖3 蝕變與未蝕變層IA柱狀圖Fig.3 Histogram of IA from altered and unaltered layers

巖芯歸位后，薄片鑒定發(fā)現(xiàn)A、B、C 3口井均存在蝕變現(xiàn)象，根據(jù)測(cè)井曲線與全巖黏土礦物含量對(duì)比分析發(fā)現(xiàn)，A井多為綠泥石化，B井和C井中碳酸鹽化與黏土化交織存在，這二者對(duì)測(cè)井曲線的影響是相反的，難以區(qū)分，因此本文提取A、B、C井中出現(xiàn)蝕變的部分巖樣的測(cè)井?dāng)?shù)據(jù)，將自然伽馬、電阻率、聲波、密度和中子5條曲線與實(shí)際的黏土礦物總量做多元回歸計(jì)算，定量計(jì)算出黏土含量：

Vcl=f(GR,RLLD,AC,DEN,CNL)

(2)

根據(jù)式(2)，可以在已知測(cè)井?dāng)?shù)據(jù)的情況下計(jì)算出黏土礦物總量。作出預(yù)測(cè)值與實(shí)際黏土分析值交會(huì)圖(圖4)，相關(guān)系數(shù)>90%，表明回歸效果良好。

圖4 預(yù)測(cè)黏土礦物Fig.4 Prediction of clay mineral content

由式(2)計(jì)算A、B、C 3口井全井段的預(yù)測(cè)黏土含量，通過(guò)預(yù)測(cè)黏土含量的高低，可以初步判斷發(fā)生蝕變的層段。計(jì)算研究區(qū)全部巖樣的IA，根據(jù)IA的大小與計(jì)算得到的預(yù)測(cè)黏土含量，再與巖芯薄片和掃描電鏡等資料結(jié)合，更準(zhǔn)確地劃分出研究?jī)?chǔ)層的蝕變層段。

3 研究區(qū)蝕變火山巖儲(chǔ)層飽和度評(píng)價(jià)

采用阿爾奇公式計(jì)算飽和度:

(3)

式中：SW為含水飽和度；RW為地層水電阻率；Rt為地層電阻率；φ為孔隙度；a為巖性參數(shù)；m為膠結(jié)指數(shù)；n為飽和度指數(shù)。

孔隙度采用聲波測(cè)井?dāng)?shù)據(jù)計(jì)算得到:

(4)

式中：AC、ACma、ACf、ACsh分別為目的層聲波時(shí)差測(cè)井值、巖石骨架聲波時(shí)差、地層流體聲波時(shí)差和黏土聲波時(shí)差；Vcl為黏土含量，由式(2)計(jì)算得到。

阿爾奇參數(shù)a、m均采用巖電實(shí)驗(yàn)測(cè)量結(jié)果。在確定飽和度指數(shù)n時(shí)，由于研究區(qū)火山巖蝕變多為黏土化蝕變，產(chǎn)生的蝕變礦物如綠泥石、伊利石等潤(rùn)濕性為油濕，而油濕礦物的存在，會(huì)導(dǎo)致儲(chǔ)層的飽和度指數(shù)n增大。因此，在計(jì)算飽和度時(shí)，蝕變區(qū)與未蝕變區(qū)的n取值不同。

根據(jù)巖電實(shí)驗(yàn)數(shù)據(jù)，做含水飽和度SW與電阻率增大系數(shù)I的交會(huì)圖(圖5)，擬合后得到蝕變層段(a)的飽和度指數(shù)n=2.936，未蝕變層段(b)的n=1.623。根據(jù)式(3)，分別計(jì)算得到蝕變層段與未蝕變層段的含水飽和度。

圖5 蝕變(a)與未蝕變(b)層段n值確定Fig.5 Determination of n in altered (a) and unaltered(b) layers

4 結(jié)論

(1)研究區(qū)存在火山巖蝕變的現(xiàn)象。蝕變后的火山巖對(duì)測(cè)井響應(yīng)有著明顯的影響?；鹕綆r蝕變產(chǎn)生的綠泥石和伊利石等黏土礦物內(nèi)部有大量結(jié)晶水存在，導(dǎo)致中子孔隙度增大，聲波時(shí)差增大，電阻率減小，密度減小。

(2)通過(guò)自然伽馬、電阻率、聲波、密度和中子5條曲線與實(shí)際的黏土礦物總量做回歸計(jì)算預(yù)測(cè)蝕變黏土含量，得到的預(yù)測(cè)黏土含量與實(shí)際含量相關(guān)性>90%。通過(guò)分析蝕變作用對(duì)測(cè)井資料的影響，定義了蝕變指數(shù)IA，該指數(shù)能判斷蝕變程度的大小，蝕變指數(shù)IA越大，表明蝕變程度越高，黏土含量越高；蝕變指數(shù)IA越小，則說(shuō)明儲(chǔ)層蝕變程度越低，黏土含量越低。

(3)蝕變后的儲(chǔ)層飽和度指數(shù)n>未發(fā)生蝕變時(shí)的n，這是由于火山巖儲(chǔ)層發(fā)生蝕變時(shí)生成了濕潤(rùn)性為油濕的蝕變產(chǎn)物，導(dǎo)致n值增大。通過(guò)巖石物理數(shù)據(jù)得到蝕變與未蝕變時(shí)的n值，利用阿爾奇公式評(píng)價(jià)研究區(qū)儲(chǔ)層的飽和度，通過(guò)對(duì)比飽和度曲線發(fā)現(xiàn)，將蝕變對(duì)n值的影響考慮到計(jì)算飽和度中是十分必要的，若忽略其影響，所獲得的飽和度誤差較大。

圖6 A井未蝕變層段測(cè)井解釋圖Fig.6 Logging interpretation of unaltered layer in Well A

圖7 A井蝕變層段測(cè)井解釋圖Fig.7 Logging interpretation of alteration layer in Well A