基于光滑约束正则化的电磁波层析反演

林芳鹏 戴前伟 雷轶

层析反演的本质是解病态方程组,且系数矩阵通常为大型稀疏矩阵。正则化可提高解的稳定性,LSQR法可减小存储空间、加快计算效率,因此将二者结合提出了基于光滑约束正则化的层析反演算法。正则化参数在反演中有着至关重要的作用,采用L曲线法选择最佳正则化参数。列举了基于走时层析成像的数值计算实例,由最佳正则化参数下的速度剖面图可清晰地识别异常体特征,证明了基于光滑约束正则化的层析反演算法的有效性。