[发明专利]基于转换波地震数据求取高精度反射系数的方法

摘要

本发明提供了一种基于转换波地震数据求取高精度反射系数的方法，所述方法包括：采集转换波地震数据，进行预处理和叠前偏移后得到转换波偏移剖面；在偏移剖面上建立稀疏约束和波阻抗约束的多薄层目标函数；对多薄层目标函数利用频域映射与L1范数联合优化算法求取高精度反射系数。本发明的方法能够基于转换波地震数据构建更加合理的多薄层目标函数并获得高精度反射系数，提高了目标函数求解精度和抗噪性，而且所获得的高精度反射系数能够为后续高分辨率剖面重构奠定坚实基础。

主权项

一种基于转换波地震数据求取高精度反射系数的方法，所述方法包括：采集转换波地震数据，进行预处理和叠前偏移后得到转换波偏移剖面，在转换波偏移剖面上取n道地震记录x(x,t)，x＝1,2,…,n；在偏移剖面上建立稀疏约束和波阻抗约束的多薄层目标函数；对多薄层目标函数利用频域映射与L1范数联合优化算法求取高精度反射系数；其中，所述目标函数为下式(1)：$O\left(r_e,r_o\right)=\left|\right|\left[\begin{array}{c}{a}_e(b_e-A_e\×r_e)\\ a_o(b_o-A_o\×r_o)\end{array}\right]|{|}_2^2+\mathrm{\λ}||r_e+r_o|{|}_1+\∂\left|\right|\frac{{(Cr-{\mathrm{\ξ}}_p)}^T(Cr-{\mathrm{\ξ}}_p)}{2}|{|}_1,$其中，r表示反射系数，a_e代表反射系数的偶分量比例，a_o代表反射系数的奇分量比例，r_e为反射系数的偶分量，r_o为反射系数的奇分量，A为变换矩阵，b为地震数据和子波相关矩阵，λ为稀疏因子，为阻抗因子，C为积分算子矩阵，ξ_p为纵波波阻抗矩阵。