[发明专利]一种自适应小波阈值求取方法

摘要

本发明提供了一种自适应小波阈值求取方法，属于地震勘探资料处理等数字信号处理领域。本方法包括(1)输入含噪信号道集fi(t)；(2)求取自适应加权叠加(3)迭代逐点求取自适应方差σi，new(t)，然后判断是否σi，new(t)‑σi‑1，new(t)＜ξ，如果是，则转入步骤(4)，如果否，则用新的σi，new(t)替代原方差，然后返同步骤(2)；(4)将自适应方差代入求取阈值公式得到阈值λ；(5)利用阈值λ对wj，k进行小波域去噪；(6)将去噪后的进行小波重构得到小波重构后生成的地震数据；(7)输出步骤(6)得到小波重构后生成的地震数据。通过理论模型及实际地震数据测试，本方法去噪效果明显，具有较强的针对性。

主权项

一种自适应小波阈值求取方法，其特征在于：所述方法包括：(1)输入含噪信号道集fi(t)；(2)求取自适应加权叠加(3)迭代逐点求取自适应方差σi,new(t)，然后判断是否σi,new(t)‑σi‑1,new(t)＜ξ，如果是，则转入步骤(4)，如果否，则用新的σi,new(t)替代原方差，然后返回步骤(2)；(4)将自适应方差代入求取阈值公式得到阈值λ；(5)利用阈值λ对wj,k进行小波域去噪；(6)将去噪后的进行小波重构得到小波重构后生成的地震数据；(7)输出步骤(6)得到小波重构后生成的地震数据，所述步骤(2)是使用下面的自适应加权叠加公式得到自适应加权叠加的：s‾(t)=1Σi=nn+m1σi(t)×Σi=nn+m1σi(t)fi(t)]]>其中，fi(t)为原始地震记录，σi(t)为方差，n表示所选自适应窗口从第n道开始，m为窗口长度，i表示第i道，t表示t时刻。