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

权利要求书

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）得到其所对应的最佳地震子波，然后再采用多子波平均得到井点位置附近的目标地震子波。