李 稳,刘伊克,陈 颙,徐善辉,谭雅丽,何银娟,花鑫升,贾 佳

(1.中国地震局地球物理勘探中心,河南郑州450002;2.南京大学地球科学与工程学院,江苏南京210023;3.中国科学院地质与地球物理研究所,北京100029;4.中国地震局地球物理研究所,北京100081)

我国地震勘探工作越来越多地涉及到复杂地表问题。一般而言,复杂地表地区是指包含沙漠、戈壁、山地、黄土塬、水网等地貌,以及风化层速度、厚度、地表高程等剧烈变化的地区[1]。在复杂地表地区,既可能存在高速地层、岩层直接出露,又可能存在低速、疏松的冲洪积物、风积物堆积。对于地震勘探而言复杂地表问题的影响主要表现在近地表地震波速度剧烈变化,强散射波、面波等干扰波发育,引起地震波场、反射波同相轴畸变等[2-3]。

在我国西部某复杂地表工区,受崎岖地表和近地表多尺度非均匀体等因素的影响,地震勘探采集的地震数据散射干扰波发育,信噪比较低。地震散射波是当散射障碍体的尺度与入射波波长相比很小或相近时,由于干涉效应而产生的与散射障碍体几何形状有关的相干或不相干的波列,且往往以众多散射波互相叠加的复杂形态出现[4]。散射干扰波在时间分布、频率特征等方面均与有效反射波组接近,采用傅氏变换、小波变换[5-6]等传统的时频域滤波方法难以达到理想的去噪效果。

目前公认的好的去噪方法是分析噪声与有效地震信号在不同“域”中的表现形式和特征差异,进而采取最有针对性的去噪方法[7]。……