GeoReveal M RIL-P核磁解譜技術在低孔滲儲層評價中的應用

羅利 王勇軍 楊嘉 謝剛 劉文龍

川慶鉆探工程有限公司測井公司

GeoReveal M RIL-P核磁解譜技術在低孔滲儲層評價中的應用

羅利 王勇軍 楊嘉 謝剛 劉文龍

川慶鉆探工程有限公司測井公司

從MRIL-P核磁共振測井儀器硬件優勢出發,深入闡述了MRIL-P核磁 T2解譜3個最核心的問題——解譜目標函數的確立、解向量的非負處理以及正則化方法的選擇。在對MRIL-P型核磁解譜技術方法研究的基礎上,展示了不同信噪比、不同巖性解譜的效果。JLS構造下侏羅統珍珠沖組礫巖儲層具有巖性復雜、儲集空間類型多樣、低孔滲等特點,準確求取地層孔隙度是測井解釋評價難點之一,結合MRIL-P型核磁共振解譜方法研究成果,提出采用求取最優正則化因子的方法來進行處理。處理結果證實MRIL-P型核磁共振測井適用于低孔滲儲層孔隙度測井需求。

GeoReveal解釋平臺 核磁測井 解譜 約束法 最優化 準則 低滲透儲集層

核磁共振測井儀器MRIL-P型相對于MRIL-C型而言,增加了更多組工作頻率,去掉了對地層FID信號進行處理電路,對核磁共振信號的分析轉而由測后數據處理完成。改進使得MRIL-P型核磁共振測井具有以下優勢:①大幅提高測井時效,M RIL-P型有9組不同的工作頻率,然后又被分為5組不同的工作頻帶,因此在一個單周期內,可以同時采集到5組不同測井參數的回波信號,而M RIL-C型核磁只能采集兩組;②能提高信噪比,在相同井眼條件下采集等量的測井信息,MRIL-P型核磁很容易將雙頻采集量增大4倍,FID信號幅度協調地增加,而噪聲成分只按均方根增加,因而,在理想情況下噪聲成分會被抑制2倍[1],在后續數據處理過程中又增加信號疊加、相位校正、信噪分離等技術進一步抑制噪聲;③針對不同的測井目的提供了近百種測井采集模式,并在測井參數設置上提供了更長的等待時間和更短的回波間隔(TWmax= 15 s,TEmin=0.6 m s),應用范圍更廣。但是,這使得后續數據處理復雜程度大為增加,必須深入研究MRIL-P型核磁共振測井數據的解譜方法。

1 解譜方法

因為核磁共振測量信號弛豫過程均按照指數規律衰減,采用多指數擬合,總的弛豫為多項弛豫的疊加[2-3]:

式中i為回波信號的序數,無量綱;j為橫向弛豫分量的序數,無量綱;gi為ti時刻觀測到的回波幅度,m V; T2j為第 j個弛豫分量的橫向弛豫時間,m s;fj為第 j個橫向弛豫分量零時刻信號的大小,%;m為橫向弛豫分量的個數,無量綱;εi為第i回波的噪聲信號,m V。且

式中 TE為回波間隔,m s。

核磁解譜的目的就是求解所有橫向弛豫分量零時刻信號的大小 f={fi,…,fj,…,fm},將式(2)代入式(1)中,假設回波數為n,信噪已完全分離,可以不考慮噪聲項εi,可寫成為:

根據最小二乘原理,解譜即求[3]:

因為回波數n遠大于橫向弛豫分量的個數m,因此,求解該方程是一個“病態”問題,即測量數據 gi即使有一個細小的波動,都會導致 fj解的結果有很大的差異。為了得到穩定的解,以正則化的思想為基礎,在式(5)中增加約束項,使用改進奇異值分解算法求解。增加約束項的方法比較多,以約束項采用解的曲率為例,式(5)改為:

式中<為目標函數;α為正則化參數,無量綱。

因為該方法約束項采用的是解的曲率,故又稱為曲率約束法。為了更方便地求解 f,需要對系數矩陣A進行奇異值分解[3]:

根據矩陣理論,當減小A的條件數,解的相對誤差變化范圍會減小,矩陣A的條件數與奇異值有關, ω1為A的最大奇異值,ωk為A的最小奇異值,則A的條件數為:

現在考慮去掉某一個奇異值ωp后,對解產生的影響,可以推導解的平方誤差增加了D2p,Dp被定義為:

要想得到穩定的解,A的條件數要小,把小的非零奇異值截去,則會增加A的條件數,使解更穩定。要去掉小的非零奇異值,同時也增加了平方誤差。均衡考慮后,在計算過程中做如下處理,將式(8)變為:

根據核磁解譜的物理意義可知,解向量{fj}必須是非負的,通常處理方法是:刪除矩陣A中與fj最小負值的分量所對應的列,同時將 fj置0,然后重新進行計算,得到一個新的解,再判斷正負。如果fj還有負值,則刪除fj最小的負值分量對應的列,同時將 fj置0,重復進行下去,直到所有的解全為非負時結束運算[4]。

另外,在實施正則化方法過程中,確定合適的正則化參數α是核心問題之一[4],其基本要求是對于待求解施加定性的(如光滑性要求)或定量的信息,以便使正則參數與原始數據資料的誤差水平相匹配,通常有先驗的和后驗的兩種策略。

當誤差水平已知時,廣為采用的后驗策略是所謂的偏差原理和廣義偏差原理[5],前者由 Mo rozov于1966年提出,而后者是由 Tikhonov的學生 Goncharsky等人對于Morozov偏差原理進行改造和推廣而得到的[5]。用上述方法確定正則化參數時,都需要預先對于原始數據的誤差水平下做出估計;但在有些情況下是難以做到的,或者即使勉強做出了估計也是不可靠的。許多先驗性的策略便于理論分析,但一般來說難于驗證,故采用后驗的方法居多。現在有了在誤差水平未知情況下的各種后驗準則,如 Tikhonov的擬最優準則[5],Hansen的L曲線準則[6]和 Engl的誤差極小化準則[7]等。

2 解譜效果對比

根據上述計算方法,對川渝地區數口井進行了解譜處理,再經時深轉換和插值處理就能得到 T2譜,同時計算出總孔隙度、有效孔隙度、滲透率等巖石物理參數。這里,選其中3口具有代表性的井,將處理得到的T2譜和巖石物理參數結果與國外DPP軟件處理結果進行對比(圖1)。

圖1 通過方法研究處理的結果與DPP工作站處理的結果對比圖

圖1中本文方法處理的結果(標有STU)與DPP工作站處理的結果(標有DPP)對比可見,兩種方法得到的 T2解譜分布范圍、峰值和波形形態一致性好。部分巖石物理參數計算結果對比見表1。

表1 各井STU、DPP計算結果對比表

結合圖1和從表1可知,對于高信噪比砂泥巖地層(E2井)、低信噪比碳酸巖地層(L 18 I井)、低信噪比砂泥巖地層(PL2井)的MRIL-P型核磁測井數據,筆者所述方法與DPP工作站解譜結果一致性好,有效孔隙度絕對誤差均在1.0個孔隙度單位以內,相對誤差不超過5%,滲透率絕對誤差在0.1 mD以內。

3 在低孔滲儲層評價中的應用

JLS構造下侏羅統珍珠沖組礫巖儲層是較為典型的低孔滲儲層,其測井解釋有以下難點。

1)巖性復雜。取心資料表明,珍珠沖組礫巖,以石英礫為主,夾燧石礫,磨圓好;礫巖間隙夾砂,見含砂礫巖,含礫砂巖;部分井礫巖還含少許灰質及黃鐵礦。

2)儲集空間類型多樣。珍珠沖組礫巖儲集空間類型涵蓋了網狀、垂直、高角度裂縫,溶孔、溶洞,粒間孔、晶間孔等多種類型。

3)低孔低滲。通過統計分析巖心孔隙度、滲透率數據表明,珍珠沖組礫巖儲層孔隙度均小于5%,一般在2%～3%;巖心滲透率一般低于0.01 mD。

因為巖性復雜,常規三孔隙度測井巖石骨架值不能確定,準確求取地層孔隙度難度大。為此,針對該礫巖段的MRIL-P型核磁測井數據,將式(11)寫成:

迭代計算方法求最優正則化因子:

圖2 501L井巖心孔隙度與常規、核磁孔隙度對比圖

圖2中井段3 272～3 274 m所示,常規處理孔隙度與巖心孔隙度對比,因為GR測井值相對上下圍巖較高,所以通過數字處理得到的泥質含量較高,進而計算得到的孔隙度較低(幾乎為0),但實際巖心孔隙度約3%。這一差異產生的原因可能是因為高 GR測井值并非完全來自于泥質的貢獻,而是巖體自身放射性較高所致。在求最優正則化因子α前,先對α取值范圍進行搜索,α分別取值SM04、SM1、SM5時與巖心孔隙度的重疊情況如第6道所示,據此確定正則化因子α取值范圍為SM04≤α≤SM5,在該取值范圍內,按式(12)進行迭代求α的最優解。當α=SM08時,式(12)達到最小,所求核磁有效孔隙度經含烴校正后與巖心孔隙度達到最佳匹配。

4 結論

1)M RIL-P型核磁 T2譜解譜3個核心問題分別為:①解譜目標函數的確立;②解向量的非負處理;③正則化方法的選擇。筆者介紹方法能有效解決這3方面的問題,與國外DPP軟件處理得到的 T2譜和巖石物理參數結果一致性好。

2)以解譜方法研究成果為基礎,發展的求最優正則化因子α方法切實可行。

3)MRIL-P型核磁能夠適應低孔滲儲層孔隙度測井需求,能真實地反映地層的孔隙度。

[1]程雪蓮,李浩,張松揚,等.新一代多頻帶核磁共振測井儀[J].石油物探譯叢,2001,6(1):87-92.

[2]王為民,李培,葉朝輝.核磁共振弛豫信號的多指數反演[J].中國科學:A輯,2001,31(8):730-736.

[3]王才志,尚衛忠.應用奇異值分解算法的核磁共振測井解譜方法[J].石油地球物理勘探,2003,38(1):91-94.

[4]DUNN K J,LA TORRACA G A,WARNER J L,et al.On the calculation and interp retation of NMR relaxation time distributions[C]∥paper 28367-MS p resented at the SPE Annual Technical Conference and Exhibition,25-28 September 1994,New Orleans,Louisiana,USA.New York: SPE,1994.

[5]袁振超,沈云中,周澤波.病態總體最小二乘模型的正則化算法[J].大地測量與地球力學,2009,29(2):131-134.

[6]徐天河,楊元喜.均方誤差意義下正則化解優于最小二乘解的條件[J].武漢大學學報:信息科學版,2004,29(3): 223-226.

[7]李功勝,馬逸塵.應用正則化子建立求解不適定問題的正則化方法的探討[J].數學進展,2000,29(6):531-541.

Application of the GeoReveal M RIL-P T2 spectrum inversion method to the evaluation of low porosity&permeability reservoirs

Luo Li,Wang Yongjun,Yang Jia,Xie Gang,Liu Wenlong
(Logging Com pany of Chuanqing D rilling Engineering Co.,L td.,CN PC,Chongqing 400021,China)

NATUR.GAS IND.VOLUM E 31,ISSUE 2,pp.59-62,2/25/2011.(ISSN 1000-0976;In Chinese)

Beginning w ith an introduction of the hardware advantages of Magnetic Resonance Imaging Logging Prime(MRIL-P) Tool,this paper discusses three co re issues of M RIL-P T2 spectrum inversion,namely the establishment of objective function,nonnegative constraints of solution-vecto r and the selection of app rop riate regularization method.Based on a study of MRIL-P T2 spectrum inversion technology,we discuss the effectsof spectrum inversion of different lithologies at different SNRs.The conglomerate reservoirs in the Lower Jurassic Zhenzhuchong Formation of JLS structure are characterized by comp lex lithology,multiple typesof reservoir space,low porosity and low permeability,challenging the accurate calculation of logging porosity.Based on the research of M RIL-P T2 spectrum inversion,we p resent a method of op timizing regularization parameters.The p rocessing results show that MRIL-P is capable of calculating accurate porosity of low permeability&porosity reservoirs.

GeoReveal Interp retation Platform,NMR logging,spectrum inversion,constraint methods,op timization,criteria,low po rosity&permeability reservoir

羅利,1971年生,博士研究生;主要從事測井解釋方法研究工作。地址:(400021)重慶市江北區大石壩大慶村400號。電話:(023)67352049。E-mail:luoli3988@163.com

羅利等.GeoReveal M RIL-P核磁解譜技術在低孔滲儲層評價中的應用.天然氣工業,2011,31(2):59-62.

10.3787/j.issn.1000-0976.2011.02.014

2010-12-03 編輯 韓曉渝)

DO I:10.3787/j.issn.1000-0976.2011.02.014

Luo L i,senior engineer,born in 1971,has long been engaged in research of log interp retation methods.

Add:No.400,Daqingcun,Dashiba,Jiangbei District,Chongqing 400021,P.R.China

Tel:+86-23-6735 2049 E-mail:luoliI3988@163.com