纯球坐标系内各项重力校正值计算方案和计算过程.pdfVIP

第34卷第6期 物 探 与 化 探 Vo1．34．No．6 2010年 12月 GEOPHYSICAL ＆GEOCHEMICALEXPLORATION Dee．．2010 纯球坐标系 内各项重力校正值计算方案和计算过程 安玉林 ，张明华 ，黄金明 ，乔计花 (1．中国地质大学 地球物理-9信息技术学院，北京 100083；2．中国地质调查局 发展研究中心，北京 looo37) 摘 要：在国内，首次阐述了 “纯球坐标系内各项重力校正值的计算方案和计算过程”。该计算方案有如下优点：① 以观测点P为中tk，以球面距离o =n +Aa·b一 (m=2，3，…，)为半径将计算区划分为 个环带，再以△A。 为由P点出发射线的角度问隔将所有环带等分为Ⅳ(整数)份的计算区划分方案是最严密拼接，因而校正值的计算 精度相当高。②在地形复杂和相对高差大的情况下 ，通过减小 Aa、b和 △A。数值，可进一步提高计算精度。③可以 直接采用带有经、纬度坐标的高程数据，进行大面积测区的各项重力校正值计算，而不受6。带和标准图幅限制。 关键词：球坐标系；重力校正值；计算方案和计算过程 中图分类号：P631．1 文献标识码：A 文章编号：1000—8918(2010)06—0697—09 本论文中，涉及笔者 2002～2008年长约7年的 则性区别 。 研究工作o。文 中重要公式是笔者独立地、严谨地 1．1 校正值计算的积分表达式 推导出来的，在我国 《重力勘探》文献 (如教科书和 球坐标系内陆地重力观测值校正，一般情况下 历年科技论文)中难于见到，故全部列举 出来。在 包括正常场校正、高度校正、中间层校正、地形校正 国内(据说 同样在 国外)，作者至今还未见到关于 和均衡校正。对重力观测值完成所有校正之后，得 “纯球坐标系内各项重力校正值计算方案和计算过 到的就是 “均衡重力异常”。 程”的论文发表。实用的计算方案依然是 “传统式 球坐标系内重力地形校正和中间层校正，就是 的”(图1左)，在相对起伏大的地形条件下误差大。 要消除水准面到观测面之间的质量在测点处产生的 因此，作者认为有必要把 自己研究成果发表 出 重力值。重力高度校正与正常校正合二为一，则是 来，目的一是要让这些公式及其应用与国内从事重 在观测点完成正常场校正。假设地球上地幔到地表 力勘探教学和科研工作 的同行见面 ；二是希望填补 观测面之间的密度均匀 ，那么在上述 4项校正完成 以后的重力教科书的空 白；三是希望按照 “纯球坐 后，测点上的异常重力值应该等于0。然而大量重 标系内各项重力校正值计算方案”设计的计算机程 力测量工程实践表明，经上述4项校正后的重力值 序能在重力野外测量中得到推广。 中，明显地存在一个与观测点高程呈负相关关系的 笔者详细阐述的计算过程和程序内的参数说明 大区域负值背景异常。这就告诉重力学者们，在地 合起来就是程序使用手册。 壳中一定存在相当大的质量亏损，使原始重力观测 值大大降低。于是普拉特和艾黎先后提出了两种均 1 纯球坐标系内重力各项校正值计算公式 衡补偿理论，并提出对重力值还应进行均衡校正。 “纯球坐标系内重力各项校正值计算 ”是指 ，所 深度起伏变化的莫霍面的存在 ，使艾黎理论被大多 有计算公式均是 以(r，A， 为坐标的球坐标系公 数重力解释人员所接受。于是，基于艾黎理论的均 式，其计算方案如后文图1中右图所示。 “纯”字在 衡校正方法，就成了消除大区