查看“︁正常重力位”︁的源代码

'''正常重力位'''（{{Lang-en|Normal gravity potential}}）是[[大地测量学]]中用于对[[地球]]的真实[[重力位]]进行近似的数学工具，是一个规则的、较为简单的重力位[[函数]]。<ref name=":0">{{Cite book|title=地球形状及外部重力场|author=宁津生|publisher=测绘出版社|year=1981|isbn=|location=|authorlink=宁津生|editor=管泽霖|first=}}</ref>{{Rp|190}}[[正常重力]]是正常重力位的[[梯度]]。<ref name=":2">{{Cite book|chapter=|url=http://archive.org/details/HeiskanenMoritz1967PhysicalGeodesy|date=1967|last=San Francisco W. H. Freeman and Company|title=Heiskanen Moritz 1967 Physical Geodesy|first=|publisher=W. H. Freeman and Company|year=1967|isbn=|location=San Francisco|pages=|language=en}}</ref>{{Rp|68}}

地球的真实重力位在研究地球形状及其外部重力场的过程中是待解的未知量，且对其直接计算需要了解地球内部的质量分布，理论上被无法精确求得，且反演过程计算复杂。因此，研究中常选用一个规则的[[正常椭球体]]对地球的形状进行近似，将其产生的重力场中的重力位称为正常重力位<ref name=":1">{{Cite book|chapter=Potential Theory and Static Gravity Field of the Earth|title=Treatise on Geophysics|url=https://linkinghub.elsevier.com/retrieve/pii/B9780444527486000547|publisher=Elsevier|date=2007|isbn=978-0-444-52748-6|pages=11–42|doi=10.1016/b978-044452748-6.00054-7|language=en|first=C.|last=Jekeli|access-date=2020-04-15|archive-date=2018-07-01|archive-url=https://web.archive.org/web/20180701075456/https://linkinghub.elsevier.com/retrieve/pii/B9780444527486000547|dead-url=no}}</ref>{{Rp|15}}，而把真实重力位与正常重力位的差异称作[[扰动位]]。选择适当的正常重力位，可以使扰动位成为微小量，便于以[[线性近似]]的方式对其进行求解。<ref name=":2" />{{Rp|64}}在地球重力场中，正常重力位可以占到到真实重力位的 <math>99.999 5\%</math><ref name=":1" />{{Rp|15}}。

== 确定要求 ==
正常重力位应当具有以下特性：<ref name=":3">{{Cite book|title=Physical Geodesy|last=Sneeuw|first=Nico|publisher=Institute of Geodesy Universität Stuttgart|year=2006|isbn=|location=|url=https://www.gis.uni-stuttgart.de/lehre/campus-docs/LNErdm.pdf|access-date=2020-04-16|archive-url=https://web.archive.org/web/20200413023934/https://www.gis.uni-stuttgart.de/lehre/campus-docs/LNErdm.pdf|archive-date=2020-04-13|dead-url=yes}}</ref>{{Rp|94}}

* 旋转对称性，即与[[经度]]无关
* 相对于[[赤道|赤道面]]对称
* 在[[椭球|椭球面]]上是[[常数]]

这些特性保证了正常重力位是规则分布的，利用这些性质能够简化复杂的计算过程。最后一个条件保证了还正常椭球体的表面是一个重力[[等位面]]。

== 数学表达 ==
正常重力位  <math>U</math> 包含两部分，一部分是因正常椭球体的[[质量]]而产生的[[引力位]] <math>V</math>，另一部分是因正常椭球体的绕轴自转而产生的离心力位 <math>\Phi</math>：<ref name=":2" />{{Rp|64}}

: <math>
U = V + \Phi = V + {1\over2}\omega^2(x^2+y^2)
</math>

上式中，<math>
\omega
</math> 表示[[地球自转]]的[[角速度]]，<math>
(x,y,z)
</math> 是正常椭球体外部空间中某点的[[笛卡尔坐标]]，该坐标系的 <math>
Z
</math> 轴与椭球体的[[自转轴]]平行或重合。其中的离心力位 <math>\Phi</math> 可以由点的坐标直接求得，不必展开为级数，只需对引力位 <math>V</math> 进行展开即可。<ref name=":0" />{{Rp|197}}

=== 斯托克斯方法 ===
采用球面作为正常椭球体的近似，将椭球体外部（<math>
r>a
</math>）的引力位展开为边界面（ <math>
r = a
</math> ）上的球谐级数，则其表达式为：<ref name=":2" />{{Rp|59}}<ref name=":3" />{{Rp|79}}

:<math>
V =
{\sum_{n=0}^{\infin}}\frac{1}{r^{n+1}}
{\sum_{m=0}^n}P_{nm}(\cos\theta)\left[A_{nm}\cos m\lambda+B_{nm}\sin m\lambda\right]
</math>
:或
:<math>
V = \frac{GM}{r}
{\sum_{n=0}^{\infin}}\left(\frac{a}{r}\right)^n
{\sum_{m=0}^n}P_{nm}(\cos\theta)\left[C_{nm}\cos m\lambda+S_{nm}\sin m\lambda\right]
</math>

上式中各项符号的物理或数学意义如下：

*<math>
(r, \theta, \lambda)
</math>是空间中某特定点的[[球座標系|球坐标]]，<math>r </math> 是空间中某点的地心距离， <math>\theta </math> 和 <math>\lambda </math> 分别是该点的[[极距]]和[[经度]]

*<math>
GM
</math> 为正常椭球体的[[標準重力參數|地心引力常数]]
*<math>a</math> 为椭球的長[[半径]]
*<math>P_n</math> 是 <math>n\ </math>阶[[勒让德多项式]]，<math>P_{nm}</math> 是 <math>n\ </math>阶 <math>m\ </math>次[[缔合勒让德多项式]]
*<math>
A_{nm}
</math> 和 <math>
B_{nm}
</math> 、<math>
C_{nm}
</math> 和 <math>
S_{nm}
</math> 是由[[积分]]求得的球谐系数，且积分表达式为：<ref name=":2" />{{Rp|59}}<ref name=":3" />{{Rp|90}}

:<math>
\begin{Bmatrix}  A_{nm} \\ B_{nm} \end{Bmatrix} = 
(2-\delta_{m,0}){(n-m)!\over(n+m)!}G
\iiint\limits_{\text{earth}} r'^n P_{nm}(\cos\theta')
\begin{Bmatrix}  \cos m\lambda' \\ \sin m\lambda' \end{Bmatrix}\operatorname{d}\!M'
\quad
</math>
:或
:<math>
\begin{Bmatrix}  C_{nm} \\ S_{nm} \end{Bmatrix} = 
(2-\delta_{m,0}){(n-m)!\over(n+m)!}{1 \over M}
\iiint\limits_{\text{earth}} {\left({r' \over a}\right)}^n P_{nm}(\cos\theta')
\begin{Bmatrix}  \cos m\lambda' \\ \sin m\lambda' \end{Bmatrix}\operatorname{d}\!M'
\quad
</math>

其中 <math>
\delta_{m,0}
</math> 为[[克罗内克δ函数]]，当 <math>
m=0
</math> 时为零，其他情况下为一。

展开后的球谐级数具有以下性质：

# 正常重力场具有旋转对称性，即其产生的正常重力位与[[经度]]无关，因此上式中[[面谐函数]]的部分可被简化作 <math>
A_{n,0}P_n(\cos\theta)
</math> 或 <math>
C_{n,0}P_n(\cos\theta)
</math>
# 当坐标系的原点与椭球体的质心重合时，可证明该级数中的一阶系数为零，即 <math>
A_{1,m} = B_{1,m} = 0
</math> 或 '''<math>
C_{1,m} = S_{1,m} = 0
</math>'''
# 保留的阶数根据观测资料的精度以及对正常重力位要求的精度确定<ref name=":0" />{{Rp|207-208}}，一般只需保留至8阶项<ref name=":3" />{{Rp|95}}

根据上述性质，该表达式可进一步简化为：

:<math>
V = {\sum_{n=0}^{8}}\frac{1}{r^{n+1}}A_{n,0}P_n(\cos\theta)
</math>
:或
:<math>
V = \frac{GM}{r}{\sum_{n=0}^{8}}\left(\frac{a}{r}\right)^n C_{n,0}P_n(\cos\theta)
</math>

=== 拉普拉斯方法 ===
利用引力位  <math>V</math>  在边界面外部满足[[拉普拉斯方程]]的性质，将地球的真实引力位可展开为[[球谐函数|球谐级数]]，保留其中的头几项作为正常重力位的引力位部分从而确定正常重力位的方法被称为拉普拉斯方法。<ref name=":0" />{{Rp|190-191}}通过选取不同的大地水准面重力位值，可以得到不同的正常重力位等位面，从中选取一个最接近于大地水准面的，这一曲面即为产生正常重力位的质体的表面。

==== 重力位的椭球谐级数 ====
将正常重力位直接展开成[[椭球面]] <math>S_0</math> 上的级数，称为[[椭球谐函数|椭球谐级数]]，形式较球谐函数更为复杂：<ref name=":2" />{{Rp|65}}

:<math>
U = V + \Phi 
= \sum_{n}^{\infin}{Q_n\left(i{u \over E}\right) \over Q_n\left(i{b \over E}\right)}A_nP_n(\sin\beta) + {1\over2}\omega^2(u^2+E^2)\cos^2\beta
</math>

上式中各项符号的物理或数学意义与球谐函数有所不同：

*<math>
(u, \beta, \lambda)
</math> 是空间中某特定点的[[椭球坐标]]，<math>u </math> 是该坐标所在椭球的[[半短軸|短轴半长]]， <math>\beta </math> 和 <math>\lambda\ </math>分别是该点的[[归化纬度]]和[[经度]]
*<math>a </math>、<math>b </math> 分别是正常椭球面 <math>S_0</math>  的[[半長軸]]和半短轴，<math>E = \sqrt{a^2-b^2} </math> 是椭球的[[线性偏心率]]（即[[半焦距]]）
*<math>Q_n </math>是 <math>n </math> 阶[[第二类勒让德多项式]]

当 <math>u = b </math> 时，由重力位 <math>U_0 </math> 所决定的等位面应当与正常椭球面 <math>S_0</math> 相重合，此时有：

:<math>
U_0 = V_0 + \Phi_0 = \sum_{n}^{\infin}A_nP_n(\sin\beta) + {1\over2}\omega^2a^2\cos^2\beta
</math>

==== 引力位的椭球谐级数 ====
当且仅当所有含 <math>
P_n(\sin\beta)
</math> 的项均为零时，对于任意 <math>
\beta
</math> 值该公式都成立。对其头三项进行展开：<ref name=":2" />{{Rp|66}}

:<math>
\left(A_0+{1\over3}\omega^2a^2-U_0\right)P_0(\sin\beta)
+ A_1P_1(\sin\beta)
+ \left(A_2-{1\over3}\omega^2a^2\right)P_2(\sin\beta)
+ \sum_{n=3}^{\infin}A_nP_n(\sin\beta) = 0
</math>

因此，有

* <math>
A_0=U_0-{1\over3}\omega^2a^2
</math>
* <math>
A_1=0
</math>
* <math>
A_2={1\over3}\omega^2a^2
</math>
* <math>
A_3=A_4=\cdots=0
</math>

得到引力位部分的椭球谐级数表达式为：<ref name=":2" />{{Rp|66}}

:<math>
V 
= \left(U_0-{1\over3}\omega^2a^2\right){Q_0\left(i{u \over E}\right) \over Q_0\left(i{b \over E}\right)}
+ {1\over3}\omega^2a^2{Q_2\left(i{u \over E}\right) \over Q_2\left(i{b \over E}\right)}P_2(\sin\beta)
= \left(U_0-{1\over3}\omega^2a^2\right){\tan^{-1}{E \over u} \over \tan^{-1}{E \over b}}
+ {1\over3}\omega^2a^2{q \over q_0}P_2(\sin\beta)
</math>

其中

* <math>
q={1\over2}\left[\left(1+3\frac{u^2}{E^2}\right)\tan^{-1}\frac{E}{u}-3\frac{u}{E}\right]
</math>
* <math>
q_0={1\over2}\left[\left(1+3\frac{b^2}{E^2}\right)\tan^{-1}\frac{E}{b}-3\frac{b}{E}\right]
</math>

==== 大地水准面重力位 ====
考虑到[[椭圆坐标]] <math>
(u,\beta)
</math> 与向径 <math>
r
</math> 存在如下转换关系：

:<math>
r^2=u^2+E^2\cos^2\beta
</math>

利用该关系对下式进行[[线性化]]，得：

:<math>
\frac{1}{u}=\frac{1}{r}+\Omicron(\frac{1}{r^3}) 
</math>
:<math>
\tan^{-1}\frac{E}{u}=\frac{E}{u}+\Omicron(\frac{1}{u^3})=\frac{E}{r}+\Omicron(\frac{1}{r^3}) 
</math>

利用这两项关系式，可以得到线性化后的引力位函数：

:<math>
V 
= \left(U_0-{1\over3}\omega^2a^2\right){E \over \tan^{-1}(E/b)}\frac{1}{r}
+ \Omicron(\frac{1}{r^3})
</math>

比较与 

:<math>
V = \frac{GM}{r} + \Omicron(\frac{1}{r^3})
</math>

得出大地水准面重力位的表达式：<ref name=":2" />{{Rp|67}}

:<math>
U_0=\frac{GM}{E}\tan^{-1}\frac{E}{b}+{1\over3}\omega^2a^2
</math>

==== 正常重力位计算公式 ====
将大地水准面重力位的表达式代入原正常重力位的计算公式中，得：<ref name=":2" />{{Rp|67}}

:<math>
U(u,\beta)=\frac{GM}{E}\tan^{-1}\frac{E}{u}
+{1\over2}\omega^2a^2\frac{q}{q_0}\left(\sin^2\beta-{1\over3}\right)
+{1\over2}\omega^2(u^2+E^2)\cos^2\beta
</math>

== 物理性质 ==

=== 球谐系数 ===
球谐系数的积分公式中包含了正常椭球体内质量的分布关系，且积分范围是整个正常椭球体，因此球谐系数与正常椭球体的某些物理性质相关。

==== 零阶项 ====
当 <math>n=0</math> 时，球谐系数只有一项：<ref name=":2" />{{Rp|61}}

:<math>
A_{0,0} = G
\iiint\limits_{\text{earth}}
\operatorname{d}\!M' = GM
</math> 
:或
:<math>
C_{0,0} = {1 \over M}
\iiint\limits_{\text{earth}}
\operatorname{d}\!M' = 1
</math>

即球谐系数的零阶项反映了正常椭球体的[[地心引力常数]]或[[质量|总质量]]。

==== 一阶项 ====
当 <math>n=1</math> 时，球谐系数有三项：<ref name=":2" />{{Rp|61}}

:<math>
A_{1,0} = G
\iiint\limits_{\text{earth}} z'
\operatorname{d}\!M' = Gz_M
</math>
:<math>
A_{1,1} = G
\iiint\limits_{\text{earth}} x'
\operatorname{d}\!M' = Gx_M
</math>
:<math>
B_{1,1} = G
\iiint\limits_{\text{earth}} y'
\operatorname{d}\!M' = Gy_M
</math>
:或
:<math>
C_{1,0} = {1 \over Ma}
\iiint\limits_{\text{earth}} z'
\operatorname{d}\!M' = {1 \over a}z_M
</math>
:<math>
C_{1,1} = {1 \over Ma}
\iiint\limits_{\text{earth}} x'
\operatorname{d}\!M' = {1 \over a}x_M
</math>
:<math>
S_{1,1} = {1 \over Ma}
\iiint\limits_{\text{earth}} y'
\operatorname{d}\!M' = {1 \over a}y_M
</math>

其中，<math>
\vec{r}_M = \left(x_M,y_M,z_M\right)
</math> 表示椭球的[[质心]]坐标。当坐标系的原点与椭球质心重合时，<math>
x_M=y_M=z_M=0
</math>，所以球谐系数的一阶项一般都为零。<ref name=":2" />{{Rp|62}}

==== 二阶项 ====
当 <math>n=2</math> 时，球谐系数有五项（仅以 <math>
A_{m,n}
</math> 和 <math>
B_{m,n}
</math> 为例）：<ref name=":2" />{{Rp|62}}

:<math>
A_{2,0} = G
\iiint\limits_{\text{earth}} {1\over2}(2z'^2-x'^2-y'^2)
\operatorname{d}\!M'
</math>
:<math>
A_{2,1} = G
\iiint\limits_{\text{earth}} x'z'
\operatorname{d}\!M'
</math>
:<math>
B_{2,1} = G
\iiint\limits_{\text{earth}} y'z'
\operatorname{d}\!M'
</math>
:<math>
A_{2,2} = G
\iiint\limits_{\text{earth}} {1\over4}(x'^2-y'^2)
\operatorname{d}\!M'
</math>
:<math>
B_{2,2} = G
\iiint\limits_{\text{earth}} {1\over2}x'y'
\operatorname{d}\!M'
</math>

其中的五项积分，既可以[[四极子|四极矩张量]] <math>
\mathbf{Q}
</math> 或[[惯性张量]] <math>
\mathbf{I}
</math> 表达：<ref name=":3" />{{Rp|91}}

:<math>
\iiint\limits_{\text{earth}} {1\over2}(2z'^2-x'^2-y'^2)\operatorname{d}\!M'
= {1\over2}Q_{zz} = {1\over2}(I_{xx}+I_{yy})-I_{zz}
</math>
:<math>
\iiint\limits_{\text{earth}} x'z'\operatorname{d}\!M'
= {1\over3}Q_{xz}
= -I_{xz}
</math>
:<math>
\iiint\limits_{\text{earth}} y'z'\operatorname{d}\!M'
= {1\over3}Q_{yz}
= -I_{yz}
</math>
:<math>
\iiint\limits_{\text{earth}} {1\over4}(x'^2-y'^2)\operatorname{d}\!M'
= {1\over12}(Q_{xx}-Q_{yy})
= {1\over4}(I_{yy}-I_{xx})
</math>
:<math>
\iiint\limits_{\text{earth}} {1\over2}x'y'\operatorname{d}\!M'
= {1\over6}Q_{xy}
= -{1\over2}I_{xy}
</math>

当且仅当坐标系的各坐标轴与地球的[[主惯性轴]]重合时，<math>
I_{xz}=I_{xy}=I_{yz}=0
</math>，因此亦有 <math>
A_{2,1}=B_{2,1}=B_{2,2}=0
</math>。反过来，坐标轴的选择又决定了这三个二阶项系数的值：当坐标系的 <math>
Z
</math> 轴指向[[协议地球极]]时，受[[极移]]等因素的影响，这一指向与地球的瞬时主惯性轴并不重合，因此 <math>
A_{2,1}
</math> 与 <math>
B_{2,1}
</math> 的并不为零；而 <math>
X
</math> 轴的指向（通常是[[本初子午线]]）则决定了 <math>
B_{2,2}
</math> 的数值。<ref name=":3" />{{Rp|91}} 另外，<math>
A_{2,2}
</math> 由[[赤道]]的形状决定：当正常椭球为对称的旋转体时，赤道是圆形，此时 <math>
I_{xx}=I_{yy}
</math>，即有<math>
A_{2,2}=0
</math>。<ref name=":0" />{{Rp|205}}{{NoteFoot}}

== 参考文献 ==
{{Reflist}}{{物理大地测量学}}
[[Category:大地测量学]]
[[Category:地球物理学]]