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

有效

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

本申請提供了一種深度域反偏移方法及裝置，該方法包括：獲取地震觀測數據；構建包括下行波場及上行波場的模擬地震波場并獲取其深度偏移成像剖面；根據深度偏移成像剖面獲取底層波場的頻率波數域波場；根據每一層的頻率波數域波場及上行波場的相移因子，獲取上一層的頻率波數域上行波場并將其變換成對應的頻率空間域上行波場；將每一層的頻率空間域上行波場與頻率空間域反射波場疊加，獲得每一層的頻率空間域地震波場；將頂層波場的頻率空間域地震波場作頻率域傅里葉逆變換，獲得時間空間域地震數據；根據地震觀測數據與時間空間域地震數據的殘差優化深度偏移成像剖面，直至殘差滿足預設的目標函數最小化要求。本申請可提高深度域反偏移的精度。

技術領域

本申請涉及勘探地震數據處理技術領域，尤其是涉及一種深度域反偏移方法及裝置。

背景技術

地震數據處理的重要目的之一是通過深度偏移(Depth Migration)得到高質量的深度域圖像。深度域圖像的質量依賴于原始地震單炮數據的質量，其對于提高地震數據的信噪比及在整個地震數據處理中起著非常重要的作用。深度域反偏移(Depth DomainDemigration)是一種數字波場成像方法，在反偏移過程中如果速度模型不變，則反偏移就是偏移的逆運算。

目前最小二乘偏移技術作為一種深度域反偏移方法，其依賴于精確的速度模型，迭代過程中正演模擬波場和實際觀測波場存在周波跳躍現象，以致于迭代過程無法收斂到可靠的結果，因此，目前現有的深度域反偏移方法的精度亟待提高。

發明內容

本申請實施例的目的在于提供一種深度域反偏移方法及裝置，以提高深度域反偏移的精度。

為達到上述目的，一方面，本申請實施例提供了一種深度域反偏移方法，包括：

獲取工區的地震觀測數據；

構建所述工區的模擬地震波場，并根據所述地震觀測數據獲取所述模擬地震波場的深度偏移成像剖面；所述模擬地震波場包括下行波場及上行波場；

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

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

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

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

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

較佳的，所述下行波場包括：

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

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

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

較佳的，所述根據所述地震觀測數據獲取所述模擬地震波場的深度偏移成像剖面，包括：