cauchy稀疏约束 matlab,基于Cauchy稀疏约束的同时迭代反演反射系数与地震子波的研究...

基于Cauchy稀疏约束的同时迭代反演反射系数与地震子波的研究

第 1 3卷第 1期 重 庆科 技学 院学报 ( 自然科学 版 ) 2 0 1 1年 2月 基于 C a u c h y稀疏约束的同时迭代反演反射系数与 地震子波的研究 李力民 姚泳 汐 ( 1 ． 西 南石 油大 学， 成都6 1 0 5 0 0 ；2 ． 中石 油 西南油 气 田分公 司重庆 气矿 ， 重庆 4 0 0 0 2 1 ) 摘要 ： 在二阶统计量实现地震反褶积基础 1 2 ， 作反射系数为高斯分布的随机序列 、 子波为最 小相位 的假设 ． 并在稀 疏反 演中以 C a u c h y约束为准则 ， 在规 则化 过程中 ， 微妙 的引入稀疏反褶积 的原 始问题 ． 作 预条件处理 。 采用预条件 共轭梯度法迭代反演 同时得到地震子波和反射系数序列 。在具体方法实现中 。 以传统预测反褶积结果作为反射系数 的迭代初始值 ， 初始子波用稀疏约束反演方法得到。利用理论模型进行试算 ， 表明该算法精度高 、 稳定 、 收敛速度快。 关键词 ： 迭代反演； C a u c h y 稀疏约束 ； 预条件共轭梯 度法 ： 预测反褶积 中图分类号： P 6 3 1 文献标 识码 ： A 文章编号 ： 1 6 7 3 — 1 9 8 0 ( 2 0 1 1 ) 0 1 - 0 l 8 4 — 0 5 地震勘 探 中基于 褶积模 型 的地震记 录可 以简单 的认 为是震 源子 波与地 下反射 系数 序列 相褶 积的结 果 ． 而地震反 褶积 技术 就是对 子波进 行处 理 ， 从 地震 记录中消除子波的影响， 得到实际反射系数序列． 提 高地震资料分辨率和抑制多次波 的有效技术I ” 。地 下 反 射系 数序 列 可 以认 为是 一 系列 脉 冲序 列 ， 在 时 间上 ， 它 的位置 具体 体现 了反射层 面确 切 的位置 ， 它 的幅值 能体 现 出反 射界 面 反射 波 能 量 的强 弱 ， 此 时 的反褶 积可 称 为脉 冲反 褶 积或 者 称 为稀 疏 反褶 积 ． 统称 为最／ b 2乘 反褶积 。 地 震 反 褶 积 实 际 上 是 一 个 盲 过 程I 2 1 ， 其 中反 射 系数 和子 波 都是 未知 量 。 因此 ， 它 实 质上 属 于 盲信 号处理过程 。 并且盲信号处理也是统计信号处理 的 范畴 ． 需 要进 行统计 性假设 I 。一类是 基 于二阶 统计 量 自相关 方法 ， 由地震 道 反 演 出地震 子 波 和 反射 系 数 ， 它必 须作地 震子 波为最 小相 位子 波 ， 反射 系数 序 列 为高斯 白噪 序列 时才能取 得好 的结 果 。另一类 是 基 于高 阶统 计 量 的非线 性 反演 实 现地 震 盲 反褶 积 。 高 阶统 计量 中包 含 了子 波 相位 的 信息 ， 可估 计 出反 射 系数序 列和非 最小相 位予 波I 。 在 地 震 资料数 字 处 理 中 ． 用 的最 多 的仍 为 基于 二 阶统 计量 的反褶 积方 法 ． 如预测反 褶积 ， 稀疏 脉 冲 反褶 积 。本文是 在 常规 的反射 系数 为随机 的高斯 分 布序列 ，子波 为最小 相位 的前 提下 ，仅从地 震道 出 发 ， 应用 C a u c h y准则 ， 对反 射 系 数序 列进 行 稀 疏 约 束 ，并采用预条件共轭梯度法对反射系数和子波同 时进行迭代反演 。本文将传统的预测反褶积结果 作 为反射 系数 序列 迭代初 值 ，与脉 冲反 褶积 的结果 相 比 ， 反演 算 法实 现 中迭 代 次数 少 ， 反 演精 度 较 高 。 子 波初 值 由 C a u c h y准 则稀 疏 反 演预 条 件 共轭 梯 度 法求 得 ， 保 证 了在计 算过程 中两者能 同步稳 定 精 度 较高且 收敛 速度 快 。 1 反演算法与原理 1 ． 1 基 于二 阶统计 量 的最小 二乘 反演原 理 地 震记 录的褶 积模 型可 以表示 为 ： S + N k ( 1 ) 式 中： S k ( = l ， 2 ⋯， n ) 一地震记 录； Q ( k = l ， 2 ⋯， h a ) 一 反 射 系 数 序列 ； W ( k = l ， 2 ⋯ ， n ) 一 子 波 序 列 ； ( k = l ， 2 ⋯ ， r ／ n ) 一 噪声 序列 。当噪音 比较微弱 时 ， 我 们 近似 值为 0 ， 式 ( 1 ) 写成 褶积矩 阵 的形式 为 ： C S ． ( 2 ) 式 中 ： 一 地 震 记 录 向量 ； 一 子 波 的褶 积矩 阵 ； q 一 反射系数序列列向量。 在最小二乘准则的约束下。 可 收 稿 日期 ： 2 0 1 0 — 0 8 — 2 5 作者简介 ： 李力 民( 1 9 8 5 一 ) ， 男 ， 满族 ， 辽宁葫芦 岛人 ， 西南石油 大学在读硕士研究生 ， 研究方向为地震信号处理与地震 反演方法 研 究 。 · 1 8 4· 李 力 民 ． 姚 泳 汐 ： 基 于 Ca u c h y稀 疏 约 束 的 同时迭代 反 演反 射 系数 与地 震子 波 的研 究 得 到 ： = s 一∑ 辱 k = l ， 2 ⋯， 目标 函数 形式 为 ： p ， ( ) = 1 1 n 1 + 1 ) ( 1 o ) ( 3 ) 式 ( 9 ) 写 成矩 阵 形式 为： ．，： ∑ p 1 l J c4 式 中： p ( M ) = 1 u ； 一每个记 录值的标准方差 。作 简 化处 理 ， 取 同值 o r ， 其矩 阵 形式 表示 为 ： T l g — S I l 2 ( 5) 2 0 - 可 得正 则化 方 程： g = S ( 6) 式 中 ： 一 子 波 褶积 矩 阵 的转 置 。 式 ( 6 ) 为 一 个 线 性 形 式 的 T o e p l i t z矩 阵 ， 可 应 用 L e v i s o n递 推 、 矩 阵 求 逆 、 共 轭 梯 度 法 等 方 法 求 解 [7 ,8 1 。 但 是 ， 方 程 中参 数 的 微 小 变 动 都 会 对 最 终 方 程 的求 解 造 成 很 大 的影 响 ． 此 时 ， 需 要做 一 些 规 则化 处理 。 1 _ 2规则 化 策略 在式 ( 6 ) 的 自相 关 矩 阵 的 对 角线 元 素上 加 入 一 个 正 则 化 因子 ， 其 目的是 使 求 得 的算 子 能 够 达 到 数 值 稳定 ， 设 计 目标 函数 为 ： = pl[ 卜 mZ P l [ 】 ( 7 ) 对 式 ( 7 ) 进 行 极 小 化 即 为规 则 化 策 略 。得 正 则 方 程 为 ： ( C w + tx I ) q = C w S ( 8) 式 中 ： ／ x为 正 则 化 因子 ， ： 一 『O “ n； 和 分 别 表 示 反 射 系数 和噪 音 的标 准 差 。式 ( 8 ) 为 阻 尼 最 小 二乘 反 演 方法 ， 与式 ( 6 ) 相 比 ， 算 法 的求 解 稳 定 性 较 高 。式 ( 8 ) 同样可以用矩阵求逆 、 L e v i