專利名稱：核磁共振測井t2譜t2截止值的確定方法
技術領域：
本發明涉及一種核磁共振測井T2譜T2截止值的確定方法，尤其是一種對核磁測井T2譜采用形態特征分類并逐點確定出每個采樣點上測井T2譜的T2截止值的方法。
背景技術：
核磁共振測井的原始數據是隨時間衰減的自旋-回波串，自旋-回波串包含了儲層物性、孔隙類型、孔徑大小、流體類型及其分布等十分豐富的信息。自旋-回波串可以表示為多個衰減指數的迭加Echo(t)=Σi=1Nφiexp(-t/T2i)---(1)]]>式中Echo(t)----核磁測井測量得到的隨時間t變化的回波幅度；T2i---第i個組分的橫向弛豫時間，i＝1，…，N；Фi---對應于T2i組分的核磁共振孔隙度，i＝1，…，N；通過對核磁測井儀器測量得到的自旋-回波串(即式(1))進行多指數擬合，可以求出每一測量點的T2i--Ф1的分布曲線，即所謂的核磁共振測井T2分布譜。該T2分布譜的橫坐標為T2值，縱坐標為信號相對幅度，總的信號幅度與巖石孔隙度有關。在單相流體條件下，T2分布定性反映孔隙大小分布。即大孔徑孔隙對應較大的T2值，小孔徑孔隙對應較小的T2值。
核磁共振測井是目前唯一一種可以區分孔隙中流體賦存狀態的測井方法，因此具有獨特的技術優勢。通過對核磁測井T2譜進行處理，可以得到巖石孔隙大小及其分布，從而更加直接的確定儲層的束縛水飽和度和滲透率等。然而與其它測井新技術(MDT、FMI等)相比，核磁測井技術無論從應用的深度和發揮的地質效果上都還存在很大的差距，其中的一個重要的因素就是現行的T2截止值的確定方法還存在應用上的困難，核磁測井T2截止值的確定技術存在缺陷和局限性。
T2截止值是核磁測井資料分析與資料應用的一個重要的基礎參數，它是計算儲層可動流體或束縛流體含量的關鍵參數，也是定量評價儲層孔隙結構、滲透率等的基礎參數。所謂T2截止值，是指T2譜上存在一個確定界限值，大于這個數值所對應的孔隙中的流體稱之為自由流體，在儲層條件下是可以流動的。小于這個數值對應的孔隙流體稱之為(毛管)束縛流體，是不能流動的。即BVI=ΣT2≤T2cutoffφi---(2)]]>FFI=ΣT2>T2cutoffφi---(3)]]>式中，T2cutoff——T2截止值；BVI——毛管束縛水孔隙體積；FFI——自由流體孔隙體積(或稱自由流體指數)。
T2截止值作為表征儲層毛管孔隙中流體特性的一個基礎參數，可以通過實驗室測量確定。國外的研究人員通過大量巖心(主要來自北美地區)實驗測量，認為碎屑巖儲層巖石的T2截止值為33ms，碳酸鹽巖儲層巖石為92ms(Morriss等，1993)，即T2截止值在已知巖性的情況下基本上是一個常數數值。事實上，目前國內外基本上也是采用這兩個截止值對核磁測井資料進行處理和應用的。
自從上世紀90年代中期核磁測井技術引入到中國以來，關于巖心T2截止值的核磁共振實驗室測量結果卻不支持國外同行普遍認同的T2截止值是某個常數數值的結論。特別是在碎屑巖儲層中，T2截止值的變化范圍相當大，根本不是一個定值。T2截止值隨地區或儲層不同而差異很大。圖1為我國部分油田砂巖儲層根據巖心核磁共振實驗得到的T2截止值數據匯總，從圖1中可以明顯看出T2截止值并不是一個固定的值，也看不出數據有趨向于33ms(碎屑巖地層)的跡象，對國外適用的33ms(碎屑巖地層)在我國的實際地層中變得不適用了。
目前在國內，大家基本上認可了砂巖碎屑巖儲層的T2截止值不是33ms這樣一個定值的客觀事實，但是關于砂巖儲層的T2截止值為什么會發生如此大的變化、是什么因素控制了砂巖儲層的T2截止值的變化規律仍沒有明確的認識。因此導致了核磁測井資料處理和應用的一個難題和尷尬局面盡管已經認識到砂巖儲層的T2截止值不是一個定值，但由于無法找到一個確定T2截止值的適用方法，在核磁測井的資料處理時，仍然無奈地沿用國外的做法，即仍然使用33ms作為砂巖儲層統一的T2截止值。這樣做的一個嚴重后果就是，技術人員對核磁測井與T2截止值有關的處理結果持懷疑態度，因此制約了核磁測井技術的應用。
目前已發表的文獻中有關求取T2截止值的方法有以下幾種①.巖心實驗室核磁共振實驗法(譜面積法)。如圖2所示，先對飽和水的巖心進行核磁共振測量，得到橫向弛豫時間的飽和T2分布譜；然后在一定的壓力條件下將飽和水的巖心在離心機上甩出其中的自由水，只剩下束縛水后，再對巖心進行核磁共振測量得到離心后的T2分布譜，兩者相比較。在測井T2譜中從低端(小的T2值)開始進行面積積分S1，當積分面積S1＝S2(離心譜總面積)時，測井T2譜中對應T2值即為T2截止值(記為T2cutoff)。這種方法是根據T2截止值的定義所確定的，也是準確的確定T2截止值的方法。離心譜面積與測井T2譜面積的比值為束縛水飽和度(用符號Swi表示)。
②.毛管壓力曲線刻度法。首先測量巖心的孔隙度和毛管壓力曲線，由毛管壓力曲線確定近似于油藏的束縛水飽和度Swir，而巖心的束縛流體體積為φSwir，然后按圖3的方法刻度T2截止值。按圖4所示的迭代圖解法優化T2截止值。圖4中，橫坐標為油層測井深度范圍巖心刻度T2截止值，縱坐標為同一油層深度范圍的T2譜和用常規巖心NMR分析確定的T2截止值所求的BVI之間的累計誤差。最小累計誤差所對應的T2截止值為地層最佳T2截止值。進行巖樣刻度時，應該收集地區有代表性取心資料，并將取心資料進行深度歸位處理。因取心資料在深度上可能存在誤差，所以應與經過校正的測井資料進行仔細對比。
③.滲透率交會圖法。在遇到取心困難，取心和巖心分析時間長以及不易獲得有代表性的巖心樣品時，可利用滲透率交會圖技術確定地層的T2截止值。由NMR測井深度范圍內的測井曲線獲得一系列T2截止值。通過改變T2截止值并計算其相應的Timur-Coates滲透率值(KNMR)，然后與同一深度范圍MDT(模塊式地層測試儀)壓降測量得出的滲透率值(KMDT)繪制交會圖，選擇交會圖上線性相關性最好的NMR測井滲透率相應的T2截止值為最佳T2截止值。
④.其他數學統計方法。參見2004年2月《石油地球物理勘探》雜志中，高楚橋，何宗斌等著作的《核磁共振T2截止值與毛管壓力的關系》的文章，其中記載了T2截止值與毛管壓力Pc之間存在指數關系，得出T2c＝A/PcB。參見2004年9月《地球物理學進展》雜志中，汪中浩，章成廣，肖承文等著作的《低滲透儲層T2截止值實驗研究》的文章，其中記載了T2截止值與巖石物性參數 之間存在一定的關系，得出T2c=a×(kφ)2+b×(kφ)+c.]]>上一部分所論述的確定T2截止值各種方法具有各自的優勢和局限性。前兩種方法的優點是可以獲得準確的T2截止值。但其應用上得局限性是需要巖心樣品，都需要巖心實驗分析或要用到其它的輔助測井方法，成本高、周期長。由于層段的限制，當有限的巖心實驗數據不存在良好的規律(如一批實驗數據變化規律很復雜或根本沒有什么變化規律可尋)，即使得到了這些實驗數據也無法應用到測井資料的連續處理。方法3和方法4都是基于地區特征得出的經驗公式，難以應用于其他地區，而且公式中的參數也不是唯一的，視地層情況不同而變化?？傊?，現行確定T2截止值各種方法的一個明顯的局限性就是需要巖心實驗或其他來源的信息進行標定。況且不論很多時候根本找不出某井或某儲層T2截止值實驗或標定結果的統計規律，即使能做到這點，也是地區規律并且往往存在很大誤差，無法普遍推廣與使用，而且其可靠程度與精度也得不到保證。
為了解決核磁測井應用中T2截止值準確確定的實際難題，應該采用全新的研究思路與方法，即避開巖心實驗或其他信息標定的傳統思路與方法，研究直接根據核磁共振測井T2譜(飽和T2譜)的形態特征確定每個測量點的T2截止值。提供確定核磁共振測井T2截止值的一種可靠方法和技術，克服現行方法的技術局限性。推動核磁測井技術的應用，更大程度地發揮其技術優勢，使這項技術在油田勘探開發中發揮重要的作用。通過研究，我們認為在只有核磁測井T2譜的情況下，是能夠根據測井T2譜的形態來確定每個深度點的T2截止值的。
發明內容
本發明的目的是直接利用核磁共振測井得到的T2譜而確定每個深度點的T2截止值。
為實現上述目的，本發明提供了一種核磁共振測井T2譜T2截止值的確定方法，包括以下步驟步驟1、對核磁共振測井T2譜按形態特征進行分類；步驟2、根據不同類型的測井T2譜擬合得出其離心譜；步驟3、將離心譜進行面積積分，得出離心譜的面積；步驟4、計算測井T2譜的面積；步驟5、將測井T2譜的面積進行累加，當累加面積等于或與離心譜面積最相近時，此時對應的T2值即為T2截止值。
所述步驟5之后還具有步驟6、將離心譜的面積除以測井T2譜的面積，即為束縛水飽和度的值。所述核磁共振井T2譜分為雙峰型、單峰型和三峰型。所述雙峰型分為雙峰小孔型、雙峰大孔型，單峰型分為單峰小孔型、單峰大孔型。所述步驟2中單峰大孔型測井T2譜T2截止值的計算方法為利用數學函數對大孔段進行擬合，其數學形式為f(x)=A×4×(1-ex-BC)×ex-BC]]>其中A-------擬合T2譜中的幅值極大值，A值取T2分布譜中的幅值極大值；B-------某個T2時間值，B值取幅度值減小到某個極小時的T2值；C-------某個待定的參數；x-----各組份的T2時間值與測井T2譜的T2時間值相同；f(x)-----擬合得到的幅度值；所述步驟2中單峰小孔型測井T2譜T2截止值的計算方法為，假設孔隙中的流體都是束縛流體，將孔隙賦予一個固定的T2截止值。所述固定的T2截止值為33ms。
所述步驟2中雙峰型測井T2譜T2截止值的計算方法為，利用正態函數在測井T2譜上擬合出一個離心譜，其數學形式為f(x)=A×e-(x-μ)22σ2]]>其中A為擬合T2譜幅度值，A值取T2分布譜中第一個峰的幅值，即使擬合的函數圖形和T2分布譜有相同的幅度；μ為擬合T2譜的期望值，μ值取第一個峰值所對應的T2時間值；σ為擬合T2譜的方差；
x為各組份的T2時間值與測井譜的T2時間值相同；f(x)為擬合得到的幅度值；所述步驟2中三峰型測井T2譜T2截止值的計算方法為利用正態函數在測井T2譜上擬合出一個離心譜，其數學形式為f(x)=A×e-(x-μ)22σ2]]>其中A為擬合T2譜幅度值，A值取T2分布譜中第一個峰的幅值，即使擬合的函數圖形和T2分布譜有相同的幅度；μ為擬合T2譜的期望值，μ值取第一個峰值所對應的T2時間值；σ為擬合T2譜的方差；x為各組份的T2時間值與測井譜的T2時間值相同；f(x)為擬合得到的幅度值；所述面積的計算方法為譜龍貝格求積法進行面積積分。
所述步驟5中T2截止值的計算方法為將測井T2譜的面積進行累加，當累加面積等于或與離心譜面積最相近時，此時對應的T2值即為T2截止值。上述各種類型均采用此種方法計算T2截止值。
因此，本發明實現了直接利用核磁共振測井得到的T2譜而確定每個深度點的T2截止值。
下面結合附圖和實施例，對本發明的技術方案做進一步的詳細描述。
圖1為新疆、華北以及西北部分地區的巖心實驗的T2截止值統計圖。
圖2為離心法確定T2截止值示意圖。
圖3為毛管壓力曲線刻度法確定T2截止值示意圖。
圖4為BVI累計誤差與T2截止值關系曲線圖。
圖5為本發明核磁共振測井T2譜T2截止值的確定方法雙峰小孔型巖心核磁共振測量結果示意圖。
圖6為本發明核磁共振測井T2譜T2截止值的確定方法雙峰小孔型實際核磁測井儀器測量結果示意圖。
圖7為本發明核磁共振測井T2譜T2截止值的確定方法雙峰大孔型巖心核磁共振測量結果示意圖。
圖8為本發明核磁共振測井T2譜T2截止值的確定方法雙峰大孔型實際核磁測井儀器測量結果示意圖。
圖9為本發明核磁共振測井T2譜T2截止值的確定方法單峰大孔型巖心核磁共振測量結果示意圖。
圖10為本發明核磁共振測井T2譜T2截止值的確定方法單峰大孔型實際核磁測井儀器測量結果示意圖。
圖11為本發明核磁共振測井T2譜T2截止值的確定方法單峰小孔型巖心核磁共振測量結果示意圖。
圖12為本發明核磁共振測井T2譜T2截止值的確定方法單峰小孔型實際核磁測井儀器測量結果示意圖。
圖13為本發明核磁共振測井T2譜T2截止值的確定方法三峰型巖心核磁共振測量結果示意圖。
圖14為本發明核磁共振測井T2譜T2截止值的確定方法三峰型實際核磁測井儀器測量結果示意圖。
圖15為本發明核磁共振測井T2譜T2截止值的確定方法實施例雙峰型示意圖1。
圖16為本發明核磁共振測井T2譜T2截止值的確定方法實施例雙峰型示意圖2。
圖17為本發明核磁共振測井T2譜T2截止值的確定方法實施例單峰型示意圖。
圖18為本發明核磁共振測井T2譜T2截止值的確定方法實施例三峰型示意圖。
圖19為本發明核磁共振測井T2譜T2截止值的確定方法新疆油田某井實際處理的效果圖。
圖20為本發明核磁共振測井T2譜T2截止值的確定方法西北地區某井實際處理的效果圖。
圖21為本發明核磁共振測井T2譜T2截止值的確定方法新疆等四個地區計算的T2截止值與巖心分析的T2截止值結果對比圖。
圖22為本發明核磁共振測井T2譜T2截止值的確定方法新疆等四個地區計算的束縛水飽和度與巖心分析的束縛水飽和度結果的對比圖。
圖23為本發明核磁共振測井T2譜T2截止值的確定方法計算機計算流程示意圖。
圖24為本發明核磁共振測井T2譜T2截止值的確定方法流程示意圖。
具體實施例方式
通過我們多年的研究，在大量實驗數據和實測資料的基礎上，提出了一種基于核磁測井T2譜形態特征來確定T2截止值的方法，簡稱形態法。形態法確定T2截止值的流程示意圖如圖24所示第一步，對核磁共振測井的T2譜按形態特征進行分類。
由于不同孔隙結構的巖心飽和譜與離心譜特征差異較大，經過對比分析，將T2飽和譜分成以下三大類五小類，分別為a.雙峰型(又分為雙峰大孔型和雙峰小孔型)；b.單峰型(又分為單峰大孔型和單峰小孔型)；c.三峰型。對應不同的孔隙結構類型，計算T2截止值的方法也有差別。對于低孔低滲儲層巖心，T2分布譜形態主要為單峰小孔型、雙峰小孔型、雙峰大孔型，對于高孔、中-高滲儲層巖石T2分布譜則通常以單峰大孔型為主。
各類譜的特征和形態描述如下a.雙峰型(亞類雙峰小孔型、雙峰大孔型)雙峰小孔型如圖5和圖6所示，為雙峰小孔型巖心核磁共振測量示意圖和實際核磁測井儀器測量示意圖。T2分布譜呈雙峰分布，在整個T2分布譜上有兩個幅值極大值，且第一個峰的幅值要大于第二個峰的幅值，由于小孔隙所對應的T2時間(第一個峰值點的T2值)要小于大孔隙的T2時間(第二個峰值點的T2值)，所以該類型被稱為雙峰小孔型。
雙峰大孔型如圖7和圖8所示，為雙峰大孔型巖心核磁共振測量示意圖和實際核磁測井儀器測量示意圖。T2分布譜呈雙峰分布，在整個T2分布譜上有兩個幅值極大值，且第一個峰的幅值要小于第二個峰的幅值，由于小孔隙所對應的T2時間(第一個峰值點的T2值)要小于大孔隙的T2時間(第二個峰值點的T2值)，所以該類型被稱為雙峰大孔型。
b.單峰型(亞類單峰大孔型、單峰小孔型)單峰大孔型如圖9和圖10所示，為單峰大孔型巖心核磁共振測量示意圖和實際核磁測井儀器測量示意圖。T2分布譜呈單峰分布，在整個T2分布譜上只有一個幅值極大值，且該極大值所對應的T2時間比較靠后(一般大于100ms)，也就是在大孔隙段的某個T2值達到極大，所以稱為單峰大孔型。
單峰小孔型如圖11圖12所示，為單峰小孔型巖心核磁共振測量示意圖和實際核磁測井儀器測量示意圖。T2分布譜呈單峰分布，在整個T2分布譜上只有一個幅值極大值，且該極大值所對應的T2時間比較靠前(一般大于100ms)，也就是在小孔隙段的某個T2值達到極大，所以稱為單峰小孔型。
C.三峰型如圖13和圖14所示，為三峰型巖心核磁共振測量示意圖和實際核磁測井儀器測量示意圖。
T2分布譜呈三峰分布，在整個T2分布譜上有三個幅值極大值，且最大的峰值不定，可能是第一個峰，可能是第二個峰，也可能是第三個峰。
第二步，針對各類T2分布譜，采用不同的方法確定其T2截止值，具體確定方法如下
a.雙峰型的T2截止值的確定方法利用正態函數在測井T2譜上擬合出一個離心譜，其數學形式為f(x)=A×e-(x-μ)22σ2]]>式中A-----擬合T2譜幅度值；μ-----擬合T2譜的期望值；σ-----擬合T2譜的方差；x-----各組份的T2時間值(和測井T2譜的T2時間值相同)；f(x)-----擬合得到的幅度值。
對于要準確確定T2截止值，則上述A、μ、σ三個參數的確定至關重要，我們的確定方法如下A值取T2分布譜中第一個峰的幅值，即使擬合的函數圖形和T2分布譜有相同的幅度；μ值取第一個峰值所對應的T2時間值；σ值和以第一個峰值的一半所對應的T2分布譜上兩個T2值的寬度有一定的相關性。
b.單峰大孔型的T2截止值的確定方法利用數學函數對大孔段進行擬合，其數學形式為f(x)=A×4×(1-ex-BC)×ex-BC]]>式中A-------擬合T2譜中的幅值極大值；B-------某個T2時間值；C-------某個待定的參數；x-----各組份的T2時間值(和測井T2譜的T2時間值相同)；f(x)-----擬合得到的幅度值。
如前所述，上述A、B、C三個參數的確定對確定T2截止值也是至關重要。研究發現，A值取T2分布譜中的幅值極大值，B值取幅度值減小到某個極小時的T2值，C值也和以第一個峰值的一半所對應的T2分布譜上兩個T2值的寬度有一定的相關性。
c.單峰小孔型的T2截止值的確定方法單峰小孔型所對應的孔隙結構多為微毛細管孔隙或互不連通的死孔隙，流體被完全封閉在其中，不能流動，這樣的孔隙結構中的流體是不能被采出的或采出量非常有限，在實際生產中也沒有工業價值。因此，我們可以給這類孔隙賦予一個固定的T2截止值例如33ms，認為這種孔隙中的流體都是束縛流體。
d.三峰型的T2截止值的確定方法三峰型的T2截止值的確定采用和雙峰型相同的正態函數來求取，其方法如雙峰型的求取方法。
有了以上各分類的算法后，用面積積分法來計算T2截止值以及束縛水飽和度①.按上述方法擬合得出離心譜；②.將離心譜按龍貝格(Romberg)求積法進行面積積分，得出的值相當于離心譜的面積；③.用同樣的算法計算測井T2譜的面積；④.將測井T2譜的面積進行累加，當累加面積等于或與離心譜面積最相近時，此時對應的T2值即為T2截止值；⑤.將離心譜的面積除以測井T2譜的面積，得到的值即為束縛水飽和度。
按照以上的思路，我們用Visual Fortran編寫了該方法的程序，實現了自動處理。流程圖如圖23所示。
利用形態法我們對新疆等地區的部分巖心做了計算，各分類的擬合效果圖如下雙峰型如圖15和16所示，對小孔段直接擬合，積分求出其面積即為束縛流體的面積，利用該面積除以總譜面積得到的是束縛水飽和度；對總譜從小T2值開始做面積積分，當累計面積和束縛流體面積相等時，此時的T2值即為T2截止值。
單峰型如圖17所示，對大孔段擬合后，積分求出其面積，然后用總譜面積減去大孔段面積，得到的即為束縛流體的面積，利用該面積除以總譜面積得到的是束縛水飽和度；對總譜從小T2值開始做面積積分，當累計面積和束縛流體面積相等時，此時的T2值即為T2截止值。
三峰型如圖18所示，對小孔段直接擬合，積分求出其面積即為束縛流體的面積，利用該面積除以總譜面積得到的是束縛水飽和度；對總譜從小T2值開始做面積積分，當累計面積和束縛流體面積相等時，此時的T2值即為T2截止值。
利用形態法自動確定核磁共振測井T2截止值方法以及自主開發的軟件包，我們對新疆油田某井和西北地區某井進行了處理，如圖19和20所示為新疆油田某井實際處理的效果圖和西北地區某井實際處理的效果圖，處理的效果圖如下圖中第一道為CALI(井徑曲線)和GR(自然伽馬曲線)；第二道為深度；第三道中SWICAL為利用形態法計算得到的束縛水飽和度曲線，SWI33為以33ms作為T2截止值時計算得到的束縛水飽和度曲線，CSWI為巖心分析得到的束縛水飽和度；第四道中T2_AMP為核磁共振測井得到的T2譜，T2CUTCAL為利用形態法計算得到的T2截止值曲線，T2CUT33為以33ms作為T2截止值時計算得到的T2截止值曲線，T2CUT_CORE為巖心分析得到的T2截止值。從圖上可以看出，利用形態法計算得到的T2截止值以及利用截止值計算的束縛水飽和度要好于以固定值33ms作為T2截止值的結果，形態法計算的結果跟接近巖心分析的結果，這說明，可以在沒有取心的情況下，直接利用形態法來確定T2截止值而不需依賴巖心實驗。
和現有的確定T2截止值的方法相比，形態法可以逐點計算出測井井段每個深度點上的T2截止值，而傳統方法由于是對巖心做實驗，所以只能是確定出取心點處的T2截止值。利用上述形態法以及自主開發的軟件包，我們就可以對核磁共振測井數據的每個深度點上的測井T2譜確定出其T2截止值，這樣就得到了一條連續的和測井井段相匹配的T2截止值曲線，進而利用該曲線可以計算出每個深度點上的束縛水飽和度以及滲透率等。利用形態法，我們對新疆等四個地區的巖心做了計算，圖21中的兩條邊界線分別是以±10ms為誤差的分界線，圖22中的兩條邊界線分別是以±5％為誤差的分界線，從圖21和圖22上可以看出利用形態法的結果可以很好的代替巖心分析的結果。
最后所應說明的是，以上實施例僅用以說明本發明的技術方案而非限制，盡管參照較佳實施例對本發明進行了詳細說明，本領域的普通技術人員應當理解，可以對本發明的技術方案進行修改或者等同替換，而不脫離本發明技術方案的精神和范圍，其均應涵蓋在本發明的權利要求范圍當中。
權利要求
1.一種核磁共振測井T2譜T2截止值的確定方法，其特征在于包括以下步驟步驟1、對核磁共振測井T2譜按形態特征進行分類；步驟2、根據不同類型的測井T2譜擬合得出其離心譜；步驟3、將離心譜進行面積積分，得出離心譜的面積；步驟4、計算測井T2譜的面積；步驟5、將測井T2譜的面積進行累加，當累加面積等于或與離心譜面積最相近時，此時對應的T2值即為T2截止值。
2.根據權利要求1所述的核磁共振測井T2譜T2截止值的確定方法，其特征在于所述步驟5之后還具有步驟6、將離心譜的面積除以測井T2譜的面積，即為束縛水飽和度的值。
3.根據權利要求1所述的核磁共振測井T2譜T2截止值的確定方法，其特征在于所述核磁共振測井T2譜分為雙峰型、單峰型和三峰型。
4.根據權利要求3所述的核磁共振測井T2譜T2截止值的確定方法，其特征在于所述雙峰型分為雙峰小孔型、雙峰大孔型，單峰型分為單峰小孔型、單峰大孔型。
5.根據權利要求4所述的核磁共振測井T2譜T2截止值的確定方法，其特征在于所述步驟2中單峰大孔型測井T2譜T2截止值的計算方法為利用數學函數對大孔段進行擬合，其數學形式為f(x)=A×4×(1-ex-BC)×ex-BC]]>其中A-------擬合T2譜中的幅值極大值，A值取T2分布譜中的幅值極大值；B-------某個T2時間值，B值取幅度值減小到某個極小時的T2值；C-------某個待定的參數；x-----各組份的T2時間值與測井T2譜的T2時間值相同；f(x)-----擬合得到的幅度值。
6.根據權利要求5所述的核磁共振測井T2譜T2截止值的確定方法，其特征在于所述步驟2中單峰大孔型測井T2譜T2截止值的計算方法為，將測井T2譜的面積減去擬合譜面積即為離心譜面積。
7.根據權利要求4所述的核磁共振測井T2譜T2截止值的確定方法，其特征在于所述步驟2中單峰小孔型測井T2譜T2截止值的計算方法為，假設孔隙中的流體都是束縛流體，將孔隙賦予一個固定的T2截止值。
8.根據權利要求7所述的核磁共振測井T2譜T2截止值的確定方法，其特征在于所述固定的T2截止值為33ms。
9.根據權利要求4所述的核磁共振測井T2譜T2截止值的確定方法，其特征在于所述步驟2中雙峰型測井T2譜T2截止值的計算方法為，利用正態函數在測井T2譜上擬合出一個離心譜，其數學形式為f(x)=A×e-(x-μ)22σ2]]>其中A為擬合T2譜幅度值，A值取T2分布譜中第一個峰的幅值，即使擬合的函數圖形和T2分布譜有相同的幅度；μ為擬合T2譜的期望值，μ值取第一個峰值所對應的T2時間值；σ為擬合T2譜的方差；x為各組份的T2時間值與測井T2譜的T2時間值相同；f(x)為擬合得到的幅度值。
10.根據權利要求4所述的核磁共振測井T2譜T2截止值的確定方法，其特征在于所述步驟2中三峰型測井T2譜T2截止值的計算方法為利用正態函數在測井T2譜上擬合出一個離心譜，其數學形式為f(x)=A×e-(x-μ)22σ2]]>其中A為擬合T2譜幅度值，A值取T2分布譜中第一個峰的幅值，即使擬合的函數圖形和T2分布譜有相同的幅度；μ為擬合T2譜的期望值，μ值取第一個峰值所對應的T2時間值；σ為擬合T2譜的方差；x為各組份的T2時間值與測井T2譜的T2時間值相同；f(x)為擬合得到的幅度值。
11.根據權利要求1所述的核磁共振測井T2譜T2截止值的確定方法，其特征在于所述步驟5中的T2截止值的計算方法為，將測井T2譜的面積進行累加，當累加面積等于或與離心譜面積最相近時，此時對應的T2值即為T2截止值。
12.根據權利要求1所述的核磁共振測井T2譜T2截止值的確定方法，其特征在于所述面積的計算方法為龍貝格求積法進行面積積分。
全文摘要
本發明涉及一種核磁共振測井T2譜T2截止值的確定方法，包括以下步驟步驟1、對核磁共振測井T2譜按形態特征進行分類；步驟2、根據不同類型的測井T2譜擬合得出其離心譜；步驟3、將離心譜進行面積積分，得出離心譜的面積；步驟4、計算測井T2譜的面積；步驟5、將測井T2譜的面積進行累加，當累加面積等于或與離心譜面積最相近時，此時對應的T2值即為T2截止值。由此，本發明實現了直接利用核磁共振測井得到的T2譜而確定每個深度點的T2截止值。
文檔編號G01V3/18GK1869733SQ20051007223
公開日2006年11月29日 申請日期2005年5月27日 優先權日2005年5月27日
發明者匡立春, 毛志強, 孫中春, 王勁 申請人:中國石油天然氣股份有限公司, 石油大學(北京)
專業數控銑床產品供應商
20年生產經驗，實力工廠，精工品質，質量有保障
