重力近、中去地形改正的计算机实现

重力近中去地形改正的计算机实现

-以兴城实习基地为例

【摘要】目前区域重力测量工作正在全国范围内普遍开展。我国山区多，地形相当复杂，重力测量的地形影响不可忽略。重力地形改正值的计算是重力勘探中一项非常繁重的工作，本文用数值积分法，编制了方域地形改正程序。它能用于各种重力工作的地形改正，并能满足重力地形改正的精度要求，通过人机对话的形式，可自选择改正条件，方便灵活。

区域重力调查地形改正公式常规用锥形体公式和扇形体公式计算，利用手工的方法进行计算，在读取高程数据时认为因素读的不准。

利用量板手算地形改正的老办法，不论测网是规则还是不规则的都能适用，但是其工作量之大，效率之低根本无法适应区域重力工作的需要。分析其原因，一是计算方

法本身存在缺陷，采用圆域“电子量板”虽然在测点

近区使用较小的扇形柱能比较准确地模拟地形起

伏，但在测点远区，扇形柱底面积较大，有限的插值

点难以精确估算出扇形柱的平均高程，模拟地形变

化的精度当然不高；二是在地形变化较大的区域，特

别是在一个扇形柱中有狭长的沟谷或狭长的山脊

时，有限的固定插值点的高程难以代表该扇形柱的

平均高程。

多年来，重力中区(50或100—2 000 m)地形改

正一直是重力资料整理过程中的一项既费时又费力

的工作，一个技术人员每天也就能完成2O个点左右

的数图工作，且错误率高，需进行一定数量的抽查、

重复数图和对算等工作，效率和精度都很低。目前，

由于《区域重力调查规范(DZ／T0082-2006)》的颁布

实施，1：20万区域重力调查布格异常总精度提高

了一倍，重力中区地形改正的精度要求由±0．44

×10 m／s 提高到4-0．20×10 m／s ，这迫切需要

使用新的中区地改的方法和技术。

2 工作流程

(1)数字地形图的准备，最好是直接购买1：5

万数字地形图或DEM高程模型数据，借助GIS软件

进行格式转换。在没有1：5万数字地形图的情况

下，将1：5万纸介质地形图扫描，采用MapGIS或

ArcGIS将等高线和高程注记点等信息数字化，并赋

予高程属性数据。应尽量使数字地形图全面地反映

原地形图上与高程有关的信息。

(2)将等高线转换为三维TIN模型数据，并用

高程注记点等高程信息编辑修改。

(3)将TIN数据转换为一定间距的网格数据，

本次计算中采用了20、50、100 m三种网格数据。

(4)计算各种网格间距的中区地形改正值等。