本發(fā)明涉及地震數據處理技術領域，尤其涉及一種采集腳印壓制方法及系統(tǒng)。

背景技術：

隨著巖性油氣藏和隱蔽性油氣藏勘探的不斷深入，若想對油氣藏和隱蔽性油氣藏的分布得出準確的油藏評價結果和儲層預測方案，需要采集準確的地震數據，因此，在采集地震數據的過程中，對地震成像分辨率和地震資料的細節(jié)刻畫提出了較高的要求。采集腳印噪聲又稱采集痕跡噪聲，是一種由于人為因素產生的地震噪聲，具體是指在地震資料采集與處理過程中留下的人為痕跡，表現方式為地震剖面或切片上不同面源位置處規(guī)則性的振幅變化假象。由于采集腳印噪聲對儲層預測有一定的影響，它的存在會模糊砂體邊界和構造特征，因此，該采集腳印噪聲會影響地震成像及地震屬性分析，進而造成一些巖性體分布的假象，對油藏評價和儲層預測產生了影響。

目前，經典的壓制采集腳印方法是指在傳統(tǒng)道集疊加后的地震數據體上，應用三維錐形濾波等地震處理技術來消除地震數據中以垂直和水平條紋形式出現的采集腳印噪聲。假設三維地震數據體的炮線和檢波線是正交的，抽取一條炮線和檢波線上的所有地震道組成一個正交子集，相當于對地下局部三維地質體進行單次覆蓋地震采集，該正交子集所構建的時空域稱為十字交叉排列域，其炮線和檢波線的交點稱為交叉點。具體的，首先將三維地震數據體抽取到十字交叉排列域，其次，利用三維錐形濾波技術來壓制地震數據采集過程中出現的線性噪聲，突出有效信號，進而壓制掉一部分的采集腳印噪聲。

然而，上述三維錐形濾波技術中對于采集腳印噪聲的壓制有其局限性，無法徹底消除采集腳印噪聲的影響，而且由炮線和檢波線組成的十字交叉排列的稀疏性會影響去噪的保真度，同時會產生空間假頻，不符合高保真處理要求，所取得的結果也很難完全滿足振幅、頻率、相位等屬性在油氣識別中 的需要。

技術實現要素：

本發(fā)明提供一種采集腳印壓制方法及系統(tǒng)，以解決現有采集腳印壓制方法無法完整消除采集腳印的影響，且不符合高保真處理要求致使地震數據處理結果不能完全滿足振幅、頻率、相位等屬性在油氣識別中需要的問題。

本發(fā)明提供的一種采集腳印壓制方法，包括：

將待處理的地震數據按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內，得到N個共偏移距共方位角域道集，所述共偏移距共方位角域道集是指具有相同偏移距和相同方位角的地震道的集合，N為相同偏移距和相同方位角區(qū)域被覆蓋的次數，N為大于或等于1的整數；

分別將所述炮線之間的距離和所述檢波線之間的距離進行減半處理，根據每個所述共偏移距共方位角域道集中的每個地震道的炮點和/或檢波點的坐標，計算出新增加的炮點和/或檢波點的坐標，每個所述共偏移距共方位角域道集生成四個子共偏移距共方位角域道集；

對每個所述子共偏移距共方位角域道集進行數據插值，得到加密數據，所述加密數據包括原始地震數據和插值后的數據；

將所述加密數據按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內，得到4N個子共偏移距共方位角域道集；

對所述4N個子共偏移距共方位角域道集進行疊加前偏移處理，得到壓制采集腳印后的地震數據分析結果。

本發(fā)明還提供一種采集腳印壓制系統(tǒng)，包括：

獲取模塊，用于將待處理的地震數據按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內，得到N個共偏移距共方位角域道集，所述共偏移距共方位角域道集是指具有相同偏移距和相同方位角的地震道的集合，N為相同偏移距和相同方位角區(qū)域被覆蓋的次數，N為大于或等于1的整數；

處理模塊，用于分別將所述炮線之間的距離和所述檢波線之間的距離進行減半處理，根據每個所述共偏移距共方位角域道集中的每個地震道的炮點和/或檢波點的坐標，計算出新增加的炮點和/或檢波點的坐標，每個所述共偏移距共方位角域道集生成四個子共偏移距共方位角域道集；

插值模塊，用于對每個所述子共偏移距共方位角域道集進行數據插值，得到加密數據，所述加密數據包括原始地震數據和插值后的數據；

所述獲取模塊，還用于將所述加密數據按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內，得到4N個子共偏移距共方位角域道集；

所述處理模塊，還用于對所述4N個子共偏移距共方位角域道集進行疊加前偏移處理，得到壓制采集腳印后的地震數據分析結果。

本發(fā)明采集腳印壓制方法及系統(tǒng)，通過將待處理的地震數據分選到十字排列域內，并根據地震道與偏移距和方位角的關系，得到了多個共偏移距共方位角域道集，通過分別將炮線距和檢波線距進行減半處理，且根據每個共偏移距共方位角域道集中的每個地震道的炮點和/或檢波點的坐標，可計算出新增加的炮點和/或檢波點的坐標，進而使每個共偏移距共方位角域道集生成四個子共偏移距共方位角域道集，分別對子共偏移距共方位角域道集進行數據插值，得到包括原始地震數據和插值后的數據的加密數據，最后將該加密數據重新分選到由減半后的炮線距和檢波線距組成的十字排列域內，得到四倍于共偏移距共方位角域道集的子共偏移距共方位角域道集，對所有的子共偏移距共方位角域道集進行疊加前偏移處理，便可得到壓制采集腳印后的地震數據分析結果。本發(fā)明實施例的技術方案，利用共偏移距共方位角域道集的單次覆蓋特性，通過在同一偏移距和同一方位角內進行插值，實現了加密炮線和檢波線的目的，實現方案簡單、容易實現，對壓制采集腳印、淺層成像和高保真處理有重要的意義，符合高保真處理要求，所取得的結果也完全滿足振幅、頻率、相位等屬性在油氣識別中的需要。

附圖說明

為了更清楚地說明本發(fā)明實施例或現有技術中的技術方案，下面將對實施例或現有技術描述中所需要使用的附圖作一簡單地介紹，顯而易見地，下面描述中的附圖是本發(fā)明的一些實施例，對于本領域普通技術人員來講，在不付出創(chuàng)造性勞動性的前提下，還可以根據這些附圖獲得其他的附圖。

圖1為本發(fā)明提供的采集腳印壓制方法實施例一的流程示意圖；

圖2為某一十字排列域及該十字排列域上的共偏移距共方位角域道集；

圖3為從十字排列域中抽取的某一共偏移距共方位角域道集；

圖4為炮線距和檢波線距減半前形成的十字排列域；

圖5為炮線距和檢波線距減半后形成的十字排列域；

圖6為共偏移距共方位角域道集進行數據插值前的數據分布示意圖；

圖7為共偏移距共方位角域道集進行數據插值后的數據分布示意圖；

圖8為原始地震數據對共偏移距共方位角域道集的劃分示意圖；

圖9為加密數據對共偏移距共方位角域道集的劃分示意圖；

圖10為原始地震數據在共偏移距共方位角域道集上形成的項鏈圖；

圖11為加密數據在共偏移距共方位角域道集上形成的項鏈圖；

圖12為原始地震數據在共偏移距共方位角域道集上的偏移距范圍；

圖13為加密數據分別在相應的四個子共偏移距共方位角域道集上的偏移距范圍；

圖14為原始地震數據的共偏移距偏移處理成果切片得到的地震數據分析結果；

圖15為加密數據的共偏移距偏移處理成果切片得到的地震數據分析結果；

圖16為本發(fā)明提供的采集腳印壓制系統(tǒng)實施例一的結構示意圖。

具體實施方式

為使本發(fā)明實施例的目的、技術方案和優(yōu)點更加清楚，下面將結合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施例是本發(fā)明一部分實施例，而不是全部的實施例?；诒景l(fā)明中的實施例，本領域普通技術人員在沒有作出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都屬于本發(fā)明保護的范圍。

三維地震勘探是一種高密度面積采集技術，是三維體積勘探。它利用炮點和檢波點網格的靈活組合獲得分布均勻的地下CDP點網格和確定的覆蓋次數。在三維地震數據采集中，炮檢點有規(guī)律地變化，會造成在疊加或偏移后的地震數據上地震波振幅出現有規(guī)律性的變化，在時間和深度切片上表現為振幅周期性變化，這種周期性振幅假象就是采集腳印。采集腳印對儲層預測有一定的影響，它的存在模糊了砂體邊界和構造特征，造成一些巖性體分布 的假象，對油藏評價和儲層預測了產生影響。

觀測系統(tǒng)參數的選擇特別是炮線距和檢波線距的選擇，會影響采集腳印現象的產生，若炮線距和檢波線距過大，則會產生較嚴重的采集腳印，并且炮線距的影響大于檢波線距的影響。在實際應用中，炮線距和檢波線距一般選擇為道距的4～6倍。

野外三維地震數據采集一旦完成，炮線距、檢波線距和道距就已確定。因此，采集到的地震數據中是否存在采集腳印主要取決于炮檢線距的大小。但是，一般采集地震數據時為了節(jié)省成本，通常采用較大炮檢線距，因此會出現采集腳印。

現階段，經典的采集腳印壓制方法是在傳統(tǒng)地震道集疊加后的地震數據體上，應用三維錐形濾波等地震處理技術，消除地震數據中以垂直和水平條紋形式出現的采集腳印噪聲。實際上，三維錐形濾波地震處理技術是在采集腳印已進入成像過程中的“事后處理”思路，其對于采集腳印現象的壓制有一定的局限性，無法通過處理完整消除采集腳印的影響，而且不符合高保真處理要求，所取得的結果很難完全滿足振幅、頻率、相位等屬性在油氣識別中的需要。

針對現有采集腳印壓制方法存在的問題和缺陷，本發(fā)明提供了一種采集腳印壓制方法及系統(tǒng)，通過在共偏移距共方位角域內插炮檢線，將炮線距和檢波線距縮小一倍，通過在共偏移距共方位角域內進行高精度插值，并按照新的炮線距和檢波線距重新抽取共偏移距共方位角域道集，獲得的道集個數是原來的4倍，即覆蓋次數增大了原來的4倍，而且每個道集的偏移距范圍縮減為原來的一半。這種高密度、小線距、小偏移距范圍的共偏移距共方位角域道集在地震數據疊前偏移時能有效阻止采集腳印進入到成像，是一種“源頭控制”的“事前壓制”思路，解決了現有三維錐形濾波技術存在的消除采集腳印噪聲不徹底，且不符合高保真處理要求造成的所取得的結果難以滿足振幅、頻率、相位等屬性在油氣識別中需要的問題。

圖1為本發(fā)明提供的采集腳印壓制方法實施例一的流程示意圖。如圖1所示，本發(fā)明實施例一提供的采集腳印壓制方法，包括：

步驟101：將待處理的地震數據按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內，得到N個共偏移距共方位角域道集；

其中，共偏移距共方位角域道集是指具有相同偏移距和相同方位角的地震道的集合，N為相同偏移距和相同方位角區(qū)域被覆蓋的次數，N為大于或等于1的整數。

具體的，根據地震資料從觀測系統(tǒng)中抽取一條炮線和檢波線，炮線由一列炮點組成，檢波線由一列檢波點組成。炮線和檢波線交叉形成十字排列域，交叉的方式可以是正交，也可以是斜交，本發(fā)明實施例并不對此進行限定。因此，將待處理的地震數據按照偏移距和方位角分選到上述十字排列域內，可以得到多個共偏移距共方位角域道集，該共偏移距共方位角域道集的數量與具有相同偏移距和相同方位角的區(qū)域被覆蓋的次數相關。

實際應用中，覆蓋次數是指炮點對地下反射界面同一個點的追蹤次數，反映到十字排列域就是地震數據分選的每個網格的次數。因此，共偏移距共方位角域內的每個地震道的數據的集合組成共偏移距共方位角域道集。

由上述可知，十字排列域包含了多個共偏移距共方位角域道集，從該十字排列域中抽取任一共偏移距共方位角域道集，該第一共偏移距共方位角域道集在空間展開后具有單次覆蓋的特征。

本發(fā)明實施例主要是對具有單次覆蓋特征的共偏移距共方位角域道集進行處理來減低或消除采集腳印噪聲。

步驟102：分別將炮線之間的距離和檢波線之間的距離進行減半處理，根據每個共偏移距共方位角域道集中的地震道的炮點和/或檢波點的坐標，計算出新增加的炮點和/或檢波點的坐標，每個共偏移距共方位角域道集生成四個子共偏移距共方位角域道集；

將組成十字排列域的炮線距和檢波線距分別減半處理，并根據該共偏移距共方位角域道集中每個地震道，提取出該地震道的炮點和/或檢波點的坐標，根據該地震道的方位角及其對應的炮點和/或檢波點的坐標，可以計算出新增加的炮點和/或檢波點的坐標。處理完成后，可使每個共偏移距共方位角域道集生成四個子共偏移距共方位角域道集。

步驟103：對每個子共偏移距共方位角域道集進行數據插值，得到加密數據；

其中，該加密數據包括原始地震數據和插值后的數據。

由于同一共偏移距共方位角域道集，其具有相同的方位角，因此，在子 共偏移距共方位角域道集進行插值得到的加密數據，其一致性顯著好于直接在共偏移距共方位角域道集內進行插值得到的數據，因此，按照本發(fā)明實施例提供的方法進行數據插值的精度明顯高于常規(guī)的共偏移距域插值的精度。

步驟104：將上述加密數據按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內，得到4N個子共偏移距共方位角域道集；

本實施例將加密數據按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內，重新抽取共偏移距共方位角域道集，共獲得4N個子共偏移距共方位角域道集，該道集個數是原來的4倍，即具有相同偏移距和相同方位角區(qū)域的覆蓋次數增大了原來的4倍。這種高密度、小線距、小偏移距范圍的子共偏移距共方位角域道集，等價于炮線距和檢波線距減半的野外采集。

由步驟103可知，插值后的數據加上原始數據稱為加密數據。加密數據按新的炮線距和檢波線距重新劃分共偏移距共方位角域，再抽取共偏移距共方位角域道集。由于炮線距和檢波線距縮小為原來的一倍，原來的一個共偏移距共方位角域一化為四。

步驟105：對上述4N個子共偏移距共方位角域道集進行疊加前偏移處理，得到壓制采集腳印后的地震數據分析結果。

具體的，對上述4N個子共偏移距共方位角域道集重新進行偏移、疊加，再提取振幅切片，進而可得到壓制采集腳印后的地震數據分析結果。

本發(fā)明實施例一提供的采集腳印壓制方法，通過將待處理的地震數據分選到十字排列域內，并根據地震道與偏移距和方位角的關系，得到了多個共偏移距共方位角域道集，通過分別將炮線距和檢波線距進行減半處理，且根據每個共偏移距共方位角域道集中的每個地震道的炮點和/或檢波點的坐標，可計算出新增加的炮點和/或檢波點的坐標，進而使每個共偏移距共方位角域道集生成四個子共偏移距共方位角域道集，分別對子共偏移距共方位角域道集進行數據插值，得到包括原始地震數據和插值后的數據的加密數據，最后將該加密數據重新分選到由減半后的炮線距和檢波線距組成的十字排列域內，得到四倍于共偏移距共方位角域道集的子共偏移距共方位角域道集，對所有的子共偏移距共方位角域道集進行疊加前偏移處理，便可得到壓制采集 腳印后的地震數據分析結果。本發(fā)明實施例的技術方案，利用共偏移距共方位角域道集的單次覆蓋特性，通過在同一偏移距和同一方位角內進行插值，實現了加密炮線和檢波線的目的，實現方案簡單、容易實現，對壓制采集腳印、淺層成像和高保真處理有重要的意義，符合高保真處理要求，所取得的結果也完全滿足振幅、頻率、相位等屬性在油氣識別中的需要。

進一步的，在上述實施例一提供的采集腳印壓制方法中，上述步驟102中的，根據每個共偏移距共方位角域道集中的每個地震道的炮點和/或檢波點的坐標，計算出新增加的炮點和/或檢波點的坐標，可具體包括以下三種情況：

情況一：通過改變共偏移距共方位角域道集中每個地震道的炮點坐標，而不改變所述共偏移距共方位角域道集中每個地震道的檢波點坐標，計算出新增加的炮點的坐標。

具體的，可采用下述坐標變換公式(1)得到新增加的炮點的坐標：

其中，x_s和y_s分別為共偏移距共方位角域道集中每個地震道的炮點的橫坐標和縱坐標，x和y分別為新增加的炮點的橫坐標和縱坐標，θ為共偏移距共方位角域道集中炮線的方位角，dl_s為共偏移距共方位角域道集中的炮線距。

情況二：通過改變共偏移距共方位角域道集中每個地震道的檢波點坐標，而不改變共偏移距共方位角域道集中每個地震道的炮點坐標，計算出新增加的檢波點的坐標。

詳細的，通過采用下述坐標變換公式(2)得到新增加的檢波點的坐標：

其中，x_r和y_r分別為共偏移距共方位角域道集中每個地震道的檢波點的橫坐標和縱坐標，x′和y′分別為新增加的檢波點的橫坐標和縱坐標，α為共偏移距共方位角域道集中檢波線的方位角，dl_r為共偏移距共方位角域道集中的檢波線距。

情況三：通過改變共偏移距共方位角域道集中每個地震道的炮點坐標和 檢波點坐標，計算出新增加的炮點的坐標和檢波點的坐標。

具體的，在該情況下可具體通過情況一中的公式(1)求出新增加的炮點的坐標，通過情況二中的公式(2)求出新增加的檢波點的坐標。

進一步的，在上述實施例一提供的采集腳印壓制方法中，上述步驟103，也即，對每個子共偏移距共方位角域道集進行數據插值，得到加密數據，具體包括：

利用五維插值法對每個子共偏移距共方位角域道集進行數據插值處理，得到加密數據。

在地震資料數字處理中，為了保證采集到的地震數據函數在時間和空間上的連續(xù)性，往往需要對其進行插值處理，根據插值原理，應用較少的地震數據作為參數值，確定出任意時間和空間的函數值。由于同一共偏移距共方位角域道集內數據具有相同的方位角，因此數據的一致性顯著好于共偏移距域道集，故數據插值的精度明顯高于常規(guī)的共偏移距域插值。

地震數據體可以實現五維采集，并且許多地震處理技術也可在五維空間內實現。地震處理技術大都對地震數據空間分布的規(guī)則性有較嚴格的要求，但地震數據往往由于各種原因，空間采樣是不規(guī)則或缺失的，因而會造成地震數據處理結果較差。因此，本發(fā)明實施例采用不規(guī)則觀測系統(tǒng)五維插值算法對地震數據進行插值重建，為地震數據處理提供高質量的資料，提高了地震數據處理的效果，精度很高。

值得說明的是，具體的插值方法可視處理系統(tǒng)的需要而定，本發(fā)明并不對此進行限定。

下面以采集腳印壓制方法在實際中的應用為例進行說明。

具體的，利用的地震數據資料具體為：炮線之間的距離為500米，檢波線之間的距離為400米，炮線和檢波線組成的十字排列域的覆蓋次數為180次，炮線的方位角為88度，檢波線的方位角為358度。

本發(fā)明該實施例的具體實現步驟如下：

第一：將待處理的地震數據按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內，得到180個共偏移距共方位角域道集；

在該實施例中，該十字排列域共被180個共偏移距共方位角域道集剖分，個數與由炮線和檢波線組成的十字排列域的覆蓋次數相同。

圖2為某一十字排列域及該十字排列域上的共偏移距共方位角域道集。如圖2所示，該十字排列域上共有180個不同顏色的方塊，其代表180個共偏移距共方位角域道集，表明相同偏移距和相同方位角區(qū)域被覆蓋的次數為180次。

第二：分別將炮線距和檢波線距進行減半處理，根據每個共偏移距共方位角域道集中的地震道的炮點和/或檢波點的坐標，計算出新增加的炮點和/或檢波點的坐標，每個共偏移距共方位角域道集生成四個子共偏移距共方位角域道集；

圖3為從十字排列域中抽取的某一共偏移距共方位角域道集。如圖3所示，從上述十字排列域中抽取該共偏移距共方位角域道集在空間的展布具有單次覆蓋的特征。值得說明的是，本發(fā)明實施例可任意抽取一個共偏移距共方位角域道集，只要其在空間的展布滿足單次覆蓋特征均可，本發(fā)明實施例并不對此進行限定。

圖4為炮線距和檢波線距減半前形成的十字排列域；圖5為炮線距和檢波線距減半后形成的十字排列域。如圖4和圖5所示，十字排列域的橫線表示炮點形成的炮線，縱線表示檢波點形成的檢波線，從圖5可知，炮線距和檢波線距均縮小為原來的一半，由于原來的炮線之間的距離為500米，檢波線之間的距離為400米，所以，形成的新的炮線距為250米，新的檢線距為200米。

第三：利用五維插值法對每個子共偏移距共方位角域道集進行數據插值，得到加密數據；該加密數據包括原始地震數據和插值后的數據。

圖6為共偏移距共方位角域道集進行數據插值前的數據分布示意圖；圖7為共偏移距共方位角域道集進行數據插值后的數據分布示意圖。圖6中的空白部分為地震數據分選到共偏移距共方位角域道集后形成的空洞，也即，本步驟主要是利用五維插值法對每個子共偏移距共方位角域道集中的空洞進行數據插值。從圖6和圖7可以看出，插值后得到的加密數據與原始地震數據的振幅、頻率、相位特性一致，包括面波的特征也很一致，這表明五維插值法的插值精度很高。

第四：將上述加密數據按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內，得到4N個子共偏移距共方位角域道集；

由于插值后的數據加上原始數據稱為加密數據，對加密數據按新的炮線距和檢波線距重新劃分共偏移距共方位角域單元，再抽取共偏移距共方位角域道集，由于炮線距和檢波線距均縮小為了原來的一倍，所以，原來的一個共偏移距共方位角域單元一化為四。

圖8為原始地震數據對共偏移距共方位角域道集的劃分示意圖；圖9為加密數據對共偏移距共方位角域道集的劃分示意圖。圖10為原始地震數據在共偏移距共方位角域道集上形成的項鏈圖；圖11為加密數據在共偏移距共方位角域道集上形成的項鏈圖。

從圖8可知，圖8上共有180個完整的共偏移距共方位角域道集，個數與原始采集覆蓋次數相同，加密數據按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內，得到720個子共偏移距共方位角域道集，也即，十字排列域形成的某一網格覆蓋次數增大為原來的4倍，具體參見圖9所示。類似的，圖10和圖11分別所示的原始地震數據和加密數據在共偏移距共方位角域道集上形成的項鏈圖也同樣說明了上述結論。

圖12為原始地震數據在共偏移距共方位角域道集上的偏移距范圍。圖13為加密數據分別在相應的四個子共偏移距共方位角域道集上的偏移距范圍。

參照圖12和圖13，由于炮線距和檢波線距均縮小了一倍，因此，偏移距范圍也縮小了一倍，顯示在圖上后可知四個子共偏移距共方位角域道集的顏色更均一。所以，偏移距范圍越小，意味著偏移距分布的一致性越好，也更有利于淺層的成像和高保真處理。

第五：對上述720個子共偏移距共方位角域道集進行疊加前偏移處理，得到壓制采集腳印后的地震數據分析結果。

具體的，對上述720個子共偏移距共方位角域道集重新進行偏移、疊加，再提取振幅切片，便可得到壓制采集腳印后的地震數據分析結果。

圖14為原始地震數據的共偏移距偏移處理成果切片得到的地震數據分析結果，圖15為加密數據的共偏移距偏移處理成果切片得到的地震數據分析結果。從圖14可以看出，雖然對原始地震數據的處理已經施加了常規(guī)壓制采集腳印的手段，但采集腳印還是明顯存在。從圖15可以看出，對加密數據處理后得到的地震數據分析結果顯示，采集腳印幾乎消失，效果非常好。

上述結果表明，本發(fā)明實施例利用共偏移距共方位角域道集的單次覆蓋特性，通過在同一方位角內進行五維插值實現了加密炮線和檢波線的目的，整個技術實現過程高效、快捷、方便，且對于壓制采集腳印、淺層成像、高保真處理都具有實用的價值。

圖16為本發(fā)明提供的采集腳印壓制系統(tǒng)實施例一的結構示意圖。如圖16所示，本發(fā)明實施例一提供的采集腳印壓制系統(tǒng)，包括：

獲取模塊161，用于將待處理的地震數據按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內，得到N個共偏移距共方位角域道集；

其中，共偏移距共方位角域道集是指具有相同偏移距和相同方位角的地震道的集合，N為相同偏移距和相同方位角區(qū)域被覆蓋的次數，N為大于或等于1的整數。

處理模塊162，用于分別將炮線之間的距離和檢波線之間的距離進行減半處理，根據每個共偏移距共方位角域道集中的每個地震道的炮點和/或檢波點的坐標，計算出新增加的炮點和/或檢波點的坐標，每個共偏移距共方位角域道集生成四個子共偏移距共方位角域道集；

插值模塊163，用于分別對四個子共偏移距共方位角域道集進行數據插值，得到加密數據；

加密數據包括原始地震數據和插值后的數據。

上述獲取模塊161，還用于將加密數據按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內，得到4N個子共偏移距共方位角域道集；

上述處理模塊162，還用于對4N個子共偏移距共方位角域道集進行疊加前偏移處理，得到壓制采集腳印后的地震數據分析結果。

本發(fā)明實施例提供的采集腳印壓制系統(tǒng)，可以用于執(zhí)行如圖1所示采集腳印壓制方法實施例的技術方案，其實現原理和技術效果類似，此處不再贅述。

進一步的，在上述實施例一提供的采集腳印壓制系統(tǒng)中，上述處理模塊162，包括：第一處理單元、第二處理單元和第三處理單元；

其中，第一處理單元，用于通過改變共偏移距共方位角域道集中每個地震道的炮點坐標，而不改變共偏移距共方位角域道集中每個地震道的檢波點 坐標，計算出新增加的炮點的坐標；

具體的，該第一處理單元，具體用于采用下述坐標變換公式(1)得到新增加的炮點的坐標：

其中，x_s和y_s分別為共偏移距共方位角域道集中每個地震道的炮點的橫坐標和縱坐標，x和y分別為新增加的炮點的橫坐標和縱坐標，θ為共偏移距共方位角域道集中炮線的方位角，dl_s為共偏移距共方位角域道集中的炮線距。

第二處理單元，用于通過改變共偏移距共方位角域道集中每個地震道的檢波點坐標，而不改變共偏移距共方位角域道集中每個地震道的炮點坐標，計算出新增加的檢波點的坐標；

具體的，該第二處理單元，具體用于采用下述坐標變換公式(2)得到新增加的檢波點的坐標：

其中，x_r和y_r分別為共偏移距共方位角域道集中每個地震道的檢波點的橫坐標和縱坐標，x′和y′分別為新增加的檢波點的橫坐標和縱坐標，α為共偏移距共方位角域道集中檢波線的方位角，dl_r為共偏移距共方位角域道集中的炮線距。

第三處理單元，用于通過改變共偏移距共方位角域道集中每個地震道的炮點坐標和檢波點坐標，計算出新增加的炮點的坐標和檢波點的坐標。

進一步的，在上述實施例一提供的采集腳印壓制系統(tǒng)中，上述插值模塊163，具體用于利用五維插值法對四個子共偏移距共方位角域道集進行數據插值處理，得到加密數據。

本發(fā)明實施例提供的采集腳印壓制系統(tǒng)，通過在炮線之間和檢波線之間內插炮線和檢波線，使得組成十字排列域的炮線和檢波線加密，并利用共偏移距共方位角域道集的單次覆蓋特性，通過在同一方位角內進行五維高精度插值得到了加密數據，將該加密數據按照偏移距和方位角重新分選到加密炮線和加密檢波線組成的十字排列域內，使獲得的子共偏移距共方位角域道集 為原來的四倍，整個技術實現過程高效、快捷、方便，最后對所有的子共偏移距共方位角域道集進行處理，得到的地震數據分析結果，有效的壓制了采集腳印的產生，因此，該技術方案對壓制采集腳印、淺層成像、高保真處理均具有很高的實用價值。

最后應說明的是：以上各實施例僅用以說明本發(fā)明的技術方案，而非對其限制；盡管參照前述各實施例對本發(fā)明進行了詳細的說明，本領域的普通技術人員應當理解：其依然可以對前述各實施例所記載的技術方案進行修改，或者對其中部分或者全部技術特征進行等同替換；而這些修改或者替換，并不使相應技術方案的本質脫離本發(fā)明各實施例技術方案的范圍。