基于奇异值分解算法的瑞雷波相速度反演研究.pdfVIP

第42卷 第 2l期 人 民 长 江 Vo1．42．No．21 2011年 11月 Yangtze River NOV．。 2011 文章编号 ：1001—4179(2011)2l一0019一O3 基于奇异值分解算法的瑞雷波相速度反演研究 李 学 良，杨 长 春，王真 理 (中国科学院地质与地球物理研究所 ，北京 100029) 摘要：进行瑞雷波相速度反演时，在对其频散 曲线的线性化过程中，因误差的引入会导致系数矩阵奇异或近似 奇异。为解决此 问题，利用奇异值分解算法进行瑞雷波相速度反演，可提高低速软弱夹层的反演精度和可靠 性 ；反演时引入权重矩阵提高了数据分辨率 ；采用 自适应修改阻尼 因子以提高迭代效果并协调分辨率与解的 关系。实测资料试算结果表明，利用奇异值分解算法对瑞雷波相速度进行反演，不但具有稳定性好、精度高、 分辨能力强的特点，而且能 自动分层和反演地层参数。 关 键 词：瑞雷波；相速度；奇异值分解算法；低速软弱夹层 中图法分类号 ：P631 文献标志码 ：A 近年来 ，瑞雷波地震勘探技术在岩土工程 、环境与 阵提高数据分辨率 ；采用 自适应修改阻尼因子提高迭 工程地球物理研究中获得了广泛应用 ，尤其是最近几 代效果并协调分辨率和解 的关系，利用函数关系约束、 年 ，国内外对瑞雷波勘探技术的应用研究不断深入 ，也 区间约束和固定约束来提高解的唯一性。在反演过程 取得 了相当多的成果 I4。利用瑞雷波进行工程地质 中，正演模拟算法采用快速 6矩阵算法 ，有效地避 勘察主要是基于瑞雷波如下3个主要特性：① 在分层 免了高频数值溢出和精度丢失等问题 ，大大提高 了运 介质 中，瑞雷波具有频散特性 ，即瑞雷波的相速度是随 算速度。 频率 的变化而变化 。② 瑞雷波 的勘探深度和能量主 1 反演相速度的奇异值分解算法 要集 中在一个波长范围内，且波长不 同，其穿透深度也 不同。③ 瑞雷波传播速度与横波速度具有密切相关 1．1 反演方程 的推导 性 。 层状介质中瑞雷波频散方程的非线性隐式形式可 通过反演瑞雷波相速度 ，可有效地获取横波速度 表示为 和进行低速软弱夹层等地质结构划分 。对于瑞雷 F( ，， ，Vp，P，H)=0 (=1，2，…，m)(1) 波相速度的反演 ，与其它地球物理反演过程一样 ，由于 式中， 是第 个频率； ，是对应第_个『频率的瑞雷波 瑞雷波频散方程是一个高度非线性方程 ，在线性化 的 相速度 ； = ( ， ，…， ) 是横波速度 向量 ； 过程中会引入误差 ，同时在利用扰动偏微分法计算雅 可 比矩阵时也不可避免地会引入误差，这将会导致系 ： (Vp。， ：，…，V ) 是纵波速度 向量 ；P= (P。，P ， 数矩阵奇异或近似奇异。在处理这一类 问题 时，奇异 … ，P)是密度向量；H =(H ，H ，…，H)是厚度 向 值分解 SVD(SingularValueDecomposition)技术是最强 量 ；m为频率个数 ；n为地层层数。 有力的工具 。 前人研究结果表明，频散方程对密度和纵波速度 为此 ，本文在前人研究 的基础上利用泰勒级数展 的变化不敏感 ，故反演时假定密度和纵波速度已知 ；