第三章1重力场基本理论

第三章1.地球重力场的基本理论——引力、离心力与重力——引力位与离心力位——地球的正常重力位——正常重力公式难，只需掌握基本概念重力场理论为研究课程上一讲应掌握的内容1、测量坐标参考系统由基准和坐标系两方面要素构成。基准是指用以描述地球形状的参考椭球的参数以及参考椭球在空间中的定位及定向，甚至还包括单位长度。狭义的坐标系是指点位表示方法（3种），广义的坐标系是由坐标原点、坐标轴的指向和尺度所定义的。2、建立地固坐标系统必须解决的问题确定椭球的形状和大小（长半径a和扁率α等）；确定椭球中心的位置（椭球定位）；确定椭球短轴的指向（椭球定向）；建立大地原点。()上一讲应掌握的内容3、1954年北京坐标系的特点4、1980年国家大地坐标系特点5、新1954年北京坐标系的特点6、WGS-84世界大地坐标系7、站心坐标系以测站为原点，测站上的法线(垂线)为Z轴方向的坐标系就称为法线(或垂线)站心坐标系。常用来描述参照于测站点的相对空间位置关系。工程上在小范围内有时也直接采用站心坐标系。8、不同空间直角坐标系转换考查学生一、牛顿万有引力定律宇宙空间任意两质点，彼此相互吸引，其引力大小与他们的质量乘积成正比，与他们之间的距离平方成反比。在相对运动中，行星相对于太阳运动的相对加速度：向心加速度二、引力、离心力与重力用F及P分别表示地球引力及质点绕地球自转轴旋转而产生的离心力，这两个力的合力称地球重力用g表示离心力P在赤道达最大值，但数值比地球引力1/200还要小一些。故重力基本上由地球引力确定的。当高出地面35730km处，重力加速度将改变符号，背向地球。三、引力位和离心力位位函数的概念：设有一标量函数，它对被吸引点各坐标方向的偏导数等于引力在相应方向上的分力，则此函数称为位函数。位函数的形式为：（一）引力位：单位质点受物质M的引力作用产生的位能称为引力位，或者说将单位质点从无穷远处移动到该点引力所做的功。即：由牛顿第二定律可导出：单位质点的物体在引力场中的加速度等于引力位的导数，方向与径向方向相反。（一）引力位万有引力定律：推导如下:假设沿力的方向做功为，则有此功等于位能的减少，积分则有：因为r→∞,V=0。所以C=0，则有取m=1，（二）离心力位质点坐标可用质点向径r，地心纬度φ及经度λ表示为：（图3-2）地球自转仅仅引起经度变化，而它对时间的一阶导数等于地球自转角速度ω时，得坐标对时间的二阶导数就是质点的离心加速度，即该方向分力。（二）离心力位(续)假定一个函数（离心力位）：则有：因此，我们可把Q称为离心力位函数。（三）重力位重力是引力和离心力的合力，重力位W是引力位V和离心力位Q之和：对三坐标轴求偏导数求得重力的分力或重力加速度分量：由各分力可计算重力加速度（模）：重力位在任意方向的偏导数等于重力在该方向上的分力：

重力位(续)当g与l相垂直时，那么dＷ=0，有Ｗ＝常数当给出不同的常数值，就得到一簇曲面，称为重力等位面，也就是我们通常说的水准面。可见水准面有无穷多个。其中，我们把完全静止的海水面所形成的重力等位面，专称它为大地水准面。如果令g与l夹角等于π，则有：水准面之间既不平行，也不相交和相切。重力的单位对于某一单位质点而言，作用其上的重力在数值上等于使它产生的重力加速度的数值，所以重力即采用重力加速度的量纲，单位是：伽(Gal=cms－２)，毫伽(mGal=Gal/1000=10－５ms－２)微伽(μGal=mGal/1000=10－８ms－２)地面点重力近似值980Gal，赤道重力值978Gal，两极重力值983Gal。由于地球的极曲率及周日运动的原因，重力有从赤道向两极增大的趋势。地球上重力的大小与方向只与被吸引点的位置有关，理论上应该是常数，但重力是随时间变化而变化，即相同的点在不同的时刻所观测到的重力不相同。(微)极赤重力差约5伽四、地球的正常重力位要精确计算出地球重力位，必须知道地球表面的形状及内部物质密度，但前者正是我们要研究的，后者分布极其不规则，目前也无法知道，故根据上式不能精确地求得地球的重力位，为此引进一个与其近似的地球重力位——正常重力位。正常重力位是一个函数简单、不涉及地球形状和密度便可直接计算得到的地球重力位的近似值的辅助重力位。当知道了地球正常重力位，想法求出它同地球重力位的差异(称扰动位)，便可求出大地水准面与这已知形状(正常位水准面)的差异。最后解决确定地球重力位和地球形状的问题。重力位用球谐函数表示的地球引力位的公式勒让德多项式称为n阶主球函数(或带球函数)。称为n阶K级的勒让德缔合函数(或伴随函数)。

称为缔合球函数(其中，当k=n时称为扇球函数，当k≠n时称为田球函数)。第n阶地球引力位公式地球正常重力位当选取前3项时，将重力位W写成U地球正常重力位的公式现在需要求系数：若地球是旋转椭球体，则有转动惯量，将系数代入则有：式中：

设赤道的离心力与重力之比为：令：地球形状参数。则有：地球正常重力位的公式与大地水准面相近的正常位水准面方程注意：如果正常重力位已知，则对应的正常水准面已知，不同的正常重力位对应不同的正常位水准面，我们寻找的是与大地水准面相近的正常位水准面的形状，上式中，对r和取不同的常数值，就得到一簇正常位水准面。取，求得与大地水准面相近的正常位水准面方程：联立求解得：

将分母展开，并略去μ、q平方以上各高次项，就得到一个旋转体，其表面是一个水准面：五、正常重力公式正常重力位对于正常水准面法线（n）的变化量即为正常重力:(忽略n与向量r的区别)

顾及r与a的关系得：特例：，赤道正常重力：

，极点处正常重力：设重力扁率为β克莱罗定理（一般正常重力公式）：精确的正常重力公式：正常重力公式反映地球扁率与重力扁率的关系用不同的观测数据，导出的正常重力公式：1901～1909年赫尔默特公式：1930年卡西尼公式：1975年国际地球正常重力公式：ＷGS84坐标系中的椭球重力公式：

不推导给出高出水准椭球面H米的正常重力计算公式

大约每升高3m，重力值减少1mGal六、正常重力场参数在物理大地测量中,正常椭球重力场可用4个基本参数决定,即：

地球正常(水准)椭球的基本参数，又称地球大地基准常数是：

其中：用地球大地基准常数表示的正常重力位式中：P2n为主球谐系数，J2n为J2的闭合表达式；J2与地球扁率α有如下关系：ae为椭球长半轴。总结引力位：离心力位：重力位：正常重力位（前3项）：与大地水准面相近的正常位水准面：大地水准面正常重力：赤道正常重力：极点正常重力：总结旋转椭球体为我们提供了一个非常简单而又精确的地球几何形状的数学模型，使一些公式推导与计算很方便。为了地面上以铅垂线为依据的观测数据归算到椭球面上，还必须给这个椭球模型加上密合于实际地球的引力场。为此，我们首先把旋转椭球赋予与实际地球椭球相等的质量，同时假定它与地球一起旋转，进而用数学约束条件把椭球面定义为其本身重力场中的一个等位面，并且这个重力场中的铅垂线方向与椭球面相垂直，由此决定的旋转椭球的重力场称为正常重力场。这样的椭球称为正常椭球，也称为水准椭球。结束再见!扇谐函数田谐函数+-1954年北京坐标系的特点1954年北京坐标系坐标系源自于原苏联采用过的1942年普尔科沃坐标系，属参心坐标系。该坐标系采用的参考椭球是克拉索夫斯基椭球，与现代精确的椭球参数的差异较大。还是一个只有几何量表示的椭球。该椭球并未依据当时我国的天文观测资料进行重新定位，而是直接由前苏联西伯利亚地区的一等锁，经我国的东北地区传算过来的。椭球定向不十分明确，椭球的短半轴不指向国际通用的CIO极。参考椭球面与我国大地水准面呈西高东低的系统性倾斜，东部高程异常达60余米，最大达68米。该坐标系统的大地点坐标是经过局部分区平差得到的，即没有进行整体平差。区与区之间存在较大的隙距。1980年国家大地坐标系特点采用国际大地测量协会1975年推荐的参考椭球IAG-75国际椭球，有四个几何和物理参数。椭球的短轴平行于地球的自转轴（由地球质心指向1968.0JYD地极原点方向），起始子午面平行于格林尼治平均天文子午面。属参心坐标系。按照椭球面与似大地水准面在我国境内符合最好的约束条件进行定位，大地原点确定在我国中部。高程系统以1956年黄海平均海水面为高程起算基准。在1980年国家大地坐标系中的大地点成果是经过整体平差。新1954年北京坐标系的特点采用克拉索夫斯基椭球参数。是综合BJ54和GDZ80建立起来的参心坐标系。采用多点定位，但椭球面与大地水准面在我国境内不是最佳拟合。定向明确，坐标轴平行与GDZ80相平行，椭球短轴平行于地球质心指向1968.0JYD地极原点方向），起始子午面平行于格林尼治平均天文子午面。大地原点与GDZ80相同，但大地起算数据不同。大地高程基准采用1956年黄海系。与BJ54旧相比，所采用的椭球参数相同，其定位相近，但定向不同。BJ54旧的坐标是局部平差的结果，而BJ54新是GDZ80整体平差结果的转换值。WGS-84世界大地坐标系WGS-84坐标系是目前GPS所采用的坐标系统，GPS卫星所发布的广播星历参数就是基于此坐标系统的。WGS-84坐标系统的全称是WorldGeodicalSystem-84（世界大地坐标系-84），它是一个地心地固坐标系统。WGS-84坐标系统由美国国防部制图局建立，于1987年取代了当时GPS所采用的坐标系统―WGS-72坐标系统而成为GPS的所使用的坐标系统。WGS-84坐标系的坐标原点位于地球的质心，Z轴指向BIH1984.0定义的协议地球极方向，X轴指向BIH1984.0的启始子午面和赤道的交点，Y轴与X轴和Z轴构成右手系.WGS-84坐标系与国际地球参考系（ITRS，最精密的地心地固坐标系）很接近，差异在2cm以内。不同空间直角坐标系转换一般存在3个平移参数和3个旋转参数以及1个尺度变化参数，共计有7个参数。如5个公共点

一般利用3个以上的公共点求解转换参数引力位的基本性质位函数沿任