基于L-BFGS反演算法的ΔT精确计算磁异常分量T_(ap)方法

磁法勘探理论中,将ΔT磁异常看作磁异常矢量Ta在地磁场方向的分量T_(ap),是ΔT异常处理与解释的物理基础,然而这种近似存在误差,理论计算及实验已经证明这种近似所产生的误差将随着Ta异常强度的增大而迅速增加。当磁异常Ta远小于地磁场T0时,误差影响小,可忽略,在强磁异常情况下,误差大,ΔT异常的处理解释精度会受到很大的影响。对于高精度磁法勘探而言,必须将ΔT转换成磁异常分量T_(ap)进行处理解释。笔者提出了基于有限储存BFGS(L-BFGS)反演算法的ΔT精确计算磁异常分量方法,首先推导了T_(ap)计算ΔT的正演公式,利用ΔT与T_(ap)的差值构建反演T_(ap)的目标函数,采用L-BFGS算法由ΔT解算T_(ap)。模型实验表明该方法计算得到的T_(ap)十分接近理论值,即可将误差降低两个数量级。在存在噪声与背景场情况下该方法也都能得到很好的结果。将本方法应用于福建阳山铁矿ΔT磁测资料的处理,得到了与实际更加符合的处理解释结果。