论文部分内容阅读
为了满足日益增长的能源需求,地震勘探已经从东部转移到西北部,从平原转移到沙漠,在这些地区采集到的地震资料中包含着很强的线性干扰和面波,其中面波的存在使得地震数据的信噪比急剧降低,有效性也遭到了极大的破坏。因此,必须对叠前资料中的面波干扰进行压制,而常用的面波处理方法如f-x滤波,f-k滤波,小波变换等,已不能满足这种需要。论文首先介绍了面波特点及多尺度分析的基本理论和性质,接着利用面波与有效波在频谱上的差异,分别采用多尺度变换中的二维小波变换和Curvelet变换,将地震记录分解到不同的频带内,只在面波的优势频段内进行压制处理,以保证其他频段的有效信号不会受到损伤。主要研究内容如下:首先,对地震记录做二维小波变换,在频率波数域中确定面波的优势频段;然后,对这些区域的小波系数进行阈值处理,本文将面波压制中的阈值选取与图像处理里的阂值构造方法联系起来,引入了均值及直方图两种图像阈值方法,并将其与常用的两种阈值函数分别应用到理论模型及实际地震数据的面波压制中,对比处理后的结果发现:地震记录经直方图阈值方法处理之后,其信噪比及分辨率都得到了明显提升。最后,将本文提出的基于小波变换及直方图阈值的面波压制算法与常用的面波处理方法分别应用到实际地震资料处理中,实验结果表明:本文提出的算法在压制面波的同时有效信号也得到了增强,处理效果要优于常用处理方法。接着,采用Curvelet变换对地震数据进行处理,它是在小波分析基础上发展起来的且能够最优表示具有曲线奇异的二维信号。算法具体思路为:首先对地震资料进行Curvelet变换,在Curvelet域中确定面波的优势频带;其次根据面波视速度低的特点,选择其在Curvelet域内对应的方向系数矩阵,利用本文建立的时域角度与Curvelet域方向的对应关系,只需根据面波在地震记录中的角度范围,就可确定其在Curvelet域的系数矩阵;接着,利用本文改进的加权均值阈值方法及其他三种阈值,对已确定存在面波的方向系数矩阵进行迭代阈值处理。理论模型及实际地震记录均证明了改进阈值的有效性。最后,将本文提出的两种分频算法分别对实际地震资料进行面波压制,通过对处理结果分析发现:两种算法都能很好的压制面波,且Curvelet变换分频算法在有效波的连续性保持上要强于小波分频方法。