初至波层析影响因素分析.docVIP

初至波层析影响因素分析刘玉柱 3董良国( 同济大学海洋地质国家 初至波层析影响因素分析 刘玉柱 3 董良国 ( 同济大学海洋地质国家重点实验室) 刘玉柱 ,董良国 . 初至波层析影响因素分析 . 石油地球物理勘探 ,2007 ,42 (5) :544～553 摘要 在山地地震勘探中 ,目前广泛采用初至波走时层析方法反演起伏地表速度结构 ,但其最终反演结果受许 多因素影响 。本文通过大量数值试验对影响表层模型层析反演精度的方程解法 、先验信息利用 、观测系统 、初 始模型及平滑因子等主要因素进行了详细的分析对比 ,总结出了最佳参数的选择原则和范围 ,并建立了内部约 束方程 。将本文方法应用于实际资料处理 ,取得了较好的成像效果 。 关键词 层析反演 表层结构 初至波 先验信息 正则化 影响因素 1 引言 2 参数分析 在山地地震勘探中 ,近地表速度结构反演是一个 亟待却又很难解决的问题 。尽管国内外学者对该问题 进行了大量的研究[ 1～8 ] ,但至今仍没有很好地解决 , 从而严重影响表层校正效果及成像质量 。通过细致 调研 ,我们认为应充分利用地震波的走时和波形 (包 括振幅 、相位 、频率等) 信息 ,结合线性反演 、拟线性反 演 、完全非线性反演等多种反演手段 ,在地表调查资 料和已知地质背景知识等先验信息的约束下 ,采用多 信息 、多方法联合反演的方法可望有效地解决此问 题 。尽管联合反演是最终目标 ,但它仍依赖于单一反 演技术的提高 。目前广泛采用的单一反演方法包括 折射层析[ 9～10 ] 和初至波走时层析[ 2 ,4～7 ,11 ,12 ] 两种 。折 射层析以层状模型假设为基础 ,尽管折射层析有时可 以得到比较好的静校正效果 ,但层状假设的前提条件 并不总能得到满足 。初至波走时层析没有层状模型 假设 ,它适用于任意表层速度结构 ,反演结果能较好 地反映表层速度的低频趋势 。鉴于初至层析仍存在 反演精度低 、稳定性差 、反演结果受诸多因素影响等 问题 ,本文通过大量数值试验 ,对影响反演精度的主 要因素进行详细对比分析 ,旨在探讨各参数影响层析 反演效果的规律 ,总结出各参数优选的原则和范围 , 达到提高初至层析反演精度及反演稳定性的目的 。 在线性假设下 ,初至波走时层析成像简化为求 解大规模线性方程组 LΔs = Δt ( 1) 其中 : L 为 m ×n 维矩阵 , 矩阵元素 l kj 代表第 k 条射 线在第 j 个模型参数单元内的长度 ;Δs 是长度为 n 的列向量 , 代表模型慢度参数修正量 ;Δt 为 m 维列 向量 , 代表观测的初至波走时与理论计算的初至波 走时残差[ 13 ] 。 不考虑拾取误差与射线追踪误差 ,概括起来影 响层析反演结果的因素包括 : ①方程解法 ; ②先验信 息的利用 ; ③观测系统 ; ④初始模型 ; ⑤平滑因子 。 下面通过理论分析 、结合数值实验探讨上述各因素 对层析反演结果的影响 。 2 . 1 方程解法 求解方程组式 (1) 通常有反投影法 (B P T) 、代数 重建法 ( A R T) 、同时迭代重建法 ( SIR T) 、奇异值分 解法 ( SVD) 和投影类迭代方法等[ 14 ] 多种算法 。其 中反投影法是最传统的一种方法 ,其基本思想是将 旅行时残差按射线在模型参数单元内的覆盖权数进 行反投影 。这种方法虽然分辨率低于代数重建法和 投影法 ,但具有较高的计算效率 ,而且稳定性好 。计 算公式为 3 上海市同济大学海洋与地球科学学院 ,200092 本文于 2006 年 11 月 23 日收到 ,修改稿于 2007 年 6 月 8 日收到 。 本课题获得教育部“新世纪优秀人才支持计划”( N CE 与“国家重点基础研究发展规划项目”( 2006 CB202402) 资助 。 巨大的稀疏矩阵 ,而且矩阵的条件数有时很大 (病态问题比较严重) ,则 巨大的稀疏矩阵 ,而且矩阵的条件数有时很大 (病态 问题比较严重) ,则求解方程组时收敛速度慢且不稳 定 ;L SC G 法改变了对该矢量的计算方式 ,使其能够 处理矩阵为病态 的情 况 ; L SQ R 法 是利 用 L a nczo s 法求解最小二乘问题的一种投影法 ,它可压制数据 误差的传递 ,收敛速度较快 ,适合求解大型病态线性 代数方程组[ 14 ] 。 图 1a 与图 1 b 是基于实际数据的拾取初至采用 相同的初始模型与反演参数 ,分别利用反投影法与 L SQ R 算法得到的层析成像结果 。图 2 是图 1a 与 图 1 b 地表以下 40 m 处的速度切片对比图 。不考虑 精度 ,可以看出 L SQ R 算法确实比反投影法具有更