王翼君，崔剛，唐洪明，沈建軍，高子涵，周士琳

（1.西南石油大學(xué)地球科學(xué)與技術(shù)學(xué)院，四川 成都 610500；2.中海石油（中國(guó)）有限公司天津分公司，天津 300452；3.中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司，河北 涿州 072751）

碳酸鹽巖核磁共振實(shí)驗(yàn)研究現(xiàn)狀

王翼君1，崔剛2，唐洪明1，沈建軍1，高子涵3，周士琳1

（1.西南石油大學(xué)地球科學(xué)與技術(shù)學(xué)院，四川 成都 610500；2.中海石油（中國(guó)）有限公司天津分公司，天津 300452；3.中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司，河北 涿州 072751）

隨著海相碳酸鹽巖油氣田的大規(guī)模開發(fā)，對(duì)碳酸鹽巖儲(chǔ)層的評(píng)價(jià)日益受到重視。低磁場(chǎng)核磁共振（NMR）巖心分析技術(shù)是近20年來(lái)在我國(guó)新興起的快速測(cè)量巖石物性參數(shù)的一種新技術(shù)，被廣泛應(yīng)用于儲(chǔ)層評(píng)價(jià)及流體識(shí)別。針對(duì)NMR當(dāng)前應(yīng)用現(xiàn)狀，文中就碳酸鹽巖中存在的大、小孔喉“擴(kuò)散耦合”、T2截止值的準(zhǔn)確求取以及溫度對(duì)NMR響應(yīng)特征的影響3個(gè)問(wèn)題，總結(jié)了國(guó)內(nèi)外相關(guān)學(xué)者的研究成果。“擴(kuò)散耦合”在各類巖石中都可能存在，只不過(guò)是在碳酸鹽巖中較為突出罷了；碳酸鹽巖的T2截止值普遍比砂巖高且變化較大，這一方面說(shuō)明碳酸鹽巖的非均質(zhì)性較強(qiáng)，另一方面也可能與碳酸鹽巖的“擴(kuò)散耦合”相關(guān)；而溫度對(duì)碳酸鹽巖NMR響應(yīng)特征的影響，則因碳酸鹽巖的礦物成分、孔隙類型、流體性質(zhì)而異。

核磁共振；擴(kuò)散耦合；T2截止值；溫度；表面弛豫強(qiáng)度

0 引言

近年來(lái)，我國(guó)的海相碳酸鹽巖油氣勘探取得了重大突破，先后發(fā)現(xiàn)了普光氣田、輪南-塔河油田等大型海相碳酸鹽巖油氣田，顯示出海相碳酸鹽巖的巨大勘探開發(fā)潛力，對(duì)我國(guó)石油工業(yè)的發(fā)展和能源結(jié)構(gòu)的改善意義重大［1-2］。在此背景下，對(duì)于碳酸鹽巖儲(chǔ)層的綜合評(píng)價(jià)也成了擺在地質(zhì)學(xué)家面前的一個(gè)重要課題。其中，低磁場(chǎng)核磁共振技術(shù)（NMR）作為一種評(píng)價(jià)手段也逐漸被應(yīng)用于碳酸鹽巖儲(chǔ)層的評(píng)價(jià)。

作為一種新興技術(shù)，NMR巖心分析技術(shù)，利用儲(chǔ)層流體（油/氣/水）中的氫原子核在磁場(chǎng)中具有共振并產(chǎn)生信號(hào)的特性，探測(cè)儲(chǔ)層孔隙度、滲透率、可動(dòng)流體飽和度等巖石物性參數(shù)以及巖心中的流體在孔喉中的分布，是有效劃分和評(píng)價(jià)儲(chǔ)層，識(shí)別儲(chǔ)層流體的重要手段，具有測(cè)試方便、快捷、操作簡(jiǎn)單、無(wú)污染等特點(diǎn)［3-5］。但是，由于種種原因，目前對(duì)于碳酸鹽巖儲(chǔ)層核磁共振特征的研究較之碎屑巖儲(chǔ)層滯后。本文總結(jié)了國(guó)內(nèi)外對(duì)碳酸鹽巖中存在的大小孔喉中流體的 “擴(kuò)散耦合”、T2截止值的求取以及溫度對(duì)碳酸鹽巖的影響等3個(gè)方面的研究成果，并且在此基礎(chǔ)上提出筆者的一些認(rèn)識(shí)，希望能為今后碳酸鹽巖的NMR特征研究以及用NMR準(zhǔn)確評(píng)價(jià)碳酸鹽巖儲(chǔ)層提供一些參考。

1 NMR原理概述

NMR技術(shù)研究的是“淹沒(méi)”在恒定磁場(chǎng)中的某些帶自旋磁矩的原子核吸收特定頻率 （Larmor頻率）的電磁波，從而改變能量狀態(tài)的現(xiàn)象，其基礎(chǔ)是具有自旋性質(zhì)的原子核的磁性及其與外加磁場(chǎng)的相互作用。在對(duì)飽和含氫流體（水/油）的巖心進(jìn)行核磁共振測(cè)試之前，地層中的氫質(zhì)子隨機(jī)取向和排列。將飽和含氫流體（水/油）的巖心放入恒定靜磁場(chǎng)B0中時(shí)，帶自旋的原子核在B0磁場(chǎng)的作用下，自旋軸沿B0方向排列，且具有2種取向和能級(jí)：與B0方向相同，處于高能級(jí)；與B0方向相反，處于低能級(jí)。

置于B0中的自旋原子核，由于處于高能級(jí)和低能級(jí)的數(shù)目不同，宏觀上會(huì)產(chǎn)生一個(gè)凈磁化矢量，稱為宏觀磁化矢量 （M），M的方向與外加磁場(chǎng)B0方向一致，原子核系統(tǒng)被“極化”。當(dāng)在垂直于B0的方向再施加一個(gè)頻率與Larmor頻率相等的交變磁場(chǎng)B1，“照射”被極化了的核自旋系統(tǒng)時(shí)，系統(tǒng)將吸收B1場(chǎng)的能量，使得M被“激發(fā)”，偏離B0方向。B1場(chǎng)結(jié)束后，自旋將逐步釋放或交換能量，核自旋系統(tǒng)恢復(fù)到施加B1前的平衡狀態(tài)［6-7］。這一過(guò)程，稱之為“弛豫”（包括縱向弛豫和橫向弛豫2個(gè)過(guò)程），這一過(guò)程需要的時(shí)間稱之為“弛豫時(shí)間”（包括縱向弛豫時(shí)間T1和橫向弛豫時(shí)間T2），目前，低場(chǎng)核磁共振的測(cè)試最常用的是T2。

儲(chǔ)層原油、天然氣和水中富含氫核，是NMR的主要研究對(duì)象，M與儲(chǔ)層巖心中油氣水的氫核含量成正比，這是核磁共振通過(guò)測(cè)量巖心中的流體信號(hào)得到孔隙度、滲透率等物性參數(shù)和可動(dòng)流體飽和度等流體參數(shù)的理論基礎(chǔ)。

2 碳酸鹽巖NMR特征研究

通過(guò)NMR法評(píng)價(jià)儲(chǔ)層好壞，除了核磁孔隙度和滲透率外，主要還依據(jù)2點(diǎn)：1）儲(chǔ)層巖心的T2截止值及其平均值大?。?）儲(chǔ)層巖心的T2時(shí)間大小。此外，對(duì)于現(xiàn)場(chǎng)的核磁共振測(cè)井而言，溫度對(duì)于核磁共振T2分布的影響也至關(guān)重要。

2.1 碳酸鹽巖的T2時(shí)間和孔喉分布的關(guān)系

T2時(shí)間可以用來(lái)表征巖石的孔徑分布，其理論依據(jù)是：在NMR測(cè)量過(guò)程中，氫質(zhì)子碰撞顆粒表面造成能量衰減，這一過(guò)程稱為表面弛豫。其中：大孔隙的比表面（S/V，單位是μm-1）小，氫質(zhì)子和顆粒表面發(fā)生碰撞概率小，弛豫時(shí)間較長(zhǎng)；而小孔喉的S/V大，因而氫質(zhì)子與顆粒表面發(fā)生碰撞概率大，弛豫時(shí)間短。

巖石顆粒表面的弛豫能力往往用表面弛豫強(qiáng)度（ρ2，又稱為表面弛豫率，μm/s）來(lái)定量表征。表面弛豫速率（1/T2，單位為s-1）是表面弛豫強(qiáng)度與巖石孔隙比表面的乘積：1/T2=ρ2（S/V）。由此，可以得出，巖石的T2時(shí)間和孔喉半徑成正比，T2時(shí)間的大小可以在一定程度上反映巖石的孔喉半徑大小。

在核磁共振中，將巖心的弛豫過(guò)程表示為

式中：M（t）和M0分別為在t時(shí)刻和0時(shí)刻的磁化矢量的幅度，a.u（任意單位）；t和Tj分別為采樣時(shí)間和弛豫時(shí)間，ms；fj為Tj時(shí)的孔隙體積，mL。

這一理論，在碎屑巖儲(chǔ)層中獲得了較好的應(yīng)用效果。但是，對(duì)于碳酸鹽巖儲(chǔ)層，該觀點(diǎn)遇到了挑戰(zhàn)。這是因?yàn)椋肨2時(shí)間來(lái)表征孔喉半徑還有一個(gè)大前提，即自旋原子核擴(kuò)散距離可以忽略不計(jì)，而每一個(gè)孔喉中流體的弛豫過(guò)程是相互獨(dú)立的，彼此之間互不干擾。但是，根據(jù)國(guó)外最近10多年來(lái)的研究表明，碳酸鹽巖中存在一個(gè)“擴(kuò)散耦合”效應(yīng)，即不同半徑孔喉之間的流體通過(guò)擴(kuò)散作用相互聯(lián)結(jié)，導(dǎo)致同一個(gè)碳酸鹽巖樣品中大、小孔喉的T2時(shí)間分布很難分清。有時(shí)甚至?xí)霈F(xiàn)不同孔喉的T2譜完全合并成一個(gè)尖峰，因而不能準(zhǔn)確反映碳酸鹽巖的孔喉分布。

Roberts等［8］研究表明，砂巖的ρ2分布范圍主要在1～10 μm/s。Chang等［9］的研究認(rèn)為，碳酸鹽巖的橫向ρ2比砂巖要低，一般為0～5 μm/s，而砂巖的ρ2一般為0～15 μm/s。Ramakrishnan［10］根據(jù)蒙特卡羅隨機(jī)游走模型模擬證明，弛豫時(shí)間能否反映孔徑分布，ρ2至關(guān)重要。如果表面弛豫微弱，大、小孔隙中的流體容易擴(kuò)散，那么巖石中的“擴(kuò)散耦合”效應(yīng)就很明顯。此時(shí)，弛豫時(shí)間就不能反映巖石的孔徑分布。Toumelin等［11］使用蒙特卡羅三維隨機(jī)步進(jìn)法分析了在恒定磁場(chǎng)梯度下微晶灰?guī)r中不同類型孔隙中流體的擴(kuò)散效應(yīng)特征，并認(rèn)為“擴(kuò)散耦合”對(duì)把測(cè)量的T2分布反演成孔隙大、小分布的影響超過(guò)確定T2截止值的影響。Vivek等［12］也認(rèn)為，傳統(tǒng)的NMR測(cè)量的前提對(duì)于粒屑灰?guī)r和泥粒灰?guī)r是行不通的。他們通過(guò)模擬巖心中微孔隙和大孔隙間流體的弛豫及他們之間的流體“擴(kuò)散耦合”現(xiàn)象，將影響T2分布的幾個(gè)參數(shù)（包括微孔隙的ρ2、流體的擴(kuò)散系數(shù)、空隙的幾何形態(tài)等），合并成了一個(gè)耦合參數(shù)α，它是孔隙特征弛豫速率和微孔—大孔間擴(kuò)散速率的比值，其數(shù)學(xué)表達(dá)式為:

式中：D為擴(kuò)散系數(shù)；L1，L2分別為孔隙半長(zhǎng)和孔隙半寬，μm；T2，C為特征孔隙的弛豫時(shí)間，ms；Vtotal為總孔隙體積，mL；Sactive為微孔隙表面積，cm2；η為孔隙縱橫比，；μ為Browns數(shù)值；β為微孔百分?jǐn)?shù)，%。

圖1為μ對(duì)t=1.0等值線的影響示意。由圖1可以看出，當(dāng)μ＜1.0時(shí)，快速擴(kuò)散在整個(gè)孔隙空間形成了幾乎均勻的磁場(chǎng)，隨著μ值的增大，縱向上的梯度越來(lái)越明顯［12］。梯度的出現(xiàn)說(shuō)明了微孔隙中自旋的弛豫速率比大孔隙中的自旋大。

圖1 t=1.0時(shí)的磁場(chǎng)強(qiáng)度等值線

Vivek等選取了不同α和μ值的T2譜進(jìn)行對(duì)比后發(fā)現(xiàn)，兩者值越小，巖心中不同大小孔隙中的流體“擴(kuò)散耦合”程度越大，T2分布就越不能準(zhǔn)確地表征巖石的孔喉分布，反之亦然（見圖2）。此外，為了定量表示巖心中孔隙流體的“擴(kuò)散耦合”程度，Vivek等［12］將“擴(kuò)散耦合”程度分為3個(gè)等級(jí)：1）α＜1.0,完全耦合。在微孔隙中，擴(kuò)散比弛豫更快，微孔和大孔全部耦合，微孔和大孔的弛豫時(shí)間譜完全重合成一個(gè)尖峰，無(wú)法區(qū)分；2）1.0≤α≤50.0，中等耦合。擴(kuò)散比弛豫快，微孔和大孔部分耦合。雖然大孔和微孔各自有它們的T2分布，但是他們的尖峰信號(hào)幅度（信號(hào)強(qiáng)度f(wàn)，a.u)并不與孔隙百分?jǐn)?shù)成正比；3）α＞50.0，不耦合。微孔和大孔的弛豫完全獨(dú)立，互不相干，它們的T2譜完全分離。

圖2 與圖1中備選參數(shù)相關(guān)的α函數(shù)的T2分布

在描述粒屑灰?guī)r的大、小孔隙耦合效應(yīng)時(shí)，為了檢驗(yàn)Ramakrishnan等提出的粒屑灰?guī)r模型（詳見2.2）的有效性，研究孔隙耦合的控制參數(shù)，Vivek等分別用不同粒徑的白堊和硅膠對(duì)飽和水及束縛水狀態(tài)下的T2分布進(jìn)行了對(duì)比，結(jié)果如圖3所示。

從圖3中可以看出：對(duì)于粗粒徑的白堊，飽和水狀態(tài)下的T2分布基本呈雙峰，且離心后的T2分布和離心前的T2譜的左峰基本重合，說(shuō)明基本上不存在“擴(kuò)散耦合”現(xiàn)象；隨著白堊粒度變細(xì)，飽水狀態(tài)下的T2分布逐漸轉(zhuǎn)化為單峰，且離心后的束縛水T2譜和離心前的T2譜出現(xiàn)了明顯的不重合，或者較之離心前的左半峰信號(hào)幅度增高，或者幅度變低，但束縛水峰出現(xiàn)了明顯左移，硅膠的情況也與之類似。這說(shuō)明，隨著巖石顆粒粒徑變小，巖石中大、小孔喉中的流體的“擴(kuò)散耦合”程度增強(qiáng)。而一般情況下，隨著巖石顆粒粒徑的減小，巖石會(huì)變得致密化，物性也會(huì)變差。盡管這一結(jié)論是針對(duì)灰?guī)r提出的，但是能否據(jù)此推論：巖石越致密，物性越差，其中不同大小孔喉中的流體越容易發(fā)生“擴(kuò)散耦合”？這還需要大量的實(shí)驗(yàn)數(shù)據(jù)來(lái)證明。

肖立志等［13］在對(duì)不同表面弛豫強(qiáng)度下的NMR響應(yīng)特征進(jìn)行研究后認(rèn)為，在連通性較好的模型中，隨著表面弛豫強(qiáng)度降低，“擴(kuò)散耦合”越來(lái)越明顯。在表面弛豫強(qiáng)度較低的模型中，T2譜已經(jīng)很難直接反映模型的孔隙結(jié)構(gòu)信息。LU Chi和Zoya［14-15］認(rèn)為，Ramakrishnan提出的“擴(kuò)散耦合”效應(yīng)理論也可延伸至孔隙-裂縫雙重介質(zhì)巖石（裂縫較為發(fā)育的巖石）中，裂縫和基質(zhì)孔隙之間也存在“擴(kuò)散耦合”效應(yīng)。他們還認(rèn)為，如果在對(duì)NMR的T2分布進(jìn)行解釋時(shí)，不考慮裂縫-孔隙的“擴(kuò)散耦合”效應(yīng)，則粒間孔或裂縫的T2時(shí)間會(huì)被低估29%，從而導(dǎo)致孔隙尺寸和裂縫寬度也被低估。

圖3 作為顆粒直徑函數(shù)的多微孔白堊和硅膠T2分布

盡管上述研究結(jié)果主要是針對(duì)碳酸鹽巖，但是筆者認(rèn)為，對(duì)于傳統(tǒng)的核磁共振巖心測(cè)試分析理論基礎(chǔ)來(lái)說(shuō)，無(wú)疑是一種顛覆。這就如同儲(chǔ)集層和蓋層組合，其本質(zhì)是由于后者物性較前者差，而與其巖性差異沒(méi)有必然聯(lián)系一樣［16-17］。從上述國(guó)外學(xué)者的研究來(lái)看，“擴(kuò)散耦合”在本質(zhì)上是由于表面弛豫強(qiáng)度低導(dǎo)致大、小孔喉之間的流體通過(guò)擴(kuò)散相互聯(lián)結(jié)，與巖性并無(wú)必然聯(lián)系。這就涉及到一個(gè)問(wèn)題，即影響巖石表面弛豫速率的因素是什么？R Dastidar等［18］研究了碎屑巖的NMR響應(yīng)之后認(rèn)為，碎屑巖中較高的石英質(zhì)量分?jǐn)?shù)對(duì)應(yīng)著較低的表面弛豫速率（見圖4）。

圖4 石英質(zhì)量分?jǐn)?shù)和表面弛豫速率相關(guān)性統(tǒng)計(jì)

K Kristina［19］研究了純石英砂樣品和有著不同質(zhì)量分?jǐn)?shù)的鍍水鐵礦膜石英砂樣品的核磁響應(yīng)后發(fā)現(xiàn)，隨著樣品中鍍水鐵礦膜石英砂的質(zhì)量分?jǐn)?shù) （CM）的增大，樣品的體弛豫速率（T2B-1）和擴(kuò)散弛豫速率（T2D-1）增加不明顯，而表面弛豫速率（T2S-1）顯著增大，且直接影響了樣品總的弛豫速率T2M-1（T2M-1=T2B-1+T2D-1+T2S-1）（見圖5，回波時(shí)間間隔TE=300 μs），即此時(shí)的樣品主要受表面弛豫的影響。此外，對(duì)這些樣品進(jìn)行比表面測(cè)試后也發(fā)現(xiàn)，隨著CM增大，樣品的比表面也會(huì)增加（純石英砂為0.398 μm-1，而完全鍍了水鐵礦膜的石英砂樣品為3.93 μm-1），這應(yīng)該是順磁性物質(zhì)引起表面弛豫速率變化的重要原因。

圖5 T2M-1與CM的相關(guān)性

此外，K Kristina還發(fā)現(xiàn)：當(dāng)樣品中的CM＜25%時(shí)，樣品總體的ρ2隨著CM的增加呈指數(shù)增加；當(dāng)CM＜2%時(shí)，鍍水鐵礦膜石英砂部分的表面弛豫強(qiáng)度ρ2m顯著增加；當(dāng)2%＜CM＜25%時(shí)，ρ2m卻呈現(xiàn)指數(shù)衰減。但是，當(dāng)CM＞25%時(shí)，樣品總體的ρ2與CM呈線性正相關(guān)且增長(zhǎng)幅度較小，且ρ2m基本為一個(gè)常數(shù)（見圖6）。這說(shuō)明鐵磁性物質(zhì)的質(zhì)量分?jǐn)?shù)較小時(shí)，巖石的ρ2增長(zhǎng)很快，但是到了一定程度之后，其增長(zhǎng)反而會(huì)陷入停滯。

圖6 ρ和CM相關(guān)性統(tǒng)計(jì)

從上述研究可以看出，巖石表面ρ2受巖石內(nèi)部的順磁-鐵磁性物質(zhì)影響很大。而不同巖石所含的順磁性物質(zhì)的類型和數(shù)量是不一樣的，因而，不同種類（甚至同種類）巖石的T2S-1也一定千差萬(wàn)別。

在三大類巖石中，火成巖和正變質(zhì)巖中含有大量的暗色鐵鎂礦物（如角閃石、輝石、黑云母等）以及磁鐵礦等，它們的磁化率較高，因而在磁法勘探中，火成巖和正變質(zhì)巖的磁異常較大。相比之下，沉積巖主要由石英、方解石、長(zhǎng)石等弱磁乃至逆磁性礦物組成，順磁-鐵磁性物質(zhì)較少，因而在磁法勘探中，通常沉積巖的磁異常很低。但如果沉積巖中有較多來(lái)自火成巖或正變質(zhì)巖碎屑的磁性顆粒，則沉積巖也會(huì)有一定的磁性（碎屑剩磁）［20-21］。同樣，在核磁共振中，火成巖和正變質(zhì)巖中的這些順磁-鐵磁性物質(zhì)對(duì)火成巖的NMR信號(hào)的影響極大，T2譜往往向短弛豫方向移動(dòng)，核磁信號(hào)幅度衰減程度也變大，小孔喉中的信號(hào)不能被檢測(cè)到，導(dǎo)致核磁孔隙度偏低［22-23］。對(duì)此，可以用上述國(guó)外學(xué)者的研究結(jié)果，即巖石中的順磁-鐵磁性物質(zhì)引起的T2S-1增大，來(lái)進(jìn)行很好地解釋。

從沉積學(xué)的角度上來(lái)看，砂巖作為一種陸源碎屑巖，它也是由地球上的火成巖、變質(zhì)巖以及先形成的沉積巖經(jīng)過(guò)風(fēng)化剝蝕、搬運(yùn)、沉積和成巖作用改造后形成的［24］，如果砂巖中含有大量的火成巖或正變質(zhì)巖巖屑而具有一定的碎屑剩磁的話，砂巖的T2S-1就會(huì)較大。反之，若砂巖的成分成熟度較高（即逆磁性的石英質(zhì)量分?jǐn)?shù)很高）的話，則砂巖的T2S-1就會(huì)變低，則砂巖發(fā)生“擴(kuò)散耦合”的可能性也會(huì)增加。對(duì)于碳酸鹽巖，由于碳酸鹽巖屬于原地成因，且它的形成環(huán)境主要為溫暖、清潔的淺水環(huán)境，不可能有太多的陸源物質(zhì)注入［25］。所以，在碳酸鹽巖中，火成巖和正變質(zhì)巖巖屑極少，因而碳酸鹽巖中的順磁-鐵磁性物質(zhì)相對(duì)于砂巖較少，碎屑剩磁較低，故它的T2S-1較低，因而更容易發(fā)生“擴(kuò)散耦合”。因此，可以這樣說(shuō)，無(wú)論是砂巖還是碳酸鹽巖（甚至更多巖類），只要它們中的順磁-鐵磁性礦物質(zhì)量分?jǐn)?shù)較低，都會(huì)使它們的表面弛豫速率降低，從而導(dǎo)致大孔和微孔之間的流體都會(huì)或多或少的出現(xiàn)“擴(kuò)散耦合”（只是在碳酸鹽巖中表現(xiàn)得更顯著罷了），使其T2譜不能準(zhǔn)確反映巖石的孔喉分布。這一點(diǎn)，值得去深入探究，以便更加精確地使用核磁法來(lái)表征碳酸鹽巖儲(chǔ)層。

目前，國(guó)內(nèi)在這方面的研究相對(duì)滯后。很多學(xué)者在用核磁法求取巖石的孔喉分布 （尤其是碳酸鹽巖的孔喉分布）時(shí)仍然堅(jiān)持傳統(tǒng)的NMR測(cè)量的基本假設(shè)，仍然認(rèn)為不同的弛豫時(shí)間在任何條件下都能對(duì)應(yīng)著不同的孔喉半徑，未考慮巖石（尤其是碳酸鹽巖）中微孔和大孔的“擴(kuò)散耦合”效應(yīng)。

2.2 碳酸鹽巖的T2截止值

在核磁共振法對(duì)儲(chǔ)層進(jìn)行評(píng)價(jià)的過(guò)程中，可動(dòng)流體飽和度的確定至關(guān)重要。目前，在低場(chǎng)核磁共振中，主要運(yùn)用T2截止值法來(lái)測(cè)量可動(dòng)流體飽和度。T2截止值法假定束縛流體主要占據(jù)小孔和微孔，可動(dòng)流體占據(jù)大孔隙。假設(shè)有一個(gè)T2時(shí)間，當(dāng)孔隙中流體的弛豫時(shí)間小于該值時(shí)，孔隙流體為束縛流體，反之，為可動(dòng)流體；因此，確定T2截止值對(duì)計(jì)算可動(dòng)流體飽和度至關(guān)重要。在實(shí)驗(yàn)室，一般采用離心法標(biāo)定T2截止值。

關(guān)于確定碳酸鹽巖的T2截止值，國(guó)內(nèi)外學(xué)者已經(jīng)做了大量的工作。劉衛(wèi)等［26］選取了135塊低滲砂巖巖樣和8塊有代表性的碳酸鹽巖巖樣，通過(guò)離心法確定T2截止值的分布（低滲砂巖的離心力選取1.38 MPa，碳酸鹽巖的離心力選取為2.76 MPa）時(shí)，低滲砂巖求取的T2截止值平均為15.82 ms，碳酸鹽巖巖樣求取的平均值為86.78 ms，明顯高于砂巖。

肖前華等［27］通過(guò)選取中東某碳酸鹽巖油田具有代表性的70塊巖樣進(jìn)行可動(dòng)流體飽和度研究，用2.76 MPa的離心力標(biāo)定碳酸鹽巖的T2截止值后發(fā)現(xiàn)，中東地區(qū)碳酸鹽巖的T2截止值分布范圍很寬，在20～180 ms。這說(shuō)明，中東碳酸鹽巖儲(chǔ)層的非均質(zhì)性較強(qiáng)和縫洞較為發(fā)育，同時(shí)也說(shuō)明，部分石油公司推薦碳酸鹽巖儲(chǔ)層的可動(dòng)流體T2截止值采用90 ms并不是普遍適用的。趙輝等［28］選取磨溪?dú)馓锛味喍蔚?16塊碳酸鹽巖樣品，分別在飽水和離心2種狀態(tài)、TE為0.35 ms的情況下進(jìn)行核磁實(shí)驗(yàn)，并對(duì)他們進(jìn)行T2截止值標(biāo)定。結(jié)果表明：碳酸鹽巖的孔隙度對(duì)于其T2截止值的影響無(wú)明顯規(guī)律；碳酸鹽巖的T2截止值不是個(gè)定值，它與泥質(zhì)質(zhì)量分?jǐn)?shù)成反比。

Chang Dahai等［29］研究認(rèn)為，對(duì)于溶蝕孔洞發(fā)育的碳酸鹽巖，T2截止值選取92 ms較為適宜。此外，還有國(guó)外學(xué)者［30-33］認(rèn)為，砂巖的T2截止值采用33 ms較合理，而碳酸鹽巖的T2截止值采用100 ms較合理。

鄧克俊等［34］研究表明，砂巖的T2截止值變化范圍為20～50 ms，碳酸鹽巖為90～150 ms，并認(rèn)為碳酸鹽巖的ρ2很低是其T2截止值遠(yuǎn)大于砂巖的原因。Toumelin等［11］也認(rèn)為，當(dāng)溫度和“擴(kuò)散耦合”效應(yīng)影響很大時(shí)，需要對(duì)核磁共振測(cè)井過(guò)程中得出的T2截止值進(jìn)行校正。Vivek等在研究上述灰?guī)r中的“擴(kuò)散耦合”效應(yīng)問(wèn)題時(shí)也認(rèn)為，對(duì)于微孔和大孔之間存在“擴(kuò)散耦合”的地層，由于其T2時(shí)間和孔徑分布的相關(guān)性已經(jīng)不復(fù)存在，因而T2截止值的應(yīng)用可能會(huì)導(dǎo)致不正確的評(píng)價(jià)。Mai等［35］在研究在研究Alberta盆地中的碳酸鹽巖的NMR響應(yīng)特征后發(fā)現(xiàn)，碳酸鹽巖中的T2截止值變化很大，從小于80 ms到大于200 ms不等，并且認(rèn)為T2截止值隨巖石類型、孔隙系統(tǒng)等因素的變化而變化。

由此可見，相對(duì)于碎屑巖來(lái)說(shuō)，碳酸鹽巖的T2截止值普遍較高，且變化較大，做出的T2截止值的經(jīng)驗(yàn)值不具有普遍適用性。影響碳酸鹽巖T2截止值的因素較多，因而通過(guò)T2截止值法精確地獲得碳酸鹽巖的可動(dòng)流體飽和度還存在許多問(wèn)題。這一方面反映了碳酸鹽巖的非均質(zhì)性更強(qiáng)，溶蝕孔洞發(fā)育但分布不均；另一方面也可能與上述提到的碳酸鹽巖中存在的大孔和微孔中流體的“擴(kuò)散耦合”有一定關(guān)系。對(duì)此，值得去進(jìn)一步研究。

2.3 溫度對(duì)碳酸鹽巖NMR特征的影響

T S Ramakrishnan等［36］提出了一個(gè)用于研究碳酸鹽巖視表面弛豫強(qiáng)度的物理模型：

式中：ρa(bǔ)為視表面弛豫強(qiáng)度，μm/s；D0為孔隙流體的擴(kuò)散系數(shù)，m2/s；φμ為晶間孔孔隙度，%；Fμ為球狀微孔的地層因子；T2μ為微孔T2弛豫時(shí)間。

由于該式中第1項(xiàng)很小，可以忽略不計(jì)，因此，碳酸鹽巖的ρa(bǔ)約等于第2項(xiàng)。從該物理模型可以得到一種認(rèn)識(shí)——碳酸鹽巖的ρa(bǔ)受流體的擴(kuò)散系數(shù)的影響，而擴(kuò)散系數(shù)受溫度的影響很大，因而也可以進(jìn)一步推測(cè)，碳酸鹽巖的表面弛豫也受溫度的影響。當(dāng)碳酸鹽巖同時(shí)存在大孔和微孔時(shí)，該模型可以用于解釋碳酸鹽巖的T2分布受溫度影響的特征。

Godefroy等［37］分別對(duì)飽和水砂巖和碳酸鹽巖進(jìn)行NMR變溫實(shí)驗(yàn)后發(fā)現(xiàn)：在飽和水情況下，砂巖的弛豫時(shí)間隨著溫度的增長(zhǎng)而減小，但在飽和油的情況下則相反；而對(duì)于碳酸鹽巖，無(wú)論是飽和油還是飽和水，其弛豫時(shí)間都隨溫度的增長(zhǎng)而增加。Toumelin等［11］在研究溫度對(duì)存在“擴(kuò)散耦合”的碳酸鹽巖NMR特征的影響時(shí)發(fā)現(xiàn)，溫度對(duì)溶洞型大孔隙的耦合響應(yīng)沒(méi)有影響，但對(duì)小孔隙類型的耦合影響顯著，即隨著溫度增加，表面弛豫減小，即溫度對(duì)表面弛豫可能有相反的影響，但是，當(dāng)不存在“擴(kuò)散耦合”效應(yīng)時(shí)，則完全相反，所以，其對(duì)溫度的響應(yīng)還需作進(jìn)一步研究。Vivek等研究了溫度對(duì)于硅膠顆粒和Thamama組儲(chǔ)層中的飽水和束縛水狀態(tài)下的碳酸鹽巖樣品的NMR特征的影響，伴隨著溫度的升高，樣品的ρ值有所下降，即相同條件下弛豫時(shí)間變長(zhǎng)、譜峰發(fā)生右移（見圖7），同時(shí)樣品中大孔和微孔中流體的“擴(kuò)散耦合”程度會(huì)大幅增加，甚至?xí)?dǎo)致大孔和微孔的T2譜合二為一（見圖8）。

謝然紅等［38-39］通過(guò)對(duì)不同的人造和天然巖心（含碳酸鹽巖儲(chǔ)層巖心）進(jìn)行NMR變溫測(cè)試后發(fā)現(xiàn)，飽水的人造碳酸鈣和碳酸鹽巖的橫向ρ和1/T2隨溫度升高而減小（或T2時(shí)間隨溫度升高而增大），與貝瑞砂巖儲(chǔ)層及人造碳化硅的結(jié)果剛好相反。這和碳酸鹽巖中與分子擴(kuò)散相關(guān)的活化能Em大于與巖石顆粒表面相互作用的活化能 （即碳酸鹽巖的橫向ρ比碎屑巖小一個(gè)數(shù)量級(jí)）有關(guān)。而對(duì)于飽和油的樣品，無(wú)論是貝瑞砂巖或碳酸鹽巖儲(chǔ)層，其T2時(shí)間都隨溫度的升高而增大。

由此可見，溫度對(duì)于碳酸鹽巖的NMR響應(yīng)特征的影響，因碳酸鹽巖的礦物成分、孔隙類型以及流體性質(zhì)的變化而變化，情況較為復(fù)雜，不同學(xué)者得出的結(jié)論也大相徑庭。目前，國(guó)內(nèi)在溫度對(duì)碳酸鹽巖的NMR特征的影響方面的研究還很不夠，有待進(jìn)一步深化。

圖7 束縛水狀態(tài)下中—粗粒硅膠顆粒的T2分布

圖8 Thamama組碳酸鹽巖不同條件下的T2分布

3 結(jié)論

1）巖石中大孔和微孔間的流體存在“擴(kuò)散耦合”，本質(zhì)上是由于巖石中順磁-鐵磁性礦物質(zhì)量分?jǐn)?shù)較少，使得巖石的表面弛豫速率降低所致，與巖性并沒(méi)有必然聯(lián)系。

2）國(guó)內(nèi)外學(xué)者研究表明，和砂巖相比，碳酸鹽巖的T2截止值普遍較高，且變化較大，所選取的經(jīng)驗(yàn)值不具有普遍適用性，其影響因素較多。因而用T2截止值來(lái)計(jì)算碳酸鹽巖的可動(dòng)流體飽和度還存在很大的問(wèn)題。

3）從不同學(xué)者的研究結(jié)果來(lái)看，溫度對(duì)碳酸鹽巖的NMR響應(yīng)特征的影響，因碳酸鹽巖的礦物成分、孔隙類型以及流體性質(zhì)而異，且仍存在較大爭(zhēng)議，目前國(guó)內(nèi)在這方面的研究較為薄弱，有待進(jìn)一步深化。

［1］李國(guó)玉，陳啟林，白云來(lái)，等.再論海相沉積是中國(guó)石油工業(yè)未來(lái)的希望［J］.巖性油氣藏，2014，26（6）：1-7.

［2］符幸子，徐守余，王淑萍，等.塔河油田六區(qū)碳酸鹽巖溶洞型儲(chǔ)集體三維地質(zhì)建模［J］.斷塊油氣田，2015，22（4）：419-424.

［3］韓文學(xué)，高長(zhǎng)海，韓霞.核磁共振及微、納米CT技術(shù)在致密儲(chǔ)層研究中的應(yīng)用：以鄂爾多斯盆地長(zhǎng)7段為例［J］.斷塊油氣田，2015，22（1）：62-66.

［4］張松揚(yáng)，范宜仁.核磁共振測(cè)井技術(shù)評(píng)述［J］.勘探地球物理進(jìn)展，2002，25（4）：21-26.

［5］王為民，趙剛，谷長(zhǎng)春，等.核磁共振巖屑分析技術(shù)的實(shí)驗(yàn)及應(yīng)用研究［J］.石油勘探與開發(fā)，2005，32（1）：56-59.

［6］林振洲，潘和平.核磁共振測(cè)井進(jìn)展［J］.工程地球物理學(xué)報(bào)，2006，3（4）：295-303.

［7］趙平，周利軍，紀(jì)常杰.核磁共振測(cè)井技術(shù)的進(jìn)展［J］.國(guó)外油田工程，2007，23（10）：32-35，38.

［8］ROBERTS S P，MACDONALD P J，RICHARD T.A bulk and spatially resolved NMR relaxa-tion study of sandstone rock plugs［J］.Journal of Magnetic Resonance：Series A，1995，116（2）：189-195.

［9］TIMUR.Nuclear magnetic resonance study of carbonate rocks［J］.The Log Analysts，1972，13（5）：3-11.

［10］RAMAKRISHNAN T S，F(xiàn)ORDHAM E J，VENKATARAMANAN L，et al.New interpretation methodology based on forward models for magnetic resonance in carbonates［C］//The Society of Professional Well Log Analysts，Inc.Transactions of the SPWLA 40th Annual Logging Symposium.SPWLA 40th Annual Logging Symposium，30 May—3June，1999.Oslo，Norway.Houston，TexasSPWLA，1999：1-14.

［11］TOUMELIN E，TORES V C，CHEN S，et al.Analysis of NMR diffusion coupling effects in two-phase carbonate rocks：Comparison of measurements with Monte Carlo simulations［J］.Petrophysics，2003，44（2）：91-107.

［12］VIVEK A，GEORGE J H.Diffusional coupling between micro and macroporosityfor NMR relaxation in sandstones and grainstones［J］. Petrophysics，2007，48（4）：289-307.

［13］肖立志，謝然紅，廖廣志.中國(guó)復(fù)雜油氣藏核磁共振測(cè)井理論與方法［M］.北京：石油工業(yè)出版社，2012：267-269.

［14］LU C，ZOYA H.Diffusional couplingbetween micro-fractures and pore structure and itsimpact on NMR measurements in multiple-porosity systems［J］.Geophysics，2015，80（1）：31-42.

［15］CHI L，HEIDARI Z.Assessment of micro-fracture density using combined interpretation of NMR relaxometry and electromagnetic logs［R］.SPE 1922804，2014.

［16］李傳亮.氣水可以倒置嗎?［J］.巖性油氣藏，2010，22（2）：128-132.

［17］朱蘇陽(yáng)，李傳亮.對(duì)頂水油藏概念的質(zhì)疑［J］.巖性油氣藏，2013，25（1）：21-23.

［18］DASTIDAR R，CARL H，SONDERGELD C.NMR Desaturation and surface relaxivity measurements on clastics rocks［R］.SPE 99629，2006.

［19］KRIXTINA K.The effect of inhomogeneous surface relaxivity on nuclear magnetic resonance relaxation rates［C］//Levin S.SEG Technical Program Expanded Abstracts.SEG Annual Meeting，17-22 October，2010，Denver，Colorado.SocietyofExplorationGeophysicists，2010：3946-3950.

［20］黃仲良.石油重·磁·電法勘探［M］.東營(yíng)：中國(guó)石油大學(xué)出版社，1998：214-215.

［21］張勝業(yè)，潘玉玲.應(yīng)用地球物理學(xué)原理［M］.武漢：中國(guó)地質(zhì)大學(xué)出版社，2004：6-13.

［22］屈樂(lè)，孫衛(wèi)，章海寧，等.火成巖儲(chǔ)層核磁共振巖石物理影響特征［J］.蘭州大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，50（1）：26-30.

［23］司馬立強(qiáng)，趙輝，戴詩(shī)華.核磁共振測(cè)井在火成巖地層應(yīng)用的適應(yīng)性分析［J］.地球物理學(xué)進(jìn)展，2012，27（1）：145-152.

［24］趙澄林，朱筱敏.沉積巖石學(xué)［M］.3版.北京：石油工業(yè)出版社，2001：99-116.

［25］姜在興.沉積學(xué)［M］.北京：石油工業(yè)出版社，2003：403-404.

［26］劉衛(wèi)，邢立，孫佃慶，等.核磁共振錄井［M］.北京：石油工業(yè)出版社，2011：39-42.

［27］肖前華，張亞蒲，楊正明，等.中東碳酸鹽巖核磁共振實(shí)驗(yàn)研［J］.科學(xué)技術(shù)與工程，2013，13（22）：6415-6420.

［28］趙輝，司馬立強(qiáng)，顏其彬，等.碳酸鹽巖巖樣核磁共振實(shí)驗(yàn)特征分析［J］.測(cè)井技技術(shù)，2011，35（2）：117-121.

［29］CHANG D H，HAROLD J V，CHRIS M，et al.Effective porosity，producible fluid and permeability in carbonates from NMR logging［J］. The Log Analysts，1997，38（2）：60-72.

［30］CHRISTIAN S，DAN R，HAROLD V，et al.Core analysis by low-field NMR.［J］.The Log Analysts，1997，38（3/4）：43-56.

［31］KLEINBERG R L，BOYD A.Tapered cutoffs for magnetic resonance bound water volume［R］.SPE 38737，1997.

［32］DUNN K J，LATORRACA G A，WARNER J L，et al.On the calculation and interpretation of NMR relaxation time distributions［R］.SPE 28367，1994.

［33］CANNON D E，MINH C C，KLEINBERG R L.Quantitative NMR interpretation［R］.SPE 49010，1998.

［34］鄧克俊，謝然紅.核磁共振測(cè)井理論及應(yīng)用［M］.青島：中國(guó)石油大學(xué)出版社.2010：60-62.

［35］MAI A，KANTZAS A.On the characterization of carbonate reservoirs using low field NMR tools［R］.SPE 75687，2002.

［36］RAMAKRISHNAN T S，SCHWARTZ L M，F(xiàn)ORDHAM E J，et al. Forward models for nuclear magnetic resonance in carbonate rocks［J］. The Log Analysts，1999，40（6）：260-270.

［37］GODEFROY S，F(xiàn)LEURY M，DEFLANDRE F，et al.Temperature effect on NMR surface relaxation［R］.SPE 71700，2001.

［38］謝然紅，肖立志.儲(chǔ)層流體及其在巖石孔隙中的核磁共振弛豫溫度特性［J］.地質(zhì)學(xué)報(bào)，2007，81（2）：284-288.

［39］謝然紅，肖立志，傅少慶.飽和水巖石核磁共振表面弛豫溫度特性［J］.中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2008，32（2）：44-46，52.

（編輯 高學(xué)民）

Research status of nuclear magnetic resonance experiment in carbonate rock

WANG Yijun1,CUI Gang2,TANG Hongming1,SHEN Jianjun1,GAO Zihan3,ZHOU Shilin1
(1.School of Geoscience and Technology,Southwest Petroleum University,Chengdu 610500,China;2.Tianjin Branch of CNOOC Ltd.,Tianjin 300452,China;3.Bureau of Geophysical Prospecting Inc.,CNPC,Zhuozhou 072751,China)

Along with the large-scale development of marine carbonate oil and gas fields,the evaluation of carbonate reservoir is becoming increasingly important.Low field nuclear magnetic resonance(NMR)core analysis,a technology introduced into China just 20 years ago,has been used ever since for rapidly measuring petrophysical parameters.It has been widely applied to reservoir fluid identification and reservoir evaluation.However,the application of NMR in carbonate reservoir still has some problems at present. This paper reviews academic research results in China and abroad concerning the following three issues:big and small pores′"diffusional coupling"in carbonate rocks,the accurate calculation of T2cutoff value of carbonate rocks,and the influence of temperature on the carbonate NMR response characteristics.On this basis,it presents three findings:first,"diffusional coupling" exists in all kinds of rocks though it is more apparent in carbonate rocks than in other kinds of rocks;second,as the T2cutoff value of carbonate rocks is generally higher than that of sandstone and also fluctuates within a broader range,it follows that carbonate rocks are highly heterogeneous,and that the T2cutoff value may be associated with the"diffusional coupling"in the carbonate rocks;third,the influence of temperature on the carbonate NMR response characteristics varies with different mineral compositions,pore types and fluid properties.

nuclear magnetic resonance(NMR);diffusional coupling;T2cutoff value;temperature;surface relaxivity

國(guó)家自然科學(xué)基金項(xiàng)目“致密氣藏儲(chǔ)層干化、提高氣體滲流能力的基礎(chǔ)研究”（51534006）

TE344

A

10.6056/dkyqt201606028

2016-04-28；改回日期：2016-08-07。

王翼君，男，1989年生，在讀碩士研究生，從事儲(chǔ)層地質(zhì)、儲(chǔ)層保護(hù)相關(guān)研究。E-mail：898789370@qq.com。

王翼君，崔剛，唐洪明，等.碳酸鹽巖核磁共振實(shí)驗(yàn)研究現(xiàn)狀［J］.斷塊油氣田，2016，23（6）：818-824.

WANG Yijun，CUI Gang，TANG Hongming，et al.Research status of nuclear magnetic resonance experiment in carbonate rock［J］.Fault-Block Oil&Gas Field，2016，23（6）：818-824.