基于GOCE卫星的重力场模型反演及应用

发布时间：2018-05-21 06:28

本文选题：GOCE + 地球重力场模型　；参考：《西南交通大学》2014年博士论文

【摘要】：地球重力场一直都是大地测量学中的研究核心之一。近年来,随着CHAMP (CHAllenging Minisatellite Payload)、GRACE (Gravity Recovery And Climate Experiment)和GOCE (Gravity field and steady-state Ocean Circulation Explorer)卫星的相继发射及数据免费释放,基于低轨卫星重力测量的研究成为了大地测量学中的研究热点。GOCE卫星是首颗采用SST-HL (Satellite-to-Satellite Tracking in High-Low mode)技术和SGG(Satellite Gravity Gradiometry)技术的低轨重力卫星,其搭载的GPS (Global Position System)接收机可采集SST-HL数据,静电重力梯度仪可采集GG (Gravity Gradient)数据。SST-HL数据可恢复地球重力场的中长波信息,GG数据可恢复地球重力场的中短波信息。因此,联合SST-HL和GG数据可恢复高精度的静态地球重力场模型。本文主要内容是研究利用GOCE卫星观测数据恢复静态地球重力场模型。在前人研究成果的基础上,导出了一种联合GOCE轨道数据和重力梯度数据的时域解法,并独立研制了一款基于GOCE卫星反演地球重力场模型的软件系统。提出顾及地球重力场模型道路工程勘测的3种实施方案,并且将利用GOCE数据恢复的重力场模型应用到道路工程勘测的模拟试算中。本文的主要内容及贡献如下：1.给出了GOCE卫星重力测量涉及到坐标系统的定义及其详细转换公式。比较了勒让德系列函数的递推计算公式,优化了勒让德积分函数的递推计算,采用向后递推法解决积分计算在极区不稳定的现象。为了避免重力梯度分量在极区出现奇异,推导了广义勒让德函数PS1和PS2的递推计算公式。推导了地固系下模拟引力梯度数据的详细公式。基于GOCE卫星轨道、行星星历文件DE405和海潮模型FES2004等计算出日月引力、固体潮、海潮和极潮等摄动力。引入了具有高度优化和并行性能的Intel MKL函数库。2.对GOCE卫星PKI轨道数据进行了预处理,利用能量守恒法、短弧长积分法和平均加速度法反演出3个模型,并比较了3种方法的优缺点。对GOCE轨道数据反演能力的研究表明,GOCE轨道数据的最大恢复能力大约在120阶次左右。提出利用GOCE卫星轨道数据探测地球重力场的时变信息,并且对扇形滤波进行了改进。分析了GOCE重力场模型球谐位系数残差之间的相关性,并与GRACE RL05模型球谐位系数残差的相关性进行了比较。利用近1年的GOCE轨道数据计算出南极冰盖质量的变化情况,探测到南极某点(-75°,250°)处的质量变化约为-11.30cm/a,并用GRACE RL05重力场模型验证了探测结果的正确性。3.对GOCE卫星重力梯度数据进行了预处理,基于重力梯度的测量原理和误差特性设计出3种重力梯度数据的滤波方法,并对3种滤波方法的滤波效果进行比较,选取其中一种效果最好的方法作为重力梯度数据的固定滤波方法。分别利用空域法和时域法基于重力梯度数据恢复了重力场模型。为解决因GOCE卫星两极空白而带来的病态问题,采用了正则化处理,并比较了标准正则化和球冠正则化。结果显示,球冠正则化效果更好。推导出了用于球冠正则化法积分计算的4类三角函数的原函数。导出了一种联合GOCE轨道数据和重力梯度数据恢复重力场模型的时域解法。利用近1年的GOCE卫星观测数据反演出一个200阶次的静态重力场模型SWJTU-GO2013,该模型的空间分辨率约为100km(半波长),大地水准面误差约为±20.32cm,与ICGEM公布的GOCE卫星第3代时域法模型进行比较,具有同等级的精度。4.针对GOCE卫星观测数据,编写了一款恢复重力场模型的软件系统SWJTU-GOCE,编写该软件系统主要使用了C#语言和Fortran语言。为使耗时巨大的重力场恢复计算在一般的PC电脑上也能开展,在Fortran代码中调用了具有高度优化和并行特性的Intel MKL函数。在可视化方面,软件可自动生成Matlab或GMT等软件的脚本程序,运行脚本程序即可得到可视化结果。软件模块涵盖了从数据预处理到精度评定全部流程。5.详细分析了目前道路工程勘测所存在的问题：GPS系统的参考基准为法线,全站仪系统的参考基准为垂线,两者混合使用时存在明显的系统误差。分别基于勘测单位和设计单位的角度,在理论上提出顾及地球重力场模型的3种道路工程勘测方案。推导出利用GPS平面坐标恢复道路工程渐变坐标系的一整套公式和精度评定公式。最后,将利用GOCE数据恢复的模型SWJTU-GO2013应用到试验计算中,并用EGM2008的计算结果进行验证。
[Abstract]:The earth gravity field has always been one of the core of Geodetic Research. In recent years, with the CHAMP (CHAllenging Minisatellite Payload), GRACE (Gravity Recovery And Climate Experiment) and GOCE (Gravity) satellites are released and free of free data, based on low rail. .GOCE satellite is the first low rail gravity satellite using SST-HL (Satellite-to-Satellite Tracking in High-Low mode) technology and SGG (Satellite Gravity Gradiometry) technology. The gravity gradiometer can collect the GG (Gravity Gradient) data.SST-HL data to restore the medium and long wave information of the earth's gravity field, and the GG data can restore the medium and short wave information of the earth's gravity field. Therefore, the combined SST-HL and GG data can restore the high precision static earth gravity field model. The main content of this paper is to study the recovery of the GOCE satellite observation data. The static earth gravity field model. On the basis of the previous research results, a time-domain method of solving the combined GOCE orbit data and gravity gradient data is derived, and a software system based on the GOCE satellite inversion of the earth gravity field model is independently developed. 3 implementation schemes which take account of the earth gravity field model road engineering survey are proposed and will be put forward. The gravity field model recovered by GOCE data is applied to the simulated trial calculation of road engineering. The main contents and contributions of this paper are as follows: 1. the definition of GOCE satellite gravity measurement and its detailed conversion formula are given. The recurrence calculation formula of Legendre series function is compared, and the Legendre integral function is optimized. In order to avoid the singularity of the gravity gradient component in the polar region, the recurrence calculation formula of the generalized Legendre function PS1 and PS2 is derived. The detailed common formula of the simulated gravitational gradient data in the earth solid system is derived. Based on the GOCE satellite orbit, the ephemeris file DE40 is derived. 5 and the sea tide model FES2004 etc. calculate the perturbation of the sun and moon gravity, the solid tide, the sea tide and the extreme tide. The Intel MKL function library.2. with high optimization and parallel performance is introduced to the GOCE satellite PKI orbit data, and the energy conservation method, the short arc length integral method and the average acceleration method are used to counter the 3 models, and 3 kinds of formulas are compared. The research on the inversion ability of GOCE track data shows that the maximum recovery ability of the GOCE orbit data is about 120 order. The time variation information of the earth gravity field is detected by using GOCE satellite orbit data, and the fan filter is improved. The correlation between the residual difference of the spherical harmonic position coefficient of the GOCE heavy force field model is analyzed. The correlation between the spherical harmonic position coefficient residuals of the GRACE RL05 model is compared with that of the GRACE RL05 model. The changes in the mass of the Antarctic ice cover are calculated by using the GOCE orbital data of nearly 1 years. The mass change at a point (-75, 250) is about -11.30cm/a, and the correctness of the detection results is verified by the GRACE RL05 gravity field model,.3. to the GOCE satellite The gravity gradient data is pretreated. Based on the measurement principle and error characteristics of gravity gradient, the filtering methods of 3 gravity gradient data are designed, and the filtering effects of the 3 filtering methods are compared. One of the best effective methods is selected as the fixed filtering method of gravity gradient data. The gravity field model is restored based on the gravity gradient data. In order to solve the ill conditioned problem caused by the gap between the two poles of the GOCE satellite, the regularization treatment is adopted and the standard regularization and the canonical regularization are compared. The results show that the regularization effect of the spherical crown is better. The original function of the 4 kinds of trigonometric functions for the integral calculation of the spherical crown regularization method is derived. A time domain solution of a combined GOCE orbit data and gravity gradient data restoration gravity field model is derived. A 200 order static gravity field model SWJTU-GO2013 with a spatial resolution of about 100km (half wavelength), the geooid error is about + 20.32cm, and the ICGEM is published with ICGEM. The GOCE satellite third generation time domain method is compared, with the same class of precision.4. for GOCE satellite observation data, a software system SWJTU-GOCE for restoring gravity field model is written. The software system mainly uses the C# language and the Fortran language. In order to make the time-consuming gravity field recovery calculation on the ordinary PC computer, it can also be used. In the Fortran code, the Intel MKL function with high optimization and parallel characteristics is called. In visualization, the software can automatically generate scripts of Matlab or GMT software, and the script program can get visual results. The software module covers the detailed analysis of the whole process from data preprocessing to precision evaluation.5.. The problems existed in the survey of the former Road Engineering: the reference datum of the GPS system is the normal line, the reference datum of the total station system is perpendicular to the vertical line. There are obvious systematic errors in the mixed use of the two. Based on the angle of the survey unit and the design unit, 3 road engineering exploration schemes which take account of the earth reforce field model are proposed in theory. A complete set of formula and accuracy evaluation formula for the restoration of road engineering gradient coordinate system using GPS plane coordinates. Finally, the model SWJTU-GO2013 using the GOCE data recovery is applied to the test calculation, and the results of the EGM2008 calculation are verified.
【学位授予单位】：西南交通大学
【学位级别】：博士
【学位授予年份】：2014
【分类号】：P223.0

【参考文献】