一种用于Occam反演中搜索拉格朗日乘子的方法  

A Method Used for Searching Lagrange Multiplier in Occam Inversion

作　　者：刘俊峰[1] 邓居智[1] 陈辉[1] 汤磊[1] 张寒韬[1] 

机构地区：[1]东华理工大学 放射性地质与勘探技术国防重点学科实验室,江西抚州344000

出　　处：《工程地球物理学报》2013年第3期344-350,共7页Chinese Journal of Engineering Geophysics

基　　金：国家自然科学基金(编号:41164003);江西省井冈之星青年科学家培养计划;国家科技支撑计划(编号:2011BAB04B03);东华理工大学研究生创新基金(编号:DYCA11008)联合资助

摘　　要：Occam反演法由于其算法稳定,对初始条件要求不高,反演效果较好,在大地电磁反演中运用较多。不过其每次迭代都需要进行偏导数计算和大量的模型计算,以便搜索到最佳的拉格朗日乘子,这造成了计算量和计算时间的增加。Occam反演在每次迭代寻找最佳模型的过程中需要搜索合适的拉格朗日乘子使拟合差最小,搜索的方法一般使用进退法和扫描法,鲜少使用其他一维搜索方法。本文将牛顿迭代搜索法和二分法组合一起用于拉格朗日乘子的搜索,取得了较好的结果,减少了模型的搜索量,在一定程度上提高了计算速度。Occam inversion method is widely used in the magnetotelluric inversion because of its algorithm stability,low requirements on the initial conditions and the inversion effect.But in each iteration,partial derivative calculation and a large number of model calculation are required in order to search the optimal Lagrange multiplier,which results in an increase in the amount of calculation and the calculation time.When Occam inversion finds the best model in each iteration,the appropriate Lagrange multiplier is needed to seek the suitable Lagrange multiplier so as to minimize the fitting error.The advance-retreat method and scanning method are generally used.In this paper,the combination of Newton iterative search method and dichotomy used in search of Lagrange multiplier achieve a good result.

关 键 词：OCCAM反演 拉格朗日乘子 一维搜索 

分 类 号：P631.4[天文地球—地质矿产勘探]