本發(fā)明涉及地震數(shù)據(jù)處理技術(shù)領(lǐng)域，特別涉及一種基于地震數(shù)據(jù)的壓制氣泡干擾的方法和裝置。

背景技術(shù)：

氣槍陣列因其性能穩(wěn)定、自動化程度高、成本低等諸多優(yōu)點已被廣泛用作海水等水域地震勘探震源。在使用氣槍陣列勘探海水等水域震源時，通過激發(fā)置于水下一定深度的氣槍陣列，使其釋放高壓氣體，高壓氣體迅速進(jìn)入周圍水中，形成氣泡，經(jīng)過氣泡的首次膨脹產(chǎn)生可在水中傳播的子波脈沖，子波脈沖到達(dá)海底時，經(jīng)海水等水域的各地質(zhì)分界面透射、折射等過程產(chǎn)生一系列的地震子波，然后，工作人員使用檢波器等設(shè)備采集一系列地震子波的振幅等地震數(shù)據(jù)，從而對海水等水域進(jìn)行勘探。然而，氣泡迅速膨脹到最大限度后便會縮減，當(dāng)縮減到內(nèi)部壓力高于周圍水的壓力時，氣泡會再次經(jīng)歷膨脹、縮減的過程，如此反復(fù)，從而在子波脈沖后形成多個氣泡脈沖，產(chǎn)生次生波，由于彈性衰減作用，子波脈沖后面的氣泡脈沖能量逐漸衰減，形成的次生波不具備參考價值，而且會嚴(yán)重影響初至波的成像質(zhì)量，是一種干擾波。因此，工作人員往往需要對該干擾波進(jìn)行壓制處理。

現(xiàn)有技術(shù)中，終端獲取原始地震數(shù)據(jù)，從原始地震數(shù)據(jù)中提取多個地震子波的地震數(shù)據(jù)，計算該多個地震子波的平均子波的地震數(shù)據(jù)，根據(jù)該平均子波的地震數(shù)據(jù)，確定出反子波的地震數(shù)據(jù)；根據(jù)該反子波的地震數(shù)據(jù)，對原始地震數(shù)據(jù)進(jìn)行濾波，從而對氣泡干擾進(jìn)行壓制。

在實現(xiàn)本發(fā)明的過程中，發(fā)明人發(fā)現(xiàn)現(xiàn)有技術(shù)至少存在以下問題：

在實際生產(chǎn)中,由于水的濃度、溫度、流速等因素的變化，不同地震子波都會有一定程度的改變，從而根據(jù)多個地震子波的地震數(shù)據(jù)確定平均子波的地震數(shù)據(jù)，確定出反子波的地震數(shù)據(jù)時，忽略了不同地震子波之間的差異，壓制氣泡干擾時仍有噪聲殘余，導(dǎo)致壓制氣泡干擾的準(zhǔn)確率較低。

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

為了解決現(xiàn)有技術(shù)的問題，本發(fā)明提供了一種基于地震數(shù)據(jù)的壓制氣泡干擾的方法和裝置。技術(shù)方案如下：

第一方面，提供一種基于地震數(shù)據(jù)的壓制氣泡干擾的方法，所述方法包括：

獲取地震道集合中每個地震道的原始地震數(shù)據(jù)；

從所述每個地震道的原始地震數(shù)據(jù)中選擇炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的第一地震數(shù)據(jù)；

從選擇的第一地震數(shù)據(jù)中確定下行初至波對應(yīng)的第二地震數(shù)據(jù)；

校正所述第二地震數(shù)據(jù)，對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，得到第三地震數(shù)據(jù)；

根據(jù)預(yù)設(shè)切除函數(shù)和所述第三地震數(shù)據(jù)，確定第四地震數(shù)據(jù)；

根據(jù)所述第三地震數(shù)據(jù)和所述第四地震數(shù)據(jù)，確定氣泡壓制因子；

根據(jù)所述氣泡壓制因子，對所述每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，使得所述每個地震道的原始地震數(shù)據(jù)中氣泡干擾被壓制。

在一種可能的設(shè)計中，所述校正所述第二地震數(shù)據(jù)，對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，得到第三地震數(shù)據(jù)，包括：

確定所述每個地震道的下行初至波的旅行時間，通過以下公式一，對所述第二地震數(shù)據(jù)進(jìn)行校正，

公式一：x″i,j(t)＝x'i,j(t-ftj)，

其中，ftj為第j個地震道的下行初至波的旅行時間，x'i,j(t)為所述第二地震數(shù)據(jù)，x″i,j(t)為校正后的第二地震數(shù)據(jù)，t為采集時間，j為地震道的道號，i為炮號；

通過以下公式二，對所述校正后的第二地震數(shù)據(jù)進(jìn)行疊加，得到模型道數(shù)據(jù)，

公式二：

其中，x'i(t)為所述模型道數(shù)據(jù)，x″i,j(t)為所述校正后的第二地震數(shù)據(jù)，t為采集時間，j為地震道的道號，i為炮號，p為炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最小道號，q為所述炮檢距在所述預(yù)設(shè)范圍內(nèi)的地震道的最大道號；

根據(jù)所述模型道數(shù)據(jù)與所述校正后的第二地震數(shù)據(jù)，通過以下公式三，確定所述時差微調(diào)對應(yīng)的互相關(guān)時差，

公式三：

其中，aj(τ)為所述模型道數(shù)據(jù)與所述校正后的第二地震數(shù)據(jù)之間的互相關(guān)函數(shù)，aj為所述互相關(guān)函數(shù)的最大值，dtj為所述時差微調(diào)對應(yīng)的互相關(guān)時差，即所述互相關(guān)函數(shù)的最大值aj對應(yīng)的時差τ1，x'i(t)為所述模型道數(shù)據(jù)，x″i,j(t)為所述校正后的第二地震數(shù)據(jù)，j為地震道的道號，i為炮號，t為采集時間，m為所述互相關(guān)函數(shù)的半支長度；

根據(jù)所述互相關(guān)時差，對所述校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)；

通過以下公式四，對所述時差微調(diào)后的第二地震數(shù)據(jù)進(jìn)行疊加，得到所述第三地震數(shù)據(jù)，

公式四：

其中，x″i(t)為所述第三地震數(shù)據(jù)，dtj為所述時差微調(diào)對應(yīng)的互相關(guān)時差，p為所述炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最小道號，q為所述炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最大道號，x”'i,j(t-dtj)所述時差微調(diào)后的第二地震數(shù)據(jù)，j為地震道的道號，i為炮號，t為采集時間。

在一種可能的設(shè)計中，所述根據(jù)預(yù)設(shè)切除函數(shù)和所述第三地震數(shù)據(jù)，確定第四地震數(shù)據(jù)，包括：

確定所述預(yù)設(shè)切除函數(shù)；

根據(jù)所述預(yù)設(shè)切除函數(shù)和所述第三地震數(shù)據(jù)，通過以下公式六，確定所述第四地震數(shù)據(jù)，

公式六：y″i(t)＝x″i(t)·wi(t)，

其中，x″i(t)為所述第三地震數(shù)據(jù)，wi(t)為所述預(yù)設(shè)切除函數(shù)，y″i(t)為所述第四地震數(shù)據(jù)，i為炮號，t為采集時間。

在一種可能的設(shè)計中，所述根據(jù)所述第三地震數(shù)據(jù)和所述第四地震數(shù)據(jù)，確定氣泡壓制因子，包括：

將所述第三地震數(shù)據(jù)作為輸入數(shù)據(jù)，將所述第四地震數(shù)據(jù)作為期望的輸出數(shù)據(jù)，通過以下公式十，確定濾波結(jié)果表達(dá)式，以及所述期望的輸出數(shù)據(jù)和所述濾波結(jié)果之間的誤差表達(dá)式，計算滿足預(yù)設(shè)條件的氣泡壓制因子，所述預(yù)設(shè)條件為所述誤差表達(dá)式的值最小，

公式十：

其中，y″i(t)為所述第四地震數(shù)據(jù)，即所述期望的輸出數(shù)據(jù)，y(t)為所述濾波結(jié)果，e為所述期望的輸出數(shù)據(jù)和所述濾波結(jié)果之間的誤差表達(dá)式；x″i(t)為所述第三地震數(shù)據(jù)，即所述輸入數(shù)據(jù)，d(t)為所述氣泡壓制因子，i為炮號，t為采集時間。

在一種可能的設(shè)計中，所述根據(jù)所述氣泡壓制因子，對所述每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，包括：

根據(jù)所述氣泡壓制因子，通過以下公式七，對所述每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，得到壓制氣泡干擾后的每個地震道的地震數(shù)據(jù)，

公式七：yi,j(t)＝xi,j(t)*d(t)，

其中，xi,j(t)為所述地震道集合中第j個地震道第i個炮號的原始地震數(shù)據(jù)，d(t)為所述氣泡壓制因子，yi,j(t)為所述壓制氣泡干擾后的每個地震道的地震數(shù)據(jù)，t為采集時間。

第二方面，提供一種基于地震數(shù)據(jù)的壓制氣泡干擾的裝置，所述裝置包括：

獲取模塊，用于獲取地震道集合中每個地震道的原始地震數(shù)據(jù)；

選擇模塊，用于從所述每個地震道的原始地震數(shù)據(jù)中選擇炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的第一地震數(shù)據(jù)；

第一確定模塊，用于從選擇的第一地震數(shù)據(jù)中確定下行初至波對應(yīng)的第二地震數(shù)據(jù)；

校正模塊，用于校正所述第二地震數(shù)據(jù)，對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，得到第三地震數(shù)據(jù)；

第二確定模塊，用于根據(jù)預(yù)設(shè)切除函數(shù)和所述第三地震數(shù)據(jù)，確定第四地震數(shù)據(jù)；

第三確定模塊，用于根據(jù)所述第三地震數(shù)據(jù)和所述第四地震數(shù)據(jù)，確定氣泡壓制因子；

濾波模塊，用于根據(jù)所述氣泡壓制因子，對所述每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，使得所述每個地震道的原始地震數(shù)據(jù)中氣泡干擾被壓制。

在一種可能的設(shè)計中，所述校正模塊，包括：

校正單元，用于確定所述每個地震道的下行初至波的旅行時間，通過以下公式一，對所述第二地震數(shù)據(jù)進(jìn)行校正，

公式一：x″i,j(t)＝x'i,j(t+ftj)，

其中，ftj為第j個地震道的下行初至波的旅行時間，x'i,j(t)為所述第二地震數(shù)據(jù)，x″i,j(t)為校正后的第二地震數(shù)據(jù)，t為采集時間，j為地震道的道號，i為炮號；

第一疊加單元，用于通過以下公式二，對所述校正后的第二地震數(shù)據(jù)進(jìn)行疊加，得到模型道數(shù)據(jù)，

公式二：

其中，x'i(t)為所述模型道數(shù)據(jù)，x″i,j(t)為所述校正后的第二地震數(shù)據(jù)，t為采集時間，j為地震道的道號，i為炮號，p為炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最小道號，q為所述炮檢距在所述預(yù)設(shè)范圍內(nèi)的地震道的最大道號；

第一確定單元，用于根據(jù)所述模型道數(shù)據(jù)與所述校正后的第二地震數(shù)據(jù)，通過以下公式三，確定所述時差微調(diào)對應(yīng)的互相關(guān)時差，

公式三：

其中，aj(τ)為所述模型道數(shù)據(jù)與所述校正后的第二地震數(shù)據(jù)之間的互相關(guān)函數(shù)，aj為所述互相關(guān)函數(shù)的最大值，dtj為所述時差微調(diào)對應(yīng)的互相關(guān)時差，即所述互相關(guān)函數(shù)的最大值aj對應(yīng)的時差τ1，x'i(t)為所述模型道數(shù)據(jù)，x″i,j(t)為所述校正后的第二地震數(shù)據(jù)，j為地震道的道號，i為炮號，t為采集時間，m為所述互相關(guān)函數(shù)的半支長度；

微調(diào)單元，用于根據(jù)所述互相關(guān)時差，對所述校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)；

第二疊加單元，用于通過以下公式四，對所述時差微調(diào)后的第二地震數(shù)據(jù)進(jìn)行疊加，得到所述第三地震數(shù)據(jù)，

公式四：

其中，x″i(t)為所述第三地震數(shù)據(jù)，dtj為所述時差微調(diào)對應(yīng)的互相關(guān)時差，p為所述炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最小道號，q為所述炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最大道號，x”'i,j(t-dtj)所述時差微調(diào)后的第二地震數(shù)據(jù)，j為地震道的道號，i為炮號，t為采集時間。

在一種可能的設(shè)計中，所述第二確定模塊，包括：

第二確定單元，用于確定所述預(yù)設(shè)切除函數(shù)；

第三確定單元，用于根據(jù)所述預(yù)設(shè)切除函數(shù)和所述第三地震數(shù)據(jù)，通過以下公式六，確定所述第四地震數(shù)據(jù)，

公式六：y″i(t)＝x″i(t)·wi(t)，

其中，x″i(t)為所述第三地震數(shù)據(jù)，wi(t)為所述預(yù)設(shè)切除函數(shù)，y″i(t)為所述第四地震數(shù)據(jù)，i為炮號，t為采集時間。

在一種可能的設(shè)計中，所述第三確定模塊，還用于將所述第三地震數(shù)據(jù)作為輸入數(shù)據(jù)，將所述第四地震數(shù)據(jù)作為期望的輸出數(shù)據(jù)，通過以下公式十，確定濾波結(jié)果表達(dá)式，以及所述期望的輸出數(shù)據(jù)和所述濾波結(jié)果之間的誤差表達(dá)式，計算滿足預(yù)設(shè)條件的氣泡壓制因子，所述預(yù)設(shè)條件為所述誤差表達(dá)式的值最小，

公式十：

其中，y″i(t)為所述第四地震數(shù)據(jù)，即所述期望的輸出數(shù)據(jù)，y(t)為所述濾波結(jié)果，e為所述期望的輸出數(shù)據(jù)和所述濾波結(jié)果之間的誤差表達(dá)式；x″i(t)為所述第三地震數(shù)據(jù)，即所述輸入數(shù)據(jù)，d(t)為所述氣泡壓制因子，i為炮號，t為采集時間。

在一種可能的設(shè)計中，所述濾波模塊，還用于根據(jù)所述氣泡壓制因子，通過以下公式七，對所述每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，得到壓制氣泡干擾后的每個地震道的地震數(shù)據(jù)，

公式七：yi,j(t)＝xi,j(t)*d(t)，

其中，xi,j(t)為所述地震道集合中第j個地震道第i個炮號的原始地震數(shù)據(jù)，d(t)為所述氣泡壓制因子，yi,j(t)為所述壓制氣泡干擾后的每個地震道的地震數(shù)據(jù)，t為采集時間。

本發(fā)明實施例中，終端通過獲取炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的第一地震數(shù)據(jù)，從第一地震數(shù)據(jù)中確定下行初至波對應(yīng)的第二地震數(shù)據(jù)，該第二地震數(shù)據(jù)為信噪比較高的地震數(shù)據(jù)；終端校正該第二地震數(shù)據(jù)，并對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，得到消除時差和誤差的第三地震數(shù)據(jù)。從而終端可以根據(jù)預(yù)設(shè)切除函數(shù)和該第三地震數(shù)據(jù)，確定消除氣泡干擾的第四地震數(shù)據(jù)；然后，終端根據(jù)該第三地震數(shù)據(jù)和該第四地震數(shù)據(jù)，確定氣泡壓制因子，終端根據(jù)該氣泡壓制因子，對每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，從而使得每個地震道的原始地震數(shù)據(jù)中氣泡干擾被壓制，終端通過對每個下行初至波對應(yīng)的第二地震數(shù)據(jù)進(jìn)行校正、時差微調(diào)，根據(jù)得到的第三地震數(shù)據(jù)以及消除氣泡干擾的第四地震數(shù)據(jù)，計算出氣泡壓制因子，即，計算該氣泡壓制因子時綜合考慮到了每個下行初至波的第二地震數(shù)據(jù)，通過該氣泡壓制因子對原始地震數(shù)據(jù)進(jìn)行濾波，從而提高了對原始地震數(shù)據(jù)進(jìn)行壓制氣泡干擾的準(zhǔn)確性。

附圖說明

圖1是本發(fā)明實施例提供的一種基于地震數(shù)據(jù)的壓制氣泡干擾的方法流程圖；

圖2是本發(fā)明實施例提供的一種基于地震數(shù)據(jù)的壓制氣泡干擾的方法流程圖；

圖3是本發(fā)明實施例提供的一種切除函數(shù)示意圖；

圖4是本發(fā)明實施例提供的一種基于地震數(shù)據(jù)的壓制氣泡干擾的裝置結(jié)構(gòu)示意圖。

具體實施方式

為使本發(fā)明的目的、技術(shù)方案和優(yōu)點更加清楚，下面將結(jié)合附圖對本發(fā)明實施方式作進(jìn)一步地詳細(xì)描述。

本發(fā)明實施例提供了一種基于地震數(shù)據(jù)的壓制氣泡干擾的方法，該方法的執(zhí)行主體可以為終端，也可以為具備處理數(shù)據(jù)功能的處理器、處理模塊等，該處理器或者處理模塊可以安裝在終端中，如圖1所示，該方法包括：

步驟101：獲取地震道集合中每個地震道的原始地震數(shù)據(jù)。

步驟102：從每個地震道的原始地震數(shù)據(jù)中選擇炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的第一地震數(shù)據(jù)。

步驟103：從選擇的第一地震數(shù)據(jù)中確定下行初至波對應(yīng)的第二地震數(shù)據(jù)。

步驟104：校正該第二地震數(shù)據(jù)，對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，得到第三地震數(shù)據(jù)。

步驟105：根據(jù)預(yù)設(shè)切除函數(shù)和該第三地震數(shù)據(jù)，確定第四地震數(shù)據(jù)。

步驟106：根據(jù)該第三地震數(shù)據(jù)和該第四地震數(shù)據(jù)，確定氣泡壓制因子。

步驟107：根據(jù)該氣泡壓制因子，對每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，使得每個地震道的原始地震數(shù)據(jù)中氣泡干擾被壓制。

在一種可能的設(shè)計中，校正該第二地震數(shù)據(jù)，對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，得到第三地震數(shù)據(jù)，包括：

確定每個地震道的下行初至波的旅行時間，通過以下公式一，對該第二地震數(shù)據(jù)進(jìn)行校正，

公式一：x″i,j(t)＝x'i,j(t-ftj)，

其中，ftj為第j個地震道的下行初至波的旅行時間，x'i,j(t)為該第二地震數(shù)據(jù)，x″i,j(t)為該校正后的第二地震數(shù)據(jù)，t為采集時間，j為地震道的道號，i為炮號；

通過以下公式二，對該校正后的第二地震數(shù)據(jù)進(jìn)行疊加，得到模型道數(shù)據(jù)，

公式二：

其中，x'i(t)為該模型道數(shù)據(jù)，x″i,j(t)為該校正后的第二地震數(shù)據(jù)，t為采集時間，j為地震道的道號，i為炮號，p為炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最小道號，q為炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最大道號；

根據(jù)該模型道數(shù)據(jù)與該校正后的第二地震數(shù)據(jù)，通過以下公式三，確定該時差微調(diào)對應(yīng)的互相關(guān)時差，

公式三：

其中，aj(τ)為該模型道數(shù)據(jù)與該校正后的第二地震數(shù)據(jù)之間的互相關(guān)函數(shù)，aj為該互相關(guān)函數(shù)的最大值，dtj為該時差微調(diào)對應(yīng)的互相關(guān)時差，即該互相關(guān)函數(shù)的最大值aj對應(yīng)的時差τ1，x'i(t)為該模型道數(shù)據(jù)，x″i,j(t)為該校正后的第二地震數(shù)據(jù)，j為地震道的道號，i為炮號，t為采集時間，m為互相關(guān)函數(shù)的半支長度；

根據(jù)該互相關(guān)時差，對該校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)；

通過以下公式四，對該時差微調(diào)后的第二地震數(shù)據(jù)進(jìn)行疊加，得到該第三地震數(shù)據(jù)，

公式四：

其中，x″i(t)為該第三地震數(shù)據(jù)，dtj為該時差微調(diào)對應(yīng)的互相關(guān)時差，p為該炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最小道號，q為該炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最大道號，x”'i,j(t-dtj)該時差微調(diào)后的第二地震數(shù)據(jù)，j為地震道的道號，i為炮號，t為采集時間。

在一種可能的設(shè)計中，根據(jù)預(yù)設(shè)切除函數(shù)和該第三地震數(shù)據(jù)，確定第四地震數(shù)據(jù)，包括：

確定該預(yù)設(shè)切除函數(shù)；

根據(jù)該預(yù)設(shè)切除函數(shù)和該第三地震數(shù)據(jù)，通過以下公式六，確定該第四地震數(shù)據(jù)，

公式六：y″i(t)＝x″i(t)·wi(t)，

其中，x″i(t)為該第三地震數(shù)據(jù)，wi(t)為該預(yù)設(shè)切除函數(shù)，y″i(t)為該第四地震數(shù)據(jù)，i為炮號，t為采集時間。

在一種可能的設(shè)計中，根據(jù)該第三地震數(shù)據(jù)和該第四地震數(shù)據(jù)，確定氣泡壓制因子，包括：

將該第三地震數(shù)據(jù)作為輸入數(shù)據(jù)，將該第四地震數(shù)據(jù)作為期望的輸出數(shù)據(jù)，通過以下公式十，確定濾波結(jié)果表達(dá)式，以及該期望的輸出數(shù)據(jù)和該濾波結(jié)果之間的誤差表達(dá)式，計算滿足預(yù)設(shè)條件的氣泡壓制因子，該預(yù)設(shè)條件為該誤差表達(dá)式的值最小，

公式十：

其中，y″i(t)為該第四地震數(shù)據(jù)，即期望的輸出數(shù)據(jù)，y(t)為該濾波結(jié)果，e為該期望的輸出數(shù)據(jù)和該濾波結(jié)果之間的誤差表達(dá)式；x″i(t)為該第三地震數(shù)據(jù)，即輸入數(shù)據(jù)，d(t)為該氣泡壓制因子，i為炮號，t為采集時間。

在一種可能的設(shè)計中，根據(jù)該氣泡壓制因子，對每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，包括：

根據(jù)該氣泡壓制因子，通過以下公式七，對每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，得到壓制氣泡干擾后的每個地震道的地震數(shù)據(jù)，

公式七：yi,j(t)＝xi,j(t)*d(t)，

其中，xi,j(t)為該地震道集合中第j個地震道第i個炮號的原始地震數(shù)據(jù)，d(t)為該氣泡壓制因子，yi,j(t)為該壓制氣泡干擾后的每個地震道的地震數(shù)據(jù)，t為采集時間。

本發(fā)明實施例中，終端通過獲取炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的第一地震數(shù)據(jù)，從第一地震數(shù)據(jù)中確定下行初至波對應(yīng)的第二地震數(shù)據(jù)，該第二地震數(shù)據(jù)為信噪比較高的地震數(shù)據(jù)；終端校正該第二地震數(shù)據(jù)，并對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，得到消除時差和誤差的第三地震數(shù)據(jù)。從而終端可以根據(jù)預(yù)設(shè)切除函數(shù)和該第三地震數(shù)據(jù)，確定消除氣泡干擾的第四地震數(shù)據(jù)；然后，終端根據(jù)該第三地震數(shù)據(jù)和該第四地震數(shù)據(jù)，確定氣泡壓制因子，終端根據(jù)該氣泡壓制因子，對每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，從而使得每個地震道的原始地震數(shù)據(jù)中氣泡干擾被壓制，終端通過對每個下行初至波對應(yīng)的第二地震數(shù)據(jù)進(jìn)行校正、時差微調(diào)，根據(jù)得到的第三地震數(shù)據(jù)以及消除氣泡干擾的第四地震數(shù)據(jù)，計算出氣泡壓制因子，即，計算該氣泡壓制因子時綜合考慮到了每個下行初至波的第二地震數(shù)據(jù)，通過該氣泡壓制因子對原始地震數(shù)據(jù)進(jìn)行濾波，從而提高了對原始地震數(shù)據(jù)進(jìn)行壓制氣泡干擾的準(zhǔn)確性。

本發(fā)明實施例提供了一種基于地震數(shù)據(jù)的壓制氣泡干擾的方法，該方法的執(zhí)行主體可以為終端，也可以為具備處理數(shù)據(jù)功能的處理器、處理模塊等，該處理器或者處理模塊可以安裝在終端中，如圖2所示，該方法包括：

步驟201：終端獲取地震道集合中每個地震道的原始地震數(shù)據(jù)。

本發(fā)明實施例中，工作人員通過激發(fā)置于水下一定深度的氣槍陣列，從而在水中產(chǎn)生包括子波脈沖和后續(xù)產(chǎn)生的次生脈沖，子波脈沖經(jīng)水底的各地質(zhì)分界面透射、折射等過程，從而產(chǎn)生地震信號。工作人員使用多個采集設(shè)備接收地震信號的地震數(shù)據(jù)；在采集過程中，采集設(shè)備還有可能采集到由次生脈沖產(chǎn)生的氣泡干擾等噪聲信號，這些噪聲信號屬于干擾數(shù)據(jù)，因此，需要對這些干擾信號進(jìn)行壓制。

其中，多個采集設(shè)備按次序排列，每次激發(fā)氣槍陣列時，將每個采集設(shè)備采集的地震信號的地震數(shù)據(jù)稱為一個地震道的地震數(shù)據(jù)。

本步驟中，終端獲取采集設(shè)備采集的每個地震道的地震信號的地震數(shù)據(jù)，作為該地震道的原始地震數(shù)據(jù)。

其中，采集設(shè)備采集的地震道的原始地震數(shù)據(jù)可以根據(jù)用戶需要設(shè)置并更改，本發(fā)明實施例對地震道的原始地震數(shù)據(jù)不做具體限定，本發(fā)明實施例對地震道的原始地震數(shù)據(jù)進(jìn)行時間域處理。地震道的原始地震數(shù)據(jù)可以為地震信號的振幅。

該采集設(shè)備可以為檢波器，終端與檢波器相連接，從而終端直接獲取檢波器采集的每個地震道的原始地震數(shù)據(jù)。

步驟202：終端從每個地震道的原始地震數(shù)據(jù)中選擇炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的第一地震數(shù)據(jù)。

本發(fā)明實施例中，實際作業(yè)時，工作人員持續(xù)激發(fā)氣槍陣列，氣槍陣列周期性的向周圍的水中發(fā)射高壓氣體，我們將氣槍陣列每次向水中發(fā)射高壓氣體的行為稱為一炮，將氣槍陣列在水中的位置稱為炮點，將采集地震信號的原始地震數(shù)據(jù)的采集設(shè)備所在的位置稱為檢波點，由于實際作業(yè)時多個采集設(shè)備按次序排列同時采集地震信號的原始地震數(shù)據(jù)，相鄰兩個采集設(shè)備之間相距預(yù)設(shè)距離，因此，多個檢波點按次序排列、兩兩之間相距預(yù)設(shè)距離。

炮檢距即檢波點與炮點之間的距離，每個地震道的炮檢距即檢測到該地震道的原始地震數(shù)據(jù)的檢波點與炮點之間的距離。

由于地震道的炮檢距越大，采集的地震道的原始地震數(shù)據(jù)噪聲越大、信噪比越低，因此，本步驟中，終端從每個地震道的原始地震數(shù)據(jù)中選擇炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的原始地震數(shù)據(jù)，將選擇的地震道的原始地震數(shù)據(jù)作為第一地震數(shù)據(jù)。

具體的，每個地震道的原始地震數(shù)據(jù)可以表示為xi,j(t)，i＝1，2，3，…n，其中i為炮號，即氣槍陣列第i次向周圍水中發(fā)射高壓氣體；j＝1,2,3,…m，j為道號，即第j個地震道，也即檢測到該地震道的原始地震數(shù)據(jù)的檢波器的編號為j；t為采集時間，即采集設(shè)備采集到地震道的原始地震數(shù)據(jù)的時間；其中，n為所有地震道的原始地震數(shù)據(jù)中的最大炮號，m為所有地震道的原始地震數(shù)據(jù)中最大的道號。

因此，每個地震道的炮檢距可用以下公式八表示為：

公式八：

其中，offj為道號為j的地震道的炮檢距。offmin為最小炮檢距，doff為道間距，即，相鄰兩個地震道之間的預(yù)設(shè)距離，也即采集到該相鄰兩個地震道的原始地震數(shù)據(jù)的相鄰兩個采集設(shè)備之間的預(yù)設(shè)距離，offmax為最大炮檢距。

其中，炮檢距的預(yù)設(shè)范圍可以根據(jù)用戶需要設(shè)置并更改，本發(fā)明實施例對此不作具體限定。例如，該預(yù)設(shè)范圍可以為800m～1000m。

例如，若該炮檢距的預(yù)設(shè)范圍為offp～offq，其中，offp、offq需滿足：offmin≤offp≤offq,offp≤offq≤offmax。則預(yù)設(shè)范圍內(nèi)的最小炮檢距對應(yīng)的地震道的道號最小，為p；最大炮檢距對應(yīng)的地震道的道號最大，為q。

因此，地震道的第一地震數(shù)據(jù)為道號在p、q之間的地震道的原始地震數(shù)據(jù)。

其中，道間距doff可以根據(jù)用戶需要設(shè)置并更改，本發(fā)明實施例對此不作具體限定。

步驟203：終端從選擇的第一地震數(shù)據(jù)中確定下行初至波對應(yīng)的第二地震數(shù)據(jù)。

本發(fā)明實施例中，該第一地震數(shù)據(jù)包括從炮點直接傳播到達(dá)檢波點的地震信號，即下行初至波的地震數(shù)據(jù)；第一地震數(shù)據(jù)還包括經(jīng)水底地質(zhì)分界面一次反射的地震信號的地震數(shù)據(jù)，即下行一次反射的有效波地震數(shù)據(jù)；第一地震數(shù)據(jù)還包括經(jīng)地質(zhì)分界面反射到達(dá)水面后又被水面反射的多次反射的地震數(shù)據(jù)；其中，經(jīng)多次反射的地震信號的地震數(shù)據(jù)噪聲較大，沒有經(jīng)過反射的下行初至波噪聲較小，可以更準(zhǔn)確的反映出每一炮激發(fā)的地震信號的狀況。

本步驟中，終端從第一地震數(shù)據(jù)中選擇下行初至波的地震數(shù)據(jù)，選擇的下行初至波的地震數(shù)據(jù)作為第二地震數(shù)據(jù)。

具體的，第一地震數(shù)據(jù)可以表示為xi,j(t)，對應(yīng)的，可以將第二地震數(shù)據(jù)表示為x'i,j(t)。

步驟204：終端校正該第二地震數(shù)據(jù)，對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，得到第三地震數(shù)據(jù)。

終端獲取第二地震數(shù)據(jù)后，還需對該第二地震數(shù)據(jù)進(jìn)行校正，以消除不同采集設(shè)備采集的不同地震道之間的時差；以及對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，以消除因海水速度等外部因素引起的誤差。

具體的，本步驟可以通過以下步驟2041-2045實現(xiàn)。

步驟2041：終端確定每個地震道的下行初至波的旅行時間，通過以下公式一，對該第二地震數(shù)據(jù)進(jìn)行校正，

公式一：x″i,j(t)＝x'i,j(t-ftj)，

其中，ftj為第j個地震道的下行初至波的旅行時間，x'i,j(t)為該第二地震數(shù)據(jù)，x″i,j(t)為該校正后的第二地震數(shù)據(jù)，t為采集時間，j為地震道的道號，i為炮號。

本步驟中，下行初至波的旅行時間為下行初至波從炮點傳播到檢波點所用的時間，終端確定每個地震道的下行初至波的旅行時間的步驟可以為：終端獲取海水速度、檢波點水深、炮點水深以及該第二地震數(shù)據(jù)對應(yīng)的每個地震道的炮檢距，通過以下公式五，確定每個地震道的下行初至波的旅行時間，

公式五：

其中，vwater為海水速度，offj為第j個地震道對應(yīng)的炮檢距，dep為檢波點水深depr與炮點水深deps的差，ftj為該第j個地震道的下行初至波的旅行時間。

終端對第二地震數(shù)據(jù)進(jìn)行校正后，將該校正后的第二地震數(shù)據(jù)初步疊加為一個初步的模型道數(shù)據(jù)。

步驟2042：終端通過以下公式二，對該校正后的第二地震數(shù)據(jù)進(jìn)行疊加，得到模型道數(shù)據(jù)，

公式二：

其中，x'i(t)為該模型道數(shù)據(jù)，x″i,j(t)為該校正后的第二地震數(shù)據(jù)，t為采集時間，j為地震道的道號，i為炮號，p為炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最小道號，q為炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最大道號。

需要說明的是，本步驟中通過對每一炮對應(yīng)的多個地震道的下行初至波地震數(shù)據(jù)進(jìn)行疊加，從而得到每一炮對應(yīng)的疊加的地震數(shù)據(jù)，將每一炮對應(yīng)的疊加地震數(shù)據(jù)作為模型道數(shù)據(jù)。

步驟2043：終端根據(jù)該模型道數(shù)據(jù)與該校正后的第二地震數(shù)據(jù)，通過以下公式三，確定該時差微調(diào)對應(yīng)的互相關(guān)時差，

公式三：

其中，aj(τ)為該模型道數(shù)據(jù)與該校正后的第二地震數(shù)據(jù)之間的互相關(guān)函數(shù)，aj為該互相關(guān)函數(shù)的最大值，dtj為該時差微調(diào)對應(yīng)的互相關(guān)時差，即該互相關(guān)函數(shù)的最大值aj對應(yīng)的時差τ1，x'i(t)為該模型道數(shù)據(jù)，x″i,j(t)為該校正后的第二地震數(shù)據(jù)，j為地震道的道號，i為炮號，t為采集時間，m為互相關(guān)函數(shù)的半支長度。

本步驟中，終端通過模型道數(shù)據(jù)與校正后的第二地震數(shù)據(jù)之間的互相關(guān)函數(shù)，計算得到該互相關(guān)函數(shù)的值最大的時差，即模型道數(shù)據(jù)與校正后的第二地震數(shù)據(jù)之間的互相關(guān)程度最大時對應(yīng)的時差，將該時差作為該時差微調(diào)對應(yīng)的互相關(guān)時差。

步驟2044：終端根據(jù)該互相關(guān)時差，對該校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)。

具體的，終端根據(jù)該互相關(guān)時差，通過以下公式九，對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，

公式九：x″′i,j(t)＝x″i,j(t-τ1)，

其中，x″′i,j(t)為時差微調(diào)后的第二地震數(shù)據(jù)，x″i,j(t-τ1)為對校正后的第二地震數(shù)據(jù)進(jìn)行的時差微調(diào)，τ1為互相關(guān)函數(shù)的最大值aj對應(yīng)的時差，即互相關(guān)時差。

需要說明的是，為了盡可能的消除校正后的第二地震數(shù)據(jù)中存在的誤差，步驟2043-2044可以多次重復(fù)執(zhí)行，直至計算出的互相關(guān)時差小于預(yù)設(shè)閾值，該預(yù)設(shè)閾值可以根據(jù)用戶需要設(shè)置并更改，本發(fā)明實施例對此不做具體限定。例如，該預(yù)設(shè)閾值可以為記錄第二地震數(shù)據(jù)的單位時間，即當(dāng)前記錄的第二地震數(shù)據(jù)最小時間間隔。

步驟2045：終端通過以下公式四，對該時差微調(diào)后的第二地震數(shù)據(jù)進(jìn)行疊加，得到該第三地震數(shù)據(jù)，

公式四：

其中，x″i(t)為該第三地震數(shù)據(jù)，dtj為該時差微調(diào)對應(yīng)的互相關(guān)時差，p為該炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最小道號，q為該炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最大道號，x”'i,j(t-dtj)該時差微調(diào)后的第二地震數(shù)據(jù)，j為地震道的道號，i為炮號，t為采集時間。

步驟205：終端根據(jù)預(yù)設(shè)切除函數(shù)和該第三地震數(shù)據(jù)，確定第四地震數(shù)據(jù)。

本步驟通過預(yù)設(shè)切除函數(shù)，對第三地震數(shù)據(jù)中存在氣泡干擾的數(shù)據(jù)進(jìn)行切除，得到?jīng)]有氣泡干擾的第四地震數(shù)據(jù)。

具體的，本步驟可以通過以下步驟2051-2052實現(xiàn)。

步驟2051：終端確定該預(yù)設(shè)切除函數(shù)。

本步驟中，該預(yù)設(shè)切除函數(shù)由通過區(qū)、斜坡區(qū)和截斷區(qū)三部分組成，具體的，該切除函數(shù)可以為：

其中，wi(t)為切除函數(shù)，i為炮號，t＝tk為地震信號的結(jié)束時間，tk表示地震信號延續(xù)的時長，即地震信號對應(yīng)的地震子波的波長為tk，t＝tw為第三地震數(shù)據(jù)的結(jié)束時間，即第三地震數(shù)據(jù)持續(xù)的總時長為tw，tb為氣泡發(fā)育時間，即采集設(shè)備開始采集到氣泡進(jìn)行第二次膨脹壓縮時產(chǎn)生的干擾波的時間。0≤t≤tk為切除函數(shù)的通過區(qū)，此時，wi(t)＝1，通過區(qū)的長度不能大于地震信號對應(yīng)的地震子波的波長tk；tk＜t＜tb為切除函數(shù)的斜坡區(qū)，此時wi(t)＝f(t-tk)，通過區(qū)和斜坡區(qū)的長度之和不能大于氣泡發(fā)育時間tb；tb≤t≤tw為切除函數(shù)的截斷區(qū)，截斷區(qū)的長度為從氣泡發(fā)育時間至第三地震數(shù)據(jù)的結(jié)束時間，

需要說明的是，通過區(qū)用于保留原始地震數(shù)據(jù)中地震信號的地震數(shù)據(jù)，截斷區(qū)用于消除原始地震數(shù)據(jù)中被氣泡干擾的地震數(shù)據(jù)，為了避免由通過區(qū)直接進(jìn)入截斷區(qū)時對原始地震數(shù)據(jù)造成的邊界抖動，此處，還定義了斜坡區(qū)。切除函數(shù)中第三地震數(shù)據(jù)的總時長tw、地震信號對應(yīng)的地震子波的波長tk、氣泡發(fā)育時間tb的關(guān)系如圖3所示。其中，斜坡區(qū)對應(yīng)的函數(shù)f(t)可以根據(jù)用戶需要設(shè)置并更改，本發(fā)明實施例對此不做具體限定。例如，該斜坡區(qū)對應(yīng)的函數(shù)f(t)可以為三角函數(shù)、漢明時窗函數(shù)等。

步驟2052：終端根據(jù)該預(yù)設(shè)切除函數(shù)和該第三地震數(shù)據(jù)，通過以下公式六，確定該第四地震數(shù)據(jù)，

公式六：y″i(t)＝x″i(t)·wi(t)，

其中，x″i(t)為該第三地震數(shù)據(jù)，wi(t)為該預(yù)設(shè)切除函數(shù)，y″i(t)為該第四地震數(shù)據(jù)，i為炮號，t為采集時間。

本步驟中，終端利用該切除函數(shù)，對第三地震數(shù)據(jù)中的氣泡干擾對應(yīng)的地震數(shù)據(jù)進(jìn)行切除，從而得到消除氣泡干擾的地震數(shù)據(jù)，終端將該消除氣泡干擾的地震數(shù)據(jù)作為第四地震數(shù)據(jù)。

步驟206：終端根據(jù)該第三地震數(shù)據(jù)和該第四地震數(shù)據(jù)，確定氣泡壓制因子。

本發(fā)明實施例中，終端將該第三地震數(shù)據(jù)作為輸入數(shù)據(jù)，將消除氣泡干擾的第四地震數(shù)據(jù)作為期望的輸出數(shù)據(jù)，通過預(yù)先定義的公式十，計算氣泡壓制因子。

本步驟可以為：終端將第三地震數(shù)據(jù)作為輸入數(shù)據(jù)，將第四地震數(shù)據(jù)作為期望的輸出數(shù)據(jù)，通過以下公式十，確定濾波結(jié)果表達(dá)式，以及期望的輸出數(shù)據(jù)和濾波結(jié)果之間的誤差表達(dá)式，計算滿足預(yù)設(shè)條件的氣泡壓制因子，預(yù)設(shè)條件為誤差表達(dá)式的值最小，

公式十：

其中，y″i(t)為第四地震數(shù)據(jù)，即期望的輸出數(shù)據(jù)，y(t)為濾波結(jié)果，e為期望的輸出數(shù)據(jù)和濾波結(jié)果之間的誤差表達(dá)式；x″i(t)為第三地震數(shù)據(jù)，即輸入數(shù)據(jù)，d(t)為氣泡壓制因子，d(t)＝(d0,d1,…dn-1)，n為第i個炮號對應(yīng)的所有地震道的總數(shù)目，i為炮號，t為采集時間。

需要說明的是，該濾波結(jié)果表達(dá)式用于表示將輸入數(shù)據(jù)經(jīng)過氣泡壓制因子壓制氣泡干擾后輸出的結(jié)果；終端通過該期望的輸出數(shù)據(jù)和該濾波結(jié)果之間的誤差表達(dá)式，計算出使得該誤差表達(dá)式的值最小的氣泡壓制因子，

本發(fā)明實施例中，終端可利用最小二乘法計算氣泡壓制因子d(t)＝(d0,d1,...dn-1)，使誤差表達(dá)式e的值最小。

具體的，首先將平方項展開，將濾波結(jié)果表達(dá)式代入公式十中的誤差表達(dá)式中，則公式十中的誤差表達(dá)式變形為：

對變形后的誤差表達(dá)式取偏微分，并且e取值為0，可將該變形后的誤差表達(dá)式表示為：

其中，(r0,r1,…rn-1)為輸入數(shù)據(jù)x″i(t)的自相關(guān)函數(shù)，(g0,g1,…gn-1)為第四地震數(shù)據(jù)，即期望的輸出數(shù)據(jù)y″i(t)與輸入數(shù)據(jù)x″i(t)之間的互相關(guān)函數(shù)。利用levinson(萊文森)遞歸算法，計算出氣泡壓制因子。

步驟207：終端根據(jù)該氣泡壓制因子，對每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，使得每個地震道的原始地震數(shù)據(jù)中氣泡干擾被壓制。

具體的，本步驟可以為：終端根據(jù)該氣泡壓制因子，通過以下公式七，對每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，得到壓制氣泡干擾后的每個地震道的地震數(shù)據(jù)，

公式七：yi,j(t)＝xi,j(t)*d(t)，

其中，xi,j(t)為該地震道集合中第j個地震道第i個炮號的原始地震數(shù)據(jù)，d(t)為該氣泡壓制因子，yi,j(t)為該壓制氣泡干擾后的每個地震道的地震數(shù)據(jù)，t為采集時間。

本步驟中，終端根據(jù)該氣泡壓制因子對該炮號對應(yīng)的地震道的原始地震數(shù)據(jù)進(jìn)行濾波，即，該每個地震道的原始數(shù)據(jù)為該該炮號對應(yīng)的地震道的原始地震數(shù)據(jù)。

需要說明的是，終端利用該氣泡壓制因子，對第i個炮號對應(yīng)的所有地震道的原始地震數(shù)據(jù)進(jìn)行濾波，得到第i個炮號對應(yīng)的地震道集合的原始地震數(shù)據(jù)中的氣泡干擾被壓制后，終端可以繼續(xù)對其它炮號中每炮對應(yīng)的地震道集合的原始地震數(shù)據(jù)均執(zhí)行步驟201-207，消除每炮對應(yīng)的地震道集合的原始地震數(shù)據(jù)中被氣泡干擾的原始地震數(shù)據(jù)，使得每炮對應(yīng)的地震道集合的原始地震數(shù)據(jù)中氣泡干擾被壓制。

本發(fā)明實施例中，終端通過獲取炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的第一地震數(shù)據(jù)，從第一地震數(shù)據(jù)中確定下行初至波對應(yīng)的第二地震數(shù)據(jù)，該第二地震數(shù)據(jù)為信噪比較高的地震數(shù)據(jù)；終端校正該第二地震數(shù)據(jù)，并對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，得到消除時差和誤差的第三地震數(shù)據(jù)。從而終端可以根據(jù)預(yù)設(shè)切除函數(shù)和該第三地震數(shù)據(jù)，確定消除氣泡干擾的第四地震數(shù)據(jù)；然后，終端根據(jù)該第三地震數(shù)據(jù)和該第四地震數(shù)據(jù)，確定氣泡壓制因子，終端根據(jù)該氣泡壓制因子，對每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，從而使得每個地震道的原始地震數(shù)據(jù)中氣泡干擾被壓制，終端通過對每個下行初至波對應(yīng)的第二地震數(shù)據(jù)進(jìn)行校正、時差微調(diào)，根據(jù)得到的第三地震數(shù)據(jù)以及消除氣泡干擾的第四地震數(shù)據(jù)，計算出氣泡壓制因子，即，計算該氣泡壓制因子時綜合考慮到了每個下行初至波的第二地震數(shù)據(jù)，通過該氣泡壓制因子對原始地震數(shù)據(jù)進(jìn)行濾波，從而提高了對原始地震數(shù)據(jù)進(jìn)行壓制氣泡干擾的準(zhǔn)確性。

本發(fā)明實施例提供了一種基于地震數(shù)據(jù)的壓制氣泡干擾的裝置，該裝置可以應(yīng)用在終端中，如圖4所示，該裝置包括：

獲取模塊401，用于獲取地震道集合中每個地震道的原始地震數(shù)據(jù)；

選擇模塊402，用于從該每個地震道的原始地震數(shù)據(jù)中選擇炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的第一地震數(shù)據(jù)；

第一確定模塊403，用于從選擇的第一地震數(shù)據(jù)中確定下行初至波對應(yīng)的第二地震數(shù)據(jù)；

校正模塊404，用于校正該第二地震數(shù)據(jù)，對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，得到第三地震數(shù)據(jù)；

第二確定模塊405，用于根據(jù)預(yù)設(shè)切除函數(shù)和該第三地震數(shù)據(jù)，確定第四地震數(shù)據(jù)；

第三確定模塊406，用于根據(jù)該第三地震數(shù)據(jù)和該第四地震數(shù)據(jù)，確定氣泡壓制因子；

濾波模塊407，用于根據(jù)該氣泡壓制因子，對每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，使得該每個地震道的原始地震數(shù)據(jù)中氣泡干擾被壓制。

在一種可能的設(shè)計中，該校正模塊404，包括：

校正單元，用于確定該每個地震道的下行初至波的旅行時間，通過以下公式一，對該第二地震數(shù)據(jù)進(jìn)行校正，

公式一：x″i,j(t)＝x'i,j(t+ftj)，

其中，ftj為第j個地震道的下行初至波的旅行時間，x'i,j(t)為該第二地震數(shù)據(jù)，x″i,j(t)為校正后的第二地震數(shù)據(jù)，t為采集時間，j為地震道的道號，i為炮號；

第一疊加單元，用于通過以下公式二，對該校正后的第二地震數(shù)據(jù)進(jìn)行疊加，得到模型道數(shù)據(jù)，

公式二：

其中，x'i(t)為該模型道數(shù)據(jù)，x″i,j(t)為該校正后的第二地震數(shù)據(jù)，t為采集時間，j為地震道的道號，i為炮號，p為炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最小道號，q為該炮檢距在該預(yù)設(shè)范圍內(nèi)的地震道的最大道號；

第一確定單元，用于根據(jù)該模型道數(shù)據(jù)與該校正后的第二地震數(shù)據(jù)，通過以下公式三，確定該時差微調(diào)對應(yīng)的互相關(guān)時差，

公式三：

其中，aj(τ)為該模型道數(shù)據(jù)與該校正后的第二地震數(shù)據(jù)之間的互相關(guān)函數(shù)，aj為該互相關(guān)函數(shù)的最大值，dtj為該時差微調(diào)對應(yīng)的互相關(guān)時差，即該互相關(guān)函數(shù)的最大值aj對應(yīng)的時差τ1，x'i(t)為該模型道數(shù)據(jù)，x″i,j(t)為該校正后的第二地震數(shù)據(jù)，j為地震道的道號，i為炮號，t為采集時間，m為該互相關(guān)函數(shù)的半支長度；

微調(diào)單元，用于根據(jù)該互相關(guān)時差，對該校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)；

第二疊加單元，用于通過以下公式四，對該時差微調(diào)后的第二地震數(shù)據(jù)進(jìn)行疊加，得到該第三地震數(shù)據(jù)，

公式四：

其中，x″i(t)為該第三地震數(shù)據(jù)，dtj為該時差微調(diào)對應(yīng)的互相關(guān)時差，p為該炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最小道號，q為該炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的最大道號，x”'i,j(t-dtj)該時差微調(diào)后的第二地震數(shù)據(jù)，j為地震道的道號，i為炮號，t為采集時間。

在一種可能的設(shè)計中，該第二確定模塊405，包括：

第二確定單元，用于確定該預(yù)設(shè)切除函數(shù)；

第三確定單元，用于根據(jù)該預(yù)設(shè)切除函數(shù)和該第三地震數(shù)據(jù)，通過以下公式六，確定該第四地震數(shù)據(jù)，

公式六：y″i(t)＝x″i(t)·wi(t)，

其中，x″i(t)為該第三地震數(shù)據(jù)，wi(t)為該預(yù)設(shè)切除函數(shù)，y″i(t)為該第四地震數(shù)據(jù)，i為炮號，t為采集時間。

在一種可能的設(shè)計中，該第三確定模塊406，還用于將該第三地震數(shù)據(jù)作為輸入數(shù)據(jù)，將該第四地震數(shù)據(jù)作為期望的輸出數(shù)據(jù)，通過以下公式十，確定濾波結(jié)果表達(dá)式，以及該期望的輸出數(shù)據(jù)和該濾波結(jié)果之間的誤差表達(dá)式，計算滿足預(yù)設(shè)條件的氣泡壓制因子，該預(yù)設(shè)條件為該誤差表達(dá)式的值最小，

公式十：

其中，y″i(t)為該第四地震數(shù)據(jù)，即該期望的輸出數(shù)據(jù)，y(t)為該濾波結(jié)果，e為該期望的輸出數(shù)據(jù)和該濾波結(jié)果之間的誤差表達(dá)式；x″i(t)為該第三地震數(shù)據(jù)，即該輸入數(shù)據(jù)，d(t)為該氣泡壓制因子，i為炮號，t為采集時間。

在一種可能的設(shè)計中，該濾波模塊407，還用于根據(jù)該氣泡壓制因子，通過以下公式七，對該每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，得到壓制氣泡干擾后的每個地震道的地震數(shù)據(jù)，

公式七：yi,j(t)＝xi,j(t)*d(t)，

其中，xi,j(t)為該地震道集合中第j個地震道第i個炮號的原始地震數(shù)據(jù)，d(t)為該氣泡壓制因子，yi,j(t)為該壓制氣泡干擾后的每個地震道的地震數(shù)據(jù)，t為采集時間。

本發(fā)明實施例中，終端通過獲取炮檢距在預(yù)設(shè)范圍內(nèi)的地震道的第一地震數(shù)據(jù)，從第一地震數(shù)據(jù)中確定下行初至波對應(yīng)的第二地震數(shù)據(jù)，該第二地震數(shù)據(jù)為信噪比較高的地震數(shù)據(jù)；終端校正該第二地震數(shù)據(jù)，并對校正后的第二地震數(shù)據(jù)進(jìn)行時差微調(diào)，得到消除時差和誤差的第三地震數(shù)據(jù)。從而終端可以根據(jù)預(yù)設(shè)切除函數(shù)和該第三地震數(shù)據(jù)，確定消除氣泡干擾的第四地震數(shù)據(jù)；然后，終端根據(jù)該第三地震數(shù)據(jù)和該第四地震數(shù)據(jù)，確定氣泡壓制因子，終端根據(jù)該氣泡壓制因子，對每個地震道的原始地震數(shù)據(jù)進(jìn)行濾波，從而使得每個地震道的原始地震數(shù)據(jù)中氣泡干擾被壓制，終端通過對每個下行初至波對應(yīng)的第二地震數(shù)據(jù)進(jìn)行校正、時差微調(diào)，根據(jù)得到的第三地震數(shù)據(jù)以及消除氣泡干擾的第四地震數(shù)據(jù)，計算出氣泡壓制因子，即，計算該氣泡壓制因子時綜合考慮到了每個下行初至波的第二地震數(shù)據(jù)，通過該氣泡壓制因子對原始地震數(shù)據(jù)進(jìn)行濾波，從而提高了對原始地震數(shù)據(jù)進(jìn)行壓制氣泡干擾的準(zhǔn)確性。

需要說明的是：上述實施例提供的基于地震數(shù)據(jù)的壓制氣泡干擾的裝置在基于地震數(shù)據(jù)的壓制氣泡干擾時，僅以上述各功能模塊的劃分進(jìn)行舉例說明，實際應(yīng)用中，可以根據(jù)需要而將上述功能分配由不同的功能模塊完成，即將裝置的內(nèi)部結(jié)構(gòu)劃分成不同的功能模塊，以完成以上描述的全部或者部分功能。另外，上述實施例提供的基于地震數(shù)據(jù)的壓制氣泡干擾的裝置與基于地震數(shù)據(jù)的壓制氣泡干擾的方法實施例屬于同一構(gòu)思，其具體實現(xiàn)過程詳見方法實施例，這里不再贅述。

本領(lǐng)域普通技術(shù)人員可以理解實現(xiàn)上述實施例的全部或部分步驟可以通過硬件來完成，也可以通過程序來指令相關(guān)的硬件完成，所述的程序可以存儲于一種計算機(jī)可讀存儲介質(zhì)中，上述提到的存儲介質(zhì)可以是只讀存儲器，磁盤或光盤等。

以上所述僅為本發(fā)明的較佳實施例，并不用以限制本發(fā)明，凡在本發(fā)明的精神和原則之內(nèi)，所作的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。