[发明专利]一种确定地震子波的方法

摘要

本发明是石油物探地震资料处理中确定地震子波的方法，在原有的目标函数中增加了确定地震子波能量最小的约束项，考虑了地震道中的噪音对确定地震子波精度的影响，通过寻找地震道残差能量与地震子波能量的最佳“平衡点”，确定出反射系数褶积矩阵奇异值的最佳截断阶数，通过对反射系数褶积矩阵奇异值截断来减弱地震道的噪音影响，确定地震子波。本发明不影响到较大奇异值对子波确定的贡献，与常用方法相比更有利于减弱地震道中噪音对子波确定精度的影响，有利于地震子波确定精度的提高。

主权项

1.一种确定地震子波的方法，特点是通过以下步骤实现：1）采集测井数据和地震道数据以及先验地震子波；2）利用先验地震子波的时间域反射系数制作合成记录，在时窗范围内对时间域反射系数进行重新标定；3）给定目标地震子波的长度，确定地震道数据时窗和反射系数时窗范围；4）构建反射系数褶积矩阵；5）构建地震道列向量X；6）按照以下公式确定目标地震子波目标函数；E={||RW(k)-X||22+||W(k)||22}→min---(9)]]>式中：R为反射系数褶积矩阵；W^(k)表示对反射系数褶积矩阵进行k阶截断情况下地震子波；X为地震道的列向量；k为反射系数褶积矩阵奇异值的截断阶数；7）按照以下公式计算不同阶数k奇异值截断情况下的地震子波W^(k)：W(k)=Σi=0kuiTxiδivi---(10)]]>式中：δ_i为反射系数褶积矩阵R的奇异值；为反射系数褶积矩阵左奇异矩阵的列向量转置；v_i为反射系数褶积矩阵右奇异矩阵列向量；x_i为地震道的列向量的各个元素；8）计算地震道残差的能量和地震子波W^(k)能量，确定最佳k值：地震道残差能量Es(k)=||RW(k)-X||22---(11),]]>和地震子波能量Ew(k)=||W(k)||22;---(12)]]>式中：为k阶截奇异值断情况下震道残差能量；式中：为k阶截奇异值截断情况地震子波能量：对不同子波对应的地震子波能量和地震道残差能量做交会图，选取地震子波能量和地震道残差能量交会拐点作为两者之间的最佳k值；9）利用最佳k值和地震子波W^(k)公式确定最佳地震子波；10）对多个地震道重复步骤2）～9）得到其所对应的最佳地震子波，然后再采用多子波平均得到井点位置附近的目标地震子波。