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

摘要

本发明提供了一种基于小波变换的新阈值函数地震资料去噪方法，该方法针对油气地震勘探的低信噪比地区地震资料品质达不到处理、解释要求的问题，通过建立一种克服硬软阈值缺点的新阈值函数，以小波变换为核心思想，以达到地震资料去噪为目标，确定新的阈值函数。该方法与目前的同类技术相比较而言，具有去噪结果更加逼近原始有用信号，地震资料品质有较大改善的优点。

主权项

一种基于小波变换的新阈值函数地震资料去噪方法，其特征在于所述的方法包括如下步骤：步骤1：输入含噪信号道集fi(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)＝[fi(t)‑s(t)]2；步骤5：将自适应方差代入求取阈值λ，其中λ是实数，N是信号的长度；步骤6：对含噪信号道fi(t)做小波变换，得到小波系数Wj,k，其中j为小波变换后的层数，k为第k个样点；步骤7：建立新阈值函数：其中λ为阈值，Wj,k为小波系数，e为实数；步骤8：利用新建阈值函数对Wj,k进行阈值处理，得到估计系数步骤9：利用所得分别进行小波重构后生成的地震数据；步骤10：输出步骤9得到小波重构后生成的地震数据。