吳 豐，司馬立強，令狐松，姚軍朋，田蔥蔥

(1.西南石油大學，四川 成都 610500;2.中油測井有限公司，陜西 西安 710021)

1 地質概況

昆北斷階帶位于柴達木盆地西部地區(qū)昆北逆沖帶南側祁漫塔格山前，包括鐵木里克凸起、切克里克凹陷、東柴山以及黃石凸起南側。由于昆侖山的抬升，導致構造活動比較強烈［1］。該區(qū)從早第三紀開始接受沉積，勘探開發(fā)目的層主要為路樂河組(E1+2l)及下干柴溝組下段(E3g1)。E1+2l—E3g1時期主要沉積有辮狀河三角洲平原亞相、辮狀河三角洲前緣亞相、濱淺湖亞相?；鶐r以花崗巖為主，其次為變質巖;E1+2l、E3g1地層主要發(fā)育棕紅色泥巖、砂質泥巖，棕紅色或灰黃色粉砂巖、泥質粉砂巖、礫巖等，巖石中含大量巖屑成分。巖石壓實成巖作用中等、溶蝕和膠結成巖作用較弱，孔隙類型以原生粒間孔為主。

2 核磁共振特征

對昆北地區(qū)E1+2l—E3g1儲層進行大量核磁共振測井［2-3］。通過與巖心氦孔隙度對比，發(fā)現(xiàn)核磁共振測井孔隙度普遍偏低。為了進一步證實該現(xiàn)象，挑選了10個樣品進行巖心核磁共振測量。實驗采用紐邁科技的全直徑巖心核磁共振成像分析系統(tǒng)(AniMR-150型)，該儀器可測量全直徑巖心、標準巖樣、巖屑，除了根據(jù)橫向弛豫時間(T2)分布得到孔隙度、滲透率、流體飽和度、孔徑分布等重要參數(shù)以外，還可進行三維成像，系統(tǒng)升級后回波間隔TE最低可調至0.1 ms。

(1)核磁共振孔隙度特征。10個樣品的巖心核磁孔隙度平均值為8.3%，與之對應的核磁共振測井總孔隙度平均值為7.9%，而巖心氦孔隙度平均值為12.2%，核磁共振孔隙度與巖心氦孔隙度的相對誤差達到32.0%。

(2)核磁共振T2分布特征。由圖1可以看出，巖心核磁共振的T2曲線分布大部分呈雙峰形態(tài)，少部分呈單峰形態(tài)，T2偏小。隨著回波間隔增大，T2分布逐漸丟失部分信號，尤其是短弛豫時間部分的信號丟失較多(圖1a)。對飽和水巖樣進行離心前后測量，1.138 MPa離心力得到的T2截止值為1.5～3.0 ms(圖1b)，該值與砂巖的理論T2截止值 33 ms［4］相比低了很多。

一般認為核磁共振測井直接反映巖石孔隙中的流體，測量結果幾乎不受巖石骨架礦物的影響，也有部分研究成果指出火成巖和綠泥石含量高的巖石受骨架影響較大［5-6］。

圖1 不同回波間隔與不同離心力的巖心核磁共振T2分布

3 孔隙度偏低原因詳析

3.1 綠泥石含量

對昆北地區(qū)的全巖X衍射及黏土X衍射分析表明:E1+2l—E3g1地層綠泥石含量較低(砂礫巖黏土含量分布范圍為5.0% ～30.0%，綠泥石占黏土總量的12.0%左右，即綠泥石僅占巖石骨架的0.6%～3.6%)。而一般綠泥石含量需要達到20.0%～30.0%才能對T2分布造成較明顯的影響［7］，因此可知，昆北地區(qū)核磁共振受影響的主要原因并非綠泥石的存在。

3.2 順磁性礦物

昆北地區(qū)E1+2l—E3g1時期巖心薄片分析表明，巖石組分中除長石顆粒、石英顆粒、泥質外，還有大量巖屑顆粒。其巖屑顆粒主要來源于花崗巖和變質巖(千枚巖、板巖、石英片巖等)，巖屑含量分布范圍為10.0% ～29.0%，平均含量為19.8%。錄井重礦物分析表明，E1+2l—E3g1時期儲層重礦物含量總體較高，平均可達500顆/m，而柴達木盆地其他地區(qū)平均為100顆/m左右。同時，E1+2l—E3g1儲層地層磁鐵礦、黃鐵礦、赤鐵礦的含量很高，某些層段超過了1000顆/m，而柴達木盆地其他地區(qū)磁鐵礦、黃鐵礦、赤鐵礦的平均含量非常低，如NBX地區(qū)約10顆/m左右(圖2)。

圖2 昆北地區(qū)磁鐵礦、赤鐵礦、黃鐵礦分布

磁鐵礦、黃鐵礦、赤鐵礦具有高磁化率，尤其是磁鐵礦最為明顯。磁鐵礦、黃鐵礦、赤鐵礦等含有順磁性礦物，一方面使得橫向表面弛豫強度增大，另一方面使得巖石產(chǎn)生內部磁場梯度［7］。這2種情況會減小地層流體孔隙的橫向弛豫時間T2，一些 較小T2的組分衰減很快，儀器測量不到這部分信號，導致核磁共振測量孔隙度降低。

4 校正核磁共振孔隙度

4.1 孔徑分布轉換

由于核磁共振測量受巖石內部梯度磁場的影響機理非常復雜，目前由順磁物質引起的核磁共振孔隙度減小還無法校正［8］，因此從其他方向入手研究。核磁共振T2分布曲線與毛管壓力曲線具有一定相似性，都可用來表征巖石孔隙結構。毛管壓力與毛管孔徑之間的關系為［9-11］:

式中:pc為毛管壓力，MPa;σ為流體界面張力，N/m;θ為潤濕接觸角，(°);Rc為毛管半徑，μm。

核磁共振觀測到的橫向弛豫時間T2可表示為［12-15］:

式中:T2B為流體體積弛豫時間，ms;T2D為流體擴散弛豫時間，ms;T2S為流體表面弛豫時間，ms;D為擴散系數(shù)，μs2/ms;γ為旋磁比，rad/(s·T);G為磁場梯度，Gs/cm;TE為回波間隔，ms;S為孔隙表面積，cm2;V為孔隙體積，cm3;ρ2為巖石橫向表面弛豫強度，μm/ms。

一般情況下，T2B的數(shù)值通常大于3 000 ms，要比T2大得多，式(2)右邊第1項可以忽略。當磁場很均勻(對應G很小)，且TE足夠小時，式(2)右邊第2項也可以忽略。此時，核磁共振T2與孔徑Rc的關系可表示為:

式中:Fs為孔隙形狀因子。

將T2分布轉換為毛管孔徑Rc的基礎，是建立在式(2)右邊第1項和第2項均可忽略的情況下。然而，昆北地區(qū)E1+2l—E3g1儲層含有大量順磁性礦物磁鐵礦、赤鐵礦、黃鐵礦，導致巖石內部磁場梯度G變大，式(2)右邊第2項不能忽略，因此直接由核磁共振轉換得到的視毛管半徑Rc與毛管壓力曲線轉換的毛管半徑Rc之間存在一定差異，下面對二者之間差異進行分析。

4.2 核磁共振與壓汞孔徑分布對比

將巖心核磁共振T2分布的孔隙度分量數(shù)據(jù)從高端向低端累計，得到一條累計曲線;將其與壓汞飽和度累計曲線進行擬合，將T2分布轉換為視毛管壓力曲線，再進一步轉換為視孔徑分布。將核磁共振視孔徑分布與壓汞孔徑分布放在同一張圖中，孔喉半徑與相對應的T2弛豫時間作為橫坐標軸，縱坐標軸為孔隙度分量(壓汞孔徑分布采用常規(guī)物性分析孔隙度進行刻度)，包絡線所包圍的面積可代表各自孔隙度的大小，二者的差值部分即為核磁共振丟失的孔隙信號。

以B31號樣品為例(圖3)，巖心氦孔隙度為10.14%，巖心核磁孔隙度為6.69%，二者孔隙度相差3.45個百分點。在孔徑分布對比圖上，壓汞孔徑分布與核磁共振視孔徑分布的差值主要集中在2個峰值附近，尤其是左峰的差值比較大。左峰孔徑較小，約為0.007～0.037 μm，對應的核磁共振T2弛豫時間為1.3～8.7 ms;右峰孔徑相對較大，在0.76 μm附近，對應的核磁共振T2弛豫時間為81 ms左右。

4.3 核磁共振孔隙度校正

昆北地區(qū)E1+2l—E3g1儲層核磁共振測井以哈里伯頓MRIL-P型儀器“雙TW/單TE”模式為主，核磁共振總孔隙度可表示為:

式中:PPR0.5、PPR1、PPR2、PPR4依次為由回波串反演得到的 T2(分別等于 0.5、1.0、2.0、4.0 ms 時)的孔隙度，%;PA8、PA16、PA32、PA64、PA128、PA256、PA512、PA1024、PA2108依次為回波串反演得到的T2(分別等于8、16、32、64、128、256、512、1 024、2 048 ms時)的孔隙度，%。

圖3 巖心核磁共振視孔徑分布與壓汞孔徑分布對比(樣品B31)

圖4 核磁共振測井孔隙度校正實例(B井)

根據(jù)文中方法校正后核磁共振總孔隙度可表 示為:

式中:n1～n13為孔隙度分量加權系數(shù)。

通過核磁共振與壓汞孔徑分布對比，確定出加權系數(shù)分別為:n2、n3、n4取值范圍為 1.5～2.0，n8取2.5，其他系數(shù)取1.0。

5 實例分析

圖4為柴達木盆地昆北地區(qū)B井的MRIL-P型核磁共振測井圖，在1 814～1 831 m井段，原始核磁總孔隙度和巖心核磁孔隙度明顯低于巖心分析氦孔隙度，利用式(5)對孔隙度分量進行補償后，重新組構的核磁總孔隙度曲線與巖心分析氦孔隙度更吻合。

6 結論

(1)柴達木盆地昆北地區(qū)核磁共振測井孔隙度較巖心孔隙度偏低，T2分布弛豫時間偏小，多表現(xiàn)為雙峰狀態(tài)，離心法測得的T2截止值為1.5～3.0 ms，較砂巖理論T2截止值33 ms偏低。

(2)柴達木盆地昆北地區(qū)巖屑組分主要來源于花崗巖和變質巖，重礦物含量較高，尤其是磁鐵礦、黃鐵礦、赤鐵礦的含量高于柴達木盆地其他地區(qū)，導致昆北地區(qū)核磁共振孔隙度降低。

(3)壓汞曲線與核磁共振T2分布均可轉化為巖石孔徑分布。當巖石含順磁性礦物時，核磁共振轉化的視孔徑分布發(fā)生畸變，將其與壓汞的孔徑分布進行對比，找出二者之間的差異，并對差異部分進行補償校正，即可重構核磁共振孔隙度曲線。

［1］石亞軍，曹正林，張小軍，等.大型高原內陸咸化湖盆油氣特殊成藏條件分析及勘探意義——以柴達木盆地柴西地區(qū)為例［J］.石油與天然氣地質，2011，32(4):577-583.

［2］羅婷婷，周立發(fā)，劉志武.核磁共振技術在渭河盆地奧陶系儲層評價中的應用［J］.特種油氣藏，2011，18(2):37-40.

［3］李太和，郭何坤，李海波，等.應用核磁共振技術研究頁巖氣儲層可動流體［J］.特種油氣藏，2012，19(1):107-109，123.

［4］蘇俊磊，王艷，孫建孟.應用可變T2截止值確定束縛水飽和度［J］.吉林大學學報:地球科學版，2010，40(6):1491-1520.

［5］司馬立強，趙輝，戴詩華.核磁共振測井在火成巖地層應用的適應性分析［J］.地球物理學進展，2012，27(1):145-152.

［6］廖廣志，肖立志，謝然紅，等.內部磁場梯度對火山巖核磁共振特性的影響及其探測方法［J］.中國石油大學學報:自然科學版，2009，33(5):55-60.

［7］肖立志，謝然紅，廖廣志.中國復雜油氣藏核磁共振測井理論與方法［M］.北京:科學出版社，2012:44-45.

［8］謝然紅，肖立志.核磁共振測井探測巖石內部磁場梯度的方法［J］.地球物理學報，2009，52(5):1341-1347.

［9］劉堂宴，馬在田，傅容珊，等.核磁共振譜的巖石孔喉結構分析［J］.石油地球物理勘探，2003，38(3):328-333.

［10］劉紅現(xiàn)，許長福，胡志明.用核磁共振技術研究剩余油微觀分布［J］. 特種油氣藏，2011，18(1):96-97，125-126.

［11］尹邦堂，李相方，李騫，等.利用核磁共振技術測量含油巖樣滲透率［J］.大慶石油地質與開發(fā)，2011，30(6):117-121.

［12］運華云，譚茂金.核磁共振測井雙等待時間觀測方式及其分析方法［J］.油氣地質與采收率，2006，13(4):96-98.

［13］田素月，孫靈芬.核磁共振測井在東濮凹陷深層氣評價中的應用［J］.油氣地質與采收率，2003，10(3):34-35.

［14］王志戰(zhàn)，許小瓊，周寶潔.孔隙流體核磁共振弛豫特征及油水層識別方法［J］.油氣地質與采收率，2011，18(2):41-44.

［15］趙文杰.利用核磁共振測井資料計算平均孔喉半徑［J］. 油氣地質與采收率，2009，16(2):43-45.