[發明專利]基于全波波場VSP/RVSP地震資料的逆時偏移方法

有效

• 申請號: 201711181298.0
• 申請日: 2017-11-23
• 公開（公告）號: CN108037526B
• 公開（公告）日: 2018-12-07
• 發明（設計）人: 孫小東;賈延睿;宋煜;葛中慧
• 申請（專利權）人: 中國石油大學（華東）
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京志霖恒遠知識產權代理事務所(普通合伙) 11435
• 代理人: 劉子成

1.一種基于全波波場VSP/RVSP地震資料的逆時偏移方法，其特征在于通過以下過程來實現：

(1)VSP/RVSP地震資料的觀測系統變換；

(2)有限差分算法求解波動方程，具體如下：

1)時間-空間域的波動方程，經過差分后，利用有限差分法進行數值求解，其中，任意階差分精度的全波波動方程如下：

2)對上式進行傅里葉變換，得到頻空域的形式：

3)確定波場外推計算的時間步長，即需要滿足的穩定性條件：

Dt<h/(2.0*Vmax)

其中，h表示網格剖分后最大的網格間距，Vmax表示速度模型中最大速度值，Dt表示波場外推計算的最大時間步長；

(3)依據速度場分布情況對VSP/RVSP逆時偏移的計算區域進行變網格剖分；

(4)利用高階有限差分算法對網格化的計算區域進行波動方程求解，實現對震源波場的正向延拓和檢波波場的反向延拓；

(5)對延拓后的波場應用成像條件，該成像條件是考慮照明度的影響，根據下式

對震源波場求取其能量分布，以此對成像的相干值進行規則化；

(6)在所有數據處理完后，進行最終的偏移疊加，并利用拉普拉斯濾波或帶通濾波方法做疊后去噪，得出最終的VSP/RVSP成像剖面。

2.根據權利要求1所述的基于全波波場VSP/RVSP地震資料的逆時偏移方法，其特征在于所述的過程(1)具體為：在逆時偏移前需要對數據進行道頭編輯；基于射線路徑可逆原理，進行炮檢置換和道集重排，得到適用于逆時偏移的共炮點道集。

3.根據權利要求1所述的基于全波波場VSP/RVSP地震資料的逆時偏移方法，其特征在于所述的過程(3)中根據地下地質構造的形態，偏移速度場的分布情況，合理進行計算區域的網格化，保證網格尺寸滿足以下公式，從而波場計算不會產生頻散：

h<Vmin/(10.0*Fmax)

其中，h表示網格剖分后最大的網格間距；Vmin表示速度模型中最小速度值；Fmax表示地震子波的最大頻率。