[发明专利]一种地震资料偏移成像方法及系统

摘要

本发明公开一种地震资料偏移成像方法及系统。该方法包括：获取待进行偏移成像的地震资料数据、观测系统参数、偏移速度模型参数、偏移密度模型参数和偏移参数；获取每炮对应的震源波场；获取当前每炮对应的模拟地震记录数据，获取每炮的每一个记录时刻的伴随波场；获得当前每炮对应的单炮梯度剖面；将所有炮的第一梯度剖面叠加，获得第二梯度剖面；获得下降方向剖面、最优化步长、偏移剖面、目标函数值；判断当前目标函数值是否满足收敛标准，如果是，将当前偏移剖面确定为最终的偏移剖面；否则，再次获取每炮对应的模拟地震记录数据。本发明方法及系统，可以获得高精度、高分辨率、高信噪比、振幅保真的叠前深度偏移剖面，且成像精度和分辨率高。

主权项

1.一种地震资料偏移成像方法，其特征在于，所述方法包括：获取待进行偏移成像的地震资料数据、观测系统参数、偏移速度模型参数、偏移密度模型参数和偏移参数；所述地震资料数据为多炮的观测地震记录数据；获取每炮对应的震源波场；所述每炮对应的震源波场包括每炮的所有记录时刻的震源波场；获取当前每炮对应的模拟地震记录数据，所述模拟地震记录数据根据前一次每炮对应的模拟地震记录数据利用非线性反偏移算子获得；根据当前每炮对应的观测地震记录数据和对应的模拟地震记录数据，获取每炮的每一个记录时刻的伴随波场；根据每炮对应的震源波场和伴随波场，获得当前每炮对应的第一梯度剖面；所述第一梯度剖面为每炮对应的单炮梯度剖面；根据所述观测系统参数和所述偏移参数，将所有炮的第一梯度剖面叠加，获得当前第二梯度剖面；利用最优化反演算法获得当前下降方向剖面；具体包括：利用δg_i＝g_i‑g_i‑1获得第i次的第二梯度剖面的增量δg_i，其中g_i为第i次的第二梯度剖面，g_i‑1为第i‑1次的第二梯度剖面；利用获得特征参数A_i；其中α_i‑1为第i‑1次的最优化步长，r_i‑1为第i‑1次的下降方向剖面，g_i为第i次的剖面梯度；利用获得特征参数B_i；利用δm_i＝ξ_iα_i‑1r_i‑1获得第i次的偏移剖面增量δm_i；其中利用获得特征参数C_i，τ为大于1的正数；利用获得第i次的下降方向剖面r_i；利用最优化步长计算方法获取当前最优化步长；根据所述当前最优化步长和所述当前下降方向剖面获得当前偏移剖面；根据获得当前目标函数值E；其中δd表示当前地震记录数据的残差；T表示转置运算；判断所述当前目标函数值是否满足收敛标准，得到第一判断结果；当所述第一判断结果表示所述当前目标函数值满足收敛标准时，将所述当前偏移剖面确定为偏移成像最终的偏移剖面；当所述第一判断结果表示所述当前目标函数值不满足收敛标准时，返回获取每炮对应的模拟地震记录数据步骤，再次获取每炮对应的模拟地震记录数据。