[发明专利]一种基于小波变换的新阈值函数地震资料去噪方法

权利要求书

1.一种基于小波变换的新阈值函数地震资料去噪方法，其特征在于所述的方法包括如下步骤：

步骤1：输入含噪信号道集f_i(t)；

步骤2：求取对应地震记录的方差σ_i(t)；

步骤3：求取加权叠加s(t)，其中σ_i(t)是方差，n表示所选自适应窗口从第n道开始，m为窗口长度，i表示第i道，t表示t时刻；

步骤4：迭代逐点求取自适应方差σ_i,new(t)，然后判断是否σ_i,new(t)-σ_i-1,new(t)<ξ，其中ξ为自己设定的一个极小的数；如果是，则转入步骤5，此时σ＝σ_i,new(t)；如果否，则用新的σ_i,new(t)替代原方差σ_i(t)，然后返回步骤3，其中σ_i,new(t)＝[f_i(t)-s(t)]²；

步骤5：将自适应方差代入求取阈值λ，其中λ是实数，N是信号的长度；

步骤6：对含噪信号道f_i(t)做小波变换，得到小波系数W_j,k，其中j为小波变换后的层数，k为第k个样点；

步骤7：建立新阈值函数：

其中λ为阈值，W_j,k为小波系数，e为实数；

步骤8：利用新建阈值函数对W_j,k进行阈值处理，得到估计系数

步骤9：利用所得分别进行小波重构后生成的地震数据；

步骤10：输出步骤9得到小波重构后生成的地震数据。