共炮點道集多級面波壓制方法和裝置與流程

2023-10-23 14:51:17

本發明涉及地質勘探技術領域，特別涉及一種共炮點道集多級面波壓制方法和裝置。

背景技術：

面波是陸上勘探中常見的幹擾波，主要分布在近炮檢距，呈「掃帚狀」，具有低頻、低速、強振幅和頻散等特點。地震記錄中面波與反射波的同相軸交織在一起，嚴重影響著地震記錄的信噪比，給疊前振幅信息的有效利用帶來極大的幹擾，從而影響了地震資料處理的質量。合理地壓制面波，有利於改善地震資料的品質，提高地震數據的信噪比，對於後續的疊前屬性分析、波形反演以及改善深部地質目標的成像質量都是至關重要的。目前，主要有以下幾種面波壓制方法：

1)基於傳統濾波的面波壓制方法

在數據採集階段可以通過檢波器的組合來壓制面波，但這種方法使得反射波的接收頻帶變窄，而且不能對已經接收到的數據進行處理，缺乏靈活性。隨著地震數據處理技術的發展，出現了多種濾波方法來衰減面波幹擾，主要包括基於頻域濾波類的高通濾波、帶通濾波以及基於視速度特徵差異的FK濾波等。在頻域中反射波與面波的頻譜不是完全分開的，部分區域存在重疊。因而這些方法都存在缺陷，例如，高通濾波在去除低頻面波時也壓制了反射波的低頻成分，而且截頻濾波會引起波形的失真，造成視覺上的假象。為了克服這些缺點，許多濾波方法利用反射波與面波在幾何特性上的差異進行面波壓制。例如：Randon變換將時空域的數據變換到τ-p域中，以直線形式分布的面波在τ-p域中是以一系列能量點的形式存在，切除這些能量點，反變換到時空域就能夠達到壓制面波的目的。但由於「端點效應」的影響，使得面波在τ-p域中不能完全集中，從而導致面波壓制不徹底。

2)基於多尺度分析的面波壓制方法

地震信號是非平穩信號，各種多尺度分析工具(例如：S變換，Wavelet、Ridgelet、Curvelet變換等)也廣泛應用於地震數據的面波衰減中。這類方法基本上都能取得一定的處理效果，但又存在各自的不足。例如：在曲波域中可以綜合利用面波與反射波在空間、尺度、方向上的差異來壓制面波，但隨著尺度的降低，波原子的局域性和方向性都降低，對於面波中低頻低波數成分，曲波低尺度波原子難以提供信噪分離所需的精度，給曲波係數的處理帶來困難；Askari等基於S變換對FK濾波法進行擴展，通過設計隨時間或炮檢距變化的視速度濾波器來壓制面波，對於面波受到地表非均質體散射或者頻散嚴重的地震記錄，能夠取得有效的濾波效果，但該方法運算量較大，非常耗時。

3)基於稀疏表徵策略的面波壓制方法

隨著稀疏表徵技術的發展出現了多種地震信號建模的策略，Karsli等利用線性或非線性的掃頻信號對面波進行參數化建模，通過反褶積的思路來預測並消除面波幹擾。Trad等將單炮記錄描述為具有線性同相軸的面波和具有雙曲同向軸的反射波的疊加，利用線性Radon變換和雙曲Radon變換，結合稀疏約束可以實現面波與反射波波場的分離。此外，利用地震記錄中面波與體波信號波形結構的差異性，利用平穩小波變換和局部離散餘弦變換作稀疏表示字典，將面波與體波分離問題轉化為稀疏優化問題來求解。由於近地表結構變化多樣，該類方法的處理精度依賴於所建模型特徵是否與實際情況相匹配。

由上述描述可知，現有的面波壓制方法主要存在以下問題：

1)由於近地表特徵複雜，依據視速度、頻率、線性特徵等難以完全地區分面波和反射波；

2)多尺度幾何分析工具具有局域性表徵能力，但隨著頻率的降低，局域性越來越差，難以精細地刻畫低頻面波與反射波的特徵，給變換係數的處理帶來困難，導致處理精度不高，並且計算效率低。

3)在現有三維地震勘探觀測系統採樣間隔條件下，當面波視速度較低時，往往存在假頻，現有方法不能有效的壓制面波假頻能量。

針對上述問題，目前尚未提出有效的解決方案。

技術實現要素：

本發明實施例提供了一種共炮點道集多級面波壓制方法，以解決現有的面波壓制方法中所存在的處理精度低、計算效率不高以及面波假頻難以消除的技術問題，該方法包括：

在共炮點道集上確定出面波區域和不含面波的反射波區域；

確定所述面波區域中面波的起始頻率、截至頻率和最低視速度；

在面波的起始頻率至截至頻率範圍內計算所述反射波區域中的每個頻率上振幅的多道統計值；

在面波的起始頻率至截至頻率範圍內對所述面波區域的數據，以所述多道統計值為期望振幅譜，進行零相位濾波處理，得到初步壓制面波的結果；

利用面波的截至頻率、最低視速度、空間採樣間隔和炮檢距計算線性相移算子；

利用所述線性相移算子對所述初步壓制面波的結果進行線性相移、相位濾波和反線性相移處理，得到壓制面波的結果。

在一個實施方式中，在面波的起始頻率至截至頻率範圍內計算所述反射波區域中的每個頻率上振幅的多道統計值，包括：

將所述反射波區域的數據變換至頻率域，得到頻率空間域的數據體；

在面波的起始頻率至截至頻率範圍內計算所述頻率空間域數據體的每個頻率上振幅的多道統計幾何平均值和標準差；

將所述幾何平均值和所述標準差相加，得到面波的起始頻率至截至頻率範圍內的每個頻率上振幅的多道統計值。

在一個實施方式中，在面波的起始頻率至截至頻率範圍內計算所述頻率空間域數據體的每個頻率上振幅的多道統計幾何平均值，包括：

對所述頻率空間域數據體取自然對數；

在面波的起始頻率至截至頻率範圍內的每個頻率上，對取自然對數後的結果計算幾何平均值；

對所述幾何平均值取做e指數處理，得到在面波的起始頻率至截至頻率範圍內的所述頻率空間域數據體的每個頻率上振幅的多道統計幾何平均值。

在一個實施方式中，利用面波的截至頻率、最小視速度、空間採樣間隔和炮檢距計算線性相移算子，包括：

利用面波的截至頻率和空間採樣間隔計算滿足無面波假頻的最小視速度；

確定所述滿足無面波假頻的最小視速度與面波的最小視速度之間的大小關係；

如果所述滿足無面波假頻的最小視速度小於面波的最小視速度，則確定線性相移係數為零；

如果所述滿足無面波假頻的最小視速度大於面波的最小視速度，則將所述滿足無面波假頻的最小視速度的倒數與面波最小視速度的倒數之間的差值作為線性相移係數；

將所述線性相移係數、頻率和炮檢距三者相乘得到線性相移算子。

在一個實施方式中，利用所述線性相移算子將所述初步壓制面波的結果進行線性相移、相位濾波和反線性相移處理，得到壓制面波的結果，包括：

計算所述初步壓制面波結果的相位譜；

對所述初步壓制面波結果的相位譜加上所述線性相移算子得到線性相移後的相位譜；

對所述線性相移後的相位譜，取e的復指數運算得到包含相位信息的信號；

對所述包含相位信息的信號進行FKK域的濾波處理，消除相位信息中的面波信息，並做去e的復指數處理，得到相位濾波後的相位譜；

對所述相位濾波後的相位譜減去線性相移因子，得到反線性相移後的相位譜；

對所述反線性相移後的相位譜與所述初步壓制面波後的振幅譜結合反變換到時間域得到壓制面波的結果。

在一個實施方式中，將所述初步壓制面波的結果進行線性相移、相位濾波和反線性相移處理，得到壓制面波的結果之後，所述方法還包括：

確定得到的所述壓制面波後的結果中面波是否被完全壓制；

如果沒有，則重複執行以下操作，直至面波被完全壓制：將所述壓制面波的結果進行線性相移、相位濾波和反線性相移處理，得到再次面波壓制的結果。

本發明實施例提供了一種共炮點道集多級面波壓制裝置，以解決現有的面波壓制方法中所存在的處理精度低、計算效率不高以及面波假頻難以消除的技術問題，該裝置包括：

第一確定模塊，用於在共炮點道集上確定出面波區域和不含面波的反射波區域；

第二確定模塊，用於確定所述面波區域中面波的起始頻率、截至頻率和最低視速度；

第一計算模塊，用於在面波的起始頻率至截至頻率範圍內計算所述反射波區域中的每個頻率上振幅的多道統計值；

濾波模塊，用於在面波的起始頻率至截至頻率範圍內對所述面波區域的數據，以所述多道統計值為期望振幅譜，進行零相位濾波處理，得到初步壓制面波的結果；

第二計算模塊，用於利用面波的截至頻率、最低視速度、空間採樣間隔和炮檢距計算線性相移算子；

相移濾波模塊，用於利用所述線性相移算子對所述初步壓制面波的結果進行線性相移、相位濾波和反線性相移處理，得到壓制面波的結果。

在一個實施方式中，所述第一計算模塊包括：

變換單元，用於將所述反射波區域的數據變換至頻率域，得到頻率空間域的數據體；

第一計算單元，用於在面波的起始頻率至截至頻率範圍內計算所述頻率空間域數據體的每個頻率上振幅的多道統計幾何平均值和標準差；

相加單元，用於將所述幾何平均值和所述標準差相加，得到面波的起始頻率至截至頻率範圍內的每個頻率上振幅的多道統計值。

在一個實施方式中，所述計算單元包括：

提取子單元，用於對所述頻率空間域數據體取自然對數；

計算子單元，用於在面波的起始頻率至截至頻率範圍內的每個頻率上，對取自然對數後的結果計算幾何平均值；

指出處理子單元，用於對所述幾何平均值取做e指數處理，得到在面波的起始頻率至截至頻率範圍內的所述頻率空間域數據體的每個頻率上振幅的多道統計幾何平均值。

在一個實施方式中，所述第二計算模塊包括：

第二計算單元，用於利用面波的截至頻率和空間採樣間隔計算滿足無面波假頻的最小視速度；

確定單元，用於確定所述滿足無面波假頻的最小視速度與面波的最小視速度之間的大小關係；如果所述滿足無面波假頻的最小視速度小於面波的最小視速度，則確定線性相移係數為零；如果所述滿足無面波假頻的最小視速度大於面波的最小視速度，則將所述滿足無面波假頻的最小視速度的倒數與面波最小視速度的倒數之間的差值作為線性相移係數；

相乘單元，用於將所述線性相移係數、頻率和炮檢距三者相乘得到線性相移算子。

在本發明實施例中，先在共炮點道集上確定出面波區域和不含面波的反射波區域，並確定所述面波區域中面波的起始頻率、截至頻率和最低視速度；然後，在面波的起始頻率至截至頻率範圍內計算所述反射波區域中的每個頻率上振幅的多道統計值；在面波的起始頻率至截至頻率範圍內對所述面波區域的數據，以所述多道統計值為期望振幅譜，進行零相位濾波處理，得到初步壓制面波的結果；利用面波的截至頻率、最低視速度、空間採樣間隔和炮檢距計算線性相移算子；最後，利用所述線性相移算子對所述初步壓制面波的結果進行線性相移、相位濾波和反線性相移處理，得到壓制面波的結果。在上例中，通過基于振幅統計和相移濾波進行有效信號與面波的分離，從而解決了現有的面波壓制方法中所存在的處理精度低、計算效率不高以及面波假頻難以消除的技術問題，達到了方便有效地實現共炮點道集多級面波壓制的目的，且實現了假頻面波的壓制。

附圖說明

此處所說明的附圖用來提供對本發明的進一步理解，構成本申請的一部分，並不構成對本發明的限定。在附圖中：

圖1是根據本發明實施例的共炮點道集多級面波壓制方法的方法流程圖；

圖2是根據本發明實施例的地震數據開時窗示意圖；

圖3是根據本發明實施例的地震數據FK譜分析圖；

圖4是根據本發明實施例的振幅統計濾波壓制面波後的剖面示意圖；

圖5是根據本發明實施例的振幅統計濾波壓制面波後的FK譜；

圖6是根據本發明實施例的相位信息的FK譜示意圖；

圖7是根據本發明實施例的相位信息相移後的結果示意圖；

圖8是根據本發明實施例的面波壓制後的剖面示意圖；

圖9是根據本發明實施例的面波壓制後的FK譜示意圖；

圖10是根據本發明實施例的共炮點道集多級面波壓制裝置的結構框圖。

具體實施方式

為使本發明的目的、技術方案和優點更加清楚明白，下面結合實施方式和附圖，對本發明做進一步詳細說明。在此，本發明的示意性實施方式及其說明用於解釋本發明，但並不作為對本發明的限定。

針對現有的面波壓制方法中所存在的問題，以減少有效信號損失，獲取有益于波形反演等所需的有效低頻信息。在本發明實施例中，綜合利用統計振幅濾波、相移濾波方法進行有效信號與面波的分離研究。主要可以包括：共炮點道集反射波振幅譜統計、零相位濾波、線性相移算子計算以及相移濾波等。通過上述方式實現方便有效的共炮點道集多級面波壓制方法，並實現假頻面波的壓制。

如圖1所示，本例所提供的共炮點道集多級面波壓制方法可以包括以下步驟：

步驟101：在共炮點道集上確定出面波區域和不含面波的反射波區域；

即，可以確定共炮點道集上的兩個時窗，一個是面波所在區域窗，一個是不含面波的反射波區域窗。這主要是考慮到面波區域內也存在反射波，但該區域的反射波被面波汙染了，因此需要確定「不含面波的」反射波區域。

具體地，可以在共炮點道集上對面波進行視速度分析，依據面波視速度範圍開時窗選取存在面波的區域，然後，開時窗選取面波幹擾區外的反射波區域，並計算共炮點道集的FK譜，觀測面波所在的頻帶範圍。

步驟102：確定所述面波區域中面波的起始頻率、截至頻率和最低視速度；

其中，面波的「起始頻率」和「截至頻率」用於進行初步壓制面波的濾波處理，面波的「最低視速度」和「截至頻率」用於計算線性相移算子。

步驟103：在面波的起始頻率至截至頻率範圍內計算所述反射波區域中的每個頻率上振幅的多道統計值；

即。對反射波窗內的地震記錄，統計頻域下每個頻率上的振幅譜值。

具體地，在面波的起始頻率至截至頻率範圍內計算所述反射波區域中的每個頻率上振幅的多道統計值，可以包括以下步驟：

S1：將所述反射波區域的數據變換至頻率域，得到頻率空間域的數據體；

S2：在面波的起始頻率至截至頻率範圍內計算所述頻率空間域數據體的每個頻率上振幅的多道統計幾何平均值和標準差；

其中，該多道統計幾何平均值可以按照以下方式確定，包括以下步驟：

S2-1：對所述頻率空間域數據體取自然對數；

S2-2：在面波的起始頻率至截至頻率範圍內的每個頻率上，對取自然對數後的結果計算幾何平均值；

S2-3：對所述幾何平均值取做e指數處理，得到在面波的起始頻率至截至頻率範圍內的所述頻率空間域數據體的每個頻率上振幅的多道統計幾何平均值。

S3：將所述幾何平均值和所述標準差相加，得到面波的起始頻率至截至頻率範圍內的每個頻率上振幅的多道統計值。

步驟104：在面波的起始頻率至截至頻率範圍內對所述面波區域的數據，以所述多道統計值為期望振幅譜，進行零相位濾波處理，得到初步壓制面波的結果；

即，在對所述的面波窗內的數據在面波所在頻帶的範圍上，利用統計值進行濾波處理。

步驟105：利用面波的截至頻率、最低視速度、空間採樣間隔和炮檢距計算線性相移算子；

在一個實施方式中，可以按照以下方式計算線性相移算子，包括以下步驟：

S1：利用面波的截至頻率和空間採樣間隔計算滿足無面波假頻的最小視速度；

S2：確定所述滿足無面波假頻的最小視速度與面波的最小視速度之間的大小關係；

S3：如果所述滿足無面波假頻的最小視速度小於面波的最小視速度，則確定線性相移係數為零；

S4：如果所述滿足無面波假頻的最小視速度大於面波的最小視速度，則將所述滿足無面波假頻的最小視速度的倒數與面波最小視速度的倒數之間的差值作為線性相移係數；

S5：將所述線性相移係數、頻率和炮檢距三者相乘得到線性相移算子。

步驟106：根據所述線性相移算子對所述初步壓制面波的結果進行線性相移、相位濾波和反線性相移處理，得到壓制面波的結果。

具體地，可以計算所述初步壓制面波結果的相位譜；對所述初步壓制面波結果的相位譜加上所述線性相移算子得到線性相移後的相位譜；對所述線性相移後的相位譜，取e的復指數運算得到包含相位信息的信號；對所述包含相位信息的信號進行FKK域的濾波處理，消除相位信息中的面波信息，並做去e的復指數處理，得到相位濾波後的相位譜；對所述相位濾波後的相位譜減去線性相移因子，得到反線性相移後的相位譜；對所述反線性相移後的相位譜與所述初步壓制面波後的振幅譜結合反變換到時間域得到壓制面波的結果。

在步驟S106之後，還可以確定得到的所述壓制面波後的結果中面波是否被完全壓制；如果沒有，則重複執行以下操作，直至面波被完全壓制：將所述壓制面波的結果進行線性相移、相位濾波和反線性相移處理，得到再次面波壓制的結果。

在上例中，綜合利用統計振幅約束和相移濾波技術，可以獲得比傳統面波壓制方法更有效的結果，為後續的疊前屬性分析、波形反演、深層目標成像以及寬頻解釋等提供可靠的數據基礎，從而為後續振幅信息的有效利用提供了一種更為精確的面波壓制手段。

在上述的計算過程中，利用反射波振幅的統計特性壓制面波振幅以及採用相移濾波來進一步衰減面波及其假頻的能量。同常規地震數據面波壓制方法相比，基于振幅統計和相移濾波的方法能夠更好地利用數據本身的信息來抑制噪聲，並對面波假頻具有很好的壓制效果，因此選擇利用反射波的振幅譜統計和相位譜相移濾波來綜合衰減面波及其假頻能量。利用該方法可以有效壓制面波能量，並保持反射信號的低頻有效信息，同時能夠自適應三維探區的激發能量的空間變化，具有較高信噪分離的精度，對面波假頻能量也具有較好的壓制能力。

下面結合一具體實施例對共炮點道集多級面波壓制方法進行說明，然而值得注意的是，該具體實施例僅是為了更好地說明本發明，並不構成對本發明的不當限定。

在本發明實施例中提供了一種基於共炮點道集振幅統計和相移濾波的面波壓制方法，可以包以下步驟：

S1：在共炮點道集上對面波進行視速度分析，依據面波視速度範圍開時窗選取存在面波的區域，然後，開時窗選取面波幹擾區外的反射波區域，並計算共炮點道集的FK譜，觀測面波所在的頻帶範圍。

如圖2所示是共炮點道集，圖2中白色虛線與底線圍成的區域是依據面波視速度範圍選出的面波幹擾區，其中，白色實線、白色虛線與底線圍成的區域是選出的面波幹擾區域外的反射波區域；對該共炮點道集變換到頻率波數域，得到FK譜，其結果可以如圖3所示，從FK譜分析的結果中可以得到存在面波幹擾的頻帶為5-15Hz。

S2：將面波區域外的反射波數據變換到頻率域，得到頻率空間域的數據體，計算頻率空間域數據體的面波所在頻帶範圍內的每個頻率上振幅的多道統計幾何平均值和標準差，將統計幾何平均值和標準差相加得到面波所在頻帶範圍內的每個頻率上的反射波振幅譜的統計值。

將圖2中白色實線、白色虛線與底線圍成的反射波區域中的數據變換到頻率域，得到頻率空間域的數據體，該數據體上5-15Hz範圍內的每一個頻率片都對應著一組隨空間位置變化的振幅值，計算這組數據取自然對數後的幾何平均值，然後取e指數處理得到統計幾何平均值，由該平均值計算數組的標準差，將幾何平均值和標準差相加，便得到了該頻率下反射波振幅譜的統計值，同理可以計算出其他頻率下的反射波振幅譜的統計值。

S3：對面波區域的數據進行傅立葉變換，並對變換後得到的每一道數據記錄振幅譜，在S1確定的面波頻帶範圍內，以S2中的反射波振幅譜統計值作為期望，設計零相位濾波算子進行濾波處理，並反變換到時間域以得到振幅譜上壓制面波的結果。

將圖2中白色虛線與底線圍成的面波區域中的數據進行傅立葉變換，並對每一道數據記錄振幅譜，在5-15Hz的頻帶範圍內，以S2中反射波振幅譜的統計值作為期望，設計零相位濾波算子進行濾波處理，得到振幅譜上壓制面波的結果，結果如圖4所示。從振幅譜上壓制面波的結果上可以看出經過振幅譜上壓制面波處理後，壓制掉了大部分的面波能量。計算振幅譜上壓制面波結果的FK譜可以如圖5所示。圖3和圖5的對比，表明振幅統計濾波處理對面波及其假頻能量都起到了有效的壓製作用。

S4：對S3的輸出數據進行相位特徵的FK譜分析，依據面波的截至頻率、最低視速度、空間採樣間隔和炮檢距計算線性相移算子θm(f,x)，對S3的輸出數據進行相移處理，計算相移後的相位特徵的FK譜，依據面波特徵，設計濾波算子在FK譜上消除相位信息中的面波。

如圖6所示對S3的輸出數據進行相位特徵的FK譜分析。對於面波而言，由於空間採樣不足，圖中FK譜上面波存在假頻，假頻能量與有效信號能量重疊在一起。依據面波的截至頻率、最低視速度、空間採樣間隔和炮檢距計算線性相移算子θm(f,x)，對S3的輸出數據的相位信息進行相移，相移後的相位特徵的FK譜可以如圖7所示。圖中結果表明相移後的相位信息的假頻特徵消失，這有利於面波的壓制。設計濾波算子在FK譜上消除相位信息中的面波特徵，得到了壓制面波後的相位信息結果。

S5：對步驟S4中消除面波特徵後的相位信息利用三角函數進行-θm(f,x)相移，並計算S3的輸出結果的振幅譜，將計算得到的振幅譜與相移後的相位譜反變換到時間域，得到如圖8所示的壓制面波的結果輸出。

S6：針對不同假頻特徵情況，考察S5輸出的面波壓制結果及其FK譜中面波特徵是否被完全消除，如果仍有面波殘留，則重複進行S4和S5，直至獲得滿意的輸出結果為止。

對比圖2和圖8可以發現面波及其假頻能量被極大地壓制，反射信號的連續性明顯增強，提高了信噪比，並對圖8的數據變換到頻率波數域，得到如圖9所示的FK譜。對比圖3和圖9可以發現頻譜上的面波及其假頻的強能量基本上完全衰減掉，有效信號的低頻能量並沒有明顯損失。

在本例中，利用反射波振幅的統計特性壓制面波振幅以及採用相移濾波來進一步衰減面波及其假頻的能量。同常規地震數據面波壓制方法相比，基于振幅統計和相移濾波的方法能夠更好地利用數據本身的信息來抑制噪聲，並對面波假頻具有很好的壓制效果，該方法具有較高的處理精度和較廣的應用範圍。

基於同一發明構思，本發明實施例中還提供了一種共炮點道集多級面波壓制裝置，如下面的實施例所述。由於共炮點道集多級面波壓制裝置解決問題的原理與共炮點道集多級面波壓制方法相似，因此共炮點道集多級面波壓制裝置的實施可以參見共炮點道集多級面波壓制方法的實施，重複之處不再贅述。以下所使用的，術語「單元」或者「模塊」可以實現預定功能的軟體和/或硬體的組合。儘管以下實施例所描述的裝置較佳地以軟體來實現，但是硬體，或者軟體和硬體的組合的實現也是可能並被構想的。圖10是本發明實施例的共炮點道集多級面波壓制裝置的一種結構框圖，如圖10所示，可以包括：第一確定模塊1001、第二確定模塊1002、第一計算模塊1003、濾波模塊1004、第二計算模塊1005和相移濾波模塊1006，下面對該結構進行說明。

第一確定模塊1001，用於在共炮點道集上確定出面波區域和不含面波的反射波區域；

第二確定模塊1002，用於確定所述面波區域中面波的起始頻率、截至頻率和最低視速度；

第一計算模塊1003，用於在面波的起始頻率至截至頻率範圍內計算所述反射波區域中的每個頻率上振幅的多道統計值；

濾波模塊1004，用於在面波的起始頻率至截至頻率範圍內對所述面波區域的數據，以所述多道統計值為期望振幅譜，進行零相位濾波處理，得到初步壓制面波的結果；

第二計算模塊1005，用於利用面波的截至頻率、最低視速度、空間採樣間隔和炮檢距計算線性相移算子；

相移濾波模塊1006，用於利用所述線性相移算子對所述初步壓制面波的結果進行線性相移、相位濾波和反線性相移處理，得到壓制面波的結果。

在一個實施方式中，第一計算模塊1003可以包括：變換單元，用於將所述反射波區域的數據變換至頻率域，得到頻率空間域的數據體；第一計算單元，用於在面波的起始頻率至截至頻率範圍內計算所述頻率空間域數據體的每個頻率上振幅的多道統計幾何平均值和標準差；相加單元，用於將所述幾何平均值和所述標準差相加，得到面波的起始頻率至截至頻率範圍內的每個頻率上振幅的多道統計值。

在一個實施方式中，上述計算單元可以包括：提取子單元，用於對所述頻率空間域數據體取自然對數；計算子單元，用於在面波的起始頻率至截至頻率範圍內的每個頻率上，對取自然對數後的結果計算幾何平均值；指出處理子單元，用於對所述幾何平均值取做e指數處理，得到在面波的起始頻率至截至頻率範圍內的所述頻率空間域數據體的每個頻率上振幅的多道統計幾何平均值。

在一個實施方式中，上述第二計算模塊可以包括：第二計算單元，用於利用面波的截至頻率和空間採樣間隔計算滿足無面波假頻的最小視速度；確定單元，用於確定所述滿足無面波假頻的最小視速度與面波的最小視速度之間的大小關係；如果所述滿足無面波假頻的最小視速度小於面波的最小視速度，則確定線性相移係數為零；如果所述滿足無面波假頻的最小視速度大於面波的最小視速度，則將所述滿足無面波假頻的最小視速度的倒數與面波最小視速度的倒數之間的差值作為線性相移係數；相乘單元，用於將所述線性相移係數、頻率和炮檢距三者相乘得到線性相移算子。

從以上的描述中，可以看出，本發明實施例實現了如下技術效果：先在共炮點道集上確定出面波區域和不含面波的反射波區域，並確定所述面波區域中面波的起始頻率、截至頻率和最低視速度；然後，在面波的起始頻率至截至頻率範圍內計算所述反射波區域中的每個頻率上振幅的多道統計值；在面波的起始頻率至截至頻率範圍內對所述面波區域的數據，以所述多道統計值為期望振幅譜，進行零相位濾波處理，得到初步壓制面波的結果；利用面波的截至頻率、最低視速度、空間採樣間隔和炮檢距計算線性相移算子；最後，利用所述線性相移算子對所述初步壓制面波的結果進行線性相移、相位濾波和反線性相移處理，得到壓制面波的結果。在上例中，通過基于振幅統計和相移濾波進行有效信號與面波的分離，從而解決了現有的面波壓制方法中所存在的處理精度低、計算效率不高以及面波假頻難以消除的技術問題，達到了方便有效地實現共炮點道集多級面波壓制的目的，且實現了假頻面波的壓制。

顯然，本領域的技術人員應該明白，上述的本發明實施例的各模塊或各步驟可以用通用的計算裝置來實現，它們可以集中在單個的計算裝置上，或者分布在多個計算裝置所組成的網絡上，可選地，它們可以用計算裝置可執行的程序代碼來實現，從而，可以將它們存儲在存儲裝置中由計算裝置來執行，並且在某些情況下，可以以不同於此處的順序執行所示出或描述的步驟，或者將它們分別製作成各個集成電路模塊，或者將它們中的多個模塊或步驟製作成單個集成電路模塊來實現。這樣，本發明實施例不限制於任何特定的硬體和軟體結合。

以上所述僅為本發明的優選實施例而已，並不用於限制本發明，對於本領域的技術人員來說，本發明實施例可以有各種更改和變化。凡在本發明的精神和原則之內，所作的任何修改、等同替換、改進等，均應包含在本發明的保護範圍之內。