[发明专利]高分辨率转换波裂缝预测方法

摘要

本发明提供了一种高分辨率转换波裂缝预测方法，所述方法包括：采集转换波地震数据，进行预处理和叠前偏移后得到转换波偏移剖面；在偏移剖面上建立稀疏约束和波阻抗约束的多薄层目标函数；对多薄层目标函数利用频域映射与L1范数联合优化算法求取高精度反射系数；构建高频子波同时保持该高频子波具有确定的频率宽度，用高频子波与高精度反射系数得到高分辨率转换波偏移剖面；利用改进的第三代特征值相干体计算法，对高分辨率转换波偏移剖面进行相干数据体切片提取，从而实现裂缝预测。本发明能够使得提取的裂缝预测剖面分辨率更高、裂缝走向更加清晰、细节刻画更加明显。

主权项

1.一种高分辨率转换波裂缝预测方法，所述方法包括：采集转换波地震数据，进行预处理和叠前偏移后得到转换波偏移剖面，在转换波偏移剖面上取n道地震记录x(x,t)，x＝1,2,…,n,t表示时间；在偏移剖面上建立稀疏约束和波阻抗约束的多薄层目标函数；对多薄层目标函数利用频域映射与L1范数联合优化算法求取高精度反射系数；构建高频子波同时保持该高频子波具有确定的频率宽度，用高频子波与高精度反射系数得到高分辨率转换波偏移剖面；利用改进的第三代特征值相干体计算法，对高分辨率转换波偏移剖面进行相干数据体切片提取，从而实现裂缝预测；其中，所述目标函数为下式(1)：其中，r表示反射系数，a_e代表反射系数的偶分量比例，a_o代表反射系数的奇分量比例，r_e为反射系数的偶分量，r_o为反射系数的奇分量，A为变换矩阵，b为地震数据和子波相关矩阵，λ为稀疏因子，为阻抗因子，C为积分算子矩阵，ξ_p为纵波波阻抗矩阵，其中，其中，t为时间，t1、t2、……、tN分别表示第1、2、……、N个时间，f为频率，f1、f2、……、fM分别表示第1、2、……、M个分析频率范围内的频率，T1、T2……、TN/2分别表示1、2、……、N/2个时间厚度，Δt表示时移量，Re表示频率域下反射系数R(f)的实部，Im表示频率域下反射系数R(f)的虚部，S表示频率域下地震记录x(x,t)，W表示频率域下地震子波；其中，t₀表示初始时间，t_M‑1表示结束时间，纵波相对波阻抗ε_t为：I(t)为纵波波阻抗值，I(t₀)表示纵波初始波阻抗值，r为反射系数，η表示在时间t₀与时间t中取值符合；C的离散形式表示为