地球重力场,大地测量.docVIP

1位理论基础 地球重力场反映了地球物质的空间分布及地球的旋转运动，它不仅决定了地球的形状和大小，而且反映了地球表面、内部以及大气和海洋的物质分布、运动和变化引力位函数对任意方向的导数等于引力在该方向上的分力；引力位是一个在无穷远处的正则函数；质面引力位是连续的、有限的和唯一的，而其一阶导数在经过层面时是不连续的；质体的引力位及其一阶导数是处处连续的、有限的和唯一的，而其二阶导数在密度发生突变时是不连续的；质体引力位在吸引质量外部满足拉普拉斯方程；质体引力位在质体内部满足泊松方程必须知道地球内部各点的密度，而后进行体积分。只要知道了水准面σ上的重力值，就可以计算地球外部任意一点的重力位或引力位这正是解引力位边值问题的理论基础。第一类边值问题——已知边界上的调和函数值的解唯一，第二类边值问题——已知边界上的调和函数的导数值的解唯一，第三类边值问题——已知边界上的调和函数与调和函数导数的线性组合的值，如果线性组合的系数异号，那么解唯一。借助于球坐标系，利用分离变量解的方法确定了拉普拉斯方程的通解的形式，其中含有勒让德多项式，我们顺带研究了勒让德多式的性质和递推公式。根据质体引力位公式，将引力位展开成球函数级数展开式，并研究了低阶位系数的物理意义：零阶位二阶系数与转动惯量和惯性积分有关。最后从重力等位面的概念及其曲率的定义以及铅垂线的曲率的推导中，我们可以得出：水准面是不平行的、铅垂线是一条空间曲线。 地球外部重力位和重力完全由形状已知的重力位水准面σ、质量M、角速度ω确定然后我们就根据水准椭球和引力的性质来推导正常位的解析表达式。 第一步，在椭球坐标系下，得到调和函数（水准椭球外部正常引力位）的椭球谐函数展开式。 第二步，由于水准椭球的旋转对称性，水准椭球外部正常引力位与经度无关，进行化简。 第三步，结合离心力位的椭球坐标表示令u=b得到关于水准椭球上的正常重力位的方程，而且水准椭球是重力等位面，那么这个方程就与纬度无关，便可将无穷级数的个数缩简为有限个数。第四步，利用数学关系式，以及地球质量和水准椭球上的正常重力位的关系，可以得到水准椭球上的正常重力和正常重力位，水准椭球外部正常引力位的椭球谐函数展开式。 仿照引力位的球函数展开式，可以将正常引力位展开成球函数，结合水准椭球的性质（旋转对称）以及椭球谐函数展开式（两式在极点的坐标、引力位相等），对其进行化简，得到可以进行应用的正常引力位的球谐函数展开式，以及系数的计算式。 由于扰动位的边值问题的解算中经常会用到大地水准面上的正常位，我们就需要根据水准椭球面上正常重力线的曲率得到正常重力向上延拓的方向和数值改正。 斯托克斯边值问题，又称大地测量边值问题，它将地表重力观测值归算到大地水准面上，计算出重力异常，然后假定大地水准面外部没有质量，解算大地水准面外部的扰动位，从而确定地球外部重力场。地球周边空间的每一点都有确定的正常重力位定义扰动位为同一点上的真实重力位与正常力位之差。因为正常力位是已知的，所以如果确定了扰动位，重力位也就确定了。另一方面，由于我们规定平均椭球体的质心与地球的质心重合，自转角速度与地球的自转角速度相等，所以正常重力位中的离心力位部分就等于地球重力位中的离心力位，因此扰动位实际上等于两者引力位之差。定义扰动重力为同一点上真实重力与正常重力之差，等于扰动位的梯度。 Dirichlet（第一）边值问题：在陆地用GPS水准、海洋通过卫星测高，得到大地水准面差距，进而根据布隆斯公式得到大地水准面上的扰动位，假定大地水准面外没有质量，则大地水准面外的扰动位就由唯一的满足在大地水准面上等于扰动位的调和函数确定（调和函数第一边值问题的解的唯一性）。 Neumann（第二）边值问题：重力测量得到地面重力g，GPS测定地面重力观测点的地心坐标，转换为大地坐标，水准测量获得地面重力观测点的正高，由大地高和正高可得到大地水准面差距。根据正高将地面重力归算到大地水准面上，根据大地水准面差距将水准椭球面上的正常重力归算成大地水准面上的正常重力，进而求得大地水准面上的扰动重力。假定大地水准面外没有质量，则大地水准面外的扰动位就由唯一的满足在大地水准面上的数值等于扰动重力的调和函数确定（调和函数第二边值问题解的唯一性） Stokes（第三）边值问题：由重力测量基本微分方程可知，重力异常是扰动位在水准椭球上的值与扰动重力在大地水准面上的值的线性组合，假定大地水准面外没有质量，则大地水准面外的扰动位就可以由满足在大地水准面上的函数值与其导数值的线性组合等于重力异常，而且具有相同线性系数的调和函数来描述（调和函数的第三边值问题）。 扰动位的三类边值问题的球近似解：因扰动位是小量，当大地水准面界面作球近似时，所引