重力异常与垂线偏差资料

重力异常与垂线偏差重力相关资料1.相关坐标系地球上任何一个质点都同时受到地心引力和由于地球自转产生的离心力的作用，两个力的合力称为重力。

离心力与引力之比约为1：300，所以重力中起主要作用的还是地心引力。

重力的作用线称为铅垂线，重力线方向就是铅垂线方向。

1.1 水准面与大地水准面当液体处于静止状态时，其表面必处处与重力方向正交，否则液体就要流动。

这个液体静止的表面就称为水准面。

水准面是一个客观存在的、处处与铅垂线正交的面。

通过不同高度的点，都有一个水准面，所以水准面有无穷多个。

为了使测量结果有一个共同的基准面，可以选择一个十分接近地球表面又能代表地球形状和大小的水准面作为共同标准。

设想海洋处于静止平衡状态，并将它延伸到大陆内部且保持处处与铅垂线正交的水准面，来表示地球的形状是最理想的，这个面称为大地水准面。

它是一个光滑的闭合曲面，又称为地球的物理表面。

由它包围的形状是地球的真实形体，称为大地体。

地球自然表面的起伏不平、地壳内部物质密度分布不均，使得引力方向产生不规则的变化。

因而引力方向除总的变化趋势外，还会出现局部变化，这就引起铅垂线方向发生不规则的变化。

由于大地水准面处处与铅垂线正交，所以它是一个略有起伏的不规则的表面。

图1 椭球面与大地水准面1.2 参考椭球面从整体上看，大地体接近于一个具有微小扁率的旋转椭球，与大地体吻合的最好的旋转椭球称为总地球椭球，也叫总椭球或平均椭球。

要确定总椭球，必须在整个地球表面上布设连成一体的天文大地网和进行全球性的重力测量。

为了大地测量工作的实际需要，各个国家和地区只有根据局部的天文、大地和重力测量资料，研究局部大地水准面的情况，确定一个于总椭球相近的椭球，以表示地球的大小，作为处理大地测量成果的依据。

这样的椭球只能较好的接近局部地区的大地水准面，不能反映整个大地体的情况，所以叫做参考椭球面。

1.3大地坐标系与天文坐标系表1 大地坐标系与天文坐标系的对比由于大地水准面起伏，导致同一点的法线和垂线不一致，两者之间的微小夹角称为垂线偏差；导致天地高和海拔高(正高)不一致，两者之间的差距称为大地水准面差距。

设垂线偏差在子午面上的分量(即南北分量)以ξ来表示，在卯酉面上的分量(即东西分量)以η来表示，N表示大地水准面差距。

考虑垂线偏差和大地水准面差距，大地坐标与天文坐标数学转换公式为()cos tan BL A N H H ξϕηλϕαηϕ=-=-=-=-正大(1)2.重力场地球重力场的确定就是通过求解某种形式的大地测量边值问题得到一种表达扰动位或其泛函的数学模型，包括解析表达模型和数值模型两种形式。

其中，解析表达模型利用斯托克司(Stokes)公式或莫洛金斯基(Molodensky)级数给出积分表达式，或利用球谐函数技术给出谱展开式；数值模型则是重力场参数的一定分辨率的格网数值，包括扰动位、大地水准面、重力异常和垂线偏差等。

例如EGM 系列地球重力场模型计算软件给出的最终计算结果便是大地水准面差距(geoid undulations)，一般用N 表示。

研究地球重力场的理论基础即是研究大地测量边值问题的解算，这一边值问题的描述为：在大地水准面或地球自然表面给定边值条件和相应的边值，确定该边界面及其外部的引力位，并满足边值条件，同时在无限空间内是调和函数。

位理论自1785年由法国数学家勒让德(Legendre)提出后，经过格林(Green)和高斯(Gauss)等数学家和大地测量学家的进一步研究，不仅解释了表征位场的基本数学关系，同时将引力或重力的3个分力通过引入一个位函数在3个分力方向上的偏导数来表达。

同时，引入等位面和大地水准面的概念，从而将地球形状和重力场的研究统一起来，即地球重力场可用地球重力位来表达，地球形状可用大地水准面的形状来代表。

由于地球表面形状的不规则和内部质量分布不均匀，通常将地球重力位分为两部分：正常场和扰动场，分别对应于正常位和扰动位。

前者可用4个大地测量基本参数确定，是接近真地球位场的一种理想化位场；后者是对应的异常质量分布产生的异常位场。

因此，研究地球形状及外部重力场的关键在于确定扰动位。

2.1重力位函数设有一个标量函数，它对各坐标轴的偏导数等于力在相应坐标轴上的分量，此函数定义为位函数。

由此定义可得重力位函数为：W V =+Φ(2)其中V 为引力产生的位，Φ为地球自转产生的位。

若ω为角速度，X 和Y 为给定点的地心坐标，则221()2X Y ωΦ=+(3)引力位函数的定义为max 20(,,)1()(cos sin )(cos )n n n nm nm nm n m GMa V r C m S m P rr θλλλθ==⎡⎤=++⎢⎥⎣⎦∑∑ (4)其中，r — 到地球质心的距离θ — 地心纬度λ — 地心经度a — 参考椭球长半轴n ，m — 阶次nm C ，nm S — 正常化的地球引力系数(cos )nm P θ — 为n 阶m 次第一类完全正常化缔合Legendre 函数1/20()!(21)(2)(cos )(cos )()!m nm nm n m n P P n m δθθ⎡⎤-+-=⎢⎥+⎣⎦，1, 0, i j i ji jδ=⎧=⎨≠⎩(5)i j δ为克罗内克(kroneker)符号，(cos )nm P θ为蒂合Legendre 函数2/22/221()(1)()(1)(1)2!m n mm m nnm n m n n md d P x x P x x x dx n dx++=-=-- (6)这些公式对于r a ≥，即在地球表面或近地空间理论上是有效的，公式计算误差比较小，但这些公式不能用于r 小于地球半径的情况。

2.2重力加速度和重力梯度在地球重力场中，重力加速度为()TTxyz W W W g grad W g g g xyz ⎡⎤∂∂∂⎡⎤===⎢⎥⎣⎦∂∂∂⎣⎦r(7)重力加速度矢量的空间梯度即是重力位的二阶空间导数，称为重力梯度，可由下述张量表示：222222222222x x x xx xy xz y y y yxyyyz zx zyzz z z z g g g WWW x y z x x y x z g g g WW W x y z y x y y z g g g WW W xyz z xz yz ⎡⎤∂∂∂∂∂∂⎡⎤⎢⎥⎢⎥∂∂∂∂∂∂∂∂⎢⎥⎢⎥⎡⎤ΓΓΓ⎢⎥∂∂∂⎢⎥∂∂∂⎢⎥Γ===ΓΓΓ⎢⎥⎢⎥⎢⎥∂∂∂∂∂∂∂∂⎢⎥⎢⎥⎢⎥ΓΓΓ⎣⎦⎢⎥⎢⎥∂∂∂∂∂∂⎢⎥⎢⎥∂∂∂∂∂∂∂∂⎢⎥⎣⎦⎣⎦(8)在不考虑地球离心力影响的情况下，地球重力场是一个保守场，重力场的旋度为零，说明重力梯度张量具有对称性；重力场的散度为零，意味着梯度张量的迹为零。

因此,,0xy yx xz zx yz zy xx yy zz Γ=ΓΓ=ΓΓ=ΓΓ+Γ+Γ= (9)重力梯度张量中只有五个独立项。

3. WGS84世界大地坐标系与EGM20083.1 WGS84坐标系概述及有关常数WGS84世界大地坐标系是由美国国家影像制图局NIMA 和其前身美国国防部测绘局DMA 从初始的世界大地坐标系WGS60开始发展，并在随后的WGS66、WGS72基础上不断改进形成的。

WGS84坐标系使用协议地面参考系(CTRS)。

这个坐标系根据国际地球自转服务技术注解概要中的标准来定义。

这些标准如下：(1)它是一个地心坐标系，地心定义为包括海洋和大气的整个地球的质量中心。

(2)范围局限于地球，即在引力相对论意义下的局部地球框架内。

(3)它的定向由国际时间局(BIH1984.0)给出。

(4)定向的时间演化相对于地壳不产生残余的全球性旋转。

WGS84坐标系是右手地固直角坐标系，坐标原点和坐标轴定义如下： 原点：地球质心。

Z 轴：指向国际地球参考系IERS 极的方向(IRP)。

这个方向与BIH 协议地面极的指向(在历元1984.0)相差±0.005″。

X轴：指向国际地球参考系IERS首子午面且垂直于Z轴的方向，这个方向与BIH协议零子午面(在历元1984.0)相差±0.005″。

Y轴与其他两轴构成右手地心地固ECEF(Earth-Centered Earth-Fixed)直角坐标系。

表2 WGS84椭球基本常数表3 WGS84椭球导出的几何常数表4 WGS84椭球导出的物理常数3.2 WGS84水准椭球正常重力公式Simigliana闭合公式给出了椭球面理论正常重力γ，即正常位函数U的梯度大小公式：2eγγ=(10)式中1p eb k a γγ=-，a 、b 表示椭球长半轴和短半轴，e γ、p γ表示赤道和两极理论重力，e 表示第一离心率，φ表示大地纬度。

这个正常重力公式就是WGS84椭球重力公式。

这个等位椭球面不仅作为地球水准面和铅垂面或地球的几何参考面，而且也作为地球正常重力的参考面。

若以椭球面法线向下为正方向，一个经常被使用的正常重力泰勒展开式为222231(12sin )h f m f h h a a γγφ⎡⎤=-++-⋅+⎢⎥⎣⎦(11)当计算点的大地高比较高时，采用上式计算得到的结果精度要小于预期的精度。

解决办法详见《大地坐标系与大地基准》P116-117。

3.3 EGM2008模型官方版地球重力场模型EGM2008由美国国家地理空间情报局地球重力场研发小组发布，该模型的球谐展开阶数达到了2159，并且还提供扩展到2190阶的扩展参数，此模型目的是希望在WGS84坐标系下解算地球水准面的起伏。

EGM2008地球重力场模型数据来源主要为地面重力、卫星测高、卫星重力等，地面数据覆盖率达83.8%，部分重力数据空白区主要集中在南极，用卫星重力数据补充。

目前能够获取到的文件包括：引力位函数的球谐系数文件EGM2008_to2190_TideFree.gz扩展到2160阶的大地水准面差距计算参数文件Zeta-to-N_to2160_egm2008.gz1′ x 1′的大地水准面网格文件Und_min1x1_egm2008_isw=82_WGS84_TideFree.gz 2.5′ x2.5′的大地水准面网格文件Und_min2.5x2.5_egm2008_isw=82_WGS84_TideFree.gz 2.5′ x 2.5′的垂线偏差及重力异常网格文件：(1) 重力异常g ∆文件：Dg01_cnt2.5x2.5_EGM08_to2190_WGS84_ell_nh (2) 垂线偏差南北分量ξ文件：xi_cnt2.5x2.5_EGM08_to2190_WGS84_ell_nh (3) 垂线偏差东西分量η文件：eta_cnt2.5x2.5_EGM08_to2190_WGS84_ell_nh 5′ x 5′的垂线偏差及重力异常网格文件：(1) 重力异常g ∆文件：Dg01_cnt5x5_EGM08_to2190_WGS84_ell_nh (2) 垂线偏差南北分量ξ文件：xi_cnt5x5_EGM08_to2190_WGS84_ell_nh (3) 垂线偏差东西分量η文件：eta_cnt5x5_EGM08_to2190_WGS84_ell_nh3.3.1 EGM2008文件说明所有以下数据均为无潮汐系统下的数据。