地球扰动引力场的计算.docVIP

地球扰动引（重）力场的计算方法 一、研究扰动引（重）力场的意义 物理大地测量的主要目标是确定地球外部重力场和大地水准面这两项任务紧密相关是一个问题的两个方面 二、扰动引（重）力场的定义 由于地球形状的不规则和内部质量分布的不均匀，地球重力场是很复杂的。它可以分解为两部分。一部分是规则的，可以看作是一个形状和质量分布都很规则的质体所产生的重力场，而且它占了地球重力场的绝大部分，通常称之为正常重力场；剩余的微小部分极不规则，称之为扰动重力场。因正常重力场是己知的，只要求得了扰动重力场，加上己知的正常重力场,实际的地球重力场就可求得。 引入正常椭球作为地球形状及其重力场的一个近似体，它所产生的重力位称为正常重力位，记为U，相应的正常重力记为γ。真正的地球重力位W和正常重力位U之差就是扰动重力场中的重力位值，称为扰动位，记为T，即 扰动位在地球表面的外部空间中是调和函数，它满足Laplace方程 由地面上测量的重力与正常重力之差的数值提供了求解上面方程的边界条件。因为重 力与正常重力γ分别是重力位与正常重力位对法线反方向的偏导，即有 和分别为沿大地水准面的外法线方向（正高方向）和沿参考椭球面的外法线方向的偏导，上式对沿这两种方向的偏导做了近似相等的处理。通常称 为纯重力异常或扰动重力，g由重力测量得到γ则通过正常重力公式计算得到。 三、扰动重力与重力异常之间的关系 扰动重力与重力异常具有近乎相同的频谱特征，尤其在高频段，两者可以说是相等的，它们在各类问题的实际应用中能够发挥相同的作用。不过，扰动重力较重力异常更易于重力场计算。对比以扰动重力为输入量的Hotine积分和以重力异常为输入量的Stokes积分，前者的核函数较更简单，因此更便于计算。在平面近似时两者完全相同。目前在重力场计算中引入的各种快速算法（如FFT）也同时适用于以扰动重力为输入量的Hotine积分。此处不加证明，直接给出两者之间的关系式如下 或 上式中，Δζ为高程异常四地球扰动引力场的计算确定地球扰动引力场的基本思路是根据一种或多种重力数据作为边值条件建立关于扰动位的相应重力的大地测量边值问题Dirichlet边值问题Neumann边值问题边值问题传统的作业方法是通过测量地球表面固定点的重力值、水平坐标和高程得到重力异常来求解混合边值问题。随着三维定位技术的发展，一些新的重力数据类型，例如地球表面的扰动重力数据、航空重力数据以及卫星重力数据等相继出现，数据精度计算方法都有了很大的进步需要解决的边值问题类型也得以扩充。利用这些数据来求定高精度的地球外部重力场显得越来越重要。以下将对物理大地测量中涉及到的主要边值问题一一解释Neumann边值 故若计算得到扰动重力，便可得到第二类边值问题具体描述为 其中，第一式为泛定方程（Laplace方程），第二式为边界条件，它决定了边值问题的类型，第三式正则条件本质上也是一种边界条件（无穷远处为地球表面外部空间的另一边界），但该边界条件并不决定边值问题的类型。 扰动重力由前文公式 可计算得到，其中，地面点A的重力由重力测量得到，地面点的正常重力可由正常重力向上延拓计算得到 式中为地面点 式中为为单位球面为待定点和积分面元之间的球面角距具体形式为 Stokes问题 该问题的数学描述为 若将大地水准面近似为正球体，则Stokes问题的边界条件可以球近似为 式中，R为正球体的半径，表示沿径向的偏导 Stokes问题有多种解算方法，Stokes解出的结果为（1894） 式中为其 Molodensky问题 Stokes问题存在一些缺陷。首先是Stokes的边界面（大地水准面）是未知的，因此Stokes问题实际上是一个自由边值问题。其次，正高不能精确地确定。地面点到平均椭球体表面的高度（即大地高H）可分为大地水准面高N和正高h两部分。N通过解算Stokes问题得到，而h定义为地面点的位基数C与沿地面点的垂线方向至大地水准面这一路径上的平均重力gm的比值，无法精确确定。其次是只能通过水准测量近似得到。此外，在进行重力归算时要对重力梯度及地形质量的密度作假设，并且要移动质量，使得大地水准面外部不再存在质量。 因此，Molodensky以正常高为基础，提出一种以（似）地球表面为边界面的边值问题来研究地球的形状和外部扰动重力场。具体来说，Molodensky理论已知边界面上所有点的重力位（水准测量结合重力测量得到）和重力向量（重力测量得到），从而确定其表面及其外部位。 其数学描述为 上式中的边值条件已经利用了球近似，并忽略了正常重力长计算中的椭球扁率f项，这对解算高程异常ζ的影响可达cm级。例如，美国大地水准面的椭球改正范围为-5.6~3.1cm，精密解算则需要加椭球改正。 其解为 一般取前2项加和 式中，和分别为地球表面层的密度与万有引力常数的乘积的