摘 要

摘　要：火成巖儲層具有巖性復(fù)雜、儲集空間復(fù)雜、孔隙結(jié)構(gòu)復(fù)雜和蝕變程度高等特征，使得其聲波、中子、密度、核磁和地層元素測井的響應(yīng)與沉積巖有較大的差別，通常適用于沉積巖基

摘　要：火成巖儲層具有巖性復(fù)雜、儲集空間復(fù)雜、孔隙結(jié)構(gòu)復(fù)雜和蝕變程度高等特征，使得其聲波、中子、密度、核磁和地層元素測井的響應(yīng)與沉積巖有較大的差別，通常適用于沉積巖基質(zhì)孔隙度計算的測井曲線并不完全適應(yīng)于火成巖儲層，為此進行了深入研究。結(jié)果表明：火成巖儲層具有孔徑差異大，孔隙分布不均勻、各向異性強的顯著特征，理論分析得知該類儲層的聲波測井值往往偏低，同時，大量實際測井資料也揭示聲波對火成巖巖性反應(yīng)不敏感；火成巖蝕變程度高，蝕變后的一些礦物含氫指數(shù)高，使補償中子測井值顯著增加；對理論模擬、測井資料和實驗資料進行綜合分析后認為，火成巖的強磁化率特征使巖石孔隙內(nèi)部產(chǎn)生強磁場梯度，使T2譜前移明顯、幅度降低，核磁孔隙度明顯偏低；地層元素測井孔隙度在部分井誤差大，反映該方法在火成巖地層中使用的局限性；而密度測井卻能有效反映地層的總孔隙度，實際資料顯示密度孔隙度與巖心分析孔隙度誤差較小。因此，火成巖儲層基質(zhì)孔隙度的計算應(yīng)首選密度測井曲線。

關(guān)鍵詞：火成巖  儲集層  孔隙度  密度測井  測井曲線  優(yōu)選

The preferred logging curve for the calculation of matrix porosity of volcanic reservoirs

Abstract：Volcanic rocks are characterized by complex lithologies，reservoir space and pore structure as well as strong alteration，making their responses of sonic logging，neutron logging，density logging，NMR logging and strata element logging quite different from those of sedimentary rocks．Such logs commonly suitable for the calculation of matrix porosity of sedimentary rocks are not completely suitable for volcanic reservoirs．Volcanic reservoirs have large variations of pore diameters，uneven distribution of pores and strong heterogeneity．Theoretical analysis shows that their sonic logging porosity is generally lower than the normal value．In addition，a large amount of real logging data also shows that sonic wavc is not sensitive to the lithology of volcanic rocks．Volcanic rocks are high in alteration and some altered minerals have a high hydrogen index(HI)，causing the significant increase of compensated neutrons．Theoretical simulation and comprehensive analyses of logging data and experimental data reveal that the strong magnetie susceptibility of volcanic rocks generates a strong magnetic field gradient within pores，resulting in significant forward moving of T2 spectrum and lowering of T：amplitude，thus their NMR porosity is significantly lower than the normal value．Great errors of the strata element logging porosity were found in some wells，indicating this method¢s limitation in application to volcanic rocks．In contrast，the density logging can effectively reflect the gross porosity of reservoirs，and real data show that the density logging porosity matches well with the porosity data by core analysis．Therefore，the density log should be the first choice for matrix porosity calculation of volcanic rocks．

Keywords：volcanic rock，reservoir，porosity，density log，logging curve，optimization

火成巖基質(zhì)孔隙度的計算不能完全照搬沉積巖的計算方法，適應(yīng)于計算沉積巖孔隙度的測井曲線并不一定完全適用于火成巖，兩者之間存在較大的差異。優(yōu)選計算火成巖基質(zhì)孔隙度的有利曲線是可靠評價火成巖基質(zhì)孔隙度的關(guān)鍵。研究表明：在火成巖儲層中，密度曲線計算火成巖孔隙度要明顯的優(yōu)于核磁^[1-2]、聲波^[3-5]、中子和地層元素測井(ECS) ^[6]曲線。

1　各類測井方法在火成巖中的響應(yīng)特征

1．1　聲波測井

聲波測井能有效反映孔隙均勻分布的均質(zhì)地層的孔隙度，而火成巖儲層孔隙孔徑差別較大，孔隙分布不均勻，孔隙分布各向異性強，連通性差(圖1)。如準噶爾盆地陸東地區(qū)和克拉美麗氣田的火成巖儲層孔隙分布不均勻，如黃驊坳陷北堡地區(qū)火成巖儲層中普遍發(fā)育有連通性較差的氣孔和杏仁體內(nèi)孔^[7]，松遼盆地慶深氣田儲層中發(fā)育有孔喉比大、孔喉分選差的儲層，可見火成巖儲層的孔隙分布往往不均勻^[8]。聲波測井記錄的是聲波傳播最快的首波，對于這種孔隙分布不均勻的儲層，聲波測井記錄的是繞過部分孔隙傳播速度最快的波，因為沿骨架傳播的聲波速度最快，故聲波測井在這類儲層中測得的速度往往偏低；此外，實測資料也顯示，從酸性巖到基性巖類，隨著鐵鎂礦物含量的增加，密度測井值和補償中子值是逐漸增大，且不同巖類間差別大，但聲波總體變化不大，差別較小。圖2是某地區(qū)基性熔巖、酸性熔巖和中性火山角礫巖的聲波時差對比圖，3種不同的巖性且孔隙度存在一定差異，但3種巖性的聲波時差差別較小。準噶爾盆地l97口井不同巖性的火成巖聲波時差統(tǒng)計也顯示，不同巖性火成巖的聲波時差無明顯的差異。大慶油田深層不同火山熔巖的聲波測井值也無明顯的差別。因此，不宜優(yōu)選聲波測井資料來計算火成巖地層的孔隙度。

 

 

1．2　中子測井

火成巖的蝕變程度高，蝕變后常常會形成高含結(jié)晶水的綠泥石、高嶺石、云母、蛇紋石和方沸石等，如綠泥石的含氫指數(shù)可達52％，高嶺石的含氫指數(shù)可達36％，云母的含氫指數(shù)在20％以上，因而蝕變后的一些礦物對補償中子測井值影響很大，會使補償中子顯著增加，一般是蝕變越嚴重，補償中子也越大。因此，在不能有效定量評價火成巖蝕變程度的情況下，利用補償中子計算的孔隙度往往不可靠，故不建議用補償中子來計算。

1．3　核磁共振測井

火成巖巖性復(fù)雜，巖性識別較難，這給孔隙度的計算帶來了很大困難，而“核磁共振測井不受巖性影響”這一廣泛認識，使得一些研究人員想到利用這一優(yōu)勢來計算火成巖地層的孔隙度，因為這樣可以避開火成巖復(fù)雜巖性識別這一難題，已發(fā)表了一些不考慮火成巖巖性直接使用核磁來計算火成巖儲層孔隙度的文獻。結(jié)合核磁共振測井響應(yīng)理論與實際測井資料研究表明，火成巖巖性會對核磁共振測井信號產(chǎn)生顯著影響，核磁共振測井在火成巖儲層中使用具有很大的局限性，尤其是在基性巖和低孔火成巖儲層中。核磁共振測井T2譜受自由弛豫、表面弛豫和擴散弛豫3種弛豫機制的作用。自由弛豫(體弛豫)與巖石空隙中的流體體積、流體類型、流體性質(zhì)等因素有關(guān)，表面弛豫是表面弛豫強度與孔隙比表面的乘積，其值與巖石顆粒比表面大小、膠結(jié)物性質(zhì)等有關(guān)，但不受溫度及壓力的影響。可見自由弛豫和表面弛豫都與巖性關(guān)系小。擴散弛豫是由梯度場中分子擴散引起的，而梯度場一是來自于測井儀器建立的，稱外部梯度磁場，另一個來源是巖石骨架顆粒與孔隙流體之間磁化率差(Dc)引起的內(nèi)部梯度磁場。若巖石顆粒(骨架)呈逆磁性時，油和水呈弱逆磁性，則巖石骨架顆粒與孔隙流體之間磁化率差便很小，孔隙內(nèi)部梯度磁場便很小，對于這類巖石，擴散弛豫與巖性關(guān)系不大，巖性對核磁共振測井T2譜影響不大。若巖石顆粒呈強順磁性時，巖石骨架顆粒與孔隙流體之間磁化率差便很大，會在孔隙內(nèi)部產(chǎn)生強內(nèi)部梯度磁場，將對核磁共振測井產(chǎn)生很大影響，對于這類巖石，巖性會對核磁共振測井產(chǎn)生較大影響。

而火成巖大多呈強順磁性，會使火成巖儲層孔隙內(nèi)部產(chǎn)生強梯度磁場，這是因為火成巖中含大量磁鐵礦、角閃石和黑云母等強順磁性礦物和鐵磁性礦物。從酸性火成巖到基性巖，鐵磁性礦物和強順磁性礦物的含量是逐漸增大，巖石的磁化率也逐漸增大，因此，基性巖的磁化率明顯的大于酸性巖?；鸪蓭r具有的強磁化率特征會使火成巖儲層孔隙內(nèi)部產(chǎn)生強梯度磁場，在孔徑一定的情況下，巖石的磁化率越大，孔隙內(nèi)部梯度磁場越強。

1．3．1強順磁性巖石對核磁共振測井T2譜影響的理論分析

式(1)與式(2)相結(jié)合理論模擬強內(nèi)部梯度磁場對T2譜的影響，模擬兩種大小的內(nèi)部磁場強度，模擬假定外部磁場梯度一定、不考慮不同尺寸孔徑的孔隙其內(nèi)部磁場梯度的差異性及增大回波間隔引起的短弛豫信號的損失，即

 

式中T2B為體積弛豫，ms；D為擴散系數(shù)，cm²／s；G為觀測時的梯度磁場；g為旋磁比；TE為回波間隔；r2為橫向表面弛豫強度；S／V為巖石比表面；B0為外加磁場強度；r為磁場變化的距離；Dc為巖石骨架顆粒與孔隙流體之間磁化率的差。

對比模擬的兩種大小的內(nèi)部梯度磁場對T2譜影響的成果圖可看出(圖3)：相同回波間隔下，內(nèi)部磁場梯度越強，相同弛豫組分對應(yīng)的T2譜時間越短，意味著T2譜前移越明顯，使其向弛豫時間變短的方向移動。

 

1．3．2強順磁性巖石對核磁共振測井T2譜影響的實際資料分析

實際測井資料與理論模擬結(jié)果有相同的認識。圖4所示井段為酸性火成巖類，巖性是花崗斑巖，孔隙度最大為l2％，平均孔隙度為9％，圖5所示井段為基性火成巖類，巖性是玄武巖，取心段最大孔隙度為18.6％，平均孔隙度為l2％，基性玄武巖的孔隙度要高于酸性花崗斑巖的孔隙度。

 

 

從孔喉大小看，圖5所示的基性玄武巖的孔喉半徑(圖5-b)要明顯的大于酸性花崗斑巖(圖4-b)?；孕鋷r孔喉半徑大于0.148mm的孔喉體積要明顯的大于酸性花崗斑巖，孔喉半徑大于0.148mm的孔喉體積兩口井分別是48％和l8.2％，孔喉半徑大小差異明顯。

根據(jù)核磁共振測井理論知，核磁共振測井T2譜時間的長短反映孔喉半徑的大小，一般孔喉半徑越大，T2譜弛豫時間越長，T2譜越靠后^[9-12]，但圖5的基性玄武巖盡管其孔隙度大、孔喉半徑也大，但其T2譜相對酸性花崗斑巖明顯偏向短弛豫時間，基性玄武巖的T2譜普遍分布在l0ms的左邊(圖5-a)，而酸性花崗斑巖盡管其孔喉半徑小，但其T2譜卻普遍分布在10ms的右邊(圖4-a)，差異顯著，這與一般的核磁理論認識不一致，這是因為基性玄武巖具有強順磁性，會使儲層孔隙內(nèi)部產(chǎn)生強梯度磁場，使擴散弛豫明顯增強，使T2譜發(fā)生明顯前移的現(xiàn)象，這與前面的關(guān)于強順磁性巖石對核磁共振測井T2譜影響的理論分析相一致。

1．3．3內(nèi)部梯度磁場的大小對核磁孔隙度的影響

火成巖強磁化率特征會使儲層孔隙內(nèi)部產(chǎn)生強梯度磁場，使擴散弛豫明顯增強，除了使T2譜發(fā)生前移外，還會使核磁T2譜發(fā)生衰減，譜面積減小，核磁孔隙度降低，核磁信號的衰減程度與孔隙內(nèi)部的梯度磁場大小成正比，孔隙內(nèi)部梯度磁場越大，核磁信號衰減越大，核磁孔隙度越小。

6塊巖心核磁實驗孔隙度與常規(guī)分析孔隙度對比看出，火成巖核磁實驗孔隙度明顯低于常規(guī)分析的孔隙度(表1)，并且基性巖(21號巖樣)減小的程度要明顯的大于酸性巖(如表l中的12、22、23號巖樣)。這是因為，基性玄武巖的磁化率高，孔隙內(nèi)部梯度磁場強，核磁信號衰減幅度大。

 

實際測井資料也具有相同的規(guī)律。圖5中基性玄武巖取心段巖心分析平均孔隙度ll.93％，核磁分析平均孔隙度8.15％，平均誤差31.7％。圖4中酸性花崗斑巖巖心分析平均孔隙度8.45％，核磁分析平均孔隙度7.97％，平均誤差5.7％。表明強磁化率巖石核磁分析的孔隙度較常規(guī)分析的偏低，巖石磁化率越大，核磁分析的孔隙度偏低幅度越大。

由式(3)可知，巖石孔隙內(nèi)部磁場梯度的強弱除與巖石磁化率有關(guān)外，還與巖石孔喉尺寸大小有關(guān)，與孔喉尺寸大小成反比，孔喉半徑越小，孔隙內(nèi)部梯度磁場越強，T2譜前移越明顯，核磁測井孔隙度也越低。

圖5中3641～3645m段核磁孔隙度與常規(guī)巖心孔隙度之間的誤差要明顯的大于3636～364lm段(圖5-a)，這是因為3641～3645m段的巖心孔喉半徑要明顯的小于3636～3641m段。

因此，利用核磁共振測井計算火成巖地層的孔隙度具有很大的局限性，尤其在中基性火成巖地層和低孔地層中^[13-14]。

1．4　密度測井

密度測井反映的是儀器探測范圍內(nèi)地層的總密度，能有效反映地層的總孔隙度，對各種類型的孔隙度都能有效反映，如火成巖地層中大的氣孔、溶蝕孔洞空隙都能反映。采用巖心刻度測井的方法，分區(qū)塊、分巖性建立了不同區(qū)塊、不同巖性的密度測井與巖心分析孔隙度間的線性回歸模型用于計算火成巖儲層的孔隙度，實際應(yīng)用表明，利用密度曲線計算的孔隙度與巖心分析孔隙度之間的誤差較小(表2)。

 

1．5　ECS測并

地層元素測井(ECS)可以獲得主要造巖元素Si、Fe、Ti、Ca、AL、S、Cl、Cr、Gd等的質(zhì)量百分含量，這些元素的含量與巖石的骨架密度直接相關(guān)，中國石油大慶油田有限責任公司、中國石油新疆油田公司(以下簡稱新疆油田公司)^[6]等建立了利用ECS測井測得的元素含量直接計算巖石體積密度的模型，可以不考慮火成巖復(fù)雜巖性識別的問題。新疆油田公司應(yīng)用準噶爾盆地18口井319組實驗數(shù)據(jù)來計算火成巖的骨架密度^[6]：

rma＝2.53553+0.10462Wsi+0.40365WFe+13.619977WTi      (4)

jD＝(rma-rb)／(rma-rf)        (5)

式中Wsi WFe、WTi分別為ECS硅元素、鐵元素和鈦元素的質(zhì)量百分含量；jD為密度孔隙度；rma、rb、rf分別為骨架密度、測井密度、流體密度值。

利用式(4)可動態(tài)計算不同深度點的巖石骨架密度，然后將骨架密度值帶入式(5)逐深度點計算火成巖地層的孔隙度。從8口井的實際計算效果看，ECS測井方法確定的巖石骨架密度帶入孔隙度模型計算的孔隙度在個別井誤差大(表2)，顯示該方法還不能完全適應(yīng)火成巖地層，該方法在確定復(fù)雜火成巖巖石骨架方面有許多不完善的地方，這也是因火成巖礦物組分多、礦物組分復(fù)雜、巖性復(fù)雜引起的。

2　結(jié)束語

綜合研究成果表明：用于火成巖孔隙度計算的聲波、中子、密度、核磁和ECS曲線中，應(yīng)首先優(yōu)選密度測曲線來計算火成巖儲層的孔隙度。

 

參考文獻

[1]張玲，王志強，張麗艷，等．火成巖油氣藏儲量計算有關(guān)問題探討[J]．石油與天然氣地質(zhì)，2009，30(2)：223-229．

ZHANG Ling，WANG Zhiqiang，ZHANG Liyan，et al．A discussion on certain issues concerning reserve calculation of igneous reservoirs[J]．Oil and Gas Geology，2009，30(2)：223-229．

[2]趙杰，雷茂盛，楊興旺，等．火山巖地層測井評價新技術(shù)[J]．大慶石油地質(zhì)與開發(fā)，2007，26(6)：134-137．

ZHAO Jie，LEI Maosheng，YANG Xingwang，et al．New logging evaluation technique for volcanic rock formation[J]．Petroleum Geology&Oilfield Development in Daqing．2007，26(6)：134-137．

[3]李寧，付有升，楊曉玲，等．大慶深層流紋巖全直徑巖心實驗數(shù)據(jù)分析[J]．測井技術(shù)，2005，29(6)：480-483．

LI Ning，FU Yousheng，YANG Xiaoling，et al．Experimental data analysis on full diameter rhyolite cores from Daqing deep layers[J]．Well Logging Technology，2005，29(6)：480-483．

[4]劉俊田，焦立新，張日供，等．三塘湖盆地牛東區(qū)塊石炭系火山巖儲集層測井評價[J]．天然氣地球科學(xué)，2009，20(4)：518-521．

LIU J untian，JIAO Lixin，ZHANG Rigong，et al．Well logging evaluation of Carboniferous volcanic reservoirs in Niudong Block，Santanghu basin[J]．Natural Gas Geoscience，2009，20(4)：518-521．

[5]陳鋼花，范宜仁，代詩華．火山巖儲層測井評價技術(shù)[J]．中國海上油氣：地質(zhì)，2000，14(6)：422-427．

CHEN Ganghua，FAN Yiren，DAI Shihua．An evaluation technique for volcanic reservoir using well logging data[J]．China Offshore Oil and Gas(Geology)，2000，14(6)：422-427．

[6]匡立春，董政，孫中春，等．基于元素俘獲譜測井計算火山巖儲集層孔隙度的方法[J]．新疆石油地質(zhì)，2009，30(3)：287-288．

KUANG Lichun，DONG Zhen9，SUN Zhongchun，et al．Method for calculation of volcanic reservoir porosity basedon ECS log[J]．Xingjiang Petroleum Geology，2009，30(3)：287-288．

[7]馬乾，鄂俊杰，李文華，等．黃驊坳陷北堡地區(qū)深層火成巖儲層評價[J]．石油與天然氣地質(zhì)，2000，21(4)：337-340．

MA Qian，E Junjie，LI Wenhua，et al．Reservoir evaluation of deep seated igneous rocks in Beipu Region，Huanghua Depression[J]．Oil＆Gas Geology，2000，21(4)：337-340．

[8]曲延明，舒萍，紀學(xué)雁，等．松遼盆地慶深氣田火山巖儲層的微觀結(jié)構(gòu)研究[J]．吉林大學(xué)學(xué)報：地球科學(xué)版，2007，37(4)：721-725．

QU Yanming，SHU Ping，JI Xueyan，et al．Micro fabrics of reservoir volcanic rocks in the Qingshen Gas Field of the Songliao Basin[J]．Journal of Jilin University：Earth Science Edition，2007，37(4)：721-725．

[9]趙輝，戴詩華，司馬直強，等．磨溪地區(qū)嘉二段儲層有效性評價方法[J]．西南石油大學(xué)學(xué)報：自然科學(xué)版，2009，31(5)：13-16．

ZHAO Hui，DAI Shihua，SIMA Liqiang，et al．Reservoir effectiveness evaluation for Tlj2 strata in Moxi Area[J]．Journal of Southwest Petroleum University：Science&Technology Edition．2009，31(5)：13-16．

[l0]顏其彬，趙輝，司馬立強，等．碳酸鹽巖核磁共振實驗與認識[J]．天然氣工業(yè)，2010，30(1)：36-38．

YAN Qibin，ZHAO Hui，SIMA Liqiang，et al．A studv of NMR experiments of carbonates[J]．Natural Gas Industry，2010，30(1)：36-38．

[11]司馬立強，吳豐，繆祥禧．龍崗地區(qū)礁灘儲層有效性評價[J]．西南石油大學(xué)學(xué)報：自然科學(xué)版，2012，34(1)：176-182．

SIMA Liqiang，WU Feng，MIAO Xiangxi．Bank reef reservoir effectiveness eValuation in Longgang District[J]．Journal of Southwest Petroleum Universitv：Science&Technology Edition，2012，34(1)：l76-l82．

[12]陳歡慶，胡永樂，冉啟全，等．徐東地區(qū)營城組一段火山噴發(fā)模式特征[J]．西南石油大學(xué)學(xué)報：自然科學(xué)版，2012，34(1)：41-48．

CHEN Huanqing，HU Yongle，RAN Qiquan，et al．Mode characteristics of volcanic eruption of the member l of Yingcheng Formation in Xudong Area in Songliao Basin[J]．Journal of Southwest Petroleum University：Science＆Technology Edition，2012，34(1)：41-48．

[13]曹寶軍，李-卡H7y，姚約東．火山巖氣藏開發(fā)難點與對策[J]．天然氣工業(yè)，2007，27(8)：82-84．

CAO Baojun，LI Xiangfang，YAO Yuedong．Challenges and strategies in volcanic gas reservoir development[J]．Natural Gas Industry，2007，27(8)：82-84．

[14]陳曉紅，何文淵，馮子輝．松遼盆地徐家圈子斷陷主要斷裂對氣藏的控制作用[J]．天然氣工業(yè)，2012，32(3)：53-56．

CHEN Xiaohong，HE Wenyuan，FENG Zihui．Controlling effect of major faults on the gas reservoirs in tbe Xujia weizi fault depression，Songliao Basin[J]．Natural Gas Industry，2012，32(3)：53-56．

 

本文作者：戴詩華  趙輝  姜淑云

作者單位：中國石油西部鉆探工程公司

  中國石油川慶鉆探工程公司地質(zhì)勘探開發(fā)研究院

  中國石油長城鉆探工程公司服務(wù)公司