本申請涉及石油地球物理勘探地震數(shù)據(jù)處理技術領域，特別涉及一種疊前深度偏移方法及裝置。

背景技術：

地下地層介質(zhì)可能存在非完全彈性介質(zhì)，其對地震波具有吸收衰減作用，因此，在疊前深度偏移過程中考慮這種介質(zhì)的吸收衰減作用，可以提高疊前深度偏移的成像結(jié)果的準確性、精度或分辨率。在疊前深度偏移過程中考慮非完全彈性介質(zhì)對地震波的吸收衰減作用的方法常稱為粘聲介質(zhì)疊前深度偏移方法。

現(xiàn)有技術中常用的粘聲介質(zhì)疊前深度偏移方法的主要過程是：獲取目的工區(qū)的疊前地震數(shù)據(jù)；根據(jù)疊前地震數(shù)據(jù)確定目的工區(qū)的速度場數(shù)據(jù)和品質(zhì)因子場數(shù)據(jù)；根據(jù)品質(zhì)因子場數(shù)據(jù)，計算目的工區(qū)中各個成像點處射線路徑補償旅行時；獲取疊前地震數(shù)據(jù)對應的旅行時表和補償旅行時表，以及疊前地震數(shù)據(jù)中的第一地震道數(shù)據(jù)；根據(jù)旅行時表和補償旅行時表，對第一地震道數(shù)據(jù)進行補償處理；根據(jù)補償處理后的第一地震道數(shù)據(jù)，確定目的工區(qū)中第一成像點的第一成像結(jié)果；對于第一成像點對應的其他地震道數(shù)據(jù)重復上述步驟，以確定目的工區(qū)中該成像點對應的由不同地震道數(shù)據(jù)對應的成像結(jié)果生成的目標成像結(jié)果。

對于一個成像點，其可以對應多個炮點或檢波點，那么不同炮點或檢波點對應的旅行時不相同，同時，不同炮點或檢波點對應的補償旅行時間也不同。對于疊前地震數(shù)據(jù)中的某一地震道數(shù)據(jù)，該地震道數(shù)據(jù)中不同采樣時間對應的相位和振幅補償因子也不同，因此，現(xiàn)有技術中常用的粘聲介質(zhì)疊前深度偏移方法需要先將第一地震道數(shù)據(jù)從時間域變換到頻率域，然后計算第一地震道數(shù)據(jù)對應的所有頻點的補償結(jié)果，再將這些補償結(jié)果進行疊加處理，得到第一成像點的第一成像結(jié)果，通過對疊前地震數(shù)據(jù)中所有地震道數(shù)據(jù)循環(huán)進行上述補償過程，得到第一成像點的目標成像結(jié)果。

發(fā)明人發(fā)現(xiàn)現(xiàn)有技術中至少存在如下問題：由于需要計算疊前地震數(shù)據(jù)中每一個地震道數(shù)據(jù)對應的各個頻點的補償結(jié)果，以及將這些補償結(jié)果進行疊加處理，導致現(xiàn)有技術中常用的粘聲介質(zhì)疊前深度偏移方法的計算量較大，從而導致采用該方法進行疊前深度偏移的運算效率較低。

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

本申請實施例的目的是提供一種疊前深度偏移方法及裝置，以提高疊前深度偏移的運算效率。

為解決上述技術問題，本申請實施例提供一種疊前深度偏移方法及裝置是這樣實現(xiàn)的：

一種疊前深度偏移方法，包括：

獲取目的工區(qū)的時間域的疊前地震數(shù)據(jù)，基于所述疊前地震數(shù)據(jù)，分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息；

根據(jù)所述補償旅行時信息和預設間隔時間，對所述疊前地震數(shù)據(jù)中第一地震道數(shù)據(jù)進行補償預處理，得到第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)；

根據(jù)所述旅行時信息和補償旅行時信息，分別確定第一地震道數(shù)據(jù)對應的炮點位置到所述目的工區(qū)中第一成像點位置的總旅行時和總補償旅行時，以及分別確定第一地震道數(shù)據(jù)對應的檢波點位置到所述目的工區(qū)中第一成像點位置的總旅行時和總補償旅行時；

基于所述總旅行時和所述總補償旅行時，對所述初始補償?shù)卣鸬罃?shù)據(jù)進行插值處理，得到第一地震道數(shù)據(jù)對應的目標補償?shù)卣鸬罃?shù)據(jù)。

優(yōu)選方案中，所述基于疊前地震數(shù)據(jù)，分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息，包括：

基于所述疊前地震數(shù)據(jù)，確定所述目的工區(qū)的速度場數(shù)據(jù)和品質(zhì)因子場數(shù)據(jù)；

根據(jù)所述速度場數(shù)據(jù)和所述品質(zhì)因子場數(shù)據(jù)，分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息。

優(yōu)選方案中，所述根據(jù)速度場數(shù)據(jù)和品質(zhì)因子場數(shù)據(jù)，分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息，包括：

根據(jù)所述速度場數(shù)據(jù)，分別計算所述目的工區(qū)中各個炮點位置到所述目的工區(qū)中各個成像點位置的旅行時，以及分別計算所述目的工區(qū)中各個檢波點位置到所述目的工區(qū)中各個成像點位置的旅行時；

根據(jù)所述品質(zhì)因子場數(shù)據(jù)，分別計算所述目的工區(qū)中各個炮點位置到所述目的工區(qū)中各個成像點位置的補償旅行時，以及分別計算所述目的工區(qū)中各個檢波點位置到所述目的工區(qū)中各個成像點位置的補償旅行時；

基于所述目的工區(qū)中各個炮點位置到所述目的工區(qū)中各個成像點位置的旅行時和補償旅行時，以及所述目的工區(qū)中各個檢波點位置到所述目的工區(qū)中各個成像點位置的旅行時和補償旅行時，分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息。

優(yōu)選方案中，采用下述公式分別計算所述目的工區(qū)中各個炮點位置到所述目的工區(qū)中各個成像點位置的補償旅行時，以及分別計算所述目的工區(qū)中各個檢波點位置到所述目的工區(qū)中各個成像點位置的補償旅行時：

公式中，T^*(x)表示所述目的工區(qū)中一個炮點位置到所述目的工區(qū)中一個成像點位置x的補償旅行時，或所述目的工區(qū)中一個檢波點位置到所述目的工區(qū)中一個成像點位置x的補償旅行時，v₀表示地震波在所述目的工區(qū)中的傳播速度，Q表示所述品質(zhì)因子場中一個成像點位置處的品質(zhì)因子，x表示所述目的工區(qū)中一個成像點位置，s表示所述目的工區(qū)中一個炮點位置分別到所述目的工區(qū)中一個成像點位置的射線路徑，或所述目的工區(qū)中一個檢波點位置分別到所述目的工區(qū)中一個成像點位置的射線路徑。

優(yōu)選方案中，所述根據(jù)補償旅行時信息和預設間隔時間，對所述疊前地震數(shù)據(jù)中第一地震道數(shù)據(jù)進行補償預處理，得到第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)，包括：

通過傅立葉變換將第一地震道數(shù)據(jù)從時間域轉(zhuǎn)換到頻率域，得到頻率域的第一地震道數(shù)據(jù)；

獲取所述補償旅行時信息中的最大補償旅行時；

利用所述預設間隔時間將所述最大補償旅行時劃分為等間隔的補償旅行時；

基于所述頻率域的第一地震道數(shù)據(jù)和所述等間隔的補償旅行時，確定第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)。

優(yōu)選方案中，采用下述公式確定第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)：

公式中，表示第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)，表示預設補償濾波函數(shù)，i表示虛數(shù)單位，ω表示第一地震道數(shù)據(jù)的角頻率，u(ω)表示所述頻率域的第一地震道數(shù)據(jù)，為所述等間隔的補償旅行時，t表示所述時間域中第一地震道數(shù)據(jù)的采樣時間。

優(yōu)選方案中，所述根據(jù)旅行時信息和補償旅行時信息，分別確定第一地震道數(shù)據(jù)對應的炮點位置到所述目的工區(qū)中第一成像點位置的總旅行時和總補償旅行時，以及分別確定第一地震道數(shù)據(jù)對應的檢波點位置到所述目的工區(qū)中第一成像點位置的總旅行時和總補償旅行時，包括：

根據(jù)所述旅行時信息和補償旅行時信息，分別獲取第一地震道數(shù)據(jù)對應的炮點位置和檢波點位置到所述目的工區(qū)中第一成像點位置的旅行時，以及分別獲取第一地震道數(shù)據(jù)對應的炮點位置和檢波點位置到所述目的工區(qū)中第一成像點位置的補償旅行時；

將所述第一地震道數(shù)據(jù)對應的炮點位置到第一成像點位置的旅行時與第一地震道數(shù)據(jù)對應的檢波點位置到第一成像點位置的旅行時相加，得到第一地震道數(shù)據(jù)對應的炮點位置和檢波點位置到所述目的工區(qū)中第一成像點位置的總旅行時；

將所述第一地震道數(shù)據(jù)對應的炮點位置到第一成像點位置的補償旅行時與第一地震道數(shù)據(jù)對應的檢波點位置到第一成像點位置的補償旅行時相加，得到第一地震道數(shù)據(jù)對應的炮點位置和檢波點位置到所述目的工區(qū)中第一成像點位置的總補償旅行時。

優(yōu)選方案中，采用下述公式進行所述插值處理，得到第一地震道數(shù)據(jù)對應的目標補償?shù)卣鸬罃?shù)據(jù)：

公式中，T和T^*分別表示所述總旅行時和所述總補償旅行時，和分別表示與所述總補償旅行時T^*緊鄰的前后兩個等間隔的補償旅行時，表示第一地震道數(shù)據(jù)對應的目標補償?shù)卣鸬罃?shù)據(jù)，表示當旅行時為T且等間隔的補償旅行時為時所對應的所述初始補償?shù)卣鸬罃?shù)據(jù)，表示當旅行時為T且等間隔的補償旅行時為時所對應的所述初始補償?shù)卣鸬罃?shù)據(jù)。

一種疊前深度偏移裝置，所述裝置包括：旅行時信息確定模塊、初始補償數(shù)據(jù)確定模塊、總旅行時確定模塊和目標補償?shù)卣鸬罃?shù)據(jù)確定模塊；其中，

所述旅行時信息確定模塊，用于獲取目的工區(qū)的時間域的疊前地震數(shù)據(jù)，基于所述疊前地震數(shù)據(jù)，分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息；

所述初始補償數(shù)據(jù)確定模塊，用于根據(jù)所述補償旅行時信息和預設間隔時間，對所述疊前地震數(shù)據(jù)中第一地震道數(shù)據(jù)進行補償預處理，得到第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)；

所述總旅行時確定模塊，用于根據(jù)所述旅行時信息和補償旅行時信息，分別確定第一地震道數(shù)據(jù)對應的炮點位置到所述目的工區(qū)中第一成像點位置的總旅行時和總補償旅行時，以及分別確定第一地震道數(shù)據(jù)對應的檢波點位置到所述目的工區(qū)中第一成像點位置的總旅行時和總補償旅行時；

所述目標補償?shù)卣鸬罃?shù)據(jù)確定模塊，用于基于所述總旅行時和所述總補償旅行時，對所述初始補償?shù)卣鸬罃?shù)據(jù)進行插值處理，得到第一地震道數(shù)據(jù)對應的目標補償?shù)卣鸬罃?shù)據(jù)。

優(yōu)選方案中，所述旅行時信息確定模塊，包括：場數(shù)據(jù)確定模塊和成像點旅行時信息確定模塊；其中，

所述場數(shù)據(jù)確定模塊，用于基于所述疊前地震數(shù)據(jù)，確定所述目的工區(qū)的速度場數(shù)據(jù)和品質(zhì)因子場數(shù)據(jù)；

所述成像點旅行時信息確定模塊，用于根據(jù)所述速度場數(shù)據(jù)和所述品質(zhì)因子場數(shù)據(jù)，分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息。

優(yōu)選方案中，所述初始補償數(shù)據(jù)確定模塊，包括：頻率域數(shù)據(jù)確定模塊、最大補償旅行時獲取模塊、等間隔補償旅行時確定模塊和初始補償?shù)卣鸬罃?shù)據(jù)確定模塊；其中，

所述頻率域數(shù)據(jù)確定模塊，用于通過傅立葉變換將第一地震道數(shù)據(jù)從時間域轉(zhuǎn)換到頻率域，得到頻率域的第一地震道數(shù)據(jù)；

所述最大補償旅行時獲取模塊，用于獲取所述補償旅行時信息中的最大補償旅行時；

所述等間隔補償旅行時確定模塊，用于利用所述預設間隔時間將所述最大補償旅行時劃分為等間隔的補償旅行時；

所述初始補償?shù)卣鸬罃?shù)據(jù)確定模塊，用于基于所述頻率域的第一地震道數(shù)據(jù)和所述等間隔的補償旅行時，確定第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)。

本申請實施例提供了一種方法及裝置，在疊前深度偏移過程中，對目的工區(qū)的疊前地震數(shù)據(jù)中某個地震道數(shù)據(jù)進行補償預處理，再對補償預處理后的地震道數(shù)據(jù)進行插值處理，得到該地震道數(shù)據(jù)對應的成像點位置處的成像結(jié)果。因此，采用插值處理的方式代替了現(xiàn)有技術中計算每一個地震道數(shù)據(jù)對應的各個頻點的補償結(jié)果以及將這些補償結(jié)果進行疊加處理的處理過程，可以減少疊前深度偏移過程中的計算量，進而可以提高疊前深度偏移的運算效率。

附圖說明

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

圖1是本申請一種疊前深度偏移方法實施例的流程圖；

圖2是本申請實施例中目的工區(qū)的速度場數(shù)據(jù)的剖面示意圖；

圖3是本申請實施例中目的工區(qū)的品質(zhì)因子場數(shù)據(jù)的剖面示意圖；

圖4是本申請實施例中針對目的工區(qū)的地震道數(shù)據(jù)分別采用現(xiàn)有技術和本申請方法得到的成像結(jié)果的對比示意圖；

圖5是本申請實施例中采用現(xiàn)有技術得到的疊前深度偏移結(jié)果；

圖6是本申請實施例中采用本申請方法得到的疊前深度偏移結(jié)果；

圖7是本申請實施例中分別采用現(xiàn)有技術和本申請方法得到的疊前深度偏移結(jié)果對應的波數(shù)譜的對比示意圖；

圖8是本申請疊前深度偏移裝置實施例的組成結(jié)構圖；

圖9是本申請疊前深度偏移裝置實施例中旅行時信息確定模塊的組成結(jié)構圖；

圖10是本申請疊前深度偏移裝置實施例中初始補償數(shù)據(jù)確定模塊的組成結(jié)構圖。

具體實施方式

本申請實施例提供一種疊前深度偏移方法及裝置。

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

圖1是本申請一種疊前深度偏移方法實施例的流程圖。如圖1所示，所述疊前深度偏移方法，包括以下步驟。

步驟S101：獲取目的工區(qū)的時間域的疊前地震數(shù)據(jù)，基于所述疊前地震數(shù)據(jù)，分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息。

品質(zhì)因子是用于描述地層巖石的非完全彈性吸收對地震波傳播能量衰減特性強弱的物理量之一，通常為當?shù)卣鸩▊鞑ヒ粋€波長或一個周期距離時，地震波儲能與耗散能的比率。

所述補償旅行時信息反映了地震波的振幅衰減和相位變化?？梢垣@取目的工區(qū)的時間域的疊前地震數(shù)據(jù)?；谒霪B前地震數(shù)據(jù)，可以確定所述目的工區(qū)的速度場數(shù)據(jù)和品質(zhì)因子場數(shù)據(jù)。根據(jù)所述速度場數(shù)據(jù)和所述品質(zhì)因子場數(shù)據(jù)，可以分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息。

對所述目的工區(qū)進行地震勘探和數(shù)據(jù)采集，可以獲取所述目的工區(qū)的時間域的疊前地震數(shù)據(jù)。

所述基于疊前地震數(shù)據(jù)，確定所述目的工區(qū)的速度場數(shù)據(jù)，可以包括：基于所述疊前地震數(shù)據(jù)，可以得到目的工區(qū)的初始均方根速度場數(shù)據(jù)；基于所述初始均方根速度場數(shù)據(jù)，對所述疊前地震數(shù)據(jù)進行層析反演處理，可以得到所述目的工區(qū)的速度場數(shù)據(jù)。

在一種實施方式中，所述基于疊前地震數(shù)據(jù)，得到目的工區(qū)的初始均方根速度場數(shù)據(jù)，可以包括：基于所述疊前地震數(shù)據(jù)，可以確定所述目的工區(qū)的疊加速度譜；對疊加速度譜進行插值處理和平滑處理，可以得到目的工區(qū)的初始均方根速度場數(shù)據(jù)。

所述基于疊前地震數(shù)據(jù)，確定所述目的工區(qū)的品質(zhì)因子場數(shù)據(jù)，可以包括：可以獲取所述目的工區(qū)的垂直地震剖面數(shù)據(jù)；基于所述垂直地震剖面數(shù)據(jù)和所述疊前地震數(shù)據(jù)，可以提取所述目的工區(qū)的等效品質(zhì)因子場數(shù)據(jù)；對所述等效品質(zhì)因子數(shù)據(jù)進行時深轉(zhuǎn)換處理和層析反演處理，可以得到所述目的工區(qū)的品質(zhì)因子場數(shù)據(jù)。

例如，圖2是本申請實施例中目的工區(qū)的速度場數(shù)據(jù)的剖面示意圖。圖3是本申請實施例中目的工區(qū)的品質(zhì)因子場數(shù)據(jù)的剖面示意圖。圖2和圖3中CMP表示所述疊前地震數(shù)據(jù)中的共中心點，time表示所述疊前地震數(shù)據(jù)中的采樣時間。

所述根據(jù)速度場數(shù)據(jù)和所述品質(zhì)因子場數(shù)據(jù)，分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息，可以包括：根據(jù)所述速度場數(shù)據(jù)，通過射線追蹤方法，可以分別計算所述目的工區(qū)中各個炮點位置到所述目的工區(qū)中各個成像點位置的旅行時，以及分別計算所述目的工區(qū)中各個檢波點位置到所述目的工區(qū)中各個成像點位置的旅行時；根據(jù)所述品質(zhì)因子場數(shù)據(jù)，可以分別計算所述目的工區(qū)中各個炮點位置到所述目的工區(qū)中各個成像點位置的補償旅行時，以及分別計算所述目的工區(qū)中各個檢波點位置到所述目的工區(qū)中各個成像點位置的補償旅行時；基于所述目的工區(qū)中各個炮點位置到所述目的工區(qū)中各個成像點位置的旅行時和補償旅行時，以及所述目的工區(qū)中各個檢波點位置到所述目的工區(qū)中各個成像點位置的旅行時和補償旅行時，可以分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息。所述旅行時信息表示所述目的工區(qū)中各個炮點位置和各個檢波點位置到所述目的工區(qū)中各個成像點位置的旅行時。所述補償旅行時信息表示所述目的工區(qū)中各個炮點位置和各個檢波點位置到所述目的工區(qū)中各個成像點位置的補償旅行時。

進一步地，可以采用下述公式分別計算所述目的工區(qū)中各個炮點位置到所述目的工區(qū)中各個成像點位置的補償旅行時，以及分別計算所述目的工區(qū)中各個檢波點位置到所述目的工區(qū)中各個成像點位置的補償旅行時：

公式(1)中，T^*(x)表示所述目的工區(qū)中一個炮點位置到所述目的工區(qū)中一個成像點位置x的補償旅行時，或所述目的工區(qū)中一個檢波點位置到所述目的工區(qū)中各個成像點位置x的補償旅行時，v₀表示地震波在所述目的工區(qū)中的傳播速度，Q表示所述品質(zhì)因子場中一個成像點位置處的品質(zhì)因子，x表示所述目的工區(qū)中一個成像點位置，s表示所述目的工區(qū)中一個炮點位置分別到所述目的工區(qū)中一個成像點位置的射線路徑，或所述目的工區(qū)中一個檢波點位置分別到所述目的工區(qū)中一個成像點位置的射線路徑。

步驟S102：根據(jù)所述補償旅行時信息和預設間隔時間，對所述疊前地震數(shù)據(jù)中第一地震道數(shù)據(jù)進行補償預處理，得到第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)。

具體地，可以通過傅立葉變換將第一地震道數(shù)據(jù)從時間域轉(zhuǎn)換到頻率域，得到頻率域的第一地震道數(shù)據(jù)。可以獲取所述補償旅行時信息中的最大補償旅行時。可以利用所述預設間隔時間將所述最大補償旅行時劃分為等間隔的補償旅行時?；谒鲱l率域的第一地震道數(shù)據(jù)和所述等間隔的補償旅行時，可以確定第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)。所述第一地震道數(shù)據(jù)可以表示所述疊前地震數(shù)據(jù)中任意一個地震道。

進一步地，可以采用下述公式確定第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)：

公式(2)中，σ＝exp[-(0.23G+1.63)]。

其中，表示第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)，表示預設補償濾波函數(shù)，表示預設振幅補償因子，i表示虛數(shù)單位，ω表示第一地震道數(shù)據(jù)的角頻率，ω₀表示第一地震道數(shù)據(jù)的主頻，u(ω)表示所述頻率域的第一地震道數(shù)據(jù)，為所述等間隔的補償旅行時，t表示所述時間域中第一地震道數(shù)據(jù)的采樣時間，G表示預設增益限制量，用于控制高頻成分振幅補償?shù)姆€(wěn)定性。

步驟S103：根據(jù)所述旅行時信息和補償旅行時信息，分別確定第一地震道數(shù)據(jù)對應的炮點位置到所述目的工區(qū)中第一成像點位置的總旅行時和總補償旅行時，以及分別確定第一地震道數(shù)據(jù)對應的檢波點位置到所述目的工區(qū)中第一成像點位置的總旅行時和總補償旅行時。

具體地，根據(jù)所述旅行時信息和補償旅行時信息，可以分別獲取第一地震道數(shù)據(jù)對應的炮點位置和檢波點位置到所述目的工區(qū)中第一成像點位置的旅行時，以及分別獲取第一地震道數(shù)據(jù)對應的炮點位置和檢波點位置到所述目的工區(qū)中第一成像點位置的補償旅行時。將所述第一地震道數(shù)據(jù)對應的炮點位置到第一成像點位置的旅行時與第一地震道數(shù)據(jù)對應的檢波點位置到第一成像點位置的旅行時相加，可以得到第一地震道數(shù)據(jù)對應的炮點位置和檢波點位置到所述目的工區(qū)中第一成像點位置的總旅行時。將所述第一地震道數(shù)據(jù)對應的炮點位置到第一成像點位置的補償旅行時與第一地震道數(shù)據(jù)對應的檢波點位置到第一成像點位置的補償旅行時相加，可以得到第一地震道數(shù)據(jù)對應的炮點位置和檢波點位置到所述目的工區(qū)中第一成像點位置的總補償旅行時。

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

步驟S104：基于所述初始補償?shù)卣鸬罃?shù)據(jù)、所述總旅行時和所述總補償旅行時，確定第一地震道數(shù)據(jù)對應的第一成像點位置處的成像結(jié)果。

具體地，基于所述總旅行時和所述總補償旅行時，可以對所述初始補償?shù)卣鸬罃?shù)據(jù)進行插值處理，得到第一地震道數(shù)據(jù)對應的目標補償?shù)卣鸬罃?shù)據(jù)?？梢詫⑺瞿繕搜a償?shù)卣鸬罃?shù)據(jù)作為第一地震道數(shù)據(jù)對應的第一成像點位置處的成像結(jié)果。

進一步地，可以采用下述公式進行所述插值處理，得到第一地震道數(shù)據(jù)對應的目標補償?shù)卣鸬罃?shù)據(jù)：

公式(3)中，T和T^*分別表示所述總旅行時和所述總補償旅行時，和分別表示與所述總補償旅行時T^*緊鄰的前后兩個等間隔的補償旅行時，u'(T,T^*)表示第一地震道數(shù)據(jù)對應的目標補償?shù)卣鸬罃?shù)據(jù)，表示當旅行時為T且等間隔的補償旅行時為時所對應的所述初始補償?shù)卣鸬罃?shù)據(jù)，表示當旅行時為T且等間隔的補償旅行時為時所對應的所述初始補償?shù)卣鸬罃?shù)據(jù)。

例如，圖4是本申請實施例中針對目的工區(qū)的地震道數(shù)據(jù)分別采用現(xiàn)有技術和本申請方法得到的補償結(jié)果的對比示意圖?？v坐標表示補償后地震數(shù)據(jù)的振幅能量，橫坐標表示地震道數(shù)據(jù)的采樣時間，星點表示現(xiàn)有技術的補償結(jié)果，實線表示本申請方法得到的補償結(jié)果。如圖4所示，本申請方法得到的補償結(jié)果保持了現(xiàn)有技術補償結(jié)果中的波形特征，波形和振幅能量都非常接近，采用現(xiàn)有技術和本申請方法得到的補償結(jié)果一致性較好。

所述疊前地震數(shù)據(jù)可以包括多個地震道數(shù)據(jù)。所述目的工區(qū)可以包括多個成像點。按照確定第一地震道數(shù)據(jù)對應的第一成像點位置處的成像結(jié)果的步驟進行循環(huán)處理，可以分別確定所述疊前地震數(shù)據(jù)中每一個地震道數(shù)據(jù)對應的所述目的工區(qū)中第一成像點位置處的成像結(jié)果；基于所述每一個地震道數(shù)據(jù)分別對應的所述目的工區(qū)中第一成像點位置處的成像結(jié)果，可以確定所述目的工區(qū)中第一成像點位置處的目標成像結(jié)果。

具體地，所述基于每一個地震道數(shù)據(jù)分別對應的所述目的工區(qū)中第一成像點位置處的成像結(jié)果，確定所述目的工區(qū)中第一成像點位置處的目標成像結(jié)果，可以包括：可以將所述每一個地震道數(shù)據(jù)分別對應的所述目的工區(qū)中第一成像點位置處的成像結(jié)果進行加權系數(shù)疊加處理，確定所述目的工區(qū)中第一成像點位置處的目標成像結(jié)果。

進一步地，所述將每一個地震道數(shù)據(jù)分別對應的所述目的工區(qū)中第一成像點位置處的成像結(jié)果進行加權系數(shù)疊加處理，確定所述目的工區(qū)中第一成像點位置處的目標成像結(jié)果，可以包括：可以獲取所述每一個地震道數(shù)據(jù)分別對應的炮檢距；根據(jù)所述每一個地震道數(shù)據(jù)分別對應的炮檢距，可以確定所述每一個地震道數(shù)據(jù)分別對應的預設加權系數(shù)；可以將所述每一個地震道數(shù)據(jù)對應的第一成像點位置處的成像結(jié)果分別與所述每一個地震道數(shù)據(jù)對應的預設加權系數(shù)相乘，分別得到所述每一個地震道數(shù)據(jù)對應的第一成像點位置處的加權成像結(jié)果；將所述每一個地震道數(shù)據(jù)對應的第一成像點位置處的加權成像結(jié)果相加，可以確定所述目的工區(qū)中第一成像點位置處的目標成像結(jié)果。

按照確定第一成像點位置處的目標成像結(jié)果的步驟進行循環(huán)處理，可以確定所述目的工區(qū)中每一個成像點位置處的目標成像結(jié)果，進而可以得到疊前深度偏移結(jié)果。

例如，圖5是本申請實施例中采用現(xiàn)有技術得到的疊前深度偏移結(jié)果。圖6是本申請實施例中采用本申請方法得到的疊前深度偏移結(jié)果。圖5和圖6中CMP表示所述疊前地震數(shù)據(jù)中的共中心點，time表示所述疊前地震數(shù)據(jù)中的采樣時間。圖7是本申請實施例中分別采用現(xiàn)有技術和本申請方法得到的疊前深度偏移結(jié)果對應的波數(shù)譜的對比示意圖。圖7中的橫坐標和縱坐標分別表示波數(shù)譜的頻率和振幅，其中Frequency表示波數(shù)譜的頻率，Amplitude表示波數(shù)譜的振幅。圖7中頻率域結(jié)果表示采用現(xiàn)有技術得到的疊前深度偏移結(jié)果對應的波數(shù)譜，時間域結(jié)果表示采用本申請方法得到的疊前深度偏移結(jié)果對應的波數(shù)譜。如圖5和圖6所示，采用本申請方法得到的疊前深度偏移結(jié)果與采用現(xiàn)有技術得到的疊前深度偏移結(jié)果基本一致。如圖7所示，采用現(xiàn)有技術和本申請方法得到的疊前深度偏移結(jié)果對應的波數(shù)譜比較相近。

表1為分別采用現(xiàn)有技術的方法和本申請方法對同一目的工區(qū)的疊前地震數(shù)據(jù)進行疊前深度偏移時的運行時間、所用的節(jié)點數(shù)和每個節(jié)點使用的CPU個數(shù)。如表1所示，采用現(xiàn)有技術方法進行疊前深度偏移時，使用了20個節(jié)點，且每個節(jié)點使用20個CPU，進行運算，最終完成疊前深度偏移耗時5201秒。而采用本申請方法進行疊前深度偏移時，使用了1個節(jié)點，且每個節(jié)點使用5個CPU，進行運算，最終完成疊前深度偏移耗時210秒。由此可見，相比現(xiàn)有技術的方法，本申請方法的運行效率約提高了倍。

表1

所述疊前深度偏移方法實施例，在疊前深度偏移過程中，對目的工區(qū)的疊前地震數(shù)據(jù)中某個地震道數(shù)據(jù)進行補償預處理，再對補償預處理后的地震道數(shù)據(jù)進行插值處理，得到該地震道數(shù)據(jù)對應的成像點位置處的成像結(jié)果。因此，采用插值處理的方式代替了現(xiàn)有技術中計算每一個地震道數(shù)據(jù)對應的各個頻點的補償結(jié)果以及將這些補償結(jié)果進行疊加處理的處理過程，可以減少疊前深度偏移過程中的計算量，進而可以提高疊前深度偏移的運算效率。

圖8是本申請疊前深度偏移裝置實施例的組成結(jié)構圖。如圖8所示，所述疊前深度偏移轉(zhuǎn)置，可以包括：旅行時信息確定模塊100、初始補償數(shù)據(jù)確定模塊200、總旅行時確定模塊300和目標補償?shù)卣鸬罃?shù)據(jù)確定模塊400。

所述旅行時信息確定模塊100，可以用于獲取目的工區(qū)的時間域的疊前地震數(shù)據(jù)，基于所述疊前地震數(shù)據(jù)，可以分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息。

所述初始補償數(shù)據(jù)確定模塊200，可以用于根據(jù)所述補償旅行時信息和預設間隔時間，對所述疊前地震數(shù)據(jù)中第一地震道數(shù)據(jù)進行補償預處理，得到第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)。

所述總旅行時確定模塊300，可以用于根據(jù)所述旅行時信息和補償旅行時信息，分別確定第一地震道數(shù)據(jù)對應的炮點位置到所述目的工區(qū)中第一成像點位置的總旅行時和總補償旅行時，以及分別確定第一地震道數(shù)據(jù)對應的檢波點位置到所述目的工區(qū)中第一成像點位置的總旅行時和總補償旅行時。

所述目標補償?shù)卣鸬罃?shù)據(jù)確定模塊400，可以用于基于所述總旅行時和所述總補償旅行時，對所述初始補償?shù)卣鸬罃?shù)據(jù)進行插值處理，可以得到第一地震道數(shù)據(jù)對應的目標補償?shù)卣鸬罃?shù)據(jù)。

圖9是本申請疊前深度偏移裝置實施例中旅行時信息確定模塊的組成結(jié)構圖。如圖9所示，圖8中旅行時信息確定模塊100，包括：場數(shù)據(jù)確定模塊110和成像點旅行時信息確定模塊120。

所述場數(shù)據(jù)確定模塊110，可以用于基于所述疊前地震數(shù)據(jù)，確定所述目的工區(qū)的速度場數(shù)據(jù)和品質(zhì)因子場數(shù)據(jù)。

所述成像點旅行時信息確定模塊120，可以用于根據(jù)所述速度場數(shù)據(jù)和所述品質(zhì)因子場數(shù)據(jù)，分別確定所述目的工區(qū)中各個成像點位置的旅行時信息和補償旅行時信息。

圖10是本申請疊前深度偏移裝置實施例中初始補償數(shù)據(jù)確定模塊的組成結(jié)構圖。如圖10所示，圖8中初始補償數(shù)據(jù)確定模塊200，可以包括：頻率域數(shù)據(jù)確定模塊210、最大補償旅行時獲取模塊220、等間隔補償旅行時確定模塊230和初始補償?shù)卣鸬罃?shù)據(jù)確定模塊240。

所述頻率域數(shù)據(jù)確定模塊210，可以用于通過傅立葉變換將第一地震道數(shù)據(jù)從時間域轉(zhuǎn)換到頻率域，得到頻率域的第一地震道數(shù)據(jù)。

所述最大補償旅行時獲取模塊220，可以用于獲取所述補償旅行時信息中的最大補償旅行時。

所述等間隔補償旅行時確定模塊230，可以用于利用所述預設間隔時間將所述最大補償旅行時劃分為等間隔的補償旅行時。

所述初始補償?shù)卣鸬罃?shù)據(jù)確定模塊240，可以用于基于所述頻率域的第一地震道數(shù)據(jù)和所述等間隔的補償旅行時，確定第一地震道數(shù)據(jù)對應的初始補償?shù)卣鸬罃?shù)據(jù)。

所述疊前深度偏移裝置實施例與所述疊前深度偏移方法實施例相對應，可以減少疊前深度偏移過程中的計算量，進而可以提高疊前深度偏移的運行效率。

在20世紀90年代，對于一個技術的改進可以很明顯地區(qū)分是硬件上的改進(例如，對二極管、晶體管、開關等電路結(jié)構的改進)還是軟件上的改進(對于方法流程的改進)。然而，隨著技術的發(fā)展，當今的很多方法流程的改進已經(jīng)可以視為硬件電路結(jié)構的直接改進。設計人員幾乎都通過將改進的方法流程編程到硬件電路中來得到相應的硬件電路結(jié)構。因此，不能說一個方法流程的改進就不能用硬件實體模塊來實現(xiàn)。例如，可編程邏輯器件(Programmable Logic Device,PLD)(例如現(xiàn)場可編程門陣列(Field Programmable Gate Array，F(xiàn)PGA))就是這樣一種集成電路，其邏輯功能由用戶對器件編程來確定。由設計人員自行編程來把一個數(shù)字系統(tǒng)“集成”在一片PLD上，而不需要請芯片制造廠商來設計和制作專用的集成電路芯片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。本領域技術人員也應該清楚，只需要將方法流程用上述幾種硬件描述語言稍作邏輯編程并編程到集成電路中，就可以很容易得到實現(xiàn)該邏輯方法流程的硬件電路。

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

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

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

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

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

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

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

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

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