脉冲反褶积研究及效果分析.doc

目 录 0 前言 1 1 脉冲反褶积的基本原理和实现方法 2 1.1 反褶积的基本原理 2 1.2 脉冲反褶积的原理 3 2 实际资料处理中脉冲反褶积的参数分析 7 2.1 反褶积因子长度的影响 8 2.2白噪系数大小的影响 11 3 结语 15 参考文献 16 致 谢 17 0 前言 地震勘探是进行油气探测的主要方法之一，如何获取高分辨率的地震剖面是地震勘探中一个重要的研究方向[4]。地震资料处理就是利用数字计算机对野外地震勘探所获得的原始资料进行加工、改造，以期得到高质量的、可靠的地震信息，为下一步资料解释提供直观的、可靠的依据和有关的地质信息[1]。 地震勘探数据处理包括许多技术方法[2-3]，例如动校正、叠加、数字滤波、反滤波、偏移成像、速度参数提取和分析等等。其中，反褶积是一种非常重要的技术，应用于岩性地震勘探和探测薄互油层及细地质结构时尤其重要[5]。反褶积可以压缩地震信号的脉冲宽度，分解复合波形，提高地震记录的纵向分辨率[9]。基于反褶积在处理地震资料中的作用，目前发展起来的反褶积方法很多[3]，但每种反褶积方法都受一定的假设条件制约，因此，每种反褶积方法只能在一定条件下使用[6]。在实际地震资料处理中，目前使用最多的反褶积方法有最小平方反褶积、预测反褶积、子波反褶积、同态反褶积[10]和最大（最小）熵反褶积[11]等。 最小平方反褶积是目前地震勘探中常用的反褶积方法，它旨在把地震记录中的地震子波压缩成为尖脉冲，从地震记录得到反射系数序列，或使地震记录接近反射系数序列[7]。最小平方反褶积的目的在于把已知的输入信号转换为与给定的期望输出信号在最小平方误差的意义下最佳接近的输出，脉冲反褶积则是期望输出为零延迟尖脉冲的最小平方反褶积。 本文依托PROMAX及MBP系统实现对地震资料原始炮集的处理，并就脉冲反褶积的特点及分辨率与信噪比的关系做一简单研究。 1 脉冲反褶积的基本原理和实现方法 1.1 反褶积的基本原理 因为大地滤波器的作用，地震脉冲变成了有一定持续时间的子波，从而使本来可以清晰反映地层层序的尖脉冲序列“模糊化”了，降低了地震记录对地层的纵向分辨能力。显然，要想利用地震记录划分岩层(主要是能区分薄层)，我们希望在所得到的地震记录上，每个界面的反射都表现为一个窄脉冲——其强弱与反射界面的反射系数大小成正比，而脉冲的极性则反映界面反射系数的符号[8]。如果这个希望能成为现实，我们就可以在常规地震勘探中划分薄层，而且可以利用地震资料研究岩性，为油气预测提供可靠的依据[2]。 这就是说，为了提高纵向分辨率，必须去掉大地滤波器的作用，把延续几十至l00ms的地震子波压缩成原来的震源脉冲形式，地震记录变为反映反射系数序列的窄脉冲组合，这就是反滤波所要完成的工作。由此可知，反褶积的目的就是为了把地震子波压缩成尖脉冲，使实际的地震记录变成反射系数序列。 假设地震记录为 （1—1） 其中为有效信号，为干扰波。 首先假设不存在干扰波，即： （1—2） 对两边求傅氏变换，则得到频率域的地震记录表示式： （1—3） 式中，、和分别为地震频谱、子波频谱和反射系数的频谱。 显然： （1—4） 如果令： （1—5） 则有： （1—6） 再对（1—6）式做反傅氏变换至时间域，就可得到： （1—7） 式中，为的时间函数。 根据（1—7）式知： （1—8） 因为为地震子波，而和之间又存在着频谱互为倒数的关系（即），所以把称为反子波，又叫做反褶积因子（deconvolution operator）。 由此可知，如已知地震子波，利用数学方法求出，再利用（1—7）式让反子波与地震记录做褶积，就可以求出反射系数序列，即 （1—9） 这样一个过程就叫做反褶积（或反滤波）。经过这样的处理，就可以达到把地震子波压缩成尖脉冲，从而达到提高地震记录纵向分辨能力的目的。 1.2 脉冲反褶积的原理 脉冲反褶积的基