專利名稱：壓低系統(tǒng)性振幅變化的方法
技術(shù)領(lǐng)域：
本發(fā)明涉及地球物理勘探領(lǐng)域，特別是涉及處理地震數(shù)據(jù)的方法。
背景技術(shù)：
對地震數(shù)據(jù)的系統(tǒng)性污染，例如源于數(shù)據(jù)獲取努力欠佳的系統(tǒng)性污染，能把自己表現(xiàn)為信號振幅的變化和表現(xiàn)為殘余噪聲振幅的變化。這兩種形式的污染橫向上與它們的來源相關(guān)，最普通的來源是獲取系統(tǒng)的幾何布局。對于放炮幾何布局(shooting geometry)造成的人為結(jié)果，它們被稱作獲取腳印(acquisition footprint)或印記。當前的算法使用各種技術(shù)來減小獲取腳印。在Morse，P.F.和Hildebramdt，G.F.(1989)的“借助疊加陣列壓低地滾波”(Geophysics，54，No.3，pp290-301)中公開的一種途徑通過校正振幅來解決信號變化問題，例如在地震處理中創(chuàng)建一個均一的疊加算子。然而，如Hampson(1994)在“波場樣本和相干噪聲衰減之間的關(guān)系”(歐洲地球物理學家和工程師協(xié)會第56屆年會，摘要集，文章H054)中指出的那樣，處理過程中形成的陣列能降低分辨力，所以更好的辦法是以適當?shù)姆绞将@取數(shù)據(jù)，但這增大獲取成本。Ronen(1994)在“處置不規(guī)則幾何布局等值化DMO及其他”(勘探地球物理學會第64屆年會擴展摘要，pp1545-1548)中的解決方案是使用特殊的DMO算法均勻地振幅變化。
Meunier，J.和Belissent，M.(1992)在“降低三維幾何分布產(chǎn)生的人為結(jié)果”(第6屆委內(nèi)瑞拉地球物理大會摘要，pp388-394)中提出的途徑依賴于污染的周期性。這一組中的算法在更適當?shù)挠蛑胁僮?，如波?shù)(K)域或最通用的頻率—波數(shù)(FK)域，例如，Necati Gulunay的“去掉獲取幾何部局腳印”(勘探地球物理學會第69屆年會擴展摘要，pp637-640)。然而，在FK域中操作的算法有一個缺點，即去掉印記只對有單一的和不變的取向或傾斜的事件才有效。

發(fā)明內(nèi)容
處理地震數(shù)據(jù)以從地震數(shù)據(jù)中去掉畸變的一種方法包含對分解后的地震數(shù)據(jù)變換分量應(yīng)用一個酉變換，以得到修正的變換分量。修正的變換分量被濾波，濾波后的數(shù)據(jù)被反演以得到時間—距離(T-X)或深度—距離(Z-X)地震數(shù)據(jù)。


通過參考下文的詳細描述和附圖，本發(fā)明及其優(yōu)點將被更好地理解。
附圖中圖1A顯示時間—距離域中的示意地震剖面；圖1B顯示圖1A數(shù)據(jù)到F-K域的付立葉變換；圖2A顯示以欠佳的獲取參數(shù)獲取的時間—距離域示意地震剖面；圖2B顯示圖2A數(shù)據(jù)到F-K域的付立葉變換；圖3顯示本發(fā)明處理技術(shù)的流程圖；圖4顯示根據(jù)圖3流程圖的數(shù)據(jù)流示意圖；圖5顯示以無畸變地震數(shù)據(jù)產(chǎn)生的合成地震剖面；圖6顯示以被畸變的地震數(shù)據(jù)產(chǎn)生的合成地震剖面；圖7顯示以標準的先有技術(shù)應(yīng)用對圖6的地震數(shù)據(jù)濾波后產(chǎn)生的合成地震剖面；圖8顯示對圖6的地震數(shù)據(jù)應(yīng)用本發(fā)明后產(chǎn)生的合成地震剖面。
盡管將結(jié)合優(yōu)選實施例描述本發(fā)明，但就該理解，本發(fā)明不限于此。相反，它要覆蓋可能落入由所附權(quán)利要求定義的本發(fā)明精神和范圍的所有替代物、修改和等效物。
具體實施例方式
本發(fā)明是一種壓低可能由獲取特性(Signature)造成的系統(tǒng)性振幅變化或其他人為信號的方法，它克服了當前方法的上述缺陷?；谙挛牡脑敿毭枋?，本發(fā)明的其他優(yōu)點對于本領(lǐng)域技術(shù)人員將會顯而易見。至于下文的詳細描述是針對本發(fā)明的特定實施例或特定用途，但這只是作為舉例說明，不能被認為是對本發(fā)明范圍的限定。
本發(fā)明的目的是提供一個不依賴于地震數(shù)據(jù)中的結(jié)構(gòu)信息的解決方案，即該算法不應(yīng)需要遭侵害的反射的先驗信息。該方法通過使用反射事件的重復(fù)性來解決它們的振幅變化問題，并以不依賴于傾角的方式完成這件事。然而，該方法對于去掉殘余噪聲印記不那么有效。
該方法對相位變化也不那么有效。如果可能的話，相位變化應(yīng)以確定性的或自適應(yīng)的方式來去掉。例如，在海上拖纜深度變化的情況中，能使用沿拖纜各位置處記錄到的實際深度確定性地去掉相應(yīng)的相位差。
使用由一個簡單模型得到的合成數(shù)據(jù)能適當?shù)孛枋鲈摲椒?。將在過程的各個階段展現(xiàn)這些合成數(shù)據(jù)，以解釋和以圖說明本發(fā)明的方法。
該優(yōu)選實施例使用付立葉變換，它是具有對這一方法應(yīng)用有利的一些性質(zhì)的酉變換。然而，也可使用許多其他的酉變換。酉變換是形如A’＝UAU*的變換，這里U*代表代表一個共軛算子(即，它是一個矩陣的轉(zhuǎn)置，該矩陣中的每個元素被其余子式所代替)。其他酉變換的例子在本領(lǐng)域是公知的。本發(fā)明的另一個實施例不是使用付立葉變換把地震數(shù)據(jù)變換到變換域，而是使用Haar變換。本發(fā)明的又一個實施例使用Haademard變換。這些變換是本領(lǐng)域有經(jīng)驗的人們都會知道的酉變換。本發(fā)明的又一個實施例使用小波變換。小波變換與Haar變換有同樣的結(jié)構(gòu)，但如YvesMeyer說明的那樣，小波變換還有一個額外的性質(zhì)，即相應(yīng)的展開是可以逐項微分的，如果地震波形道是光滑的話。利用這些酉變換中的任何一個，可以有效地進行變換、濾波和壓低不希望的人為信號。
圖1A顯示時間-距離域(T-X域)中的示意地震剖面，代表穿過地球的垂直切片。為簡化該示意圖示，未畫出地震波形。一個地震剖面包含大量依賴于時間(T)的波形道，被一個挨一個地沿水平(X)方向排列。該地震剖面含有來自一個水平地震反射體7的響應(yīng)、來自傾斜反射體8的響應(yīng)以及來自地球表面的散射障礙體9的響應(yīng)。
為了說明，數(shù)據(jù)首先從T-X域變換到F-K域，例如通過使用付立葉變換算法。這些數(shù)據(jù)付立葉變換到F-K域產(chǎn)生的結(jié)果如圖1B中所示。這里，來自水平反射表面的能量被垂直投射到K＝0軸10，而來自傾斜表面的能量被投射到與K＝0軸成一夾角的線11上。來自地震散射體的能量投射為與K＝0軸成另一夾角的線12。這種情況說明把數(shù)據(jù)變換到F-K域的好處之一；盡管在T-X域中散射體干擾所希望的反射，但在F-K域中，如果能避免假頻的話，則存在完全的分離。不希望的能量能從F-K域中切除掉，然后能執(zhí)行逆變換，從而恢復(fù)沒有干擾能量的原始數(shù)據(jù)。
現(xiàn)在讓我們假定使用減少的獲取努力產(chǎn)生欠佳數(shù)據(jù)集，由此激發(fā)出這同一個剖面。于是，其響應(yīng)看來如同圖2A中所示剖面。這些振幅被以系統(tǒng)方式5橫向調(diào)制，如同在數(shù)據(jù)獲取過程中被例如太粗糙的“逐點前進步長(roll along step)”所致。如本領(lǐng)域有經(jīng)驗的人們所知，“逐點前進步長”是指從一個震源位置到下一個震源位置的距離。在這個示意性舉例中，變化被簡化以顯示一個矩形類調(diào)制。在變換到F-K域之后的數(shù)據(jù)示于圖2B中。在這一情況中，三個事件的干擾存在于F-K域中，甚至比在原來的時間—距離域中還要差。具體地說，現(xiàn)在水平事件不僅將其能量投射到主線10上，而且還投射到若干個重復(fù)線上，圖中只顯示其中的兩個重復(fù)20和30。傾斜反射體將其能量投射到線11、21和31上，而散射體映射到線12、22和32上。對于傾斜事件和散射，還是存在更多的再發(fā)生，但只顯示兩個重復(fù)。當把圖2B與圖1B比較時，能得出如下結(jié)論污染已使再發(fā)生事件20、21、22、30、31和32與投射到線10、11和12上的未被污染的主能量交叉。在那些交叉點上，屬于污染的重復(fù)能量不能與所希望的能量分開。在實踐中，交叉能位于F-K域中的任何地方，因為地震反射體能有從-90度到+90度范圍的所有可能傾斜。所以，在標準的F-K域中，如果不是同時也去掉屬于所希望事件的能量的話，則去掉屬于散射體12、22和32的能量是不現(xiàn)實的。要孤立和去掉屬于系統(tǒng)性污染的能量，即在波數(shù)域(沿K方向)重復(fù)的所有能量，也是不可能的。這是因為重復(fù)的能量干擾主能量10、11和12。
本發(fā)明利用標準的快速付立葉變換以及其他方法將數(shù)據(jù)從時間域變換到頻率域和從距離域變換到波數(shù)域。本發(fā)明還對復(fù)付立葉變換后的數(shù)據(jù)采用數(shù)據(jù)分解和重構(gòu)技術(shù)。
本發(fā)明的前向過程稱作前向修正的F-K變換，而反向過程稱作反向修正的F-K變換。修正的F-K變換過程使能孤立與系統(tǒng)性污染關(guān)聯(lián)的能量。
本發(fā)明處理技術(shù)的優(yōu)選實施例使用修正的前向F-K變換、壓低系統(tǒng)性污染和反向修正的F-K變換將數(shù)據(jù)轉(zhuǎn)換回到原來的T-X域。該過程示于圖3和圖4。圖3中的流程圖顯示一套處理步驟，而圖4補充顯示通過這些過程的數(shù)據(jù)流。圖4的數(shù)據(jù)板，即示意性顯示的輸入板39和輸出板60，與圖2a中顯示的響應(yīng)關(guān)聯(lián)。
圖3和圖4中顯示的前向修正的F-K變換首先使用垂直酉變換，即一維前向付立葉變換或垂直FFT40，以把輸入數(shù)據(jù)從T-X域轉(zhuǎn)換到F-K域。這一變換也稱作垂直前向變換，因為它沿著時間軸起作用，而時間軸在標準的地震剖面上是垂直顯示的。在前向付立葉變換中，數(shù)據(jù)被從實數(shù)轉(zhuǎn)換成復(fù)數(shù)，在圖4中顯示為具有實部41和虛部42。使用下列公式將復(fù)數(shù)據(jù)分解(43)為振幅和相位分量A=Re2+Im2-----(1)]]>φ=tan-1(ImRe)----(2)]]>在公式(1)和(2)中，‘A’是振幅，‘φ’是相位，‘Re’是復(fù)數(shù)據(jù)的實部，‘Im’是復(fù)數(shù)據(jù)的虛部。以與付立葉變換相似的方式，其他酉變換有與特定變換的基函數(shù)對應(yīng)的分解。
參考圖4，振幅分量44是輸入波形道的振幅譜。這一振幅譜含有不希望的污染，所以被用于進一步變換。相位分量45通常沒有受到系統(tǒng)性污染的影響或只受輕微影響，被存儲供其后在反演序列中用于數(shù)據(jù)重構(gòu)。
該過程中的下一階段是振幅分量從F-X域到F-K域的正向付立葉變換(46)。因為這一酉變換是沿X方向進行的，故稱作水平變換，在該優(yōu)選實施例中還稱作水平前向付立葉變換。由于被剝掉了相位屬性，這些數(shù)據(jù)失掉了它們的傾斜關(guān)聯(lián)(dip association)，從而使所有事件成為虛擬水平。然而，該過程是完全可逆的，即振幅數(shù)據(jù)能與保存的相位數(shù)據(jù)45組合以恢復(fù)原始輸入數(shù)據(jù)。
現(xiàn)在，圖4中所示地震剖面已被變換到修正的F-K域，并再次包括實部47和虛部48。在實部47和虛部48二者中，所希望的響應(yīng)49與源自污染的不希望的重復(fù)50分開了。因為所有傾斜事件現(xiàn)在的表現(xiàn)如同它們是水平的，所以它們?nèi)纪渡涞終＝0軸上，而且所有污染能量完全分離地投射到具有離散K值的平行線50上。
當適當獲取的數(shù)據(jù)的地動振幅譜(即F-X變換)在X方向沒有顯著變化時，在本發(fā)明中描述的過程工作得很好。例如，當波形道含有許多反射隨時間T隨機分布時就能滿足這一條件。于是，相應(yīng)的波形道頻譜是寬和平坦的(也稱作白譜)。假定記錄的時間窗足夠長。作為一級近似，地動譜是白的。然而，一般地說，即使其反射的分布不那么理想，橫向振幅連續(xù)這個條件也足以能夠得到滿足。
當處理系統(tǒng)性污染時，這構(gòu)成F-X空間中的主要橫向變化。在獲取印記的情況中，作為現(xiàn)場逐點前進(roll along)過程的結(jié)果，振幅周期性變化。所以，相關(guān)的能量將集中在在修正的F-K域中的常數(shù)基本波數(shù)及其諧波值50。相關(guān)的關(guān)鍵條件是在F-X域中數(shù)據(jù)的相位值被置為零。這還造成包括系統(tǒng)性污染在內(nèi)的所有事件的能量對稱地位于K＝0軸附近。
在這一階段，能得出如下結(jié)論修正的變換方法和標準的F-K變換方法有兩個基本差別。首先，在應(yīng)用前向修正的變換得到修正變換過的地震數(shù)據(jù)之后，所有事件被剝掉了與它們關(guān)聯(lián)的傾斜信息；所有事件變?yōu)樘摂M水平的，而且它們的能量投射到K＝0軸上。其次，在前向修正的F-K變換之后，污染投射到對稱位于K＝0附近的常數(shù)K值處。
作為前向修正的F-K變換的結(jié)果，被完全孤立的污染與所希望的能量分開投射，并能簡單地通過刪除所有相應(yīng)的K值(51)來去掉。壓低步驟51產(chǎn)生處理后的復(fù)數(shù)據(jù)，有實部52和虛部53，其能量主要在K＝0軸上。
其后的處理序列使數(shù)據(jù)通過一個反向修正的F-K變換。反向修正的F-K變換使用從F-K空間到F-X空間的逆付立葉變換(IFT)(54)，產(chǎn)生處理后的振幅分量55。反向修正的F-K變換進一步使用重構(gòu)過程56將保存的相位數(shù)據(jù)45與處理過的振幅數(shù)據(jù)55組合。造成的復(fù)頻譜包括一個實部57和一個虛部58。重構(gòu)過程(‘合成’步驟)的公式是Re＝Acos(φ) (3)Im＝Asin(φ) (4)在反向修正的F-K變換中的最后步驟涉及從重構(gòu)的F-X譜(即反向合成的數(shù)據(jù))到原來T-X域60的逆付立葉變換59。過程的結(jié)果60是存在于輸入39上的系統(tǒng)性污染已被壓低或消除。
本發(fā)明不限于應(yīng)用在時間——距離陣列或地震時間剖面，而是它可應(yīng)用于垂直深度—距離陣列或深度剖面，那里反射時間T已被轉(zhuǎn)換成深度Z，例如使用深度遷移技術(shù)。在那種情況中，該過程保持完全相同，唯獨垂直向前變換是從深度域(Z域)轉(zhuǎn)換成垂直波數(shù)域(Kz域)而反向垂直變換將數(shù)據(jù)從垂直波數(shù)域轉(zhuǎn)換回到深度域。
這里公開的實施例的應(yīng)用示于圖5至圖8中的合成T-X地震數(shù)據(jù)。在這些圖中，垂直軸時間參照是以樣本數(shù)來描述，而水平軸由任意波形道號來描述。
圖5顯示的合成數(shù)據(jù)集類似于圖1A的合成數(shù)據(jù)。圖5顯示在時間—距離域(T-X域)中的合成地震剖面，代表穿過地球的一個垂直切片。一個地震剖面含有大量依賴于時間的波形道，一個挨一個地沿水平方向排列。該地震剖面含有來自水平地震反射體7’的響應(yīng)、來自傾斜反射體8’的響應(yīng)以及來自由于地面散射障礙造成的反射體9’的響應(yīng)。圖5代表沒有由獲取參數(shù)引入畸變的地震反射體響應(yīng)。
圖6顯示與圖5類似的合成數(shù)據(jù)，只是輸入數(shù)據(jù)被畸變，類似于圖2A的示例合成數(shù)據(jù)，這里不均一的源獲取努力造成畸變的輸入數(shù)據(jù)集。圖6顯示一個時間-距離域(T-X域)中的合成地震剖面，代表穿過地球的一個垂直切片。該地震剖面含有來自水平反射體7’的響應(yīng)、來自傾斜反射體8’的響應(yīng)以及來自地球表面的散射障礙的不希望的響應(yīng)9’。圖6代表有由獲取參數(shù)引入的畸變的地震反射體響應(yīng)。
圖6輸入的合成數(shù)據(jù)是由仿真的逐點前進步長源努力系統(tǒng)性污染的，該步長等于圖5、圖6、圖7和圖8中由T8所示8個波形道之間的距離。圖6地震數(shù)據(jù)板顯示代表圖4中的輸入數(shù)據(jù)39。在FX域數(shù)據(jù)的振幅分量中顯而易見由這一周期性調(diào)制造成的畸變。在分解(圖3和圖4中的43)之后，數(shù)據(jù)的相位分量被存儲供其后根據(jù)本發(fā)明的方法對振幅分量濾波之后再使用。與此同時，F(xiàn)X域數(shù)據(jù)的相位分量將不包含系統(tǒng)性污染。
在使用修正的正的FK變換將圖6的數(shù)據(jù)變換到FK域之后，系統(tǒng)性污染將投射為FK空間中的垂直線。然后，這系統(tǒng)性污染將由一個對應(yīng)于該周期性的K濾波器壓低或消除。
為了比較，圖7顯示先有技術(shù)以標準FK變換進行的濾波方法對圖6的畸變數(shù)據(jù)的應(yīng)用，這是不均一的源獲取努力已造成畸變的輸入數(shù)據(jù)集。圖7顯示時間—距離域(T-X域)中的一個合成地震剖面，代表穿過地球的一個垂直切片。該地震剖面包含來自水平地震反射體7’的濾波后的地震響應(yīng)。標準的先有技術(shù)濾波已去掉了水平地震反射體7’的畸變后的數(shù)據(jù)集響應(yīng)，造成穿過地震剖面的均一反射體響應(yīng)的所希望的橫向連續(xù)反射振幅序列。然而，非水平反射體的響應(yīng)，即來自傾斜反射體8’和地球表面散射反射體9’的響應(yīng)，保持不變。
圖8顯示以本發(fā)明的修正FK變換進行的濾波對圖6的畸變數(shù)據(jù)的應(yīng)用，這里不均一的源獲取努力已造成畸變的輸入數(shù)據(jù)集。圖8顯示時間—時間域(T-X域)中的一個地震剖面，代表應(yīng)用本發(fā)明后的一個穿過地球的垂直切片。圖8對應(yīng)于最終數(shù)據(jù)板60，它代表通過圖流程圖中所示序列處理后的數(shù)據(jù)。圖8的地震剖面含有來自水平地震反射體7’的濾波后的地震響應(yīng)，這里標準的先有技術(shù)濾波已去掉了畸變的數(shù)據(jù)集響應(yīng)，造成穿過地震剖面的均一的反射體響應(yīng)。非水平反射體，即來自傾斜反射體8’和地球表面散射反射體9’的響應(yīng)，也已經(jīng)被校正了畸變。與先有技術(shù)的方法相反，該校正也造成了地震剖面傾斜反射體的均一反射體響應(yīng)。
概括根據(jù)圖3的本發(fā)明優(yōu)選實施例，通過在FK域中對地震數(shù)據(jù)振幅分量應(yīng)用一個K濾波器，地震數(shù)據(jù)能被校正由獲取參數(shù)引入的周期性畸變。如圖3中所示，這是通過對有周期性畸變的T-X地震輸入數(shù)據(jù)應(yīng)用垂直FFT(40)來實現(xiàn)的。所造成的F-X數(shù)據(jù)被分解(43)為振幅分量和相位響應(yīng)分量。相位響應(yīng)分量被“保存”(43至56)，用于其后要進行的反演過程。水平FFT46應(yīng)用于F-X數(shù)據(jù)的振幅分量以得到F-K數(shù)據(jù)的振幅分量。因為這一變換是沿X方向進行的，所以這里也稱作水平向前付立葉變換。以要壓低或消除的能量周期為依據(jù)的一個K濾波器應(yīng)用于(51)這一F-K數(shù)據(jù)。例如，該K濾波器可以是一個陷波濾波器或任何其他適當?shù)臑V波器。
濾波之后數(shù)據(jù)被反演(54、56和59)。第一個反演是水平FFT54，這里得到一個振幅分量，以與分解步驟43得到的相位響應(yīng)分量合成(56)。然后，合成步驟56的反演輸出成為垂直IFT59的輸入，由此得到一個T-X地震剖面，其數(shù)據(jù)已被濾波從而去掉了畸變。
再次概括根據(jù)圖4的數(shù)據(jù)流程圖公開說明的本發(fā)明的優(yōu)選實施例通過對修正的FK域數(shù)據(jù)應(yīng)用一個K濾波器，T-X地震數(shù)據(jù)可被校正由獲取參數(shù)引入的周期性畸變。在該優(yōu)選實施例中，一個垂直FFT40應(yīng)用于T-X輸入地震數(shù)據(jù)39，造成實響應(yīng)41和虛響應(yīng)42。使用上述公式1和2。這些響應(yīng)被分解(43)為振幅響應(yīng)44和相位響應(yīng)45。相位響應(yīng)分量被“保存”(45至46)供其后進行的反演過程使用。一個水平FFT46應(yīng)用于F-X數(shù)據(jù)的振幅分量，從而得到振幅響應(yīng)44的實數(shù)據(jù)響應(yīng)47和虛響應(yīng)48。因為這一變換46是沿X方向進行的，所以這里也稱作水平向前付立葉變換。以要壓低或消除的能量周期為依據(jù)的一個K濾波器應(yīng)用于(51)修正的F-K數(shù)據(jù)47和48。例如，該K濾波器可以是一個陷波濾波器或任何其他適當?shù)臑V波器。應(yīng)用濾波器之后的數(shù)據(jù)是修正FK數(shù)據(jù)的實響應(yīng)52和虛響應(yīng)53。
濾波之后，數(shù)據(jù)被反演回到T-X數(shù)據(jù)。第一個反演是水平FFT54，這里得到一個振幅分量55，以與分解步驟43得到的相位響應(yīng)分量45合成。然后，合成步驟56的反演輸出57、58成為垂直IFT59的輸入，由此得到一個T-X地震剖面60，所造成的數(shù)據(jù)已被濾波從而去掉了畸變。圖8是以本發(fā)明處理后的T-X數(shù)據(jù)舉例。
本領(lǐng)域技術(shù)人員將會理解，這里描述的壓低系統(tǒng)性振幅變化的方法可在任何類型數(shù)據(jù)上實踐，包括但不限于諸如地震數(shù)據(jù)和勢場數(shù)據(jù)等地球物理數(shù)據(jù)。再有，應(yīng)該理解，本發(fā)明不是要過度地限定于為說明目的已提出的前述內(nèi)容。對于本領(lǐng)域技術(shù)人員，顯然有各種修改和替代物不脫離如所附權(quán)利要求中定義的本發(fā)明的范圍。
權(quán)利要求
1.一種處理地震數(shù)據(jù)的方法，包含(a)對分解后的地震數(shù)據(jù)應(yīng)用一種水平酉變換，以給出變換后第一分量和變換后第二分量，以及(2)對變換后的第一分量和第二分量中的至少一個進行濾波，該濾波不依賴于變換后的第一分量和第二分量中的另一個。
2.權(quán)利要求1的方法，其中所述水平酉變換包含付立葉變換。
3.權(quán)利要求1的方法，其中所述分解后的地震數(shù)據(jù)是在頻率—空間域中。
4.權(quán)利要求1的方法，其中所述濾波進一步包含應(yīng)用一濾波器以在頻率—波數(shù)域中形成過濾后的分量數(shù)據(jù)。
5.權(quán)利要求1的方法，其中所述濾波進一步包含應(yīng)用一個逆變換形成逆變換后的數(shù)據(jù)，并將所述逆變換后的數(shù)據(jù)與變換后的第一和第二分量中的另一個分量合成，從而形成逆合成數(shù)據(jù)。
6.權(quán)利要求5的方法，進一步包含對所述逆合成數(shù)據(jù)應(yīng)用一個逆變換以形成T-X數(shù)據(jù)。
7.一種處理地震數(shù)據(jù)的方法，包含(a)對分解后的變換后地震數(shù)據(jù)的一個分量應(yīng)用水平酉變換，該變換不依賴于所述分解后的變換后地震數(shù)據(jù)的至少一個其它分量，以得到修正的變換后地震數(shù)據(jù)；以及(b)對所述修正的變換后地震數(shù)據(jù)進行濾波。
8.權(quán)利要求7的方法，進一步包含反演所述濾波后的修正的變換后地震數(shù)據(jù)，以得到T-X地震數(shù)據(jù)。
9.權(quán)利要求7的方法，其中所述濾波進一步包含應(yīng)用一個逆酉變換，以給出反演后的修正的變換后地震數(shù)據(jù)。
10.權(quán)利要求7的方法，進一步包含對所述濾波后的修正的變換后地震數(shù)據(jù)應(yīng)用一個逆水平酉變換以形成逆濾波后分量，并將所述逆濾波后分量數(shù)據(jù)與分解后的變換數(shù)據(jù)的所述至少一個其它分量合成，從而形成合成變換數(shù)據(jù)。
11.權(quán)利要求10的方法，進一步包含對合成變換數(shù)據(jù)應(yīng)用一個逆變換以形成T-X地震數(shù)據(jù)。
12.一種處理地震數(shù)據(jù)的方法，包含(a)對地震數(shù)據(jù)應(yīng)用一個垂直酉變換以給出變換后第一分量和變換后第二分量；(b)對所述變換后第一分量和變換后第二分量應(yīng)用一個水平酉變換；以及(c)對變換后第一分量和第二分量中的至少一個進行濾波，該濾波不依賴于變換后第一分量和第二分量中的另一個。
13.權(quán)利要求12的方法，其中垂直酉變換包含付立葉變換。
14.權(quán)利要求13的方法，其中(i)所述地震數(shù)據(jù)是在時間—空間域中，(ii)應(yīng)用所述付立葉變換進一步包含將所述數(shù)據(jù)變換到頻率—空間域，(iii)所述變換后第一分量包含頻率—空間域中的振幅，以及(iv)應(yīng)用所述水平酉變換進一步包含對所述變換后第一分量應(yīng)用一付立葉變換以給出頻率—波數(shù)域中的第三分量。
15.權(quán)利要求12的方法，其中所述濾波進一步包含應(yīng)用一個K濾波器以形成頻率—波數(shù)域中的濾波后的分量數(shù)據(jù)。
16.權(quán)利要求15的方法，進一步包含對所述濾波后的分量數(shù)據(jù)應(yīng)用一個逆水平變換，以形成逆水平變換數(shù)據(jù)，并將所述逆水平變換數(shù)據(jù)與變換后第一分量和第二分量中的另一個合成以形成逆合成數(shù)據(jù)。
17.權(quán)利要求16的方法，其中逆合成數(shù)據(jù)被逆變換以形成T-X數(shù)據(jù)。
18.一種處理地震數(shù)據(jù)的方法，包含(a)對分解后的地震數(shù)據(jù)振幅分量應(yīng)用一個水平酉變換以得到修正的酉變換數(shù)據(jù)；(b)對所述修正的酉變換數(shù)據(jù)進行濾波；以及(c)反演所述修正的酉變換數(shù)據(jù)以得到T-X地震數(shù)據(jù)。
19.權(quán)利要求18的方法，其中反演所述修正的酉變換數(shù)據(jù)進一步包含應(yīng)用下列變換中的至少一個i)逆水平酉變換，和ii)逆垂直酉變換。
20.權(quán)利要求18的方法，其中反演所述修正的酉變換數(shù)據(jù)進一步包含對所述修正的酉變換數(shù)據(jù)應(yīng)用一個逆水平酉變換以得到逆變換后的修正數(shù)據(jù)，并將所述逆變換后的修正數(shù)據(jù)與地震數(shù)據(jù)的分解后的相位分量合成以形成T-X地震數(shù)據(jù)。
全文摘要
通過對修正的酉變換數(shù)據(jù)進行濾波糾正地震數(shù)據(jù)中由獲取參數(shù)引入的周期性畸變的一種方法。一個垂直變換應(yīng)用于地震數(shù)據(jù)以得到頻率—空間地震數(shù)據(jù)。這些數(shù)據(jù)被分解成振幅分量和相位響應(yīng)分量。相位響應(yīng)分量被保存供反演過程使用。一個水平變換應(yīng)用于該F－X數(shù)據(jù)的振幅分量以得到修正的變換數(shù)據(jù)。因為這一變換是沿X方向進行的，所以這里也稱作水平向前付立葉變換。基于要壓低或消除的能量周期的一個K濾波器應(yīng)用于這個F－K數(shù)據(jù)。然后，這些數(shù)據(jù)被反演以得到其畸變被壓低或消除的T－X數(shù)據(jù)。
文檔編號G01V1/00GK1646940SQ03804108
公開日2005年7月27日 申請日期2003年4月23日 優(yōu)先權(quán)日2002年4月23日
發(fā)明者羅博爾·J·德庫克 申請人:維斯特恩格科地震控股有限公司