1、1 第三章第三章 地球重力场及形状的基本理论地球重力场及形状的基本理论2 地球重力场状基本理论地球重力场状基本理论3.1.1 地球的概说（略）地球的概说（略）3.1.2 地球运动概说地球运动概说 地球是太阳系中的一颗行星，它有自转和公转运动。地球是太阳系中的一颗行星，它有自转和公转运动。 1、地球的自转、地球的自转 地球的自转地球的自转即地球绕地轴由西向东旋转。即地球绕地轴由西向东旋转。 地球的地球的绕地轴旋转绕地轴旋转360度的时间：太阳日、恒星日。度的时间：太阳日、恒星日。 地球的自转速度：地球的自转速度： 2cos)RhVT （2T 3地球重力场状基本理论地球重力场状基本理论2、地球的公
2、转、地球的公转 地球的公转满足开普勒三大行星运动定律地球的公转满足开普勒三大行星运动定律 (1) 行星运动轨迹是椭圆，太阳位于其行星运动轨迹是椭圆，太阳位于其 椭圆的椭圆的一个焦点上一个焦点上 直角坐标方程：直角坐标方程： 极坐标方程：极坐标方程： f 真近点角，真近点角，p 为焦参数（半通径）为焦参数（半通径） 22abea 22(1)bpaea 22221xyab 1cospref 4地球重力场状基本理论地球重力场状基本理论(2)行星运动在单位时间内扫过的面积相等；行星运动在单位时间内扫过的面积相等； 在时间在时间 t 内扫过的面积内扫过的面积 s 相等，则面速度相等，则面速度 可根据能量
3、守恒定律导出。可根据能量守恒定律导出。 (3) 行星运动的周期的平方与轨道的长半轴的立方的比为行星运动的周期的平方与轨道的长半轴的立方的比为常数。常数。 设设a 和和a1 , T 和和 T1分别表示两行星轨道的长半径与轨道分别表示两行星轨道的长半径与轨道运行周期。运行周期。221sabaetTT ABCDEFABCDEFVVV5地球重力场状基本理论地球重力场状基本理论则第三定律表达为：则第三定律表达为：一般可以用来计算行星或卫星的质量。一般可以用来计算行星或卫星的质量。牛顿万有引力定律：牛顿万有引力定律： 开普勒定律是牛顿万有引力定律的基础。开普勒定律是牛顿万有引力定律的基础。 天体力学天体力
4、学221331TTaa 322()4afMmT 2323111)TMmaTMma 6地球重力场状基本理论地球重力场状基本理论222M mM mFkfrr 22FMakmr 22222()()MmMmakkrrr 22224, vravrarTT 322() 4af M mT 宇宙空间任意两质点，彼此相互吸引，其引力大小与宇宙空间任意两质点，彼此相互吸引，其引力大小与他们的质量成积成正比，与他们之间的距离平方成反他们的质量成积成正比，与他们之间的距离平方成反比。比。在相对运动中，行星相对于太阳运动的相对加速度：在相对运动中，行星相对于太阳运动的相对加速度：7地球重力场状基本理论地球重力场状基本理
5、论考虑到Mm 注意： f 、 G、 k2 在不同的教材都表示引力常数。32 fMnTa 3223()2 n=4Taf MmfMTa 82rmMfF2mP PFg地球重力场的基本原理地球重力场的基本原理3.2.1 引力与离心力引力与离心力其它作用力（太阳、月亮）大多数情况下可忽略。其它作用力（太阳、月亮）大多数情况下可忽略。9地球重力场的基本原理地球重力场的基本原理3.2.2 引力位和离心力位引力位和离心力位 由理论力学可知，如果某一空间（有限或无限）的由理论力学可知，如果某一空间（有限或无限）的任意一点都有一定力的作用，而力的大小与方向只与该任意一点都有一定力的作用，而力的大小与方向只与该点的
6、位置有关，则这一空间称为力场。就力场而言，具点的位置有关，则这一空间称为力场。就力场而言，具有共同的特性，即力场所做的功与路径无关，只与起点有共同的特性，即力场所做的功与路径无关，只与起点与终点有关。这样的力称为保守力。引力与离心力都是与终点有关。这样的力称为保守力。引力与离心力都是保守力。保守力。 引力位：引力位：单位质点受物质单位质点受物质M的引力作用产生的位能称为的引力作用产生的位能称为引力位，或者说引力位，或者说将单位质点从无穷远处移动到该点引力将单位质点从无穷远处移动到该点引力所做的功。即所做的功。即： rMfVdrdVa10地球重力场的基本原理地球重力场的基本原理dA2rMmfFd
7、rrMmfdA2drrMmfdAdV2CrMmfVrMfVrMmfV万有引力定律：万有引力定律：推导如下推导如下:假设沿力线方向做功为假设沿力线方向做功为，则有，则有此功等于位能的减少，此功等于位能的减少，积分则有：积分则有：因为因为r, V=0。所以。所以 C=0 ，则有，则有取取 m=1，11地球重力场的基本原理地球重力场的基本原理rdmfdVVM )(2 rMmfFmaF2rMfa2rMfdrdVdrdVa地球总体的位函数：地球总体的位函数：1、由牛顿第二定律可知：、由牛顿第二定律可知：2、对位函数求导：、对位函数求导：, 则有则有12地球重力场的基本原理地球重力场的基本原理 结论：结论
8、： 单位质点的物体在引力场中的加速度等于引力位的单位质点的物体在引力场中的加速度等于引力位的导数，方向与径向方向相反。导数，方向与径向方向相反。 推论：推论： 位对被吸引点各坐标轴的偏导数等于相应坐标轴上位对被吸引点各坐标轴的偏导数等于相应坐标轴上的加速度的加速度( (或引力或引力) )向量的负值。向量的负值。 zVayVaxVazyx,13地球重力场的基本原地球重力场的基本原理理 离心力位 在离心力场中， PdldQ222222Q2ldlldldQ222222sin2)(2ryxQ14QVW)(2222yxrdmfW)()()(zQzVzWgyQyVyWgxQxVxWgzyx地球重力场的基本
9、原理地球重力场的基本原理3.2.3 重力位重力位重力是引力和离心力的合力，重力位重力是引力和离心力的合力，重力位W是引力位是引力位V和离和离心力位心力位Q之和：之和：对三坐标轴求偏导数求得重力的分力或重力加速度分量对三坐标轴求偏导数求得重力的分力或重力加速度分量：15各分力的模：各分力的模：方向余弦：方向余弦： 重力位在任意方向的偏导数等于重力在该方向上的重力位在任意方向的偏导数等于重力在该方向上的分力：分力： 222zyxggggggzgggygggxgzyx),cos(,),cos(,),cos(),cos(lggglWl地球重力场的基本原理地球重力场的基本原理16gdWdl地球重力场的基
10、本原理地球重力场的基本原理u当当g与与l相垂直时，那么相垂直时，那么d=0，常数，常数当给出不同的常数值，就得到一簇曲面，称为重力等当给出不同的常数值，就得到一簇曲面，称为重力等位面，也就是我们通常说的水准面。可见水准面有无位面，也就是我们通常说的水准面。可见水准面有无穷多个。其中，我们把完全静止的海水面所形成的重穷多个。其中，我们把完全静止的海水面所形成的重力等位面，专称它为力等位面，专称它为大地水准面大地水准面。u如果令如果令g与与l夹角等于夹角等于，则有：则有：u水准面之间既不平行，也不相交和相切。水准面之间既不平行，也不相交和相切。17 对于某一单位质点而言，作用其上的重力在数值上等于
11、对于某一单位质点而言，作用其上的重力在数值上等于使它产生的重力加速度的数值，所以重力即采用重力加速度使它产生的重力加速度的数值，所以重力即采用重力加速度的量纲，单位是：的量纲，单位是： 伽伽(Gal=cms)， 毫伽毫伽(mGal= Gal/1000=10ms) 微伽微伽(Gal= mGal/1000=10m s) 1 1、地面点重力近似值、地面点重力近似值 980Gal，赤道重力值赤道重力值 978Gal，两两极重力值极重力值 983Gal。由于地球的极曲率及周日运动的原因，重由于地球的极曲率及周日运动的原因，重力有从赤道向两极增大的趋势。力有从赤道向两极增大的趋势。 2 2、地球上重力的大
12、小与方向只与被吸引点的位置有关，、地球上重力的大小与方向只与被吸引点的位置有关，理论上应该是常数，但重力是随时间变化而变化，即相同的理论上应该是常数，但重力是随时间变化而变化，即相同的点在不同的时刻所观测到的重力不相同点在不同的时刻所观测到的重力不相同。地球重力场的基本原理地球重力场的基本原理183.2.4 地球的正常重力位和正常重力地球的正常重力位和正常重力 要精确计算出地球重力位，必须知道地球表面的形要精确计算出地球重力位，必须知道地球表面的形状及内部物质密度，但前者正是我们要研究的，后者分状及内部物质密度，但前者正是我们要研究的，后者分布极其不规则，目前也无法知道，故根据上式不能精确布极
13、其不规则，目前也无法知道，故根据上式不能精确地求得地球的重力位，为此引进一个与其近似的地球重地求得地球的重力位，为此引进一个与其近似的地球重力位力位正常重力位正常重力位。 )(2222yxrdmfWM地球重力场的基本原理地球重力场的基本原理19地球重力场的基本原理地球重力场的基本原理( , , )x y zMdmfV( , , )r (,)mmmxyz(, )mmR 正常重力位是一个函数简单、不涉及地球形状和密正常重力位是一个函数简单、不涉及地球形状和密度便可直接计算得到的地球重力位的近似值的辅助重力度便可直接计算得到的地球重力位的近似值的辅助重力位。当知道了地球正常重力位，想法求出它同地球重
14、力位。当知道了地球正常重力位，想法求出它同地球重力位的差异位的差异(称扰动位称扰动位)，便可求出大地水准面与这已知形状，便可求出大地水准面与这已知形状(正常位水准面正常位水准面)的差异。最后解决确定地球重力位和地球的差异。最后解决确定地球重力位和地球形状的问题。形状的问题。1 地球引力位的数学表达式地球引力位的数学表达式 地球惯性矩表达引力位地球惯性矩表达引力位 (方法方法1)设地球上的点坐标为设地球上的点坐标为:地球表面点坐标为地球表面点坐标为:与与与与20 建立空间直角坐标系与球面极坐标系建立空间直角坐标系与球面极坐标系cos2)(1 cos222222rRrRrRrRrcos2)(2rR
15、rRl21)1 (11lrdmlllrfV)16583211 (32niivvvvV0210地球重力场的基本原理地球重力场的基本原理21地球重力场的基本原地球重力场的基本原理理由于由于rMfdmrfvM00cos1MdmrRrfvMdmrRrfv)21cos23()(222MdmrRrfv)cos23cos25()(333MmmMmmMmmdmyxCdmzxBdmzyA)()()(222222MmmMmmMmmdmyxFdmzxEdmzyD)()()(rRzzyyxxmmmcos22地球重力场的基本原理地球重力场的基本原理理论力学可知：物体的重心为理论力学可知：物体的重心为定义坐标系：定义坐标
16、系： ，则有：，则有： 0000zyxByzxAxzyrfv)2()2(222222252666)2(222xyFxzEyzDCzyxMmMmMmdmMzdmMydmxMxz1, y1, 1000Mrfv00)(31dmzzdmyydmxxrfvMmMmMm0000zyx23 用球谐函数表达地球引力位（方法2） 勒让德多项式勒让德多项式nnnnndxxdnxP)1(!21)(2)(1)(112)(11xPnnxxPnnxPnnn)()(01xxPxP地球重力场的基本原理地球重力场的基本原理24地球重力场的基本原理地球重力场的基本原理cos23cos25)(cos21cos23)(coscos)
17、(cos1)(cos332210PPPPdmPrRrfVnnn)(cos)(25地球重力场的基本原理地球重力场的基本原理 勒让德多项式中：勒让德多项式中： 称为称为n阶主球函数阶主球函数(或带球或带球函数函数)， 称为称为n阶阶K级的勒让德缔合函数级的勒让德缔合函数(或伴随函或伴随函数数) )。 称为缔合球函数称为缔合球函数( (其中，当其中，当k=nk=n时称为扇球函数，当时称为扇球函数，当knkn时称时称为田球函数为田球函数) ) 090 nKKnKnKnnnnnPKBKAPArV11)(cos)sincos()(cos1)(cosnP)(cosKnP)(cossin ),(coscosK
18、nKnPKPK(cos )(cos )sin(cos )kkknnkd pPd26地球重力场的基本原理地球重力场的基本原理 用球谐函数表示的地球引力位的公式用球谐函数表示的地球引力位的公式 2 地球正常重力位 )(cos1010nnnnnnPArVV)(cos)sincos(1KnKnKnnKPKBKA222sin2rVW00kn10)(cos)sincos(1nnkknknnnnPkBkArVV27地球重力场的基本原理地球重力场的基本原理 当选取前当选取前3项时，将重力位项时，将重力位W W写成写成U MnnndmPRfA)cos(00)sincos()(cos121201KKnKnnnnn
19、KBKAPArU222sin2)(cosrPKn()!2(cos)sindm ,1,()!knknnmmMnkBfR Pkknnk00(cos)nnnmMAfR Pdm()!2(cos)cosdm()!knknnmmMnkAfR Pknk28地球重力场的基本原理地球重力场的基本原理现在需要求系数：现在需要求系数：若地球是旋转椭球体，则有转动惯量若地球是旋转椭球体，则有转动惯量 ，将系数代入，将系数代入则有：则有：式中：式中： 001101122011122222, , , , ,AAABAABABsin2)cos31 (21 23222fMrrKrMfUKMA C00AfM22= () 4BA
20、Af0111110AAB01122222() , =02ABAfCABBAB29地球重力场的基本原理地球重力场的基本原理设赤道的离心力与重力之比为：设赤道的离心力与重力之比为：令：令：则有：则有：22232eaaaqagfMfMsin2)cos31 (31 22qrMfU232Ka30地球重力场的基本原理地球重力场的基本原理 注意：如果正常重力位已知，则对应的正常水准面已知，注意：如果正常重力位已知，则对应的正常水准面已知，不同的正常重力位对应不同的正常位水准面，我们寻找不同的正常重力位对应不同的正常位水准面，我们寻找的是与大地水准面相近的正常位水准面的形状，的是与大地水准面相近的正常位水准面
21、的形状，上式中，上式中，对对r r和和 取不同的常数值，就得到一簇正常位水准面，取不同的常数值，就得到一簇正常位水准面，取取 ，求得与大地水准面相近的正常位水，求得与大地水准面相近的正常位水准面方程：准面方程： 取：取： ，则有，则有ar ,90ar ,900132MqUfa2201(13cos)sin32MqUfUr221(1 3cos)sin/(1)3232qqra31地球重力场的基本原理地球重力场的基本原理 另外，旋转椭球面的方程：另外，旋转椭球面的方程： 则有：则有： 4.4.3 4.4.3正常重力公式正常重力公式 因为：因为： ()2qcos)2(1 2qardrdU2(1cos)r
22、a222(1(13cos)sin)fMqrcos)2(1 2qar32地球重力场的基本原地球重力场的基本原理理 特例：特例： ，赤道正常重力：赤道正常重力： ，极点处正常重力：，极点处正常重力： 令：令： 则有：则有： 上述正常重力公式称为上述正常重力公式称为克莱罗定理。克莱罗定理。9020235(1()cos)22fMqqa)231 (2qafMe0)1 (2qafMp55 , +=22peeqq20(1sin)e33地球重力场的基本原理地球重力场的基本原理 顾及到扁率的二次项的正常重力公式51 7(1)23 5q2111()84 2201(1sinsin2)e式中：34 19011909年
23、赫尔默特公式：年赫尔默特公式： 19301930年卡西尼公式：年卡西尼公式： 19751975年国际地球正常重力公式：年国际地球正常重力公式： GS84GS84坐标系中的椭球重力公式：坐标系中的椭球重力公式： )2sin000007. 0sin005302. 01 (030.978220)2sin000059. 0sin0052882. 01 (049.978220)2sin0000058. 0sin005302. 01 (032.978220/ )sin86390019318513. 01 (03267714.9782B122(1 0.00669437999013sin)B地球重力场的基本原
24、理地球重力场的基本原理35 高出水准椭球面高出水准椭球面H米的正常重力计算公式米的正常重力计算公式2)(HRMfg20RMfg)(11(2201HRRfMggg)1 (11 (22RHRfM2200220132)321 (1 RHRHRHRHg2711072. 03086. 0HHgH3086. 00地球重力场的基本原理地球重力场的基本原理36 4 正常重力场参数 在物理大地测量中在物理大地测量中, ,正常椭球重力场可用正常椭球重力场可用4 4个基本参数决定个基本参数决定, ,即：即： 地球正常地球正常(水准水准)椭球的基本参数，又称椭球的基本参数，又称地球大地基准常数地球大地基准常数是：是：
25、 其中：其中： 002, , (), UAfMAf ACfKM ,2fMJa222JfMaA2223322222AqqqJa fM fMaq32)231 (0qafMU223aK地球重力场的基本原理地球重力场的基本原理373.2.5 正常椭球、水准椭球、总地球椭球与参考椭球正常椭球、水准椭球、总地球椭球与参考椭球 正常椭球面正常椭球面 是大地水准面的规则形状（一般指旋转椭球是大地水准面的规则形状（一般指旋转椭球面）。因此引入正常椭球后，地球重力位被分成正常重力面）。因此引入正常椭球后，地球重力位被分成正常重力位和扰动位两部分，实际重力也被分成正常重力和重力异位和扰动位两部分，实际重力也被分成正
26、常重力和重力异常两部分。常两部分。 正常椭球的确定：正常椭球的确定： 1 1、除了确定其、除了确定其M M和和值外，其规则形状可以任意选择。值外，其规则形状可以任意选择。但考虑到实际使用的方便，又顾及几何大地测量中采用旋但考虑到实际使用的方便，又顾及几何大地测量中采用旋转椭球的实际情况，目前都采用水准椭球作为正常椭球。转椭球的实际情况，目前都采用水准椭球作为正常椭球。 2 2、对于正常椭球，除了确定其对于正常椭球，除了确定其4个基本参数：个基本参数：a, J，fM和和外，也要定位和定向。正常椭球的定位是使其中心和外，也要定位和定向。正常椭球的定位是使其中心和地球质心重合，正常椭球的定向是使其短
27、轴与地轴重合，地球质心重合，正常椭球的定向是使其短轴与地轴重合，起始子午面与起始天文子午面重合。起始子午面与起始天文子午面重合。 地球重力场的基本原理地球重力场的基本原理38地球重力场的基本原理地球重力场的基本原理 总的地球椭球：总的地球椭球： 一个和整个大地体最为密合的。总地球椭球中心和地球一个和整个大地体最为密合的。总地球椭球中心和地球质心重合质心重合，总的地球椭球的短轴与地球地轴相重合，起始总的地球椭球的短轴与地球地轴相重合，起始大地子午面和起始天文子午面重合，总地球椭球和大地体大地子午面和起始天文子午面重合，总地球椭球和大地体最为密合。最为密合。 从几何和物理两个方面来研究全球性问题，
28、我们可把从几何和物理两个方面来研究全球性问题，我们可把总地球椭球定义为最密合于大地体的正常椭球。正常椭球总地球椭球定义为最密合于大地体的正常椭球。正常椭球参数是根据天文大地测量，重力测量及人卫观测资料一起参数是根据天文大地测量，重力测量及人卫观测资料一起处理确定的，并由国际组织发布。处理确定的，并由国际组织发布。 参考椭球：参考椭球： 其大小及定位定向最接近于本国或本地区的地球椭球。其大小及定位定向最接近于本国或本地区的地球椭球。这种最接近，表现在两个面最接近及同点的法线和垂线最这种最接近，表现在两个面最接近及同点的法线和垂线最接近。接近。393.3.1一般说明一般说明 大地高由两部分组成：地
29、形高部分大地高由两部分组成：地形高部分(含含H正或正或H正常正常)及大地水准面及大地水准面(或似大地水准面或似大地水准面)高部分。地形高基本上高部分。地形高基本上确定着地球自然表面的地貌，大地水准面高度又称确定着地球自然表面的地貌，大地水准面高度又称大地大地水准面差距水准面差距 N，似大地水准面高度又称，似大地水准面高度又称高程异常高程异常，它，它们基本上确定着大地水准面或似大地水准面的起伏。因们基本上确定着大地水准面或似大地水准面的起伏。因此，大地高可表示为：此，大地高可表示为： HHN大正高H+H大正 常3.3 3.3 高高 程系程系 统统40设由设由OAB路线水准测量得到路线水准测量得到
30、B点的高程点的高程由由ONB线路得到线路得到B点高程点高程由于水准面不平行，对应的由于水准面不平行，对应的和和不相等，水准不相等，水准环线高程闭合差也不等于零，称为环线高程闭合差也不等于零，称为 理论闭合差理论闭合差。 hHBhHB高高 程系程系 统统413.3.2 正高系统正高系统 正高系统正高系统是以大地水准面为高程基准面，地面上任一是以大地水准面为高程基准面，地面上任一点的正高是该点沿垂线方向至大地水准面的距离。点的正高是该点沿垂线方向至大地水准面的距离。 因为无限接近两水准面其位能差可以写为因为无限接近两水准面其位能差可以写为 CBCBBdHHH正dHggdhBdhggdHBOABBC
31、BBdhggdHH正OABBmBgdhgH1正高高 程系程系 统统423.3.3 正常高系统正常高系统 将正高系统中不能精确测定的将正高系统中不能精确测定的 用正常重力代替，便得用正常重力代替，便得到另一种系统的高程，称其为正常高。到另一种系统的高程，称其为正常高。我国规定采用正我国规定采用正常高高程系统作为我国高程的统一系统。常高高程系统作为我国高程的统一系统。 正常高高差的实际计算公式正常高高差的实际计算公式 BmggdhHBmB1常BABAABABdhHH常常mHm2sin0000015395.0()mmgH高高 程系程系 统统43 说明：说明： 1 1、正常高与正高不同，它不是地面点到
32、大地水准面的、正常高与正高不同，它不是地面点到大地水准面的距离，而是地面点到一个与大地水准面极为接近的基距离，而是地面点到一个与大地水准面极为接近的基准面的距离，这个基准面称为似大地水准面。因此，准面的距离，这个基准面称为似大地水准面。因此，似大地水准面是由地面沿垂线向下量取正常高所得的似大地水准面是由地面沿垂线向下量取正常高所得的点形成的连续曲面，它不是水准面，只是用以计算的点形成的连续曲面，它不是水准面，只是用以计算的辅助面。因此，我们可以把正常高定义为以似大地水辅助面。因此，我们可以把正常高定义为以似大地水准面为基准面的高程。准面为基准面的高程。 2 2、正常高和正高之差、正常高和正高之
33、差, ,在高山地区可达在高山地区可达4 4米，在平原地米，在平原地区数厘米，在海水面上相等，大地水准面的高程原点区数厘米，在海水面上相等，大地水准面的高程原点对似大地水准面也是适用的。对似大地水准面也是适用的。 高高 程系程系 统统443.3.4 力高和地区力高高程系统力高和地区力高高程系统 同一个重力位水准面上两点的正高或正常高是不相等同一个重力位水准面上两点的正高或正常高是不相等的。对于大型水库等工程项目，它的静止水面是一个重力的。对于大型水库等工程项目，它的静止水面是一个重力等位面，在设计、施工、放样等工作中，通常要求这个水等位面，在设计、施工、放样等工作中，通常要求这个水面是一个等高面
34、。这时若继续采用正常高或正高显然是不面是一个等高面。这时若继续采用正常高或正高显然是不合适的，为了解决这个矛盾，可以采用所谓力高系统，它合适的，为了解决这个矛盾，可以采用所谓力高系统，它按下式定义：按下式定义： AOAgdhH451力AOAgdhH1力高高 程系程系 统统45 注意：注意：说明力高是区域性的，主要用于大型水库等工说明力高是区域性的，主要用于大型水库等工程建设中。它不能作为国家统一高程系统。在工程测程建设中。它不能作为国家统一高程系统。在工程测量中，应根据测量范围大小，测量任务的性质和目的量中，应根据测量范围大小，测量任务的性质和目的等因素，合理地选择正常高，力高或区域力高作为工
35、等因素，合理地选择正常高，力高或区域力高作为工程的高程系统。程的高程系统。 高高 程系程系 统统463.3.5 国家高程基准国家高程基准1 1、高程基准面、高程基准面 高程基准面：高程基准面：就是地面点高程的统一起算面，由于就是地面点高程的统一起算面，由于大地水准面所形成的体形大地水准面所形成的体形大地体是与整个地球最大地体是与整个地球最为接近的体形，因此通常采用大地水准面作为高程基为接近的体形，因此通常采用大地水准面作为高程基准面。准面。 高程基准面的确定高程基准面的确定：在海洋近岸的一点处竖立水位在海洋近岸的一点处竖立水位标尺，成年累月地观测海水面的水位升降，根据长期标尺，成年累月地观测海
36、水面的水位升降，根据长期观测的结果可以求出该点处海洋水面的平均位置，假观测的结果可以求出该点处海洋水面的平均位置，假定大地水准面就是通过这点处实测的平均海水面。定大地水准面就是通过这点处实测的平均海水面。 验潮、验潮站验潮、验潮站 高高 程系程系 统统47 19561956年黄海高程系统年黄海高程系统：19501950年至年至19561956年年7 7年间青岛年间青岛验潮站的潮汐资料推求的平均海水面作为我国的高程验潮站的潮汐资料推求的平均海水面作为我国的高程基准面。基准面。 19851985国家高程基准国家高程基准：根据青岛验潮站根据青岛验潮站 1952195219791979年中取年中取19
37、19年的验潮资料计算确定，并从年的验潮资料计算确定，并从19881988年年1 1月月1 1日日开始启用。开始启用。 高高 程系程系 统统482、水准原点水准原点 为了长期、牢固地表示为了长期、牢固地表示出高程基准面的位置，作出高程基准面的位置，作为传递高程的起算点，必为传递高程的起算点，必须建立稳固的水准起算点，须建立稳固的水准起算点，用精密水准测量方法将它用精密水准测量方法将它与验潮站的水准标尺进行与验潮站的水准标尺进行联测，以高程基准面为零联测，以高程基准面为零推求水准原点的高程。推求水准原点的高程。高高 程系程系 统统49 19561956年黄海高程系统中，我国水准原点的高程为年黄海高
38、程系统中，我国水准原点的高程为72.289m 1985国家高程基准系统中，我国水准原点的高程为国家高程基准系统中，我国水准原点的高程为72.260m。 地面上的点相对于高程基准面的高度，通常称为绝对地面上的点相对于高程基准面的高度，通常称为绝对高程或海拔高程，也简称为标高或高程。海洋的深度高程或海拔高程，也简称为标高或高程。海洋的深度也是相对于高程基准面而言的，例如太平洋的平均深也是相对于高程基准面而言的，例如太平洋的平均深度为度为40004000m m，就是说在高程基准面以下就是说在高程基准面以下40004000m m。 mHH029. 05685高高 程系程系 统统503.4 测定垂线偏差
39、和大地水准面差距的基本概念测定垂线偏差和大地水准面差距的基本概念 大地坐标同天文坐标的区别主要是由同一点的法大地坐标同天文坐标的区别主要是由同一点的法线和垂线不一致，亦即由垂线偏差引起的。线和垂线不一致，亦即由垂线偏差引起的。 地面一点上的重力向量地面一点上的重力向量g g和相应椭球面上的法线向和相应椭球面上的法线向量量 n n之间的夹角定义为该点的之间的夹角定义为该点的垂线偏差垂线偏差。很显然，根。很显然，根据所采用的椭球不同可分为据所采用的椭球不同可分为绝对垂线偏差绝对垂线偏差及及相对垂线相对垂线偏差偏差，垂线同总地球椭球，垂线同总地球椭球( (或参考椭球或参考椭球) )法线构成的角法线构
40、成的角度称为绝对度称为绝对( (或相对或相对) )垂线偏差，它们统称为天文大地垂线偏差，它们统称为天文大地垂线偏差垂线偏差。 51 图中，图中，u是垂线偏差，是垂线偏差，、分别是分别是u在子午圈和卯酉在子午圈和卯酉圈上的分量圈上的分量 AAuAsincos222uBB)90(90cos)(LsecLB天ZAALAcot)cossin(sin)(sin)(LAAAZZsincos0垂线偏差垂线偏差52垂线偏垂线偏差差1 1、天文大地测量方法、天文大地测量方法 在天文大地点上，既进行大地测量取得大地在天文大地点上，既进行大地测量取得大地坐标坐标(B,L)，又进行天文测量取得天文坐标又进行天文测量取
41、得天文坐标(，）。）。 2 2、重力重力测量方法测量方法 建立扰动位与垂线偏差的关系，即扰动位与观建立扰动位与垂线偏差的关系，即扰动位与观测量测量( (重力异常重力异常) )的函数的函数000()TF g0011TT XY Bcos)(L53垂线偏差垂线偏差00gg020cos)(21dAAdgQ020sin)(21dAAdgQ2222cos3( )csc12sin32sin12sinln(sinsin)21 sinQ2维宁.曼尼兹公式 此公式是在假定大地水准面之外没有扰动物质及全此公式是在假定大地水准面之外没有扰动物质及全球重力异常球重力异常都已知的情况下推导的。然而这两个条都已知的情况下推
42、导的。然而这两个条件都还不能实现，所以重力方法至今也没有得到独立的件都还不能实现，所以重力方法至今也没有得到独立的应用。应用。543 3、天文重力方法、天文重力方法 综合利用天文大地方法和重力测量方法来确定垂线综合利用天文大地方法和重力测量方法来确定垂线偏差偏差 4、GPS测量方法测量方法 在在GPS相对定位中，只要测出基线长相对定位中，只要测出基线长D，大地方位角大地方位角A及高程异常差及高程异常差，便可求得垂线偏差。但这种方法应便可求得垂线偏差。但这种方法应用是有条件的，比如，地形平坦，基线不长，精度要求用是有条件的，比如，地形平坦，基线不长，精度要求较低等。较低等。 D),2, 1(si
43、ncosniAAiiiii垂线偏差垂线偏差55 3.4.2 测定大地水准面差距的基本概念测定大地水准面差距的基本概念 1、用地球重力场模型法计算大地水准面差距用地球重力场模型法计算大地水准面差距 大地水准面上一点大地水准面上一点P P的实际重力位的实际重力位 与相应于点与相应于点P P的正常重力位的正常重力位 U U 之差，称之为该点的扰动位之差，称之为该点的扰动位T T，用下式用下式表示表示 由于在选择正常重力位时总是使地球离心力位对和由于在选择正常重力位时总是使地球离心力位对和U U的影响相同，因此扰动位具有引力位的性质。的影响相同，因此扰动位具有引力位的性质。 UWT大地水准面差距大地水
44、准面差距56大地水准面差距大地水准面差距57CTUW00000000000000UWTUTWUUdUN000TTNnmmnmnmnnnPmSmCrarGMrT0,20)(cos)sincos()(),(,20( )(cossin)(cos )nnn mn mn mnmGMaNCmSmPrr大地水准面差距大地水准面差距58大地水准面差距大地水准面差距 dSgRT)()(4100000020)()(4dAdSgRN2( )csc()6sin1 5cos3cosln(sinsin)2222S 000TTN2 2、利用斯托克司积分公式计算、利用斯托克司积分公式计算已知：已知：则有则有：59大地水准面差
45、距大地水准面差距 3、卫星无线电测高方法研究大地水准面卫星无线电测高方法研究大地水准面0rrh604、利用利用GPS高程拟合法研究似大地水准面高程拟合法研究似大地水准面 如果在测区中选择一定的如果在测区中选择一定的GPSGPS点同时联测几何水准点同时联测几何水准测量，求出这些点的正常高测量，求出这些点的正常高h h，于是在这些点上便可求于是在这些点上便可求出高程异常出高程异常: :代入下面数学拟合方程中用最小二乘求解代入下面数学拟合方程中用最小二乘求解： hH yaxaa210 xyayaxaa321025243210yaxaxyayaxaa大地水准面差距大地水准面差距615、利用利用最小二乘
46、配置法最小二乘配置法研究大地水准面研究大地水准面0)( LXFeLL01nSlAX)(XFLl0nBSlAX0000000000011lXKSnAABIBCICTTSn)()()()()()()()(11111111111XAlCCCnXAlCCBCSXAlCCKlCCAACCAXSnnSnTSSnSnTSnT大地水准面差距大地水准面差距623.5 地球形状的基本概念地球形状的基本概念1、天文大地测量方法天文大地测量方法 弧线法弧线法面积法面积法 现代推求新的椭球元素是在原有旧的椭球元素基础现代推求新的椭球元素是在原有旧的椭球元素基础上，综合利用天文、大地、重力及空间测量等资料，同上，综合利用
47、天文、大地、重力及空间测量等资料，同椭球定向、定位等一起实现的。椭球定向、定位等一起实现的。 43222111,(),(BBafSBBafS地球形状地球形状63000sinsincoscoscos)(cos)(sinsin)(cossin0)(cos)(sinZYXBLBLBHMBHMLBHMLBHNLHNLN旧新新新mBeNBBeMNLBBNeLBBNeLLBLBLBzyx)sin1 (cossin00coscossinsincossin0cossincossinsincossin22222旧地球形状地球形状64地球形状地球形状旧NBBLcos)(aBBeMBeaNBBHMBeMBBeaHMN旧22222222sin)sin1 (1)sin1 (cossin)1)()sin2(cossin)(0065)sin1()sin1(222121BgBgeeq25eaq2地球形状地球形状2 2、重力测量方法、重力测量方法根据克莱罗定理：根据克莱罗定理：3 3、 空间大地测量方法空间大地测量方法
