最小二乘法求解三类卫星重力梯度边值问题的研究

研究了最小二乘法求解3类卫星重力梯度边值问题的理论和方法，给出了3类梯度观测值{Гzz}、{Гxz、Гyz}和{Гxx-Гyy,2Гxy}对应边值问题解的核函数严密表达式。模拟试算结果表明，最小二乘法求解的卫星重力梯度积分公式用于恢复地球重力场是有效而严密的。     

Tikhonov正则化方法在GOCE重力场求解中的模拟研究

本文在阐述Tikhonov正则化方法基本原理的基础上,给出了四类可用于重力场解算的正则化矩阵(零次、一次、二次和Kaula),以及用于确定正则化参数的L曲线法和GCV方法的数学模型。基于SA方法利用模拟数据分析讨论了零次、一次以及Kaula正则化矩阵应用于GOCE全球重力场模型确定的有效性,并由Kaula正则化矩阵分析了L曲线法和GCV方法确定正则化参数的可行性。数值结果表明三类正则化矩阵获得的最优解(以大地水准面MSE最小为准则确定)的精度水平相近,关键在于相应正则化参数的确定,数值结果同时说明了GCV方法和L曲线法可用于确定正则化参数,且前者较后者具有更好的稳定性。     

卫星重力资料分析的同解法研究及其仿真

对同解法的基本思想和流程进行归纳,研制相应的软件系统.利用模拟的GPS和星间测距数据进行同时确定重力场模型和卫星精密轨道的仿真模拟,验证本文设计思路的可行性和软件系统的可靠性.     

一种新型重力测量卫星系统确定全球重力场的性能分析

本文设计了一种高-低卫星跟踪卫星、低-低卫星跟踪卫星和卫星重力梯度测量相结合的新型重力测量卫星系统,其可在一定程度上发挥卫星重力梯度和低低卫星跟踪卫星两种测量模式各自的优势.基于重力卫星系统指标设计的半解析法,深入分析了不同重力测量卫星系统配置和不同观测量及其不同白噪声水平情况下,新型卫星重力测量模式反演重力场模型的能力.数值模拟分析结果表明：在观测值精度和星间距离相同的条件下,轨道高度是影响重力场反演精度的关键因素;随着星间距离的增大,高频重力场信号反演精度会先提高后降低,轨道高度在200~350 km之间时,星间距离在150~180 km之间时反演精度最优;星间距离变率和卫星重力梯度两类观测值仅在某些精度配置时可达到优势互补,如果某一类观测值精度很高,则另一类观测值在联合解算时贡献非常小或者没有贡献.在300 km轨道高度,若以GRACE和GOCE任务的设计指标1 μm·s^-1/√Hz和5 mE/√Hz来配置新型重力测量卫星系统中星间距离变率和引力梯度观测值的精度,联合两类观测值解算200阶次模型大地水准面的精度比独立解算分别提高1.2倍和2.8倍.如果以实现100 km空间分辨率1~2 cm精度大地水准面为科学目标,考虑卫星在轨寿命,建议轨道高度选择300 km,星间距离变率和卫星重力梯度的精度分别为0.1 μm·s^-1/√Hz和1 mE/√Hz.本文的研究成果可为中国研制自主的重力测量卫星系统提供参考依据.     

GOCE卫星引力梯度数据滤波方法研究

针对GOCE卫星引力梯度观测值中低精度分量和低频有色噪声的处理策略问题,该文采用模型模拟值代替低精度分量Vxy和Vyz,以减弱低精度分量在坐标系转换中对高精度分量的影响。深入分析比较了多种滤波方法处理GOCE卫星引力梯度观测值中有色噪声的效果,提出采用Butterworth零相移滤波方法加移去-恢复技术的思路,实测数据的处理效果验证了该方法的有效性。     

极空白对三类重力梯度边值问题球谐分析解的影响分析

就极空白对三类重力梯度边值问题球谐分析解的影响进行了定性探讨和定量分析,结果表明,基于球面边界的垂直-垂直重力梯度边值问题球谐分析解的零次项、垂直-水平重力梯度边值问题球谐分析解的一次项、水平-水平重力梯度边值问题球谐分析解的二次项受极空白问题的影响最为显著。     

区域与全球高程基准差异的确定

全球高程基准统一是继全球大地测量坐标系及其参考基准统一之后,大地测量学科面临和亟待解决的一个重要问题,也是全球空间信息共享与交换的基础。本文针对区域高程基准与全球高程基准间基准差异确定的理论、方法及实际问题开展研究。利用物理大地测量高程系统的经典理论方法,给出了高程基准差异的定义,并推导了计算基准差异的严密公式,该公式可将高程基准差异确定的现有3种方法统一起来。在此基础上,分析顾及了不同椭球参数对于计算基准差异的影响及量级,同时,高程异常差法还需考虑全球高程基准重力位与模型计算大地水准面位值不一致引起的零阶项改正。利用青岛原点附近152个GPS水准点数据,分别选择GRS80、WGS-84、CGCS2000参考椭球以及EGM2008、EIGEN-6C4、SGG-UGM-1模型,采用位差法和高程异常差法,确定了我国1985高程基准与全球高程基准的差异。其中,EIGEN-6C4模型计算的我国高程基准与WGS-84参考椭球正常重力位U0定义的全球高程基准之间的差异约为-23.1cm。也就是说,我国高程基准低于采用WGS-84参考椭球正常重力位U0定义的全球高程基准,当选取基于平均海面确定的Gauss-Listing大地水准面作为全球高程基准时,我国1985高程基准高于全球基准约21.0cm。从计算结果还可看出,当前重力场模型在青岛周边不同GPS/水准点的精度差别依然较大,这会导致选择不同数据对确定我国85国家高程基准与全球基准之间的差异影响较大,因此,若要实现厘米级精度区域高程基准与全球高程基准的统一,全球重力场模型的精度和可靠性还需要进一步提高。     

OpenMP并行计算在卫星重力数据处理中的应用

本文对新一代卫星重力数据处理中的计算密集型任务进行了分析,总结出需要采用并行计算技术提高效率的几个关键任务。对不同的并行手段进行了比较,采用了OpenMP并行方法,并通过算例验证了并行设计方法的有效性。结果表明,并行计算能明显提高GRACE/GOCE任务的数据处理效率。