论文部分内容阅读
本文进行了强剩磁条件下的三维磁异常反演研究,通过同时反演磁化率与剩磁矢量丰富磁测反演结果、提高磁异常解释的精度。首先,采用非结构化网格有限单元法进行同时考虑磁导率与剩磁矢量的正演,非结构化网格有限元算法更适用于复杂模型计算,同时具有较高的计算速度,为反演提供了保障。在有限元正演的基础上,分别开展了磁导率反演、同时反演磁导率与剩磁矢量的反演,建立了三维磁测数据反演目标函数,通过高斯-牛顿法进行目标函数迭代求解;为减少内存的使用,采用了有限内存技术,通过正演计算灵敏度矩阵及其转置与向量的乘积,进而避免灵敏度矩与海森矩阵的存储。最后通过理论数据与实测数据反演,验证了本文算法的有效性。