刘卓, 曾昭发, 李静，何荣钦，霍祉君，林景宜

1.吉林大学 地球探测科学与技术学院，长春 130026;2.国土资源部 应用地球物理重点实验室，长春 130026

0 引言

居里面是地壳中岩石铁磁性矿物在温度升高到居里点变为顺磁性物质时的深度界面，为磁法勘探的磁性层底界面，表征地下温度场的分布，对地热田评价开发、油气资源的预测、地震火山灾害防治及原生热液矿产勘查等具有重要指导意义[1]。居里面深度计算方法包括谱分析方法和界面反演方法等。谱分析方法[2--3]是通过计算磁异常数据对数功率谱来确定磁源深度，但对数功率谱的复杂性和混乱性难以有效地划分磁源平均深度，具有较大的误差。Parker在位场界面正演计算中引入快速傅里叶变换(FFT)[4]。Oldenburg根据Parker公式，提出频率域深度界面迭代反演方法[5]，改变了过去采用棱柱或长方体模型反演界面深度不连续性，同时计算速度快[6]。Parker--Oldenburg反演算法存在3个问题：①该方法属于单磁性界面反演方法，磁性体通常是由双界面或多界面组成的,利用单界面计算与实际情况不符。王万银[7]和相鹏[8]采用双磁性界面模型来提高反演结果精度。②在反演过程中假设磁化率参数为常数，而实际地质介质磁化率参数在横向和纵向上均具有较复杂的变化。冯娟等[9]采用了垂向变物性参数模型来计算基底界面，提高了反演精度。③在Parker--Oldenburg公式中包含向下延拓因子中的不收敛项，影响计算的稳定性。尽管可以通过添加滤波器来确保迭代的收敛，但将会使信号的高频信息损失，降低反演结果的分辨率和精度。……