顾及测站点上重力场信息的大地水准面高的拟合方法

利用测站点上的重力测量信息,根据Bruns公式研究了大地水准面高的变化特性,得出垂线偏差是大地水准面相对水准椭球面倾斜的线性改正,重力异常是大地水准面相对水准椭球面弯曲的线性改正,以及地形起伏效应构成大地水准面高的二阶变化等结论。在此基础上,提出了顾及测站点上重力场信息的大地水准面高的拟合方法,并分析了该方法相对于二次曲面函数拟合的优越性。     

整体大地测量的确定性参数模型

本文将整体大地测量各类观测方程中作为随机信号的重力场参数（ΔN，ξ，η）均表达为重力异常Δg的泛函，又引入Δg的B样条函数并以样条系数C_(11)作待定参数，提出一种整体大地测量的确定性参数模型，给出了相应系数矩阵的表达式和算法，最后对该模型的求解问题作了说明。     

不同扰动位泛函间的积分变换广义核函数

基于物理大地测量边值问题的解,利用一阶边界算子定义,推导重力异常Δg、单层密度μ、大地水准面高N,垂线偏差ε、扰动重力δg等扰动场元的解。利用球谐函数的正交特性,通过对核函数的算子运算,可以得到上述扰动场元的有关逆变换公式。相对经典物理大地测量公式应用的边界面条件,笔者将含有因子r的对应扰动场元反演关系的公式称为广义积分公式。针对常用的重力异常Δg、大地水准面高N,垂线偏差ε、扰动重力δg计算,重点分析它们之间的变换关系,给出利用某个选定扰动场元计算其他扰动场元的广义积分公式。同时,通过对积分边界面的讨论,分析经典公式与广义积分公式的差异和联系。最后,给出所有外部扰动场元与核函数映射的关系表。     

地球重力场参数计算程序的设计与应用

根据地球重力场参数的计算模型,借助全球地球重力场模型EGM2008,利用Microsoft Visual Studio 2010面向对象的功能,设计和开发了可以计算任意点重力异常、高程异常、垂线偏差分量的计算程序,以德国地学中心(GFZ)研制的相关软件对程序的正确性进行了验证,以吉林省各城市的经纬度数据计算了对应的重力异常、高程异常、垂线偏差分量。     

An evaluation of some systematic error sources affecting terrestrial gravity anomalies

Terrestrial free-air gravity anomalies form a most essential data source in the framework of gravity field determination. Gravity anomalies depend on the datums of the gravity, vertical, and horizontal networks as well as on the definition of a normal gravity field; thus gravity anomaly data are affected in a systematic way by inconsistencies of the local datums with respect to a global datum, by the use of a simplified free-air reduction procedure and of different kinds of height system. These systematic errors in free-air gravity anomaly data cause systematic effects in gravity field related quantities like e.g. absolute and relative geoidal heights or height anomalies calculated from gravity anomaly data. In detail it is shown that the effects of horizontal datum inconsistencies have been underestimated in the past. The corresponding systematic errors in gravity anomalies are maximum in mid-latitudes and can be as large as the errors induced by gravity and vertical datum and height system inconsistencies. As an example the situation in Australia is evaluated in more detail: The deviations between the national Australian horizontal datum and a global datum produce a systematic error in the free-air gravity anomalies of about −0.10 mgal which value is nearly constant over the continent     

基于逆Stokes公式的测高重力反演中央区效应计算

为提高利用逆Stokes公式反演测高重力的精度,将中央区大地水准面高表示成双三次多项式插值形式,引入了非奇异变换,推导出了重力异常的计算公式。大地水准面高理论模型下的分析表明,该公式有较高的精度。以分辨率为2′×2′的大地水准面高数据为背景场进行了实际计算,结果说明中央区对反演重力异常有不容忽视的贡献。本文导出的公式可为高精度重力异常的反演提供理论依据。     

Inverse Vening Meinesz formula and deflection-geoid formula: applications to the predictions of gravity and geoid over the South China Sea

Using the spherical harmonic representations of the earth's disturbing potential and its functionals, we derive the inverse Vening Meinesz formula, which converts deflection of the vertical to gravity anomaly using the gradient of the H function. The deflection-geoid formula is also derived that converts deflection to geoidal undulation using the gradient of the C function. The two formulae are implemented by the 1D FFT and the 2D FFT methods. The innermost zone effect is derived. The inverse Vening Meinesz formula is employed to compute gravity anomalies and geoidal undulations over the South China Sea using deflections from Seasat, Geosat, ERS-1 and TOPEX//POSEIDON satellite altimetry. The 1D FFT yields the best result of 9.9-mgal rms difference with the shipborne gravity anomalies. Using the simulated deflections from EGM96, the deflection-geoid formula yields a 4-cm rms difference with the EGM96-generated geoid. The predicted gravity anomalies and geoidal undulations can be used to study the tectonic structure and the ocean circulations of the South China Sea. Received: 7 April 1997 / Accepted: 7 January 1998     

基于球谐函数的重力异常和垂线偏差误差匹配关系

重力异常和垂线偏差是测高卫星非常重要的产品。二者的精度指标对于未来的测高卫星方案设计至关重要。本文利用球谐函数来对重力异常和垂线偏差的精度指标进行讨论,首先从理论上推导了重力异常和垂线偏差误差的近似匹配关系,然后通过6个超高阶重力场模型验证了有关结论的正确性。数值试验表明:垂线偏差误差和重力异常误差满足近似的比例关系,即若垂线偏差各方位向等精度测量,且假定精度均为1μrad,则所对应的重力异常精度约为1.4mGal;反之,若重力异常的精度为1mGal,则所对应的垂线偏差的精度约为0.7μrad。     

Relation between geoidal undulation, deflection of the vertical and vertical gravity gradient revisited

The vertical gradients of gravity anomaly and gravity disturbance can be related to horizontal first derivatives of deflection of the vertical or second derivatives of geoidal undulations. These are simplified relations of which different variations have found application in satellite altimetry with the implicit assumption that the neglected terms—using remove-restore—are sufficiently small. In this paper, the different simplified relations are rigorously connected and the neglected terms are made explicit. The main neglected terms are a curvilinear term that accounts for the difference between second derivatives in a Cartesian system and on a spherical surface, and a small circle term that stems from the difference between second derivatives on a great and small circle. The neglected terms were compared with the dynamic ocean topography (DOT) and the requirements on the GOCE gravity gradients. In addition, the signal root-mean-square (RMS) of the neglected terms and vertical gravity gradient were compared, and the effect of a remove-restore procedure was studied. These analyses show that both neglected terms have the same order of magnitude as the DOT gradient signal and may be above the GOCE requirements, and should be accounted for when combining altimetry derived and GOCE measured gradients. The signal RMS of both neglected terms is in general small when compared with the signal RMS of the vertical gravity gradient, but they may introduce gradient errors above the spherical approximation error. Remove-restore with gravity field models reduces the errors in the vertical gravity gradient, but it appears that errors above the spherical approximation error cannot be avoided at individual locations. When computing the vertical gradient of gravity anomaly from satellite altimeter data using deflections of the vertical, the small circle term is readily available and can be included. The direct computation of the vertical gradient of gravity disturbance from satellite altimeter data is more difficult than the computation of the vertical gradient of gravity anomaly because in the former case the curvilinear term is needed, which is not readily available.     

EGM2008地球重力场模型的可视化探究

利用EGM2008地球重力场模型可以将地球重力场解析化,由此我们可以十分方便快速地表示和计算重力异常、高程异常、垂线偏差等其它重力场参数。本文运用Surfer软件将EGM2008重力场模型可视化,获取了全球和局部地区重力场相关参数的等值线图、线框图、3D表面图、影像图。