联合高低卫-卫跟踪和卫星重力梯度数据恢复地球重力场的谱组合法

GOCE采用的高低卫-卫跟踪和卫星重力梯度测量技术在恢复重力场方面各有所长并互为补充,如何有效利用这两类观测数据最优确定地球重力场是GOCE重力场反演的关键问题。本文研究了联合高低卫-卫跟踪和卫星重力梯度数据恢复地球重力场的最小二乘谱组合法,基于球谐分析方法推导并建立了卫星轨道面扰动位T和径向重力梯度Tzz、以及扰动位T和重力梯度分量组合{Tzz-Txx-Tyy}的谱组合计算模型与误差估计公式。数值模拟结果表明,谱组合计算模型可以有效顾及各类数据的精度和频谱特性进行最优联合求解。采用61天GOCE实测数据反演的两个180阶次地球重力场模型WHU_GOCE_SC01S（扰动位和径向重力梯度数据求解）和WHU_GOCE_SC02S（扰动位和重力梯度分量组合数据求解）,结果显示后者精度优于前者,并且它们的整体精度优于GOCE时域解,而与GOCE空域解的精度接近,验证了谱组合法的可行性与有效性。     

利用GOCE卫星轨道反演地球重力场模型

根据积分方程法反演地球重力场的数学模型,利用GOCE卫星2009-11-02～2010-01-02共61d的精密轨道数据反演了几组地球重力场模型。结果表明,GOCE卫星轨道能有效提取地球重力场的长波信息,弥补了GOCE卫星重力梯度带宽的限制,在106阶次的大地水准面误差为±9.6cm,该阶次精度优于EIGEN-CHAMP03S及GRACE卫星两个月轨道反演地球重力场的精度,但由于两极空白,反演的带谐位系数精度偏低。联合GOCE及GRACE卫星轨道反演的模型在106阶次的大地水准面误差为±6.9cm,弥补了GOCE卫星轨道的缺陷。     

GOCE实测数据反演高阶重力场模型的Torus方法

不同于当前广泛使用的空域法、时域法、直接解法,本文尝试采用Torus方法处理GOCE实测数据,利用71 d的GOCE卫星引力梯度数据反演了200阶次GOCE地球重力场模型,实现了对参考模型的精化。首先,采用Butterworth零相移滤波方法加移去—恢复技术,处理引力梯度观测值中的有色噪声,并利用泰勒级数展开和Kriging方法对GOCE卫星引力梯度数据进行归算和格网化,计算得到了名义轨道上格网点处的引力梯度数据。然后,利用2D-FFT技术和块对角最小二乘方法处理名义轨道上数据,获得了200阶次的GOCE地球重力场模型GOCE_Torus。利用中国和美国的GPS/水准数据进行外部检核结果说明,GOCE_Torus与ESA发布的同期模型的精度相当;GOCE_Torus模型与200阶次的EGM2008模型相比,在美国区域精度相当,但在中国区域精度提高了4.6 cm,这充分体现了GOCE卫星观测数据对地面重力稀疏区的贡献。Torus方法拥有快速高精度反演卫星重力场模型的优势,可以在重力梯度卫星的设计、误差分析及在轨快速评估等方面得到充分应用。     

利用最小二乘直接法反演卫星重力场模型的MPI并行算法

针对海量卫星重力数据反演高阶次地球重力场模型的密集型计算任务与高内存耗用问题,基于MPI实现了最小二乘直接法恢复高阶次位系数的并行算法。引入并行读写、分块存储与分块计算等方式完成了设计矩阵的构建、法方程的形成与求解等密集型计算任务的并行算法,数值计算结果表明三者的并行相对效率峰值可分别达到95%、68%、63%。利用GOCE轨道跟踪和径向扰动重力梯度数据(共518 400个历元)分别反演了120、240阶次地球重力场模型,计算时间仅为40 min、7 h,内存耗用峰值仅为290 MB、1.57 GB;采用与GOCE同等噪声水平的观测数据恢复的重力场模型精度与GOCE已发布模型的解算精度相一致,联合GRACE和GOCE的解算模型能够实现二者独立信息的频谱互补,表明本文方法可高效稳定地恢复高阶次地球重力场模型。     

利用GOCE卫星轨道数据恢复地球重力场模型的方法

欧空局早期公布的时域法和空域法解算的GOCE模型均采用能量守恒法处理轨道数据,但恢复的长波重力场信号精度较低,而且GOCE卫星在两极存在数据空白,利用其观测数据恢复重力场模型是一个不适定问题,导致解算的模型带谐项精度较低,需进行正则化处理。本文分析了基于轨道数据恢复重力场模型的方法用于处理GOCE数据的精度,对最优正则化方法和参数的选择进行了研究。利用GOCE卫星2009-11-01—2010-01-31共92d的精密轨道数据,采用不依赖先验信息的能量守恒法、短弧积分法和平均加速度法恢复GOCE重力场模型,利用Tikhonov正则化技术处理病态问题。结果表明,平均加速度法恢复模型的精度最高,能量守恒法的精度最低,短弧积分法的精度稍差于平均加速度法。未来联合处理轨道和梯度数据时,建议采用平均加速度法或短弧积分法处理轨道数据,并且轨道数据可有效恢复120阶次左右的模型。Kaula正则化和SOT处理GOCE病态问题的效果最好,并且两者对应的最优正则化参数基本一致,但利用正则化技术不能完全抑制极空白问题的影响,需要联合GRACE等其他数据才能获得理想的结果。     

利用GOCE卫星轨道数据恢复地球重力场模型方法的分析

欧空局早期公布的时域法和空域法解算的GOCE模型均采用能量守恒法处理轨道数据, 但恢复的长波重力场信号精度较低, 而且GOCE卫星在两极存在数据空白, 利用其观测数据恢复重力场模型是一个不适定问题, 导致解算的模型带谐项精度较低, 需进行正则化处理。本文分析了基于轨道数据恢复重力场模型的方法用于处理GOCE数据的精度, 对最优正则化方法和参数的选择进行研究。利用GOCE卫星2009-11-01—2010-01-31共92 d的精密轨道数据, 采用不依赖先验信息的能量守恒法、短弧积分法和平均加速度法恢复GOCE重力场模型, 利用Tikhonov正则化技术处理病态问题。结果表明, 平均加速度法恢复模型的精度最高, 能量守恒法的精度最低, 短弧积分法的精度稍差于平均加速度法。未来联合处理轨道和梯度数据时, 建议采用平均加速度法或短弧积分法处理轨道数据, 并且轨道数据可有效恢复120阶次左右的模型。Kaula正则化和SOT处理GOCE病态问题的效果最好, 并且两者对应的最优正则化参数基本一致, 但利用正则化技术不能完全抑制极空白问题的影响, 需要联合GRACE等其他数据才能获得理想的结果。     

星载原子干涉技术用于地球重力场测量及其精度评估

重力梯度卫星GOCE通过搭载静电式重力梯度仪,将全球静态重力场恢复至200阶以上。目前GOCE卫星已结束寿命,亟须发展下一代更高分辨率的卫星重力梯度测量来完善200~360阶的全球静态重力场模型。原子干涉型的重力梯度测量在空间微重力环境下可获得较长的干涉时间,因此具有很高的星载测量精度,是下一代卫星重力梯度测量的候选技术之一。本文针对未来更高分辨率全球重力场测量的科学需求,提出了一种适用于空间微重力环境下的原子干涉重力梯度测量方案,其梯度测量噪声可低至0.85mE/Hz1/2。文中对不同类型的卫星重力梯度测量方案进行了重力场反演精度的对比评估,仿真结果表明,相比于现有静电式卫星重力梯度测量,原子干涉型的卫星重力梯度测量有望将重力场的恢复阶数提升至252~290阶,对应的累积大地水准面误差7~8cm,累积重力异常误差3×10-5 m/s2。     

利用能量守恒法和GOCE卫星轨道数据反演地球重力场模型

利用傅立叶级数拟合GOCE卫星的耗散能,解决了基于能量守恒法恢复GOCE重力场模型时耗散能的计算问题。采用Helmert-Wolf参数估计法统一求解位系数、能量常数和耗散能的傅立叶级数拟合参数,并采用消局部参数的最小二乘法求解位系数。该方法不需要任何初始值或参考模型,不需要采用差分方法处理能量常数,也不需要进行迭代计算。利用GOCE卫星2009-11-01~2010-02-12共103d的精密轨道数据反演了三组100阶次的重力场模型GOCE-ECP01S、GOCE-ECP02S和GOCE-ECP03S,并与EIGEN-5C、EIGENCHAMP05S和GOCO03S模型进行比较。结果表明,采用一阶傅立叶级数拟合GOCE卫星的耗散能效果最好,反演的GOCE-ECP01S模型精度最高,整体精度优于EIGEN-CHAMP05S,但较GOCO03S模型的精度偏低;在100阶次的大地水准面误差为±3.2cm,但由于极空白的影响,恢复模型的带谐项位系数精度偏低。     

GOCE卫星重力梯度数据反演重力场的滤波器设计与比较分析

地球重力场和海洋环流探测（gravity field and steady-state ocean circulation explorer,GOCE）卫星重力梯度数据有色噪声和低频系统误差的滤波处理是反演高精度地球重力场的一个关键问题。针对GOCE卫星重力梯度数据的滤波处理,基于移动平均（moving average,MA）方法和CPR(circle per revolution)经验参数方法设计了两类低频系统误差滤波器,并分别将这两类滤波器与基于自回归移动平均（auto-regressive and moving average,ARMA）模型设计的有色噪声滤波器组合起来形成级联滤波器。为了分析滤波器处理的实际效果,基于空域最小二乘法采用70 d的GOCE观测数据,并联合重力恢复与气候实验（gravity recovery and climate experiment,GRACE）数据分别反演了224阶次的重力场模型GOGR-MA(MA+ARMA级联滤波)和GOGR-CPR(CPR+ARMA级联滤波)。将反演模型与采用同期数据求解的第一代GOCE系列模型及GOCE和GRACE联合模...     

联合GOCE卫星数据和GRACE法方程确定SWJTU-GOGR01S全球重力场模型

联合地球重力场和海洋环流探测器（Gravity Field and Steady-State Ocean Circulation Explorer,GOCE）和重力恢复与气候实验（Gravity Recovery and Climate Experiment,GRACE）卫星观测数据确定全球静态重力场模型是当前大地测量学的研究热点之一。联合近3 a的GOCE卫星梯度数据和7 a左右的GRACE星间距离变率数据计算的ITG-GRACE2010S模型的法方程恢复了210阶次的重力场模型SWJTU-GOGR01S。采用带通数字滤波方法处理GOCE卫星的4个高精度梯度观测分量,利用梯度数据恢复重力场模型的观测方程直接建立在梯度仪坐标系中,可以避免坐标转换过程中高精度的梯度观测分量受低精度分量的影响;联合法方程解的最优权采用方差分量估计迭代计算,GOCE数据的两极空白引起的病态问题采用Kaula正则化方法进行约束。基于EIGEN-6C2模型和北美地区的GPS水准网观测数据,对SWJTU-GOGR01S模型进行内外符合精度分析,结果表明,SWJTU-GOGR01S模型在210阶次的大地水准面误差和累计误差分别为1.3 cm和5.7 cm,精度与欧洲空间局公布的第四代时域法模型相当,略优于GOCO02S和GOCO03S模型的精度。     

GOCE卫星重力梯度观测值的时变重力场变化改正及影响分析

GOCE卫星重力梯度观测值为高阶静态重力场反演提供了重要的数据支撑,但其在使用前需考虑扣除时变重力场变化的影响.本文研究了GOCE卫星重力梯度观测值的时变重力场变化改正方法,更新了ESA标准和背景模型,以更好地扣除时变重力场变化的影响,自主实现了由GOCE卫星Level1b重力梯度数据直接进行重力场反演.本文通过3种时...     

利用GOCE模拟观测反演重力场的Torus法

在介绍Torus方法反演地球重力场模型的基本原理和方法的基础上,基于圆环面上均匀分布的卫星引力梯度模拟观测值解算了200阶次的地球重力场模型,在无误差情况下,Torus方法解算模型的阶误差RMS小于10-16,验证了该方法的严密性。利用61dGOCE卫星轨道上无误差的模拟引力梯度观测值解算了200阶次的地球重力场模型,分析了格网化误差、极空白对解算精度的影响,迭代3次后,在不考虑低次系数情况下,模型的大地水准面阶误差和累积误差均较小,最大值仅为0.022mm和0.099mm。在沿轨卫星引力梯度模拟数据中加入5mE/Hz1/2的白噪声,基于Torus方法和空域最小二乘法解算了200阶次的地球重力场模型,Torus方法的精度略低于空域最小二乘法的精度,在不考虑低次项的情况下,两种方法解算模型的大地水准面阶误差最大值分别为1.58cm和1.45cm,累积误差最大值分别为6.37cm和5.55cm。但由于采用了二维快速傅里叶技术和块对角最小二乘法,极大地提高了计算效率。本文数值结果说明Torus方法是一种独立有效的方法,可用于GOCE任务海量卫星引力梯度观测值反演重力场的快速解算。     

联合GOCE卫星轨道和重力梯度数据严密求解重力场的模拟研究

论述了联合卫星轨道和重力梯度数据严密求解重力场的方法及数据处理方案,研究了GOCE重力场反演中有色噪声的AR去相关滤波、病态法方程的Kaula正则化和观测值最优加权的方差分量估计等关键问题。模拟结果表明:①极空白问题会降低法方程求解的稳定性,导致低次位系数的求解精度较低,而Kaula正则化可有效用于GOCE病态法方程的求解,并得到合理稳定的解;②重力梯度有色噪声会降低GOCE重力场求解的整体精度,特别是对低阶位系数的影响最为明显,而AR去相关滤波法可有效处理有色噪声,但解算结果仍含有低频误差;③方差分量估计可有效确定SST和SGG两类观测值的最优权比,并且有色噪声造成的低频误差经过联合求解后得到了抑制;④利用30d、5s采样的GOCE模拟数据恢复200阶次的重力场模型,其大地水准面和重力异常精度在纬度±83°范围内分别为±3.81cm和±1.056mGal。     

一种计算GOCE卫星加速度的实用算法

引入了一种由卫星轨道数据计算瞬时加速度的实用数值微分算法,分别采用模拟卫星轨道数据、地球重力场和稳态海洋环流探测(GOCE)卫星任务发布的几何法轨道数据验证该算法的有效性。数值模拟结果表明,采用8阶9点10s采样间隔的移动窗口多项式数值微分算法计算加速度的精度最佳,71天实测几何法轨道数据采用相同算法恢复的100阶重力场模型GOCE-PKI-AA与EIGEN-5C模型在80阶以前具有良好的一致性,GOCE-PKI-AA模型整体精度优于EIGEN-CHAMP03S,低阶(小于15)部分较欧空局发布的第一代GOCE重力场模型GO-CONS-GCF-2-TIM-R1更接近于EIGEN-5C模型。     

GOCE卫星测量恢复地球重力场模型的理论与方法

研究了GOCE卫星测量恢复地球重力场模型的理论与方法。论文的主要工作和创新点有： (1) 建立了扰动重力梯度张量各分量没有奇异性的详细计算模型,解决了重力梯度张量Txx分量在两极地区计算的奇异性难题。 (2) 系统研究了卫星重力梯度数据向下延拓的解析法、泊松积分迭代法和卫星重力梯度数据格网化的移动平均法、反距离加权法、普通克里金法,建立了相应的数学模型,导出了相应的计算公式,并采用“直接法”和“移去-恢复法”两种方案对其向下延拓和格网化效果进行了测试。 (3) 分析了能量守恒方程中各项误差对沿轨扰动位计算结果的影响,建立了利用GOCE模拟数据确定地球重力场的最小二乘直接法、调和分析法、最小二乘配置法的实用数学模型,并做了大量的模拟计算。 (4) 建立了利用扰动引力梯度张量各单分量和组合分量确定地球重力场的最小二乘直接法去奇异性计算模型;推导了利用扰动引力梯度张量单分量和组合分量解算地球重力场的调和分析法模型;进一步推导了扰动引力梯度张量各个分量之间的自协方差和互协方差函数及其与引力位系数之间协方差函数的具体计算公式。 (5) 推导了利用不同类型重力测量数据确定地球重力场的联合平差法数学模型,介绍并分析了模型中各类数据最优定权的参数协方差法和方差分量估计法。 (6) 论述了谱组合法的基本原理,给出了多种类型重力测量数据联合处理的谱权及谱组合的通用表达式,基于调和分析方法推导了SST+SGG、SST+SGG＋Δg和SST+SGG＋Δg+N恢复地球重力场模型的谱组合公式及对应谱权的具体形式。 (7) 推导了利用迭代法联合不同类型重力测量数据反演地球重力场模型的基本原理公式,并给出了其具体实现步骤。 (8) 分析并计算了重力卫星轨道高度、卫星星间距离和卫星轨道倾角的设计指标;讨论了双星轨道长半轴的一致性要求、双星姿态俯仰角的控制要求以及双星编队保持机动的时间间隔要求。 (9) 确定了KBR系统的星间距离、星间距离变化率和星间加速度的精度指标;设计了星载GPS系统的卫星轨道位置和速度以及加速度计测量的精度指标;计算了加速度计检验质量质心到卫星质心的调整距离精度指标;分析了恒星敏感器的姿态角测量精度和稳定度;计算了参考重力场模型对于累计大地水准面精度和积分卫星轨道的影响。 (10) 研制了一套利用卫星重力测量数据反演地球重力场模型的软件平台,可对卫星重力测量数据处理及其精度评估提供一些基本方法,并为我国卫星重力测量系统的总体战技指标和主要有效载荷技术指标的量化分析、论证提供理论和技术支持,为我国未来的卫星重力测量系统提供可能的积累和参考。     

由GOCE引力梯度张量不变量确定卫星重力模型的半解析法

提出了利用半解析方法解算重力场与稳态海洋环流探测器(gravity field and steady-state ocean circulation explorer,GOCE)引力梯度张量不变量观测值来恢复全球卫星重力模型的方法,该方法比最小二乘方法效率高,同样可给出模型的验后方差。推导了半解析法快速解算引力梯度张量不变量的理论公式,给出了该方法利用卫星重力梯度观测数据快速求解重力场模型的基本步骤。利用89.5°倾角圆形严格重复轨道上的模拟梯度观测值验证了半解析方法用于张量不变量解算的理论严密性,并利用GOCE任务61d梯度仪坐标系(gradiometer reference frame,GRF)下无噪声和有噪声的I2和Vzz模拟观测值进行了模型反演和结果分析。结果表明,利用I_2观测值求解模型的精度略优于仅用Vzz分量解算模型的精度。     

利用先验重力场模型对GOCE卫星重力梯度观测值进行校准分析

GOCE卫星任务搭载了高灵敏度的重力梯度仪,其观测值用于恢复高精度高分辨率的地球重力场。本文利用EIGEN-5C、EGM2008、GOTIM3、GGM03S高精度全球重力场模型,确定了GOCE引力梯度张量的对角分量观测值(Vxx、Vyy、Vzz)的校准参数,分析了比例因子的稳定性,并讨论了相同模型不同阶次、同阶次不同模型以及是否估计漂移参数对比例因子、偏差参数及校准观测值的影响。研究表明比例因子的稳定性在10-4的量级,利用250阶的EIGEN-5C模型和EGM2008模型校准得到观测值的差异小于10-4 E,远远小于观测误差,以1d为周期估计校准参数时,是否估计漂移对校准结果的影响达到0.4E。同时,校准前后观测值差异的频谱说明校准过程主要影响Vxx、Vyy、Vzz观测值的低频部分,即来自先验重力场模型的中低(150)阶次,考虑到GOCE引力梯度的观测频带,校准后的观测值可用于恢复中高频的重力场信号。     

扰动重力梯度的球冠谐分析建模

从球冠谐理论出发,详细推导了球冠坐标系下扰动重力梯度的无奇异性计算公式。基于Tikhonov正则化方法,利用GOCE卫星实际观测数据解算局部重力场球冠谐模型。数值计算表明,基于扰动重力梯度的球冠谐分析建模方法能够有效地恢复局部重力场中的短波信号,与GO_CONS_GCF_2_DIR_R5模型的差异在±0.3×10~(-5) m/s~2水平。     

利用引力梯度数据计算径向梯度的优化方法

根据重力梯度观测各分量的方差及协方差信息,提出了利用GOCE梯度数据计算径向重力梯度的优化方法。首先给出了径向重力梯度的计算方法,并深入分析了误差传播规律,通过建立相应的条件极值问题,给出了计算径向重力梯度最优组合因子的方法;通过模拟数据验证了本文所提出的优化因子的优越性。实际数据计算表明:相对于传统方法,采用优化组合因子可使反演所得引力位模型的累积大地水准面精度在250阶时提高约2 cm。由于径向重力梯度不仅可以用于地球引力场模型的求解,也可直接应用于地球物理问题的讨论,因此本文所提出的优化方法也可对部分地球动力学问题的讨论提供方便。     

利用先验重力场模型求定GOCE卫星重力梯度仪校准参数

重力梯度仪校准参数的确定是GOCE重力梯度观测数据处理的关键环节。本文对GOCE卫星重力梯度观测值中的时变信号与粗差进行了分析,利用高精度全球重力场模型,确定了GOCE重力梯度观测值各分量的尺度因子与偏差,并对校准结果进行了精度评定。结果表明,在测量带宽内,海潮对重力梯度观测值影响在mE量级,与重力梯度仪的精度水平相当,陆地水等非潮汐重力场时变信号略小于海潮,量级约为10^-4E;各分量重力梯度观测值的粗差比例均大于0.2%;除EGM96模型外的其他模型对GOCE重力梯度仪进行校准后,V_xx、V_yy、V_zz、V_yz分量上尺度因子的稳定性均在10^-4量级,V_xz分量能达到10^-5量级,V_xy分量为10^-2量级,这与梯度观测值各分量的精度水平一致。