[发明专利]一种基于反比例模型的地震数据同时插值与去噪方法

摘要

本发明涉及一种基于反比例模型的地震数据同时插值与去噪方法，其步骤：输入初始不规则含噪声地震数据并进行正Seislet变换，得到其稀疏表示；采用阈值算子对Seislet系数进行阈值处理；对阈值处理后的Seislet系数进行反Seislet变换；根据反比例模型确定加权因子；选取POCS算法作为迭代算法，根据迭代公式和加权因子的取值，将反变换得到的时间域数据根据加权因子选择保留的数据，原观测数据按照加权因子的权重的比例回加到反Seislet变换后的数据中，得到本次迭代的重建结果；经过多次迭代后输出最终重建结果。本发明能够在不同的噪声强度下达到较好的重建效果并保证计算效率。

主权项

1.一种基于反比例模型的地震数据同时插值与去噪方法，其特征在于包括以下步骤：1)输入初始不规则含噪声地震数据，对输入的数据进行正Seislet变换，将输入数据变换到Seislet域后得到了稀疏的表示；2)通过阈值算子保留大于阈值的Seislet系数，将小于阈值的Seislet系数置零；3)对阈值处理后的Seislet系数进行反Seislet变换得到时间域的数据；4)根据反比例模型确定加权因子；所述步骤4)中，加权因子α的取值由如下反比例模型确定：其中，k为迭代次数，N为预先设定的最大迭代次数，q为控制因子，控制因子q用来决定加权因子的下降速率，所述控制因子q根据地震数据的噪声强度来确定；5)根据加权POCS算法的迭代公式和加权因子的取值，将步骤3)中得到的时间域数据根据确定的加权因子按照(I‑αR)的比例保留，原观测数据按照加权因子α的权重的比例回加到反Seislet变换后的数据中，得到本次迭代的重建结果；并判断k是否小于预先设定的最大迭代次数N，小于则返回步骤1)，再将重建结果作为下次迭代的输入，继续迭代，反之则输出最终重建结果，使缺失的数据被重建，噪声也得到压制；其中，I是单位矩阵；R是采样矩阵。