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

说明书

本申请实施例公开了一种地震数据处理方法及装置。所述方法包括：对原始地震道数据中的振幅值进行统计处理，得到概率密度函数；确定符合概率密度函数的目标随机数集合；确定多个子波数据，将多个子波数据进行第一叠加，得到拟合地震道数据；确定残差地震道数据；当残差地震道数据中的振幅的绝对值的最大值大于或等于预设振幅阈值时，将残差地震道数据乘以指定特征参数得到的结果作为新的原始地震道数据，生成新的拟合地震道数据和新的残差地震道数据，直至新的残差地震道数据中振幅的绝对值的最大值小于预设振幅阈值，将新的拟合地震道数据以及在其之前生成的拟合地震道数据进行第二叠加，得到目标地震道数据。可以提高地震道数据的分辨率。

技术领域

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

背景技术

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

发明内容

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

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

一种地震数据处理方法，提供有原始地震道数据；所述方法包括：

对所述原始地震道数据中的振幅值进行统计处理，得到概率密度函数；

基于所述原始地震道数据中的振幅值范围，生成初始随机数集合，从所述初始随机数集合中确定符合所述概率密度函数的目标随机数集合；

基于所述目标随机数集合，确定所述原始地震道数据对应的多个子波数据，并将所述多个子波数据进行第一叠加，将第一叠加的结果作为拟合地震道数据；

将所述原始地震道数据减去所述拟合地震道数据，得到残差地震道数据；

当所述残差地震道数据中的振幅的绝对值的最大值大于或等于预设振幅阈值时，将所述残差地震道数据乘以指定特征参数得到的结果作为新的原始地震道数据，并生成所述新的原始地震道数据对应的新的拟合地震道数据和新的残差地震道数据，直至所述新的残差地震道数据中振幅的绝对值的最大值小于所述预设振幅阈值，并将所述新的拟合地震道数据以及在所述新的拟合地震道数据之前生成的拟合地震道数据进行第二叠加，并将第二叠加的结果作为目标地震道数据；其中，所述指定特征参数与迭代次数相关联。

优选方案中，所述对所述原始地震道数据中的振幅值进行统计处理，得到概率密度函数，包括：

将所述原始地震道数据中数值为负数的振幅值替换为零，得到第一地震道数据；

当所述第一地震道数据中的指定采样点的振幅值大于预设振幅阈值时，将所述预设振幅阈值作为所述指定采样点的振幅值，得到第二地震道数据；

对所述第二地震道数据中的振幅值进行统计处理，得到所述概率密度函数。

优选方案中，所述基于所述原始地震道数据中的振幅值范围，生成初始随机数集合，包括：

采用线性同余法生成多个均匀分布的随机数；其中，所述多个均匀分布的随机数的值在所述原始地震道数据中的振幅值范围内；

将所述多个均匀分布的随机数构成的集合作为所述初始随机数集合。

优选方案中，所述初始随机数集合中随机数的数量的取值范围为：5000～10000。