基于离散小波变换的地震资料自适应高频噪声压制-物探与化探.pdfVIP

　第３７卷第１期 物　探　与　化　探 Ｖｏｌ．３７，Ｎｏ．１　 　２０１３年２月 ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ Ｆｅｂ．，２０１３　 基于离散小波变换的地震资料自适应高频噪声压制 １ １ ２ １ 李雪英 ，张晶 ，孔祥琦 ，侯相辉 （１．东北石油大学地球科学学院，黑龙江 大庆　１６３３１８；２．大庆钻探工程公司地球物理勘探一公 司，黑龙江 大庆　１６３３５７） 摘要：针对现有的小波变换阈值去噪方法无法适用于地震资料高频噪声压制的缺点，笔者分别对阈值函数和阈值 选取方案进行了改进，提出了连续硬阈值函数与自适应阈值相结合的地震资料高频噪声压制方法。连续硬阈值函 数兼具软、硬阈值函数的优点，可提高重构地震信号的保真度，减少人为噪声误差；自适应阈值方案可根据非平稳 地震数据中的能量时变、空变分布特点，通过引入不同子带小波系数标准差、代数平均值以及几何平均值等统计参 数，使阈值能够随不同子带的小波系数能量变化而自动调整，以适应地震资料高频噪声压制的要求。实际地震数 据处理结果表明，笔者提出的方法在提高信噪比的同时，可保护陡倾角反射界面信号，提高噪声压制后地震数据的 保真度。 关键词：离散小波变换；地震资料；高频噪声压制；阈值函数；自适应阈值 中图分类号：Ｐ６３１．４　　　文献标识码：Ａ　　　文章编号：１０００－８９１８（２０１３）０１－０１６５－０６ 　　基于离散小波变换阈值去噪方法由于具有算法 值函数虽然保证了阈值处连续性，但改变了小波系 简洁、计算速度快、适应性广、去噪效果显著等特点 数大小，整体消减了小波能量；改进后的非线性阈值 成为图像处理领域应用最为广泛的一种信号降噪处 函数可保证阈值处连续性，并随着小波系数增大而 理技术［１－５］。此类方法可以得到信号的最优近似估 逐渐趋近于硬阈值函数。它介于软、硬阈值函数之 计，主要适用于含有高斯白噪声且噪声能量不是很 间，是一种折衷策略，其阈值附近的小波系数也有较 高的二维信号的降噪处理。其中阈值函数的选取和 大的改变。与软阈值函数一样，非线性阈值函数也 自适应阈值的确定是该类方法中的技术关键，二者 将会导致重构信号缺乏必要的保真度。这就需要寻 选择正确与否直接关系到高频噪声的压制效果。然 找一种全新的阈值函数，使其能够完全继承硬阈值 而，现有的阈值函数和阈值确定方法并不能够完全 函数的保真特性，同时在阈值处又具有软阈值函数 适用于地震资料的降噪处理分析。因此，有必要针 连续特性。 对地震数据噪声分布特性提出一种新的阈值方案以 第二，地震资料是一种时变、空变的非平稳数 适应地震数据处理的要求。 据，随着传播时间的增加，地震波的能量逐渐被吸 收，导致深层地震子波的能量与浅层相比逐渐减弱， １　研究意义 远偏移距与近偏移距相比，能量也逐渐减弱。根据 第一，地震资料高频噪声压制只是地震数据处 巴塞瓦定理，小波系数的平方具有能量的量纲，因 理流程中的一个环节，其处理准则是在提高信噪比 此，小波系数分布也具有上述时变、空变规律。实际 的同时保持有效信号的保真度，从而为后续处理提 地震数据研究表明，不同分解层次的小波系数能量 供必要的前提与基础，这与图像处理主要为了改善 随着分解层次的增加是逐渐增大的，且 ＨＨ、ＬＨ、ＨＬ 视觉效果的处理准则有着本质的区别。现有的小波 三个子带在不同分解层次的能量分布规律也不尽相