本發明屬于海洋環境以及內陸湖泊的氣槍震源地震勘探采集、地質調查等
技術領域：
，尤其涉及一種實現氣槍震源質控的處理方法及裝置。
背景技術：
：氣槍震源地震勘探采集技術在海洋環境得到了廣泛的應用，隨著勘探技術水平的不斷發展，對氣槍震源的質量控制要求越來越高，體現在例如氣槍壓力、容量、沉放深度、不同步、自激和氣泡壓制等方面。氣槍震源類型主要為空氣槍組合模式，總槍數一般在20-40支之間。組合槍陣的槍總數越多，槍的不同步、自激和氣泡壓制等不良現象越容易出現，而氣泡壓制效果的好壞，對單炮地震采集的數據結果影響重大。槍的不同步故障可以通過槍控系統及其log文件可以快速識別。當氣槍陣列氣泡壓制效果差時，氣槍陣列就會產生嚴重的氣泡效應，對應的地震單炮就會出現異常初至，類似重復沖擊現象，直接影響獲取的地震資料品質。一般來說，氣槍震源陣列設計中采用相干槍組合以及優化陣列的沉放深度可以將震源激發產生的氣泡振幅變小，達到壓制氣泡的目的。但是，目前現有公開資料還沒有關于氣槍陣列氣泡壓制方面的質控判別方法，以實現定性定量來判別氣槍陣列是否存在故障。技術實現要素：本申請目的在于提供一種實現氣槍震源質控的處理方法及裝置，可以利用近場子波信息快速有效地識別氣槍陣列是否存在氣泡壓制效果差的故障，從而有針對性去排除相關氣槍故障，進而保證地震采集資料品質，實現氣槍震源質控監控和處理。本申請提供的一種實現氣槍震源質控的處理方法及裝置是這樣實現的：一種實現氣槍震源質控的處理方法，所述方法包括：將單炮地震記錄對應的近場子波記錄在距離氣槍震源陣列中心最近的檢波器接收點處進行線性疊加，形成疊加后近場子波記錄；計算所述疊加后近場子波記錄的氣泡比值，所述氣泡比值包括所述疊加后近場子波記錄的第一個壓力脈沖的最大振幅值與疊加后近場子波記錄的第一個氣泡脈沖的最大振幅值之比；取相同水深范圍內的連續k個正常單炮對應的近場子波記錄的氣泡比值作的平均值，作為正常單炮的標準氣泡比值，k≥3；基于所述標準氣泡比值來確定所述氣槍震源地震采集數據的品質信息。優選的實施例中，利用下述公式實現所述線性疊加：上式中，f(t)表示為疊加子波，f(t)為單槍子波，g(t)為虛反射子波，di為第i條氣槍距所述最近的檢波器接收點的距離，di'為第i條氣槍的虛像距所述最近的檢波器接收點的距離，n為陣列中氣槍個數，c為聲波在水中的速度。優選的實施例中，所述k的取值設置為：k≥10。優選的實施例中，所述k的取值設置為10。優選的實施例中，所述基于所述標準氣泡比值來確定所述氣槍震源地震采集數據的品質信息，包括：當有出現連續m個以上單炮對應的氣泡比值小于所述標準氣泡比值時，則確定不符合氣槍陣列氣泡壓制要求，m≥1。優選的實施例中，所述m的取值為5。優選的實施例中，當有出現連續m個以上單炮對應的氣泡比值小于所述標準氣壓比值的預定百分比時，則確定不符合氣槍陣列氣泡壓制要求，m≥1。一種實現氣槍震源質控的處理裝置，所述裝置包括處理器以及存儲所述處理器可執行指令的存儲器，所述處理器執行指令時可以實現下：將單炮地震記錄對應的近場子波記錄在距離氣槍震源陣列中心最近的檢波器接收點處進行線性疊加，形成疊加后近場子波記錄；計算所述疊加后近場子波記錄的氣泡比值，所述氣泡比值包括所述疊加后近場子波記錄的第一個壓力脈沖的最大振幅值與疊加后近場子波記錄的第一個氣泡脈沖的最大振幅值之比；取相同水深范圍內的連續k個正常單炮對應的近場子波記錄的氣泡比值作的平均值，作為正常單炮的標準氣泡比值，k≥3；基于所述標準氣泡比值來確定所述氣槍震源地震采集數據的品質信息。優選的實施例中，所述存儲器存儲的指令被設置成包括使得處理器執行時實現利用下述公式實現所述線性疊加：上式中，f(t)表示為疊加子波，f(t)為單槍子波，g(t)為虛反射子波，di為第i條氣槍距所述最近的檢波器接收點的距離，di'為第i條氣槍的虛像距所述最近的檢波器接收點的距離，n為陣列中氣槍個數，c為聲波在水中的速度。優選的實施例中，所述k的取值設置為：k≥10。優選的實施例中，所述基于所述標準氣泡比值來確定所述氣槍震源地震采集數據的品質信息，包括：當有出現連續m個以上單炮對應的氣泡比值小于所述標準氣泡比值時，則確定不符合氣槍陣列氣泡壓制要求，m≥1。或者，當有出現連續m個以上單炮對應的氣泡比值小于所述標準氣壓比值的預定百分比時，則確定不符合氣槍陣列氣泡壓制要求，m≥1。優選的實施例中，所述m的取值為5。優選的實施例中，所述預定百分比為50％。本申請提供的一種實現氣槍震源質控的處理方法及裝置，可以通過計算近場子波的氣泡比并結合地震單炮記錄，快速有效識別氣槍陣列是否存在氣泡壓制效果差而導致地震單炮記錄出現異常初至的現象，保證地震采集資料品質，填補了一項氣槍震源地震采集質控技術的空白。附圖說明為了更清楚地說明本申請實施例或現有技術中的技術方案，下面將對實施例或現有技術描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本申請中記載的一些實施例，對于本領域普通技術人員來講，在不付出創造性勞動性的前提下，還可以根據這些附圖獲得其他的附圖。圖1是本申請提供的一種實現氣槍震源質控的處理方法一種實施例的方法流程圖；圖2是本申請提供的一種實現氣槍震源質控的處理裝置一種實施例的模塊結構示意圖；圖3是本申請一個示例中連續10炮正常單炮記錄示意圖；圖4是本申請一個示例中10炮正常單炮記錄對應的線性疊加近場子波記錄示意圖；圖5是本申請一個示例中連續5炮異常單炮記錄的示意圖；圖6是本申請一個示例中多個疊加后近場子波記錄的示意圖，圖7是本申請一個示例中項目異常單炮及其對應的近場子波示意圖圖8是本申請一個示例中正常單炮及其對應的近場子波的示意圖；圖9是本申請一個示例中異常單炮及其對應的近場子波的示意圖。具體實施方式為了使本
技術領域：
的人員更好地理解本申請中的技術方案，下面將結合本申請實施例中的附圖，對本申請實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是本申請一部分實施例，而不是全部的實施例。基于本申請中的實施例，本領域普通技術人員在沒有作出創造性勞動前提下所獲得的所有其他實施例，都應當屬于本申請保護的范圍。氣槍震源的槍陣上通常會安裝了近場子波檢波器得到近場子波信息，一般來說，近場子波主要有兩大用途，一是模擬遠場子波，二是監控震源工作狀態。目前實際應用中，通常都是利用近場子波監控震源主要是檢查單槍是否掛錯，判斷定深繩的深度是否有錯，觀測其所在位置的槍是否不點火、是否存在延遲錯和自激現象(定性)等，而沒有利用近場子波所獲得的信息進行氣泡壓制方面的質控，沒有形成一種快速有效的質控方法來定性定量來判別氣槍陣列是否存在故障。現有技術中也有人提出些氣泡比的計算方法，主要用于氣槍自激等故障現象，并且沒有同地震記錄結合。因此，本申請提供的一種實現氣槍震源質控的處理方法對實際地震資料采集工作具有非常重要的意義。本申請提供的方法主要利用在定量分析槍陣近場子波信息的基礎上，結合近場子波上的異常信息與單炮的地震波對應關系，快速有效地判別氣槍陣列是否存在氣泡壓制效果差而導致單炮出現異常初至的故障。圖1是本申請所述一種實現氣槍震源質控的處理方法一種實施例的方法流程圖。雖然本申請提供了如下述實施例或附圖所示的方法操作步驟或裝置結構，但基于常規或者無需創造性的勞動在所述方法或裝置中可以包括更多或者更少的操作步驟或模塊單元。在邏輯性上不存在必要因果關系的步驟或結構中，這些步驟的執行順序或裝置的模塊結構不限于本申請實施例或附圖所示的執行順序或模塊結構。所述的方法或模塊結構的在實際中的裝置或終端產品應用時，可以按照實施例或者附圖所示的方法或模塊結構進行順序執行或者并行執行(例如并行處理器或者多線程處理的環境、甚至包括分布式處理的實施環境)。本申請一個實施例提供的方法，可以通過計算近場子波的氣泡比并結合地震單炮記錄，快速有效識別氣槍陣列是否存在氣泡壓制效果差而導致地震單炮記錄出現異常初至的現象。氣槍陣列可以由單個或者多個子陣列組成，而子陣列是由多支單槍或者多組相干槍組成。一般情況下，在單槍或者相干槍的震源中心上方1m處掛接有近場子波傳感器。氣槍震源每激發一次，這些近場子波傳感器能夠接收到氣槍震源產生的近場子波信息，并通過bigshot槍控系統記錄到hyd文件中(類似seg-y文件)，同時地震儀器能接收到相應的地震采集信息并形成單炮記錄。在本實施例中，可以獲取一些計算所需的已知參數，例如氣槍陣列中每個單槍的坐標位置、近場子波傳感器的坐標位置、離陣列震源中心最近距離的檢波點坐標位置、海面反射系數、聲波在水中的速度、地震單炮記錄及氣槍陣列實際接收的近場子波等。這些所述已知參數對于本領域技術人員來說可以在實際現場應用中讀取、測繪或計算得到，當然，不排除還可能使用到其他參數，或者與本發明已知參數等同或變形、變換后的其他參數，在此不做詳細的一一論述。具體的如圖1所述，本申請提供的一種實現氣槍震源質控的處理方法的一種實施例可以包括：s1：將單炮地震記錄對應的近場子波記錄在距離氣槍震源陣列中心最近的檢波器接收點處進行線性疊加，形成疊加后近場子波記錄。在本實施例實時場景中，氣槍震源激發時，可以實時獲取每一炮的單炮地震記錄所對應的所有近場子波記錄，具體的如前述通過近場子波傳感器接收到的氣槍震源產生的近場子波信息。然后可以將每一炮的地震單炮記錄對應的近場子波記錄進行線性疊加，疊加的方式可以采集在距離氣槍震源陣列中心最近的檢波器接收點的位置處。所述的線性疊加，主要是指指根據激發點和接收點的位置關系所計算的時移量，在接收點位置將所有近場子波道進行時移求和，從而得到一道新的疊加后的近場子波，該近場子波同時考慮了氣槍激發和接收兩個環節對地震資料的影響，其近場子波特征值(氣泡比，等)為資料評價提供了依據。經過線性疊加后，可以形成一道新的疊加后的近場子波記錄。本申請提供的一種實施例中，可以采用下述方式實現近場子波記錄的線性疊加。具體的，本申請所述方法的另一種實施例中，可以利用下述公式實現所述線性疊加：上式中，f(t)表示為疊加子波(即經過疊加后形成的疊加后近場子波記錄)，f(t)為單槍子波，g(t)為虛反射子波，di為第i條氣槍距所述最近的檢波器接收點的到距離，di'為第i條氣槍的虛像距所述最近的檢波器接收點的距離，n為陣列中氣槍個數，c為聲波在水中的速度。虛反射子波數據的獲取，可以利用虛反射子波與直達波在海平面位置極性相反，振幅相同，其能量與傳播距離成反比這些信息確定得到。基于虛像和震源關于海平面對稱這一信息可以獲取虛像距離接收點的距離。當前，上述提供的公式僅是實現本實施例線性疊加的方式之一，通過參數的變形、公式的變換等還可以存在其他的疊加公式。應當理解的是，基于本實施例單槍子波和虛反射子波在所述最近的檢波器的求和相減后得到的結果均應屬于本申請的技術實施范圍，上述公式僅是給出的實際可操作的一種實施例的計算實現方式。本申請實施例中，可以采集獲取氣槍震源產生的近場子波記錄，然后將單炮地震記錄對應的近場子波記錄在距離氣槍震源陣列中心最近的檢波器接收點處進行線性疊加，形成疊加后近場子波記錄。s2：計算所述疊加后近場子波記錄的氣泡比值，所述氣泡比值包括所述疊加后近場子波記錄的第一個壓力脈沖的最大振幅值與疊加后近場子波記錄的第一個氣泡脈沖的最大振幅值之比。其中，所述的第一個壓力脈沖與第一個氣泡脈沖所在位置可以通過氣槍陣列的設計氣泡周期獲取，通過選定時間窗口，可以取最大真值為最大振幅值。s3：取相同水深范圍內的連續k個正常單炮對應的近場子波記錄的氣泡比值作的平均值，作為正常單炮的標準氣泡比值，k≥3。具體的處理過程可以包括在相同水深范圍內(可以允許存在一定誤差，例如20米±2米，此時應當理解符號本實施例所述的相同水深范圍)，取連續k個正常單炮對應的近場子波記錄進行疊加后計算的氣泡比值，并計算著k個氣泡比值的平均值，以這個平均值作為正常單炮的標準氣泡比值，這里的正常單炮可以是指槍陣的壓力、容量、沉放深度相同，且單炮上沒有不同步、自激、異常初至現象。當前，也可以根據實際現場作業情況或質控要求設置其他的正常單炮條件。在本實施例中，所述的k的取值范圍可以為至少3個正常單炮的氣泡比值(一個無意義，兩個幾乎無實用價值)。在本申請實施方案中，一般情況下，如果計算標準氣泡比值所使用的正常單炮的取值個數較少，則可靠性較差，對于評判氣槍震源的質控結果參考意義和準確性較低。基于此，本申請的一種實施方式中，所述的k的取值為大于等于3。本申請的另一個實施例中，結合實際多數場景下的質控預期需求，可以取相同水深范圍內的連續10炮或者以上正常單炮對應的氣泡比值作平均，所述的k的取值可以設置為：k≥10。當然，優選的另一個實施例中，k的取值可以為10，這樣，技能滿足一定的質控要求，由于是10的倍數，數據處理和質控的意義表征等都易于作業人員理解和計算。需要說明的是，這個水深范圍可以有作業人員根據實際現場施工情況和質控要求進行選取和設置。另外，當出現多處連續k個符合要求的氣泡比值時，一種實施方式可以任意選取其中一段作為的平均值作為標準氣泡比值。例如，當有編號1010～1020的連續k個正常單炮的氣泡比值和1050～1060的連續k個正常單炮的氣泡比值時，可以任意選擇其中一段的平均值作為質控的標準氣泡比值。s4：基于所述標準氣泡比值來確定所述氣槍震源地震采集數據的品質信息。通過上述方式確定了用于評判地震資料品質的標準氣泡比值，則可以進一步利用該標準氣泡比值來實現氣槍震源質控。具體的利用所述標準氣泡比值實現質控要求的處理方式可以根據實際的質控要求采用不同的處理方式。例如直接將采集的單炮進行計算后與標準氣泡比值進行比較，如果滿足一定的預置要求，則可以認為該單炮符合質控要求。也可以在聯系出現多個單炮的氣泡比值過小時判定出現異常初至現場等。本申請實施例提供的一種實現氣槍震源質控的處理方法，能夠及時發現氣槍震源激發作業過程中因氣泡壓制效果差而導致單炮出現異常初至的現象，保證地震采集資料品質，填補了一項氣槍震源地震采集質控技術的空白。所述方法的另一種實施例提供了一個具體的量化輸出質控處理方法。具體的，所述方法的另一種實施例中，所述基于所述標準氣泡比值來確定所述氣槍震源地震采集數據的品質信息可以包括：s401：當有出現連續m個以上單炮對應的氣泡比值小于所述標準氣泡比值時，則確定不符合氣槍陣列氣泡壓制要求，m≥1。此時可以輸出告警或提示信息。一種實施方式中，可以認為當出現地震單炮對應的氣泡比小于標準氣泡比值時則氣泡壓制效果差。其他的實施方式中，一般的，可以在聯系出現多個氣泡比值小于標準氣泡比值時才認為出現因氣泡壓制效果差而出現單炮異常初至的現場。一般在海上地震勘探中，不允許有連續5炮單炮出現質量問題，如果出現，則要求將這5炮進行重新補炮。因此，本申請結合實際應用效果提供所述方法的另一種實施例中，所述m的取值為5。當有超過連續5炮地震單炮記錄對應的氣泡比值小于標準氣泡比值時，可以將地震單炮與其對應的疊加后的近場子波進行記錄時間上的對應，若近場子波氣泡脈沖與單炮上出現的異常初至能夠對應上，說明整個槍陣壓制氣泡效果差，需要進行槍陣檢修。或者，另一種實施例中，所述基于所述標準氣泡比值來確定所述氣槍震源地震采集數據的品質信息可以包括：s402：當有出現連續m個以上單炮對應的氣泡比值小于所述標準氣壓比值的預定百分比時，則確定不符合氣槍陣列氣泡壓制要求。此時可以輸出告警或提示信息。這個實施例的實施場景下，可以直接與標準氣泡比值進行比較，而是與其一定的百分比的數值比較。本申請具體的一個實施例中，所述預定百分比為50％，即可以將連續出現m個以上的單炮的氣泡比值均小于標準氣泡比值的50％時，可以將地震單炮與其對應的疊加后的近場子波進行記錄時間上的對應，若近場子波氣泡脈沖與單炮上出現的異常初至能夠對應上，說明整個槍陣壓制氣泡效果差，需要進行槍陣檢修。本申請上述實施例提供的時序氣槍震源質控的處理方法，明能夠及時發現氣槍震源激發作業過程中因氣泡壓制效果差而導致單炮出現異常初至的現象，保證地震采集資料品質，填補了一項氣槍震源地震采集質控技術的空白。基于上述所述的實現氣槍震源質控的處理方法，本申請還提供一種實現氣槍震源質控的處理裝置。所述的裝置可以包括使用了本申請所述方法的系統(包括分布式系統)、軟件(應用)、模塊、組件、服務器等并結合必要的實施硬件的裝置。基于同一創新構思，本申請提供的一種實施例中的裝置如下面的實施例所述。由于裝置解決問題的實現方案與方法相似，因此本申請具體的裝置的實施可以參見前述方法的實施，重復之處不再贅述。以下所使用的，術語“單元”或者“模塊”可以實現預定功能的軟件和/或硬件的組合。盡管以下實施例所描述的裝置較佳地以軟件來實現，但是硬件，或者軟件和硬件的組合的實現也是可能并被構想的具體的。圖2是本申請提供的一種實現氣槍震源質控的處理裝置一種實施例的模塊結構示意圖，如圖2所示，所述裝置可以包括處理器以及存儲所述處理器可執行指令的存儲器，所述處理器執行指令時可以實現下：將單炮地震記錄對應的近場子波記錄在距離氣槍震源陣列中心最近的檢波器接收點處進行線性疊加，形成疊加后近場子波記錄；計算所述疊加后近場子波記錄的氣泡比值，所述氣泡比值包括所述疊加后近場子波記錄的第一個壓力脈沖的最大振幅值與疊加后近場子波記錄的第一個氣泡脈沖的最大振幅值之比；取相同水深范圍內的連續k個正常單炮對應的近場子波記錄的氣泡比值作的平均值，作為正常單炮的標準氣泡比值，k≥3；基于所述標準氣泡比值來確定所述氣槍震源地震采集數據的品質信息。一種實施例中，所述存儲器存儲的指令被設置成包括使得處理器執行時實現利用下述公式實現所述線性疊加：上式中，f(t)表示為疊加子波(即經過疊加后形成的疊加后近場子波記錄)，f(t)為單槍子波，g(t)為虛反射子波，di為第i條氣槍距所述最近的檢波器接收點的距離，di'為第i條氣槍的虛像距所述最近的檢波器接收點的距離，n為陣列中氣槍個數，c為聲波在水中的速度。所述裝置優選的一種實施例中，所述k的取值設置為：k≥10。如前所述，所述裝置的另一種實施例中，所述基于所述標準氣泡比值來確定所述氣槍震源地震采集數據的品質信息，包括：當有出現連續m個以上單炮對應的氣泡比值小于所述標準氣泡比值時，則確定不符合氣槍陣列氣泡壓制要求，m≥1。或者，當有出現連續m個以上單炮對應的氣泡比值小于所述標準氣壓比值的預定百分比時，則確定不符合氣槍陣列氣泡壓制要求，m≥1。上述兩種方式中的m，取值均可以如前述方法所述的為5，所述預定百分比可以取值為50％。即當有超過連續5炮地震單炮記錄對應的氣泡比值小于或者等于標準氣泡比值的50％時,將地震單炮與其對應的疊加后的近場子波進行記錄時間上的對應，若近場子波氣泡脈沖與單炮上出現的異常初至能夠對應上，說明整個槍陣壓制氣泡效果差，需要進行槍陣檢修。上述實施例所述的裝置具體的實現方式可以參考前述方法實施例的相關描述，在此不做一一贅述。本申請實施例提供的方法或裝置，利用近場子波信息快速有效地識別氣槍陣列是否存在氣泡壓制效果差的故障，從而有針對性去排除相關氣槍故障，進而保證地震采集資料品質，實現氣槍震源質控監控和處理。下述是利用本申請方法或裝置的一個具體應用，以此來進一步說明本申請的創新思想及實際操作應用。本申請實施例在某底電纜三維地震資料采集項目中得到了充分應用。工區位于渤海灣南部，該項目使用的震源類型為空氣槍陣列組合，組合方式為3排線性面積組合，共有29支槍，掛接有17個近場子波傳感器。具體實施包括如下：1、取連續10炮正常單炮記錄(如圖3所示，圖3是，圖4是圖3中10炮正常單炮對應的線性疊加近場子波示意圖)對應的近場子波記錄分別進行線性疊加。陣列中單槍、近場子波傳感器的位置坐標如表1所示。每一單炮最近的接收檢波點位置坐標如表2所示，海水聲速為1500m/s，海面反射系數為-1。表1陣列中單槍、近場子波傳感器的位置坐標表210炮正常單炮記錄最近的檢波點接收位置坐標最近接收點x坐標y坐標z坐標單炮112.634.722.0單炮215.823.922.0單炮321.027.822.0單炮419.225.521.9單炮515.226.521.8單炮618.027.721.9單炮716.028.222.0單炮815.826.721.9單炮919.823.221.9單炮1019.826.021.8通過線性疊加公式，將每一單炮對應的17道近場子波疊加為1道近場子波，如圖4所示，圖4中左側表示氣槍震源每激發一次，所有近場子波傳感器接收到的近場子波信息。圖右側表示所有的近場子波進行線性疊加后得到一道新的近場子波記錄。圖5是多個疊加后近場子波記錄的示意圖。2、計算線性疊加后的近場子波的氣泡比(p/bration)值，見表3。表310炮正常單炮對應的疊加近場子波的氣泡比值將10炮正常單炮對應的疊加近場子波的氣泡比值作平均，得到正常單炮的標準氣泡比值為18.46。3、在施工163線束時，通過計算，發現連續5炮單炮記錄對應的疊加近場子波氣泡比在4.20-5.49之間，低于標準氣泡比值的50％。每一單炮最近的接收檢波點位置坐標如表4所示，疊加近場子波氣泡比計算結果如表5所示。通過單炮與近場子波對應，發現近場子波的氣泡脈沖與單炮上出現的異常初至能夠對應上(如圖6、圖7所示，圖6是連續5炮異常單炮記錄示意圖，圖7是項目異常單炮及其對應的近場子波示意圖)，得出氣槍陣列存在故障，通過檢修后，槍陣工作狀態正常。表45炮異常單炮記錄最近的檢波點接收位置坐標表55炮異常單炮對應的疊加近場子波的氣泡比值單炮記錄單炮1單炮2單炮3單炮4單炮5平均值氣泡比5.495.264.384.204.354.74圖8是本申請一個示例中正常單炮及其對應的近場子波的示意圖，圖8中左側為正常地震單炮，右側為正常單炮對應的近場子波，其氣泡比值為19.97。圖9是本申請一個示例中異常單炮及其對應的近場子波的示意圖，圖9中，最左側為異常地震單炮對應的所有近場子波，圖中間為線性疊加后的近場子波，其氣泡比為4.73，最右側為對應的異常單炮。在圖9中，因近場子波是在震源激發前50ms開始接收信號，而地震單炮是在震源激發同時接收的，所以近場子波50ms時刻對應地震單炮的0ms時刻，從圖9中可以看出近場子波的氣泡脈沖與地震單炮異常初至對應(圖中虛線標志線)。該項目在將近4個月的作業時間里共計生產炮數為132772炮，通過本申請提供的質控方法，資料合格率達到了100％，單炮中沒有因氣泡壓制效果差而導致的異常初至現象，實際效果非常明顯，用戶認可度非常高。本發明不僅保證了cfd項目最佳地震資料的取得，同時也為項目高效運行保駕護航，也為海洋地震資料采集提供了很好的借鑒參考經驗，值得大力推廣和應用，將成為氣槍震源地震采集質量控制不可或缺的技術手段。本申請的一個實施例中，采取了包括：a、取相同水深范圍內的連續10炮或者以上正常單炮對應的近場子波的氣泡比值作平均，得到正常單炮的標準氣泡比值。這里的正常單炮是指槍陣的壓力、容量、沉放深度相同，且單炮上沒有不同步、自激、異常初至現象；b、結合氣槍震源地震采集空廢炮的質控要求，當有超過連續5炮地震單炮記錄對應的氣泡比值小于或者等于標準氣泡比值的50％時，將地震單炮與其對應的疊加后的近場子波進行記錄時間上的對應，若近場子波氣泡脈沖與單炮上出現的異常初至能夠對應上，說明整個槍陣壓制氣泡效果差，需要進行槍陣檢修。這兩步的理解是將地震記錄與近場子波的氣泡比值結合起來，得到氣槍陣列是否存在壓制氣泡效果差的現象。現有技術中的一些方法中并沒有提到地震記錄。也就是說地震記錄存在異常初至現象，并沒有人提出這個是槍陣壓制氣泡差引起的。槍陣壓制氣泡差而導致地震記錄存在異常初至現象，是個認知的問題，之前是沒有的，本申請首先提出來了，并且有實例證明。在實際應用中具有明顯的預測效果。盡管本申請內容中提到線性疊加的計算公式、標準氣泡比值計算方法、判斷氣泡壓制效果差的參數取值等的公式定義、取值、判斷、交互、計算等描述，但是，本申請并不局限于必須是符合標準地震勘探采集及數據處理方法、標準數學公式計算、質控條件的設置等以及實施例所描述的情況等，某些行業標準或者使用自定義方式或實施例描述的實施基礎上略加修改后的實施方案也可以實現上述實施例相同、等同或相近、或變形后可預料的實施效果。應用這些修改或變形后的數據獲取、定義、判斷、取值方式等獲取的實施例，仍然可以屬于本申請的可選實施方案范圍之內。雖然本申請提供了如實施例或流程圖所述的方法操作步驟，但基于常規或者無創造性的手段可以包括更多或者更少的操作步驟。實施例中列舉的步驟順序僅僅為眾多步驟執行順序中的一種方式，不代表唯一的執行順序。在實際中的裝置或客戶端產品執行時，可以按照實施例或者附圖所示的方法順序執行或者并行執行(例如并行處理器或者多線程處理的環境，甚至為分布式數據處理環境)。術語“包括”、“包含”或者其任何其他變體意在涵蓋非排他性的包含，從而使得包括一系列要素的過程、方法、產品或者設備不僅包括那些要素，而且還包括沒有明確列出的其他要素，或者是還包括為這種過程、方法、產品或者設備所固有的要素。在沒有更多限制的情況下，并不排除在包括所述要素的過程、方法、產品或者設備中還存在另外的相同或等同要素。第一，第二等詞語用來表示名稱，而并不表示任何特定的順序。上述實施例闡明的單元、裝置或模塊等，具體可以由計算機芯片或實體實現，或者由具有某種功能的產品來實現。為了描述的方便，描述以上裝置時以功能分為各種模塊分別描述。當然，在實施本申請時可以把各模塊的功能在同一個或多個軟件和/或硬件中實現，也可以將實現同一功能的模塊由多個子模塊或子單元的組合實現等。以上所描述的裝置實施例僅僅是示意性的，例如，所述單元的劃分，僅僅為一種邏輯功能劃分，實際實現時可以有另外的劃分方式，例如多個單元或組件可以結合或者可以集成到另一個系統，或一些特征可以忽略，或不執行。另一點，所顯示或討論的相互之間的耦合或直接耦合或通信連接可以是通過一些接口，裝置或單元的間接耦合或通信連接，可以是電性，機械或其它的形式。本領域技術人員也知道，除了以純計算機可讀程序代碼方式實現控制器以外，完全可以通過將方法步驟進行邏輯編程來使得控制器以邏輯門、開關、專用集成電路、可編程邏輯控制器和嵌入微控制器等的形式來實現相同功能。因此這種控制器可以被認為是一種硬件部件，而對其內部包括的用于實現各種功能的裝置也可以視為硬件部件內的結構。或者甚至，可以將用于實現各種功能的裝置視為既可以是實現方法的軟件模塊又可以是硬件部件內的結構。本申請可以在由計算機執行的計算機可執行指令的一般上下文中描述，例如程序模塊。一般地，程序模塊包括執行特定任務或實現特定抽象數據類型的例程、程序、對象、組件、數據結構、類等等。也可以在分布式計算環境中實踐本申請，在這些分布式計算環境中，由通過通信網絡而被連接的遠程處理設備來執行任務。在分布式計算環境中，程序模塊可以位于包括存儲設備在內的本地和遠程計算機存儲介質中。通過以上的實施方式的描述可知，本領域的技術人員可以清楚地了解到本申請可借助軟件加必需的通用硬件平臺的方式來實現。基于這樣的理解，本申請的技術方案本質上或者說對現有技術做出貢獻的部分可以以軟件產品的形式體現出來，該計算機軟件產品可以存儲在存儲介質中，如rom/ram、磁碟、光盤等，包括若干指令用以使得一臺計算機設備(可以是個人計算機，移動終端，服務器，或者網絡設備等)執行本申請各個實施例或者實施例的某些部分所述的方法。本說明書中的各個實施例采用遞進的方式描述，各個實施例之間相同或相似的部分互相參見即可，每個實施例重點說明的都是與其他實施例的不同之處。本申請可用于眾多通用或專用的計算機系統環境或配置中。例如：個人計算機、服務器計算機、手持設備或便攜式設備、平板型設備、多處理器系統、基于微處理器的系統、置頂盒、可編程的電子設備、網絡pc、小型計算機、大型計算機、包括以上任何系統或設備的分布式計算環境等等。雖然通過實施例描繪了本申請，本領域普通技術人員知道，本申請有許多變形和變化而不脫離本申請的精神，希望所附的權利要求包括這些變形和變化而不脫離本申請的精神。當前第1頁12