摘要：在频率域弹性波有限元正演方程的基础上,依据匹配函数(也就是观测数据和正演数据残差的二次范数)最小的准则,用矩阵压缩存储与LU分解技术来存储和求解频率域正演方程中的大型稀疏复系数矩阵、用可调阻尼因子的Levenberg-Marquard方法求解反演方程组,直接求取地下介质的弹性波速度,导出了频率域弹性波有限元最小二乘反演算法.为了利用地下地质体的分布规律,减少反演所求的未知数个数,本文又提出了规则地质块体建模方法引入到反演中来.经数值模型验证,在噪声干扰很大(噪声达到50%)或初始模型与真实模型相差很大的情况下,反演也能取得很满意的效果,证明本方法具有很好的抗噪性与"强壮性".

关键词：频率域 弹性波 有限元 反演 矩阵压缩存储 

单位：中国科学院地质与地球物理研究所，北京100029