1、大地测量学与测量工程专业优秀论文 应用低轨卫星跟踪数据反演地球重力场模型关键词：CHAMP 卫星 GRACE 卫星 地球重力场模型 动力学法 卫星跟踪数据反演 大地水准面 卫星重力探测摘要：低轨卫星跟踪测量技术(CHAMP 卫星和 GRACE 卫星)可提供覆盖全球的、高精度的空间重力观测资料：精密轨道数据、GRACE 卫星的 K 波段观测数据等，这些观测数据包含了丰富的地球重力场信息，尤其是 GRACE 卫星的 K 波段观测数据。在当代低轨卫星的定轨精度达到比较高的程度(约 2cm)，K 波段观测数据将是地球重力场信息提取的主要载体，如何充分利用此类信息进行反演地球重力场模型的研究是当代大地测
2、量研究的前沿和热点之一。其研究成果将对现代大地测量、地球物理、海洋科学以及地球动力学等学科的发展起到巨大的推动作用。 本文围绕应用低轨卫星跟踪数据反演地球重力场模型这一主题展开讨论，主要研究内容和创新点如下： (1)在查阅大量国内外相关文献的基础上，对卫星重力探测技术的发展历程、国内外研究现状以及研究意义进行了阐述。讨论了目前常用的几种地球重力场模型反演算法，并分析了它们的优缺点。 (2)详细研究了卫星跟踪卫星重力测量所涉及到的时间系统、坐标系统以及空间运行卫星所受到的各类摄动力模型，给出了卫星重力测量的基本原理、常用的算法以及各类偏导数的计算。定性分析了卫星初始位置误差对应用卫星精密轨道数据
3、反演地球重力场模型的影响，若不顾及其它测量误差(如加速度计测量误差)的影响，以 2 小时反演弧长为宜。 (3)详细讨论了 CHAMP 卫星和 GRACE 卫星相关观测数据的预处理方法，分析了姿态数据间断的三种情况，给出了利用二次插值法对间断的姿态数据进行填补，算例表明：在姿态数据间断 40 历元的情况下，二次插值法完全能满足精度要求。 (4)利用动力学法对 CHAMP 卫星的加速度计观测数据进行了重新标定，结果表明：加速度计的标定参数并不是一成不变的，一天标定结果的变化范围比 1 小时的要小，更接近 GFZ 给出的结果。加速度计重新标定后的重力场模型的精度要高于未标定的结果，因此对加速度计观测
4、数据进行重新标定是十分必要的。 (5)详细介绍了笔者独立自主开发的卫星重力数据处理软件-GRASTAR，该软件是利用 Visual C+语言平台开发完成的，可对 CHAMP 卫星和 GRACE 卫星观测数据进行处理，具有很好的可移植性和可拓展性。 (6)基于应用低轨卫星精密轨道数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始状态向量与地球重力场位系数一起求解的数学模型，利用 120 天 CHAMP 卫星精密轨道和加速度计观测数据解算出 50 阶次的地球重力场模型 TJCHAMP01S，检核结果表明：TJCHAMP01S 模型精度优于相同阶次的 EGM96 和
5、EIGEN-1S 模型，并且更接近 EIGEN-2 模型。在我国区域，TJCHAMP01S 和 EIGEN-2 模型与 EIGEN-CG01C 模型前 30 阶所表示的大地水准面差值的中误差分别为 4.4cm 和 3.3cm。 (7)基于动力学法，利用 26 天 GRACE 卫星精密轨道和加速度计观测数据解算出 36 阶次的地球重力场模型 TJGRACE01S-OR，检核结果表明：该模型的精度要优于相同阶次的 EGM96 和 EIGEN-2 模型，前 20 阶模型位系数与 EIGEN-CHAMP03S 精度相当。 (8)基于应用 GRACE 卫星星间距离变率数据反演地球重力场模型的动力学法，推
6、导了将加速度计观测数据的尺度和偏差以及卫星初始相对速度向量与地球重力场位系数一起求解的数学模型，利用 31 天 GRACE 卫星星间距离变率和加速度计观测数据解算出 60 阶次的地球重力场模型 TJGRACE01S，检核结果表明：TJGRACE01S 模型的精度要明显优于相同阶次的 EIGEN-CHAMP03S 模型，但要稍逊于 EIGEN-GRACE01S。利用 EGM96、EIGEN-CHAMP03S、TJGRACE01S 和EIGEN-GRACE01S 模型的前 60 阶分别计算全球 1#176;1#176;格网点的大地水准面高度，并与 EIGEN-GL04C 模型的大地水准面高度做比较
7、，相应差值的中误差分别为：0.446m、0.144m、0.027m 和 0.008m。 (9)利用 67 天GRACE 卫星星间距离变率和加速度计观测数据解算出 80 阶次的地球重力场模型TJGRACE02S，检核结果表明：TJGRACE02S 模型的精度要明显优于相同阶次的EIGEN-CHAMP03S 模型，但要稍逊于 GGM01S。利用 EGM96、EIGEN-CHAMP03S、GGM01S 和 TJGRACE02S 模型的前 72 阶分别计算全球2.5#176;2.5#176;格网点的大地水准面高度和重力异常，并与EIGEN-GL04C 模型的计算结果进行比较，大地水准面高度差值的中误差
8、分别为：0.524m、0.305m、0.021m 和 0.052m，而重力异常差值的中误差分别为：4.00mGal、3.07 mGal、0.21 mGal 和 0.50 mGal。 (10)对本文的工作进行总结，对未来的工作进行展望和规划。正文内容低轨卫星跟踪测量技术(CHAMP 卫星和 GRACE 卫星)可提供覆盖全球的、高精度的空间重力观测资料：精密轨道数据、GRACE 卫星的 K 波段观测数据等，这些观测数据包含了丰富的地球重力场信息，尤其是 GRACE 卫星的 K 波段观测数据。在当代低轨卫星的定轨精度达到比较高的程度(约 2cm)，K 波段观测数据将是地球重力场信息提取的主要载体，如
9、何充分利用此类信息进行反演地球重力场模型的研究是当代大地测量研究的前沿和热点之一。其研究成果将对现代大地测量、地球物理、海洋科学以及地球动力学等学科的发展起到巨大的推动作用。 本文围绕应用低轨卫星跟踪数据反演地球重力场模型这一主题展开讨论，主要研究内容和创新点如下： (1)在查阅大量国内外相关文献的基础上，对卫星重力探测技术的发展历程、国内外研究现状以及研究意义进行了阐述。讨论了目前常用的几种地球重力场模型反演算法，并分析了它们的优缺点。 (2)详细研究了卫星跟踪卫星重力测量所涉及到的时间系统、坐标系统以及空间运行卫星所受到的各类摄动力模型，给出了卫星重力测量的基本原理、常用的算法以及各类偏导
10、数的计算。定性分析了卫星初始位置误差对应用卫星精密轨道数据反演地球重力场模型的影响，若不顾及其它测量误差(如加速度计测量误差)的影响，以 2 小时反演弧长为宜。 (3)详细讨论了 CHAMP 卫星和 GRACE 卫星相关观测数据的预处理方法，分析了姿态数据间断的三种情况，给出了利用二次插值法对间断的姿态数据进行填补，算例表明：在姿态数据间断 40 历元的情况下，二次插值法完全能满足精度要求。 (4)利用动力学法对 CHAMP 卫星的加速度计观测数据进行了重新标定，结果表明：加速度计的标定参数并不是一成不变的，一天标定结果的变化范围比 1 小时的要小，更接近 GFZ 给出的结果。加速度计重新标定
11、后的重力场模型的精度要高于未标定的结果，因此对加速度计观测数据进行重新标定是十分必要的。 (5)详细介绍了笔者独立自主开发的卫星重力数据处理软件-GRASTAR，该软件是利用 Visual C+语言平台开发完成的，可对 CHAMP 卫星和 GRACE 卫星观测数据进行处理，具有很好的可移植性和可拓展性。 (6)基于应用低轨卫星精密轨道数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始状态向量与地球重力场位系数一起求解的数学模型，利用 120 天 CHAMP 卫星精密轨道和加速度计观测数据解算出 50 阶次的地球重力场模型 TJCHAMP01S，检核结果表明：TJ
12、CHAMP01S 模型精度优于相同阶次的 EGM96 和 EIGEN-1S 模型，并且更接近 EIGEN-2 模型。在我国区域，TJCHAMP01S 和 EIGEN-2 模型与 EIGEN-CG01C 模型前 30 阶所表示的大地水准面差值的中误差分别为 4.4cm 和 3.3cm。 (7)基于动力学法，利用 26 天 GRACE 卫星精密轨道和加速度计观测数据解算出 36 阶次的地球重力场模型 TJGRACE01S-OR，检核结果表明：该模型的精度要优于相同阶次的 EGM96 和 EIGEN-2 模型，前 20 阶模型位系数与 EIGEN-CHAMP03S 精度相当。 (8)基于应用 GRA
13、CE 卫星星间距离变率数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始相对速度向量与地球重力场位系数一起求解的数学模型，利用 31 天 GRACE 卫星星间距离变率和加速度计观测数据解算出 60 阶次的地球重力场模型 TJGRACE01S，检核结果表明：TJGRACE01S 模型的精度要明显优于相同阶次的 EIGEN-CHAMP03S 模型，但要稍逊于 EIGEN-GRACE01S。利用 EGM96、EIGEN-CHAMP03S、TJGRACE01S 和EIGEN-GRACE01S 模型的前 60 阶分别计算全球 1#176;1#176;格网点的大地水准面高度
14、，并与 EIGEN-GL04C 模型的大地水准面高度做比较，相应差值的中误差分别为：0.446m、0.144m、0.027m 和 0.008m。 (9)利用 67 天GRACE 卫星星间距离变率和加速度计观测数据解算出 80 阶次的地球重力场模型TJGRACE02S，检核结果表明：TJGRACE02S 模型的精度要明显优于相同阶次的EIGEN-CHAMP03S 模型，但要稍逊于 GGM01S。利用 EGM96、EIGEN-CHAMP03S、GGM01S 和 TJGRACE02S 模型的前 72 阶分别计算全球2.5#176;2.5#176;格网点的大地水准面高度和重力异常，并与EIGEN-GL
15、04C 模型的计算结果进行比较，大地水准面高度差值的中误差分别为：0.524m、0.305m、0.021m 和 0.052m，而重力异常差值的中误差分别为：4.00mGal、3.07 mGal、0.21 mGal 和 0.50 mGal。 (10)对本文的工作进行总结，对未来的工作进行展望和规划。低轨卫星跟踪测量技术(CHAMP 卫星和 GRACE 卫星)可提供覆盖全球的、高精度的空间重力观测资料：精密轨道数据、GRACE 卫星的 K 波段观测数据等，这些观测数据包含了丰富的地球重力场信息，尤其是 GRACE 卫星的 K 波段观测数据。在当代低轨卫星的定轨精度达到比较高的程度(约 2cm)，K
16、 波段观测数据将是地球重力场信息提取的主要载体，如何充分利用此类信息进行反演地球重力场模型的研究是当代大地测量研究的前沿和热点之一。其研究成果将对现代大地测量、地球物理、海洋科学以及地球动力学等学科的发展起到巨大的推动作用。 本文围绕应用低轨卫星跟踪数据反演地球重力场模型这一主题展开讨论，主要研究内容和创新点如下： (1)在查阅大量国内外相关文献的基础上，对卫星重力探测技术的发展历程、国内外研究现状以及研究意义进行了阐述。讨论了目前常用的几种地球重力场模型反演算法，并分析了它们的优缺点。 (2)详细研究了卫星跟踪卫星重力测量所涉及到的时间系统、坐标系统以及空间运行卫星所受到的各类摄动力模型，给
17、出了卫星重力测量的基本原理、常用的算法以及各类偏导数的计算。定性分析了卫星初始位置误差对应用卫星精密轨道数据反演地球重力场模型的影响，若不顾及其它测量误差(如加速度计测量误差)的影响，以 2 小时反演弧长为宜。 (3)详细讨论了 CHAMP 卫星和 GRACE 卫星相关观测数据的预处理方法，分析了姿态数据间断的三种情况，给出了利用二次插值法对间断的姿态数据进行填补，算例表明：在姿态数据间断 40 历元的情况下，二次插值法完全能满足精度要求。 (4)利用动力学法对 CHAMP 卫星的加速度计观测数据进行了重新标定，结果表明：加速度计的标定参数并不是一成不变的，一天标定结果的变化范围比 1 小时的
18、要小，更接近 GFZ 给出的结果。加速度计重新标定后的重力场模型的精度要高于未标定的结果，因此对加速度计观测数据进行重新标定是十分必要的。 (5)详细介绍了笔者独立自主开发的卫星重力数据处理软件-GRASTAR，该软件是利用 Visual C+语言平台开发完成的，可对 CHAMP 卫星和 GRACE 卫星观测数据进行处理，具有很好的可移植性和可拓展性。 (6)基于应用低轨卫星精密轨道数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始状态向量与地球重力场位系数一起求解的数学模型，利用 120 天 CHAMP 卫星精密轨道和加速度计观测数据解算出 50 阶次的地球重
19、力场模型 TJCHAMP01S，检核结果表明：TJCHAMP01S 模型精度优于相同阶次的 EGM96 和 EIGEN-1S 模型，并且更接近EIGEN-2 模型。在我国区域，TJCHAMP01S 和 EIGEN-2 模型与 EIGEN-CG01C 模型前 30 阶所表示的大地水准面差值的中误差分别为 4.4cm 和 3.3cm。 (7)基于动力学法，利用 26 天 GRACE 卫星精密轨道和加速度计观测数据解算出 36 阶次的地球重力场模型 TJGRACE01S-OR，检核结果表明：该模型的精度要优于相同阶次的 EGM96 和 EIGEN-2 模型，前 20 阶模型位系数与 EIGEN-CH
20、AMP03S 精度相当。 (8)基于应用 GRACE 卫星星间距离变率数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始相对速度向量与地球重力场位系数一起求解的数学模型，利用 31 天 GRACE 卫星星间距离变率和加速度计观测数据解算出 60 阶次的地球重力场模型 TJGRACE01S，检核结果表明：TJGRACE01S 模型的精度要明显优于相同阶次的 EIGEN-CHAMP03S 模型，但要稍逊于 EIGEN-GRACE01S。利用 EGM96、EIGEN-CHAMP03S、TJGRACE01S 和EIGEN-GRACE01S 模型的前 60 阶分别计算全球
21、 1#176;1#176;格网点的大地水准面高度，并与 EIGEN-GL04C 模型的大地水准面高度做比较，相应差值的中误差分别为：0.446m、0.144m、0.027m 和 0.008m。 (9)利用 67 天GRACE 卫星星间距离变率和加速度计观测数据解算出 80 阶次的地球重力场模型TJGRACE02S，检核结果表明：TJGRACE02S 模型的精度要明显优于相同阶次的EIGEN-CHAMP03S 模型，但要稍逊于 GGM01S。利用 EGM96、EIGEN-CHAMP03S、GGM01S 和 TJGRACE02S 模型的前 72 阶分别计算全球2.5#176;2.5#176;格网点
22、的大地水准面高度和重力异常，并与EIGEN-GL04C 模型的计算结果进行比较，大地水准面高度差值的中误差分别为：0.524m、0.305m、0.021m 和 0.052m，而重力异常差值的中误差分别为：4.00mGal、3.07 mGal、0.21 mGal 和 0.50 mGal。 (10)对本文的工作进行总结，对未来的工作进行展望和规划。低轨卫星跟踪测量技术(CHAMP 卫星和 GRACE 卫星)可提供覆盖全球的、高精度的空间重力观测资料：精密轨道数据、GRACE 卫星的 K 波段观测数据等，这些观测数据包含了丰富的地球重力场信息，尤其是 GRACE 卫星的 K 波段观测数据。在当代低轨
23、卫星的定轨精度达到比较高的程度(约 2cm)，K 波段观测数据将是地球重力场信息提取的主要载体，如何充分利用此类信息进行反演地球重力场模型的研究是当代大地测量研究的前沿和热点之一。其研究成果将对现代大地测量、地球物理、海洋科学以及地球动力学等学科的发展起到巨大的推动作用。 本文围绕应用低轨卫星跟踪数据反演地球重力场模型这一主题展开讨论，主要研究内容和创新点如下： (1)在查阅大量国内外相关文献的基础上，对卫星重力探测技术的发展历程、国内外研究现状以及研究意义进行了阐述。讨论了目前常用的几种地球重力场模型反演算法，并分析了它们的优缺点。 (2)详细研究了卫星跟踪卫星重力测量所涉及到的时间系统、坐
24、标系统以及空间运行卫星所受到的各类摄动力模型，给出了卫星重力测量的基本原理、常用的算法以及各类偏导数的计算。定性分析了卫星初始位置误差对应用卫星精密轨道数据反演地球重力场模型的影响，若不顾及其它测量误差(如加速度计测量误差)的影响，以 2 小时反演弧长为宜。 (3)详细讨论了 CHAMP 卫星和 GRACE 卫星相关观测数据的预处理方法，分析了姿态数据间断的三种情况，给出了利用二次插值法对间断的姿态数据进行填补，算例表明：在姿态数据间断 40 历元的情况下，二次插值法完全能满足精度要求。 (4)利用动力学法对 CHAMP 卫星的加速度计观测数据进行了重新标定，结果表明：加速度计的标定参数并不是
25、一成不变的，一天标定结果的变化范围比 1 小时的要小，更接近 GFZ 给出的结果。加速度计重新标定后的重力场模型的精度要高于未标定的结果，因此对加速度计观测数据进行重新标定是十分必要的。 (5)详细介绍了笔者独立自主开发的卫星重力数据处理软件-GRASTAR，该软件是利用 Visual C+语言平台开发完成的，可对 CHAMP 卫星和 GRACE 卫星观测数据进行处理，具有很好的可移植性和可拓展性。 (6)基于应用低轨卫星精密轨道数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始状态向量与地球重力场位系数一起求解的数学模型，利用 120 天 CHAMP 卫星精密
26、轨道和加速度计观测数据解算出 50 阶次的地球重力场模型 TJCHAMP01S，检核结果表明：TJCHAMP01S 模型精度优于相同阶次的 EGM96 和 EIGEN-1S 模型，并且更接近EIGEN-2 模型。在我国区域，TJCHAMP01S 和 EIGEN-2 模型与 EIGEN-CG01C 模型前 30 阶所表示的大地水准面差值的中误差分别为 4.4cm 和 3.3cm。 (7)基于动力学法，利用 26 天 GRACE 卫星精密轨道和加速度计观测数据解算出 36 阶次的地球重力场模型 TJGRACE01S-OR，检核结果表明：该模型的精度要优于相同阶次的 EGM96 和 EIGEN-2
27、模型，前 20 阶模型位系数与 EIGEN-CHAMP03S 精度相当。 (8)基于应用 GRACE 卫星星间距离变率数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始相对速度向量与地球重力场位系数一起求解的数学模型，利用 31 天 GRACE 卫星星间距离变率和加速度计观测数据解算出 60 阶次的地球重力场模型 TJGRACE01S，检核结果表明：TJGRACE01S 模型的精度要明显优于相同阶次的 EIGEN-CHAMP03S 模型，但要稍逊于 EIGEN-GRACE01S。利用 EGM96、EIGEN-CHAMP03S、TJGRACE01S 和EIGEN-
28、GRACE01S 模型的前 60 阶分别计算全球 1#176;1#176;格网点的大地水准面高度，并与 EIGEN-GL04C 模型的大地水准面高度做比较，相应差值的中误差分别为：0.446m、0.144m、0.027m 和 0.008m。 (9)利用 67 天GRACE 卫星星间距离变率和加速度计观测数据解算出 80 阶次的地球重力场模型TJGRACE02S，检核结果表明：TJGRACE02S 模型的精度要明显优于相同阶次的EIGEN-CHAMP03S 模型，但要稍逊于 GGM01S。利用 EGM96、EIGEN-CHAMP03S、GGM01S 和 TJGRACE02S 模型的前 72 阶分
29、别计算全球2.5#176;2.5#176;格网点的大地水准面高度和重力异常，并与EIGEN-GL04C 模型的计算结果进行比较，大地水准面高度差值的中误差分别为：0.524m、0.305m、0.021m 和 0.052m，而重力异常差值的中误差分别为：4.00mGal、3.07 mGal、0.21 mGal 和 0.50 mGal。 (10)对本文的工作进行总结，对未来的工作进行展望和规划。低轨卫星跟踪测量技术(CHAMP 卫星和 GRACE 卫星)可提供覆盖全球的、高精度的空间重力观测资料：精密轨道数据、GRACE 卫星的 K 波段观测数据等，这些观测数据包含了丰富的地球重力场信息，尤其是
30、GRACE 卫星的 K 波段观测数据。在当代低轨卫星的定轨精度达到比较高的程度(约 2cm)，K 波段观测数据将是地球重力场信息提取的主要载体，如何充分利用此类信息进行反演地球重力场模型的研究是当代大地测量研究的前沿和热点之一。其研究成果将对现代大地测量、地球物理、海洋科学以及地球动力学等学科的发展起到巨大的推动作用。 本文围绕应用低轨卫星跟踪数据反演地球重力场模型这一主题展开讨论，主要研究内容和创新点如下： (1)在查阅大量国内外相关文献的基础上，对卫星重力探测技术的发展历程、国内外研究现状以及研究意义进行了阐述。讨论了目前常用的几种地球重力场模型反演算法，并分析了它们的优缺点。 (2)详细
31、研究了卫星跟踪卫星重力测量所涉及到的时间系统、坐标系统以及空间运行卫星所受到的各类摄动力模型，给出了卫星重力测量的基本原理、常用的算法以及各类偏导数的计算。定性分析了卫星初始位置误差对应用卫星精密轨道数据反演地球重力场模型的影响，若不顾及其它测量误差(如加速度计测量误差)的影响，以 2 小时反演弧长为宜。 (3)详细讨论了 CHAMP 卫星和 GRACE 卫星相关观测数据的预处理方法，分析了姿态数据间断的三种情况，给出了利用二次插值法对间断的姿态数据进行填补，算例表明：在姿态数据间断 40 历元的情况下，二次插值法完全能满足精度要求。 (4)利用动力学法对 CHAMP 卫星的加速度计观测数据进
32、行了重新标定，结果表明：加速度计的标定参数并不是一成不变的，一天标定结果的变化范围比 1 小时的要小，更接近 GFZ 给出的结果。加速度计重新标定后的重力场模型的精度要高于未标定的结果，因此对加速度计观测数据进行重新标定是十分必要的。 (5)详细介绍了笔者独立自主开发的卫星重力数据处理软件-GRASTAR，该软件是利用 Visual C+语言平台开发完成的，可对 CHAMP 卫星和 GRACE 卫星观测数据进行处理，具有很好的可移植性和可拓展性。 (6)基于应用低轨卫星精密轨道数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始状态向量与地球重力场位系数一起求解的
33、数学模型，利用 120 天 CHAMP 卫星精密轨道和加速度计观测数据解算出 50 阶次的地球重力场模型 TJCHAMP01S，检核结果表明：TJCHAMP01S 模型精度优于相同阶次的 EGM96 和 EIGEN-1S 模型，并且更接近EIGEN-2 模型。在我国区域，TJCHAMP01S 和 EIGEN-2 模型与 EIGEN-CG01C 模型前 30 阶所表示的大地水准面差值的中误差分别为 4.4cm 和 3.3cm。 (7)基于动力学法，利用 26 天 GRACE 卫星精密轨道和加速度计观测数据解算出 36 阶次的地球重力场模型 TJGRACE01S-OR，检核结果表明：该模型的精度要
34、优于相同阶次的 EGM96 和 EIGEN-2 模型，前 20 阶模型位系数与 EIGEN-CHAMP03S 精度相当。 (8)基于应用 GRACE 卫星星间距离变率数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始相对速度向量与地球重力场位系数一起求解的数学模型，利用 31 天 GRACE 卫星星间距离变率和加速度计观测数据解算出 60 阶次的地球重力场模型 TJGRACE01S，检核结果表明：TJGRACE01S 模型的精度要明显优于相同阶次的 EIGEN-CHAMP03S 模型，但要稍逊于 EIGEN-GRACE01S。利用 EGM96、EIGEN-CHA
35、MP03S、TJGRACE01S 和EIGEN-GRACE01S 模型的前 60 阶分别计算全球 1#176;1#176;格网点的大地水准面高度，并与 EIGEN-GL04C 模型的大地水准面高度做比较，相应差值的中误差分别为：0.446m、0.144m、0.027m 和 0.008m。 (9)利用 67 天GRACE 卫星星间距离变率和加速度计观测数据解算出 80 阶次的地球重力场模型TJGRACE02S，检核结果表明：TJGRACE02S 模型的精度要明显优于相同阶次的EIGEN-CHAMP03S 模型，但要稍逊于 GGM01S。利用 EGM96、EIGEN-CHAMP03S、GGM01S
36、 和 TJGRACE02S 模型的前 72 阶分别计算全球2.5#176;2.5#176;格网点的大地水准面高度和重力异常，并与EIGEN-GL04C 模型的计算结果进行比较，大地水准面高度差值的中误差分别为：0.524m、0.305m、0.021m 和 0.052m，而重力异常差值的中误差分别为：4.00mGal、3.07 mGal、0.21 mGal 和 0.50 mGal。 (10)对本文的工作进行总结，对未来的工作进行展望和规划。低轨卫星跟踪测量技术(CHAMP 卫星和 GRACE 卫星)可提供覆盖全球的、高精度的空间重力观测资料：精密轨道数据、GRACE 卫星的 K 波段观测数据等，
37、这些观测数据包含了丰富的地球重力场信息，尤其是 GRACE 卫星的 K 波段观测数据。在当代低轨卫星的定轨精度达到比较高的程度(约 2cm)，K 波段观测数据将是地球重力场信息提取的主要载体，如何充分利用此类信息进行反演地球重力场模型的研究是当代大地测量研究的前沿和热点之一。其研究成果将对现代大地测量、地球物理、海洋科学以及地球动力学等学科的发展起到巨大的推动作用。 本文围绕应用低轨卫星跟踪数据反演地球重力场模型这一主题展开讨论，主要研究内容和创新点如下： (1)在查阅大量国内外相关文献的基础上，对卫星重力探测技术的发展历程、国内外研究现状以及研究意义进行了阐述。讨论了目前常用的几种地球重力场
38、模型反演算法，并分析了它们的优缺点。 (2)详细研究了卫星跟踪卫星重力测量所涉及到的时间系统、坐标系统以及空间运行卫星所受到的各类摄动力模型，给出了卫星重力测量的基本原理、常用的算法以及各类偏导数的计算。定性分析了卫星初始位置误差对应用卫星精密轨道数据反演地球重力场模型的影响，若不顾及其它测量误差(如加速度计测量误差)的影响，以 2 小时反演弧长为宜。 (3)详细讨论了 CHAMP 卫星和 GRACE 卫星相关观测数据的预处理方法，分析了姿态数据间断的三种情况，给出了利用二次插值法对间断的姿态数据进行填补，算例表明：在姿态数据间断 40 历元的情况下，二次插值法完全能满足精度要求。 (4)利用
39、动力学法对 CHAMP 卫星的加速度计观测数据进行了重新标定，结果表明：加速度计的标定参数并不是一成不变的，一天标定结果的变化范围比 1 小时的要小，更接近 GFZ 给出的结果。加速度计重新标定后的重力场模型的精度要高于未标定的结果，因此对加速度计观测数据进行重新标定是十分必要的。 (5)详细介绍了笔者独立自主开发的卫星重力数据处理软件-GRASTAR，该软件是利用 Visual C+语言平台开发完成的，可对 CHAMP 卫星和 GRACE 卫星观测数据进行处理，具有很好的可移植性和可拓展性。 (6)基于应用低轨卫星精密轨道数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差
40、以及卫星初始状态向量与地球重力场位系数一起求解的数学模型，利用 120 天 CHAMP 卫星精密轨道和加速度计观测数据解算出 50 阶次的地球重力场模型 TJCHAMP01S，检核结果表明：TJCHAMP01S 模型精度优于相同阶次的 EGM96 和 EIGEN-1S 模型，并且更接近EIGEN-2 模型。在我国区域，TJCHAMP01S 和 EIGEN-2 模型与 EIGEN-CG01C 模型前 30 阶所表示的大地水准面差值的中误差分别为 4.4cm 和 3.3cm。 (7)基于动力学法，利用 26 天 GRACE 卫星精密轨道和加速度计观测数据解算出 36 阶次的地球重力场模型 TJGR
41、ACE01S-OR，检核结果表明：该模型的精度要优于相同阶次的 EGM96 和 EIGEN-2 模型，前 20 阶模型位系数与 EIGEN-CHAMP03S 精度相当。 (8)基于应用 GRACE 卫星星间距离变率数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始相对速度向量与地球重力场位系数一起求解的数学模型，利用 31 天 GRACE 卫星星间距离变率和加速度计观测数据解算出 60 阶次的地球重力场模型 TJGRACE01S，检核结果表明：TJGRACE01S 模型的精度要明显优于相同阶次的 EIGEN-CHAMP03S 模型，但要稍逊于 EIGEN-GRA
42、CE01S。利用 EGM96、EIGEN-CHAMP03S、TJGRACE01S 和EIGEN-GRACE01S 模型的前 60 阶分别计算全球 1#176;1#176;格网点的大地水准面高度，并与 EIGEN-GL04C 模型的大地水准面高度做比较，相应差值的中误差分别为：0.446m、0.144m、0.027m 和 0.008m。 (9)利用 67 天GRACE 卫星星间距离变率和加速度计观测数据解算出 80 阶次的地球重力场模型TJGRACE02S，检核结果表明：TJGRACE02S 模型的精度要明显优于相同阶次的EIGEN-CHAMP03S 模型，但要稍逊于 GGM01S。利用 EGM
43、96、EIGEN-CHAMP03S、GGM01S 和 TJGRACE02S 模型的前 72 阶分别计算全球2.5#176;2.5#176;格网点的大地水准面高度和重力异常，并与EIGEN-GL04C 模型的计算结果进行比较，大地水准面高度差值的中误差分别为：0.524m、0.305m、0.021m 和 0.052m，而重力异常差值的中误差分别为：4.00mGal、3.07 mGal、0.21 mGal 和 0.50 mGal。 (10)对本文的工作进行总结，对未来的工作进行展望和规划。低轨卫星跟踪测量技术(CHAMP 卫星和 GRACE 卫星)可提供覆盖全球的、高精度的空间重力观测资料：精密轨
44、道数据、GRACE 卫星的 K 波段观测数据等，这些观测数据包含了丰富的地球重力场信息，尤其是 GRACE 卫星的 K 波段观测数据。在当代低轨卫星的定轨精度达到比较高的程度(约 2cm)，K 波段观测数据将是地球重力场信息提取的主要载体，如何充分利用此类信息进行反演地球重力场模型的研究是当代大地测量研究的前沿和热点之一。其研究成果将对现代大地测量、地球物理、海洋科学以及地球动力学等学科的发展起到巨大的推动作用。 本文围绕应用低轨卫星跟踪数据反演地球重力场模型这一主题展开讨论，主要研究内容和创新点如下： (1)在查阅大量国内外相关文献的基础上，对卫星重力探测技术的发展历程、国内外研究现状以及研
45、究意义进行了阐述。讨论了目前常用的几种地球重力场模型反演算法，并分析了它们的优缺点。 (2)详细研究了卫星跟踪卫星重力测量所涉及到的时间系统、坐标系统以及空间运行卫星所受到的各类摄动力模型，给出了卫星重力测量的基本原理、常用的算法以及各类偏导数的计算。定性分析了卫星初始位置误差对应用卫星精密轨道数据反演地球重力场模型的影响，若不顾及其它测量误差(如加速度计测量误差)的影响，以 2 小时反演弧长为宜。 (3)详细讨论了 CHAMP 卫星和 GRACE 卫星相关观测数据的预处理方法，分析了姿态数据间断的三种情况，给出了利用二次插值法对间断的姿态数据进行填补，算例表明：在姿态数据间断 40 历元的情
46、况下，二次插值法完全能满足精度要求。 (4)利用动力学法对 CHAMP 卫星的加速度计观测数据进行了重新标定，结果表明：加速度计的标定参数并不是一成不变的，一天标定结果的变化范围比 1 小时的要小，更接近 GFZ 给出的结果。加速度计重新标定后的重力场模型的精度要高于未标定的结果，因此对加速度计观测数据进行重新标定是十分必要的。 (5)详细介绍了笔者独立自主开发的卫星重力数据处理软件-GRASTAR，该软件是利用 Visual C+语言平台开发完成的，可对 CHAMP 卫星和 GRACE 卫星观测数据进行处理，具有很好的可移植性和可拓展性。 (6)基于应用低轨卫星精密轨道数据反演地球重力场模型
47、的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始状态向量与地球重力场位系数一起求解的数学模型，利用 120 天 CHAMP 卫星精密轨道和加速度计观测数据解算出 50 阶次的地球重力场模型 TJCHAMP01S，检核结果表明：TJCHAMP01S 模型精度优于相同阶次的 EGM96 和 EIGEN-1S 模型，并且更接近EIGEN-2 模型。在我国区域，TJCHAMP01S 和 EIGEN-2 模型与 EIGEN-CG01C 模型前 30 阶所表示的大地水准面差值的中误差分别为 4.4cm 和 3.3cm。 (7)基于动力学法，利用 26 天 GRACE 卫星精密轨道和加速度计观测
48、数据解算出 36 阶次的地球重力场模型 TJGRACE01S-OR，检核结果表明：该模型的精度要优于相同阶次的 EGM96 和 EIGEN-2 模型，前 20 阶模型位系数与 EIGEN-CHAMP03S 精度相当。 (8)基于应用 GRACE 卫星星间距离变率数据反演地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始相对速度向量与地球重力场位系数一起求解的数学模型，利用 31 天 GRACE 卫星星间距离变率和加速度计观测数据解算出 60 阶次的地球重力场模型 TJGRACE01S，检核结果表明：TJGRACE01S 模型的精度要明显优于相同阶次的 EIGEN-CHA
49、MP03S 模型，但要稍逊于 EIGEN-GRACE01S。利用 EGM96、EIGEN-CHAMP03S、TJGRACE01S 和EIGEN-GRACE01S 模型的前 60 阶分别计算全球 1#176;1#176;格网点的大地水准面高度，并与 EIGEN-GL04C 模型的大地水准面高度做比较，相应差值的中误差分别为：0.446m、0.144m、0.027m 和 0.008m。 (9)利用 67 天GRACE 卫星星间距离变率和加速度计观测数据解算出 80 阶次的地球重力场模型TJGRACE02S，检核结果表明：TJGRACE02S 模型的精度要明显优于相同阶次的EIGEN-CHAMP03S 模型，但要稍逊于 GGM01S。利用 EGM96、EIGEN-CHAMP03S、GGM01S 和 TJGRACE02S 模型的前 72 阶分别计算全球2.5#176;2.5#176;格网点的大地水准面高度和重力异常，并与EIGEN-GL04C 模型的计算结果进行比较，大地水准面高度差值的中误差分别为：0.524m、0.305m、0.021m 和 0.052m，而重力异常差值的中误差分别为：4.00m
