论文部分内容阅读
提高信噪比,是伴随整个地震勘探资料采集处理解释流程的一项重要工作。在高分辨率地震资料的处理工作中,信噪比的提高是根本和前提。而纵观各种已成熟或发展中的地震数据去噪方法,我们可以发现,大部分的地震去噪方法和技术都是建立在通用的数字技术基础上的,这些数字技术往往也可以用在类似的其他领域,如图像图形的去噪增强或压缩算法,语音及机械或其他振动信号的加强或分析等等,并且每种数字技术往往都有多方面的应用。先进的数字信号处理技术有利于在地震反射数据上进行复杂地质目标的成像和分析。经验模态分解(EMD)和Hilbert-Huang变换(HHT)是一个处理非线性非平稳信号的方法。经验模态分解先首先将信号分解为一系列的子信号,这些子信号在经过hilbert变换后会落在相同的时频区域中。HHT为随机信号或伪随机信号提供了一种时频域瞬时属性分析的方法。在本文中,EMD和HHT在地震数据上的应用研究主要包括:1)评价EMD和HHT在时间域及时频域确定有意义的地质信息的能力,2)在地震数据的f-x域内使用EMD分解来压制随机噪声,3)使用HHT的瞬时属性来设计一个好的滤波器,提高地震数据的信噪比。在滤波器设计和实际滤波过程中,HHT所表现出来的经验特性,可以比其他的滤波器更好的提高信噪比。本文使用模拟数据和实际地震数据证明了HHT和EMD可以成功应用在地震反射数据的处理上。在其中的一个实验中,强烈的工业交流电噪声可以通过HHT的陷频滤波进行消除。在另一实验中,使用HHT方法进行面波的压制,与其它方法相比,显示其具有较好的去噪效果。本篇文章的主要研究工作是如何更好的将EMD和HHT技术应用在地震反射数据上,尝试将EMD,HHT在其他领域中的新特点使用在地震勘探上,以提高反射地震数据的质量。这些研究不是要取代传统的地震数据处理方法,而是要将EMD和HHT补充到其中去。EMD是一种非常有效的时间域滤波器,HHT又是一种可以同时在时间和频率方面分析信号的优秀方法。同时,HHT的最大优点是它可以不通过任何先验信息就可以获取信号的瞬时振幅,瞬时频率和瞬时相位信息。先验信息通常是指采样率,噪声尖峰及步函数等时间序列的局部分布。本文首先介绍了EMD分解的有关基本概念及实施过程,以一个由多频率简谐波合成的调频信号为例来说明EMD分解的具体做法,相关概念,分解后的结果及原始信号的重构方法等。随后讨论了EMD分解应该注意的问题及EMD分解的不当会产生的问题,EMD分解的一个缺陷是容易产生频率混叠及出现原信号中不包含的虚假IMF分量,本文提出使用可变阈值的方法来改善这一缺陷,同时,阈值方法的改进也大幅加快了EMD分解的效率。为了讨论f-x域内的地震数据去噪方法,叙述了f-x反褶积的基本原理及发展过程,并讨论了其局限性,由此引出了在f-x域内使用EMD滤波方法来消除地震噪声的观点。使用模拟数据加入噪声前后的f-x域谱等频信息的对比来说明EMD滤波压制噪声的可行性。提出了f-x域EMD与小波阈值法的联合去噪的方法。使用简单模型,复杂模型的模拟数据来验证这种方法的有效性并与其他去除随机噪声的方法比较说明其优缺点。随后使用叠前及叠后的实际地震数据进行去噪实验,并对去噪结果进行详细分析。本文所提出的f-x域EMD与小波阈值法联合噪声衰减的方法是Battista B.M的f-x域EMD去噪方法进行了改进,通过分析EMD分解的过程指出直接删除IMF1分量的去噪方法是不严格的,而本文所提出的f-x域EMD与小波阈值法联合噪声衰减的方法在有效衰减噪声的同时可以更大程度的保留有效信号成分,整体效果优于f-x反褶积等其他去噪方法。本文系统的将HHT滤波技术引入到地震数据规则噪声的衰减过程中来。具体包括:基于HHT的工频噪声衰减,基于HHT的面波衰减方法研究和基于HHT的海上地震数据涌浪噪声衰减方法研究。首先分析了这几种规则噪声的特点,它们在EMD分解及Hilbert-Huang变换的过程中是易于分离的,从中得出结论,适合使用HHT滤波的方法来衰减这几种地震数据中的规则噪声。对模拟和实际地震数据使用所设计的方法进行去噪测试,结果表明基于HHT的方法与基于FFT的方法相比可以更好的压制噪声并保留有效信号。本文简要介绍了去除地震数据中工业交流电及面波干扰的常规方法,并讨论了这些方法的优缺点。随后,简要介绍了Hilbert-Huang变换的基本过程,并随即给出了基于HHT的时频域滤波的概念及实施过程。EMD和HHT对工频噪声具有突出的分离和时频分析的能力,因此可以设计一种使用HHT进行陷频滤波以消除地震数据中的工频噪声干扰的方法。对含工频噪声的地震数据进行EMD分解及HHT的计算,结果表明EMD分解可以将工频噪声较为完整的分离出来,进而在HHT谱上有很好的对应于工频噪声的时频显示。对单道含工频噪声的地震信号进行HHT陷频滤波去噪实验,结果表明该方法可以有效的压制噪声并保留有效信号。与其他基于FFT的滤波方法比较,该方法的去噪效果是最好的。随后对实际数据进行去噪实验,并与普通陷频滤波方法进行对比。本文进行了HHT时频滤波对面波干扰的压制研究。对合成含面波数据及实际数据都进行了HHT分析,所得时频谱中可以清楚的得到面波所对应的时频信息,对实际数据使用HHT滤波方法及FK滤波两种方法来压制噪声,对所得结果进行对比分析,本文方法优于传统的FK面波压制方法。而对于海上地震数据通常包含的低频强振幅的涌浪噪声,进行了HHT滤波方法压制研究。通过对实际数据的分析指出,直接使用EMD分解的方法来剔除涌浪噪声是不合适的,针对实际数据的涌浪噪声的特点,本文设计了HHT低切滤波的方法,有效的滤波除了海上数据的涌浪噪声,与FFT低切滤波方法的去噪效果进行对比分析,结果表明本文方法在有效压制噪声的同时更大程度的保留了有效反射信号。本文对基于BEMD的二维去噪技术及地震属性提取技术进行了研究,首先简要介绍了BEMD分解的基本过程,包括二维数据局部极值点求取方法,构造包络曲面的插值计算方法等。随后介绍了目前BEMD分解技术还存在的问题及改进方向。使用BEMD技术对二维双曲面数据进行滤波,与中值滤波,均值滤波,维纳滤波及小波阈值滤波方法比较,BEMD滤波对该数据的噪声消除效果更好。将BEMD技术引入到地震属性提取及相干数据体上断层信息的提取。BEMD对二维信号高低频分离的特性可以用于识别时间切片上的隐伏圈闭等构造。BEMD同时具有很好的边缘检测和纹理识别的能力,对相干数据的切片进行BEMD处理可以压制非构造的噪声信号,将其用于提高在地震相干数据切片上的断层识别能力,使断层信号更清晰连续,效果很好。对数谱在常规地震勘探数据处理中的使用方法主要是反褶积及子波提取,本文进行了基于对数谱的薄层估计研究。在数字信号处理的其他领域,如音频信号的处理,对数谱(倒谱)的定义方法将信号频谱的自然对数作为信号再求其频谱。本文分析了这种对数谱的特点,首先研究如何使用对数谱来快速测试地震仪正弦波信号的畸变。对数谱应用在地震仪信号畸变程度的检测当中时,可以更快速直观的判断正弦波信号的畸变程度。使用对数谱可以进行薄层厚度预测,所得的关于薄层厚度的分辨率要高于直接在时间域及频率域的厚度识别能力,本文使用楔状模型试验来检验对数谱在薄层识别上的效果,信号主频在30Hz时,薄层的识别能力可以达到10ms以内。本文提出了一种使用对数谱进行低频信号检测的方法并与其他谱分解方法进行比较,结果表明对数谱方法在对低频信号检测过程中,所取得的时间和频率分辨能力要优于其他的谱分解的方法。