基于MATLAB语言的二维大地电磁OCCAM快速反演

以TM模式为例,基于MATLAB语言对二维大地电磁快速正反演算法进行了研究.针对有限元正演计算中多次矩阵赋值和多频率计算的特点,引入了矢量化编程和并行计算两种策略,编制了实现快速正演的MATLAB程序.以此为基础,开展了基于OCCAM理论的反演算法研究.首先,使用拟正演手段对偏导数矩阵进行快速求取.其次,为平衡模型参数与拟合误差的关系,采用了固定步长减小的拉格朗日乘子,简化了拉格朗日乘子的求解过程,提高了反演计算的效率.最后,对模型数据进行反演,结果显示反演算法收敛快速、稳定性好.