[发明专利]一种地震数据分解方法、系统及存储介质与终端

说明书

本公开提供了地震数据分解方法、系统及存储介质与终端，通过确定地震波形中绝对值最大的振幅值，该振幅值对应的反射时间即为零相位雷克子波的中心时间；以最大绝对值振幅对应的反射时间为中心点，确定雷克子波的主周期；主周期计算雷克子波的主频与幅度值，根据中心时间、主频和幅度值计算零相位雷克子波波形；从地震波形中减去所得的子波分量，得到残差信号；将残差信号作为输入信号，继续循环迭代，直到残差信号的能量小于给定的门槛值为止。不再需要对地震数据进行希尔伯特变换，不需要将地震数据映射到复数域。

技术领域

本公开涉及一种地震数据分解方法、系统及存储介质与终端。

背景技术

本部分的陈述仅仅是提供了与本公开相关的背景技术信息，不必然构成在先技术。

地震数据稀疏表示和时频谱分析技术在油气储层表征和预测中应用广泛，其中匹配追踪时频谱具有极高的分辨率，能更精确地表征小尺度地质目标，例如厚度较小的碎屑岩岩储层或碳酸盐岩孔洞发育带。

匹配追踪方法存在的主要问题是计算量庞大，目前发展的地震数据快速匹配追踪方法已被用于上述类型的储层预测中，一般是利用地震数据的瞬时振幅为先验信息来确定子波的中心时间，利用瞬时频率为先验信息来确定子波的主频，加快了地震数据匹配追踪分解过程。但为了利用这些先验信息，需要首先对地震数据进行希尔伯特变换，将其映射到复数域中，然后在复数域中完成分解过程，计算仍然很复杂，在确保地震数据分解与重构精度的前提下，匹配追踪分解速度有望进一步提高。

发明内容

本公开为了解决上述问题，提出了一种地震数据分解方法、系统及存储介质与终端，本公开利用实数域零相位雷克子波对地震数据体进行稀疏分解，提供在实数域直接对地震波形进行自适应子波分解的过程，不需要地震数据的瞬时振幅和瞬时频率信息，也不需要将地震数据进行希尔伯特变换，计算过程简洁，能够快速完成大型地震数据体的稀疏分解与精确重构。

根据一些实施例，本公开采用如下技术方案：

一种地震数据分解方法，包括以下步骤：

步骤一、确定地震波形中绝对值最大的振幅值，该振幅值对应的反射时间即为零相位雷克子波的中心时间；

步骤二、以步骤一的最大绝对值振幅对应的反射时间为中心点，确定雷克子波的主周期；

步骤三、根据步骤二的主周期计算雷克子波的主频与幅度值；

步骤四、根据中心时间、主频和幅度值计算零相位雷克子波波形；

步骤五、从地震波形中减去步骤四所得的子波分量，得到残差信号；将残差信号作为输入信号，转至步骤一继续循环迭代，直到残差信号的能量小于给定的门槛值为止。

作为进一步的限定，所述步骤一中，对地震波形逐点求绝对值，确定出其中绝对值最大的振幅值，将其对应的反射时间设为零相位雷克子波的中心时间。

作为进一步的限定，所述步骤二中，以步骤一的最大绝对值振幅对应的反射时间为中心点，向其两侧寻找地震波形中的相邻极值点，确定雷克子波的主周期。

作为进一步的限定，所述步骤二中具体过程包括：

(1)若最大绝对值振幅为波峰，则寻找两侧相邻的波谷，此两个波谷之间的时间跨度为雷克子波的主周期；

(2)若最大绝对值振幅为波谷，则寻找两侧相邻的波峰，此两个波峰之间的时间跨度为雷克子波的主周期。

作为进一步的限定，所述步骤三中具体过程为对主周期取倒数，对主周期倒数值乘以第一因子，得到雷克子波的主频。

作为进一步的限定，所述步骤三中根据地震波形与雷克子波之间的相关系数计算雷克子波的幅度。