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

权利要求书

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

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

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

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

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

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

2.根据权利要求1所述的一种地震数据处理方法，其特征在于，所述对所述原始地震道数据中的振幅值进行统计处理，得到概率密度函数，包括：

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

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

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

3.根据权利要求1所述的一种地震数据处理方法，其特征在于，所述基于所述原始地震道数据中的振幅值范围，生成初始随机数集合，包括：

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

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

4.根据权利要求3所述的一种地震数据处理方法，其特征在于，所述初始随机数集合中随机数的数量的取值范围为：5000～10000。

5.根据权利要求1所述的一种地震数据处理方法，其特征在于，所述从所述初始随机数集合中确定符合所述概率密度函数的目标随机数集合，包括：

采用舍选法从所述初始随机数集合中确定符合所述概率密度函数的目标随机数集合。

6.根据权利要求1所述的一种地震数据处理方法，其特征在于，所述基于所述目标随机数集合，确定所述原始地震道数据对应的多个子波数据，包括：

基于预设子波模型和所述目标随机数集合，确定目标似然函数模型；

采用期望最大化方法确定所述目标似然函数模型中的子波参数集合中子波参数的值；其中，所述子波参数集合包括多组子波参数，所述子波参数包括子波的振幅系数、中心时间、主频和相位；

根据所述预设子波模型和所述子波参数集合中子波参数的值，确定所述原始地震道数据对应的多个子波数据；其中，所述子波参数与所述子波数据一一对应，以使得一组子波参数对应一个子波数据。

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

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

以及采用下述公式确定目标似然函数模型：

其中，L(Θ)表示目标似然函数，表示所述多个子波数据中第j个子波数据，x_z表示所述目标随机数集合中第z个随机数，P表示所述目标随机数集合中随机数的数量；表示所述子波参数集合，α_j表示第j个子波数据中的子波的振幅系数，u_j表示第j个子波数据中的子波的中心时间，ζ_j表示第j个子波数据中的子波的主频，表示第j个子波数据中的子波的相位，M表示所述多个子波数据中子波数据的数量。