MATLAB 共炮点地震合成,基于小波阈值的地震信号去噪处理

《基于小波阈值的地震信号去噪处理》由会员分享，可在线阅读，更多相关《基于小波阈值的地震信号去噪处理(4页珍藏版)》请在技术文库上搜索。

1、采 样 点 1 501、 道 采 样 间 隔 2 m s)中 提 取 了 一 道 地 震信 号 数 据 ,信 噪 比 为 利 用 常 规 小 波 阈 值 和 改 进 的 方法 进 行 去 噪 处 理 。 图 7是 从 实 际 地 震 数 据 中 提 取 出的 一 道 地 震 数 据 ,图 8是 常 规 小 波 阈 值 去 噪 后 的 波形 图 ,图 9是 用 前 面 提 到 的 方 法 去 噪 后 的 波 形 图 。图 8 常 规 小 波 阈 值 去 噪 后 波 形 图图 9 改 进 方 法 去 噪 后 波 形 图从 图 8和 图 9中 可 以 看 出 ,常 规 小 波 阈 值 和 改进 的 。

2、方 法 都 能 很 好 地 抑 制 噪 声 ,但 是 常 规 小 波 阈 值去 噪 不 仅 去 除 随 机 噪 声 还 去 除 部 分 有 效 信 号 ,改 进方 法 具 有 更 好 的 视 觉 效 果 ,从 而 证 明 了 该 方 法 的 有效 性 。3 结 论本 文 将 常 规 小 波 阈 值 去 噪 方 法 和 由 Heursure选 取 阈 值 进 行 二 次 小 波 分 层 阈 值 去 噪 方 法 应 用 到测 试 信 号 、 合 成 地 震 信 号 和 实 际 地 震 资 料 中 ,从 信噪 比 的 提 高 和 仿 真 结 果 可 以 看 出 ,对 含 有 高 频 分 量的 信 号。

3、 采 用 改 进 的 方 法 去 噪 ,不 仅 抑 制 了 噪 声 ,而且 很 好 地 保 留 了 高 频 有 效 信 号 ,具 有 更 好 的 视 觉 效果 ,从 而 证 明 了 该 方 法 的 有 效 性 。参 考 文 献1 何 樵 登 , 熊 维 纲 . 应 用 地 球 物 理 教 程 地 震 勘 探 . 北 京 :地 质 出 版社 , 19912 Jonesand L F, Levy S. Signal2to2noise ratio enhancement in multichan2nel seism ic data via the Karhunen2Loeve transform. 。

4、Geophysical Pros2pecting, 1987； 35: 12 323 冷 军 发 ,荆 双 喜 ,李 臻 . 基 于 小 波 多 分 辨 分 析 的 阈 值 去 噪 . 焦作 工 学 院 学 报 (自 然 科 学 版 ) , 2002； 21 (6) : 443 4464 胡 昌 华 ,李 国 华 ,刘 涛 ,等 . 基 于 MATLAB 6. X的 系 统 分 析 与 设计 小 波 分 析 . 西 安 :西 安 电 子 科 技 科 技 大 学 出 版 社 , 2004Se ism ic Da ta D e2no ising Ba sed on W avelet Thresho。

5、ld ingL IU Xia, PAN Hong2p ing3 , GAO Xiao2chun1 , WU W ei2long( School of Electrical & Information Engineering, Northeast Petroleum University, Daqing 163318, P. R. China；Hongwei Thermo2Electricity Power Plant1 , Daqing 163411, P. R. China) Abstract De2noising of seism ic data is a main task of geo。

6、physical exp loration, imp roving thesignal2to2noise ra2dio is turning to more and more important accordingly. The threshold is selected by Heursure method and softthreshold function, launches p rocess of second2decomposition on wavelet layer quantification to test signal, synthet2ic seism ic signal。

7、s and p ractical seism ic data. Compared with conventionalwavelet threshold de2noising, thismethodhas more advantage on de2noising to the signal containing high frequency components. Key words seism ic data wavelet transform threshold de2noising4527 科 学 技 术 与 工 程 10卷 1994-2011 China Academic Journal。

8、 Electronic Publishing House. All rights reserved. http:/www.cnki.net法 去 噪 后 ,信 号 的 高 频 有 用 信 息 相 对 较 好 地 保 留 下来 ,信 噪 比 比 常 规 小 波 阈 值 也 高 ,具 有 很 好 的 去 噪效 果 。2. 2 合 成 地 震 信 号 去 噪下 面 利 用 MATLAB对 地 震 资 料 进 行 模 拟 ,用 一个 低 频 ricker子 波 和 一 个 高 频 ricker子 波 叠 加 合 成共 炮 点 道 集 剖 面 图 ,具 体 参 数 :深 度 为 100 m, 反 射层 。

9、上 层 速 度 为 2 000 m / s,下 层 速 度 为 3 000 m / s,最小 偏 移 距 为 10 m,道 数 为 30,记 录 长 度 为 400 m s。图 3为 合 成 的 共 炮 点 道 集 剖 面 图 ,图 4为 加 入 5 db随 机 噪 声 后 的 共 炮 点 道 集 ,图 5是 常 规 小 波 阈 值 去噪 后 的 共 炮 点 道 集 ,图 6是 用 本 文 改 进 的 方 法 去 噪后 的 共 炮 点 道 集 。图 3 合 成 的 共 炮 点 道 集图 4 加 入 噪 声 后 的 共 炮 点 道 集常 规 小 波 阈 值 去 噪 和 改 进 方 法 去 噪 后。

10、 的 信 噪比 分 别 为 10. 804 8 db和 12. 071 2 db,信 噪 比 都 有所 提 高 ,从 图 5和 图 6中 也 可 以 看 出 ,当 有 效 信 号 含有 高 频 分 量 时 ,用 改 进 的 方 法 去 噪 视 觉 效 果 比 常 规小 波 阈 值 去 噪 的 视 觉 效 果 好 。图 5 常 规 小 波 阈 值 去 噪 后 的 共 炮 点 道 集图 6 改 进 方 法 去 噪 后 的 共 炮 点 道 集图 7 一 道 实 际 地 震 数 据2. 3 实 际 地 震 信 号 去 噪实 际 地 震 数 据 一 般 都 是 以 segy格 式 存 储 ,需 要进 行。

11、 转 换 ,去 掉 道 头 和 卷 头 ,对 真 正 的 数 据 进 行 处理 。 下 面 是 从 真 实 的 地 震 信 号 数 据 (201道 数 据 、 道352729期 刘 霞 ,等 :基 于 小 波 阈 值 的 地 震 信 号 去 噪 处 理 1994-2011 China Academic Journal Electronic Publishing House. All rights reserved. http:/www.cnki.nets ( k)中 恢 复 出 真 实 信 号 f ( k) 。图 1 分 解 过 程1. 2 二 次 小 波 分 层 阈 值 去 噪 方 法 的 。

12、实 现 步 骤在 阈 值 去 噪 算 法 中 ,硬 阈 值 和 软 阈 值 函 数 是 两种 最 常 用 的 阈 值 函 数 ,但 是 硬 阈 值 函 数 的 不 连 续 性使 去 噪 后 的 信 号 仍 然 含 有 明 显 的 噪 声 ,软 阈 值 方 法连 续 性 较 好 。 当 有 效 信 号 含 有 高 频 分 量 时 ,进 行 二次 小 波 分 解 ,本 文 采 用 软 阈 值 函 数 进 行 分 层 量 化 处理 。 此 方 法 的 实 现 步 骤 :1) 选 择 小 波 和 小 波 分 解 层 次 对 信 号 进 行 小 波分 解 ,得 到 低 频 小 波 系 数 和 高 频 小。

13、 波 系 数 ,其 中 高频 部 分 由 噪 声 和 信 号 组 成 ；2) 对 每 层 高 频 小 波 系 数 进 行 二 次 小 波 分 解 ,二次 分 解 后 低 尺 度 上 的 高 频 部 分 主 要 由 噪 声 控 制 ；3) 由 Heursure选 取 阈 值 和 软 阈 值 函 数 对 二 次分 解 后 的 高 频 系 数 进 行 处 理 ,然 后 重 构 ,重 构 后 作为 第 一 次 分 解 的 高 频 系 数 ；4) 将 处 理 后 的 高 频 系 数 与 第 一 次 分 解 得 到 的低 频 系 数 进 行 重 构 ,即 得 到 去 噪 后 的 信 号 。软 阈 值 函 。

14、数 为 = sgn ( j, k ) ( | j, k | - ) , | j, k | 0, | j, k | 。其 中 , j, k为 分 解 后 小 波 系 数 ； 为 施 加 阈 值 后 小 波系 数 ； 为 阈 值 。2 仿 真 实 验2. 1 测 试 信 号 去 噪信 号 B lock、 Bump s、 Heavy sine、 Dopp ler和 Quad2chirp 5种 信 号 几 乎 含 盖 了 图 像 、 光 谱 、 地 震 及 其 他科 学 领 域 所 遇 到 的 信 号 特 征 。 它 们 已 经 成 为 检 验小 波 去 噪 性 能 的 标 准 测 试 函 数 。下 。

15、面 对 含 有 噪 声 的 典 型 测 试 信 号 Quadchirp信号 进 行 试 验 ,信 号 的 长 度 为 n = 1 000。 图 2是 对 含噪 声 的 Quadchirp信 号 处 理 对 照 图 。 a) 原 始 Quad2chirp信 号 ； b) 含 噪 声 的 Quadchirp信 号 ,信 噪 比 为7 db； c) 常 规 小 波 阈 值 去 噪 波 形 ； d) 改 进 方 法 去 噪波 形 的 对 比 。 其 中 ,所 使 用 的 是 sym4小 波 。图 2 对 含 噪 声 的 Quadchirp信 号 去 噪常 规 小 波 阈 值 去 噪 后 的 信 噪 比。

16、 为 14. 950 9 db,改 进 方 法 去 噪 后 的 信 噪 比 为 17. 731 2 db,可 以 看出 ,两 种 方 法 去 噪 后 信 噪 比 都 提 高 了 ,大 部 分 噪 声都 得 到 了 抑 制 ,但 是 使 用 二 次 小 波 分 解 分 层 阈 值 方2527 科 学 技 术 与 工 程 10卷 1994-2011 China Academic Journal Electronic Publishing House. All rights reserved. http:/www.cnki.net第 10卷 第 29期 2010年 10月1671 1815 (201。

17、0) 2927251204科 学 技 术 与 工 程Science Technology and EngineeringVol110 No129 Oct12010Z 2010 Sci1Tech1Engng1石 油 技 术基 于 小 波 阈 值 的 地 震 信 号 去 噪 处 理刘 霞 潘 洪 屏 3 高 晓 春 1 吴 伟 龙(东 北 石 油 大 学 电 气 信 息 工 程 学 院 ,大 庆 163318；宏 伟 热 电 厂 1 ,大 庆 163411)摘 要 地 震 信 号 去 噪 处 理 是 地 震 勘 探 所 要 解 决 的 一 项 主 要 任 务 ,提 高 地 震 信 号 信 噪 比 。

18、也 相 应 显 得 尤 为 重 要 。 提 出 由 Heur2sure选 取 阈 值 和 软 阈 值 函 数 对 测 试 信 号 、 合 成 地 震 信 号 及 实 际 地 震 数 据 进 行 二 次 小 波 分 解 分 层 量 化 处 理 ,与 常 规 小 波 阈 值 去噪 相 比 ,对 含 有 高 频 分 量 的 信 号 ,该 方 法 具 有 很 好 的 去 噪 效 果 。关 键 词 地 震 信 号 小 波 变 换 阈 值 去 噪中 图 法 分 类 号 TE132. 14； 文 献 标 志 码 A2010年 7月 23日 收 到第 一 作 者 简 介 :刘 霞 ( 1970 )女 ,硕 士。

19、 ,副 教 授 ,研 究 方 向 :智 能控 制 与 数 字 信 号 处 理 。3 通 信 作 者 简 介 :潘 洪 屏 (1980 ) ,女 ,硕 士 研 究 生 ,研 究 方 向 :智 能信 息 与 控 制 。 E2mail: garychong sohu. com。野 外 采 集 的 地 震 资 料 中 不 可 避 免 包 含 一 些 规则 的 或 不 规 则 的 干 扰 噪 声 1 ,这 些 干 扰 对 地 质 资 料解 释 将 产 生 很 大 影 响 。 因 此 ,降 噪 成 为 地 震 信 号 处理 的 重 要 环 节 。 而 小 波 变 换 具 有 很 好 的 时 频 局 部化 。

20、能 力 ,将 其 应 用 于 地 震 信 号 去 噪 2 处 理 ,有 着 良好 的 应 用 价 值 和 研 究 前 景 。对 地 震 信 号 去 噪 的 方 法 有 多 种 ,在 小 波 领 域 有模 极 大 值 去 噪 、 小 波 阈 值 去 噪 和 小 波 包 去 噪 ,其 中阈 值 去 噪 最 为 常 用 ,方 法 简 单 。 在 阈 值 去 噪 方 法中 ,如 何 选 取 阈 值 和 阈 值 函 数 是 至 关 重 要 的 。 阈 值选 取 的 合 适 与 否 直 接 影 响 信 号 去 噪 的 效 果 和 重 构信 号 的 失 真 程 度 ,如 果 阈 值 选 取 过 大 ,虽 然。

21、 能 够 减少 重 构 信 号 中 残 留 的 噪 声 成 分 ,但 会 使 信 号 有 较 大的 失 真 ；反 之 ,降 低 阈 值 能 减 小 重 构 信 号 的 失 真 ,但恢 复 信 号 中 残 留 的 噪 声 增 多 。本 文 主 要 介 绍 对 含 有 高 频 有 效 信 号 进 行 二 次小 波 分 层 阈 值 去 噪 的 基 本 思 想 和 实 现 步 骤 ,并 由Heursure选 取 阈 值 和 软 阈 值 函 数 对 测 试 信 号 、 合 成地 震 信 号 及 实 际 地 震 数 据 进 行 二 次 小 波 分 解 分 层量 化 处 理 。1 二 次 小 波 分 层 阈。

22、 值 去 噪 方 法 的 基 本 思 想和 实 现 步 骤1. 1 二 次 小 波 分 层 阈 值 去 噪 方 法 的 基 本 思 想二 次 小 波 分 解 分 层 阈 值 去 噪 3, 4 主 要 是 对 含 有高 频 有 效 信 号 进 行 二 次 小 波 分 解 ,利 用 有 效 信 号 和噪 声 在 小 波 变 换 下 奇 异 性 截 然 不 同 的 表 现 特 征 来分 层 去 除 噪 声 保 留 有 效 信 号 。 一 个 含 噪 信 号 的 一维 信 号 模 型 可 以 表 示 成 如 下 形 式 ,即 s ( k) = f ( k) + e ( k) ； k = 0, 1, ,。

23、 n - 1,其 中 , f ( k)为 有 效 信 号 ,e ( k)为 噪 声 , s ( k)为 含 噪 声 的 信 号 。 当 有 效 信 号 含有 高 频 分 量 时 ,对 第 一 次 分 解 后 的 高 频 系 数 进 行 二次 小 波 分 解 ,采 用 分 层 阈 值 处 理 。 在 实 际 的 工 程 应用 中 ,有 用 信 号 通 常 表 现 为 低 频 信 号 或 是 一 些 比 较平 稳 的 信 号 ,但 也 可 能 含 有 高 频 信 号 ,而 噪 声 信 号通 常 表 现 为 高 频 信 号 ,所 以 消 噪 过 程 可 按 如 下 方 法进 行 处 理 :首 先 对。

24、 信 号 进 行 一 次 小 波 分 解 (如 进 行二 层 分 解 ) ,则 噪 声 部 分 主 要 包 含 在 高 频 小 波 系 数CD1 , CD2 中 ,其 中 包 含 部 分 有 效 信 号 ,对 高 频 系 数进 行 二 次 小 波 分 解 (如 进 行 二 层 分 解 ) ,分 解 过 程 如图 1所 示 ,然 后 由 Heursure选 取 阈 值 进 行 分 层 阈 值处 理 ,重 构 信 号 即 可 达 到 消 噪 的 目 的 。 对 信 号 s ( k)消 噪 的 目 的 就 是 要 抑 制 信 号 中 的 噪 声 部 分 ,从 而 在 1994-2011 China Academic Journal Electronic Publishing House. All rights reserved. http:/www.cnki.net。