[发明专利]一种利用单程波算子的双程波叠前深度偏移方法

摘要

本发明涉及一种利用单程波算子的双程波叠前深度偏移方法，包括如下步骤：1，采用上下双检波器地震数据采集系统，形成检波器波场；2，提取地震子波信息，形成震源波场；3，读取单炮地震数据、速度模型和震源波场；4，对震源波场和检波器波场在时间域作傅里叶变换得到频率‑空间域波场数据；5，对频率‑空间域波场数据进行深度延拓计算；6，利用互相关成像原理或反射系数成像原理进行偏移成像；7，重复步骤4～6延拓计算到最大深度，保存偏移结果，计算下一炮地震数据；8，重复步骤3～7至最后一炮，输出偏移剖面。本发明，利用单程波算子进行全波动方程深度偏移，实现了快速、准确的双程波深度偏移算法，更有工业利用价值。

主权项

1.一种利用单程波算子的双程波叠前深度偏移方法，其特征在于，包括如下步骤：步骤1，通过上下双检波器地震数据采集系统获取双检波器地震数据，形成检波器波场；步骤2，从双检波器地震数据中提取地震子波信息，形成震源波场；步骤3，读取单炮地震数据、速度模型和震源波场；步骤4，对震源波场和检波器波场在时间域作傅里叶变换得到频率‑空间域波场数据；步骤5，对频率‑空间域波场数据进行深度延拓计算；步骤6，利用互相关成像原理或反射系数成像原理进行偏移成像；步骤7，重复步骤4～6延拓计算到最大深度，保存偏移结果，计算下一炮地震数据；步骤8，重复步骤3～7至最后一炮，输出偏移剖面；步骤5具体包括：输入z和z+Δz的波场值和根据下式计算z+2Δz的波场值式(5)中，为单程波算子；采用不同的近似算法来计算kz，则得到不同的双检波器叠前深度偏移算法；具体计算单程波算子中的k_z时，采用裂步傅里叶方法、傅里叶有限差分方法或高阶广义屏方法；分别形成：双检波器裂步傅里叶叠前深度偏移、双检波器傅里叶有限差分叠前深度偏移和双检波器高阶广义屏叠前深度偏移；以裂步傅里叶方法来近似单程波算子时，其中kz可以表示为：其中：v_r为介质参考速度；k_x和k_y分别表示x方向和y方向波数；A^SSF表示利用裂步傅里叶方法近似的单程波算子；表示利用裂步傅里叶方法近似计算的k_z；以高阶广义屏方法来近似单程波算子时，其中kz可以表示为：其中：A^GSP表示利用高阶广义屏方法近似的单程波算子；表示利用高阶广义屏方法近似计算的k_z；其中，ω为角频率；v为模型速度；Δz为z方向网格剖分间距，也是延拓步长；kZ为z方向的波数。