简介： 本发明公开了，包括以下步骤：步骤1、将三维地质体剖分成多个四面体单元；分别对观测点和三维地质体建立对应的八叉树结构；步骤2、将三维地质体划分为近区与远区，步骤3、分别计算近区与远区重力参数；基于近区所包含的四面体单元的重力积分无奇异性解析表达式直接计算近区的积分，得到近区重力参数；并基于观测点八叉树和源八叉树，采用自适应快速多极展开算法计算远区的积分，得到远区重力参数；步骤4、将步骤3中计算得到的近区重力参数与远区重力参数求和，得到观测点重力参数。本发明能够基于非结构化网格剖分，对任意分布的观测点进行重力正演计算，大幅度地提高了重力正演计算的效率。