秦 熠，杨宇山，刘天佑，朱 丹

(1.中国地质大学 地球物理与空间信息学院，湖北 武汉 430074；2.四川省煤田地质局 地质工程勘察设计研究院，四川 成都 610000)

1 引 言

视磁化强度填图方法是一种利用磁异常反演地下地质体磁性强弱的方法，利用航磁资料进行视磁化强度填图可以有效确定磁性体边界，划分磁性地层，特别对于存在大规模隐伏岩体地区具有良好的应用效果，作为研究区域构造的辅助手段，已被广泛地应用于地质构造解释与研究之中[1]。

1973年Parker[2]将傅里叶变化引入位场正演计算，提出了连续单界面模型快速正演方法，即Parker法。随后，Oldenburg[3](1974)在Parker法基础上推导出了频率域密度界面迭代反演方法，由于该方法解决了过去空间域离散模型界面反演方法计算速度慢的问题，因此被广泛地应用于起伏界面的计算。王万银等[4](1993)在单界面重力场正反演基础之上研究了更加符合实际情况的双界面重力场正反演问题，并提出了保证反演快速收敛的措施。

视磁化强度填图兴起于20世纪80年代，Pilkington[5](1989)根据Parker界面正反演方法推导出了磁化强度反演公式，并应用于Athabasca盆地。早期视磁化强度填图，尤其在化极处理之后存在负值的问题，难以解释。管志宁等[6](1990)研究了视磁化强度填图出现负值的原因，并提出一种采用立体角滑动窗口迭代反演方法的视磁化强度填图，应用于秦巴地区。这种方法适用于地形起伏不大，观测平面水平的航磁资料处理。张贵宾等[7](1996)从理论与模型算例两方面对Parker公式为基础的频率域位场反演方法做了深入研究，提出新的快速正则化反演方法，并将其推广为单界面无限延深视密度和视磁化强度填图方法。……