基于匹配追踪算法子波分解技术在薄互层储层预测中的应用

１ 子波分解原理及过 程实现
随着 勘探 开发 的不 断深入 ， 勘探 开发 目标 逐步 转 向薄互 层 。因此 ， 薄互层 储层 砂体 预测 技术 已成 为 国 内、 各 油 田的重 要研 究 课 题 。多年 来 ， 球 外 地 物理学家 对利 用地 震 资 料解 决 薄 层 问题 开 展 了广 泛的研究 ， 并取 得 了显 著成 就 ｌ ， 利 用 地震 资 － ｌ 为 料对 薄层 的研 究 奠定 了理 论 和 方法 基 础 。但 多 数 利 用地震 资料进 行 薄 互层 预 测 技 术 均是 基 于单 一
布规律 ， 别薄互层 储层 ， 识 为油藏 的地质 及地 球物 理 解释 提供 有 力的 支持 。
关 键词 ：匹配追踪 ； 波 ； 互层 ； 解与 重构 ； 震数 据体 子 薄 分 地
中图分 类号 ：Ｐ６ １４ ４ ３ ． ３ 文 献标识 码 ：Ａ
布进行 预测 。
０ 前 言
第３ 卷 第６ ２ 期
物 探 化 探 计 算技 术
２１ 年 １ 月 ０ ０ １
文章编 号 ：１０ — １４ （０ ０ ０ － ０ ４ —０ ｏ １ ７ ９ ２ １ ）６ ＿ ６ １ ４
基 于 匹 配 追 踪 算 法 子 波 分 解 技 术 在 薄 互 层 储 层 预 测 中 的 应 用
个 多子波 地震 道 模 型 【 Ｊ １ 。在 该 模 型 中 ， 博 地震 子
变化 。如储层 和 非储层 ， 油气储 层段 和不 含油气 含 储层段 的地震 响应 不 同。 因此 ， 种单 子波 褶积模 这 型并不 符合 实 际情 况 。作 者 主 要是 利 用 基 于 匹 配 追踪 （ ＭＰ算 法 ） Ａ 区块地 震 道 数 据进 行 子 波分 对
波 的每 一次反射 ， 其形 状都 可能会 发生 改变 。在每 次反 射 时 ， 形状 都 可 能不 同。 因此 ， 个 地 震道 其 每
就是 所有这 些反 射 的对 应 组合叠 加 ， ： 即
Ⅳ
解， 并在 此基 础上 ， 该 地 区 砂 泥岩 薄 互 层 储 层展 对
收稿 日期 ：２ １ ０ ０—０ ５—１ ７ 改 回 日期 ：２ ０一ｌ ０１ ０—１ ３
何胡军 王秋语２ ， ，程会明２
（． １ 成都理 工 大学 能 源学 院 ，四川 成都 ６０５ ； １０ ９
２ 中国石 化胜 利 油田 地质科 学研 究 院 ，山 东 东营 ．
摘
２７ １ ） ５０５
要 ：薄互层储 层砂 体预 测技 术 的研 究 已成 为 国 内、 外各 油田的重要 研 究课题 ， 已取得较 多 并
同物理特 性 的地 层 时 ， 形 状 会 发生 不 同 的变 化 ， 其
导致 这个假 设 与实 际情 况 尚有 一 定 的距 离 。 因为 在现 实 中 ， 由于地 层 在 横 向上 有 变化 ， 震 子 波 在 地 穿过 不 同物理特 性 的地 层后 ， 由于不 同地层 对地 震
子波 的响应 不 同 ， 会使 地震 子波 的形状 发生 不 同的
研 究成果 。但 多数 利 用地震 资料 进行 薄互层 预 测技 术 ， 均基 于单 一地 震 子 波与 地震 反 射 系数 函 数 的褶积假 设 ， 实际上 地震 子波在 穿过 不 同物理 特 性地 层 时 ， 其形 状会 发 生不 同变化 ， 导致 该 假
设与 实际情 况 尚有一 定距 离， 致使 薄 互层储 层预 测 的精 度 降低 。这里 主要 论 述 的是 基 于 匹 配追 踪算 法子 波分 解技 术识 别 薄互层储 层 。通过 研 究表 明， 该技 术 能较 好 地 分辨 出储 层 在 空 间的 展
层 和非 储层 ， 油 气层 和不 含 油 气地 层 ， 地 震信 含 对 号（ 子波 ） 响应 也 不 同 ， 的 与之 相应 的地震 信号 （ 子
波） 的形 状或频 率 也不 同 ， 而导 致这 一 假 设 与实 从 际情 况 尚有 一 定 的 距 离 ¨ 。为 此 ， 人 设计 了一 前
地震 子波 的假设 ， 而在 实 际 中， 震 子波 在 穿 过 不 地
１ １ 多子波 地震道 模型 ．
地 震资 料解 释 的基本 模 型 ， 是单 一地震 子波 与
地 层反 射系数 序 列 的褶 积 。常规 的地震 信号 处理
和解 释中所用 的 褶 积 、 褶 积 ， 反 以及 地震 波阻 抗 反 演 等 ， 是基 于单 一地 震 子 波 的假 设 （ ｈ ｒ 、 ｏ — 都 Ｓ ｅ ｆ Ｒ ｂ ｉ ｅｔ １９ ｒ ９ ９年 ） Ｅ， 。而实 际地 震 波 的频 率 会 随着 深度 的增 加 而降低 ， 且 同一 地 震 数 据 ， 同类 型 的储 而 不
为对应 的反射 系数序列 ； （）为噪音项 。 ６ｆ
与常规单一 子波地震道 模型相 比较 ， 该模 型 中 的地震 子波形状 是随 传播 时 间 和空 间 的变化 而 变 化 的， 即多子波地震道 模型 中的地震子 波形状 不是 固定不 变的 。而是 随实 际地层 情况 ， 以及 含油性 变 化而变化 的。因此 ， 多子 波地震 道模 型更 符合地震 波传播 的实 际情况 。 １２ 子波分解 过程实现—— Ｍ ． Ｐ算 法 子 波 分解是 把 一个地 震道 分解 成 不 同能 量 的 地震子 波 的集合 ， 体实现过 程是将输 人 的地 震数 具 据体 中所 给定 的数据 段 ， 分解成 不 同能量 的子 波分 量， 而不 同能量 的子波 分量是基 于输入地 震数据 段 通 过统计计算 而得到 的。 作者在本 文的子波 分解过程 的实现 ， 主要 采用
ｔ ）＝
（） ｒ（）＋ｂ ￡ ｆ ￡ （）
（） １
６２ ４
物探 化探 计 算技 术
３ ２卷
式 中
（）为第 ｎ个 反 射层 的地震 子 波 ； ｔ ｔ ｒ（ ）
ＲＳｔ （）＝ ［ Ｓ ｔ ， ｔ］ ｔ （）ｇ （）ｇ （）＋
Ｒ ．ｔ ｓ ） （
，一 ｖ １