利用地质规则块体建模方法的频率域有限元弹性波速度反演

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