叶益信，李泽林，付 宸，丁尚见

（1.东华理工大学 放射性地质与勘探技术国防重点学科实验室，南昌 330013；2.中国地质大学 地球内部多尺度成像湖北省重点实验室，武汉 430074）

利用阻尼型高斯牛顿法的激发极化数据聚焦反演

叶益信1，2，李泽林1，付 宸1，丁尚见1

（1.东华理工大学 放射性地质与勘探技术国防重点学科实验室，南昌 330013；2.中国地质大学 地球内部多尺度成像湖北省重点实验室，武汉 430074）

利用阻尼型高斯牛顿法求解反演目标函数，并将粗糙度矩阵约束和最小支撑泛函约束以及参考模型和模型参数界限约束引入到目标函数中，同时结合预条件共轭梯度法进行反演迭代，实现了激发极化数据三维聚焦反演快速计算。通过对几个典型模型的试算，并与传统光滑模型约束的反演结果进行比较，表明该反演算法的可靠性和稳定性较好，并且算法速度快，占用内存低且易于并行化。

激发极化；高斯牛顿法；聚焦反演；最小支撑泛函

0 引言

激发极化（IP）法的应用范围非常广泛，无论是在金属矿和非金属矿产勘查［1－2］，还是在寻找地下水资源和地热田方面［3－4］，都获得了成功地应用，近年来又拓展到环境监测领域［5－6］。国内、外有关学者对激发极化法正反演也有较多的研究，Pelton等［7］利用参考数据库中数据插值计算偏导数，实现了激发极化二维最小二乘反演方法；日本的Sasaki博士［8］在前人的基础上发展了有限元电阻率／极化率数据的二维最小二乘反演方法；Ellis等［9］也将解的先验信息加入目标函数中，使反演结果更符合实际；……