局部大地水准面精化系统软件的设计与实现

给出了局部大地水准面精化系统的总体结构以及各分系统的详细设计与实现方法,开发了旨在对相关重力测量数据实现自动化处理的局部大地水准面精化系统工具软件,并应用于某区域实际大地水准面的精确计算。实践证明,该系统能有效提高相关数据的处理效率,有助于局部大地水准面精化过程的规范化。     

利用空中平均重力异常确定区域大地水准面

提出了直接利用空中平均重力异常计算区域大地水准面的方法。模拟计算的结果表明, 该方法与传统的利用地面平均空间重力异常确定的大地水准面精度相当, 但其显著优点是勿需空中重力异常的向下解析延拓, 从而可以避免延拓误差对大地水准面精化的影响。     

区域似大地水准面再精化精度影响因素分析

介绍区域似大地水准面再精化方法,以江苏省域内不同精度的重力似大地水准面及试验区全球导航卫星系统(global navigation satellite system,GNSS)/水准数据为例,计算分析了GNSS/水准点数、初始重力似大地水准面精度和内插拟合方法对再精化水准面的影响。结果表明,初始重力似大地水准面精度优于10 cm时,其再精化水准面结果精度与初始重力水准面精度关系不大,但其离散程度和初始重力水准面精度相关;平面/二次拟合与反距离平方加权内插结合使用效果更优,三次拟合与Shepard内插结合使用效果最优。结论对今后区域似大地水准面再精化工作中方案的选择、成本的节约等方面有一定的指导意义。     

似大地水准面的误差分析与抑制技术

大地水准面误差分析与精度评定是局部重力场逼近技术的重要组成部分,是大地水准面精化工程外业方案优化、算法设计和工程质量评价的基本依据。本文分别从地面重力数据误差和局部重力场算法两个方面,分析cm级大地水准面误差的影响特性,提出重力数据误差与大地水准面精度之间普遍适用的规律,推荐一种GPS水准和地面重力数据联合平差的精度评定方法,结合实例和模拟计算分析,介绍大地水准面误差分析与误差抑制方法。     

GPS水准和重力数据联合平差技术及其性能分析

本文推荐一种适合高精度大地水准面精化的GPS水准和重力数据联合平差技术。这种技术以重力场积分公式为约束条件,在统计框架中,使GPS水准高程异常和其他重力场元(平差值)满足重力场边值条件,进而利用条件平差技术构造大地水准面的统计解。实测和模拟计算结果表明,利用本文推荐的GPS水准和重力数据联合平差技术,可以较好地解决目前高精度大地水准面精化中地面重力数据精度分布存在较大差异,重力数据、GPS水准数据存在多种系统偏差,以及局部重力场积分的“远区效应”等一系列难题。     

最新中国陆地数字高程基准模型:重力似大地水准面CNGG2011

本文回顾了近20年国内外国家局部大地水准面模型研究的概况和发展背景,采用Stokes-Helmert方法,计算了一个新的2′×2′中国重力和1985国家高程基准似大地水准面数值模型（CNGG2011）,采用了1百万余陆地重力数据和SRTM 7″.5×7″.5地形高数据,以及649个B级GPS水准点数据。CNGG2011平均精度为±0.13m,东部地区±0.07m,西部地区±0.14m。各省区局部似大地水准面平均精度为±0.06m,东部为±0.05m,西部为±0.11m。西藏精度为±0.22m。本文还讨论了重力大地水准面与GPS水准的关系,提出了今后进一步精化我国高程基准大地水准面模型的构想。     

地球及其形状非对称性的重力学研究

依据卫星重力、大地水准面观测结果,对表示地球形状的大地水准面异常进行了不同阶次的球谐函数计算和分析,发现2~6阶大地水准面异常表现出了地球双重非对称的基本形状,而高于6阶的大地水准面异常只表示了地球形态的局部变化特征。另外,采用阻尼最小二乘方法,利用全球地震层析成像资料和大地水准面异常资料联合反演了三维全球地幔密度异常。结果显示了地幔密度异常不仅在横向上,而且在纵向上也存在着明显的不均匀性。对比分析低阶大地水准面异常和地幔密度异常结果表明,地球形状非对称性主要是由下地幔的物质密度不均匀引起的。     

水体对精化局部似大地水准面的影响

在精化局部似大地水准面时，由于对大地水准面精度的要求不高或者精化区域的水下地形数据缺乏等因素，往往很少考虑湖泊水体对局部似大地水准面的影响。本文将通过分析水体对局部地形改正和地形均衡改正的影响，然后通过实例计算得出了湖泊水体对精化局部似大地水准面的影响。     

基于EGM2008重力场模型的松原灌区大地水准面精化

以松原灌区GPS控制网为例,采用高程异常直接拟合和插值法,以及基于EGM2008重力场模型和移去-计算-恢复技术的残差二次曲面拟合和插值法,精化了测区大地水准面。首先简要介绍了研究区数据,进而探讨了高程异常的获取,在此基础上,详细介绍了二次曲面拟合的大地水准面精化和空间插值法大地水准面精化方法,结果表明,EGM2008重力场模型精度较高,对精化局域和国家大地水准面具有重要意义。     

精化区域似大地水准面应顾及的有关问题

简述精化区域似大地水准面的计算方法;分析误差的主要来源;讨论已知平均重力异常栅格分辨率,求高程异常与布测GPS水准点距的关系;论述在精化区域似大地水准面项目设计中应采用统一的测绘基准、现势性好的高程起算数据,要全面规划和建设好地方测绘基础控制网等有关方面的问题。     

应用双输入单输出法确定局部似大地水准面的实例研究

利用了双输入单输出法，融合处理了我国某地区的重力异常和地形资料两类数据，结合WDM94地球重力场模型和63个高精度GPS水准数据，计算了该区域的似大地水准面。     

利用重力异常数据构建重力梯度场的方法比较

为实现大范围、高精度基准重力梯度数据库的构建,考虑到重力梯度场对地形质量的敏感效应,一般利用恒密度数字高程模型来求取重力梯度值,从而忽略了地形密度变化以及水准面以下密度异常对重力梯度的影响。根据重力位理论中求解边值问题的数值应用方法,直接利用重力异常数据求取重力梯度场,弥补了密度变化和密度异常在重力梯度上的反映。根据模型算例和实测重力异常数据求取了剖面重力梯度值,结果表明,限于重力数据空间分辨率的影响,利用重力异常数据可恢复中长波段重力梯度场。该方法与地形数据求取重力梯度和卫星重力梯度测量等方法技术相结合,对重力梯度数据库的建设具有实际应用价值。     

海洋重力似大地水准面与区域测高似大地水准面的拟合问题

研究了将陆地重力似大地水准面与GPS水；住似大地水准面拟合的处理方法推广到海洋的问题．首先从理论上证明了当存在海面地形．则海洋大地水准面与似大地水准面不重合．导出了在海洋上大地水；住面差距与高程异常之间差值的公式．由此给出了求定平均海面相对于区域高程基准的正常高以及测高似大地水准面的计算公式。由于测高平均海面与GPS大地高有相近的精度．提出了将海洋重力似大地水准面与区域测高似大地水准面拟合的处理方法．并利用当前最新的海面地形模型和测高平均海面模型做了数值估计。     

组合法精化胶州市大地水准面

本文利用全球重力位模型、胶州市地面重力观测数据、胶州市GPS水准数据和数字地面模型(DTM),采用组合法应用移去-恢复技术计算剩余大地水准面,并与地球位模型计算的高程异常进行拟合,得到该地区重力似大地水准面,再和布测、计算得到的GPS/水准所构成的几何大地水准面拟合,利用多项式拟合完成系统改正,获得最终的大地水准面结果及相关的精度信息。     

Geoid and high resolution sea surface topography modelling in the mediterranean from gravimetry,altimetry and GOCE data: evaluation by simulation

The determination of local geoid models has traditionally been carried out on land and at sea using gravity anomaly and satellite altimetry data, while it will be aided by the data expected from satellite missions such as those from the Gravity field and steady-state ocean circulation explorer (GOCE). To assess the performance of heterogeneous data combination to local geoid determination, simulated data for the central Mediterranean Sea are analyzed. These data include marine and land gravity anomalies, altimetric sea surface heights, and GOCE observations processed with the space-wise approach. A spectral analysis of the aforementioned data shows their complementary character. GOCE data cover long wavelengths and account for the lack of such information from gravity anomalies. This is exploited for the estimation of local covariance function models, where it is seen that models computed with GOCE data and gravity anomaly empirical covariance functions perform better than models computed without GOCE data. The geoid is estimated by different data combinations and the results show that GOCE data improve the solutions for areas covered poorly with other data types, while also accounting for any long wavelength errors of the adopted reference model that exist even when the ground gravity data are dense. At sea, the altimetric data provide the dominant geoid information. However, the geoid accuracy is sensitive to orbit calibration errors and unmodeled sea surface topography (SST) effects. If such effects are present, the combination of GOCE and gravity anomaly data can improve the geoid accuracy. The present work also presents results from simulations for the recovery of the stationary SST, which show that the combination of geoid heights obtained from a spherical harmonic geopotential model derived from GOCE with satellite altimetry data can provide SST models with some centimeters of error. However, combining data from GOCE with gravity anomalies in a collocation approach can result in the estimation of a higher resolution geoid, more suitable for high resolution mean dynamic SST modeling. Such simulations can be performed toward the development and evaluation of SST recovery methods.     

利用Abel-Poisson径向基函数模型化局部重力场

多种类型高分辨率重力场数据的不断增加,使得在局部范围内精化重力场模型成为了可能。本文采用Abel-Poisson核将重力场量表示成有限个径向基函数线性求和的形式,对局部区域的多种重力场数据进行联合建模。为了提高运算速度,运用了基于自适应精化格网算法的最小均方根误差准则(RMS)来求解径向基函数平均带宽。以南海核心地区为例,联合两种不同类型、不同分辨率的重力场资料(大地水准面起伏6'×6'、重力异常2'×2'),构建了局部区域高分辨率的重力场模型。所建模型表示的重力场参量达到了2'×2'的分辨率,对原始的重力异常数据(2'×2')拟合的符合程度达到±0.8×10-5m/s2。结果表明,利用径向基函数方法进行局部重力场建模,避免了球谐函数建模收敛慢的问题,有效提高了模型表示重力场的分辨率。     

Residual geoid and free-air gravity over the indian offshore from ers-1 high resolution altimeter data

In this study, ERS-1 altimeter data over the Indian offshore have been processed for deriving marine geoid and gravity. Processing of altimeter data involves corrections for various atmospheric and oceanographic effects, stacking and averaging of repeat passes, cross-over correction, removal of deeper earth and bathymetric effects, spectral analyses and conversion of geoid into free-air gravity anomaly. Methods for generation of residual geoid and free-air gravity anomaly using high resolution ERS-1 168 day repeat altimeter data were developed. High resolution detailed geoid maps, gravity anomaly and their spectral components have been generated over the Indian offshore using ERS-I altimeter data and ARCGIS system. A number of known megastructures over the study area have been successfully interpreted e.g. Bombay High, Saurastra platform, 90° east ridge etc. from these maps.     

Ellipsoidal geoid computation

Modern geoid computation uses a global gravity model, such as EGM96, as a third component in a remove–restore process. The classical approach uses only two: the reference ellipsoid and a geometrical model representing the topography. The rationale for all three components is reviewed, drawing attention to the much smaller precision now needed when transforming residual gravity anomalies. It is shown that all ellipsoidal effects needed for geoid computation with millimetric accuracy are automatically included provided that the free air anomaly and geoid are calculated correctly from the global model. Both must be consistent with an ellipsoidal Earth and with the treatment of observed gravity data. Further ellipsoidal corrections are then negligible. Precise formulae are developed for the geoid height and the free air anomaly using a global gravity model, given as spherical harmonic coefficients. Although only linear in the anomalous potential, these formulae are otherwise exact for an ellipsoidal reference Earth—they involve closed analytical functions of the eccentricity (and the Earths spin rate), rather than a truncated power series in e². They are evaluated using EGM96 and give ellipsoidal corrections to the conventional free air anomaly ranging from –0.84 to +1.14 mGal, both extremes occurring in Tibet. The geoid error corresponding to these differences is dominated by longer wavelengths, so extrema occur elsewhere, rising to +766 mm south of India and falling to –594 mm over New Guinea. At short wavelengths, the difference between ellipsoidal corrections based only on EGM96 and those derived from detailed local gravity data for the North Sea geoid GEONZ97 has a standard deviation of only 3.3 mm. However, the long-wavelength components missed by the local computation reach 300 mm and have a significant slope. In Australia, for example, such a slope would amount to a 600-mm rise from Perth to Cairns.