[发明专利]分相位地震数据的处理方法、装置和服务器

权利要求书

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

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

对所述反射地震数据进行基于Gabor变换的稀疏约束反演谱分解处理，得到对应的时频振幅谱数据和时频相位谱数据；

在频率域上，根据预设的相位误差限值、所述时频振幅谱数据和所述时频相位谱数据，进行分相位重构，以得到分相位重构后的数据体；

对所述分相位重构后的数据体进行基于Gabor变换的稀疏约束反演谱分解的逆处理，得到分相位地震数据；其中，所述分相位地震数据包含有不同相位分量的地震数据；

其中，对所述反射地震数据进行基于Gabor变换的稀疏约束反演谱分解处理，包括：构建子波库；其中，所述子波库包括多个频率的Gabor变换基函数；根据所述反射地震数据、所述子波库，构建基于Gabor变换的稀疏约束反演的目标函数；求解所述目标函数，以得到对应的时频数据；根据所述时频数据，确定出所述时频振幅谱数据和所述时频相位谱数据。

2.根据权利要求1所述的方法，其特征在于，在得到分相位地震数据之后，所述方法还包括：

从所述分相位地震数据中，筛选出指定相位分量的地震数据；

根据所述指定相位分量的地震数据，确定出目标区域的储层特征。

3.根据权利要求1所述的方法，其特征在于，在得到分相位地震数据之后，所述方法还包括：

对所述分相位地震数据进行全相位叠加处理，以重构地震记录。

4.根据权利要求1所述的方法，其特征在于，按照以下算式构建基于Gabor变换的稀疏约束反演的目标函数：

其中，J为目标函数的函数值，G为基于所述子波库所构建的褶积矩阵，m为时频数据，s为反射地震数据的矩阵形式，γ为正则化参数，表示L2范数的平方，||·||₁表示L1范数。

5.根据权利要求4所述的方法，其特征在于，求解所述目标函数，包括：

获取噪声数据；

根据所述噪声数据，调整所述目标函数中的正则化参数，以得到调整后的目标函数；

通过迭代软阈值算法，求解所述调整后的目标函数。

6.根据权利要求5所述的方法，其特征在于，根据所述噪声数据，调整所述目标函数中的正则化参数，包括：

在确定所述噪声数据大于第一预设阈值的情况下，增大所述正则化参数的数值；

在确定所述噪声数据小于第二预设阈值的情况下，减小所述正则化参数的数值。

7.根据权利要求1所述的方法，其特征在于，在频率域上，根据预设的相位误差限值、所述时频振幅谱数据和所述时频相位谱数据，进行分相位重构，以得到分相位重构后的数据体，包括：

根据时频相位谱数据、预设的相位值，以及预设的相位误差限值，确定出符合要求的相位；

保留所述时频振幅谱数据中的符合要求的相位处的时频振幅数据，以得到分相位重构后的数据体。

8.一种分相位地震数据的处理装置，其特征在于，包括：

获取模块，用于获取目标区域的反射地震数据；

第一处理模块，用于对所述反射地震数据进行基于Gabor变换的稀疏约束反演谱分解处理，得到对应的时频振幅谱数据和时频相位谱数据；

重构模块，用于在频率域上，根据预设的相位误差限值、所述时频振幅谱数据和所述时频相位谱数据，进行分相位重构，以得到分相位重构后的数据体；

第二处理模块，用于对所述分相位重构后的数据体进行基于Gabor变换的稀疏约束反演谱分解的逆处理，得到分相位地震数据；其中，所述分相位地震数据包含有不同相位分量的地震数据；

其中，所述第一处理模块具体用于构建子波库；其中，所述子波库包括多个频率的Gabor变换基函数；根据所述反射地震数据、所述子波库，构建基于Gabor变换的稀疏约束反演的目标函数；求解所述目标函数，以得到对应的时频数据；根据所述时频数据，确定出所述时频振幅谱数据和所述时频相位谱数据。

9.一种服务器，其特征在于，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现权利要求1至7中任一项所述方法的步骤。