[发明专利]一种最小二乘偏移成像方法及系统

权利要求书

1.一种最小二乘偏移成像方法，其特征在于，包括以下步骤：

(1)采集野外的观测数据，计算得到炮记录数据，预设残差阈值为设定低值；

(2)利用子波估计方法计算震源子波函数；

(3)将深度域偏移速度场变换为变换域偏移速度场计算得到偏移成像结果以及反偏移模拟数据；

(4)计算所述反偏移模拟数据与所述炮记录数据之间的残差；

(5)计算梯度，根据梯度迭代地更新偏移成像结果，直到所述残差值小于所述设定低值，输出最终成像结果。

2.根据权利要求1所述的最小二乘偏移成像方法，其特征在于，所述最终成像结果为变换域偏移速度场的成像结果，或者为将变换域偏移速度场的成像结果变换回深度域偏移速度场的成像结果。

3.根据权利要求1或2所述的最小二乘偏移成像方法，其特征在于，所述变换域偏移速度场为伪时间域偏移速度场，将深度域偏移速度场变换为伪时间域偏移速度场的变换公式为：

其中，ν为偏移速度，x为横向坐标，z为垂向深度，τ为单程垂直旅行时间。

4.根据权利要求1所述的最小二乘偏移成像方法，其特征在于，所述设定低值为1％。

5.根据权利要求3所述的最小二乘偏移成像方法，其特征在于，计算伪时间域偏移速度场的偏移成像结果过程中，伪时间域反偏移过程的波场延拓公式为：

其中，ν为偏移速度，ξ＝x，η＝τ，x为横向坐标，τ为单程垂直旅行时间，Δs²为慢度场平方的扰动，为背景慢度场的平方，p₀为背景波场，p_s为扰动波场，t为延拓时间。

6.根据权利要求5所述的最小二乘偏移成像方法，其特征在于，在计算梯度过程中，伪时间域偏移过程的波场延拓公式为：

其中，ν为偏移速度，ξ＝x，η＝τ，x为横向坐标，τ为单程垂直旅行时间，p为延拓波场，t为延拓时间。

7.根据权利要求6所述的最小二乘偏移成像方法，其特征在于，将伪时间域偏移速度场的成像结果变换回深度域偏移速度场的成像结果的变换公式为：

其中，ν为偏移速度，x为横向坐标，z为垂向深度，τ为单程垂直旅行时间。

8.一种最小二乘偏移成像系统，包括控制器，所述控制器包括处理器与存储器，其特征在于，所述处理器运行存储在所述存储器中的计算机程序，以实现权利要求1-7中任意一项所述的方法。