有限元矩量法三维位场曲面延拓

摘要

针对三维曲面延拓在高精度数据处理与解释中的重要性和实现的困难性,本文提出了有限元矩量法(FEMM)曲面延拓算法。首先采用等参元有限单元剖分最小二乘矩量法将延拓方程离散化成大型实满矛盾方程组,然后用共轭梯度法求解。这种算法可直接利用起伏观测面上的不规则网数据,既可作向上延拓,也可以完成向下延拓曲化平。其精度、分辨率和保真性大大高于现有的其它方法。按照本算法编制的程序使用方便,能适应各种复杂观测条件,在各类计算机上运行。对模型数据与野外资料的处理结果,证明了这种算法的有效性与实用性。