用于地震數(shù)據(jù)處理的迭代傾角導(dǎo)向中值濾波器的制造方法
【專利摘要】提供了用于對在地震數(shù)據(jù)中示出有沖突傾角的地震數(shù)據(jù)進行隨機噪聲衰減的迭代傾角導(dǎo)向中值濾波器。所述數(shù)據(jù)的處理窗口或樣本內(nèi)的若干支配傾角通過傅里葉徑向變換在頻率-波數(shù)域中被識別。之后沿著所述支配傾角施加中值濾波器以去除噪聲，濾波之后的剩余信號被保留用于進一步的中值濾波器迭代。重復(fù)迭代以在所述剩余數(shù)據(jù)中沿著主要支配傾角施加中值濾波器。在隨后的迭代中繼續(xù)處理直到所有選定的傾角被處理完。每次迭代的剩余信號之后被累加用于最終輸出。
【專利說明】用于地震數(shù)據(jù)處理的迭代傾角導(dǎo)向中值濾波器

【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及用于地下地層結(jié)構(gòu)分析的地震數(shù)據(jù)處理，特別是涉及用于這種處理的迭代傾角導(dǎo)向中值濾波器。

【背景技術(shù)】
[0002]在地震數(shù)據(jù)處理中，噪聲衰減對于產(chǎn)生相關(guān)區(qū)域中的地震勘測結(jié)果的更加精確的表征非常重要。地震數(shù)據(jù)的疊前道集中的信噪比的優(yōu)化處理可產(chǎn)生更佳的后續(xù)處理、成像和解釋。
[0003]在努力抑制噪聲方面已采用若干的信號處理技術(shù)。這些技術(shù)被分類成三個主要組:頻-空(f-x)域預(yù)測濾波、奇異值分解(SVD)法和中值濾波。
[0004]在若干年前，在f-x域中采取通過預(yù)測反褶積來進行隨機噪聲衰減。F-X反褶積是基于假設(shè)在每個單一頻率上的空間信號是由有限數(shù)量的復(fù)諧波所構(gòu)成的。在噪聲存在的情況下，自回歸模型適合于預(yù)測諧波的迭加。F-X反褶積對衰減隨機噪聲有效并且能夠處理被稱為沖突傾角(conflicting dip)的事件。沖突傾角指的是地震數(shù)據(jù)示出可能會出現(xiàn)多于一個可能的傾角的情形。然而，f_x反褶積已知會在存在極強噪聲時使信號電平顯著失真。由此建議采取投影濾波器來取代預(yù)測濾波器。
[0005]奇異值分解(SVD)是優(yōu)化地震道集中的橫向相關(guān)事件的另一工具。在采用特征值分解以提取相關(guān)事件前，奇異值分解形成數(shù)據(jù)協(xié)方差矩陣。SVC能通過只累加表示橫向相關(guān)信號的最大奇異值的貢獻來有效地抑制隨機噪聲。擴展的SVD應(yīng)用已用于f-χ域中的地震數(shù)據(jù)。
[0006]中值濾波同樣已經(jīng)在石油和天然氣工業(yè)中被廣泛接受。在Bednar, J.B.于1983年發(fā)表在 Geophysics 第 48 期第 1598-1610 頁上的文章“Applicat1ns of median filteringto deconvolut1n, pulse estimat1n, and statistical editing of seismic data，，中討論了中值濾波在地震勘探中的一些應(yīng)用。在Duncan, G.和G.Beresford于1995年發(fā)表在Geophysics 第 60 期第 1157-1168 頁上的文章“Some analyses of2D median f-k filters”中介紹了 2D中值f-k濾波器，其使用f一k濾波器的截短的脈沖響應(yīng)系數(shù)作為加權(quán)系數(shù)用于加權(quán)的中值處理。在Zhang, R.和T.J.Ulrych于2003年發(fā)表在SEG擴充文摘第1949-1952頁的文章“Multiple suppress1n based on the migrat1n operator and a hyperbolicmedian filter”中探討了使用雙曲線中值濾波器來進行多重抑制，而在Liu, C.、Y.Liu、B.Yang、D.Wang 和 J.Sun 于 2006 年發(fā)表在 Geophysics 第 71 期 V105-V110 上的文章 “A2Dmultistage median filter to reduce random seismic noise，，中米用 2D 多級中值濾波器來抑制陸地地震數(shù)據(jù)中的隨機噪聲。在Liu, Y.、Y.Luo和Y.Wang于2009年發(fā)表在SEG 擴充文摘 29 第 3342-3346 頁的文章 “Vector median filter and its Applicat1nsin Geophysics”中提出在地質(zhì)物理中應(yīng)用矢量中值濾波器(VMF)。在Huo，S.、Luo，Y.和P.G.Kelamis 于 2009 年發(fā)表在 SEG 擴充文摘 28 第 31-35 頁的文章“Simultaneous sourcesseparat1n via mult1-direct1nal vector-median filter” 中討論了將矢量中值濾波器擴展為多方向矢量中值濾波器(MD-VMF)以分隔混合的現(xiàn)場地震數(shù)據(jù)。
[0007]中值濾波器處理假設(shè)已經(jīng)預(yù)先對相關(guān)事件進行平整化，而MD-VMF濾波假設(shè)在操作窗口中存在單一傾角。因此，目前已知的，中值濾波器或多方向矢量中值濾波器方法都不能夠處理具有沖突傾角的地震道集。
[0008]當(dāng)?shù)卣饠?shù)據(jù)示出沖突傾角時使用被稱為頻率-波數(shù)(或F-K)濾波的技術(shù)。經(jīng)驗顯示頻率-波數(shù)濾波器不是邊緣保持濾波器，因此在濾波之后具有邊緣效應(yīng)。F-K濾波器不能令人滿意地執(zhí)行尖峰或脈沖噪聲的消除。如上所記載的，中值濾波器，作為邊緣保持濾波器，一般適用于衰減尖峰噪聲，但是只對平整的事件有用。


【發(fā)明內(nèi)容】

[0009]簡而言之，本發(fā)明提供了一種新的和改進的用于衰減其中時間-空間處理結(jié)果指示有沖突傾角的地震數(shù)據(jù)中的隨機噪聲的處理地震數(shù)據(jù)的計算機實施的方法。根據(jù)本發(fā)明的計算機實施的方法，將地震數(shù)據(jù)的時間-空間處理結(jié)果在計算機中聚集以形成該地震數(shù)據(jù)的一系列重疊的時間-空間樣本窗口。將單個地震數(shù)據(jù)的樣本窗口的數(shù)據(jù)轉(zhuǎn)換為頻率-波數(shù)域數(shù)據(jù)。隨后，將通過轉(zhuǎn)化樣本窗口數(shù)據(jù)所獲得的頻率-波數(shù)域數(shù)據(jù)轉(zhuǎn)換為傅里葉徑向域數(shù)據(jù)。確定樣本窗口的傅里葉徑向域數(shù)據(jù)的峰值圖，對選定的傅里葉徑向域數(shù)據(jù)的峰值進行轉(zhuǎn)換，以指示出樣本窗口數(shù)據(jù)中的支配傾角。沿著指示出的支配傾角中的一個所選定的支配傾角將中值濾波器施加到地震數(shù)據(jù)的樣本窗口中的數(shù)據(jù)，從而對該樣本窗口中的數(shù)據(jù)進行噪聲衰減。對中值濾波后的數(shù)據(jù)的剩余信號進行存儲。
[0010]本發(fā)明還提供了一種新的和改進的用于衰減其中時間-空間處理結(jié)果指示有沖突傾角的地震數(shù)據(jù)中的隨機噪聲的地震數(shù)據(jù)的數(shù)據(jù)處理系統(tǒng)。該數(shù)據(jù)處理系統(tǒng)包括存儲地震數(shù)據(jù)的時間-空間處理結(jié)果的數(shù)據(jù)存儲器，以及處理器。該數(shù)據(jù)處理系統(tǒng)的處理器將地震數(shù)據(jù)的時間-空間處理結(jié)果在計算機中聚集以形成該地震數(shù)據(jù)的一系列重疊的時間-空間樣本窗口，將單個地震數(shù)據(jù)的樣本窗口的數(shù)據(jù)轉(zhuǎn)換為頻率-波數(shù)域數(shù)據(jù)。隨后，處理器將通過轉(zhuǎn)化樣本窗口數(shù)據(jù)所獲得的頻率-波數(shù)域數(shù)據(jù)轉(zhuǎn)換為傅里葉徑向域數(shù)據(jù)，并確定樣本窗口的傅里葉徑向域數(shù)據(jù)的峰值圖。隨后，處理器根據(jù)樣本窗口數(shù)據(jù)中指示出的支配傾角對選定的傅里葉徑向域數(shù)據(jù)的峰值進行轉(zhuǎn)換。處理器沿著所指示出的支配傾角的一個選定的支配傾角將中值濾波器施加到樣本窗口的數(shù)據(jù)，從而對該樣本窗口中的數(shù)據(jù)進行噪聲衰減。處理器然后對樣本窗口的中值濾波后的數(shù)據(jù)的剩余信號進行存儲。
[0011]本發(fā)明還提供了一種新的和改進的在計算機可讀介質(zhì)中存儲計算機可操作指令的數(shù)據(jù)存儲裝置，所述計算機可操作指令使數(shù)據(jù)處理系統(tǒng)對其中時間-空間處理結(jié)果指示有沖突傾角的地震數(shù)據(jù)進行處理以衰減該地震數(shù)據(jù)中的隨機噪聲。存儲在數(shù)據(jù)存儲裝置中的指令引發(fā)數(shù)據(jù)處理系統(tǒng)將地震數(shù)據(jù)的時間-空間處理結(jié)果在計算機中聚集以形成該地震數(shù)據(jù)的一系列重疊的時間-空間樣本窗口。該指令還引發(fā)數(shù)據(jù)處理系統(tǒng)將單個地震數(shù)據(jù)的樣本窗口的數(shù)據(jù)轉(zhuǎn)換為頻率-波數(shù)域數(shù)據(jù)，并將通過轉(zhuǎn)化樣本窗口數(shù)據(jù)所獲得的頻率-波數(shù)域數(shù)據(jù)轉(zhuǎn)換為傅里葉徑向域數(shù)據(jù)。該指令還引發(fā)數(shù)據(jù)處理系統(tǒng)確定樣本窗口的傅里葉徑向域數(shù)據(jù)的峰值圖，并根據(jù)樣本窗口數(shù)據(jù)中指示出的支配傾角對選定的傅里葉徑向域數(shù)據(jù)的峰值進行轉(zhuǎn)換。該指令引發(fā)數(shù)據(jù)處理系統(tǒng)沿著所指示出的支配傾角中的一個選定的支配傾角將中值濾波器施加到地震數(shù)據(jù)的樣本窗口中的數(shù)據(jù)，從而對樣本窗口中的數(shù)據(jù)進行噪聲衰減。隨后，該指令引發(fā)數(shù)據(jù)處理系統(tǒng)對樣本窗口的中值濾波后數(shù)據(jù)的剩余信號進行存儲。

【專利附圖】

【附圖說明】
[0012]圖1是根據(jù)本發(fā)明的在數(shù)據(jù)處理系統(tǒng)中執(zhí)行的用于迭代傾角導(dǎo)向中值濾波器的數(shù)據(jù)處理步驟組的功能框圖。
[0013]圖2是根據(jù)本發(fā)明的用于迭代傾角導(dǎo)向中值濾波器的數(shù)據(jù)處理系統(tǒng)的示意框圖。
[0014]圖3A、圖3B和圖3C是在不同域中的樣本平面波的變換示意圖。
[0015]圖4A、圖4B和圖4C是處理后的地震數(shù)據(jù)的繪圖組。
[0016]圖5A、圖5B、圖5C、圖和圖5E是處理后的地震數(shù)據(jù)的另一繪圖組。
[0017]圖6A、圖6B和圖6C是處理后的地震數(shù)據(jù)的又一繪圖組。
[0018]圖7A、圖7B、圖7C、圖7D和圖7E是處理后的地震數(shù)據(jù)的追加繪圖組。

【具體實施方式】
[0019]本發(fā)明提供了在用于相關(guān)地下結(jié)構(gòu)分析的地震數(shù)據(jù)中衰減隨機噪聲的迭代傾角導(dǎo)向濾波器。如將要闡述的，在對數(shù)據(jù)執(zhí)行傅里葉徑向變換之后通過在數(shù)據(jù)窗口中搜尋多個支配傾角來處理數(shù)據(jù)。之后以降序序列施加中值濾波。通過這種方法，可以消除具有沖突傾角的數(shù)據(jù)中的隨機噪聲。
[0020]由一組數(shù)據(jù)處理步驟所組成的流程圖F (圖1)示出了在計算機程序軟件中實現(xiàn)的本發(fā)明的邏輯結(jié)構(gòu)。流程圖F是高級邏輯流程圖，其示出了根據(jù)本發(fā)明的通過根據(jù)本發(fā)明的迭代傾角導(dǎo)向中值濾波來處理地震數(shù)據(jù)的方法。應(yīng)當(dāng)理解的是，該流程圖示出了根據(jù)本發(fā)明運行的計算機程序碼元的結(jié)構(gòu)。本發(fā)明在其必要實施例中通過計算機組件來實行，所述計算機組件以這樣一種形式使用程序代碼指令，即，指示數(shù)字數(shù)據(jù)處理系統(tǒng)D(圖2)執(zhí)行對應(yīng)于流程圖F中所示步驟的處理步驟序列。圖1的流程圖F包括針對用于隨機噪聲衰減的在數(shù)據(jù)處理系統(tǒng)D中迭代傾角導(dǎo)向中值濾波地震數(shù)據(jù)的計算機實施處理數(shù)據(jù)的優(yōu)選步驟序列。
[0021]流程圖F是高級邏輯流程圖，其示出了根據(jù)本發(fā)明的處理方法。本發(fā)明的方法在數(shù)據(jù)處理系統(tǒng)D的計算機10(圖2)中執(zhí)行并且可以使用圖1的存儲在存儲器12中且由計算機10的系統(tǒng)處理器14執(zhí)行的計算機程序步驟實施。如將要闡述的，流程圖F示出了采用數(shù)據(jù)處理系統(tǒng)D來迭代傾角導(dǎo)向中值濾波地震數(shù)據(jù)以衰減地震數(shù)據(jù)中的不需要的隨機噪聲及用于其他目的的優(yōu)選實施例。
[0022]在流程圖F的步驟20(圖1)期間，數(shù)據(jù)處理系統(tǒng)D的存儲器中的輸入地震數(shù)據(jù)被分隔并聚集到若干個時間-空間重疊的小數(shù)據(jù)樣本窗口中。應(yīng)當(dāng)理解的是，用于根據(jù)本發(fā)明進行處理的輸入地震數(shù)據(jù)可以是疊前道集或疊后剖面。輸入數(shù)據(jù)被聚集到的作為開始的第一個時間-空間窗口以及相繼的其他數(shù)據(jù)窗口中，處理如圖1所示的后續(xù)序列繼續(xù)進行。
[0023]在步驟22期間，將被處理的地震數(shù)據(jù)的選定時間-空間樣本窗口從時間-空間域變換為頻率-波數(shù)(f-k)域。在步驟24期間，對樣本窗口數(shù)據(jù)執(zhí)行步驟22所產(chǎn)生的f-k域數(shù)據(jù)被變換為傅里葉徑向域，優(yōu)選地使用公知的環(huán)繞徑向道變換。
[0024]在步驟26期間，通過累加從步驟24產(chǎn)生的傅里葉徑向變換數(shù)據(jù)中所示出的每個徑向形成峰值圖。在步驟28中，評估在步驟26期間所形成的傅里葉徑向域峰值圖中的選定峰值以示出用于樣本窗口的原始t-x域數(shù)據(jù)中的支配傾角。按照峰值圖中示出的它們各自的峰值來評估傾角的支配度。
[0025]在步驟30期間，之后沿著在步驟28期間示出的支配傾角的第一選定支配傾角的角度施加中值濾波器以衰減來自信號的噪聲。以在峰值圖中示出的不同傾角的各自的峰值的降序序列以不同的示出傾角迭代地執(zhí)行中值濾波步驟30。由于通常存在多個支配傾角，本發(fā)明使用以示出峰值的降序序列按照傾角進行濾波的迭代方法。對于中值濾波的每個處理迭代，作為步驟30的一部分所施加的中值濾波器被沿著在剩余傾角中具有最大示出峰值的剩余傾角的對應(yīng)角度施加。處理針對每個選定傾角繼續(xù)進行。在步驟30期間，在每個迭代之后，將前面預(yù)估的信號數(shù)據(jù)從輸入數(shù)據(jù)中減去，剩余或殘留的數(shù)據(jù)被用作在步驟30期間的后續(xù)中值濾波的樣本窗口數(shù)據(jù)的輸入數(shù)據(jù)。殘留數(shù)據(jù)同樣存儲在計算機存儲器中。當(dāng)沿著由峰值圖示出的每個選定支配傾角的中值濾波結(jié)束時，剩余的噪聲衰減后的信號被累加以形成用于樣本數(shù)據(jù)窗口的最終輸出值。結(jié)果是當(dāng)前被處理的樣本窗口的噪聲衰減后的信號的傾角導(dǎo)向中值濾波最終輸出。針對樣本窗口的在中值濾波步驟30期間的迭代數(shù)量取決于數(shù)據(jù)的復(fù)雜性(即，存在于樣本窗口的數(shù)據(jù)中的示出沖突傾角的數(shù)量)。
[0026]之后在步驟32期間，將步驟30產(chǎn)生的針對當(dāng)前被處理的樣本窗口的中值濾波后的數(shù)據(jù)存儲到數(shù)據(jù)處理系統(tǒng)D的存儲器中。在步驟34期間，確認輸入數(shù)據(jù)中的每個樣本窗口是否都已經(jīng)被處理。如果否，則在步驟36期間增加數(shù)據(jù)窗口序列或計數(shù)到新的相關(guān)數(shù)據(jù)窗口并且處理回歸到步驟22以繼續(xù)以上述方式進行處理。針對連續(xù)樣本窗口的處理數(shù)據(jù)的序列以這樣的方式繼續(xù)直到相關(guān)的輸入時間-空間數(shù)據(jù)的每個數(shù)據(jù)樣本窗口被處理。
[0027]當(dāng)在步驟34期間確定該處理完成時，步驟38使得每個時間-空間數(shù)據(jù)樣本窗口的信號輸出被從數(shù)據(jù)處理系統(tǒng)D的存儲器中聚合。步驟38通過顯示器26產(chǎn)生傾角導(dǎo)向中值濾波后的噪聲衰減的地震數(shù)據(jù)的輸出地震顯示用于分析、解釋和評估。
[0028]基于下面的分析，可進一步理解圖1中所示出及上面所描述的處理序列。在t-x域中平面波d(x，t)可以表達為
[0029]

【權(quán)利要求】
1.一種用于衰減其中時間-空間處理結(jié)果指示有沖突傾角的地震數(shù)據(jù)中的隨機噪聲的處理地震數(shù)據(jù)的計算機實施的方法，該方法包括該計算機實施的處理步驟: (a)將地震數(shù)據(jù)的時間-空間處理結(jié)果在計算機中聚集以形成該地震數(shù)據(jù)的一系列重疊的時間-空間樣本窗口； (b)將單個地震數(shù)據(jù)的樣本窗口的數(shù)據(jù)轉(zhuǎn)換為頻率-波數(shù)域數(shù)據(jù)； (C)將樣本窗口數(shù)據(jù)的頻率-波數(shù)域數(shù)據(jù)轉(zhuǎn)換為傅里葉徑向域數(shù)據(jù)； (d)確定樣本窗口的傅里葉徑向域數(shù)據(jù)的峰值圖； (e)對選定的樣本窗口的傅里葉徑向域數(shù)據(jù)的峰值進行轉(zhuǎn)換，以指示出樣本窗口數(shù)據(jù)中的支配傾角； (f)沿著選定的一個所指示的支配傾角將中值濾波器施加到樣本窗口中的數(shù)據(jù)，從而對該樣本窗口中的數(shù)據(jù)進行噪聲衰減；以及 (g)對樣本窗口的中值濾波后的數(shù)據(jù)中的剩余信號進行存儲。
2.如權(quán)利要求1所述的計算機實施的方法，還包括步驟: 確定所指示的樣本窗口數(shù)據(jù)中的支配傾角中的每一個是否已被選擇；以及如果否，選擇另一個指示出的支配傾角并沿著該選定的另一個指示出的支配傾角重復(fù)步驟(f)和(g);或 如果是，形成沿著選定的多個支配傾角輸出的所存儲的中值濾波后的數(shù)據(jù)的剩余信號的累加和。
3.如權(quán)利要求1所述的計算機實施的方法，還包括步驟: 確定每個所聚集的地震數(shù)據(jù)的樣本窗口是否已被選擇；以及 如果否，選擇另一個所聚集的地震數(shù)據(jù)的樣本窗口，并返回到步驟(b)對該另一個所聚集的地震數(shù)據(jù)的樣本窗口的數(shù)據(jù)執(zhí)行步驟(b)到(h);或 如果是，對聚集的所確定的針對每個重疊的樣本窗口輸出的噪聲衰減后的數(shù)據(jù)的累加和進行存儲。
4.如權(quán)利要求2所述的計算機實施的方法，還包括步驟: 顯示所存儲的聚集的所確定的針對每個重疊的樣本窗口輸出的中值濾波后的數(shù)據(jù)的累加和。
5.如權(quán)利要求1所述的計算機實施的方法，其中所處理的地震數(shù)據(jù)包括疊前道集。
6.如權(quán)利要求1所述的計算機實施的方法，其中所處理的地震數(shù)據(jù)包括疊后剖面。
7.如權(quán)利要求1所述的計算機實施的方法，其中所述確定傅里葉徑向域數(shù)據(jù)的峰值圖的步驟包括步驟: 確定傅里葉徑向域數(shù)據(jù)中每個徑向的累加和。
8.如權(quán)利要求7所述的計算機實施的方法，其中所述施加中值濾波器的步驟包括步驟: 以傾角降序順序?qū)Ω鲀A角實施中值濾波器。
9.一種用于衰減其中時間-空間處理結(jié)果指示有沖突傾角的地震數(shù)據(jù)中的隨機噪聲的地震數(shù)據(jù)的數(shù)據(jù)處理系統(tǒng)，該數(shù)據(jù)處理系統(tǒng)包括: 存儲地震數(shù)據(jù)的時間-空間處理結(jié)果的數(shù)據(jù)存儲器； 執(zhí)行以下步驟的處理器:(a)將地震數(shù)據(jù)的時間-空間處理結(jié)果在計算機中聚集以形成該地震數(shù)據(jù)的一系列重疊的時間-空間樣本窗口； (b)將單個地震數(shù)據(jù)的樣本窗口的數(shù)據(jù)轉(zhuǎn)換為頻率-波數(shù)域數(shù)據(jù)； (C)將樣本窗口數(shù)據(jù)的頻率-波數(shù)域數(shù)據(jù)轉(zhuǎn)換為傅里葉徑向域數(shù)據(jù)； (d)確定樣本窗口的傅里葉徑向域數(shù)據(jù)的峰值圖； (e)對選定的樣本窗口的傅里葉徑向域數(shù)據(jù)的峰值進行轉(zhuǎn)換，以指示出樣本窗口數(shù)據(jù)中的支配傾角； (f)沿著所指示的支配傾角將中值濾波器施加到樣本窗口中的數(shù)據(jù),從而對該樣本窗口中的數(shù)據(jù)進行噪聲衰減；以及 (g)對樣本窗口的中值濾波后的數(shù)據(jù)輸出的剩余信號進行存儲。
10.如權(quán)利要求9所述的數(shù)據(jù)處理系統(tǒng)，其中處理器還執(zhí)行步驟: 確定所指示的樣本窗口數(shù)據(jù)中的支配傾角中的每一個是否已被選擇；以及 如果否，選擇另一個指示出的支配傾角并沿著該選定的另一個指示出的支配傾角重復(fù)步驟(f)和(g);或 如果是，形成所存儲的沿著選定的多個支配傾角輸出的中值濾波后的數(shù)據(jù)的剩余信號的累加和。
11.如權(quán)利要求9所述的數(shù)據(jù)處理系統(tǒng)，其中處理器還執(zhí)行步驟: 確定每個所聚集的地震數(shù)據(jù)的樣本窗口的是否已被選擇；以及 如果否，選擇另一個所聚集的地震數(shù)據(jù)的樣本窗口，并返回到步驟(b)從而對該另一個所聚集的地震數(shù)據(jù)的樣本窗口的數(shù)據(jù)執(zhí)行步驟(b)到(h);或 如果是，對聚集的所確定的針對每個重疊的樣本窗口輸出的中值濾波后的數(shù)據(jù)的累加和進行存儲。
12.如權(quán)利要求11所述的數(shù)據(jù)處理系統(tǒng)，還包括: 顯示器，其對所存儲的聚集的所確定的針對每個重疊的樣本窗口輸出的中值濾波后的數(shù)據(jù)的累加和形成輸出記錄。
13.如權(quán)利要求11所述的數(shù)據(jù)處理系統(tǒng)，其中處理器在執(zhí)行所述確定傅里葉徑向域數(shù)據(jù)的峰值圖的步驟中執(zhí)行步驟: 確定傅里葉徑向域數(shù)據(jù)中的每個徑向的累加和。
14.如權(quán)利要求10所述的數(shù)據(jù)處理系統(tǒng)，其中處理器在執(zhí)行所述施加中值濾波器的步驟中執(zhí)行步驟: 以傾角降序順序?qū)Ω鲀A角實施中值濾波。
15.在計算機可讀介質(zhì)中存儲計算機可操作指令的數(shù)據(jù)存儲裝置，所述計算機可操作指令使數(shù)據(jù)處理系統(tǒng)對其中時間-空間處理結(jié)果指示有沖突傾角的地震數(shù)據(jù)進行處理以衰減該地震數(shù)據(jù)中的隨機噪聲，存儲在數(shù)據(jù)存儲裝置中的指令引發(fā)數(shù)據(jù)處理系統(tǒng)執(zhí)行下列步驟: (a)將地震數(shù)據(jù)的時間-空間處理結(jié)果在計算機中聚集以形成該地震數(shù)據(jù)的一系列重疊的時間-空間樣本窗口； (b)將單個地震數(shù)據(jù)的樣本窗口的數(shù)據(jù)轉(zhuǎn)換為頻率-波數(shù)域數(shù)據(jù)； (c)將樣本窗口數(shù)據(jù)的頻率-波數(shù)域數(shù)據(jù)轉(zhuǎn)換為傅里葉徑向域數(shù)據(jù)；(d)確定樣本窗口的傅里葉徑向域數(shù)據(jù)的峰值圖； (e)對選定的樣本窗口的傅里葉徑向域數(shù)據(jù)的峰值進行轉(zhuǎn)換，以指示出樣本窗口數(shù)據(jù)中的支配傾角； (f)沿著所指示的支配傾角將中值濾波器施加到樣本窗口中的數(shù)據(jù),從而對該樣本窗口中的數(shù)據(jù)進行噪聲衰減；以及 (g)對中值濾波后的數(shù)據(jù)輸出的剩余信號進行存儲。
16.如權(quán)利要求15所述的數(shù)據(jù)存儲裝置，其中指令還包括引發(fā)數(shù)據(jù)處理系統(tǒng)執(zhí)行如下步驟的指令: 確定所指示的樣本窗口數(shù)據(jù)中的支配傾角中的每一個是否已被選擇；以及 如果否，選擇另一個指示出的支配傾角并沿著該選定的另一個指示出的支配傾角重復(fù)步驟(f)和(g);或 如果是，形成所存儲的沿著選定的多個支配傾角輸出的中值濾波后的數(shù)據(jù)的剩余信號的累加和。
17.如權(quán)利要求15所述的數(shù)據(jù)存儲裝置，其中指令還包括引發(fā)數(shù)據(jù)處理系統(tǒng)執(zhí)行如下步驟的指令: 確定每個所聚集的地震數(shù)據(jù)的樣本窗口的是否已被選擇；以及 如果否，選擇另一個所聚集的地震數(shù)據(jù)的樣本窗口，并返回到步驟(b)從而對該另一個所聚集的地震數(shù)據(jù)的樣本窗口的數(shù)據(jù)執(zhí)行步驟(b)到(h);或 如果是，對聚集的所確定的針對每個重疊的樣本窗口輸出的中值濾波后的數(shù)據(jù)的累加和進行存儲。
18.如權(quán)利要求15所述的數(shù)據(jù)存儲裝置，其中指令還包括引發(fā)數(shù)據(jù)處理系統(tǒng)執(zhí)行如下步驟的指令: 對所存儲的聚集的所確定的針對每個重疊的樣本窗口輸出的中值濾波后的數(shù)據(jù)的累加和進行顯示。
19.如權(quán)利要求15所述的數(shù)據(jù)存儲裝置，其中指令還包括引發(fā)數(shù)據(jù)處理系統(tǒng)在確定傅里葉徑向域數(shù)據(jù)的峰值圖時執(zhí)行如下步驟的指令: 確定傅里葉徑向域數(shù)據(jù)中的每個徑向的累加和。
20.如權(quán)利要求15所述的數(shù)據(jù)存儲裝置，其中指令還包括引發(fā)數(shù)據(jù)處理系統(tǒng)在施加中值濾波器時執(zhí)行如下步驟的指令: 以傾角降序順序?qū)Ω鲀A角實施中值濾波。
【文檔編號】G01V1/36GK104081226SQ201280062160
【公開日】2014年10月1日 申請日期:2012年12月14日 優(yōu)先權(quán)日:2011年12月15日
【發(fā)明者】S·霍, W·朱 申請人:沙特阿拉伯石油公司