[发明专利]一种基于多源混合动态编码的最小二乘逆时偏移方法

摘要

一种基于多源混合动态编码的最小二乘逆时偏移方法,在多源编码最小二乘逆时偏移中，首先利用多源编码矩阵将单炮道集分组处理转化为超炮道集，然后利用超炮道集来定义估计反射系数的目标函数，在最小二乘逆时偏移结果中引入交叉项噪音，将每个阶段的编码矩阵选作正交矩阵，每个子问题的形成采用混合编码方式来实现。本发明基于双向波波动方程理论和反演理论的最小二乘逆时偏移方法能够提供振幅更加保真，成像更加聚焦，成像假象更少的地震偏移结果，从而降低地震解释的难度和油藏开发的风险。同时本发明的多源编码方法在最小二乘逆时偏移迭代的过程更好地实现计算效率的提升和编码的交叉项噪音的压制。

主权项

一种基于多源混合动态编码的最小二乘逆时偏移方法,其特征是：在多源编码最小二乘逆时偏移中，首先利用多源编码矩阵将单炮道集分组处理转化为超炮道集，然后利用超炮道集来定义估计反射系数的目标函数，Je=||dsupobs-dsupmod||2=||Bdobs-BLm||2,---(6)]]>其中，代表观测的多源超炮道集；代表模拟的多源超炮道集；L代表正演算子；m为模型；B代表多源编码矩阵，设定BTB≠I，在最小二乘逆时偏移结果中引入交叉项噪音，定义如下动态目标泛函，Jie＝||Bidobs‑BiLm||2             (7)其中，i代表反演过程中的阶段编号；将每个阶段的编码矩阵选作正交矩阵，它的每一个子问题定义的误差泛函为，Jje＝||bjdobs‑bjLm||2           (8)其中，j代表反演过程中的子问题编号，bj为子问题中的多源编码矩阵；每个子问题的形成采用如下混合编码方式来实现，usupjs(x,ω)=Σk=1nsigei(fk,j+ωτk,i)usk(x,ω)---(9)]]>usupjr(x,ω)=Σk=1nsigei(fk,j+ωτk,i)urk,j(x,ω)---(10)]]>其中，代表合成的超炮源波场；usk代表单炮源波场；k代表炮号；代表合成的超炮检波点残差波场；urk,j代表单炮残差波场，代表反演的每个阶段对单炮道集的随机时间延迟；代表在每个阶段内形成每个子问题所用的正交编码向量。