犚 犲 狊 犲 犪 狉 犮 犺狆 狉 狅 狉 犲 狊 狊狅 犳 狊 犲 犻 狊 犿 犻 犮狑 犪 狏 犲 犾 犲 狋 犲 狓 狋 狉 犪 犮 狋 犻 狅 狀 犵
１２ ２ １ ２ ， ， ， ＧＡＯＳ ｈ ａ ｏ ｗ ｕ ＨＡＯＢ ｏ ｈ ｅ ｎ ｈ ｕ ａ ｕ ｎ ｉ ｎ Ｚ ＨＥＺ ＭＡＯ Ｙ ｇ ，
收稿日期 ２ ０ ０ ８ １ ２ ２ ０； ０ ０ ９ ０ ３ ２ ２． 修回日期 ２ 作者简介 高少武 ， 陕西周至人 ， 高级工程师 ， 博士研究生 ， 发表论文 ２ 申请发明专利 ７ 项 ， 现在东方地球物理公司从 事 物 １ ９ ６ ５ 年生 ， ０ 余篇 ， 探方法研究与软件开发工作 ． （ ： ） Ｅ ｍ ａ ｉ ｌ ａ ｏ ｓ ｈ ａ ｏ ｗ ｕ ｎ ｃ ． ｃ ｏ ｍ． ｃ ｎ ＠ｃ ｇ ｐ
第２ ４卷 第４期 页码 ： ） ２ ０ ０ ９年 ８月（ １ ３ ８ ４～１ ３ ９ １
地 球 物 理 学 进 展
Ｐ Ｒ Ｏ Ｇ Ｒ Ｅ Ｓ Ｓ Ｉ Ｎ Ｇ Ｅ Ｏ Ｐ ＨＹ Ｓ Ｉ Ｃ Ｓ
Ｖ ｏ ｌ ． ２ ４，Ｎ ｏ ． ４ Ａ ｕ ． ２ ０ ０ ９ ｇ
高少武 ，赵 波 ， 贺振华 ， 等． 地震子波提取 方 法 研 究 进 展 ． 地 球 物 理 学 进 展， ： ： ／ ２ ０ ０ ９， ２ ４（ ４） １ ３ ８ ４～１ ３ ９ １， Ｄ Ｏ Ｉ １ ０． ３ ９ ６ ９ ． ｉ ｓ ｓ ｎ ． ｊ １ ０ ０ ４ ２ ９ ０ ３． ２ ０ ０ ９． ０ ４． ０ ２ ８． ） ， Ｇ ａ ｏＳ Ｗ， Ｚ ｈ ａ ｏＢ，Ｈ ｅＺ Ｈ， 犲 狋 犪 犾． Ｒ ｅ ｓ ｅ ａ ｒ ｃ ｈｐ ｒ ｏ ｒ ｅ ｓ ｓｏ ｆ ｓ ｅ ｉ ｓ ｍ ｉ ｃｗ ａ ｖ ｅ ｌ ｅ ｔ ｅ ｘ ｔ ｒ ａ ｃ ｔ ｉ ｏ ｎ ． 犘 狉 狅 狉 犲 狊 狊 犻 狀犌 犲 狅 犺 狊．（ ｉ ｎＣ ｈ ｉ ｎ ｅ ｓ ｅ ２ ０ ０ ９， ｇ 犵 狆 狔 （ ） ： ： ／ ２ ４ ４ １ ３ ８ ４～１ ３ ９ １， Ｄ Ｏ Ｉ １ ０． ３ ９ ６ ９ ． ｉ ｓ ｓ ｎ ． １ ０ ０ ４ ２ ９ ０ ３． ２ ０ ０ ９． ０ ４． ０ ２ ８． ｊ
， 犃 犫 狊 狋 狉 犪 犮 狋 ａ ｖ ｅ ｌ ｅ ｔ ｅ ｘ ｔ ｒ ａ ｃ ｔ ｉ ｏ ｎｏ ｆ ｓ ｅ ｉ ｓ ｍ ｉ ｃｄ ａ ｔ ａ ｉ ｓ ａｂ ａ ｓ ｉ ｃ ｔ ａ ｓ ｋ ｉ ｎ ｔ ｈ ｅｐ ｒ ｏ ｃ ｅ ｓ ｓ ｉ ｎ ｆ ｓ ｅ ｉ ｓ ｍ ｉ ｃｄ ａ ｔ ａｄ ｅ ｃ ｏ ｎ ｖ ｏ ｌ ｕ ｔ ｉ ｏ ｎ ｉ ｍ ｅ ｄ ａ ｎ ｃ ｅ Ｗ ｇｏ ｐ ， ｉ ｎ ｖ ｅ ｒ ｓ ｉ ｏ ｎａ ｎ ｄ ｆ ｏ ｒ ｗ ａ ｒ ｄｍ ｏ ｄ ｅ ｌ ｉ ｎ ． Ｆ ｏ ｒ ｔ ｈ ｅ ｓ ｅ ｉ ｓ ｍ ｉ ｃｄ ａ ｔ ａｐ ｒ ｏ ｃ ｅ ｓ ｓ ｉ ｎ ｆ ｈ ｉ ｈｒ ｅ ｓ ｏ ｌ ｕ ｔ ｉ ｏ ｎ ｈ ｉ ｈｓ ｉ ｎ ａ ｌ ｔ ｏ ｎ ｏ ｉ ｓ ｅ ｒ ａ ｔ ｉ ｏａ ｎ ｄｈ ｉ ｈ ｇ ｇｏ ｇ ｇ ｇ ｇ ， ， ｆ ｉ ｄ ｅ ｌ ｉ ｔ ｉ ｔｈ ａ ｓａ ｎｅ ｘ ｔ ｒ ｅ ｍ ｅｓ ｉ ｎ ｉ ｆ ｉ ｃ ａ ｎ ｃ ｅｔ ｏｅ ｓ ｔ ｉ ｍ ａ ｔ ｅｓ ｅ ｉ ｓ ｍ ｉ ｃｗ ａ ｖ ｅ ｌ ｅ ｔａ ｃ ｃ ｕ ｒ ａ ｔ ｅ ｌ ．Ｏ ｖ ｅ ｒｔ ｈ ｅｐ ａ ｓ ｔｔ ｗ ｅ ｎ ｔ ｅ ａ ｒ ｓ ｔ ｈ ｅ ｏ ｒ ｉ ｅ ｓ ｙ ｇ ｙ ｙｙ ， ， ｎ ｏ ｎ ｌ ｉ ｎ ｅ ａ ｒｓ ｃ ｉ ｅ ｎ ｃ ｅ ｓａ ｎ ｄｏ ｔ ｉ ｍ ｉ ｚ ａ ｔ ｉ ｏ ｎｈ ａ ｖ ｅｍ ａ ｄ ｅｇ ｒ ｅ ａ ｔｄ ｅ ｖ ｅ ｌ ｏ ｍ ｅ ｎ ｔ ａ ｎ ｄ ａ ｎ ｄｔ ｅ ｃ ｈ ｎ ｉ ｕ ｅ ｓ ｉ ｎｄ ｉ ｉ ｔ ａ ｌｓ ｉ ｎ ａ ｌｐ ｒ ｏ ｃ ｅ ｓ ｓ ｉ ｎ ｐ ｐ ｑ ｇ ｇ ｇ ， ｔ ｈ ｅ ｉ ｒｕ ｓ ｅ ｓｃ ｏ ｍ ｅｄ ｏ ｗ ｎ ｔ ｏ ｔ ｈ ｅ ｓ ｅ ｉ ｓ ｍ ｉ ｃｄ ａ ｔ ａｐ ｒ ｏ ｃ ｅ ｓ ｓ ｉ ｎ ａ ｔ ｔ ｈ ｅ ｓ ａ ｍ ｅ ｔ ｉ ｍ ｅ ｉ ｔ ｉ ｖ ｅ ｓｇ ｒ ｅ ａ ｔ ｉ ｍ ｅ ｔ ｕ ｓ ｔ ｏ ｔ ｈ ｅｄ ｅ ｖ ｅ ｌ ｏ ｍ ｅ ｎ ｔ ｏ ｆ ｔ ｈ ｅ ｇ ｇ ｐ ｐ ， ｗ ｅｆ ｉ ｒ ｓ ｔｂ ｒ ｉ ｅ ｆ ｌ ｘ ｌ ａ ｉ ｎｔ ｈ ｅｃ ｏ ｎ ｖ ｏ ｌ ｕ ｔ ｉ ｏ ｎｍ ｏ ｄ ｅ ｌ ｔ ｅ ｃ ｈ ｎ ｉ ｕ ｅａ ｎ ｄｍ ｅ ｔ ｈ ｏ ｄｏ ｆ ｔ ｈ ｅｓ ｅ ｉ ｓ ｍ ｉ ｃｗ ａ ｖ ｅ ｌ ｅ ｔ ｅ ｘ ｔ ｒ ａ ｃ ｔ ｉ ｏ ｎ ． Ｉ ｎｔ ｈ ｉ ｓｐ ａ ｅ ｒ ｙｅ ｐ ｑ ｐ ， ｗ ｉ ｔ ｈｔ ｈ ｅｃ ｌ ｏ ｓ ｅ ｓ ｔ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎｔ ｏｔ ｈ ｅｓ ｅ ｉ ｓ ｍ ｉ ｃｗ ａ ｖ ｅ ｌ ｅ ｔ ｔ ｈ ｅ ｎｗ ｅｇ ｉ ｖ ｅａｓ ｕ ｍｍ ａ ｒ ｆｓ ｅ ｉ ｓ ｍ ｉ ｃｗ ａ ｖ ｅ ｌ ｅ ｔｅ ｘ ｔ ｒ ａ ｃ ｔ ｉ ｏ ｎ ．Ｆ ｉ ｎ ａ ｌ ｌ ｅ ｙｏ ｙｗ ｓ ｉ ｍ ｌ ｉ ｓ ｃ ｕ ｓ ｓ ｔ ｈ ｅｐ ｒ ｏ ｓ ｅ ｃ ｔｏ ｆ ｔ ｈ ｅｓ ｅ ｉ ｓ ｍ ｉ ｃｗ ａ ｖ ｅ ｌ ｅ ｔ ｅ ｘ ｔ ｒ ａ ｃ ｔ ｉ ｏ ｎ ． ｐ ｙｄ ｐ ，ｗ ， ， 犓 犲 狑 狅 狉 犱 狊 ｅ ｉ ｓ ｍ ｉ ｃｗ ａ ｖ ｅ ｌ ｅ ｔ ａ ｖ ｅ ｌ ｅ ｔ ｅ ｘ ｔ ｒ ａ ｃ ｔ ｉ ｏ ｎ ｃ ｏ ｎ ｖ ｏ ｌ ｕ ｔ ｉ ｏ ｎｍ ｏ ｄ ｅ ｌ ｒ ｅ ｓ ｅ ａ ｒ ｃ ｈｐ ｒ ｏ ｒ ｅ ｓ ｓ ｓ ｇ 狔
４ 期
高少武 ， 等： 地震子波提取方法研究进展
１ ３ ８ ５
序列是最小相位的 ， 则可求得一个零相位或者最小 相位地震子波序列 ； 另外一种情况就是高阶统计量 假设 ： 假设反射系数序列为一个非高斯 、 平稳和统计 独立的随机过程 ， 假设子波为一个非因果 、 非最小相 位系统 ， 加性噪声为高斯噪声 ， 则可求得一个混合相 在利用地震记录的统计信息进行 位地震子波序列 ． 子波估计时 ， 二阶统计量中没有相位信息 ， 只有振幅 而高阶统计量不仅能保留地震子波相位信息 ， 信息 ； 而且能较好地压制高斯噪声 ．
在有测井数据的情况下 ， 利用 的关键制约因素之一 ． 测井数据计算出反射系数序列 ， 基于褶积模型 ， 利用 井旁地震道 ， 可以估算出地震子波 ． 在没有测井数据 的情况下 ， 使用统计性子波提取方法 ， 首先对反射系 数序列的分布做某 种 假 设 ， 然后利用地震数据的统 计特征进行子波估 计 ． 在没有任何先验知识的情况 下， 一种情况就是通常的地震褶积模型假设 （ 二阶统 ： 假设反 射 系 数 序 列 是 白 噪 的 ， 地震子波 计量假设 ）
（ １． 犆 犺 犲 狀 犱 狌犝 狀 犻 狏 犲 狉 狊 犻 狋 犲 犮 犺 狀 狅 犾 狅 犆 犺 犲 狀 犱 狌６ １ ０ ０ ５ ９， 犆 犺 犻 狀 犪； 犵 狔狅 犳犜 犵 狔， 犵 ２． 犌 犲 狅 犺 狊 犻 犮 犪 犾犜 犲 犮 犺 狀 犻 狌 犲犚 犲 狊 犲 犪 狉 犮 犺犆 犲 狀 狋 犲 狉， 犅 犌 犘， 犆犖犘 犆， 犣 犺 狌 狅 狕 犺 狅 狌０ ７ ２ ７ ５ ０， 犆 犺 犻 狀 犪） 狆 狔 狇
０ 引 言
随着地震勘探数据处理技术和解释技术的不断 发展 ， 从地震记录中 提 取 地 震 子 波 已 经 成 为 地 震 勘 探数字信号处理的 经 典 研 究 课 题 之 一 ， 尽管人们提 出了许多子波提取 方 法 技 术 ， 但是由于实际地震资 料往往很难满足子 波 提 取 的 一 些 基 本 假 设 ， 使得子 波提取技术成为影响地震资料处理质量进一步提高
反褶积的根本目的就是 波压缩或者消除处 理 过 程 ． 在地震反射数据已知条件下 ， 如果已知地震子波 ， 那 么根据褶积模型 ， 确定地震反射系数的过程 ， 这就是 确定子波反褶积 ； 如果不知地震子波 ， 那么根据褶积 模型 ， 对反射系数或者／和地震子波相位做出某种假 压缩地震子波 ， 提 高 地 震 数 据 分 辨 率 的 过 程， 这 设， 就是未知子波反褶积 ． 常规地震反演 ： 常规的地震反 根据褶 演就是在地震反射 数 据 和 地 震 子 波 已 知 时 ， 积模型 震 波 阻 抗 的 过程 ． 根据地震褶积模型 （ 和（ ， 一个方程三个变 １） ２） 如果知道其中两个变元 ， 那么就可以确定另外一 元， 这是一个 确 定 问 题 ， 正 演 模 拟、 确定子波反 个变元 ， 褶积 、 确定子波提取就是这种情况 ． 如果仅仅知道一 个变元 ， 且需要确 定 一 个 变 元 ， 这 是 一 个 欠 定 问 题， 那么就需要对第 三 个 变 元 和／或 者 待 确 定 变 元 做 出 未知子波 地 震 反 褶 积 和 估 计 子 波 提 取 就 某种假设 ， 是这种情况 ． 常规地 震 反 演 与 地 震 反 褶 积 的 情 况 类 似， 在反褶积的基础上确定地震波阻抗 ．
地震子波提取方法研究进展
２ 高少武１， ， 赵 波２， 贺振华１， 马玉宁２
（ 成都理工大学 ， 成都 ６ 东方地球物理公司物探技术研究中心 ， 涿州 ０ ） １． １ ０ ０ ５ ９； ７ ２ ７ ５ ０ ２．
摘 要 地震数据子波提取 ， 是地震资料反褶积处理 、 波阻抗反演以及 正 演 模 拟 的 基 础 工 作 ． 准确的地震子波估计技 术对于高分辨率 、 高信噪比 、 高保真度的地震勘 探 数 据 处 理 具 有 极 为 重 要 的 意 义 ． 在 过 去 的 二 三 十 年 中， 数字信号处 非线性理论 、 优化理论得到了很大发展 ， 这 些 为 地 震 子 波 提 取 技 术 和 方 法 注 入 了 新 的 巨 大 活 力． 本文首先简 理理论 、 要地说明了与地震子波紧密相关的褶积模型 ， 然 后 对 地 震 子 波 提 取 现 状 进 行 了 概 括 阐 述， 最后给出了地震子波提取 前景展望 ． 关键词 地震子波 ， 子波提取 ， 褶积模型 ， 研究进展 ： ／ ． ｉ ｓ ｓ ｎ ． １ ０ ０ ４ ２ ９ ０ ３． ２ ０ ０ ９． ０ ４． ０ ２ ８ 中图分类号 Ｐ Ｄ Ｏ Ｉ １ ０． ３ ９ ６ ９ ６ ３ １ 文献标识码 Ａ ｊ