[发明专利]绕射波成像方法、装置和电子设备

权利要求书

1.一种绕射波成像方法，其特征在于，包括：

获取地震数据和速度数据；

根据所述速度数据对预设的地震模型进行仿真处理，得到模拟地震数据；所述地震模型包括反射模型和绕射模型；

基于所述地震数据和所述模拟地震数据，得到波场残差目标函数和波场残差值；其中，根据所述地震数据和所述模拟地震数据，利用L2范数构建所述波场残差目标函数，并根据构建的所述波场残差目标函数计算所述波场残差值；

对所述波场残差值进行偏移处理，得到所述波场残差目标函数的梯度；所述梯度包括反射梯度和绕射梯度；

根据所述反射梯度对所述反射模型进行更新处理，得到更新后的反射模型；其中，根据所述反射梯度，得到所述反射模型的更新方向；根据所述更新方向和预设的更新步长对所述反射模型进行更新处理，得到所述更新后的反射模型；其中，所述预设的更新步长采用内插法进行计算；

根据所述绕射梯度对所述绕射模型进行更新处理，得到更新后的绕射模型；

根据所述更新后的反射模型和所述更新后的绕射模型对所述波场残差目标函数进行迭代更新处理，得到绕射波成像结果；其中，根据所述更新后的反射模型和所述更新后的绕射模型对所述波场残差目标函数进行更新，并根据更新后的波场残差目标函数，得到反射梯度和绕射梯度，再次不断对所述反射模型和所述绕射模型进行更新，以及根据预先设置的迭代次数对所述波场残差目标函数进行迭代更新处理，得到最优的反射模型和绕射模型，从而得到绕射波成像结果；

其中，所述根据所述速度数据对预设的地震模型进行仿真处理，得到模拟地震数据，包括：

根据所述速度数据分别对所述反射模型和所述绕射模型进行仿真处理，得到相对应的反射波波场数据和绕射波波场数据；

分别对所述反射波波场数据和所述绕射波波场数据进行偏移处理，得到相对应的反射波偏移数据和绕射波偏移数据；

分别对所述反射波偏移数据和所述绕射波偏移数据进行正演处理，得到相对应的反射波正演数据和绕射波正演数据；

对所述反射波正演数据和所述绕射波正演数据进行加法处理，得到所述模拟地震数据。

2.根据权利要求1所述的方法，其特征在于，所述基于所述地震数据和所述模拟地震数据，得到波场残差目标函数和波场残差值，包括：

利用下述公式对所述地震数据和所述模拟地震数据进行处理，得到所述波场残差目标函数：

其中，d_mod表示所述模拟地震数据；d_obs表示所述地震数据；m^diff表示所述绕射模型；m^ref表示所述反射模型；α表示预设的正则化参数；

利用下述公式对所述地震数据和所述模拟地震数据进行处理，得到所述波场残差值：

e(m^ref，m^diff)＝d_obs-d_mod。

3.根据权利要求2所述的方法，其特征在于，所述对所述波场残差值进行偏移处理，得到所述波场残差目标函数的梯度，包括：

利用下述公式对所述波场残差值进行偏移处理，得到所述反射梯度：

g(m^ref)＝L_refe(m^ref，m^diff)；

其中，L_ref为反射偏移算子，表示对所述波场残差值进行反射偏移；

利用下述公式对所述波场残差值进行偏移，得到所述绕射梯度：

g(m^diff)＝L_diffe(m^ref，m^diff)+α*m^diff；

其中，L_diff为绕射偏移算子，表示对所述波场残差值进行绕射偏移。