基于重力卫星几何轨道线性化的地球重力场反演方法

传统动力学法的观测方程以6个初始轨道参数和先验力模型为初值进行线性化,其线性化误差随积分弧长拉长而增大.本文直接以重力卫星的几何观测轨道为初值进行线性化,其线性化误差与轨道弧长无关,且不需要初始重力场模型和初始轨道参数.导出了基于卫星轨道观测值反演重力场模型的相关公式,利用JPL公布的RL02版本2008年全年的GRACE双星轨道数据和加速度计数据解算了90阶次的地球重力场模型TJGRACE01S,并以EGM2008模型为基准与其他模型进行了比较分析,结果表明:TJGRACE01S模型直到90阶次的大地水准面累积误差为17.6 cm,优于同阶次的EIGEN-CHAMP03S和EIGEN-CHAMP05S模型,前27阶位系数整体精度优于EIGEN-GRACE01S,前15阶位系数整体精度与EIGEN-GRACE02S模型精度大致相当.利用美国8221个GPS水准点数据的分析结果也表明,本文模型也优于同阶次的EIGEN-CHAMP03S和EIGEN-CHAMP05S模型.     

基于GRACE卫星重力数据确定地球重力场模型WHU-GM-05

基于卫星轨道运动的能量积分方程,可导出利用卫星跟踪卫星数据求解地球重力场的实用公式.本文在Jekeli给出的公式基础上导出了基于能量守恒方程利用两颗低-低卫星跟踪的扰动位差求解重力位系数的严密关系式.基于两颗GRACE卫星的观测数据,采用本文导出的严密能量积分方法求解得到120阶的GRACE地球重力场模型,命名为WHU-GM-05;将WHU-GM-05模型与国际上同类重力场模型EIGEN-GRACE系列和GGM02S分别在阶方差和大地水准面高等方面作了比较,并与美国和中国的部分地区GPS水准观测值进行了精度分析.结果表明基于本文推导的严密双星能量守恒方程得到的WHU-GM-05重力场模型精度与国际上同类重力场模型的精度相当.     

基于CHAMP短弧长动力学轨道的地球重力场模型

讨论了基于CHAMP卫星动力学轨道数据以及加速度计数据推求地球重力场模型的动力学法，推导了将加速度计观测数据的尺度和偏差以及卫星初始状态向量与地球重力场位系数一起求解的数学模型. 采用CHAMP卫星120天的动力学轨道数据和加速度数据解算出50阶次的地球重力场模型TJCHAMP01S,并利用各种方法对该模型进行了检核，结果表明：TJCHAMP01S模型精度优于相同阶次的EGM96和EIGEN_1S模型.     

基于新型能量插值法精确建立GRACE-only地球重力场模型

本文基于新型能量插值法,利用美国喷气推进实验室(JPL)公布的2008年的GRACE-Level-1B实测数据,反演了120阶GRACE地球重力场.首先,由于GPS轨道测量精度相对较低,通过将K波段测距仪高精度的星间距离观测量插值引入双星动能差中,进而建立了新型能量插值卫星观测方程.其次,详细对比分析了2点、4点、6点和8点能量插值观测方程对地球重力场反演精度的影响.研究结果表明:基于最优的信噪比,6点能量插值公式有利于提高120阶GRACE地球重力场的反演精度.最后,基于美国、欧洲和澳大利亚的GPS/水准观测数据检验了本文新建立的WHIGG-GEGM03S地球重力场模型的正确性和有效性.     

基于下一代三向车轮双星编队改善地球重力场空间分辨率研究

由于当前GRACE(Gravity Recovery and Climate Experiment)串行式编队存在"南北向条带误差"等缺陷,因此本文基于星间速度插值法开展了利用下一代三向车轮双星编队ACR(Along-Cross-Radial)-Cartwheel提高地球重力场空间分辨率的可行性研究论证.第一,采用GRACE卫星轨道参数和关键载荷精度,利用三向车轮双星编队ACR-Cartwheel-A/B反演了120阶地球重力场.结果表明:基于ACR-Cartwheel-A/B双星编队反演地球重力场的模拟精度较德国波茨坦地学研究中心(GFZ)公布的EIGEN-GRACE02S地球重力场模型的实测精度平均提高2.6倍,从而检验了基于下一代三向车轮双星编队ACR-Cartwheel-A/B反演地球重力场精度优于当前GRACE串行式双星编队的可行性.第二,通过星间速度插值法,采用卫星轨道参数(初始轨道高度350km、平均星间距离100km、初始轨道倾角89°、初始轨道离心率0.0046)、卫星关键载荷精度指标(星间速度10-7 m·s-1、轨道位置10-3 m、轨道速度10-6 m·s-1、非保守力10-11 m·s-2)、观测时间30天和采样间隔10s,基于经向车轮双星编队Lo-AR(Longitudinal-Along-Radial)-Cartwheel-A/B、纬向车轮双星编队La-AR(Latitudinal-Along-Radial)-Cartwheel-A/B和三向车轮双星编队ACR-Cartwheel-A/B,分别反演了120阶地球重力场;在120阶处,累计大地水准面精度分别为5.115×10-4 m、4.923×10-4 m和3.488×10-4 m.结果表明:(1)由于La-AR-Cartwheel-A/B编队的轨道稳定性优于Lo-AR-Cartwheel-A/B编队,因此基于La-AR-Cartwheel-A/B编队反演重力场精度高于Lo-AR-CartwheelA/B编队;(2)由于ACR-Cartwheel-A/B编队可以同时获得轨向、垂向和径向的重力场信息,卫星观测数据具有各向同性优点,因此ACR-Cartwheel-A/B编队是建立下一代高精度和高空间分辨地球重力场模型的优化选择.     

GRACE卫星关键载荷实测数据的有效处理和地球重力场的精确解算

首先对美国喷气推进实验室（JPL）公布的2007-06-01～2007-12-31时间段内的GRACE Level 1B卫星GPS轨道位置和速度、K波段系统星间速度、加速度计非保守力以及恒星敏感器姿态实测数据对应进行了轨道拼接、粗差探测、线性内插、重新标定、坐标转换、误差分析等有效处理;其次,基于改进的能量守恒法恢复了120阶GRACE地球重力场,在120阶处累计大地水准面精度为25.313 cm;最后,检验了本文地球重力场模型IGG-GRACE的可靠性,同时分析了IGG-GRACE解算精度在低频部分略优于德国波兹坦地学研究中心（GFZ）公布的EIGEN-GRACE02S地球重力场模型,而在中高频部分略低的原因.     

基于星间加速度法精确和快速确定GRACE地球重力场

本文基于星间加速度法(IRAM),利用预处理共轭梯度迭代法(PCCG)开展了地球重力场恢复的模拟研究.第一,通过9点Newton-Gregory插值法得到星问加速度,并联合星间距离、星间速度和星间加速度建立了卫星观测方程.第二,在120阶处,恢复GRACE-IRAM累计大地水准面的精度为7.215 cm.第三,分析了在...     

基于半解析法有效和快速估计GRACE全球重力场的精度

首先基于半解析法建立了新的GRACE卫星K波段测量系统星间测速、GPS接收机轨道位置和加速度计非保守力误差联合影响累计大地水准面的误差模型;其次,基于各关键载荷精度指标的匹配关系,论证了误差模型的可靠性;最后,基于美国喷气动力实验室（JPL）公布的2006年的GRACE Level 1B实测误差数据,有效和快速地估计了120阶全球重力场的精度,在120阶处累计大地水准面的精度为18.368 cm,其结果和德国地学研究中心（GFZ）公布的EIGEN-GRACE02S全球重力场模型符合较好. 本文的研究为将来国际卫星重力测量计划（如GRACE Follow-On, 360阶）中高阶全球重力场模型精度的有效和快速估计提供了理论基础和计算保证.     

联合串行式和钟摆式卫星编队精确建立下一代HIP-3S地球重力场模型

第一,基于扰动星间距离观测量对地球重力场反演精度的敏感性优于星间距离观测值的特性,本文构建了新型扰动星间距离法(DIRM).第二,有效检验了下一代HIP-3S编队的轨道稳定性,结果表明:HIP-3S编队较稳定,有利于提高地球重力场反演精度.第三,基于扰动星间距离法,分别利用当前GRACE-2S串行式双星编队和下一代HIP-3S复合式三星编队精确反演了120阶地球重力场,在120阶处累计大地水准面精度为2.271×10~(-1)m和1.923×10~(-3)m,结果表明:HIP-3S复合式三星编队有利于建立下一代高精度和高空间分辨率的地球重力场模型.     

基于GRACE Follow-On卫星重力梯度法精确反演地球重力场

由于GRACE Follow-On双星系统等效于基线长为星间距离的一维水平重力梯度仪,因此本文基于GRACE Follow-On卫星重力梯度法开展了精确和快速反演下一代地球重力场的可行性论证研究. 研究结果表明：第一,基于GRACE Follow-On卫星重力梯度法（GFO-SGGM）,利用卫星轨道参数（轨道高度250 km、星间距离50 km、轨道倾角89°、轨道离心率0.001）、关键载荷测量精度（星间距离10^-6 m、星间速度10^-7 m·s^-1、星间加速度10^-10 m·s^-2、轨道位置10^-3 m、轨道速度10^-6 m·s^-1、非保守力10^-11 m·s^-2）、观测时间30天和采样间隔10 s反演了120阶地球重力场,在120阶处累计大地水准面精度为9.331×10^-4 m. 第二,在120阶内,利用将来GRACE Follow-On双星反演地球重力场精度较现有GRACE双星平均提高61倍,因此GRACE Follow-On卫星重力梯度法是进一步提高地球重力场反演精度的优选方法. 第三,下一代GRACE Follow-On计划较当前GRACE计划的优点如下：轨道高度更低（200~300 km）、载荷精度更高（10^-7 ~10^-9 m·s^-1）和星间距离更短（50~100 km）.     

卫星重力反演的短弧长积分法研究

给出了统一求解球谐位系数、弧段边界轨道改正向量、有偏距离改正及加速度计偏差的短弧长积分法,通过对力模型梯度改正减弱了轨道误差对反演地球重力场的影响.采用GRACE卫星1个月的实测轨道及星间距离数据计算表明,短弧长积分法加了梯度改正的精度比不加梯度改正整体提高了近一倍,且该方法在高阶次位系数的精度优于动力学法.基于GRA...     

基于卫-卫跟踪观测技术利用能量守恒法恢复地球重力场的数值模拟研究

根据卫－卫跟踪观测技术的测量原理,基于能量守恒法建立了一种新的双星相互跟踪和三星相互跟踪的卫星观测方程. 通过数值模拟,采用预处理共轭梯度法恢复120阶地球重力场. 模拟结果表明：第一,双星相互跟踪恢复地球重力场的精度和美国喷气动力实验室公布的EIGEN GRACE02S的结果相符合;第二,三星相互跟踪恢复地球重力场的精度较双星提高约2倍.     

基于时空域混合法利用Kaula正则化精确和快速解算GOCE地球重力场

为了研究卫星重力梯度技术对中高频地球重力场反演精度的影响,本文基于时空域混合法,利用Kaula正则化反演了250阶GOCE地球重力场.模拟结果表明:第一,时空域混合法是精确和快速求解高阶地球重力场的有效方法;第二,Kaula正则化是降低正规阵病态性的重要方法;第三,基于改进的预处理共轭梯度迭代法可快速求解大型线性方程组...     

基于残余星间速度法精确和快速反演下一代GRACE Follow-On地球重力场

基于新型残余星间速度法（RIRM）反演了120阶GRACE Follow-On地球重力场. 第一,由于GPS定轨精度相对较低,通过将激光干涉测距仪的高精度残余星间速度（测量精度10^-7 m·s^-1）引入残余轨道速度差分矢量的视线分量构建了新型RIRM观测方程. 第二,基于2点、4点、6点和8点RIRM公式对比论证了最优的插值点数. 如果相关系数和采样间隔一定,随着插值点数的增加,卫星观测值的信号量被有效加强,而卫星观测值的误差量也同时增加. 因此,6点RIRM公式是提高下一代地球重力场精度的较优选择. 第三,相关系数对地球重力场精度的影响在不同频段表现为不同特性. 随着相关系数的逐渐增大,地球长波重力场精度逐渐降低,而地球中长波重力场精度逐渐升高. 第四,基于6点RIRM公式,通过30天观测数据和采样间隔5 s,分别利用星间速度和残余星间速度观测值,在120阶次处反演下一代GRACE Follow-On累计大地水准面精度为1.638×10^-3 m和1.396×10^-3 m. 研究结果表明：（1）残余星间速度观测量较星间速度对地球重力场反演精度更敏感;（2）GRACE Follow-On地球重力场精度较GRACE至少高10倍.