大地测量学基础第3章地球重力场及地球型状的基本理论课件

1、1第三章地球重力场及形状的基本理论(续2)1、大地水准面差距和垂线偏差2、确定大地水准面的理论和方法3、重力归算4、区域大地水准面精化5、地球形状的确定方法1第三章1、大地水准面差距和垂线偏差2 选择一个形状和大小都已知的质体(称：正常椭球)，并要求满足如下条件（物理特性），问题转化为： 1）质量和地球相同；旋转角度值与地球相同; 2）尽量不改变大地水准面形状； 3）外部重力场尽量接近地球外部重力场； 4）椭球表面是一个重力等位面（水准面），且恰好包裹全部地球质量； 5）椭球的质心与地球质心重合一、大地水准面差距和垂线偏差V |s = V0lim V = 01、正常椭球Lanplace 算子作

2、用与地球外部重力场=02 选择一个形状和大小都已知的质体(称：3 补充说明： 1、理论上除了确定其M和值外，其规则形状可以任意选择。但考虑到实际使用的方便，又顾及几何大地测量中采用旋转椭球的实际情况，目前都采用水准椭球作为正常椭球。 2、对于正常椭球，除了确定其4个基本参数：a, J，fM和外，也要定位和定向。正常椭球的定位是使其中心和地球质心重合，正常椭球的定向是使其短轴与地轴重合，起始子午面与起始天文子午面重合。 3、正常椭球面 是大地水准面的规则形状（一般指旋转椭球面）。因此引入正常椭球后，地球重力位被分成正常重力位和扰动位两部分，实际重力也被分成正常重力和重力异常两部分。一、大地水准面

3、差距和垂线偏差3 补充说明：一、大地水准面差距和垂线偏差42、正常椭球重力位 将地球重力位展开为球谐级数，然后取前面几项作为正常重力位：当选取前3项时： 并假设：正常椭球是旋转椭球体；坐标系原点在地球质心；坐标轴为地球主惯性轴；赤道上重力约等于 G M/a2:一、大地水准面差距和垂线偏差42、正常椭球重力位一、大地水准面差距和垂线偏差5一、大地水准面差距和垂线偏差5一、大地水准面差距和垂线偏差63、扰动位 大地水准面上一点P的实际重力位 与相应于点P的正常重力位 U 之差，称之为该点的扰动位T，用下式表示：由于在选择正常重力位时总是使地球离心力位对、U的影响相同，因此扰动位具有引力位的性质。一

4、、大地水准面差距和垂线偏差63、扰动位由于在选择正常重力位时总是使地球离心力位对、U74、大地水准面差距由假定：Wp = Upo布隆公式：将几何量N与物理量T联系起来一、大地水准面差距和垂线偏差74、大地水准面差距由假定：Wp = Upo布隆公式：一、大85、垂线偏差通过对Np在子午圈方向和夘酉圈方向求导，球可以得出垂线偏差分量（按弧度表示）：M：子午圈半径；N：夘酉圈半径一、大地水准面差距和垂线偏差85、垂线偏差通过对Np在子午圈方向和夘酉圈方向求导，球可以96、 边值问题线性化一、大地水准面差距和垂线偏差96、 边值问题线性化一、大地水准面差距和垂线偏差10 将p点正常重力展开为P0点的泰

5、劳级数，并代入上式。同时用球面进行近似：一、大地水准面差距和垂线偏差10 将p点正常重力展开为P0点的泰劳级数，并代入上11 如果令椭球的离心力位等于地球的离心力为，则待解的T 微分方程可以结算出来，因为它可以看作是地球外部空间两个引力位之差，所以T 是一个调和函数，T 可以球谐级数展开1、斯托克司（Stokes）解法 如果我们通过重力归算，将大地水准面以外的地形质量去掉，边值条件变为：二、确定大地水准面的理论和方法11 如果令椭球的离心力位等于地球的离12代入边值条件得“理论重力异常”的球谐表达式： 12代入边值条件得“理论重力异常”的球谐表达式： 13 斯托克司积分公式 若已知大地水准面上

6、的重力异常，可通过Stokes积分求解地球外部扰动位T或二、确定大地水准面的理论和方法13斯托克司积分公式或二、确定大地水准面的理论和方法14垂线偏差：维宁.曼尼兹公式： 此公式是在假定大地水准面之外没有扰动物质及全球重力异常都已知的情况下推导的。然而这两个条件都还不能实现，所以重力方法至今也没有得到独立的应用。二、确定大地水准面的理论和方法14垂线偏差：维宁.曼尼兹公式： 此公式是在假定大152、Molodensky方法 Stokes法假定正常椭球之外没有质量，所有观测之必须在属于大地水准面上，而实际的大地测量工作都是在地球自然表面上，这就意味着要知道大地水准面之上的物质密度。为避免这种状况

7、，Molodensky提出采用重力方法直接确定大地水准面的方法。 lim V = 0二、确定大地水准面的理论和方法152、Molodensky方法lim V = 0二、确定大16Molodensky 公式 注意：二、确定大地水准面的理论和方法16Molodensky 公式注意：二、确定大地水准面的理论173、Stokes与Molodensky关系 1） Stokes 要求大地水准面外部没有质量，这就需要将重力观测值归算到大地水准面上，从而需要对密度进行假设。 Molodensky避免了复杂的重力归算。2） Molodensky的零次解就是Stokes解； Molodensky的一次解相当于St

8、okes解+局部地形改正；3） Stokes确定的是大地水准面， Molodensky确定的是大地水准面。和 Stokes和Molodensky对应的正高和正常高。 二、确定大地水准面的理论和方法173、Stokes与Molodensky关系二、确定大地水18重力归算概述 在确定大地水准面形状的基本原理中，有两个前提，一是大地水准面外部必须没有质量；二是所用的实测重力值g 应当是大地水准面上的数值g0。但事实上大地水准面外部有大陆存在，而观测也是在地面上进行的。Monodensky也需要局部地形改正。 重力归化的三个主要目的： （1）求定大地水准面； （2）内插和外推重力值（需要先移去高频变化

9、，然后再恢复）； （3）研究地壳状态。三、重力归算18重力归算概述三、重力归算19三、重力归算重力归化包括以下步骤： 首先将大地水准面外部的地形质量全部去掉，或者移到大地水准面以下去，然后再将重力测量结果从地面降低到大地水准面上。19三、重力归算重力归化包括以下步骤：201空间改正 假设在图中，A 为地面上一点，A0 为大地水准面上相应的投影点，A点的高程为H，我们要将A点的重力加以改正归算到大地水准面上，求出A0 点的重力值。现求其改正数。三、重力归算空间重力异常201空间改正三、重力归算空间重力异常21 2层间改正 将这两个平面间的质量去掉，这样引起的重力改正叫层间改正，记为 。由于质量层

10、在重力观测点的下方，去掉它时重力的值必然减小，所以 总是为负的。两个平行平面间的质量层可被看作半径趋于无穷大的圆柱体，而此时的重力观测点可被当作位于圆柱体项面的中心，所以：布格重力异常（不完全）21 2层间改正布格重力异常22 3. 局部地形改正 由于中间层的上表面并不是平面，而是曲面。进行地形改正 ，就是要把地面凸起部分削去，把凹陷部分填平。 部分：它的引力指向上方，由于它的存在使A点的重力减小了，去掉它必然使A点的重力增大，所以 取正值； 部分：该区域内本来没有质量，填进的质量在A点的引力指向下方，所以使A点的重力增大， 也取正值。 22 3. 局部地形改正23将每一个小部分看作是两个四棱

11、柱的差，可得：i23将每一个小部分看作是两个四棱柱的差，i243gi 将每一个小部分看作是两个扇形柱的差，可得：243gi 将每一个小部分看作是两个扇形柱的差，25 4.布格重力异常（1）布格改正我们称 为完全的布格改正。（2）布格重力异常加入布格改正后计算的重力异常称为布格重力异常。即 5.法耶重力异常（1）法耶改正我们称 为法耶改正。 （2）法耶重力异常称 为法耶重力异常。25 4.布格重力异常2626277、重力均衡和均衡异常-均衡问题的产生277、重力均衡和均衡异常-均衡问题的产生28 如果大地水准面以上的质量是引起重力异常的主要原因，那么布格重力异常中去掉了大地水准面以上质量的影响，

12、它应该很小，但实际情况恰恰相反，山区的布格重力异常一般都量级很大，而且是负的，这说明山区下面的质量有亏损。 垂线偏差的测量结果也表明类似的结论，真正的垂线偏差比看得见的地形质量引起的垂线偏差小得多这也只能用山区下面存在质量亏损来解释。 现在一般比较认可的是地壳均衡理论。28 如果大地水准面以上的质量是引起重力异常的主29 1）普拉特海福特系统（均衡改正）bukan 这种补偿系统的概念由普拉特提出，后来由海福特引进数学公式，系统地应用于大地测量。如图所示：在大地水准面下方某一深度处存在一个补偿面，该面下方的密度是均匀的，上方的质量在该面上产生的压强处处相等。29 1）普拉特海福特系统（均衡改正）

13、buk30 用D表示由大地水准面到补偿面的距离，叫补偿深度，H表示地面到大地水准面的距离，即为海拔高程，大地水准面上方的地壳密度用表示，一般取2.67gcm3，大地水准面与补偿面之间地壳的密度用表示。若H=0，则大地水准面下方的质量无需补偿，。在大陆地区，设两个柱体和的底面积相等，则它们的质量也相等，所以由此可求得补偿密度为30 用D表示由大地水准面到补偿面的距离，叫补偿31在海洋地区，设海洋面与大地水准面重合，P为海洋深度，取海水密度为1.03gcm3，则这可由比较和两个柱体得到，并且可以解得补偿密度为31在海洋地区，设海洋面与大地水准面重合，P为海洋深度，取海32 均衡改正就是补偿密度为0

14、的物质对重力观测值的改正，因此，其计算方法与地形改正相似。 只需将上式中的换成 0，将z的积分限由0到H换成H到HD。均衡重力异常32 均衡改正就是补偿密度为0的物质对重力观测值33 2） 爱黎海斯卡涅系统 这种系统由爱黎提出，海斯卡涅导出了实用于大地测量的公式。其原理如图所示， 地球表层的密度是均匀的，它浮在密度较大的底层上，底层密度为 1，一般取2.67gcm3， 13.27 gcm3， T为大陆地区海拔高程H等于零时表层的厚度，高出大地水准面的大陆质量由深于T的质量亏损得到补偿，这相当于在山的下面有一山根一样，山越高，山根越深。33 2） 爱黎海斯卡涅系统高出大地水准34用t表示表层低于

15、T的深度，则在大陆地区有均衡改正计算：积分区域在 - (h+T) 到 -（h+Tt）34用t表示表层低于T的深度，则在大陆地区有均衡改正计算：积35中国区域空间异常中国区域布格异常中国区域均衡异常中国区域地形图中国区域地形、异常图35中国区域空间异常中国区域布格异常中国区域均衡异常中国区域361、GPS高程拟合法 如果在测区中选择一定的GPS点同时联测几何水准测量，求出这些点的正常高h，于是在这些点上便可求出高程异常:代入下面数学拟合方程中用最小二乘求解： 四、区域大地水准面精化或或361、GPS高程拟合法四、区域大地水准面精化或或372、移去-恢复法确定“重力似大地水准面” 主要思路是：利用

16、地球重力场模型将对大地水准面中长波信息移去，利用重力归算将大地水准面的短波信息移去，再用剩余的重力异常（随机量）进行拟合。提高剩余异常的光滑度，达到提高精度的目的。 四、区域大地水准面精化372、移去-恢复法确定“重力似大地水准面”四、区域大地水准38 第一次移去-恢复：计算出基础格网地面平均空间异常：通过重力归算将重力观测值归算到大地水准面上；移去正常重力值，获得高精度的“地形均衡异常”；拟合出区域“地形均衡异常” 的关系式；通过推估内插、形成平均“地形均衡异常”格网基础数据；再 利用高分辨率DEM，按重力归算的逆过程分别减去：层间改正、局部地形改正、均衡改正；恢复每个格网点的地面平均“空间

17、异常”。四、区域大地水准面精化38 第一次移去-恢复：计算出基础格网地面平均空间异常：39 第二次移去-恢复；计算出重力高程异常和重力似大地面：利用地球重力场模型，由位系数计算出地面格网相同分辨率的模型空间异常；将地面空间异常减去模型空间异常得到格网残差空间异常；在残差空间异常中加上局部地形改正得到残差法耶异常；采用Molodensky积分公式对残差法耶异常进行积分计算，求取每个每个格网中点的残差高程异常。然后利用位系数，由FFT技术计算模型高程异常，并加上残差高程异常得到重力似大地水准面四、区域大地水准面精化39 第二次移去-恢复；计算出重力高程异常和重力似大地面403、组合法-区域大地水准

18、面精化 GPS拟合法确定大地水准面精度高，但难以表征长波效应，分辨率较低。重力法确定大地水准面具有较高的分辨率，但精度较低（一般为分米级）。 组合法结合两者的优点： 首先采用移去-恢复法确定重力似大地水面，然后利用GPS实测高程异常对重力似大地水准面进行纠正。四、区域大地水准面精化403、组合法-区域大地水准面精化四、区域大地水准面精化41组合法-区域大地水准面精化-步骤： 由观测值和相应水准观测值的计算实测高程异常：内插观测点上重力似大地水准面的重力高程异常：；求观测点上，重力高程异常与实测高程异常之差值，并组成不符值序列；由不符值序列和各观测点的球面坐标组成多项式观测方程；按最小二乘法求解多项式系数；由拟合的多项式和格网中心坐标，对重力似大地水准面进行纠正。四、区域大地水准面精化41组合法-区域大地水准面精化-步骤：四、区域大地水准面精化424、利用最小二乘配置法研究大地水准面四、区域大地水准面精化424、利用最小二乘配置法研究大地水准面四、区域大地水准面精43 5、卫星无线电测高方法研究大地水准面四、区域大地水准面精化43 5、卫星无线电测高方法研究大地水准面四、区域大44总的地球椭球： 从几何和物理两个方面来研究全球性问题，我们可把总地球椭球定义为最密合于大地体的正常椭球。 总地球椭球中心和地球质心重合，总的地