[发明专利]基于Shearlet变换的地震信号去噪方法

摘要

本发明公开了一种基于Shearlet变换的地震信号去噪方法，其包括如下步骤：①读入二维地震剖面数据；②将二维地震剖面数据扩展为一长宽为奇数的方阵；③构造频域方向滤波器组；④将各个变换矩阵分别与信号向量相乘，再分别做二维傅立叶反变换，得到各个方向和尺度的Shearlet系数；⑤阈值处理；⑥将经过阈值处理的shearlet变换系数做shearlet反变换获得去噪之后的信号。本发明将含有噪声的地震信号做拉普拉斯分解，然后利用剪切波函数进行滤波处理，得到对应的剪切波系数，通过阈值处理滤除噪声信号，在通过逆非下采样剪切波变换恢复去噪的信号，获得较好的去噪效果，具有很好的实用价值。

主权项

1.基于Shearlet变换的地震信号去噪方法，其特征在于，包括如下步骤：步骤一：读入一组二维地震剖面数据，把数据存为矩阵S，读出该矩阵的大小为[m n]，同时设定Shearlet的分解级数k以及Shearlet分解所采用的Meyer小波；步骤二：将二维地震剖面数据S依照下式扩展为一长宽为奇数的方阵S₁：S₁=MS其中，M为一变换矩阵，大小为：a=max(m,n)若max为奇数a=max(m,n)+1若max为偶数；步骤三：构造频域方向滤波器组；步骤四：将各个变换矩阵分别与S₁信号向量相乘，再分别做二维傅立叶反变换，得到各个方向和尺度的Shearlet系数c_i,j。步骤五：阈值处理；根据每个尺度下各个方向子带的稀疏的平方差与该尺度改方向子带下的噪声方差的比值来确定该子带的阈值，公式如下：τi,j=σϵi,j2/σi,j,n2]]>其中表示第j尺度下的第i个方向子带的第n个系数的方差、表示j尺度下i方向的噪声方差，通过比值确定各个子带的阈值，对信号变换得到的shearlet系数做阈值处理，获得所有尺度下的变换稀疏；步骤六：将经过阈值处理的shearlet变换系数做shearlet反变换获得去噪之后的信号。