第四章——地球正常重力场.doc

第四章 地球的正常重力场 重力测量结果表明，地球在其表面上的重力分布是有规律的；总的说来，它由赤道向两极逐渐增加，由赤道上的978Gal逐渐增加到两极的983Gal。在大地测量中，参数合适的旋转椭球是地面点坐标的参考架，当参考椭球选定后，大地水准面相对参考椭球面的起伏不超过110m，起伏只占参考椭球赤道半径的2×10-6。因而自然想到，用质量等于地球总质量、以地球自转角速度绕其极半径旋转的旋转椭球来模拟真实地球，用这种地球模型（正常场地球模型），在其表面上和外部空间产生的重力场称为地球的正常重力场。当正常场地球模型在地球内部定位后，地球的重力场可以分解为两部分，一部分是正常场地球模型在该点产生的重力场，第二部分为真实地球与正常场地球模型的密度分布不同在该点产生的重力场；前者称为地球在该点产生的正常重力场，后者称为地球在该点产生的重力异常场。重力测量结果表明，当正常场地球模型选择合适后，大地水准面上的重力异常场不超过150 mGal，约占地球正常重力场的1×10-4~2×10-4。地球的重力异常场虽只占地球重力场的万分之一二，但它却包含了有关地球内部结构和大地水准面形状的重要信息，因而研究地球重力异常场空间分布规律以及它们与地球内部结构和大地水准面形状之间的关系已成为重力测量的重要目的之一。根据第三章的结果，本章给出正常场地球模型在旋转椭球面上产生的重力、正常重力位二次导数张量以及它在其外部空间产生的大地位球函数展开系数。 4.1 旋转椭球的几何参数 引入笛卡尔直角坐标系，坐标原点置于旋转椭球的中心，沿其极半径，在其赤道平面内，则旋转椭球面的方程为 其子午椭圆的方程为 其中、分别为旋转椭球的赤道半径和极半径，它们是决定旋转椭球形状的两个几何参数。考虑到参考椭球的赤道半径和极半径相差很小，其扁率约为3×10-3量级，因而参考椭球的子午椭圆与圆非常接近，为了讨论问题方便，对子午椭圆常引入下面几个几何参数： 子午椭圆的扁率、第一偏心率、第二偏心率有下述关系 如图4.1.1所示，与轴之间的角度为点的地心纬度，点子午椭圆的法线与轴之间的角度称为点的大地纬度，因为子午椭圆与圆非常接近，点的地心纬度和大地纬度相差很小，其差约为子午椭圆扁率的量级。在图4.1.1中，有 根据（4.1.2）式，有 因而有 将上式代入（4.1.5）式，得 大地纬度和地心纬度相差很小，根据（3.2-1.6）式可以求出它们之间的相互换算关系，与（4.1.6）式相对应的（3.2-1.6）式中的、分别为 因而有 考虑到子午椭圆的扁率约为为3×10-3量级，有时将子午椭圆的方程写成极坐标的方式比较方便。将（4.1.6）式代入（4.1.2）式，把子午椭圆的直角坐标方程（4.1.2）写成极坐标的形式，考虑到，有 因为 将（4.1.11）式代入（4.1.10）式，化简，舍去高于的项，即舍去小于30×10-9的项，得 在（4.1.12）式中用大地纬度代替地心纬度，根据（4.1.9）式，舍去含高于的项，即舍去小于30×10-9的项，有 参考椭球面上大地纬度为的子午椭圆的曲率半径和卯酉圈的曲率半径的数学表达式分别为 4.2 索米格兰纳（Somigliana）正常重力公式 正常场地球模型在其表面上产生的重力矢量是正常重力位在此表面上的梯度，考虑到旋转椭球面是正常场地球模型的一个重力等位面，因而正常场地球模型在其表面上产生的重力矢量应垂直于旋转椭球面，亦即 其中、分别为坐标的单位坐标基矢量和它的拉梅系数，根据（3.2.9）式，旋转椭球面上上的拉梅系数为 其中、分别为改化余纬和改化纬度。习惯上正常重力矢量的方向约定为，即约定它指向旋转椭球内部为正，则参考椭球在其表面上产生的正常重力等于正常重力矢量在方向上的投影，即 （4.2.3）式称为正常重力公式。在旋转椭球体的重力位表达式（3.3.25）式中用改化纬度代替改化余纬，考虑到为的二阶勒让德多项式，可以把（3.3.25）式写成 根据（3.3.25）式，有 将（4.2.2）、（4.2.4）式代入（4.2.3）式，得 用、分别表示赤道上和两级的正常重力，根据（4.2.6）式，有 将（4.2.7）式代入（4.2.6）式，得 根据（3.2-1.3）式、（4.1.7）式，可以求出大地纬度和改化纬度之间的关系，它们是 将（4.2.9）式代入（4.2.8）式，化简得到以带地纬度为变量的正常重力公式，它为 （4.2.10）式是意大利人索米格兰纳于1929年导出的，它称为索米格兰纳正常重力公式。 4.3展成级数形式的正常重力公式，克雷诺（Clairaut）定理 斯托克斯定理表明，正常场地球模型的赤道半径、扁率、总质量和旋转角速度唯一地决定了旋转椭球在其表面上和外部空间产生的重力场。正常重力公式（4.2.10）式给出了以大地纬度为变量、以赤道上和两极重力为