一种地震数据重建方法

摘要

本发明提供一种地震数据重建方法，获取地震资料后，对其进行傅里叶变换，在时间-频率域对地震数据进行奇异值分解和截断处理，通过在频率域对Hankel矩阵化的地震资料奇异值分解并进行截断处理，可以在对有效信号伤害较小的情况下压制随机噪声并对缺道的地震数据进行插值，插值后的信号与原信号的振幅差会随着迭代不断减小，应用了分步插值的方法改善地震数据的相干性，本发明在时间-频率域通过奇异值分解进行去噪和插值处理，克服了常规方法不能同时进行去噪和插值的缺陷，并能够有效地减少去噪时有效信号的损失，是一种相对保幅的数据重建方法。

主权项

一种地震数据重建方法，其特征在于，该地震数据重建方法包括以下步骤：(1)将波场变换至频率域；(2)对每个频率的波场分别构建相应的分块矩阵；(3)对每个频率分量的分块矩阵进行奇异值分解；(4)秩减得到每个频率分量对应分块矩阵的近似矩阵；(5)对秩减近似矩阵反对角线求平均值获取频率域序列；(6)对频率域序列计算得到处理过后的地震数据；(7)重复步骤2到步骤6进行迭代处理得到满足要求的频率域数据；(8)将频率域数据变换到时间域。