子波基本理论与提取方法

1、子波基本理论与提取方法1地震子波基本原理由震源激发、经地下传播并被人们在地面或井中接收到的地震波通常是一个 短的脉冲振动，称该振动为振动子波。它可以理解为有确定起始时间和有限能量， 在很短时间内衰减的一个信号。地震子波其振动的一个根本属性是振动的非周期 性。因此，它的动力学参数应有别于描述周期振动的振幅、频率、相位等参数， 而用振幅谱、相位谱等概念来描述。子波一般是物理可实现的，特别是地震子波，作为一个物理滤波器的响应函 数，自然是物理可实现的，所有必定为非零相子波，但不同子波相位延迟不同。 子波包括最小相位子波、最大相位子波、混合相位子波。子波的Z变换是一个多项式：B(z) = b + b

2、z + b z2 +. + b zn若此多项式的全部零点均在单位圆外，则为最小相位子波；在单位圆内，为 最大相位子波；零点在单位圆的内外都有，则为混合相位子波。2地震子波的数学模型实际中的地震子波是一个很复杂的问题，因为地震子波与地层岩石性质有 关，地层岩石性质本身就是一个复杂体。为了研究方便，仍需要对地震子波进行 模拟，目前普遍认为雷克提出的地震子波数学模型具有广泛的代表性，即称雷克 子波。最小相位的地震子波的数学模型为：b(t) = e - at2 sin 2兀 ft式中：f为子波的主频；口二2 f 2ln(M)为子波衰减系数；M =|mi/m2|为最 大波峰值%和最大波谷值m2之对比。其

3、波形大致如图所示：3地震子波提取的基本方法地震子波的提取方法有两大类：第一类是确定性子波提取方法；第二类是统 计性子波提取方法。确定性子波提取方法指的是利用测井资料首先计算出反射系 数序列，然后结合井旁地震道由褶积理论求出地震子波，它的优点是不需要对反 射系数序列的分布作任何假设，能得到较为准确的子波，而统计性方法的优点是 不需要测井信息也可以得到子波的估计，但缺点是需对所用的地震资料和地下的 反射系数序列的分布进行某种假设，所得子波理论上的精度不是高很。正是由于这个原因，所以本文将分别最这两种方法进行C语言的编程实现， 来对统计性子波提取方法和确定性子波进行地震信号反褶积的效果作一下对比 和

4、演示。4统计性子波提取方法原理虽然子波一般是未知的，但地震记录中包含有子波，因此，可以从地震记录 中求取子波。虽然目前求取地震子波的方法也很多，但下面还是介绍一下统计求 取子波方法。若将子波作为一般信号对待，则子波也可用s表示，假设反射系数是随机 的白噪序列，则有地震记录)的自相关和子波)的自相关相等，于是有记录 的振幅谱1 X)|和子波的振幅谱1 S(3)I相等I S ()IT X(3 )1及其对数也相等lnI S(3)I=lnI X(3)I理论已证明，当子波为最小相位时，其对数谱序列(或称复赛谱)S(n)是因果 序列1S (n) = 一j ln I S(3)Iej 兀33兀一兀由于lnI

5、S (3 )I为实偶函数，因此S (n)是实的因果序列。任何实序列都可写成奇部和偶部序列之和，故S(n)可写成：S (n) = S o(n) + S e(n)即子波对数谱序列S(n)的奇部S o(n)和偶部Se(n)有下述两个性质：第一，由于S(n)的因果性，起奇部和偶部有如下关系 TOC o 1-5 h z HYPERLINK l bookmark30 o Current Document S o(n) = sign(n) S e(n)-S e(n) = sign(n) S o(n) + S o(n)6 (n) /式中，1n 00n =0-1n 0第二，S(n)的偶部和奇部的傅里叶变换为其傅

6、里叶变换的实部和虚部。设S(n)的傅里叶变换为Sr(f) = ReSi(f)，S i(f) = ImS(f)，S(f)为子波的对数谱， 则S(f) = S r(f) + iS i(f)由傅里叶变换性质有F S e(n) = S e( f)、F S o(n) = iS o(f)故有S(f) = S r(f) + iS i(f)=S e(f) + iS o(f)即S e( f) = S r(f)-S o(f) = S i(f)于是求子波的方法可归结为：1)用多道统计的方法获得可靠的子波对数谱的实部。由子波谱S (o) =1 S (o) I e网则有S ()= lnl S (w)l= lnl S ()l +却()由若干道振幅谱的几何平均(或多道记录的相关函数平均)确定子波振幅谱 的对数谱lnl S (o )1。2)由子波振幅谱对数求子波相位谱l中(o)。计算公式为lnl S (o )l 马 S e(n)S o(n) = sign(n) S e(n)S o(n)鸟中(o)3)计算子波S。由S (o)及中(o)得S (o) =l S (o)l e mod S (t)由于干扰的影响和反射