本發(fā)明涉及地球物理信號(hào)處理領(lǐng)域，尤其涉及一種地震偏移數(shù)據(jù)融合的方法、裝置、設(shè)備及介質(zhì)。

背景技術(shù)：

1、為了對(duì)油氣勘探區(qū)進(jìn)行精細(xì)勘探，采用分年度、分區(qū)塊對(duì)勘探區(qū)進(jìn)行三維地震采集施工，在完成每個(gè)年度采集施工完成后進(jìn)行全流程地震處理，包括獲得單個(gè)區(qū)塊的疊前時(shí)間偏移、疊前深度偏移成果。

2、進(jìn)一步，需要對(duì)油氣勘探區(qū)進(jìn)行連片構(gòu)造解釋、儲(chǔ)層反演和油氣預(yù)測(cè)，以及對(duì)不同年度采集的三維資料進(jìn)行疊前連片處理或疊后數(shù)據(jù)拼接。其中，疊前連片處理數(shù)據(jù)量大、計(jì)算機(jī)資料消耗大、時(shí)效慢，難以及時(shí)滿足區(qū)域油氣勘探的要求；疊后拼接處理數(shù)據(jù)量小、時(shí)效快，是最佳的技術(shù)手段之一。

3、由于分年度、分區(qū)塊資料的地震處理流程和參數(shù)具有繼承性，其波形特征差異較小，滿足進(jìn)行疊后拼接處理的要求。

4、但是，地震偏移成果存在兩項(xiàng)差異導(dǎo)致偏移成果數(shù)據(jù)難以直接進(jìn)行拼接處理：不同區(qū)塊的激發(fā)接收條件差異、采集參數(shù)差異、地震能量的系統(tǒng)性差異造成不同區(qū)塊偏移成果的能量、頻率存在著一定的差異；由于偏移速度、偏移參數(shù)等差異造成不同區(qū)塊間存在著一定的成像差異，尤其在資料的邊界處常疊加缺失邊界方向的反射信息導(dǎo)致邊界處成像差異較大。

5、基于現(xiàn)有技術(shù)中存在的以上問題導(dǎo)致地震偏移數(shù)據(jù)無(wú)法直接進(jìn)行疊后拼接處理或者無(wú)法獲得一致性好的連片偏移成果。

技術(shù)實(shí)現(xiàn)思路

1、有鑒于此，本發(fā)明提出了一種地震偏移數(shù)據(jù)融合的方法、裝置、設(shè)備及介質(zhì)，解決了三維區(qū)塊之間因能量、頻率差異和邊界成像差異而導(dǎo)致疊后拼接處理成果存在能量、頻率和成像差異的問題，獲得一致性好的連片偏移成果的問題。

2、基于以上目的，本發(fā)明的實(shí)施例首先提供了一種地震偏移數(shù)據(jù)融合的方法，包括：獲取若干區(qū)塊分別對(duì)應(yīng)的地震偏移數(shù)據(jù)并基于目的層特征和地質(zhì)任務(wù)選定所述若干區(qū)塊中的一個(gè)區(qū)塊為基準(zhǔn)區(qū)塊；基于所述基準(zhǔn)區(qū)塊的絕對(duì)值振幅與所述若干區(qū)塊中的非基準(zhǔn)區(qū)塊分別對(duì)應(yīng)的絕對(duì)值振幅之間的能量差異構(gòu)建振幅校正序列，以基于所述振幅校正序列對(duì)所述地震偏移數(shù)據(jù)進(jìn)行能量校正；基于能量校正后的地震偏移數(shù)據(jù)提取所述若干區(qū)塊分別在所述目的層的相關(guān)子波，以基于所述相關(guān)子波和所述地質(zhì)任務(wù)構(gòu)建目標(biāo)子波；計(jì)算所述目標(biāo)子波與各個(gè)所述相關(guān)子波之間的頻率差異并基于所述頻率差異獲取所述相關(guān)子波重構(gòu)到所述目標(biāo)子波上的子波重構(gòu)因子，以基于所述子波重構(gòu)因子對(duì)對(duì)應(yīng)區(qū)塊的所述能量校正后的地震偏移數(shù)據(jù)進(jìn)行頻率校正；對(duì)所述若干區(qū)塊之間的重合區(qū)域的頻率校正后的地震偏移數(shù)據(jù)進(jìn)行時(shí)差校正后與所述若干區(qū)塊之間的未重合區(qū)域的頻率校正后的地震偏移數(shù)據(jù)進(jìn)行融合，得到對(duì)所述地震偏移數(shù)據(jù)的融合成果。

3、在一些實(shí)施例中，所述基于所述基準(zhǔn)區(qū)塊的絕對(duì)值振幅與所述若干區(qū)塊中的非基準(zhǔn)區(qū)塊分別對(duì)應(yīng)的絕對(duì)值振幅之間的能量差異構(gòu)建振幅校正序列，以基于所述振幅校正序列對(duì)所述地震偏移數(shù)據(jù)進(jìn)行能量校正的步驟包括：基于所述目的層對(duì)所述若干區(qū)塊分別對(duì)應(yīng)的地震偏移數(shù)據(jù)做層拉平處理，并按照預(yù)設(shè)步長(zhǎng)的時(shí)窗提取所述若干區(qū)塊分別在所述目的層對(duì)應(yīng)的絕對(duì)值振幅；計(jì)算所述非基準(zhǔn)區(qū)塊的絕對(duì)值振幅與所述基準(zhǔn)區(qū)塊的絕對(duì)值振幅之間的能量差異；基于所述能量差異構(gòu)建所述非基準(zhǔn)區(qū)塊的絕對(duì)值振幅相對(duì)于所述基準(zhǔn)區(qū)塊的絕對(duì)值振幅的振幅校正因子；將所述振幅校正因子按照時(shí)間順序進(jìn)行排列以及線性插值，得到振幅校正序列；基于所述振幅校正序列對(duì)所述層拉平處理后的地震偏移數(shù)據(jù)進(jìn)行反層拉平處理，以對(duì)所述地震偏移數(shù)據(jù)進(jìn)行能量校正。

4、在一些實(shí)施例中，所述基于所述振幅校正序列對(duì)所述地震偏移數(shù)據(jù)進(jìn)行能量校正的步驟包括：將層拉平處理后的所述非基準(zhǔn)區(qū)塊對(duì)應(yīng)的地震偏移數(shù)據(jù)在時(shí)間域上乘以所述振幅校正序列以及反層拉平處理；判斷所述能量差異是否在預(yù)設(shè)能量范圍內(nèi)，根據(jù)判斷結(jié)果確認(rèn)是否重構(gòu)所述振幅校正因子并相應(yīng)進(jìn)行重構(gòu)或保留；基于重構(gòu)/保留后的振幅校正因子，通過(guò)所述能量差異在所述預(yù)設(shè)能量差異范圍內(nèi)所對(duì)應(yīng)的振幅校正序列對(duì)所述非基準(zhǔn)區(qū)塊在所述目的層對(duì)應(yīng)的地震偏移數(shù)據(jù)進(jìn)行能量校正。

5、在一些實(shí)施例中，所述基于能量校正后的地震偏移數(shù)據(jù)提取所述若干區(qū)塊分別在所述目的層的相關(guān)子波，以基于所述相關(guān)子波和所述地質(zhì)任務(wù)構(gòu)建目標(biāo)子波的步驟包括：將能量校正后的地震偏移數(shù)據(jù)在所述目的層進(jìn)行層拉平處理后提取所述若干區(qū)塊分別對(duì)應(yīng)的相關(guān)子波；基于所述若干區(qū)塊分別對(duì)應(yīng)的相關(guān)子波獲取所述若干區(qū)塊分別對(duì)應(yīng)的有效頻帶范圍，以基于所述若干區(qū)塊分別對(duì)應(yīng)的有效頻帶范圍得到目標(biāo)頻帶范圍；基于最小二乘誤差最小方式以及所述地質(zhì)任務(wù)構(gòu)建所述目標(biāo)頻帶范圍對(duì)應(yīng)的目標(biāo)子波。

6、在一些實(shí)施例中，所述計(jì)算所述目標(biāo)子波與各個(gè)所述相關(guān)子波之間的頻率差異并基于所述頻率差異獲取所述相關(guān)子波重構(gòu)到所述目標(biāo)子波上的子波重構(gòu)因子，以基于所述子波重構(gòu)因子對(duì)對(duì)應(yīng)區(qū)塊的所述能量校正后的地震偏移數(shù)據(jù)進(jìn)行頻率校正的步驟包括：基于互相關(guān)算法以及所述目標(biāo)頻帶范圍構(gòu)建所述若干區(qū)塊分別對(duì)應(yīng)的相關(guān)子波重構(gòu)到所述目標(biāo)子波的子波重構(gòu)因子；基于所述子波重構(gòu)因子對(duì)所述能量校正后的地震偏移數(shù)據(jù)進(jìn)行頻率校正；判斷頻率校正后的地震偏移數(shù)據(jù)是否滿足所述目的層的構(gòu)造解釋成果；響應(yīng)于不滿足所述目的層的構(gòu)造解釋成果，返回到重構(gòu)所述子波重構(gòu)因子的步驟，直至頻率校正后的地震偏移數(shù)據(jù)滿足所述目的層的構(gòu)造解釋成果。

7、在一些實(shí)施例中，所述對(duì)所述若干區(qū)塊之間的重合區(qū)域的頻率校正后的地震偏移數(shù)據(jù)進(jìn)行時(shí)差校正后與所述若干區(qū)塊之間的未重合區(qū)域的頻率校正后的地震偏移數(shù)據(jù)進(jìn)行融合，得到對(duì)所述地震偏移數(shù)據(jù)的融合成果的步驟包括：獲取所述若干區(qū)塊之間的重合區(qū)域；基于互相關(guān)算法計(jì)算所述重合區(qū)域的時(shí)差并進(jìn)行平滑處理；基于平滑處理后的時(shí)差對(duì)頻率校正后的地震偏移數(shù)據(jù)進(jìn)行時(shí)移得到所述若干區(qū)塊之間進(jìn)行拼接的最優(yōu)位置，以基于所述最優(yōu)位置得到所述重合區(qū)域的頻率校正后的地震偏移數(shù)據(jù)在時(shí)間域上的時(shí)差校正因子；基于所述時(shí)差校正因子對(duì)所述重合區(qū)域的頻率校正后的地震偏移數(shù)據(jù)進(jìn)行時(shí)差校正以及與所述若干區(qū)塊之間的未重合區(qū)域的頻率校正后的地震偏移數(shù)據(jù)進(jìn)行融合，得到對(duì)所述地震偏移數(shù)據(jù)的融合成果。

8、在一些實(shí)施例中，所述基于所述最優(yōu)位置得到所述重合區(qū)域的頻率校正后的地震偏移數(shù)據(jù)在時(shí)間域上的時(shí)差校正因子的步驟包括：基于所述最優(yōu)位置獲取所述重合區(qū)域?qū)?yīng)的時(shí)差校正應(yīng)用程度因子；將所述時(shí)差校正應(yīng)用程度因子與所述重合區(qū)域?qū)?yīng)的平滑處理后的時(shí)差相乘，得到對(duì)應(yīng)的時(shí)差校正因子。

9、本發(fā)明實(shí)施例的另一方面，還提供了一種地震偏移數(shù)據(jù)融合的裝置，包括：第一模塊，用于獲取若干區(qū)塊分別對(duì)應(yīng)的地震偏移數(shù)據(jù)并基于目的層特征和地質(zhì)任務(wù)選定所述若干區(qū)塊中的一個(gè)區(qū)塊為基準(zhǔn)區(qū)塊；第二模塊，用于基于所述基準(zhǔn)區(qū)塊的絕對(duì)值振幅與所述若干區(qū)塊中的非基準(zhǔn)區(qū)塊分別對(duì)應(yīng)的絕對(duì)值振幅之間的能量差異構(gòu)建振幅校正序列，以基于所述振幅校正序列對(duì)所述地震偏移數(shù)據(jù)進(jìn)行能量校正；第三模塊，用于基于能量校正后的地震偏移數(shù)據(jù)提取所述若干區(qū)塊分別在所述目的層的相關(guān)子波，以基于所述相關(guān)子波和所述地質(zhì)任務(wù)構(gòu)建目標(biāo)子波；第四模塊，用于計(jì)算所述目標(biāo)子波與各個(gè)所述相關(guān)子波之間的頻率差異并基于所述頻率差異獲取所述相關(guān)子波重構(gòu)到所述目標(biāo)子波上的子波重構(gòu)因子，以基于所述子波重構(gòu)因子對(duì)對(duì)應(yīng)區(qū)塊的所述能量校正后的地震偏移數(shù)據(jù)進(jìn)行頻率校正；第五模塊，用于對(duì)所述若干區(qū)塊之間的重合區(qū)域的頻率校正后的地震偏移數(shù)據(jù)進(jìn)行時(shí)差校正后與所述若干區(qū)塊之間的未重合區(qū)域的頻率校正后的地震偏移數(shù)據(jù)進(jìn)行融合，得到對(duì)所述地震偏移數(shù)據(jù)的融合成果。

10、本發(fā)明實(shí)施例的另一方面，還提供了一種電子設(shè)備，包括至少一個(gè)處理器；以及存儲(chǔ)器，存儲(chǔ)器存儲(chǔ)有可在處理器上運(yùn)行的計(jì)算機(jī)指令，指令由處理器執(zhí)行時(shí)實(shí)現(xiàn)上述方法的步驟。

11、本發(fā)明實(shí)施例的另一方面，還提供了一種計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)，計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)存儲(chǔ)有被處理器執(zhí)行時(shí)實(shí)現(xiàn)如上述方法步驟的計(jì)算機(jī)程序。

12、本發(fā)明至少具有以下有益效果：本發(fā)明提出一種地震偏移數(shù)據(jù)融合的方法，解決了三維區(qū)塊之間的能量、頻率差異和邊界成像差異，改善了疊后拼接處理成果的能量、頻率和邊界成像差異，可以得到一致性好的連片偏移成果，沒有明顯拼接痕跡。同時(shí)，本發(fā)明提出的一種地震偏移數(shù)據(jù)融合的方法可以對(duì)不同輪的處理數(shù)據(jù)之間的能量、頻率、時(shí)差進(jìn)行校正，可操作性強(qiáng)、運(yùn)行效率高。