地震道子波分解.pdfVIP

2010年第2期 内蒙古石油化工 59 地震道子波分解 邬 静 ，贺智娟，祁子元 (成都理工大学信息工程学院．四JII成都 610059) 摘 要：根据时频分析的思想，引入全新时频分析方法——子渡分解。子波分解可将地震道分解成 不同主频、不同时间的子波集合，从而实现地震道在时、频率域内的精确分解。分解后的子波可以按特定 频率和时间深度重构出多种地震道。采用树型模型对地震道实现子渡分解，全子波重构道与原始地震道 的一致性验证 了分解方法的正确性 。 关键词：时频分析；雷克子波；子渡分解；树 中图分类号 ：P631．443 文献标识码：A 文章编号：lOo6—7981(2010)02一OO59一O2 地震信号无论在时间上还是在空间上都是典型 1．2 树 的非平稳、非线性信号 ，只有从时间和频率域同时表 树既是一种递归结构。如图2所示即是一棵树 征信号，才能全面细微地反映信号特征。时频分析能 T．它是由根结点A和两棵子树T1和T2所组成，T1 较好地描述信号的频率随时间变化的规律。目前有 和T2分别位于A结点的左下部和右下部 ，A是B、C 很 多时频 分析方法 ，如 Hilbert变换、Hilbert— 的父结点，B、C为A的孩子结点，A结点位于第一 Huang变换 、短时傅立叶变换和小波变换等。但这 层 ．D、E、F、G结点位于第三层 ，树的深度为4。 些方法在时间、频率分辨率，以及对多频率成分信号 适应能力方面不能同时兼顾。本文根据地震道中含 有多种子波信息特点，通过子波分解对地震资料作 高分辨率时频分析。 1 理论基础 1．1 地震子波 地震道可以理解为或解释为单一地震子波b(t) 和地层反射系数序列R(t)的褶积： N 图2 树 的结构 X(t)一 R ×b 一R *b (1) n= 1 递归性和方向性是树的重要特性 ：递归保证 1r 此模型被广泛应用 ，但单一地震子波的褶积模 树的产生 ；方向保证树的结点按一定顺序产生，方便 型与实际情况有很大的距离。实际地震记录具有多 对其结点的操作 。 种子波属性 (如图1)，这是因为： 2 树型分解 1．1．1 实际地层往往是尺度各异的薄互层。 2．1 理论记录的合成 1．1．2 实际子波是变化的。 任意给定理论频率组合f~123：=：{0，77，77，98， 1．1．3 各道地震波的传播距离、入射角各不相同。 86，77，98，0，0，86，98，0}。将离散后的子波按fE12] 中的频率顺序 以dt一0．002s为相邻子波平移量进 行”移位叠加”，理论合成记录道X0．如图3所示 ： I 图1 地震记录子波属性 地震道的子波属性是子波分解 的理论依据 。若 能在解释前对地震道作子波分解 ，实时分频重构 出 能最大限度地突出目标体特征的剖面，将提高处理 解释精度。本文采用不同频率 的雷克子波 (Ricker 圈3 理论合成记 录 wavelets)，其数学表达式为： 2．2 子波分解 b(t)一 (1—2 f0t2)e (2)