[發(fā)明專利]一種深度域反偏移方法及裝置

有效

• 申請?zhí)枺?/td>: 201810146031.6
• 申請日: 2018-02-12
• 公開（公告）號: CN108445532B
• 公開（公告）日: 2019-11-08
• 發(fā)明（設(shè)計）人: 章曉;戴南潯;武威;王東;李文卉;吳曉豐
• 申請（專利權(quán)）人: 中國石油天然氣集團有限公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京三友知識產(chǎn)權(quán)代理有限公司 11127
• 代理人: 李輝;劉飛

1.一種深度域反偏移方法，其特征在于，包括：

獲取工區(qū)的地震觀測數(shù)據(jù)；

構(gòu)建所述工區(qū)的模擬地震波場，并根據(jù)所述地震觀測數(shù)據(jù)獲取所述模擬地震波場的深度偏移成像剖面；所述模擬地震波場包括下行波場及上行波場；

根據(jù)所述深度偏移成像剖面獲取底層波場的頻率波數(shù)域波場；

從底層波場開始，根據(jù)所述工區(qū)內(nèi)每一層波場的頻率波數(shù)域波場及所述上行波場的相移因子，獲取對應的上一層波場的頻率波數(shù)域上行波場，并將所述頻率波數(shù)域上行波場變換成對應的頻率空間域上行波場；

將每一層波場的頻率空間域上行波場與其頻率空間域反射波場疊加，獲得每一層波場的頻率空間域地震波場；

將頂層波場的頻率空間域地震波場進行頻率域傅里葉逆變換，獲得時間空間域地震數(shù)據(jù)；

確定所述地震觀測數(shù)據(jù)與所述時間空間域地震數(shù)據(jù)的殘差，并根據(jù)所述殘差優(yōu)化所述深度偏移成像剖面，直至所述殘差滿足預設(shè)的目標函數(shù)最小化要求。

2.如權(quán)利要求1所述的深度域反偏移方法，其特征在于，所述下行波場包括：

且所述下行波場的相移因子為exp(-ik_zΔz)；

所述上行波場包括且所述上行波場的相移因子為exp(ik_zΔz)；

其中，x為地下介質(zhì)中質(zhì)點的平面位置，x_s為震源的平面位置，z為地下介質(zhì)中質(zhì)點的深度位置，ω為地震頻率，p_z為地震波在傳播過程中地下介質(zhì)中每個質(zhì)點的位移，P_D(x,z,ω)為下行波場，P_D(x,z＝0,ω)為零深度時的下行波場，δ(x-x_s)為狄拉克函數(shù)；P_U(x,z,ω)為上行波場，P_U(x,z＝0,ω)為零深度時的上行波場，為地震觀測數(shù)據(jù)，Δz為網(wǎng)格步長。

3.如權(quán)利要求2所述的深度域反偏移方法，其特征在于，所述根據(jù)所述地震觀測數(shù)據(jù)獲取所述模擬地震波場的深度偏移成像剖面，包括：

將所述地震觀測數(shù)據(jù)、所述下行波場的相移因子及預設(shè)的速度模型作為輸入并求解所述模擬地震波場；

根據(jù)預設(shè)的互相關(guān)成像條件I(x,z_j)＝P^S(x,z_j,ω_l)P^R(x,z_j,ω_l)j＝1,2,…,N，對求解后的模擬地震波場進行疊前深度偏移，獲得所述模擬地震波場的深度偏移成像剖面；

其中，I(x,z_j)為介質(zhì)反射系數(shù)，P^S(x,z_j,ω_l)為下行波場，P^R(x,z_j,ω_l)為上行波場，z_j為z方向第j個網(wǎng)格點的坐標，ω_l為頻率。

4.如權(quán)利要求2所述的深度域反偏移方法，其特征在于，所述根據(jù)所述深度偏移成像剖面獲取底層波場的頻率波數(shù)域波場，包括：

根據(jù)公式獲取所述工區(qū)的頻率空間域反射波場；

對所述頻率空間域反射波場沿x方向進行傅里葉變換，獲得底層波場的頻率波數(shù)域波場；

其中，為所述工區(qū)的頻率空間域反射波場，I(x,z_j)為介質(zhì)反射系數(shù)，P^S(x,z_j,ω_l)為下行波場，z_j為z方向第j個網(wǎng)格點的坐標，ω_l為頻率。