肖 晓 段雅婷 胡双贵 汤井田 谢 勇 刘长生

(①中南大学有色金属成矿预测与地质环境监测教育部重点实验室，湖南长沙 410083； ②有色资源与地质灾害探查湖南省重点实验室，湖南长沙 410083； ③中南大学地球科学与信息物理学院，湖南长沙 410083； ④自然资源部覆盖区深部资源勘查工程技术创新中心，湖南长沙 410083； ⑤长沙航空职业技术学院，湖南长沙 410124)

0 引言

卫星、航空、地面、井中等多维空间的磁场数据测量技术迅速发展并逐渐成熟，配套的数据处理和解释技术也在不断更新和完善，这势必拓宽磁法勘探的应用领域，提高磁法的探测能力[1-4]。磁场数据三维反演不仅能够提供磁异常体的空间位置，还能获得磁性异常体的形状、物性分布等定量信息，在磁场数据处理和解释中具有重要意义，一直是该领域的研究热点问题。

目前的磁异常反演包括磁化率和磁性界面两方面的内容。磁化率反演首先固定反演网格，假设反演单元内磁化率为常数或者线性分布，建立反演目标函数，然后最小化反演目标函数，从而获得目标的磁化率分布[5-6]。Li等[6]利用一个或多个适当的加权函数将先验信息引入目标函数；Pilkington[7]对目标函数进行平滑及深度加权;Miguel[8]通过地质统计学模型实现了岩性约束下的重磁数据联合反演；Fregoso等[9]实现了重磁数据的交叉梯度联合反演。上述方法可定量描述地质构造的磁化率分布。磁性界面反演方法一般假定已知地下异常体的磁化率，通过反演得到异常体的边界和位置。……