[发明专利]一种地震数据去噪方法、装置、电子设备及存储介质

说明书

本申请提供一种地震数据去噪方法、装置、电子设备及存储介质，方法包括：获得地震观测数据；根据所述地震观测数据构造稀疏反演模型；构造正则参数序列，所述正则参数序列中正则参数为正数，且从首个正则参数起逐渐变小；根据所述正则参数序列，采用迭代阈值法求解所述稀疏反演模型，获得正则参数对应的解，所述正则参数对应的解为去噪后的地震数据；将正则参数以及所述正则参数对应的解代入评价函数，确定评价函数的极值；确定用于计算获得所述极值的所述正则参数对应的解为最终的去噪后的地震数据。由此，在对不同的地震数据进行去噪时，无需人工选择阈值参数，同时，也提高了地震数据去噪的工作效率。

技术领域

本申请涉及地震数据处理领域，具体而言，涉及一种地震数据去噪方法、装置、电子设备及存储介质。

背景技术

地震数据在采集、处理等过程中会不可避免的产生各种噪声，噪声的存在会降低地震数据的信噪比，使得偏移成像、全波形反演等关键处理的结果受到影响，因此地震数据去噪技术是地震数据处理的重要手段。一直以来研究者对地震数据去噪的研究热情不断，并且提出了各种去噪方法。

传统的地震数据去噪方法由于地震数据在变换域能够和随机噪声分离开，其中地震数据在变换域有稀疏的表达和较大的系数，而噪声在变换域中的表现为稠密分布，系数振幅较小，通过阈值运算可以消除随机噪声，然后通过反变换到时间空间域得到去噪的结果。但传统的地震数据去噪方法需要人为的调节阈值参数，才能针对不同的地震数据获得可靠的去噪效果。

发明内容

有鉴于此，本申请实施例的目的在于提供一种地震数据去噪方法、装置、电子设备及存储介质，以改善传统的地震数据去噪方法需要人为针对不同的地震数据调节阈值参数的问题。

第一方面，本申请实施例提供了一种地震数据去噪方法，包括：获得地震观测数据；根据所述地震观测数据构造稀疏反演模型；构造正则参数序列，所述正则参数序列中正则参数为正数，且从首个正则参数起逐渐变小；根据所述正则参数序列，采用迭代阈值法求解所述稀疏反演模型，获得正则参数对应的解，所述正则参数对应的解为去噪后的地震数据；将正则参数以及所述正则参数对应的解代入评价函数，确定评价函数的极值；确定用于计算获得所述极值的所述正则参数对应的解为最终的去噪后的地震数据。

本申请实施例通过设置评价函数来对多个正则参数以及对应的稀疏反演模型的解进行评价，通过确定评价函数的极值来确定最终的去噪后的地震数据。由此，在对不同的地震数据进行去噪时，无需人工选择阈值参数，同时，也提高了地震数据去噪的工作效率。

进一步地，所述评价函数表征所述正则参数对应的解与所述地震观测数据的差的范数，与所述正则参数之间的差异。

本申请实施例通过设置评价函数来对多个稀疏反演模型的解与地震观测数据的差值以及对应正则参数的差异进行评价，使得通过评价函数可以更加准确的确定最终的去噪后的地震数据。

进一步地，所述评价函数为：

其中，P_i为与第i个正则参数对应的评价函数的值；λ_i为第i个正则参数；d_i为与第i个正则参数对应的解；d_obs为所述地震观测数据；所述确定评价函数的极值，包括：确定评价函数的极大值。

本申请实施例通过设置评价函数来对多个正则参数以及对应的稀疏反演模型的解进行评价，由此，通过确定评价函数的极大值来确定最终的去噪后的地震数据。使得在对不同的地震数据进行去噪时，无需人工选择阈值参数，同时，也提高了地震数据去噪的工作效率。

进一步地，所述构造正则参数序列，包括：对所述地震观测数据进行稀疏化，根据稀疏化的地震观测数据确定所述正则参数序列中首个正则参数的值。