正常水准面不平行性及其改正数计算

§5.7 正常水准面不平行性及其改正数计算如果假定不同高程的水准面是相互平行的，那么水准测量所测定的高差，就是水准面间的垂直距离。

这种假定在较短距离内与实际相差不大，而在较长距离时，这种假定是不正确的。

5.7.1 水准面不平行性在空间重力场中的任何物质都受到重力的作用而使其具有位能。

对于水准面上的单位质点而言，它的位能大小与质点所处高度及该点重力加速度有关。

我们把这种随着位置和重力加速度大小而变化的位能称为重力位能，并以表示，则有W（5-gh W =17）式中，为重力加速度；为单位质点所处的高度。

g h 我们知道，在同一水准面上各点的重力位能相等，因此，水准面称为重力等位面，或称重力位水准面。

如果将单位质点从一个水准面提高到相距的另一个水准面，其所做功就h ∆等于两水准面的位能差，即。

在图5-51中，设h g W ∆=∆、分别表示两个非常接近的水准面在两点的垂A h ∆B h ∆B A ,直距离，、为两点的重力加速度。

由于水准面具A gB g B A ,有重力位能相等的性质，因此两点所在水准面的位能差应有下列关系B A ,W ∆（5-B B A A h g h g W ∆=∆=∆18）我们知道，在同一水准面上的不同点重力加速度值是不同的，因此由式（5-g 18）可知，与必定不相等，也就是说，任何两邻近的水准面之间的距离在不A h ∆B h ∆同的点上是不相等的．并且与作用在这些点上的重力成反比。

以上的分析说明水准面不是相互平行的，这是水准面的一个重要特性，称为水准面不平行性。

重力加速度值是随纬度的不同而变化的，在纬度较低的赤道处有较小的值，g g 而在两极处值较大，因此，水准面是相互不平行的、且为向两极收敛的、接近椭圆g 形的曲面。

水准面的不平行性，对水准测量将产生什么影响呢？我们知道，水准测量所测定的高程是由水准路线上各测站所得高差求和而得到的。

在图5-34中，地面点的高程可以按水准路线各测站测得高差之和B OAB ,,21h h ∆∆图5-33求得，即∑∆=OABB hH 测如果沿另一条水准路线施测，则点的高程应为水准路线各测站测得高差ONB B ONB 之和，即 ,,21h h '∆'∆∑'∆='ONBB h H 测由水准面的不平行性可知，因此也必定不等，也就是说，用水∑∑'∆≠∆ONBOABh h B H 测'准测量测得两点间高差的结果随测量所循水准路线的不同而有差异。

如果将水准路线构成闭合环形，既然，可见，即使水准测量完OABNO B BH H 测测'≠全没有误差，这个水准环形路线的闭合差也不为零。

在闭合环形水准路线中，由于水准面不平行所产生的闭合差称为理论闭合差。

由于水准面的不平行性，使得两固定点间的高差沿不同的测量路线所测得的结果不一致而产生多值性，为了使点的高程有惟一确定的数值，有必要合理地定义高程系，在大地测量中定义下面三种高程系统：正高，正常高及力高高程系。

5.7.2 正高高程系正高高程系是以大地水准面为高程基准面，地面上任一点的正高高程（简称正高），即该点沿垂线方向至大地水准面的距离。

如图5-34中，点的正高，设以表示，B BH 正则有（5-⎰∑=∆=BCBCB dH H H 正19）设沿垂线的重力加速度用表示，在垂BC B g 线的不同点上，也有不同的数值。

由式BC B g (5-18)的关系可以写出gdh dH g B =或（5-dh g gdH B=20）将（5-20）式代入（5-19）式中，得图5-34（5-⎰⎰==OABBBCBdh g gdH H 正21）如果取垂线上重力加速度的平均值为，上式又可写为BC B m g（5-⎰=OABB mB gdh g H 1正22）从（5-22）式可以看出，某点的正高不随水准测量路线的不同而有差异，这是因为B 式中为常数，为过点的水准面与大地水准面之间的位能差，也不随路线而Bm g ⎰gdh B 异，因此，正高高程是惟一确定的数值，可以用来表示地面的高程。

如果沿着水准路线每隔若干距离测定重力加速度，则（5-22)式中的值是可以得g 到的。

但是由于沿垂线的重力加速度不但随深入地下深度不同而变化，而且还BC B g 与地球内部物质密度的分布有关，所以重力加速度的平均值并不能精确测定，也不B m g 能由公式推导出来，所以严格说来，地面一点的正高高程不能精确求得。

5.7.3 正常高高程系将正高系统中不能精确测定的用正常重力代替，便得到另一种系统的高程，B m g Bm γ称其为正常高，用公式表达为（5-⎰=gdh H BmBγ1常23）式中，g 由沿水准测量路线的重力测量得到；是水准测量的高差，是按正常重力dh B m γ公式算得的正常重力平均值，所以正常高可以精确求得，其数值也不随水准路线而异，是惟一确定的。

因此，我国规定采用正常高高程系统作为我国高程的统一系统。

下面推导正常高高差的实际计算公式。

首先推导高出水准椭球面的正常重力的计算公式。

在这里，我们把水准椭球m H 看成是半径为的均质圆球，则地心对地面高的点的引力为R H 2)(H R Mfg +=对大地水准面上点的引力为20R M fg =两式相减，得重力改正数))(11(2201H R R fM g g g +-=-=∆)1(11[22RH R fM+-=上式右端括号外项，可认为是地球平均正常重力；由于<<，可把2R fM0γH R 展开级数，并取至二次项，经整理得2)1(-+RH )]321(1[2201RH R H g +--=∆γ220032R H R Hγγ-=将地球平均重力及地球半径代入上式，最后得0γR 2711072.03086.0H H g -⨯-=∆这就是对高出地面点的重力改正公式，式中以m 为单位，以mGal 为单位。

H H g 1∆显然式中第一项是主项，大约每升高3m ，重力值减少1mGa1。

第二项是小项，只在特高山区才顾及它，在一般情况下可不必考虑，这样通常可把上式写成Hg 3086.01=∆于是得出地面高度处的点的正常重力计算公式H（5-H 3086.00-=γγ24）式中为水准椭球面上的正常重力值，在大地控制测量中，采用1901～1909年赫尔0γ默特正常重力公式：（5-)2sin 000007.0sin 005302.01(030.978220ϕϕγ-+=25）将重力写成下面的形式g（5-γγγγ-+-+=Bm B m g g 26）式中用（5-24）式计算。

在有限路线上，可以认为正常重力是线性变化，因此可认γ为是处的值，即，进而B m γB H γ)23086.0(0BB BmH ⋅-=γγγγγγ--+⋅--+=)3086.0()23086.0(00H H g g BB Bm （5-)2(3086.0)()(00H H g BB B m -+-+-+=γγγγ27）分项积分得到⎰⎰⎰-=-OABOABB OABB Hdhdh Hdh H H 2)2(可近似地写成：022)2(22=-=-⎰B B OABB H H dh H H 因此，有正常高计算公式：（5-⎰⎰⎰-+-+=OABB mB B mOABBdh g dh dh H )(1)(100γγγγγ常28）上式右端第一项是水准测量测得的高差，这是主项；第二项中的是沿水准0γB A O --路线上各点的正常重力值，随纬度而变化，亦即，所以第二项称为正常位水准B 00γγ≠面不平行改正数。

第一、二项之和称为概略高程。

第三项是由正常位水位面与重力等位面不一致引起的，称之为重力异常改正项。

当计算两点高差时，有式（5-⎥⎥⎦⎤⎢⎢⎣⎡---+⎥⎥⎦⎤⎢⎢⎣⎡---+=-⎰⎰⎰⎰⎰OBOA AmB m OBOA A AmB Bm AB A B dh g dh g dh dh dh H H )(1)(1)(1)(10000γγγγγγγγγγ常常29）将上式右端第二、三大项分别用和表示，则ελ（5-λε++=-⎰ABAB dh H H 常常30）上式中称为正常位水准面不平行引起的高差改正，称为由重力异常引起的高差改ελ正，经过和改正后的高差称为正常高高差。

ελ 下面推导和的计算公式。

首先推导的计算公式。

ελε由于⎰⎰---=OAA A mOBB B mdhdh )(1)(10000γγγγγγε⎰⎰⎰⎰---+-+---=OAA Am OAB A A Bm OAB Bm OB B Bm dhdh dh dh )(1)(1)(1)(10000000000γγγγγγγγγγγγγγ⎰⎰⎰⎰---+-+-=OAA Am OAA Bm OAB A Bm AB B Bm dhdh dh dh )(1)(1[)(1)(100000000γγγγγγγγγγγγ于是（5-⎰⎰-⋅-+-+-=OAA Bm A m BmA m AB m B A ABB B mdh H dh )()(1000000γγγγγγγγγγγγε31）上式中最后一项数值很小，可略去；第一项在间距不大的情况下，可认为呈线B A ,0γ性变化，可用平均值代替，亦即，则0γ)(21000B Aγγγ+=（5-2)()2(1)(1000000hdh dh Bm A B ABB BA B mABB B m∆⋅--=-+=-⎰⎰γγγγγγγγγγ32）这样 （5-m BmA B A B m A B H H h γγγγγγε)()2()(0000--=+∆--=33）式中，为两点平均高度（可用近似值代替），。

又由（5-25）式m H B A ,γγγ∆=-A B 00可知，若忽略右端第三项（即含项），并令，则把它改写ϕ2sin 2ϕϕ2cos 21212sin 2-=成)]2cos 2121(1[0ϕβγγ-+=e （5-]2cos 21211[ϕββγ-+=e 34）当时，得。

因此上式可写成 45=ϕ)211(45βγγ+=e )2cos 21(45450ϕγγβγγe⋅-=将有关数值代入，于是（5-)2cos 002644.01(9806160ϕγ-⨯=35）因此对上式取微分得ρϕϕγ''⨯⨯=22sin 2002644.09806160d 亦即（5-ϕϕγ'∆⨯=∆2sin 508344.136）当（5-33）式中的以我国平均纬度代入算得Bm γ 35=ϕ)70cos 002644.01(980616 -⨯=Bm γ979773=将以上关系式及数据代入（5-33）式，得的最后计算公式：ε（5-m m H ϕϕε'∆⋅-=2sin 0000015395.037）或 （5-m H A ⋅'∆-=ϕε38）式中，是两点平均纬度，系数可按在水准测量规范中查取，m ϕB A ,A m ϕ是两点的纬度差，以分为单位。