邱隆君，孙诚业*，杨亚斌，吴新刚

（1. 中国地质科学院地球物理地球化学勘查研究所，河北廊坊 065000；2. 国家现代地质勘查工程技术研究中心，河北廊坊 065000）

0 引言

重力勘探是一种传统的地球物理勘探方法，是对实测重力数据进行自由空气校正、中间层校正、地形校正、正常场校正等预处理获得布格重力异常［1-2］，再用场分离技术和重力反演方法获得地下密度异常体的位置信息或者密度信息［3-4］的一种间接勘探方法。

在数据预处理阶段，为了消除起伏地形密度体的干扰，需要进行地形改正。它的计算过程主要考虑两个方面：第一，对起伏地形选择恰当的密度数值。既可选择能代表研究区地层特征的密度值，也可以按照有关规范设定岩石密度、海水密度、淡水密度［5-6］。第二，结合实际的地形改正范围选择合适的地形数据库。一般近区和中区地形改正范围为0～2000 m，采用空间分辨率和精度较高的地形数据库，如应用高精度差分卫星定位测地技术（GPS-RTK）或航摄获得的高精度数字高程模型（DEM）数据［7］；远区地形改正范围为2000～166700 m，需要覆盖率较高的地形数据库，如航天飞机雷达地形测绘使命（SRTM）或高级星载热发射和反射辐射仪（ASTER）得到的DEM 数据。

地形数据库的空间分辨率、精度、覆盖率都会影响计算结果，分析某单一变量难以清晰地论述这些关系。因此，地形数据库是全面、深入地研究地形改正的各种影响因素及其影响程度的关键要素。全站仪或GPS-RTK 方式获取的高程数据精度最高，但成本高，工作效率低；……