[发明专利]一种地震随机噪声衰减方法

权利要求书

1.一种地震随机噪声衰减方法，其特征在于：所述方法包括：

(1)对含噪数据进行行和列方向的平移得到新的二维数组；

(2)对步骤(1)得到的二维数组进行Curvelet变换得到曲波系数；

(3)对步骤(2)得到的曲波系数进行阈值去噪得到新的曲波系数；

(4)对步骤(3)得到的新的曲波系数进行Curvelet逆变换得到新的二维数组；

(5)对步骤(4)得到的新的二维数组进行行和列方向的逆向平移，恢复数据原先的状态，即得到一次平移的去噪结果；

(6)判断是否达到设定的循环处理次数，如果达到，则处理完成，转入步骤(7)；如果否，则返回步骤(1)；

(7)对所有去噪结果进行线性平均得到最终去噪结果。

2.根据权利要求1所述的地震随机噪声衰减方法，其特征在于：所述步骤(1)是这样实现的：

采用循环平移方法，即通过一个平移函数对由含噪数据组成的二维矩阵沿行和列方向进行平移，形成新的二维数组。

3.根据权利要求2所述的地震随机噪声衰减方法，其特征在于：所述步骤(5)中所述一次平移的去噪结果如下：

s^i,j=S-i-j(C-1(Λ[C(Si,j(x))])]]>

其中，S为循环平移算子，下标i，j和-i，-j分别为行和列方向上的平移量，C为Curvelet变换算子，C^-1为逆Curvelet变换算子，Λ为阈值算子。

4.根据权利要求3所述的地震随机噪声衰减方法，其特征在于：所述步骤(6)中的所述循环处理次数为1到含噪数据二维数组的行数×列数之间，即最小平移次数为1，最大平移次数为行数×列数。

5.根据权利要求4所述的地震随机噪声衰减方法，其特征在于：所述步骤(7)是这样实现的：

利用下面的线性平均公式求得最终去噪结果：

s^=1K1K2Σi=0,j=0K1,K2si,j,]]>

其中，K₁，K₂分别表示行和列方向上的最大平移量。