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

权利要求书

1.一种地震数据处理方法，其特征在于，提供有目的层位的原始叠前地震道数据；所述方法包括：

对所述原始叠前地震道数据进行希尔伯特变换，得到希尔伯特数据；

确定所述原始叠前地震道数据对应的目标子波；其中，所述目标子波根据预设子波模型、中心时间和目标主频确定；所述目标主频根据所述预设子波模型、所述原始叠前地震道数据和所述中心时间确定；所述中心时间按照下述方法确定：确定所述希尔伯特数据中多个极大值数据点以及所述极大值数据点对应的采样时间，其中，所述极大值数据点与采样时间一一对应，并计算所述多个极大值数据点分别与所述目的层位之间的距离，所述极大值数据点与所述目的层位之间的距离为所述极大值数据点对应的采样时间与所述目的层位对应的预设采样时间之间的差值的绝对值，并将所述多个极大值数据点中与所述目的层位之间的距离最小的极大值数据点作为目标极大值数据点，并将所述目标极大值数据点对应的采样时间作为所述中心时间；

将所述原始地震道数据减去所述目标子波，得到残差地震道数据；

当所述残差地震道数据中的振幅的绝对值的最大值大于或等于预设振幅阈值时，将所述残差地震道数据作为新的原始叠前地震道数据，并生成所述新的原始叠前地震道数据对应的新的目标子波和新的残差地震道数据，直至所述新的残差地震道数据中振幅的绝对值的最大值小于所述预设振幅阈值，并将所述新的目标子波以及在所述新的目标子波之前生成的目标子波进行叠加，并将叠加的结果作为目标叠前地震道数据。

2.根据权利要求1所述的一种地震数据处理方法，其特征在于，所述目标主频按照下述方法确定：

根据所述预设子波模型和所述中心时间，确定所述原始叠前地震道数据对应的初始子波模型；其中，所述初始子波模型中子波的相位为零；

确定所述目标极大值数据点对应的瞬时频率，并根据所述瞬时频率确定目标主频范围；

遍历所述目标主频范围，以获取初始主频，并将所述初始主频作为所述初始子波模型的主频，得到所述原始叠前地震道数据对应的初始子波，并计算所述初始子波与所述原始叠前地震道数据的相似度；其中，所述目标主频范围中的初始主频与相似度一一对应；所述目标主频范围中的多个初始主频对应多个相似度；

将所述多个相似度中数值最大的相似度对应的主频作为所述目标主频。

3.根据权利要求2所述的一种地震数据处理方法，其特征在于，所述计算所述初始子波与所述原始叠前地震道数据的相似度，包括：

计算所述初始子波与所述原始叠前地震道数据之间的内积，并将内积结果作为所述相似度。

4.根据权利要求2所述的一种地震数据处理方法，其特征在于，采用下述公式确定目标主频范围：

(f_k-f_T)≤f≤(f_k+f_T)

其中，f表示所述目标主频范围内的任一主频，f_k表示所述瞬时频率。

5.根据权利要求1所述的一种地震数据处理方法，其特征在于，采用下述公式表征所述预设子波模型：

其中，g_γ(x)表示所述预设子波模型，u表示子波的中心时间，ζ表示子波的主频，表示子波的相位，x表示采样时间。

6.根据权利要求1所述的一种地震数据处理方法，其特征在于，所述预设振幅阈值为所述原始叠前地震道数据中的振幅的绝对值的最大值乘以5％。