地震子波的提取及其在最小平方反褶积中的应用.pdfVIP

地震子波的提取及其在 最小平方反褶积中的应用 王 超 (长安大学地测学院 陕西 ·西安 710064) 摘 要 本文通过希尔伯特变换法从实际地震资料中提取最小相位子波，将最小相位子波与通过求解托普利兹方程 组得到的反滤波因子做褶积，实现了对地震资料的最小平方反滤波，提高了地震剖面的分辨率。 关键词 希尔伯特变换 地震记录 最小相位子波 托普利兹方程组 中图分类号：TN911．7 文献标识码：A 为了精细地解释地下地质情况，在实际地震资料处理过 相位子波最关键的一步是选择地震数据质量好的时间段作为 程中，常常会对偏移(或叠后)数据体进行修饰性处理，从而提 提取时窗，利用傅里叶正变换将地震记录由时间域一一空间 高地震剖面的分辨率、改善其波组特征。对地震资料做最小 域变换到频率域一一波数域，求取地震记录的振幅谱，然后求 平方反滤波就是提高地震剖面分辨率的一种常用方法。该方 出对数谱，用希尔伯特变换法求出最小相位子波的振幅谱，最 法主要分为两个步骤：第一步为地震子波的提取，第二步为最 后通过反傅里叶变换即可求出所需要的最小相位子波。通过 小平方反滤波。准确提取地震子波是进行反褶积的前提，而 FORTRAN语言编程，输出最小相位子波序列，经过MATLAB 反褶积是提高地震资料分辨率的重要手段。地震子波是一段 成图，通过分析可得出所求出的子波为符合要求的最小相位 具有确定的起始时间和有限能量，有限连续长度的信号，它是 子波 。 地震记录中地震的基本单元；一般认为，地震震源激发时所产 1．2最小平方反褶积 生的地震波仅是一个连续时间很短的尖脉冲，随着尖脉冲在 记最小相位子波为序列为，反滤波因子为，最小平方反滤 介质中的传播，尖脉冲的高频成分很快衰减，波形随之增长，变 波就是通过将与做褶积运算，将地震子波压缩成一个窄带尖 成了一个具有有限频带宽度和一定延续时间的地震子波。实 脉冲。对地震记录做最小平方反滤波运算分为以下几步：(1) 际的地震子波是一个很复杂的问题，子波是震源波形经过地 求取最小相位子波的自相关序列；(2)构造托普利兹方程组； 层滤波和检波器耦合与仪器效应等综合因素共同作用的结果， (3)用Levmson递推算法求解Toepli~方程组，求出反滤波因 具体因素的作用会有所不同，对具体勘探 目标的影响也不相 子；(4)最小相位子波序列与滤波因子做褶积，压缩地震子波。 同。在地震子波的选取过程中，由于其能量和相位特征较难 在最小平方反褶积程序的编写过程中要注意以下几个问题： 确定，实际地震资料通常很难满足子波提取的一些假设，使得 (1)选择合适的反滤波因子长度，虽然理论上反滤波因子的长 地震子波提取技术成为制约地震资料处理质量进一步提高的 度选择是任意的，但是通常情况下，一般选择的反滤波因子长 关键因素。由于本文的最终 目的是对地震资料做最小平方反 度与最小相位子波的长度接近。(2)选择合适的自相关时窗 滤波，因此本文选择提取的子波是最小相位子波，采用的是希 长度，经过反复试算发现最佳时窗长度应该大于反滤波因子 尔伯特变换法。反褶积作为地震数据处理中的一个基本的处 长度的两倍。 理环节和提高地震分辨率的一种重要手段，主要是通过压缩 本文所采用的编程语言为FORTRAN，本文所采用的思路 地震记录中的基本地震子波 以及交混回响和周期多次波的方 是：首先给定一个理论上的一定长度的最小相位子波，通过程 式，提取反射系数序列；最小平方反滤波，其又称为脉冲反滤 序试算发现通过给定的最小相位子波与求出的反滤波因子做 波，是一种比较简单但是很实用的滤波方法；最小平方反滤波 褶积后确实变成了一个窄带的尖脉冲，达到了压缩子波的目 就是将期望输出设计为尖脉冲，构造托普利兹矩阵方程组，通 的，说明了本程序编写的正确性；然后通过对程序的进一步改 过解方程最终求解出反滤波因子，然后用求出的反滤波因子 进，实现了将从实际地震资料提取的最小相位子波与其反滤 与己求出的地震子波做褶积，从而达到对地震记录做最小平 波因子做褶积，通过剖面显示发现滤波效果很好，达到了提高 方反滤波的 目的。 地震资料分辨率的 目的。