地球重力模型-武汉大学课程中心.ppt

地球重力学 武汉大学地球物理系 第七章 重力模型的确定与应用 地球重力模型 地球重力模型及其作用　地球重力场在地球物理学、海洋学和空间技术中占有特别重要的地位。它直接反映地球内部的密度分布。从地幔产生的长波信号，到大陆岩石圈和海底地壳的局部特征等，都反映在地球重力场中。用一组重力位系数来表示相应的地球重力场，称为地球重力场模型。它的作用可简单归纳为以下几点：　　①卫星大地测量定位的精度取决于卫星定轨的精度，而全球重力场模型是精密定轨的基础。　　②通过地球重力场模型及对地球外部重力场的分析，可为地球物理学和地质学提供地球内部结构和状态的信息。　　③地球重力场模型可精确确定地球的扁率。　　④各国的区域性坐标系与全球坐标系的精确转换，需要区域性大地水准面资料，而大地水准面属地球重力场的一个等位面。　　⑤大地测量观测是在地球重力场内进行的，数据的处理和归算要知道地球重力场。　　⑥人造卫星、洲际导弹轨道的摄动与地球外部重力场密切相关。　　⑦重力勘探是重力学原理在勘探地下资源方面的应用，根据局部重力场变化规律可以反推矿藏位置和范围。 地球重力模型的发展发展历程 地球重力模型 …………… 1.重力场展开与快速计算 重力场的球谐函数展开 函数的勒让德展开 函数的勒让德展开 重力场的球谐函数展开 说明：①系数由地球重力模型给定，②求和从2开始是满足“质量相等、质心重合”，③尽管地球模型的不同，但求和教师有限的，④面球谐函数的计算稍后讨论。 勒让德函数的计算方法 勒让德函数的计算 勒让德函数的计算 正向横推法 其基本思路(见图)是利用固定列的勒让德函数递推公式计算较高阶的勒让德函数。求出部份0,1，2阶项，先求出纵向的0,1列（红线），其公式为 其中： 自第3列开始，依次用下式从左向右递推（称为正向横推方法）： 说明：后式中①若n=k→∞时首项收敛末项为0，②若k固定n→∞时首项为0末项系数为－1，③n固定k→∞不可能（k≤n）。 勒让德函数的计算 正向纵推法 正向纵推方法是先求得对角线上元素，再向下求出其它的勒让德函数(见图)。①先推导出对角线上的函数的公式：②导得利用对角线上的函数值推导出紧挨对角线下方(同列)一函数的值，即 ③再用列递推公式求得每列余下的函数值：  说明：后式中①若k固定n→∞时首项系数为2x末项系数为-1，②若n→∞且k→∞时首项系数为0末项系数为－1。（若2x1方收敛） 勒让德函数的计算 反向横推 该方法是正向纵推方法的改进(见右图)。①先推导出对角线上的函数的公式：②导得利用对角线上的函数值推导出紧挨对角线下方(同列)一函数的值，即 ③再用列递推公式求得每列余下的函数值：  说明：后式中①若k固定n→∞时首项系数为0末项系数为-1，②若n→∞且k→∞时首项系数为0末项系数为－1。 展开级数的快速计算 级数的快速计算 常规计算方法　首先，在实际计算过程当中，大地水准面高的计算公式中的地球位系数的个数不可能是无穷多个，它总是一个有限值，即写为在上式中我们不难发现，计算某上点的大地水准面高需要多次计算三角函数，可以证明需计算三角函数计算量十分庞大。为避免这种重复，交换求和的顺序，即 级数的快速计算 基本公式与分析　1979年,利兹（Rizos）描述了二维均匀间距的地理网格上用位系数计算大地水准面高或重力异常的计算机技术，常称为利兹法。他编制的程序,是为一个由经纬度所限定的区域中,计算已确定了经纬度间隔的网格点值而设计的。这个计算区,可以是局部区域,也可以是全球。大地水准面的高的计算公式为： 级数的快速计算方法之一：格网分裂法 级数的快速计算 级数的快速计算方法之一：格网分裂法 经度的分裂　将经度按等间距写为　　　　　，　是格网最左列的经度原点。由于：　　　 其中Q仅与纬度θ相关。 方法评述　方法的优点是①先由径度原点计算每一纬度点上的系数Q，再代入公式计算大地水准面高。这样对纬度而言，没有重复计算。②利用递推性质计算三角函数。具体公式由读者自选导出。 计算速度　结果表明，利兹方法比常规方法快近100倍。可见利兹方法比常规方法具有更高的效率。但必须指出，利兹方法只适用于计算等间距地理网格点的重力量，若计算单点的重力量，则这两种算法的计算速度将不会有太大的差别。 级数的快速计算 基本原理　将公式　　　　　　　　　　　　　　　　　　写成为如下形式 级数的快速计算方法之二：矩阵求和法 化简：将勒让德函数的（纵向）递推性质写为矩阵形式有 矩阵形式 方法评述　该方法比较稳定，减少了有效数字的损失。用该方法计算局部区域或全球网格点的重力量时，比常规方法快20倍。该方法的另一个优点是还可以计算重力场的其它量。 级数的快速计算 级数的快速计算方法之二：矩