胡正旺 杜劲松*③ 孙石达 陈 超

(①中国地质大学(武汉)地球物理与空间信息学院，湖北武汉 430074；②地球内部多尺度成像湖北省重点实验室，湖北武汉 430074；③地质过程与矿产资源国家重点实验室，湖北武汉 430074；④河北地质大学勘查技术与工程学院，河北石家庄 050031)

0 引言

在地球物理学方法体系中，磁法发展最早，至今仍被广泛应用于地质调查、资源勘探、军事与国防建设、工程勘察与环境监测等诸多领域[1-3]。地磁场测量的平台、方式以及物理量多种多样，目前磁法勘探普遍采用标量磁力仪，基于航空、地面与船载等移动平台测量地磁场总强度的空间分布[4]。随着地磁场测量精度和测点定位精度的大幅提高，以及计算机与计算技术的快速发展，对深部目标及小型磁性体的探测需求日益增长。以往在磁力异常的校正计算、数据预处理、正演、反演及定量解释中的许多假设条件和近似条件逐渐变得不再合理，在传统算法中被忽略的系列因素也逐步纳入考虑范围，例如自退磁效应[5]、磁化率各向异性[6]、剩余磁化强度[7]、总磁场强度异常的非线性效应[8-17]等。

总磁场强度异常(ΔT)的非线性效应，指的是实际测量得到的磁异常ΔT为地磁场总场(T)的模量|T|与主磁场(T0)的模量|T0|之差(即ΔT=|T|-|T0|)，与磁场矢量异常(Ta)以及岩石圈磁源的磁矩不满足线性关系。在实际数据预处理、正反演和定量解释时，通常将Ta在主磁场或正常地磁场T0方向上的投影分量ΔTPro近似为ΔT[10]，这就意味着ΔT与其他磁分量异常一样，与场源的磁矩呈简单的线性关系，与主磁场T0无耦合关系，并且在场源外部空间满足拉普拉斯方程从而具备调和性质，因而可以简化实测ΔT数据的预处理、正反演及定量解释过程[8-9]。……