正常重力

正常重力（英語：Normal gravity）是正常椭球体在其外部空间所产生的重力，由意大利数学物理学家卡洛·索米里安在1929年引入，在大地测量学与地球物理学的研究中常用于对真实地球所产生的重力进行近似。在正常重力场中，正常椭球所产生的重力位和能够以较为简单的函数关系表达，且与真实的地球重力位相接近，而正常重力即为这一正常重力位所对应的重力。根据不同的定义方式，真实重力与正常重力之间的差异被称为重力异常或重力扰动。正常重力与真实重力之间的比例约为 99.995 % {\displaystyle 99.995\%} 。 由于正常重力能够被精确计算，其在高程系统中也用于代替真实重力来作为正常高系统所采用的测量值。 正常重力值在两极最大，在赤道处最小，随纬度降低呈递减趋势，相对于赤道面对称而与经度无关。椭球面上几个特殊的重力值分别为： 设正常椭球体在其外部空间产生的正常重力位为 U {\displaystyle U} ，则正常重力矢量被定义为该正常重力位的梯度： γ = ∇ U {\displaystyle {\boldsymbol {\gamma }}=\nabla U} 在椭球坐标系 ( u , β , λ ) {\displaystyle (u,\beta ,\lambda )} 中，正常重力矢量的三个分量具体表示为： γ u = 1 w ∂ U ∂ u {\displaystyle \gamma _{u}={1 \over w}{\partial U \over \partial u}} γ β = 1 w u 2 + E 2 ∂ U ∂ β {\displaystyle \gamma _{\beta }={1 \over w{\sqrt {u^{2}+E^{2}}}}{\partial U \over \partial \beta }} γ λ = 1 u 2 + E 2 cos ⁡ β ∂ U ∂ λ = 0 {\displaystyle \gamma _{\lambda }={1 \over {\sqrt {u^{2}+E^{2}}}\cos \beta }{\partial U \over \partial \lambda }=0} 上式中的 w = u 2 + E 2 sin 2 ⁡ β u 2 + E 2 {\displaystyle w={\sqrt {u^{2}+E^{2}\sin ^{2}\beta \over u^{2}+E^{2}}}} 是为简化公式而引入的辅助量， E {\displaystyle E} 是椭球的半焦距。又因正常重力位 U {\displaystyle U} 与经度无关，所以正常重力矢量的经度分量为零。 由正常重力的数学表达式可以得出，正常重力的值可以根据正常重力位 U {\displaystyle U} 的偏导数，以及正常椭球体本身的几何性质得到。而正常椭球体的确定只需要四个基本参数：椭球的半长轴 a {\displaystyle a} 、几何扁率 f {\displaystyle f} 、赤道上的正常重力值 γ e {\displaystyle \gamma _{e}} ，以及地球自转的角速度 ω {\displaystyle \omega } ，其他的几何参数可以由上述基本参数确定： 椭球的半短轴 b = a ( 1 − f ) {\displaystyle b=a(1-f)} 椭球的第一偏心率 e = a 2 − b 2 / a = 2 f − f 2 {\displaystyle e={\sqrt {a^{2}-b^{2}}}/a=2f-f^{2}} 椭球的第二偏心率 e = a 2 − b 2 / b {\displaystyle e={\sqrt {a^{2}-b^{2}}}/b} 亦有一些坐标系统会选择其他的基本参数，例如GRS80椭球选用的是地心引力常数 G M {\displaystyle GM} 、地球动力学形状因子 J 2 {\displaystyle J_{2}} 、地球自转角速度 ω {\displaystyle \omega } 和椭球的半长轴 a {\displaystyle a} ，但其他的椭球参数仍能由这些基本参数计算而得。 法国数学家克莱罗在其发表于1743年的著作中给出了地球的几何扁率 f {\displaystyle f} 与重力扁率 f ∗ {\displaystyle f^{*}} 之间的对应关系，即克莱罗定理。在顾及至扁率的平方项的情况下，该定理可表述为： f + f ∗ = 5 2 ω 2 b γ e ( 1 + 9 35 e ′ 2 ) {\displaystyle f+f^{*}={5 \over 2}{\omega ^{2}b \over \gamma _{e}}(1+{9 \over 35}{e'}^{2})} 重力扁率 f ∗ {\displaystyle f^{*}} 的定义与几何扁率类似，其由椭球赤道处的重力 γ e {\displaystyle \gamma _{e}} 和椭球极点处的重力 γ p {\displaystyle \gamma _{p}} 决定 ： f ∗ = γ p − γ e γ e {\displaystyle f^{*}={\gamma _{p}-\gamma _{e} \over \gamma _{e}}} γ e = G M a 2 ( 1 + m + 3 7 e ′ 2 m ) {\displaystyle \gamma _{e}={GM \over a^{2}}\left(1+m+{3 \over 7}{e'}^{2}m\right)} γ p = G M a b ( 1 − 3 2 m − 3 14 e ′ 2 m ) {\displaystyle \gamma _{p}={GM \over ab}\left(1-{3 \over 2}m-{3 \over 14}{e'}^{2}m\right)} 其中 m = ω 2 a 2 b G M {\displaystyle m={\omega ^{2}a^{2}b \over GM}} ，且有 ω 2 b γ e = m + 3 2 m 2 {\displaystyle {\omega ^{2}b \over \gamma _{e}}=m+{3 \over 2}m^{2}} 。 克莱罗定理给出了椭球赤道处的正常重力值和极点处的正常重力值，而椭球面上其他纬度的正常重力则可由正常重力公式计算得到，这一公式由索米里安在1929年给出： γ = a γ p sin 2 ⁡ β + b γ e cos 2 ⁡ β a 2 sin 2 ⁡ β + b 2 cos 2 ⁡ β {\displaystyle \gamma ={a\gamma _{p}\sin ^{2}\beta +b\gamma _{e}\cos ^{2}\beta \over {\sqrt {a^{2}\sin ^{2}\beta +b^{2}\cos ^{2}\beta }}}} 其中 β {\displaystyle \beta } 是椭球面上某点的归化纬度，顾及到大地纬度 φ {\displaystyle \varphi } 与归化纬度 β {\displaystyle \beta } 存在如下转换关系： tan ⁡ β = b a tan ⁡ φ {\displaystyle \tan \beta ={b \over a}\tan \varphi } 则正常重力公式也可以表达成大地纬度 φ {\displaystyle \varphi } 的函数： γ = a γ e cos 2 ⁡ φ + b γ p sin 2 ⁡ φ a 2 cos 2 ⁡ φ + b 2 sin 2 ⁡ φ {\displaystyle \gamma ={a\gamma _{e}\cos ^{2}\varphi +b\gamma _{p}\sin ^{2}\varphi \over {\sqrt {a^{2}\cos ^{2}\varphi +b^{2}\sin ^{2}\varphi }}}} 正常重力公式也可以展开为几何扁率 f {\displaystyle f} 的级数，其截断形式为： γ = γ e ( 1 + f 2 sin 2 ⁡ φ + f 4 sin 4 ⁡ φ ) {\displaystyle \gamma =\gamma _{e}(1+f_{2}\sin ^{2}\varphi +f_{4}\sin ^{4}\varphi )} 其中的系数为： f 2 = − f + 5 2 m + 1 2 f 2 − 26 7 f m + 15 4 m 2 {\displaystyle f_{2}=-f+{5 \over 2}m+{1 \over 2}f^{2}-{26 \over 7}fm+{15 \over 4}m^{2}} f 4 = − 1 2 f 2 + 5 2 f m {\displaystyle f_{4}=-{1 \over 2}f^{2}+{5 \over 2}fm} 这一公式也可写为： γ = γ e ( 1 + f ∗ sin 2 ⁡ φ − 1 4 f 4 sin 4 ⁡ 2 φ ) {\displaystyle \gamma =\gamma _{e}(1+f^{*}\sin ^{2}\varphi -{1 \over 4}f_{4}\sin ^{4}2\varphi )} 其中的 f ∗ = f 2 + f 4 {\displaystyle f^{*}=f_{2}+f_{4}} 为上述提到的重力扁率。 正常重力公式还可以闭合形式表达： γ = γ e 1 + k sin 2 ⁡ φ 1 − e 2 s i n 2 φ {\displaystyle \gamma =\gamma _{e}{1+k\sin ^{2}\varphi \over {\sqrt {1-e^{2}sin^{2}\varphi }}}} 其中的系数 k {\displaystyle k} 为： k = b γ p − a γ e a γ e {\displaystyle k={b\gamma _{p}-a\gamma _{e} \over a\gamma _{e}}} 采用不同的椭球参数和不同的表达形式，正常重力公式可以有不同的数值计算形式，常用的几条公式包括： 在椭球面外部不远处，其正常重力 γ h {\displaystyle \gamma _{h}} 可以在其沿法线到椭球面上投影处展开为正常高 h {\displaystyle h} 的级数： γ h = γ + ∂ γ ∂ h h + 1 2 ∂ 2 γ ∂ h 2 h 2 + ⋯ {\displaystyle \gamma _{h}=\gamma +{\partial \gamma \over \partial h}h+{1 \over 2}{\partial ^{2}\gamma \over \partial h^{2}}h^{2}+\cdots } 由广义布隆斯方程，椭球面的外部空间的重力梯度与椭球面（水准面）的平均曲率半径 J {\displaystyle J} 的关系为： ∂ γ ∂ h = − 2 γ J − 2 ω 2 = − 2 γ a ( 1 + f + m − 2 f sin 2 ⁡ φ ) {\displaystyle {\partial \gamma \over \partial h}=-2\gamma J-2\omega ^{2}=-{2\gamma \over a}\left(1+f+m-2f\sin ^{2}\varphi \right)} 又二次导数 ∂ 2 γ / ∂ h 2 {\displaystyle \partial ^{2}\gamma /\partial h^{2}} 是微小量，可以将其近似近似于在球面外部微分（即以半长轴 a {\displaystyle a} 代替 r {\displaystyle r} ），得到： ∂ 2 γ ∂ h 2 = 6 G M a 4 = 6 γ a 2 {\displaystyle {\partial ^{2}\gamma \over \partial h^{2}}={6GM \over a^{4}}={6\gamma \over a^{2}}} 得到正常重力的向上延拓公式为： γ h = γ [ 1 − 2 a ( 1 + f + m − 2 f sin 2 ⁡ φ ) h + 3 a 2 h 2 ] {\displaystyle \gamma _{h}=\gamma \left[1-{2 \over a}\left(1+f+m-2f\sin ^{2}\varphi \right)h+{3 \over a^{2}}h^{2}\right]} 上式的数值形式近似为： γ h = γ − 0.3086 h + 0.72 × 10 − 7 h 2 {\displaystyle \gamma _{h}=\gamma -0.3086h+0.72\times 10^{-7}h^{2}} 重力异常 布隆斯公式 似大地水准面 其中 u {\displaystyle u} 表示椭球的半短轴， β {\displaystyle \beta } 表示归化纬度， λ {\displaystyle \lambda } 表示经度 GeographicLib: Normal gravity. geographiclib.sourceforge.io. [2020-04-13]. （原始内容存档于2021-05-06）.  Somigliana, Carlo. Teoria generale del campo gravitazionale dell'ellissoide di rotazione. Libr. Editr. Politecnica. 1929 （意大利语）.  宁津生. 管泽霖 , 编. 地球形状及外部重力场. 测绘出版社. 1981: 154–293.  Jekeli, C. Potential Theory and Static Gravity Field of the Earth. Treatise on Geophysics. Elsevier. 2007: 11–42 [2020-04-15]. ISBN 978-0-444-52748-6. doi:10.1016/b978-044452748-6.00054-7. （原始内容存档于2018-07-01） （英语）.  孔祥元; 郭际明; 刘宗泉. 大地测量学基础. 武汉大学出版社. 2001. ISBN 978-7-30-707562-7.  Torge, Wolfgang. Geodesy. Walter de Gruyter GmbH & Co KG. 2001. ISBN 978-3-11-017072-6 （英语）.  H. Moritz. Geodetic Reference System 1980 (pdf). （原始内容存档 (PDF)于2021-04-01） （英语）.  San Francisco W. H. Freeman and Company. Heiskanen Moritz 1967 Physical Geodesy. San Francisco: W. H. Freeman and Company. 1967 （英语）.  引文格式1维护：日期与年 (link) Alexis Claude, Clairaut. Théorie de la figure de la terre : tirée des principes de l'hydrostatique (PDF). 1743. [失效連結] Department of Defense World Geodetic System 1984 ― Its Definition and Relationships with Local Geodetic Systems (PDF). [2020-04-14]. （原始内容 (pdf)存档于2021-03-01） （英语）.  Vanícek, P.; Krakiwsky, E. J. Geodesy: The Concepts. Elsevier. 2015-06-03. ISBN 978-1-4832-9079-9 （英语）.