專利名稱：多屬性地震波形分類的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及對油、氣儲層的地震勘探，更具體地說，它涉及使用一種處理技術(shù)在高速數(shù)字計算機中處理地震數(shù)據(jù)，該處理技術(shù)包括對整個數(shù)據(jù)集(a data volume)中的完全波形形狀進行分類，并繪出由波形分類所揭示的多種屬性的區(qū)域變化趨勢，這些變化趨勢能代表類似的巖性和/或烴類效應(yīng)。
背景技術(shù)：
多年來，油、氣地震勘探已涉及使用地震能量源和用地震檢測器(通常稱作地震檢波器)陣列對其進行接收。當(dāng)在陸地上使用時，地震能量源可以是在地層上選定點處的井孔中電起爆的高爆炸力炸藥，或者說是有能力向地面發(fā)送一系列沖擊或機械振動的其他能量源。在海上，通常使用空氣槍源和水中地震檢波器接收器。由這些源在地球中產(chǎn)生的聲波從地層邊界和/或其他不連續(xù)面返回傳輸，并以不同的時間間隔達到地球表面，這時間間隔取決于穿過的距離和所經(jīng)過的地表特性。在陸地上，這些返回波被地震檢波器檢測到，這些地震檢波器的作用是把這些聲波轉(zhuǎn)換成代表聲波的電模擬信號，它們通常稱作記錄線(trace)。在陸地上使用時，一地震檢波器陣列沿著覆蓋感興趣地區(qū)的網(wǎng)格布設(shè)，以在所希望的地區(qū)內(nèi)形成一組間隔分布的觀測站，從而能構(gòu)建在廣大區(qū)域上的反射體位置的三維(3D)分布圖。與地震檢波器相隔一希望距離的地震源向地球內(nèi)注入聲波信號，且在該陣列中每個地震檢波器處檢測到的信號被記錄下來供以后用數(shù)字計算機處理，在那里模擬數(shù)據(jù)通常被量化為數(shù)字采樣點，例如每2毫秒一個采樣，從而可以對每個采樣點單獨進行操作。于是，連續(xù)記錄的地震現(xiàn)場(field)記錄線被轉(zhuǎn)換成近似于地下結(jié)構(gòu)的垂直剖面或體積表示、或水平分布圖。然后，地震檢波器陣列移到新的位置并重復(fù)這一過程以提供地震勘察結(jié)果。3D地震勘查結(jié)果是在表面收集的并展現(xiàn)為地下一部分的體積表示的數(shù)據(jù)。
在完成一個地區(qū)的勘探之后，已經(jīng)記錄了在多個地震檢波器處檢測到的與能量有關(guān)的數(shù)據(jù)，其中這些地震檢波器位于距爆點不同距離的地點。然后對數(shù)據(jù)進行重新組織，以收集在各爆點發(fā)射的和在各地震檢波器位置記錄的數(shù)據(jù)所產(chǎn)生的記錄線，其中這些記錄線被分組，使得可以假定這些反射是從地球內(nèi)一個特定點，即共中點反射的。然后，單個記錄或“記錄線”被加以不同的距離校正，該距離是地震能量從相應(yīng)爆點出發(fā)穿過地球到達共中點再向上到達各地震檢測器的距離。這一步驟包括因穿過不同類型巖石層的速度變化需要的校正以及源和接收器深度的變化需要的校正。由于爆點/地震檢波器對的間距變化所需要的校正稱作“正常時差(normal move out)”校正。在完成這一校正之后，對來自各中點的信號組進行相加。因為地震信號為正弦振動性質(zhì)，所以求和過程的作用是減小地震記錄中的噪聲，從而增大它的信號-噪聲比。這一過程稱作共中點數(shù)據(jù)的“疊加這對本領(lǐng)域技術(shù)人員是熟知的。于是地震現(xiàn)場數(shù)據(jù)受到上述校正，還可以受到偏移，這是對未經(jīng)解釋的數(shù)據(jù)進行的一種操作，涉及對地震信息的重新安排，從而使傾斜層畫在它們的真正位置。在把連續(xù)記錄的記錄線轉(zhuǎn)換成重直或水平剖面或水平分布圖之前，還可以應(yīng)用其他奇特的已知處理技術(shù)，這些技術(shù)例如能增強顯示斷層、地層特征、振幅隨與炮檢距(offset)的變化(AVO)或者某一屬性，如峰值振幅，瞬時頻率或相位、極性等。
在地震勘探過程中，可以由穿入感興趣地層的井孔建立控制點。然而，相當(dāng)經(jīng)常發(fā)生的情況是這些井孔間距很大，而且只有在這樣稀少的控制點處才能把選定的地震屬性與目標(biāo)地層的結(jié)構(gòu)(texture)和組分進行比較，以此對地震觀測進行校準(zhǔn)。地震勘查有相距較近的觀測點，它們分布在稀少的控制點之間，因此地震勘查具有潛力提供數(shù)據(jù)以估計那些遠離井孔延伸的儲層狀況。
普查地震分析經(jīng)常使用共同(common)地震波形的變化趨勢來識別有希望的地質(zhì)特征。許多當(dāng)前進行的處理使用簡單的屬性，如峰值振幅、平均頻率等，來顯示地區(qū)性趨勢，以峰值振幅應(yīng)用最普遍，因為它能代表一個“亮點”，該亮點可能是烴類的直接指示，更新近的一些方法已使用實際波形來顯示遠離井孔地方的地區(qū)性趨勢。然而，仍然需要包括利用范圍廣泛的記錄線到記錄線的相似性來進行波形形狀分析，它比較地震波形的多個屬性，以估計巖石類型或結(jié)構(gòu)的變化，即沿著給定地層的巖性變化和/或烴類效應(yīng)(如孔隙空間中的油或氣)的變化。
因此，希望允許解釋員能勘查多個屬性沿著要解釋的層面的變化。
還希望在多次地震疊加(例如近炮檢距和遠炮檢距的疊加)中分析記錄線波形形狀以及對波形分類，并繪制所發(fā)現(xiàn)的任何地震波形變化趨勢圖。
還希望允許由解釋員選擇若干個波形類別，這些類別被指定為對地震數(shù)據(jù)集中的記錄線進行分類。
還希望利用一種穩(wěn)健的方法去比較記錄線與記錄線的相似性，它對被比較波形的絕對振幅和形狀都予以考慮。
再有，還希望提取從地震數(shù)據(jù)集中能得到的附加巖性信息。

發(fā)明內(nèi)容
根據(jù)本發(fā)明，在一種方法和裝置中可實現(xiàn)上述的以及其他的目的和優(yōu)點，該方法和裝置使用選自記錄線集的地震記錄線子集來識別多個地震屬性沿一地下層面的變化趨勢。該方法包括在一預(yù)先確定的在感興趣的地下層面的上方和下方延伸的區(qū)域(或時間間隔)上將該子集中每個記錄線的全波形與該子集中的每個其他記錄線進行比較，從而確定該子集中每個記錄線的相似性因子。根據(jù)該子集中各單個記錄線的相似性因子的離散(Variance)范圍，將該子集中的記錄線分成若干波形類別。最后，根據(jù)該波形類別將該記錄線分組，從而顯示出該地下多個地震屬性的地區(qū)性變化趨勢。
在一優(yōu)選實施例中，該方法首先是選擇一被選排(gather)、疊加和偏移的地震記錄線集，其中標(biāo)識出被解釋的層面，該方法還進一步選擇一或多個其他處理過的記錄線集(例如近疊加、遠疊加、AVO截矩或梯度、相干性、P波段、S波段等)，這些記錄線集有識別出的同一層面。然后，根據(jù)本發(fā)明從要處理的每個記錄線集中選出相應(yīng)的代表性記錄線子集。一種穩(wěn)健的算法用于比較該子集中的記錄線對，其中該子集中的每個記錄線與該子集中的每個其他記錄線進行比較，該算法對被比較的兩個波形的形狀和絕對振幅都予以考慮。對于要處理的全部記錄線集的各相應(yīng)子集中的記錄線進行類似的比較。根據(jù)該子集中每個記錄線與該子集中每個其他記錄線的比較結(jié)果，確定該子集中每個記錄線的相似性因子，并對取自每個記錄線集的子集的相應(yīng)相似因子進行平均，從而形成一復(fù)合子集。可選地，可對被平均的相似性因子使用加權(quán)。
接下來，根據(jù)先前確定的對于復(fù)合子集的相似性因子離散范圍(arange of variance)，由解釋員(interpreter)選擇波形類別個數(shù)和要賦予一個類別的波形個數(shù)。形成該復(fù)合子集相似性因子矩陣，在一平均過程中確定相似性因子的最有共同性的類別和其后的逐個類別，該平均過程從該矩陣中迭代地去掉相似性因子，首先去掉的是最有共同性的波形類別的相似性因子。
然后，對于具有這些最有共同性波形的每個要處理的整個記錄線集中的每個記錄線，確定其相似性因子，并通過平均相應(yīng)的因子確定這些記錄線集的復(fù)合相似性因子，仍然是可選用加權(quán)。最后，根據(jù)最高的復(fù)合相似性因子對每個記錄線集中的每個記錄線賦予一波形類別，除非這個最高相似性因子低于一所希望的閾值。在這種情況下該記錄線被賦予一空(null)類別。為便于觀察，為整個記錄線集輸出波形類別圖或地震立方疊加圖。
按照本發(fā)明另一方面的裝置，包括一計算機，被編程以實現(xiàn)多個地震屬性沿地下層面變化趨勢的識別方法。對于大規(guī)模的3D勘探數(shù)據(jù)來說，該計算機最好是大型并行機。
本發(fā)明的方法和裝置對在一個時間(或深度)間隔在被解釋層面上方和下方延伸的各記錄線的全波形進行比較，并進一步進行比較，這些比較包括被解釋的層面所在3D體積中的全部記錄線。然后每個記錄線被分類，從而使類似的記錄線能分組在一起以供顯示，例如按照與該類中選出的第一記錄線匹配的質(zhì)量順序進行顯示。這樣，所顯示的記錄線突出了地震屬性的變化趨勢，它可以代表類似巖性和/或烴類效應(yīng)的變化趨勢。由下文的詳細描述和附圖，對于本領(lǐng)域的技術(shù)人員而言，本發(fā)明的其他目的和優(yōu)點將變得顯而易見，其中只顯示和描述本發(fā)明的一個最佳實施例。
附圖簡述

圖1是一地震記錄線集的透視圖，表示從記錄線中拾取的一被解釋的層面，這些記錄線從表面網(wǎng)格延伸并切入該層面。
圖2(a)-2(b)是本發(fā)明主要數(shù)據(jù)處理步驟的計算機軟件流程圖。
圖2(c)-2(d)為表示根據(jù)本發(fā)明在地下地層中波形分類變化趨勢圖。
圖3是流程圖，說明根據(jù)圖2中的步驟52決定波形分類的數(shù)據(jù)處理步驟。
圖4表示10個地震記錄線的子集，這些記錄線被分成兩類，并按照與該類中第一記錄線匹配的質(zhì)量順序顯示。
優(yōu)選實施例的詳細說明圖1代表一地震記錄線集，表示記錄線20按照要勘查地區(qū)上方的表面網(wǎng)格陣列排列。通常要勘查的地區(qū)將包括一個由上千個甚至更多記錄線覆蓋的區(qū)域。記錄線20已被垂直疊加，從地球表面22向下延伸切入地下層面24。該層面把地震波部分地反射回地面，在那里它們被檢測到。如圖所示，這些記錄線被安排在一個矩形網(wǎng)格中，該網(wǎng)格部分地覆蓋要勘查的地區(qū)。然而，如前所述，該網(wǎng)格陣列向前移動到一個新的位置并重復(fù)勘探過程以得到完全的勘查。如圖1中表示的地下層面24就將是由解釋員拾取的或用計算機軟件自動拾取的，這在本領(lǐng)域是眾所周知的。
現(xiàn)在參考圖2，圖中說明本發(fā)明的主要處理步驟，第一步驟30是選擇具有被解釋的層面的地震記錄線集?？赡馨ㄒ粋€被解釋的層面的幾類記錄線集在本領(lǐng)域是眾所周知的。各記錄線集可能有不同的獲取特性，例如，記錄線集可能是作為共中點、共炮檢距或共反射角范圍被聚集的。各記錄線集還可能有不同的數(shù)據(jù)處理特性，例如AVO梯度或AVO截距。其他特性可以區(qū)分記錄線集，如壓縮波或剪切波。根據(jù)本發(fā)明，如像前文提到的那樣，將處理多個地震記錄線集。
在步驟32中，選出一個有代表性的、切入層面24的、在空間上分離的記錄線子集。該子集可能有100個記錄線，它們包括感興趣的層面的延伸范圍。例如，該子集可能包括該網(wǎng)格陣列的均勻間距記錄線，例如每第5行和第5記錄線，它們是從具有上千或更多記錄線的記錄線集中選出的。
在步驟34中，確定該子集中每個記錄線與該子集中每個其他記錄線的相似性因子。波形相似性的最通用的度量是歸一化相關(guān)系數(shù)，它是兩個波形間線性關(guān)系的量度。這一量度的主要局限是它只考慮波形形狀而不考慮絕對振幅(即兩個波形能有極其不同的絕對振幅但有高相關(guān)系數(shù))。在絕對振幅也應(yīng)作為相似性測試中的一個屬性予以考慮的情況中，歸一化樣本差可與歸一化相關(guān)系數(shù)結(jié)合，從而形成一個新的度量，其中稱作RB因子。
歸一化相關(guān)系數(shù)是一個人們相當(dāng)熟知的度量。給定兩個長度為N的時間序列(Ai和Bi)，歸一化相關(guān)系數(shù)(X(a，b))按下式計算X(a,b)=Σi=1Nai·biΣi=1Nai2·Σi=1Nbi2]]>其中長度N是與包括層面24的選定時間或深度間隔相對應(yīng)的比較窗。
歸一化樣本差(D(a，b))定義為D(a,b)=Σi=1N(|ai|+|bi|)-Σi=1N|ai-bi|Σi=1N(|ai|+|bi|)]]>
這產(chǎn)生一個差值系數(shù)，它與相關(guān)系數(shù)在形式上的相似之處是兩個完全相同的波形將產(chǎn)生系數(shù)1.0。
然后，RB因子定義為這兩個度量的加權(quán)平均，即RB(a，b)＝WX(a，b)+(1-W)D(a，b)其中W是介于0和1之間的加權(quán)因子(通常把W設(shè)為0.5以成為等加權(quán))。
于是，對于每個子集的每個記錄線，根據(jù)與該子集中每個其他記錄線的比較，計算出RB因子。
編程的方法在步驟36詢問是否已處理了全部記錄線集。當(dāng)結(jié)果為“否”時，表明還可得到更多記錄線集，于是經(jīng)由塊38(它訪問下一個記錄線集)重復(fù)步驟30，并在步驟32形成相應(yīng)的子集。當(dāng)結(jié)果為“是”時，表明所有可得到的記錄線集都已處理過，于是在步驟40中對處理過的每個記錄線集的子集中相應(yīng)記錄線的RB因子進行平均，從而計算出復(fù)合RB因子。在由參考數(shù)字40表示的步驟中表明了這一處理，其中還表明，如果希望的話，可以對各因子加權(quán)。接下來程序進入步驟42，它要求限定若干波形類別以分割復(fù)合子集，這通常是由解釋員根據(jù)步驟40中計算出的多個復(fù)合RB因子的離散來確定的。
在把復(fù)合RB因子分成波形類別的準(zhǔn)備過程中，這些RB因子被排列在一正方矩陣中，如在步驟44中要求的那樣。
在步驟46中，使用一迭代過程把復(fù)合子集的復(fù)合RB因子分成所希望個數(shù)的類別，對于該迭代過程，將參考圖3做更充分的解釋。接下來在步驟48中，對整個記錄線集中的每一個記錄線計算RB因子(相似性)，在步驟50中對來自每個記錄線集的相應(yīng)單個因子進行平均(可選地采用先前使用的加權(quán))，以形成單個記錄線的復(fù)合RB因子。然后，使用該記錄線的復(fù)合RB因子將該記錄線與各診斷波形形狀(先前在步驟46中確定的)進行比較。在步驟52，根據(jù)該記錄線的最高復(fù)合RB因子對該記錄線賦予一先前確定的波形類別。如果該記錄線的最高復(fù)合RB因子低于一由使用者定義的閾值(例如0.75)，則對該記錄線賦予空值。最后，在步驟54中，為整個記錄線集輸出一波形類別圖或地震立方疊加圖。
圖2(c)以分布圖形式表示6個波形類別，根據(jù)圖上的編碼條對它們進行了編碼。圖2(d)是一個剖面圖，對應(yīng)于圖2(c)上的線AB，它表示沿著一被解釋的層面編碼的波形類別。
現(xiàn)在參考圖3，圖中說明確定波形類別的處理步驟。首先在步驟60中分別把數(shù)M和N定義為該矩陣中的RB因子數(shù)和要提供的波形類別數(shù)。于是，M對應(yīng)于為該子集選擇的記錄線數(shù)，如前所述，N是由解釋員根據(jù)該子集的RB因子離散范圍選定的。接下來在步驟62中計算出賦予每個類別的記錄線數(shù)NTR＝M/N。例如，如果在該子集中有100個記錄線，并希望有10個類別，那么在每個類別中將有10個記錄線。
然后被編程的方法進入步驟64，那里對該子集中M個記錄線中的每一個，選出(NTR-1)個其他記錄線，這些記錄線相對于所考慮的該記錄線有最高的RB因子(即最佳匹配)。接下來如在步驟66中要求的那樣，具有在先前步驟中計算出的最高平均RB因子的那個記錄線被選為第一類別中的第一記錄線。然后該方法進入步驟68，在該步驟中從矩陣中去掉先前步驟中選出的第一記錄線及其(NTR-1)個最佳匹配記錄線的RB因子。然后，在步驟70中重復(fù)上述步驟64-68，以確定其余類別。
形狀分類舉例這個例子以數(shù)字說明本方法如何用于對波形分類以能更好地顯示多個屬性變化趨勢。給定在圖2的步驟40中建立的復(fù)合相似性矩陣RB(i，j)為記錄線i對記錄線j的相似性，其中i和j的范圍是從1到M，找出最有共同性的波形形狀的過程如下對本例假定M＝10個記錄線，N＝2個類別。
首先計算每個類別中的記錄線數(shù)NTR每個類別中有10/2＝5個記錄線。
現(xiàn)在參見表I，其中以數(shù)字說明10個記錄線彼此比較所得到的復(fù)合RB因子，并以矩陣格式顯示，我們希望把這10個記錄線分成兩類。

接下來計算每個記錄線的最高(NTR-1)個(即4個)RB因子的平均值如下a.記錄線1的平均相似性＝(0.91+0.86+0.86+0.77)/4＝0.85b.記錄線2的平均相似性＝(0.90+0.86+0.86+0.75)/4＝0.8425c.記錄線3的平均相似性＝(0.91+0.87+0.86+0.83)/4＝0.8675d.記錄線4的平均相似性＝(0.90+0.87+0.86+0.79)/4＝0.855e.記錄線5的平均相似性＝(0.83+0.79+0.77+0.77)/4＝0.79f.記錄線6的平均相似性＝(0.78+0.77+0.63+0.63)/4＝0.7025g.記錄線7的平均相似性＝(0.83+0.78+0.67+0.65)/4＝0.7325
h.記錄線8的平均相似性＝(0.86+0.83+0.80+0.63)/4＝0.78i.記錄線9的平均相似性＝(0.91+0.86+0.67+0.61)/4＝0.7625j記錄線10的平均相似性＝(0.91+0.80+0.66+0.65)/4＝0.755選擇具有最高平均RB因子的記錄線，它在本例中是3號記錄線，而在這一類別中的其他(NTR-1)個記錄線是1、2、4和5。從相似性矩陣中去掉這些記錄線，我們得到下面的表II。

現(xiàn)在對第二類，我們再對每個剩下的記錄線計算最高的4個RB因子的平均值如下a.記錄線6的平均相似性＝(0.78+0.63+0.48+0.46)/4＝0.5875b.記錄線7的平均相似性＝(0.83+0.78+0.67+0.65)/4＝0.7325
c.記錄線8的平均相似性＝(0.86+0.83+0.80+0.63)/4＝0.78d.記錄線9的平均相似性＝(0.91+0.86+0.67+0.48)/4＝0.73e.記錄線10的平均相似性＝(0.91+0.80+0.65+0.46)/4＝0.705選擇記錄線8作為具有最高RB因子的記錄線，其他(NTR-1)個記錄線是6、7、9和10。
圖4顯示這10個記錄線實例，其時間軸以秒為單位。這些記錄線根據(jù)各自的類別顯示出來，其中這些記錄線按它們各自的類別分組，其排列順序是根據(jù)該類別中與第一記錄線的匹配質(zhì)量，即第一類的順序是3、1、4、2、5，第二類的順序是8、9、7、10、6。能對整個記錄線集輸出的波形類別分布圖或地震立方疊加圖，以在儲層狀況從井孔延伸到遠處時幫助識別這些儲層狀況的變化趨勢。
在所公開的內(nèi)容中只表示了本發(fā)明的最佳實施例?？梢灶A(yù)期，對這些實施例進行修改和替換，而且這些修改和替換可得到本發(fā)明的優(yōu)點和益處，這對于本領(lǐng)域那些已經(jīng)參考了這一說明及其附圖的技術(shù)人員而言是顯然的。應(yīng)該理解，這些修改和替換是在下面的權(quán)利要求書中申明的本發(fā)明的范圍之內(nèi)。
權(quán)利要求
1.一種使用一或多個地震記錄線集識別多個地震屬性沿地下層面的地區(qū)性變化趨勢的方法，所述方法包括a)從所述地震記錄線集當(dāng)中選出有代表性的記錄線子集，其中所述記錄線子集切入所述層面；b)在預(yù)先確定的在所述地下層面之上和之下的區(qū)域上將每個記錄線集的每個所述子集中的每個記錄線的全波形與所述子集中的每個其他記錄線進行比較，以確定在所述記錄線子集中每個記錄線與每個其他記錄線的相似性因子；c)根據(jù)所述子集中各記錄線的所述相似性因子的離散程度將所述子集中的記錄線分成若干波形類別；以及d)根據(jù)所述波形類別將所述子集中的記錄線分組，其中根據(jù)所述波形類別顯示的記錄線揭示出多個地震屬性的地區(qū)性變化趨勢。
2.根據(jù)權(quán)利要求1的方法，其中所述一或多個地震記錄線集是從一組地震記錄線集當(dāng)中選出的，它們在獲取特性方面彼此不同，所述組包含共中點集；共炮檢距集；共反射角集；近疊加集；以及遠疊加集。
3.根據(jù)權(quán)利要求1的方法，其中所述一或多個地震記錄線集是從一組地震記錄線集當(dāng)中選出的，它們在某種數(shù)據(jù)處理特性方面彼此不同，所述集組包含振幅對炮檢距(AVO)梯度集；以及AVO截距集。
4.根據(jù)權(quán)利要求1的方法，其中所述一或多個地震記錄線集是從一組地震記錄線集當(dāng)中選出的，它們在傳播的地震波類型方面彼此不同，所述集組包含由壓縮波產(chǎn)生的集；以及由剪切波產(chǎn)生的集。
5.根據(jù)權(quán)利要求1的方法，其中所述記錄線的代表性子集包含一組在網(wǎng)格陣列中有均勻間距的記錄線。
6.根據(jù)權(quán)利要求1的方法，其中所述相似性因子是用這樣一種算法確定的，該算法考慮了被比較的波形段的形狀和絕對振幅差兩個方面。
7.根據(jù)權(quán)利要求6的方法，其中所述算法使用如下表達式確定歸一化絕對振幅差D(a,b)=Σi=1N(|ai|+|bi|)-Σi=1N|ai-bi|Σi=1N(|ai|+|bi|)]]>而且此處所述子集的復(fù)合相似性因子是通過對所處理的全部記錄線集平均各子集中的相應(yīng)因子得到的。
8.根據(jù)權(quán)利要求1的方法，其中所述根據(jù)波形類別對所述子集中的記錄線進行分組的步驟包括(a)把所述子集的所述相似性因子排列在一正方矩陣中；(b)定義一個數(shù)，表示為M，它等于所述矩陣中的所述相似性因子的個數(shù)；(c)定義一個數(shù)，表示為N，它等于所希望的類別個數(shù)；(d)計算數(shù)M/N，表示為NTR，它等于待賦予每個類別的記錄線數(shù)；(e)對M個記錄線中的每一個，選擇(NTR-1)個其他具有最高相似性因子的記錄線并平均這(NTR-1)個相似性因子；(f)選擇具有最高平均相似性因子的記錄線作為第一類的形狀；(g)從所述矩陣中去掉有最高相似性因子的記錄線以及那其他(NTR-1)個記錄線；以及(h)重復(fù)步驟(e)、(f)和(g)以確定其余類別中的記錄線。
9.使用一或多個地震記錄線集識別多個地震屬性沿地下層面的地區(qū)性變化趨勢的裝置，所述裝置包括一計算機，被編程以實現(xiàn)下述方法步驟(a)從所述地震記錄線集當(dāng)中選出有代表性的記錄線子集，其中所述記錄線子集切入所述層面；(b)在預(yù)先確定的在所述地下層面之上和之下的區(qū)域上將每個記錄線集的每個所述子集中的每個記錄線的全波形與所述子集中的每個其他記錄線進行比較，以確定在所述記錄線子集中每個記錄線與每個其他記錄線的相似性因子；(c)根據(jù)所述子集中各記錄線的所述相似性因子的離散程度將所述子集中的記錄線分成若干波形類別；以及(d)根據(jù)所述波形類別將所述子集中的記錄線分組，其中根據(jù)所述波形類別顯示的記錄線揭示出多個地震屬性的地區(qū)性變化趨勢。
10.根據(jù)權(quán)利要求9的裝置，其中所述一或多個地震記錄線集是從一組地震記錄線集當(dāng)中選出的，它們在獲取特性方面彼此不同，所述組包含共中點集；共炮檢距集；共反射角集；近疊加集；以及遠疊加集。
11.根據(jù)權(quán)利要求9的裝置，其中所述一或多個地震記錄線集是從一組地震記錄線集當(dāng)中選出的，它們在某種數(shù)據(jù)處理特性方面彼此不同，所述集組包含振幅對炮檢距(AVO)梯度集；以及AVO截距集。
12.根據(jù)權(quán)利要求9的裝置，其中所述一或多個地震記錄線集是從一組地震記錄線集當(dāng)中選出的，它們在傳播的地震波類型方面彼此不同，所述集組包含由壓縮波產(chǎn)生的集；以及由剪切波產(chǎn)生的集。
13.根據(jù)權(quán)利要求9的裝置，其中所述記錄線的代表性子集包含一組在網(wǎng)格陣列中有均勻間距的記錄線。
14.根據(jù)權(quán)利要求9的裝置，其中所述根據(jù)波形類別對所述子集中的記錄線進行分組的步驟包括(a)把所述子集的所述相似性因子排列在一正方矩陣中；(b)定義一個數(shù)，表示為M，它等于所述矩陣中所述相似性因子的個數(shù)；(c)定義一個數(shù)，表示為N，它等于所希望的類別個數(shù)；(d)計算數(shù)M/N，表示為NTR，它等于待賦予每個類別的記錄線數(shù)；(e)對M個記錄線中的每一個，選擇(NTR-1)個其他具有最高相似性因子的記錄線并平均這(NTR-1)個相似性因子；(f)選擇具有最高平均相似性因子的記錄線作為第一類的形狀；(g)從所述矩陣中去掉有最高相似性因子的記錄線以及那其他(NTR-1)個記錄線；以及(h)重復(fù)步驟(e)、(f)和(g)以確定其余類別中的記錄線。
全文摘要
公開了一種由計算機實現(xiàn)的方法(圖2a、2b和3),用于分析地震數(shù)據(jù)以發(fā)現(xiàn)地震數(shù)據(jù)集中多個屬性的地區(qū)性變化趨勢。根據(jù)所公開的方法和裝置,在一被解釋的平面(24)上方和下方延伸的區(qū)域中對多個記錄線的完整波形(圖4)進行比較(34),這種比較是使用一種相關(guān)算法進行的,該算法考慮被比較波形的形狀和絕對振幅兩者?；诙啻伪容^,根據(jù)波形形狀對記錄線進行分類(40,52),于是共同波形類別分布圖顯示出具有相似巖性和/或烴類效應(yīng)的地區(qū)。
文檔編號G01V1/28GK1382262SQ00814563
公開日2002年11月27日 申請日期2000年10月11日 優(yōu)先權(quán)日1999年10月20日
發(fā)明者丹尼斯·B·奈夫, 斯科特·A·朗納斯特蘭, 埃德加·L·布特勒 申請人:菲利浦石油公司