基于非结构化单元的磁场及梯度模拟计算

在地球物理磁法勘探中,磁场及其梯度张量数据是磁法正反演及解释的重要基础。本文结合磁梯度张量解析解表达式,实现了一种基于非结构化四面体单元的磁场及梯度模拟计算方法,该方法弥补了传统方法难以处理任意复杂模型的缺陷,同时大幅提高了模拟计算精度。首先,采用非结构化单元剖分异常体可以显著地减小几何误差,更好地模拟带地形复杂模型。其次,采用新的磁梯度张量解析解表达式,其精度比传统的有限元、有限差分方法更高。本文设计了单个球体模型和组合体模型,模拟计算结果体现了该方法的适用性和准确性。