[发明专利]地震偏移成像方法及装置

权利要求书

1.一种地震偏移成像方法，其特征在于，包括：

基于待处理区域的地震观测数据和待求解模拟数据建立波场残差目标函数；

结合当前反射系数模型和所述波场残差目标函数确定当前波场残差值，其中，所述当前反射系数模型用于确定所述待求解模拟数据；

基于时移成像条件的偏移方法对所述当前波场残差值进行偏移，以确定得到目标函数梯度；

基于所述目标函数梯度采用共轭梯度法确定更新后的反射系数模型；

将所述更新后的反射系数模型作为所述当前反射系数模型进行更新迭代，直至达到预设迭代次数，得到目标反射系数模型，进而得到偏移成像结果。

2.根据权利要求1所述的方法，其特征在于，基于待处理区域的地震观测数据和待求解模拟数据建立波场残差目标函数包括：

获取所述待处理区域的地震观测数据和偏移速度体；

通过L2范数建立所述地震观测数据与所述待求解模拟数据之间的波场残差目标函数。

3.根据权利要求1所述的方法，其特征在于，结合当前反射系数模型和所述波场残差目标函数确定当前波场残差值包括：

根据模拟数据计算算式d_mod＝Lm_k计算所述待求解模拟数据，其中，d_mod表示所述待求解模拟数据，L表示Born近似条件下依赖于速度信息的地震波场正演算子，m_k表示所述当前反射系数模型；

根据所述波场残差目标函数的计算算式计算所述当前波场残差值，其中，E(m_k)表示所述当前波场残差值，d_mod表示所述待求解模拟数据，d_obs表示所述地震观测数据。

4.根据权利要求2所述的方法，其特征在于，基于时移成像条件的偏移方法对所述当前波场残差值进行偏移，以确定得到目标函数梯度包括：

通过所述偏移速度体将检波点波场进行反向延拓，并通过所述偏移速度体将震源波场进行正向延拓；

基于时移成像条件对反向延拓后的检波点波场和正向延拓后的震源波场进行成像，得到关于时移量的偏移剖面；

在所述关于时移量的偏移剖面中将零时刻的偏移剖面作为所述目标函数梯度。

5.根据权利要求4所述的方法，其特征在于，基于时移成像条件对所述反向延拓后的检波点波场和所述正向延拓后的震源波场进行成像包括：

根据基于时移成像条件的成像算式R(m_k,t,τ)＝P_r(m_k,t+τ)*P_s(m_k,t-τ)进行成像，得到所述关于时移量的偏移剖面，其中，R(m,t,τ)表示所述关于时移量的偏移剖面，*表示时间域互相关，t表示时间，τ表示时移成像条件中的时移量，P_r表示所述反向延拓后的检波点波场，P_s表示所述正向延拓后的震源波场，m_k表示所述当前反射系数模型。

6.根据权利要求1所述的方法，其特征在于，基于所述目标函数梯度采用共轭梯度法确定更新后的反射系数模型包括：

基于所述目标函数梯度确定更新步长和更新方向；

根据反射系数模型更新算式m_k+1＝m_k+α_kd_k计算所述更新后的反射系数模型，其中，m_k+1表示所述更新后的反射系数模型，m_k表示当前反射系数模型，α_k表示所述更新步长，d_k表示所述更新方向。