李金朋 范红波 刘 利 张英堂 李志宁 武炳阳

(①93114部队，北京100195；②陆军工程大学石家庄校区，河北石家庄050003；③北京市遥感信息研究所，北京 100192)

0 引言

磁测数据反演主要应用于军事侦察(未爆弹、潜艇和水雷等)、水文、能源及工程地质勘探等领域[1-4]。磁性目标体的三维反演是利用观测面上磁测数据，对其空间形状、位置、磁化强度及磁化方向等参数进行求解的过程，为下一步地质解释提供可靠的数据支持。

磁性体反演主要包括物性反演和形态反演两方面[5]。物性反演的主要思路是将观测面对应的地下区域离散化为规则的长方体网格，通过反演获得这些网格的磁性参数，并最终获得目标体的磁性参数空间分布[6-8]。物性反演过程中，由于目标函数固有的欠定性特征，计算结果存在多解性；同时，由于需要进行多次模型正演计算，核函数矩阵会大大增加计算时间，影响计算效率。形态反演的主要思路是利用一定的计算准则、基于观测数据对地下多面体进行拟合，利用多面体的形态对待测模型的三维分布进行求解[9]。Uieda等[10]利用L2范数定义数据泛函约束函数对待测目标进行形态反演；曹书锦等[11]基于L1范数利用种子反演方法，提高了重力梯度张量反演精度；Uieda等[12]利用重力梯度张量，对不同待测目标的“种子”分配不同密度值，实现特定目标的反演，有效排除了非目标场源的干扰。常规的形态反演方法能够克服物性反演方法在迭代过程中对核函数矩阵进行多次计算以及计算结果存在多解性的问题，但前提是需要获得待测目标的先验信息。……