基于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