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

说明书

本申请实施例公开了一种地震数据处理方法及装置。所述方法包括：对原始叠前地震道数据进行希尔伯特变换，得到希尔伯特数据；确定原始叠前地震道数据对应的目标子波；将原始地震道数据减去目标子波，得到残差地震道数据；当残差地震道数据中的振幅的绝对值的最大值大于或等于预设振幅阈值时，将残差地震道数据作为新的原始叠前地震道数据，并生成新的原始叠前地震道数据对应的新的目标子波和新的残差地震道数据，直至新的残差地震道数据中振幅的绝对值的最大值小于预设振幅阈值，并将新的目标子波以及在新的目标子波之前生成的目标子波进行叠加，将叠加的结果作为目标叠前地震道数据。本申请实施例提供的技术方案可以提高叠前地震数据的分辨率。

技术领域

本申请涉及地球物理勘探数据处理技术领域，特别涉及一种地震数据处理方法及装置。

背景技术

在地球物理勘探处理和解释过程中，叠前地震道集数据通常需经过叠加处理后形成叠后数据，然后交付给解释软件(或模块)使用。在实际生产实践中，叠前地震道集数据在叠加处理之前，不对叠前地震道集中的叠前地震道数据进行任何处理与分析。针对分辨率较低的叠前地震道集，经过叠加处理后，往往难以满足地震资料的高精度处理和储层精细刻画的要求，也不利于解释人员对目的层位的地震数据进行更加细致的解释。

发明内容

本申请实施例的目的是提供一种地震数据处理方法及装置，以提高叠前地震数据的分辨率。

为解决上述技术问题，本申请实施例提供一种地震数据处理方法及装置是这样实现的：

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

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

确定所述原始叠前地震道数据对应的目标子波；其中，所述目标子波根据预设子波模型、中心时间和目标主频确定；所述中心时间根据所述希尔伯特数据确定；所述目标主频根据所述预设子波模型、所述原始叠前地震道数据和所述中心时间确定；

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

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

优选方案中，所述中心时间按照下述方法确定：

确定所述希尔伯特数据中多个极大值数据点以及所述极大值数据点对应的采样时间；其中，所述极大值数据点与采样时间一一对应；

计算所述多个极大值数据点分别与所述目的层位之间的距离；具体地，所述极大值数据点与所述目的层位之间的距离为所述极大值数据点对应的采样时间与所述目的层位对应的预设采样时间之间的差值的绝对值；

将所述多个极大值数据点中与所述目的层位之间的距离最小的极大值数据点作为目标极大值数据点，并将所述目标极大值数据点对应的采样时间作为所述中心时间。

优选方案中，所述目标主频按照下述方法确定：

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

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