本發(fā)明涉及地質(zhì)勘探技術(shù)領(lǐng)域，特別涉及一種地震資料振幅調(diào)整方法和裝置。

背景技術(shù)：

目前，由于各個(gè)區(qū)塊間的施工因素、激發(fā)以及接收條件和信號(hào)采集時(shí)間均不相同，因而區(qū)塊間采集到的參數(shù)差異較大，例如：覆蓋次數(shù)、單炮能量強(qiáng)弱等參數(shù)。當(dāng)利用各個(gè)區(qū)塊的地震資料進(jìn)行連片成像處理時(shí)，在區(qū)塊的拼接區(qū)，由于是由兩塊或者兩塊以上的地震資料疊加得到的，疊加處地震資料的覆蓋次數(shù)明顯高于其他地方，如果不能進(jìn)行有效的振幅調(diào)整，必然會(huì)導(dǎo)致劃弧現(xiàn)象；在區(qū)塊的邊緣區(qū)域，由于采集原因，地震道集的覆蓋次數(shù)明顯低于滿覆蓋地區(qū)，也必將導(dǎo)致其他地區(qū)的強(qiáng)能量往區(qū)塊的邊緣劃弧。

通常，對(duì)上述連片區(qū)塊進(jìn)行成像處理時(shí)，可以先采用幾何擴(kuò)散補(bǔ)償和地表一致性振幅補(bǔ)償?shù)确椒▽?duì)連片的地震資料進(jìn)行振幅調(diào)整。然而，采用上述方法時(shí)，由于各個(gè)區(qū)塊內(nèi)以及各個(gè)區(qū)塊間原始地震信息的振幅能量差異過(guò)大，遠(yuǎn)、近偏移距能量較弱，因而進(jìn)行調(diào)整后的地震資料效果較差。

針對(duì)上述問(wèn)題，目前尚未提出有效的解決方案。

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

本發(fā)明提供了一種地震資料振幅調(diào)整方法和裝置，以實(shí)現(xiàn)在對(duì)連片區(qū)塊進(jìn)行成像處理時(shí)，在保證地震資料信息不丟失的情況下，各個(gè)區(qū)塊間能量一致的目的。

本發(fā)明實(shí)施例提供了一種地震資料振幅調(diào)整方法，可以包括：獲取待測(cè)道集的信噪比；對(duì)所述待測(cè)道集進(jìn)行動(dòng)校正拉伸切除處理，并確定所述待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)；根據(jù)所述待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)，確定出多個(gè)時(shí)間段；統(tǒng)計(jì)切除后的待測(cè)道集落入所述多個(gè)時(shí)間段中各個(gè)時(shí)間段的地震道數(shù)；根據(jù)所述信噪比和所述各個(gè)時(shí)間段的地震道數(shù)，計(jì)算得到所述待測(cè)道集在所述各個(gè)時(shí)間段的能量調(diào)整因子；基于所述能量調(diào)整因子，對(duì)所述待測(cè)道集進(jìn)行振幅調(diào)整，得到調(diào)整后的地震信號(hào)。

在一個(gè)實(shí)施例中，可以按照以下公式根據(jù)所述信噪比和所述各個(gè)時(shí)間段的地震道數(shù)，計(jì)算得到所述待測(cè)道集在所述各個(gè)時(shí)間段的能量調(diào)整因子：

cof(t)＝p*cov(t)^-1/k+(1-p)*snr^-1/l

其中，cof(t)表示所述待測(cè)道集在時(shí)間段t的能量調(diào)整因子，cov(t)表示所述待測(cè)道集在時(shí)間段t的地震道數(shù)，snr表示所述信噪比，0≤p≤1，k＞1，l＞0。

在一個(gè)實(shí)施例中，在獲取待測(cè)道集的信噪比之前，所述方法還可以包括：獲取待測(cè)工區(qū)的地震數(shù)據(jù)；對(duì)所述地震數(shù)據(jù)進(jìn)行預(yù)處理，得到處理后的地震數(shù)據(jù)，其中，所述預(yù)處理可以包括但不限于以下至少之一：疊前去噪、地表一致性振幅補(bǔ)償、幾何擴(kuò)散吸收補(bǔ)償和道分選；根據(jù)所述處理后的地震數(shù)據(jù)，確定多個(gè)所述待測(cè)道集。

在一個(gè)實(shí)施例中，在根據(jù)所述信噪比和所述各個(gè)時(shí)間段的地震道數(shù)，計(jì)算得到所述待測(cè)道集在所述各個(gè)時(shí)間段的能量調(diào)整因子之后，所述方法還可以包括：對(duì)所述待測(cè)道集的能量調(diào)整因子進(jìn)行平滑處理，得到所述待測(cè)道集平滑處理后的能量調(diào)整因子。

在一個(gè)實(shí)施例中，可以按照以下公式基于所述能量調(diào)整因子，對(duì)所述待測(cè)道集進(jìn)行振幅調(diào)整，得到調(diào)整后的地震信號(hào)：

s'(t)＝s(t)·pcof(t)

其中，s'(t)表示調(diào)整后的在時(shí)間段t的地震信號(hào)，s(t)表示所述待測(cè)道集調(diào)整前在時(shí)間段t的地震信號(hào)，pcof(t)表示所述待測(cè)道集在時(shí)間段t進(jìn)行平滑處理后的能量調(diào)整因子。

在一個(gè)實(shí)施例中，對(duì)所述待測(cè)道集的能量調(diào)整因子進(jìn)行平滑處理，得到所述待測(cè)道集平滑處理后的能量調(diào)整因子，可以包括：按照預(yù)定的順序依次對(duì)所述待測(cè)工區(qū)中各個(gè)待測(cè)道集的位置進(jìn)行標(biāo)定，得到進(jìn)行位置標(biāo)定后的各個(gè)待測(cè)道集；按照以下方式獲取各個(gè)待測(cè)道集平滑處理后的能量調(diào)整因子：以當(dāng)前待測(cè)道集的位置為中心，預(yù)定的平滑半徑為半徑，得到平滑區(qū)域；根據(jù)所述平滑區(qū)域中各個(gè)待測(cè)道集的能量調(diào)整因子，計(jì)算得到所述當(dāng)前待測(cè)道集平滑處理后的能量調(diào)整因子。

在一個(gè)實(shí)施例中，可以按照以下公式計(jì)算得到所述當(dāng)前待測(cè)道集平滑處理后的能量調(diào)整因子：

其中，當(dāng)時(shí)，wij＝0；

當(dāng)時(shí)，

其中，pcofkl表示位置為(k,l)的所述當(dāng)前待測(cè)道集進(jìn)行平滑處理后的能量調(diào)整因子，cofij表示位置為(i,j)的待測(cè)道集進(jìn)行平滑處理前的能量調(diào)整因子，wij表示位置為(i,j)的待測(cè)道集的平滑系數(shù)，r表示所述平滑半徑，r≥1，i＝k-r,k-r+1,…,k+r，j＝l-r,l-r+1,…,l+r。

在一個(gè)實(shí)施例中，所述平滑半徑可以是根據(jù)所述待測(cè)工區(qū)中待測(cè)道集的個(gè)數(shù)確定的。

在一個(gè)實(shí)施例中，可以按照以下公式獲取待測(cè)道集的信噪比：

其中，

其中，snr表示所述待測(cè)道集的信噪比，r(xc,xd)表示所述待測(cè)道集中第c個(gè)地震道與第d個(gè)地震道的互相關(guān)值，c＝1,2…,n，d＝1,2…,n，n表示所述第c個(gè)地震道分析時(shí)窗內(nèi)的總樣點(diǎn)數(shù)。

在一個(gè)實(shí)施例中，可以按照以下方式確定所述待測(cè)道集中第c個(gè)地震道與第d個(gè)地震道的互相關(guān)值：按照以下公式確定所述待測(cè)道集中第c個(gè)地震道與第d個(gè)地震道對(duì)應(yīng)于地震道分析時(shí)窗內(nèi)的各個(gè)樣點(diǎn)序號(hào)的互相關(guān)值，得到與所述各個(gè)樣點(diǎn)序號(hào)相對(duì)應(yīng)的多個(gè)互相關(guān)結(jié)果：

其中，t＝0,1,…,l，l表示所述第c個(gè)地震道分析時(shí)窗內(nèi)的總樣點(diǎn)數(shù)，τ表示地震道分析時(shí)窗內(nèi)的樣點(diǎn)序號(hào)，τ＝0,1,…,l；

選取所述多個(gè)互相關(guān)結(jié)果中的最大值作為所述第c個(gè)地震道與第d個(gè)地震道的互相關(guān)值。

本發(fā)明實(shí)施例還提供了一種地震資料振幅調(diào)整裝置，可以包括：信噪比獲取模塊，用于獲取待測(cè)道集的信噪比；道集切除模塊，用于對(duì)所述待測(cè)道集進(jìn)行動(dòng)校正拉伸切除處理，并確定所述待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)；時(shí)間段確定模塊，用于根據(jù)所述待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)，確定出多個(gè)時(shí)間段；地震道數(shù)統(tǒng)計(jì)模塊，用于統(tǒng)計(jì)切除后的待測(cè)道集落入所述多個(gè)時(shí)間段中各個(gè)時(shí)間段的地震道數(shù)；調(diào)整因子計(jì)算模塊，用于根據(jù)所述信噪比和所述各個(gè)時(shí)間段的地震道數(shù)，計(jì)算得到所述待測(cè)道集在所述各個(gè)時(shí)間段的能量調(diào)整因子；振幅調(diào)整模塊，用于基于所述能量調(diào)整因子，對(duì)所述待測(cè)道集進(jìn)行振幅調(diào)整，得到調(diào)整后的地震信號(hào)。

在本發(fā)明實(shí)施例中，由于動(dòng)校正拉伸切除處理可以有效切除待測(cè)道集中由于偏移距變化，而失去波動(dòng)特征的噪聲信號(hào)，因此，可以先對(duì)待測(cè)道集進(jìn)行動(dòng)校正拉伸切除處理，并確定待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)，通過(guò)統(tǒng)計(jì)切除后的待測(cè)道集落入由切除時(shí)間點(diǎn)所確定的多個(gè)時(shí)間段中各個(gè)時(shí)間段的地震道數(shù)，從而考慮了由于待測(cè)道集的偏移距變化引起的信號(hào)能量分布不均的問(wèn)題。進(jìn)一步地，利用上述多個(gè)地震道數(shù)以及信噪比所確定的待測(cè)道集的能量調(diào)整因子進(jìn)行振幅調(diào)整，可以在連片區(qū)塊的拼接區(qū)和邊緣區(qū)得到更加符合地質(zhì)規(guī)律的勘探信號(hào)，應(yīng)用該勘探信號(hào)對(duì)連片區(qū)塊進(jìn)行成像時(shí)，可以得到質(zhì)量更高的地質(zhì)成像結(jié)果，從而實(shí)現(xiàn)在保證地震資料信息不丟失的情況下，各個(gè)區(qū)塊間能量的一致。

附圖說(shuō)明

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

圖1是本申請(qǐng)?zhí)峁┑囊环N地震資料振幅調(diào)整方法流程圖；

圖2是基于常規(guī)振幅調(diào)整方法對(duì)某地區(qū)進(jìn)行振幅調(diào)整后地震信號(hào)的偏移成像示意圖；

圖3是基于本申請(qǐng)?zhí)峁┑恼穹{(diào)整方法對(duì)某地區(qū)進(jìn)行振幅調(diào)整后地震信號(hào)的偏移成像示意圖；

圖4是本申請(qǐng)?zhí)峁┑囊环N地震資料振幅調(diào)整裝置的結(jié)構(gòu)框圖。

具體實(shí)施方式

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

考慮到現(xiàn)有技術(shù)中各個(gè)地震道集內(nèi)不同偏移距變化所引起的信號(hào)能量分布不均的問(wèn)題，而動(dòng)校正拉伸切除處理可以切除待測(cè)道集中由于偏移距變化所產(chǎn)生的噪聲信號(hào)，發(fā)明人提出了將動(dòng)校正拉伸切除處理和多次統(tǒng)計(jì)地震道數(shù)相結(jié)合的方式，先統(tǒng)計(jì)切除后的待測(cè)道集在由切除時(shí)間點(diǎn)所確定的多個(gè)時(shí)間段中各個(gè)時(shí)間段的地震道數(shù)，再利用各個(gè)地震道數(shù)確定待測(cè)道集的能量調(diào)整因子，最后根據(jù)該能量調(diào)整因子對(duì)待測(cè)道集進(jìn)行振幅調(diào)整?；诖耍岢隽艘环N地震資料振幅調(diào)整方法，如圖1所示，可以包括以下步驟：

s101：獲取待測(cè)道集的信噪比。

其中，待測(cè)道集可以是共中心點(diǎn)道集、共炮點(diǎn)道集、共接收點(diǎn)道集、共反射點(diǎn)道集、共深度點(diǎn)道集等。

在本申請(qǐng)的一個(gè)實(shí)施例中，在獲取待測(cè)道集的信噪比之前，還可以包括：獲取待測(cè)工區(qū)的地震數(shù)據(jù)；對(duì)地震數(shù)據(jù)進(jìn)行預(yù)處理，得到處理后的地震數(shù)據(jù)，其中，預(yù)處理可以包括但不限于以下至少之一：疊前去噪、地表一致性振幅補(bǔ)償、幾何擴(kuò)散吸收補(bǔ)償和道分選；根據(jù)處理后的地震數(shù)據(jù)，確定多個(gè)待測(cè)道集。

采用上述預(yù)處理方式對(duì)待測(cè)道集進(jìn)行去噪處理之后，可以去除待測(cè)道集中較為常見(jiàn)的噪聲信號(hào)，從而可以得到更加符合地質(zhì)特征的待測(cè)道集。

在本申請(qǐng)的一個(gè)實(shí)施例中，可以按照以下公式獲取待測(cè)道集的信噪比：

其中，

其中，snr表示待測(cè)道集的信噪比，r(xc,xd)表示待測(cè)道集中第c個(gè)地震道與第d個(gè)地震道的互相關(guān)值，c＝1,2…,n，d＝1,2…,n，n表示第c個(gè)地震道分析時(shí)窗內(nèi)的總樣點(diǎn)數(shù)。

在確定待測(cè)道集的信噪比的過(guò)程中，可以先通過(guò)以下方式確定待測(cè)道集中第c個(gè)地震道與第d個(gè)地震道的互相關(guān)值：

s1-1：按照以下公式確定待測(cè)道集中第c個(gè)地震道與第d個(gè)地震道對(duì)應(yīng)于地震道分析時(shí)窗內(nèi)的各個(gè)樣點(diǎn)序號(hào)的互相關(guān)值，得到與各個(gè)樣點(diǎn)序號(hào)相對(duì)應(yīng)的多個(gè)互相關(guān)結(jié)果：

其中，t＝0,1,…,l，l表示第c個(gè)地震道分析時(shí)窗內(nèi)的總樣點(diǎn)數(shù)，τ表示地震道分析時(shí)窗內(nèi)的樣點(diǎn)序號(hào)，τ＝0,1,…,l；

s1-2：選取多個(gè)互相關(guān)結(jié)果中的最大值作為第c個(gè)地震道與第d個(gè)地震道的互相關(guān)值。

將得到的多個(gè)互相關(guān)值代入信噪比計(jì)算公式中，從而可以計(jì)算得到該待測(cè)道集的信噪比。

s102：對(duì)待測(cè)道集進(jìn)行動(dòng)校正拉伸切除處理，并確定所述待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)。

s103：根據(jù)所述待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)，確定出多個(gè)時(shí)間段。

s104：統(tǒng)計(jì)切除后的待測(cè)道集落入所述多個(gè)時(shí)間段中各個(gè)時(shí)間段的地震道數(shù)。

地震資料動(dòng)校正關(guān)系到地震資料處理結(jié)果的精度和應(yīng)用效果。動(dòng)校正過(guò)程中，信號(hào)時(shí)間會(huì)拉長(zhǎng)、頻率向低頻移動(dòng)且會(huì)產(chǎn)生波形變異的現(xiàn)象，拉伸畸變嚴(yán)重時(shí)，不僅會(huì)破壞疊加效果，還會(huì)使得地震信號(hào)失去應(yīng)有的波動(dòng)力學(xué)特征，降低分辨率，從而會(huì)影響地震勘探的精度和地質(zhì)效果。因此，我們可以對(duì)動(dòng)校正拉伸之后的波形進(jìn)行切除處理，從而可以實(shí)現(xiàn)消除上述影響的目的。

應(yīng)用到本申請(qǐng)中，即，可以通過(guò)對(duì)待測(cè)道集進(jìn)行動(dòng)校正拉伸切除處理，并確定待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)，根據(jù)待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)，確定出多個(gè)時(shí)間段，統(tǒng)計(jì)切除后的待測(cè)道集落入所述多個(gè)時(shí)間段中各個(gè)時(shí)間段的地震道數(shù)。

在實(shí)際應(yīng)用的過(guò)程中，可以通過(guò)讀取每個(gè)地震道的道頭中的信息，從而獲取各個(gè)地震道的切除時(shí)間，再對(duì)應(yīng)到待測(cè)道集中進(jìn)行地震道數(shù)的統(tǒng)計(jì)。通過(guò)統(tǒng)計(jì)由切除時(shí)間點(diǎn)確定的多個(gè)時(shí)間段中各個(gè)時(shí)間段的切除處理后的地震道數(shù)，從而考慮了待測(cè)道集中由于偏移距變化對(duì)待測(cè)道集的影響，提高了連片區(qū)塊進(jìn)行成像處理時(shí)的效果。

進(jìn)一步的，在本申請(qǐng)中，待測(cè)道集中所包含的地震道數(shù)即為待測(cè)道集的覆蓋次數(shù)。因而，上述統(tǒng)計(jì)地震道數(shù)的方式即為統(tǒng)計(jì)待測(cè)道集的覆蓋次數(shù)。因而，對(duì)于單個(gè)待測(cè)道集而言，采用上述方式統(tǒng)計(jì)到的覆蓋次數(shù)是隨切除時(shí)間所確定的時(shí)間段變化的。例如：某待測(cè)道集具有10個(gè)地震道，進(jìn)行動(dòng)校正拉伸切除處理后，在0-5s，10個(gè)地震道均沒(méi)有被切除；在5s-7s，8個(gè)地震道沒(méi)有被切除；在7s-8s，5個(gè)地震道沒(méi)有被切除；在8s-9s，3個(gè)地震道沒(méi)有被切除；在9s-10s，1個(gè)地震道沒(méi)有被切除。對(duì)于該例子中的待測(cè)道集而言，在0-5s，覆蓋次數(shù)為10；在5s-7s，覆蓋次數(shù)為8；在7s-8s，覆蓋次數(shù)為5；在8s-9s，覆蓋次數(shù)為3；在9s-10s，覆蓋次數(shù)為1。

s105：根據(jù)所述信噪比和所述各個(gè)時(shí)間段的地震道數(shù)，計(jì)算得到所述待測(cè)道集在所述各個(gè)時(shí)間段的能量調(diào)整因子。

在采用s101至s105所述方法分布得到待測(cè)道集的信噪比以及隨時(shí)間段變化的地震道數(shù)之后，可以按照以下公式根據(jù)上述信噪比和上述各個(gè)時(shí)間段的地震道數(shù)，計(jì)算得到待測(cè)道集在各個(gè)時(shí)間段的能量調(diào)整因子：

cof(t)＝p*cov(t)^-1/k+(1-p)*snr^-1/l

其中，cof(t)表示待測(cè)道集在時(shí)間段t的能量調(diào)整因子，cov(t)表示待測(cè)道集在時(shí)間段t的地震道數(shù)，snr表示信噪比，0≤p≤1，k＞1，l＞0。

對(duì)于一個(gè)待測(cè)道集而言，覆蓋次數(shù)是隨時(shí)間段變化的，但是信噪比是不變的。

在本申請(qǐng)的一個(gè)實(shí)施例中，在得到上述能量調(diào)整因子之后，可以對(duì)待測(cè)道集的能量調(diào)整因子進(jìn)行平滑處理，得到待測(cè)道集平滑處理后的能量調(diào)整因子。進(jìn)行平滑處理之前，待測(cè)工區(qū)中各個(gè)待測(cè)道集已經(jīng)先計(jì)算得到了與各個(gè)待測(cè)道集對(duì)應(yīng)的能量調(diào)整因子。當(dāng)進(jìn)行平滑處理時(shí)，是對(duì)與各個(gè)待測(cè)道集對(duì)應(yīng)的能量調(diào)整因子進(jìn)行平滑運(yùn)算的。具體的，可以包括以下步驟：

s5-1：按照預(yù)定的順序依次對(duì)待測(cè)工區(qū)中各個(gè)待測(cè)道集的位置進(jìn)行標(biāo)定，得到進(jìn)行位置標(biāo)定后的各個(gè)待測(cè)道集；

s5-2：按照以下方式獲取各個(gè)待測(cè)道集平滑處理后的能量調(diào)整因子：以當(dāng)前待測(cè)道集的位置為中心，預(yù)定的平滑半徑為半徑，得到平滑區(qū)域；根據(jù)平滑區(qū)域中各個(gè)待測(cè)道集的能量調(diào)整因子，計(jì)算得到當(dāng)前待測(cè)道集平滑處理后的能量調(diào)整因子。

即，可以按照以下公式計(jì)算得到上述當(dāng)前待測(cè)道集平滑處理后的能量調(diào)整因子：

其中，當(dāng)時(shí)，wij＝0；

當(dāng)時(shí)，

其中，pcofkl表示位置為(k,l)的當(dāng)前待測(cè)道集進(jìn)行平滑處理后的能量調(diào)整因子，cofij表示位置為(i,j)的待測(cè)道集進(jìn)行平滑處理前的能量調(diào)整因子，wij表示位置為(i,j)的待測(cè)道集的平滑系數(shù)，r表示平滑半徑，i＝k-r,k-r+1,…,k+r，j＝l-r,l-r+1,…,l+r。

其中，平滑半徑可以是根據(jù)待測(cè)工區(qū)中待測(cè)道集的個(gè)數(shù)確定的。即，當(dāng)待測(cè)工區(qū)中有10個(gè)待測(cè)道集時(shí)，平滑半徑可以是10。

s106：基于所述能量調(diào)整因子，對(duì)所述待測(cè)道集進(jìn)行振幅調(diào)整，得到調(diào)整后的地震信號(hào)。

在本申請(qǐng)的一個(gè)實(shí)施例中，在得到能量調(diào)整因子或者平滑處理后的能量調(diào)整因子之后，可以按照以下公式基于上述能量調(diào)整因子或者平滑處理后的能量調(diào)整因子，對(duì)待測(cè)道集進(jìn)行振幅調(diào)整，得到調(diào)整后的地震信號(hào)：

s'(t)＝s(t)·pcof(t)

其中，s'(t)表示調(diào)整后的在時(shí)間段t的地震信號(hào)，s(t)表示待測(cè)道集調(diào)整前在時(shí)間段t的地震信號(hào)，pcof(t)表示待測(cè)道集在時(shí)間段t進(jìn)行平滑處理時(shí)的能量調(diào)整因子。

通過(guò)對(duì)地震資料進(jìn)行振幅調(diào)整，使得面元道集總體能量大致保持在同一水平，進(jìn)而獲得覆蓋次數(shù)較均勻的疊前道集。

如圖2所示為基于常規(guī)振幅調(diào)整方法對(duì)某地區(qū)進(jìn)行振幅調(diào)整后地震信號(hào)的偏移成像示意圖，圖3所示為本發(fā)明實(shí)施例的振幅調(diào)整后地震信號(hào)的偏移成像示意圖，對(duì)比圖2和圖3可以發(fā)現(xiàn)：在圖3中的6個(gè)白色箭頭處可以明顯看到6個(gè)斷面，即，圖3可以較好地顯示地層的細(xì)節(jié)。由此可見(jiàn)，通過(guò)本例所提供的振幅調(diào)整方法可以有效解決疊前偏移畫弧及邊界效應(yīng)問(wèn)題，有效地提高了偏移剖面信噪比，使地震數(shù)據(jù)解釋人員能夠更加準(zhǔn)確地判斷和解釋地下地質(zhì)構(gòu)造情況，提高了尋找有利油氣圈閉以及確定井位的精度。

基于同一發(fā)明構(gòu)思，本發(fā)明實(shí)施例中還提供了一種地震資料振幅調(diào)整裝置，如下面的實(shí)施例所述。由于地震資料振幅調(diào)整裝置解決問(wèn)題的原理與地震資料振幅調(diào)整方法相似，因此地震資料振幅調(diào)整裝置的實(shí)施可以參見(jiàn)地震資料振幅調(diào)整方法的實(shí)施，重復(fù)之處不再贅述。以下所使用的，術(shù)語(yǔ)“單元”或者“模塊”可以實(shí)現(xiàn)預(yù)定功能的軟件和/或硬件的組合。盡管以下實(shí)施例所描述的裝置較佳地以軟件來(lái)實(shí)現(xiàn)，但是硬件，或者軟件和硬件的組合的實(shí)現(xiàn)也是可能并被構(gòu)想的。圖4是本發(fā)明實(shí)施例的地震資料振幅調(diào)整裝置的一種結(jié)構(gòu)框圖，如圖4所示，可以包括：信噪比獲取模塊401、道集切除模塊402、時(shí)間段確定模塊403、地震道數(shù)統(tǒng)計(jì)模塊404、調(diào)整因子計(jì)算模塊405、振幅調(diào)整模塊406，下面對(duì)該結(jié)構(gòu)進(jìn)行說(shuō)明。

信噪比獲取模塊401，可以用于獲取待測(cè)道集的信噪比；

道集切除模塊402，可以用于對(duì)所述待測(cè)道集進(jìn)行動(dòng)校正拉伸切除處理，并確定所述待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)；

時(shí)間段確定模塊403，可以用于根據(jù)所述待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)，確定出多個(gè)時(shí)間段；

地震道數(shù)統(tǒng)計(jì)模塊404，可以用于統(tǒng)計(jì)切除后的待測(cè)道集落入所述多個(gè)時(shí)間段中各個(gè)時(shí)間段的地震道數(shù)；

調(diào)整因子計(jì)算模塊405，可以用于根據(jù)所述信噪比和所述各個(gè)時(shí)間段的地震道數(shù)，計(jì)算得到所述待測(cè)道集在所述各個(gè)時(shí)間段的能量調(diào)整因子；

振幅調(diào)整模塊406，可以用于基于所述能量調(diào)整因子，對(duì)所述待測(cè)道集進(jìn)行振幅調(diào)整，得到調(diào)整后的地震信號(hào)。

在本申請(qǐng)的一個(gè)實(shí)施例中，所述調(diào)整因子計(jì)算模塊具體可以用于按照以下公式根據(jù)所述信噪比和所述各個(gè)時(shí)間段的地震道數(shù)，計(jì)算得到所述待測(cè)道集在所述各個(gè)時(shí)間段的能量調(diào)整因子：

cof(t)＝p*cov(t)^-1/k+(1-p)*snr^-1/l

其中，cof(t)表示所述待測(cè)道集在時(shí)間段t的能量調(diào)整因子，cov(t)表示所述待測(cè)道集在時(shí)間段t的地震道數(shù)，snr表示所述信噪比，0≤p≤1，k＞1，l＞0。

在本申請(qǐng)的一個(gè)實(shí)施例中，所述信噪比獲取模塊可以包括：地震數(shù)據(jù)獲取單元，可以用于在獲取待測(cè)道集的信噪比之前，獲取待測(cè)工區(qū)的地震數(shù)據(jù)；預(yù)處理單元，可以用于對(duì)所述地震數(shù)據(jù)進(jìn)行預(yù)處理，得到處理后的地震數(shù)據(jù)，其中，所述預(yù)處理可以包括但不限于以下至少之一：疊前去噪、地表一致性振幅補(bǔ)償、幾何擴(kuò)散吸收補(bǔ)償和道分選；道集確定單元，可以用于根據(jù)所述處理后的地震數(shù)據(jù)，確定多個(gè)所述待測(cè)道集。

在本申請(qǐng)的一個(gè)實(shí)施例中，所述調(diào)整因子計(jì)算模塊在根據(jù)所述信噪比和所述各個(gè)時(shí)間段的地震道數(shù)，計(jì)算得到所述待測(cè)道集在所述各個(gè)時(shí)間段的能量調(diào)整因子之后，還可以用于對(duì)所述待測(cè)道集的能量調(diào)整因子進(jìn)行平滑處理，得到所述待測(cè)道集平滑處理后的能量調(diào)整因子。

在本申請(qǐng)的一個(gè)實(shí)施例中，所述振幅調(diào)整模塊具體可以用于按照以下公式基于所述能量調(diào)整因子，對(duì)所述待測(cè)道集進(jìn)行振幅調(diào)整，得到調(diào)整后的地震信號(hào)：

s'(t)＝s(t)·pcof(t)

其中，s'(t)表示調(diào)整后的在時(shí)間段t的地震信號(hào)，s(t)表示所述待測(cè)道集調(diào)整前在時(shí)間段t的地震信號(hào)，pcof(t)表示所述待測(cè)道集在時(shí)間段t進(jìn)行平滑處理后的能量調(diào)整因子。

在本申請(qǐng)的一個(gè)實(shí)施例中，所述調(diào)整因子計(jì)算模塊可以包括：位置標(biāo)定單元，可以用于按照預(yù)定的順序依次對(duì)所述待測(cè)工區(qū)中各個(gè)待測(cè)道集的位置進(jìn)行標(biāo)定，得到進(jìn)行位置標(biāo)定后的各個(gè)待測(cè)道集；因子計(jì)算單元，可以用于按照以下方式獲取各個(gè)待測(cè)道集平滑處理后的能量調(diào)整因子：以當(dāng)前待測(cè)道集的位置為中心，預(yù)定的平滑半徑為半徑，得到平滑區(qū)域；根據(jù)所述平滑區(qū)域中各個(gè)待測(cè)道集的能量調(diào)整因子，計(jì)算得到所述當(dāng)前待測(cè)道集平滑處理后的能量調(diào)整因子。

在本申請(qǐng)的一個(gè)實(shí)施例中，所述因子計(jì)算單元具體可以用于按照以下公式計(jì)算得到所述當(dāng)前待測(cè)道集平滑處理后的能量調(diào)整因子：

其中，當(dāng)時(shí)，wij＝0；

當(dāng)時(shí)，

其中，pcofkl表示位置為(k,l)的所述當(dāng)前待測(cè)道集進(jìn)行平滑處理后的能量調(diào)整因子，cofij表示位置為(i,j)的待測(cè)道集進(jìn)行平滑處理前的能量調(diào)整因子，wij表示位置為(i,j)的待測(cè)道集的平滑系數(shù)，r表示所述平滑半徑，i＝k-r,k-r+1,…,k+r，j＝l-r,l-r+1,…,l+r。

在本申請(qǐng)的一個(gè)實(shí)施例中，所述平滑半徑可以是根據(jù)所述待測(cè)工區(qū)中待測(cè)道集的個(gè)數(shù)確定的。

在本申請(qǐng)的一個(gè)實(shí)施例中，所述信噪比獲取模塊具體可以用于按照以下公式獲取待測(cè)道集的信噪比：

其中，

其中，snr表示所述待測(cè)道集的信噪比，r(xc,xd)表示所述待測(cè)道集中第c個(gè)地震道與第d個(gè)地震道的互相關(guān)值，c＝1,2…,n，d＝1,2…,n，n表示所述第c個(gè)地震道分析時(shí)窗內(nèi)的總樣點(diǎn)數(shù)。

在本申請(qǐng)的一個(gè)實(shí)施例中，所述信噪比獲取模塊具體可以用于按照以下方式確定所述待測(cè)道集中第c個(gè)地震道與第d個(gè)地震道的互相關(guān)值：按照以下公式確定所述待測(cè)道集中第c個(gè)地震道與第d個(gè)地震道對(duì)應(yīng)于地震道分析時(shí)窗內(nèi)的各個(gè)樣點(diǎn)序號(hào)的互相關(guān)值，得到與所述各個(gè)樣點(diǎn)序號(hào)相對(duì)應(yīng)的多個(gè)互相關(guān)結(jié)果：

其中，t＝0,1,…,l，l表示所述第c個(gè)地震道分析時(shí)窗內(nèi)的總樣點(diǎn)數(shù)，τ表示地震道分析時(shí)窗內(nèi)的樣點(diǎn)序號(hào)，τ＝0,1,…,l；

選取所述多個(gè)互相關(guān)結(jié)果中的最大值作為所述第c個(gè)地震道與第d個(gè)地震道的互相關(guān)值。

利用上述各實(shí)施例所提供的地震資料振幅調(diào)整裝置的實(shí)施方式，可以自動(dòng)實(shí)施所述地震資料振幅調(diào)整方法，對(duì)所述待測(cè)道集進(jìn)行振幅調(diào)整，可以不需要實(shí)施人員的具體參與，可以直接輸出振幅調(diào)整結(jié)果，操作簡(jiǎn)單快捷，有效提高了用戶體驗(yàn)。

所述地震資料振幅調(diào)整裝置中，所述能量調(diào)整因子計(jì)算、平滑處理、信噪比獲取的實(shí)施方式的擴(kuò)展可以參照前述方法的相關(guān)描述。

從以上的描述中，可以看出，本發(fā)明實(shí)施例實(shí)現(xiàn)了如下技術(shù)效果：由于動(dòng)校正拉伸切除處理可以有效切除待測(cè)道集中由于偏移距變化，而失去波動(dòng)特征的噪聲信號(hào)，因此，可以先對(duì)待測(cè)道集進(jìn)行動(dòng)校正拉伸切除處理，并確定待測(cè)道集中各個(gè)地震道的切除時(shí)間點(diǎn)，通過(guò)統(tǒng)計(jì)切除后的待測(cè)道集落入由切除時(shí)間點(diǎn)所確定的多個(gè)時(shí)間段中各個(gè)時(shí)間段的地震道數(shù)，從而考慮了由于待測(cè)道集的偏移距變化引起的信號(hào)能量分布不均的問(wèn)題。進(jìn)一步地，利用上述多個(gè)地震道數(shù)以及信噪比所確定的待測(cè)道集的能量調(diào)整因子進(jìn)行振幅調(diào)整，可以在連片區(qū)塊的拼接區(qū)和邊緣區(qū)得到更加符合地質(zhì)規(guī)律的勘探信號(hào)，應(yīng)用該勘探信號(hào)對(duì)連片區(qū)塊進(jìn)行成像時(shí)，可以得到質(zhì)量更高的地質(zhì)成像結(jié)果，從而實(shí)現(xiàn)在保證地震資料信息不丟失的情況下，各個(gè)區(qū)塊間能量的一致。

盡管本申請(qǐng)內(nèi)容中提到信噪比獲取、能量調(diào)整因子計(jì)算、平滑處理等描述，但是，本申請(qǐng)并不局限于必須是本申請(qǐng)實(shí)施例所描述的情況。某些行業(yè)標(biāo)準(zhǔn)或者使用自定義方式或?qū)嵤├枋龅膶?shí)施基礎(chǔ)上略加修改后的實(shí)施方案也可以實(shí)現(xiàn)上述實(shí)施例相同、等同或相近、或變形后可預(yù)料的實(shí)施效果。應(yīng)用這些修改或變形后的數(shù)據(jù)計(jì)算、處理等獲取的實(shí)施例，仍然可以屬于本申請(qǐng)的可選實(shí)施方案范圍之內(nèi)。

雖然本申請(qǐng)?zhí)峁┝巳鐚?shí)施例或流程圖所述的方法操作步驟，但基于常規(guī)或者無(wú)創(chuàng)造性的手段可以包括更多或者更少的操作步驟。實(shí)施例中列舉的步驟順序僅僅為眾多步驟執(zhí)行順序中的一種方式，不代表唯一的執(zhí)行順序。在實(shí)際中的裝置或終端產(chǎn)品執(zhí)行時(shí)，可以按照實(shí)施例或者附圖所示的方法順序執(zhí)行或者并行執(zhí)行(例如并行處理器或者多線程處理的環(huán)境，甚至為分布式數(shù)據(jù)處理環(huán)境)。術(shù)語(yǔ)“包括”、“包含”或者其任何其他變體意在涵蓋非排他性的包含，從而使得包括一系列要素的過(guò)程、方法、產(chǎn)品或者設(shè)備不僅包括那些要素，而且還包括沒(méi)有明確列出的其他要素，或者是還包括為這種過(guò)程、方法、產(chǎn)品或者設(shè)備所固有的要素。在沒(méi)有更多限制的情況下，并不排除在包括所述要素的過(guò)程、方法、產(chǎn)品或者設(shè)備中還存在另外的相同或等同要素。

上述實(shí)施例闡明的單元、裝置或模塊等，具體可以由計(jì)算機(jī)芯片或?qū)嶓w實(shí)現(xiàn)，或者由具有某種功能的產(chǎn)品來(lái)實(shí)現(xiàn)。為了描述的方便，描述以上裝置時(shí)以功能分為各種模塊分別描述。當(dāng)然，在實(shí)施本申請(qǐng)時(shí)可以把各模塊的功能在同一個(gè)或多個(gè)軟件和/或硬件中實(shí)現(xiàn)，也可以將實(shí)現(xiàn)同一功能的模塊由多個(gè)子模塊或子單元的組合實(shí)現(xiàn)等。以上所描述的裝置實(shí)施例僅僅是示意性的，例如，所述單元的劃分，僅僅為一種邏輯功能劃分，實(shí)際實(shí)現(xiàn)時(shí)可以有另外的劃分方式，例如多個(gè)單元或組件可以結(jié)合或者可以集成到另一個(gè)系統(tǒng)，或一些特征可以忽略，或不執(zhí)行。另一點(diǎn)，所顯示或討論的相互之間的耦合或直接耦合或通信連接可以是通過(guò)一些接口，裝置或單元的間接耦合或通信連接，可以是電性，機(jī)械或其它的形式。

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

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

通過(guò)以上的實(shí)施方式的描述可知，本領(lǐng)域的技術(shù)人員可以清楚地了解到本申請(qǐng)可借助軟件加必需的通用硬件平臺(tái)的方式來(lái)實(shí)現(xiàn)。基于這樣的理解，本申請(qǐng)的技術(shù)方案本質(zhì)上或者說(shuō)對(duì)現(xiàn)有技術(shù)做出貢獻(xiàn)的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來(lái)，該計(jì)算機(jī)軟件產(chǎn)品可以存儲(chǔ)在存儲(chǔ)介質(zhì)中，如rom/ram、磁碟、光盤等，包括若干指令用以使得一臺(tái)計(jì)算機(jī)設(shè)備(可以是個(gè)人計(jì)算機(jī)，移動(dòng)終端，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本申請(qǐng)各個(gè)實(shí)施例或者實(shí)施例的某些部分所述的方法。

本說(shuō)明書中的各個(gè)實(shí)施例采用遞進(jìn)的方式描述，各個(gè)實(shí)施例之間相同或相似的部分互相參見(jiàn)即可，每個(gè)實(shí)施例重點(diǎn)說(shuō)明的都是與其他實(shí)施例的不同之處。本申請(qǐng)可用于眾多通用或?qū)Ｓ玫挠?jì)算機(jī)系統(tǒng)環(huán)境或配置中。例如：個(gè)人計(jì)算機(jī)、服務(wù)器計(jì)算機(jī)、手持設(shè)備或便攜式設(shè)備、平板型設(shè)備、多處理器系統(tǒng)、基于微處理器的系統(tǒng)、置頂盒、可編程的電子設(shè)備、網(wǎng)絡(luò)pc、小型計(jì)算機(jī)、大型計(jì)算機(jī)、包括以上任何系統(tǒng)或設(shè)備的分布式計(jì)算環(huán)境等等。

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