两种新的长方体重力异常正演公式与其理论推导.pdfVIP

两种新的长方体重力异常正演公式及其理论推导 骆遥1, 2 1 100029 中国科学院地质与地球物理研究所，北京 （ ） 2 100049 中国科学院研究生院，北京 （ ） E-mail：geo@ 摘 要： 在前人推导长方体重力场、磁场正演理论表达式工作的基础上，重新对长方体重 力场正演理论表达式进行理论推导，提出了两种全新的长方体重力异常正演公式形式，并给 出了全部的理论推导过程，对比模型正演计算结果表明，新导出长方体重力场正演理论表达 式的正确。 关键词：长方体，重力场，正演，积分 中图分类号：P631 1. 引 言 在地球物理勘探中，重力异常的正演计算在重磁异常的解释中有着重要的意义，它反映 场源与场分布特征之间的联系，是地质解释的基础。长方体是位场正演计算中最为常用的三 度体模型，可以用平行于直角坐标面的平面对任意形状的三维地质体进行分割，将三维复杂 地质形体重力异常的正演计算转化为一系列长方体体元重力异常的迭加；所以，长方体重力 异常的正演计算是重力勘探中最基本的问题。同时，随着高精度重磁资料处理解释工作的开 展，二维重磁反演已难以满足全方位确定地质构造精细展布的要求，重磁反演逐渐发展到三 维；重磁三维反演分为形态模型反演和物性模型反演，限于形态反演中存在的较大困难，近 年来物性反演成为三维重磁反演中的主要形式[1~6] 。重磁三维物性反演中要将地质场源区域 离散化成离散的长方体单元，长方体重力异常的正演计算成为重磁三维物性反演中最为基本 的问题。因此，研究长方体重力异常的正演是一项基础性工作，具有重要的理论和实际意义。 2. 长方体重力异常 P ( , , ) P (x, y , z ) 位于空间中 点 具有密度 的质量单元在 点 处产生的重力场为： z g = G v 1 ， （） R 3 R P P v G V 式中 是点 到 的距离， 质量单元的体积， 是万有引力常数。通过对体积 的积 分可得到均匀密度体的重力场： z g (x, y , z ) = G dv 2 （） R3 V 2 2 2 1/ 2 其中，R = {( x ) + ( y ) + ( z ) } 。 1 于是，如图 所示的均匀密度长方体产生的重力场可以表示为： 2 2 2 z g (x, y , z ) = G