一种多尺度地震资料随机噪音衰减方法

摘要

本发明提供了一种多尺度地震资料随机噪音衰减方法，属于地震信号处理领域。本方法包括：(1)对包含随机噪音的地震资料进行Contourlet变换，将地震数据变换到Contourlet域，获得该地震数据在Contourlet域的一个低频子带，即低频系数，以及多个方向子带，即高频系数c_j，k；(2)设定当前大小为M×N的方向子带Contourlet高频系数的基本阈值；(3)计算当前大小为M×N的方向子带Contourlet高频系数的调整因子；(4)计算当前大小为M×N的方向子带Contourlet高频系数的自适应阈值λ_j，k＝μ_j，kλ；(5)对步骤(1)获得的高频系数进行Contourlet高频系数运算，得到。

主权项

一种多尺度地震资料随机噪音衰减方法，其特征在于：所述方法包括：(1)对包含随机噪音的地震资料进行Contourlet变换，将地震数据变换到Contourlet域，获得该地震数据在Contourlet域的一个低频子带，即低频系数，以及多个方向子带，即高频系数c_j，k；(2)设定当前大小为M×N的方向子带Contourlet高频系数的基本阈值$\lambda =\sigma \sqrt{2\log \mathrm{MN}};$(3)计算当前大小为M×N的方向子带Contourlet高频系数的调整因子${\mu }_{j,k}=\sqrt{c_{j,k}^2/S_{j,k}^2};$(4)计算当前大小为M×N的方向子带Contourlet高频系数的自适应阈值λ_j，k＝μ_j，kλ；(5)对步骤(1)获得的高频系数进行Contourlet高频系数运算，得到(6)将邻域窗口按照步骤(5)得到系数逐个移动，每移动一次，重复第(2)‑(5)步，直到所有Contourlet高频系数处理完成；(7)利用步骤(1)得到的Contourlet低频系数和运算后的Contourlet高频系数进行Contourlet反变换，得到随机噪音压制后的地震资料；(8)输出随机噪音压制后的地震资料。