[发明专利]地震子波提取方法

权利要求书

1.一种地震子波提取方法，包括：

基于井资料与地震资料，获得多个井的反射系数序列与对应的井旁地震道，其中所述多个井的反射系数序列以反射系数矩阵表示，所述对应的井旁地震道以井旁地震道矩阵表示；

基于所述反射系数矩阵与所述井旁地震道矩阵，建立关于地震子波矩阵的字典学习方程；

基于所述反射系数矩阵、所述井旁地震道矩阵与所述地震子波矩阵，建立目标函数，计算使所述目标函数的值最小的地震子波矩阵，作为最终的地震子波矩阵；

对所述最终的地震子波矩阵进行提取与处理，获得平均子波；

其中，所述目标函数为：

其中，f(x)表示目标函数，X表示反射系数矩阵、Y表示井旁地震道矩阵，D表示地震子波矩阵，γ表示权重因子，||X||_p表示X的p范数。

2.根据权利要求1所述的地震子波提取方法，其中，基于井资料与地震资料，获得多个井的反射系数序列与对应的井旁地震道，其中所述多个井的反射系数序列以反射系数矩阵表示，所述对应的井旁地震道以井旁地震道矩阵表示包括：

基于所述井资料与所述地震资料，计算每口井的反射系数序列，并将所述每口井的反射系数序列转换到时间域；

对每口井进行井-震标定，获得每口井对应的井旁地震道；

将时间域的每个反射系数序列作为列向量，获得所述反射系数矩阵，将每口井对应的井旁地震道作为列向量，获得所述井旁地震道矩阵。

3.根据权利要求1所述的地震子波提取方法，其中，所述字典学习方程为：

Y＝D·X (2)

其中，Y表示井旁地震道矩阵，X表示反射系数矩阵，D表示地震子波矩阵。

4.根据权利要求1所述的地震子波提取方法，其中，对所述最终的地震子波矩阵进行提取与处理，获得平均子波包括：

分别提取所述最终的地震子波矩阵的每一列作为一个地震子波序列，获得多个地震子波序列对应的起始位置与长度；

针对每个地震子波序列，基于其对应的起始位置与长度，对其进行傅立叶变换，获得所述长度内的所述地震子波序列的振幅谱与相位谱；

针对多个地震子波序列的振幅谱与相位谱进行平均计算，获得平均振幅谱与平均相位谱，进而获得所述平均子波。

5.根据权利要求4所述的地震子波提取方法，其中，所述平均子波为：

其中，A(ω)表示平均振幅谱，表示平均相位谱，i表示虚数单位，w(ω)表示频率域的平均子波。

6.根据权利要求5所述的地震子波提取方法，其中，还包括：将所述频率域的平均子波通过傅里叶反变换转换至时间域，获得时间域的平均子波。