地震子波提取研究.pdfVIP

地震子波提取研究 王 辉 张苗苗 胡鹏飞 (长安大学地质工程与测绘学院 陕西 ·西安 710054) 摘 要 地震子波的提取是层速度反演的基础和关键，而层速度又是地震勘探中一个极其重要的参数，它在研究地 层的压力变化、预测高压带、圈定砂体范围、静校正、确定岩性和储层描述中起着重要作用。地震子波的提取是数据分 析和处理的关键，采用不同的子波提取计算方法，提取出的子波有所差异。 关键词 地震子波 层速度 反演 子波提取 中图分类号：P315 文献标识码：A 1引言 ·o¨ 4“ o00 口12 3E．007 目前提取子波的方法主要有直接观测法 (只适用于海上 2E 007 勘探)、自相关法 (WoldKolmogorov法)、Claerbout算法、多项 1巳007 式求根法(无法处理单位圆上的根)、同态滤波法(抗噪声能力 O 弱，但没有苛刻的假设条件，只要求子波的对数谱和反射系数 1￡-007 的对数谱序列在时间上可分)。由于勘探早期往往没有满足 0 20 40 60 80 要求的井资料，研究直接从地震记录中提取子波的方法显得 图 1：合成地震记录 图2：Claerbout法提取子波 非常重要。本试验提取子波采用了Claerbout提出的算法进 从上述测试效果可以看出，Claerbout法提取的子波与原 行求解。 始雷克子波很相似，但经过测试，若窗口函数长度选取为不同 2Claerbout基本算法原理 的采样点，则提取的子波则会有明显的差异。因此实际中，对 由信号理论知道，若已知某一信号的振幅谱，可以求出该 于窗函数窗口大小的选取需要试算，这无疑增加了工作量，可 振幅谱对应的唯一的最小相位信号，且最小相位信号的相位 以设法进一步研究自适应的窗口选取办法。 谱，可以通过Hierbout变换唯一的求出，即： 4结语 9(w)=-HT(1n(IB(w)I))=-HT(1n(Ix(w)})) 反褶积的结果好坏与子波的提取有很大关系，Claerbout 由于直接使用Hierbout算子在零点处是奇异的，没有定 算法得出的子波与原始子波较为吻合，实际应用中可以考虑 义的，可以先变换到时间域，再反变换回来。但Claerbout(1992) 用这种方法来提取子波，但提取过程中对窗口的选择直接决 提出了一种可以不直接求Hierbout变换，但理论上与Hierbout 定了计算工作量的大小和计算精度。在提取地震子波之前， 变换等价的可以提取最小相位子波的方法，事实上也可以作 应首先估计地震子波的有效长度，而不能随便用一个长度参 为同态法滤波的特例，其算法流程如下： 数进行子波提取。实际计算提取的反子波和原始子波进行相 (1)计算对数，B(w)=ln(F1(w)J，F(w)为原始信号(地震记 关计算，其结果并非理想状态下的一个脉冲信号，表明从地震 录)的频谱。 记录中提取的子波不一定是最小相位的子波，这是零延迟脉 (2)计算逆傅里叶变换，B(w)一b(t) 冲信号反褶积的基础。 (3)t0时lb(t)乘0，tO时，b(t)乘2。但实际中，离散信号 的处往往只截取右端有限项处理，其具体处理办法可以为 参考文献 b(n)=j 2．b(n)0《％(