扰动重力边值问题与实际数据处理的研究

讨论了应用GPS测量重力点的大地高 ,从而获得地面点的扰动重力 (或称纯重力异常 )值这样新的重力边值条件 ,如何从理论与实际上确定重力场的问题。     

样条函数拟合GPS高程探讨

主要讨论了运用三次样条函数曲线拟合法来拟合GPS点的高程,并与多项式曲线拟合法比较,通过实例验证了此方法的可行性及其所能达到的精度,得出了一些有益的结论,并给出了一些建议。     

弹道扰动引力的有限元逼近算法

为了克服多项式逼近弹道扰动引力的缺点，根据有限元插值的原理，采用了对弹道周围空间区域进行有限元剖分的方法，利用剖分单元各顶点的扰动引力分量内插出弹道点对应的扰动引力分量值。结果表明，文中提出的逼近算法能够快速精确可靠地逼近弹道扰动引力，是一种具有应用价值的方法。     

Generation and study of satellite gravity over Gujarat,India and their possible correlation with earthquake occurrences

High-resolution satellite gravity data have been generated and utilized to infer subsurface geological structures in the area of devastating earthquake that struck the Bhuj region in Gujarat on 26 January 2001. Latitudinal gravity profiles have been generated in the Bhuj, Anjar and IBF regions across the epicentres (23.5° N, 69.8° E/M_w 7.0 in 2001; 23.2° N, 70° E/M_w 7.0 in 1956; 24.2° N, 69.2° E/M_w 7.8 in 1819). Substantial differences in gravity anomaly patterns as high as 37 mGal could be observed existing near the epicentre regions. These gravitational differences might have caused due to the plate tectonic processes and due to the changes in densities of different lithospheric zones/sedimentary layers. Temporal variations of the satellite-derived gravity and their probable relations with already occurred major earthquakes in this region have been studied. Hence we conclude that drastic changes in gravity anomalies can be considered as a precursor for occurrences of substantially large earthquakes.     

削弱残余误差对多波束测深综合影响的方法研究

提出了一种基于地形变化长波项与短波项相结合的残余误差综合削弱方法。实验验证了该方法的有效性。     

球谐函数变换快速计算扰动引力

在球谐函数变换基础上,利用新极下轨道的特殊性,在新坐标下引入Clenshaw求和计算轨道扰动引力。从理论上对比分析了传统方法、球谐函数变换方法和改进方法的计算速度和存储模型需要的物理空间。模拟试验分别采用3种方法计算了一段轨道的扰动引力,试验结果表明,改进的球谐函数变换方法比传统球谐函数变换方法计算速度可提高100倍,数据存储量仅占传统方法的3%。     

Accuracy analysis for SST gravity field model in China

Taking China as the region for test the potential of the new satellite gravity technique, satelliteto-satellite tracking for improving the accuracy of regional gravity field model is studied. With WDM94 as reference, the gravity anomaly residuals of three models, the latest two GRACE global gravity field model (EIGEN_GRACE02S, GGM02S) and EGM96, are computed and compared. The causes for the differences among the residuals of the three models are discussed. The comparison between the residuals shows that in the selected region, EIGEN_GRACE02S or GGM02S is better than EGM96 in lower degree part (less than 110 degree). Additionally, through the analysis of the model gravity anomaly residuals, it is found that some systematic errors with periodical properties exist in the higher degree part of EIGEN and GGM models, the results can also be taken as references in the validation of the SST gravity data.     

应用Slepian局部谱方法解算中国大陆重力场球谐模型

根据经典的球谐函数方法,为满足正交化要求,观测数据需要覆盖整个球面,而对于地表局部测量数据,则无法应用球谐方法解算重力场模型。针对此问题,采用Slepian局部谱分析方法解算中国大陆范围内的实测重力场变化数据,并以GOCE卫星球谐函数解作为已知模型,评估由于实际陆地重力测点的非均匀分布对球谐函数解的误差影响。通过计算多个阶次中国大陆局部范围的Slepian基函数分布;采用GOCE卫星获得重力场模型的前72阶球谐系数作为已知结果,评价实际测点非均匀分布的解算有效性,并针对中国大陆地区采用Slepian基函数进行解算,通过模型对比选择最优截段项数;针对2005—2008年中国大陆地区流动重力测量获得的重力场变化信号进行解算,获得了72阶重力场变化模型。     

多个超高阶重力场模型精度分析

分析了国际和国内公布的EGM2008、GECO、Eigen6c4、UGM2005、DQM2000d和DQM2006等超高阶重力场模型的内符合精度;利用中国境内的实测GPS/水准数据对各类模型进行了外符合精度检验;结合不同的方法对系统差进行了处理;最后基于试验区数据讨论了不同模型之间频谱组合方法对高程异常计算结果的改善情况。实验表明,Eigen-6C4模型具有更高的内符合精度,国内公布的UGM2005和DQM2006两种模型精度相当;本文提出的十参数法相比于传统四参数方法精度获得了4.6%~22.20%的提升;GECO模型和Eigen-6C4模型在实验区内外符合精度较好;不同模型之间的频谱组合可以在一定程度上提高试验区内模型确定高程异常的精度。     

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

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