[发明专利]地震数据处理方法和装置

权利要求书

1.一种地震数据处理方法，其特征在于，包括：

获取目标区域的地震数据；

从所述地震数据中提取第一地震子波；

在时间域中对所述第一地震子波进行压缩处理，以获取第二地震子波；

根据所述地震数据、所述第一地震子波、所述第二地震子波建立反演目标函数，其中，所述反演目标函数为基于反演机理建立的用于进行高分辨率地震数据处理的函数；

根据所述反演目标函数，确定处理后的地震数据；

其中，在时间域中对所述第一地震子波进行压缩处理，以获取第二地震子波，包括：

根据傅里叶变换伸缩性，按照以下算式对所述第一地震子波进行压缩处理，以获取所述第二地震子波：

w_h(t)＝w_t(at)

其中，w_h(t)为第二地震子波，w_t(t)为第一地震子波，t为时间点的编号，a为压缩因子；

其中，所述第二地震子波为高频地震子波；

其中，根据所述地震数据、所述第一地震子波、所述第二地震子波建立反演目标函数，包括：按照以下公式，建立所述反演目标函数：

其中，minimize表示求取最小值，表示求取L2范数，y表示反演目标函数的解析解，W为基于第一地震子波的第一地震子波矩阵，W_h为基于第二地震子波的第二地震子波矩阵，x为地震数据的列向量形式，μ为正则化参数。

2.根据权利要求1所述的方法，其特征在于，根据所述第二地震子波，按照以下公式确定所述基于第二地震子波的第二地震子波矩阵：

其中，W_h为基于第二地震子波的第二地震子波矩阵，w_h(i)为第二地震子波中编号为i的时间点的振幅值，n为第二地震子波时间点的总数。

3.根据权利要求1所述的方法，其特征在于，根据所述反演目标函数，确定处理后的地震数据，包括：

根据所述反演目标函数，求取所述反演目标函数的导数为0时所对应的解析解，并将所述解析解作为所述处理后的地震数据。

4.根据权利要求3所述的方法，其特征在于，按照以下公式求取所述反演目标函数的导数为0时所对应的解析解：

y′＝(W^TW+μI)^-1W^Td

其中，(·)^-1表示矩阵求逆，y′为反演目标函数的导数为0时所对应的解析解，W为基于第一地震子波的第一地震子波矩阵，W^T为基于第一地震子波的第一地震子波矩阵转置，I为单位矩阵，d为模型向量，μ为正则化参数。

5.根据权利要求4所述的方法，其特征在于，按照以下方式确定所述模型向量：

d＝W_h·x

其中，d为模型向量，W_h为基于第二地震子波的第二地震子波矩阵，x为地震数据的列向量。