[发明专利]一种基于阻抗敏感核函数的最小二乘逆时偏移成像方法

权利要求书

1.一种基于阻抗敏感核函数的最小二乘逆时偏移成像方法，包括：

获取初始输入数据，所述初始输入数据包括：偏移速度模型v₀(x)、偏移密度模型ρ₀、震源函数子波w(t)、观测数据d^obs、预设最大迭代次数N_iter和冗余误差ε；

根据所述初始输入数据和预先确定的伴随偏移算子确定反射系数模型m_k的初始值m₀；

对反射系数模型m_k进行如下迭代更新，直至满足预设的迭代终止条件：

对第k次获得的反射系数m_k，进行线性正演模拟，得到合成地震数据d^syn＝Lm_k，其中L为线性正演算子，0≤k≤N_iter；

计算数据残差r＝d^syn-d^obs，根据该数据残差和所述伴随偏移算子确定第k+1次迭代的更新梯度

根据所述更新梯度计算第k+1次的更新步长α_k+1、共轭步长β_k+1和更新方向h_k+1，对反射系数模型进行第k+1次的更新m_k+1＝m_k+α_k+1h_k+1；

当迭代停止时，输出当前的反射系数模型，作为最终成像结果；

其中，所述线性正演算子L的形式为：

其中，x为成像空间位置，x_s为震源位置，和x_r为检波点位置，t为波场传播时间，T为记录长度，v₀(x)为偏移速度模型，z₀(x)为背景阻抗模型，p(x_s,x,t)为正向传播波场，为伴随波场，为空间梯度算子，f(x_s,t)为震源函数，为虚拟震源；

所述伴随偏移算子的形式为：

其中，为伴随波场，r(x_r,t)为数据残差，δ为狄拉克δ函数，K(x_s,x_r,x)为阻抗敏感核函数。

2.如权利要求1所述的方法，第k+1次更新步长为第k+1次的共轭步长为第k+1次迭代的更新方向为h_k+1＝H^-1g_k+1+β_k+1h_k，其中，H^-1为震源照明预条件算子，表示点乘运算，g_k为第k次迭代时的更新梯度。

3.如权利要求1所述的方法，所述预设的迭代终止条件，包括：数据残差L₂模小于所述冗余误差ε，或者，当前迭代次数超过最大迭代次数N_iter。

4.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其中，所述处理器执行所述程序时实现如权利要求1至3任一所述的方法。