本申請涉及石油地震勘探中的地震數(shù)據(jù)處理技術(shù)領(lǐng)域，特別涉及一種地震采集腳印噪音壓制方法及系統(tǒng)。

背景技術(shù)：

在地震數(shù)據(jù)的采集過程中所使用的觀測系統(tǒng)會對采集到的地震數(shù)據(jù)產(chǎn)生干擾，可以形象地將這種干擾被喻為在地震數(shù)據(jù)中留下了腳印，即采集腳印。這種腳印表現(xiàn)為在疊后地震數(shù)據(jù)的時間切片上出現(xiàn)有規(guī)律的振幅變化，是調(diào)制在地層反射信號上的一種系統(tǒng)性噪聲，它嚴(yán)重地影響了地震屬性反演的精確度。

通常有三種途徑衰減采集腳印噪音，第一種途徑是通過采集系統(tǒng)參數(shù)設(shè)置使不同偏移距的道數(shù)變化最??；第二種途徑是通過疊前處理使要疊加的道集間差異最小；第三種途徑是通過疊后處理壓制采集腳印。地震數(shù)據(jù)處理實例表明，通過優(yōu)化采集參數(shù)以及疊前處理措施后，采集腳印依然發(fā)育。因此，需要對疊后地震數(shù)據(jù)進(jìn)行噪音壓制，以提高地震數(shù)據(jù)信噪比和突出有效反射信號。

目前基于疊后地震勘探數(shù)據(jù)壓制采集腳印噪音的方法主要是頻率波數(shù)(F-Kx-Ky)濾波方法，這種方法是通過二維傅里葉變換將原始地震數(shù)據(jù)由空間域變換至波數(shù)域，得到原始地震數(shù)據(jù)波數(shù)譜，根據(jù)采集腳印周期性特點設(shè)計濾波因子，利用所述濾波因子對原始地震數(shù)據(jù)波譜數(shù)進(jìn)行采集腳印噪音壓制，得到目標(biāo)地震數(shù)據(jù)波譜數(shù)，通過二維傅里葉反變換將所述目標(biāo)地震數(shù)據(jù)波譜數(shù)由波數(shù)域變換至空間域，完成采集腳印噪音壓制。

發(fā)明人發(fā)現(xiàn)現(xiàn)有技術(shù)中至少存在如下問題：

部分采集腳印，例如海上采集腳印噪音，具有隨時間和空間變化的特征，在整個地震數(shù)據(jù)中沒有明顯的周期性分布特點。由于現(xiàn)有技術(shù)中頻率波數(shù)濾波方法是根據(jù)采集腳印周期性特點設(shè)計濾波因子來完成采集腳印噪音壓制的，所以，該方法很難有效識別和壓制所述隨時間和空間變化的采集腳印，對地質(zhì)信號保持能力較差。

技術(shù)實現(xiàn)要素：

本申請實施例的目的是提供一種地震采集腳印噪音壓制方法及系統(tǒng)，以提高自動識別并壓制采集腳印噪音以及保持有效地質(zhì)信號的能力。

為解決上述技術(shù)問題，本申請實施例提供一種地震采集腳印噪音壓制方法及系統(tǒng)是這樣實現(xiàn)的：

一種地震采集腳印噪音壓制方法，包括：

獲取原始疊后地震數(shù)據(jù)；

調(diào)整所述原始疊后地震數(shù)據(jù)的排列方式，得到標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)；

獲取所述標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)中第一時間切片的地震數(shù)據(jù)，對第一時間切片的地震數(shù)據(jù)進(jìn)行拉普拉斯變換，得到拉普拉斯變換地震數(shù)據(jù)；

將所述拉普拉斯變換地震數(shù)據(jù)進(jìn)行二維傅里葉變換，得到拉普拉斯變換地震數(shù)據(jù)波數(shù)譜，根據(jù)所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜，獲取第一時間切片的濾波因子；

將第一時間切片的地震數(shù)據(jù)進(jìn)行二維傅里葉變換，得到第一時間切片的地震數(shù)據(jù)波數(shù)譜；

利用第一時間切片的濾波因子對第一時間切片的地震數(shù)據(jù)波數(shù)譜進(jìn)行采集腳印噪音壓制。

優(yōu)選方案中，在利用所述第一時間切片的濾波因子對第一時間切片的地震數(shù)據(jù)波數(shù)譜進(jìn)行采集腳印噪音壓制之后，所述方法還包括：

將采集腳印噪音壓制后的第一時間切片的地震數(shù)據(jù)波數(shù)譜從波數(shù)域轉(zhuǎn)換到空間域，得到中間疊后地震數(shù)據(jù)；

調(diào)整所述中間疊后地震數(shù)據(jù)的排列方式，得到目標(biāo)疊后地震數(shù)據(jù)。

優(yōu)選方案中，通過矩陣轉(zhuǎn)置調(diào)整所述中間疊后地震數(shù)據(jù)的排列方式。

優(yōu)選方案中，所述根據(jù)所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜，獲取第一時間切片的濾波因子，包括：

根據(jù)所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜，獲取所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中的最大波數(shù)值和最小波數(shù)值；

根據(jù)預(yù)設(shè)規(guī)則以及所述最大波數(shù)值和所述最小波數(shù)值，調(diào)整所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中的波數(shù)值；

將所述調(diào)整后的拉普拉斯變換地震數(shù)據(jù)波數(shù)譜作為第一時間切片的濾波因子。

優(yōu)選方案中，所述根據(jù)預(yù)設(shè)規(guī)則以及所述最大波數(shù)值和所述最小波數(shù)值，調(diào)整所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中的波數(shù)值，包括：

采用下述公式實現(xiàn)：

$A_i=\frac{A_{0\max }-A_{0i}}{A_{0\max }-A_{0\min }}$

公式中，A_i表示所述調(diào)整后的拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中第i個波數(shù)值；A_0max表示所述調(diào)整前的拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中的最大波數(shù)值；A_0min表示所述調(diào)整前的拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中的最小波數(shù)值；A_0i表示所述調(diào)整前的拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中第i個波數(shù)值。

優(yōu)選方案中，所述利用所述第一時間切片的濾波因子對第一時間切片的地震數(shù)據(jù)波數(shù)譜進(jìn)行采集腳印噪音壓制，包括：將第一時間切片的濾波因子與第一時間切片的地震數(shù)據(jù)波數(shù)譜在波數(shù)域相乘。

優(yōu)選方案中，所述地震數(shù)據(jù)包括：時間、振幅值、主測線號和聯(lián)絡(luò)測線號。

優(yōu)選方案中，所述調(diào)整原始疊后地震數(shù)據(jù)的排列方式，包括：調(diào)整所述原始疊后地震數(shù)據(jù)的至少兩個維度的排列方式；

具體地，將所述原始疊后地震數(shù)據(jù)的至少兩個維度中的時間維度設(shè)置為最慢維度；所述原始疊后地震數(shù)據(jù)的至少兩個維度包括：所述時間對應(yīng)的時間維度、所述主測線號對應(yīng)的主測線維度和所述聯(lián)絡(luò)測線號對應(yīng)的聯(lián)絡(luò)測線維度。

優(yōu)選方案中，通過矩陣轉(zhuǎn)置調(diào)整所述原始疊后地震數(shù)據(jù)的排列方式。

一種地震采集腳印噪音壓制系統(tǒng)，包括：原始疊后地震數(shù)據(jù)獲取單元、排列方式調(diào)整單元、拉普拉斯變換單元、濾波因子獲取單元、傅里葉變換單元和噪音壓制單元；其中，

所述原始疊后地震數(shù)據(jù)獲取單元，用于獲取原始疊后地震數(shù)據(jù)；

所述排列方式調(diào)整單元，用于調(diào)整所述原始疊后地震數(shù)據(jù)的排列方式，得到標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)；

所述拉普拉斯變換單元，用于獲取所述標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)中第一時間切片的地震數(shù)據(jù)，對第一時間切片的地震數(shù)據(jù)進(jìn)行拉普拉斯變換，得到拉普拉斯變換地震數(shù)據(jù)；

所述濾波因子獲取單元，用于將所述拉普拉斯變換地震數(shù)據(jù)進(jìn)行二維傅里葉變換，得到拉普拉斯變換地震數(shù)據(jù)波數(shù)譜，根據(jù)所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜，獲取第一時間切片的濾波因子；

所述傅里葉變換單元，用于將第一時間切片的地震數(shù)據(jù)進(jìn)行二維傅里葉變換，得到第一時間切片的地震數(shù)據(jù)波數(shù)譜；

所述噪音壓制單元，用于利用第一時間切片的濾波因子對第一時間切片的地震數(shù)據(jù)波數(shù)譜進(jìn)行采集腳印噪音壓制。

優(yōu)選方案中，所述系統(tǒng)還包括：中間疊后地震數(shù)據(jù)獲取單元和目標(biāo)疊后地震數(shù)據(jù)獲取單元；其中，

所述中間疊后地震數(shù)據(jù)獲取單元，用于將采集腳印噪音壓制后的第一時間切片的地震數(shù)據(jù)波數(shù)譜從波數(shù)域轉(zhuǎn)換到空間域，得到中間疊后地震數(shù)據(jù)；

所述目標(biāo)疊后地震數(shù)據(jù)獲取單元，用于調(diào)整所述中間疊后地震數(shù)據(jù)的排列方式，得到目標(biāo)疊后地震數(shù)據(jù)。

本申請?zhí)峁┝艘环N地震采集腳印噪音壓制方法及系統(tǒng)，基于時間切片地震數(shù)據(jù)的拉普拉斯變換和波數(shù)濾波結(jié)合的方式，所述方法對疊后地震數(shù)據(jù)進(jìn)行采集腳印壓制，能夠根據(jù)海上拖纜地震數(shù)據(jù)采集腳印隨空間和時間變化的特點，自動識別并壓制采集腳印噪音，并且能夠有效保持地質(zhì)信號。

附圖說明

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

圖1是本申請一種地震采集腳印噪音壓制方法實施例的流程圖；

圖2是本申請第1000毫秒處的時間切片的地震數(shù)據(jù)及其波數(shù)譜的示意圖；

圖3是本申請第1000毫秒處的時間切片的地震數(shù)據(jù)拉普拉斯變換后的拉普拉斯變換地震數(shù)據(jù)和拉普拉斯變換地震數(shù)據(jù)波數(shù)譜的示意圖；

圖4是本申請采集腳印噪音壓制后的第1000毫秒處的時間切片的地震數(shù)據(jù)及其波數(shù)譜的示意圖；

圖5是本申請采集腳印噪音壓制前后的第1000毫秒處的時間切片的地震數(shù)據(jù)差及其波數(shù)譜的示意圖；

圖6是本申請一種地震采集腳印噪音壓制系統(tǒng)實施例的組成結(jié)構(gòu)圖。

具體實施方式

本申請實施例提供一種地震采集腳印噪音壓制方法及系統(tǒng)。

為了使本技術(shù)領(lǐng)域的人員更好地理解本申請中的技術(shù)方案，下面將結(jié)合本申請實施例中的附圖，對本申請實施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實施例僅僅是本申請一部分實施例，而不是全部的實施例。基于本申請中的實施例，本領(lǐng)域普通技術(shù)人員在沒有作出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都應(yīng)當(dāng)屬于本申請保護(hù)的范圍。

圖1是本申請一種地震采集腳印噪音壓制方法實施例的流程圖。如圖1所示，所述地震采集腳印噪音壓制方法，包括以下步驟。

步驟S101：獲取原始疊后地震數(shù)據(jù)。

具體地，可以將原始疊后地震數(shù)據(jù)讀入計算機(jī)內(nèi)存。其中，將地表接受到的來自地下同一反射點的地震數(shù)據(jù)經(jīng)過動校正后進(jìn)行疊加，可以有效提高地震數(shù)據(jù)的信噪比，疊加后的地震數(shù)據(jù)，稱為疊后地震數(shù)據(jù)。所述原始疊后地震數(shù)據(jù)可以包括：時間、振幅值、主測線(Inline)號和聯(lián)絡(luò)測線(Xline)號。所述原始疊后地震數(shù)據(jù)為三維地震數(shù)據(jù)，可以包括：時間方向、Inline方向和Xline方向。其中，所述原始疊后地震數(shù)據(jù)中Inline方向和Xline方向相互垂直，時間方向與Inline方向和Xline方向均垂直。

例如，所述原始疊后地震數(shù)據(jù)中的振幅值可以為所述原始疊后地震數(shù)據(jù)中Inline方向上的Inline號5574和Xline方向上的Xline號3221在時間方向上的時間1000毫秒處對應(yīng)的時間樣點的振幅值。

步驟S102：調(diào)整所述原始疊后地震數(shù)據(jù)的排列方式，得到標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)。

具體地，所述原始疊后地震數(shù)據(jù)的維度可以包括：所述時間對應(yīng)的時間維度、所述主測線號對應(yīng)的主測線維度和所述聯(lián)絡(luò)測線號對應(yīng)的聯(lián)絡(luò)測線維度。所述原始疊后地震數(shù)據(jù)是按照Inline維度-Xline維度-時間維度的序列進(jìn)行排列的。由于采集腳印沿時間切片具有條帶狀等規(guī)律性特征，可以通過矩陣轉(zhuǎn)置將所述原始疊后地震數(shù)據(jù)的排列方式調(diào)整為時間維度-Inline維度-Xline維度的序列。即可以將所述時間維度設(shè)置為最慢維度?？梢詫⑺雠帕蟹绞秸{(diào)整后的原始疊后地震數(shù)據(jù)作為標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)。

例如，所述矩陣轉(zhuǎn)置的過程可以采用下述方式實現(xiàn)：

定義矩陣A轉(zhuǎn)置為n×m階矩陣B，滿足B＝a(j，i)，即b(i，j)＝a(j，i)(B的第i行第j列元素是A的第j行第i列元素)，記A^T＝B。二維矩陣轉(zhuǎn)置同理三維數(shù)組的轉(zhuǎn)置，保持一維不變，另外兩維數(shù)組按照二維矩陣轉(zhuǎn)置方式實現(xiàn)轉(zhuǎn)置。

步驟S103：獲取所述標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)中第一時間切片的地震數(shù)據(jù)，對第一時間切片的地震數(shù)據(jù)進(jìn)行拉普拉斯變換，得到拉普拉斯變換地震數(shù)據(jù)。

具體地，可以將所述標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)由時間域轉(zhuǎn)換到時間切片域，可以得到第一時間處對應(yīng)的第一時間切片的地震數(shù)據(jù)?？梢詫Λ@得的第一時間切片的地震數(shù)據(jù)進(jìn)行拉普拉斯變換，得到拉普拉斯變換地震數(shù)據(jù)，可以使所述標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)中的采集腳印噪音在拉普拉斯變換地震數(shù)據(jù)中突顯出來。其中，第一時間切片的地震數(shù)據(jù)可以是所述標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)中任意時間處對應(yīng)的任意時間切片的地震數(shù)據(jù)。

例如，如圖2所示，圖2中左圖為第1000毫秒處的時間切片的地震數(shù)據(jù)。如圖3所示，圖3中左圖為圖2中第1000毫秒處的時間切片的地震數(shù)據(jù)拉普拉斯變換后的拉普拉斯變換地震數(shù)據(jù)。從圖3中左圖中可以看出，第1000毫秒處的時間切片的地震數(shù)據(jù)經(jīng)過拉普拉斯變換以后，條帶狀的采集腳印噪音在所述拉普拉斯變化地震數(shù)據(jù)中突顯出來，從而使采集腳印噪音得到有效表征。

步驟S104：將所述拉普拉斯變換地震數(shù)據(jù)進(jìn)行二維傅里葉變換，得到拉普拉斯變換地震數(shù)據(jù)波數(shù)譜，根據(jù)所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜，獲取第一時間切片的濾波因子。

具體地，將圖3中左圖中拉普拉斯變換地震數(shù)據(jù)進(jìn)行二維傅里葉變換，得到拉普拉斯變換地震數(shù)據(jù)波譜數(shù)。根據(jù)圖3中右圖中所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜，可以獲取所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中的最大波數(shù)值和最小波數(shù)值。根據(jù)預(yù)設(shè)規(guī)則以及最大波數(shù)值和最小波數(shù)值，可以調(diào)整所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中的波數(shù)值。所述預(yù)設(shè)規(guī)則可以采用下述公式實現(xiàn)：

$A_i=\frac{A_{0\max }-A_{0i}}{A_{0\max }-A_{0\min }}$

公式中，A_i表示所述調(diào)整后的拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中第i個波數(shù)值；A_0max表示所述調(diào)整前的拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中的最大波數(shù)值；A_0min表示所述調(diào)整前的拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中的最小波數(shù)值；A_0i表示所述調(diào)整前的拉普拉斯變換地震數(shù)據(jù)波數(shù)譜中第i個波數(shù)值。按照上述公式對所述拉普拉斯變化地震數(shù)據(jù)波數(shù)譜的波數(shù)值進(jìn)行調(diào)整，將調(diào)整后的拉普拉斯變換地震數(shù)據(jù)波數(shù)譜作為第一時間切片的濾波因子。所述得到的第一時間切片的濾波因子是基于第一時間切片的地震數(shù)據(jù)得到的。因此，每一個時間切片的濾波因子可以隨著拉普拉斯變換后識別的每一個時間切片的地震數(shù)據(jù)中的采集腳印噪音的變化而變化，可以對每一個時間切片的地震數(shù)據(jù)中的采集腳印噪音進(jìn)行有效壓制。

例如：如圖3所示，圖3中右圖為圖3中左圖中所述拉普拉斯變換地震數(shù)據(jù)經(jīng)過二維傅里葉變換后，得到的拉普拉斯變換地震數(shù)據(jù)波譜數(shù)。按照所述公式對圖3中右圖中的拉普拉斯變換地震數(shù)據(jù)波譜數(shù)進(jìn)行波數(shù)值調(diào)整，將調(diào)整后的拉普拉斯變換地震數(shù)據(jù)波數(shù)譜作為第1000毫秒處的時間切片的濾波因子。

步驟S105：將第一時間切片的地震數(shù)據(jù)進(jìn)行二維傅里葉變換，得到第一時間切片的地震數(shù)據(jù)波數(shù)譜。

具體地，將所述標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)中第一時間切片的地震數(shù)據(jù)進(jìn)行二維傅里葉變換，得到第一時間切片的地震數(shù)據(jù)波數(shù)譜，以便利用上述步驟中的濾波因子在波數(shù)域?qū)Φ谝粫r間切片的地震數(shù)據(jù)進(jìn)行采集腳印壓制。

例如，將圖2中左圖中第1000毫秒處的時間切片的地震數(shù)據(jù)進(jìn)行二維傅里葉變換，得到第1000毫秒處的時間切片的地震數(shù)據(jù)波數(shù)譜。圖2中右圖為經(jīng)過二維傅里葉變換后的第1000毫秒處的時間切片的地震數(shù)據(jù)的波譜數(shù)。

需要說明的是，步驟S105可以在步驟S103和/或步驟S104之前或之后，本申請對此并不作出限定。

步驟S106：利用第一時間切片的濾波因子對第一時間切片的地震數(shù)據(jù)波數(shù)譜進(jìn)行采集腳印噪音壓制。

具體地，可以將第一時間切片的濾波因子與第一時間切片的地震數(shù)據(jù)波數(shù)譜在波數(shù)域相乘，對第一時間切片的地震數(shù)據(jù)波數(shù)譜進(jìn)行采集腳印壓制，得到所述采集腳印壓制后的第一時間切片的地震數(shù)據(jù)波數(shù)譜。

例如，圖4是本申請采集腳印噪音壓制后的第1000毫秒處的時間切片的地震數(shù)據(jù)及其波數(shù)譜的示意圖。如圖4所示，圖4中左圖為采集腳印噪音壓制后的第1000毫秒處的時間切片的地震數(shù)據(jù)，圖4中右圖為經(jīng)過二維傅里葉變換后的采集腳印噪音壓制后的第1000毫秒處的時間切片的地震數(shù)據(jù)波數(shù)譜。對比圖2中左圖和圖4中左圖，可以看出，條帶狀的采集腳印噪音可以得到有效壓制。從圖4中右圖中可以看出，低、高波數(shù)類采集腳印噪音均可以得到壓制。

圖5是本申請采集腳印噪音壓制前后的第1000毫秒處的時間切片的地震數(shù)據(jù)差及其波數(shù)譜的示意圖。如圖5所示，圖5中左圖為采集腳印噪音壓制前后的第1000毫秒處的時間切片的地震數(shù)據(jù)差，圖5中右圖為經(jīng)過二維傅里葉變換后的采集腳印噪音壓制前后的第1000毫秒處的時間切片的地震數(shù)據(jù)差波數(shù)譜。從圖5可以看出，采集腳印壓制前后地質(zhì)信息得到有效保持。

在另一個實施方式中，還可以將采集腳印噪音壓制后的第一時間切片的地震數(shù)據(jù)波數(shù)譜從波數(shù)域轉(zhuǎn)換到空間域，得到中間疊后地震數(shù)據(jù)?？梢哉{(diào)整所述中間疊后地震數(shù)據(jù)的排列方式，得到目標(biāo)疊后地震數(shù)據(jù)。

具體地，可以采用上述實施方式對所述標(biāo)準(zhǔn)疊后地震數(shù)據(jù)中除第1000毫秒處的時間切片以外其他時間處對應(yīng)的時間切片的地震數(shù)據(jù)進(jìn)行采集腳印壓制，得到采集腳印噪音壓制后的所述標(biāo)準(zhǔn)疊后地震數(shù)據(jù)中所有時間切片的地震數(shù)據(jù)波數(shù)譜。采用二維傅里葉反變換，可以將所述采集腳印噪音壓制后的所述標(biāo)準(zhǔn)疊后地震數(shù)據(jù)中所有時間切片的地震數(shù)據(jù)波數(shù)譜從波數(shù)域轉(zhuǎn)換到時間切片域，然后將其再由時間切片域轉(zhuǎn)換到時間域，得到中間疊后地震數(shù)據(jù)?？梢酝ㄟ^矩陣轉(zhuǎn)置將所述中間疊后地震數(shù)據(jù)的排列方式調(diào)整為Inline維度-Xline維度-時間維度的序列，得到目標(biāo)疊后地震數(shù)據(jù)。從而實現(xiàn)了整個原始疊后地震數(shù)據(jù)的采集腳印噪聲壓制。其中，這一個步驟中矩陣轉(zhuǎn)置的方法與步驟S102中矩陣轉(zhuǎn)置的方法相同。

所述地震采集腳印噪音壓制方法實施例，基于時間切片地震數(shù)據(jù)的拉普拉斯變換和波數(shù)濾波結(jié)合的方式，對疊后地震數(shù)據(jù)進(jìn)行采集腳印壓制，能夠根據(jù)海上拖纜地震數(shù)據(jù)采集腳印隨空間和時間變化的特點，自動識別并壓制采集腳印噪音，并且能夠有效保持地質(zhì)信號。

圖6是本申請一種地震采集腳印噪音壓制系統(tǒng)實施例的組成結(jié)構(gòu)圖。如圖6所示，所述地震采集腳印噪音壓制系統(tǒng)可以包括：原始疊后地震數(shù)據(jù)獲取單元100、排列方式調(diào)整單元200、拉普拉斯變換單元300、濾波因子獲取單元400、傅里葉變換單元500和噪音壓制單元600。

所述原始疊后地震數(shù)據(jù)獲取單元100，可以用于獲取原始疊后地震數(shù)據(jù)。

所述排列方式調(diào)整單元200，可以用于調(diào)整所述原始疊后地震數(shù)據(jù)的排列方式，得到標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)。

所述拉普拉斯變換單元300，可以用于獲取所述標(biāo)準(zhǔn)原始疊后地震數(shù)據(jù)中第一時間切片的地震數(shù)據(jù)?？梢詫Φ谝粫r間切片的地震數(shù)據(jù)進(jìn)行拉普拉斯變換，得到拉普拉斯變換地震數(shù)據(jù)。

所述濾波因子獲取單元400，可以用于將所述拉普拉斯變換地震數(shù)據(jù)進(jìn)行二維傅里葉變換，得到拉普拉斯變換地震數(shù)據(jù)波數(shù)譜。根據(jù)所述拉普拉斯變換地震數(shù)據(jù)波數(shù)譜，可以獲取第一時間切片的濾波因子。

所述第二傅里葉變換單元500，可以用于將第一時間切片的地震數(shù)據(jù)進(jìn)行二維傅里葉變換，得到第一時間切片的地震數(shù)據(jù)波數(shù)譜。

所述噪音壓制單元600，可以用于利用第一時間切片的濾波因子對第一時間切片的地震數(shù)據(jù)波數(shù)譜進(jìn)行采集腳印噪音壓制。

在另一個實施方式中，所述地震采集腳印噪音壓制系統(tǒng)還可以包括：中間疊后地震數(shù)據(jù)獲取單元700和目標(biāo)疊后地震數(shù)據(jù)獲取單元800。

所述中間疊后地震數(shù)據(jù)獲取單元700，可以用于將采集腳印噪音壓制后的第一時間切片的地震數(shù)據(jù)波數(shù)譜從波數(shù)域轉(zhuǎn)換到空間域，得到中間疊后地震數(shù)據(jù)。

所述目標(biāo)疊后地震數(shù)據(jù)獲取單元800，可以用于調(diào)整所述中間疊后地震數(shù)據(jù)的排列方式，得到目標(biāo)疊后地震數(shù)據(jù)。

所述地震采集腳印噪音壓制系統(tǒng)實施例與所述地震采集腳印噪音壓制方法實施例相對應(yīng)，可以實現(xiàn)對疊后地震數(shù)據(jù)進(jìn)行采集腳印壓制，能夠根據(jù)海上拖纜地震數(shù)據(jù)采集腳印隨空間和時間變化的特點，自動識別并壓制采集腳印噪音，并且能夠有效保持地質(zhì)信號。

在20世紀(jì)90年代，對于一個技術(shù)的改進(jìn)可以很明顯地區(qū)分是硬件上的改進(jìn)(例如，對二極管、晶體管、開關(guān)等電路結(jié)構(gòu)的改進(jìn))還是軟件上的改進(jìn)(對于方法流程的改進(jìn))。然而，隨著技術(shù)的發(fā)展，當(dāng)今的很多方法流程的改進(jìn)已經(jīng)可以視為硬件電路結(jié)構(gòu)的直接改進(jìn)。設(shè)計人員幾乎都通過將改進(jìn)的方法流程編程到硬件電路中來得到相應(yīng)的硬件電路結(jié)構(gòu)。因此，不能說一個方法流程的改進(jìn)就不能用硬件實體模塊來實現(xiàn)。例如，可編程邏輯器件(Programmable Logic Device,PLD)(例如現(xiàn)場可編程門陣列(Field Programmable Gate Array，F(xiàn)PGA))就是這樣一種集成電路，其邏輯功能由用戶對器件編程來確定。由設(shè)計人員自行編程來把一個數(shù)字系統(tǒng)“集成”在一片PLD上，而不需要請芯片制造廠商來設(shè)計和制作專用的集成電路芯片2。而且，如今，取代手工地制作集成電路芯片，這種編程也多半改用“邏輯編譯器(logic compiler)”軟件來實現(xiàn)，它與程序開發(fā)撰寫時所用的軟件編譯器相類似，而要編譯之前的原始代碼也得用特定的編程語言來撰寫，此稱之為硬件描述語言(Hardware Description Language，HDL)，而HDL也并非僅有一種，而是有許多種，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby Hardware Description Language)等，目前最普遍使用的是VHDL(Very-High-Speed Integrated Circuit Hardware Description Language)與Verilog2。本領(lǐng)域技術(shù)人員也應(yīng)該清楚，只需要將方法流程用上述幾種硬件描述語言稍作邏輯編程并編程到集成電路中，就可以很容易得到實現(xiàn)該邏輯方法流程的硬件電路。

控制器可以按任何適當(dāng)?shù)姆绞綄崿F(xiàn)，例如，控制器可以采取例如微處理器或處理器以及存儲可由該(微)處理器執(zhí)行的計算機(jī)可讀程序代碼(例如軟件或固件)的計算機(jī)可讀介質(zhì)、邏輯門、開關(guān)、專用集成電路(Application Specific Integrated Circuit，ASIC)、可編程邏輯控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320，存儲器控制器還可以被實現(xiàn)為存儲器的控制邏輯的一部分。

本領(lǐng)域技術(shù)人員也知道，除了以純計算機(jī)可讀程序代碼方式實現(xiàn)控制器以外，完全可以通過將方法步驟進(jìn)行邏輯編程來使得控制器以邏輯門、開關(guān)、專用集成電路、可編程邏輯控制器和嵌入微控制器等的形式來實現(xiàn)相同功能。因此這種控制器可以被認(rèn)為是一種硬件部件，而對其內(nèi)包括的用于實現(xiàn)各種功能的裝置也可以視為硬件部件內(nèi)的結(jié)構(gòu)?；蛘呱踔?，可以將用于實現(xiàn)各種功能的裝置視為既可以是實現(xiàn)方法的軟件模塊又可以是硬件部件內(nèi)的結(jié)構(gòu)。

上述實施例闡明的系統(tǒng)、裝置、模塊或單元，具體可以由計算機(jī)芯片或?qū)嶓w實現(xiàn)，或者由具有某種功能的產(chǎn)品來實現(xiàn)。

為了描述的方便，描述以上裝置時以功能分為各種單元分別描述。當(dāng)然，在實施本申請時可以把各單元的功能在同一個或多個軟件和/或硬件中實現(xiàn)。

通過以上的實施方式的描述可知，本領(lǐng)域的技術(shù)人員可以清楚地了解到本申請可借助軟件加必需的通用硬件平臺的方式來實現(xiàn)。基于這樣的理解，本申請的技術(shù)方案本質(zhì)上或者說對現(xiàn)有技術(shù)做出貢獻(xiàn)的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來，在一個典型的配置中，計算設(shè)備包括一個或多個處理器(CPU)、輸入/輸出接口、網(wǎng)絡(luò)接口和內(nèi)存。該計算機(jī)軟件產(chǎn)品可以包括若干指令用以使得一臺計算機(jī)設(shè)備(可以是個人計算機(jī)，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本申請各個實施例或者實施例的某些部分所述的方法。該計算機(jī)軟件產(chǎn)品可以存儲在內(nèi)存中，內(nèi)存可能包括計算機(jī)可讀介質(zhì)中的非永久性存儲器，隨機(jī)存取存儲器(RAM)和/或非易失性內(nèi)存等形式，如只讀存儲器(ROM)或閃存(flash RAM)。內(nèi)存是計算機(jī)可讀介質(zhì)的示例。計算機(jī)可讀介質(zhì)包括永久性和非永久性、可移動和非可移動媒體可以由任何方法或技術(shù)來實現(xiàn)信息存儲。信息可以是計算機(jī)可讀指令、數(shù)據(jù)結(jié)構(gòu)、程序的模塊或其他數(shù)據(jù)。計算機(jī)的存儲介質(zhì)的例子包括，但不限于相變內(nèi)存(PRAM)、靜態(tài)隨機(jī)存取存儲器(SRAM)、動態(tài)隨機(jī)存取存儲器(DRAM)、其他類型的隨機(jī)存取存儲器(RAM)、只讀存儲器(ROM)、電可擦除可編程只讀存儲器(EEPROM)、快閃記憶體或其他內(nèi)存技術(shù)、只讀光盤只讀存儲器(CD-ROM)、數(shù)字多功能光盤(DVD)或其他光學(xué)存儲、磁盒式磁帶，磁帶磁磁盤存儲或其他磁性存儲設(shè)備或任何其他非傳輸介質(zhì)，可用于存儲可以被計算設(shè)備訪問的信息。按照本文中的界定，計算機(jī)可讀介質(zhì)不包括短暫電腦可讀媒體(transitory media)，如調(diào)制的數(shù)據(jù)信號和載波。

本說明書中的各個實施例均采用遞進(jìn)的方式描述，各個實施例之間相同相似的部分互相參見即可，每個實施例重點說明的都是與其他實施例的不同之處。尤其，對于系統(tǒng)實施例而言，由于其基本相似于方法實施例，所以描述的比較簡單，相關(guān)之處參見方法實施例的部分說明即可。

本申請可用于眾多通用或?qū)Ｓ玫挠嬎銠C(jī)系統(tǒng)環(huán)境或配置中。例如：個人計算機(jī)、服務(wù)器計算機(jī)、手持設(shè)備或便攜式設(shè)備、平板型設(shè)備、多處理器系統(tǒng)、基于微處理器的系統(tǒng)、置頂盒、可編程的消費電子設(shè)備、網(wǎng)絡(luò)PC、小型計算機(jī)、大型計算機(jī)、包括以上任何系統(tǒng)或設(shè)備的分布式計算環(huán)境等等。

本申請可以在由計算機(jī)執(zhí)行的計算機(jī)可執(zhí)行指令的一般上下文中描述，例如程序模塊。一般地，程序模塊包括執(zhí)行特定任務(wù)或?qū)崿F(xiàn)特定抽象數(shù)據(jù)類型的例程、程序、對象、組件、數(shù)據(jù)結(jié)構(gòu)等等。也可以在分布式計算環(huán)境中實踐本申請，在這些分布式計算環(huán)境中，由通過通信網(wǎng)絡(luò)而被連接的遠(yuǎn)程處理設(shè)備來執(zhí)行任務(wù)。在分布式計算環(huán)境中，程序模塊可以位于包括存儲設(shè)備在內(nèi)的本地和遠(yuǎn)程計算機(jī)存儲介質(zhì)中。

雖然通過實施例描繪了本申請，本領(lǐng)域普通技術(shù)人員知道，本申請有許多變形和變化而不脫離本申請的精神，希望所附的權(quán)利要求包括這些變形和變化而不脫離本申請的精神。