本发明涉及一种基于泛克里金法的重力场插值方法，其主要技术特点是包括以下步骤：针对空间中的重力异常模型建立回归模型；针对空间中的重力异常模型建立相关模型；计算回归模型和相关模型的权重系数，代入重力异常模型，便可获得待插值位置上的重力异常空间分布表达式，再将待插值位置代入重力异常模型，便可计算出待插值位置的重力值。本发明采用逐步回归法对重力场的非平稳性及相关性分别建模描述，使用多元逐步回归法，结合拉格朗日乘数法，完整的描述了重力场的非平稳性，保证了插值结果的无偏性，具有较高的插值精度。