**正常重力位**（英語：）是大地测量学中用于对地球的真实重力位进行近似的数学工具，是一个规则的、较为简单的重力位函数。正常重力是正常重力位的梯度。

地球的真实重力位在研究地球形状及其外部重力场的过程中是待解的未知量，且对其直接计算需要了解地球内部的质量分布，理论上被无法精确求得，且反演过程计算复杂。因此，研究中常选用一个规则的正常椭球体对地球的形状进行近似，将其产生的重力场中的重力位称为正常重力位，而把真实重力位与正常重力位的差异称作扰动位。选择适当的正常重力位，可以使扰动位成为微小量，便于以线性近似的方式对其进行求解。在地球重力场中，正常重力位可以占到到真实重力位的 {\displaystyle 99.9995\%}。

确定要求
----

正常重力位应当具有以下特性：

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

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

数学表达
----

正常重力位 U 包含两部分，一部分是因正常椭球体的质量而产生的引力位 V，另一部分是因正常椭球体的绕轴自转而产生的离心力位 \Phi ：

{\displaystyle U=V+\Phi =V+{1 \over 2}\omega ^{2}(x^{2}+y^{2})}
上式中，{\displaystyle \omega } 表示地球自转的角速度，{\displaystyle (x,y,z)} 是正常椭球体外部空间中某点的笛卡尔坐标，该坐标系的 {\displaystyle Z} 轴与椭球体的自转轴平行或重合。其中的离心力位 \Phi  可以由点的坐标直接求得，不必展开为级数，只需对引力位 V 进行展开即可。

### 斯托克斯方法

采用球面作为正常椭球体的近似，将椭球体外部（{\displaystyle r>a}）的引力位展开为边界面（ {\displaystyle r=a} ）上的球谐级数，则其表达式为：

{\displaystyle V={\sum \_{n=0}^{\infty }}{\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]}
或
{\displaystyle V={\frac {GM}{r}}{\sum \_{n=0}^{\infty }}\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]}
上式中各项符号的物理或数学意义如下：

* {\displaystyle (r,\theta ,\lambda )}是空间中某特定点的球坐标，{\displaystyle r} 是空间中某点的地心距离， \theta  和 \lambda  分别是该点的极距和经度

* {\displaystyle GM} 为正常椭球体的地心引力常数
* a 为椭球的長半径
* P\_{n} 是 n\ 阶勒让德多项式，{\displaystyle P\_{nm}} 是 n\ 阶 m\ 次缔合勒让德多项式
* {\displaystyle A\_{nm}} 和 {\displaystyle B\_{nm}} 、{\displaystyle C\_{nm}} 和 {\displaystyle S\_{nm}} 是由积分求得的球谐系数，且积分表达式为：

{\displaystyle {\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 }
或
{\displaystyle {\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 }
其中 {\displaystyle \delta \_{m,0}} 为克罗内克δ函数，当 {\displaystyle m=0} 时为零，其他情况下为一。

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

1. 正常重力场具有旋转对称性，即其产生的正常重力位与经度无关，因此上式中面谐函数的部分可被简化作 {\displaystyle A\_{n,0}P\_{n}(\cos \theta )} 或 {\displaystyle C\_{n,0}P\_{n}(\cos \theta )}
2. 当坐标系的原点与椭球体的质心重合时，可证明该级数中的一阶系数为零，即 {\displaystyle A\_{1,m}=B\_{1,m}=0} 或 **{\displaystyle C\_{1,m}=S\_{1,m}=0}**
3. 保留的阶数根据观测资料的精度以及对正常重力位要求的精度确定，一般只需保留至8阶项

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

{\displaystyle V={\sum \_{n=0}^{8}}{\frac {1}{r^{n+1}}}A\_{n,0}P\_{n}(\cos \theta )}
或
{\displaystyle V={\frac {GM}{r}}{\sum \_{n=0}^{8}}\left({\frac {a}{r}}\right)^{n}C\_{n,0}P\_{n}(\cos \theta )}
### 拉普拉斯方法

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

#### 重力位的椭球谐级数

将正常重力位直接展开成椭球面 S\_{0} 上的级数，称为椭球谐级数，形式较球谐函数更为复杂：

{\displaystyle U=V+\Phi =\sum \_{n}^{\infty }{Q\_{n}\left(i{u \over E}\right) \over Q\_{n}\left(i{b \over E}\right)}A\_{n}P\_{n}(\sin \beta )+{1 \over 2}\omega ^{2}(u^{2}+E^{2})\cos ^{2}\beta }
上式中各项符号的物理或数学意义与球谐函数有所不同：

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

当 {\displaystyle u=b} 时，由重力位 {\displaystyle U\_{0}} 所决定的等位面应当与正常椭球面 S\_{0} 相重合，此时有：

{\displaystyle U\_{0}=V\_{0}+\Phi \_{0}=\sum \_{n}^{\infty }A\_{n}P\_{n}(\sin \beta )+{1 \over 2}\omega ^{2}a^{2}\cos ^{2}\beta }
#### 引力位的椭球谐级数

当且仅当所有含 {\displaystyle P\_{n}(\sin \beta )} 的项均为零时，对于任意 {\displaystyle \beta } 值该公式都成立。对其头三项进行展开：

{\displaystyle \left(A\_{0}+{1 \over 3}\omega ^{2}a^{2}-U\_{0}\right)P\_{0}(\sin \beta )+A\_{1}P\_{1}(\sin \beta )+\left(A\_{2}-{1 \over 3}\omega ^{2}a^{2}\right)P\_{2}(\sin \beta )+\sum \_{n=3}^{\infty }A\_{n}P\_{n}(\sin \beta )=0}
因此，有

* {\displaystyle A\_{0}=U\_{0}-{1 \over 3}\omega ^{2}a^{2}}
* {\displaystyle A\_{1}=0}
* {\displaystyle A\_{2}={1 \over 3}\omega ^{2}a^{2}}
* {\displaystyle A\_{3}=A\_{4}=\cdots =0}

得到引力位部分的椭球谐级数表达式为：

{\displaystyle V=\left(U\_{0}-{1 \over 3}\omega ^{2}a^{2}\right){Q\_{0}\left(i{u \over E}\right) \over Q\_{0}\left(i{b \over E}\right)}+{1 \over 3}\omega ^{2}a^{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 \over 3}\omega ^{2}a^{2}\right){\tan ^{-1}{E \over u} \over \tan ^{-1}{E \over b}}+{1 \over 3}\omega ^{2}a^{2}{q \over q\_{0}}P\_{2}(\sin \beta )}
其中

* {\displaystyle q={1 \over 2}\left[\left(1+3{\frac {u^{2}}{E^{2}}}\right)\tan ^{-1}{\frac {E}{u}}-3{\frac {u}{E}}\right]}
* {\displaystyle q\_{0}={1 \over 2}\left[\left(1+3{\frac {b^{2}}{E^{2}}}\right)\tan ^{-1}{\frac {E}{b}}-3{\frac {b}{E}}\right]}

#### 大地水准面重力位

考虑到椭圆坐标 {\displaystyle (u,\beta )} 与向径 {\displaystyle r} 存在如下转换关系：

{\displaystyle r^{2}=u^{2}+E^{2}\cos ^{2}\beta }
利用该关系对下式进行线性化，得：

{\displaystyle {\frac {1}{u}}={\frac {1}{r}}+\mathrm {O} ({\frac {1}{r^{3}}})}
{\displaystyle \tan ^{-1}{\frac {E}{u}}={\frac {E}{u}}+\mathrm {O} ({\frac {1}{u^{3}}})={\frac {E}{r}}+\mathrm {O} ({\frac {1}{r^{3}}})}
利用这两项关系式，可以得到线性化后的引力位函数：

{\displaystyle V=\left(U\_{0}-{1 \over 3}\omega ^{2}a^{2}\right){E \over \tan ^{-1}(E/b)}{\frac {1}{r}}+\mathrm {O} ({\frac {1}{r^{3}}})}
比较与 

{\displaystyle V={\frac {GM}{r}}+\mathrm {O} ({\frac {1}{r^{3}}})}
得出大地水准面重力位的表达式：

{\displaystyle U\_{0}={\frac {GM}{E}}\tan ^{-1}{\frac {E}{b}}+{1 \over 3}\omega ^{2}a^{2}}
#### 正常重力位计算公式

将大地水准面重力位的表达式代入原正常重力位的计算公式中，得：

{\displaystyle U(u,\beta )={\frac {GM}{E}}\tan ^{-1}{\frac {E}{u}}+{1 \over 2}\omega ^{2}a^{2}{\frac {q}{q\_{0}}}\left(\sin ^{2}\beta -{1 \over 3}\right)+{1 \over 2}\omega ^{2}(u^{2}+E^{2})\cos ^{2}\beta }

物理性质
----

### 球谐系数

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

#### 零阶项

当 n=0 时，球谐系数只有一项：

{\displaystyle A\_{0,0}=G\iiint \limits \_{\text{earth}}\operatorname {d} \!M'=GM}
或
{\displaystyle C\_{0,0}={1 \over M}\iiint \limits \_{\text{earth}}\operatorname {d} \!M'=1}
即球谐系数的零阶项反映了正常椭球体的地心引力常数或总质量。

#### 一阶项

当 n=1 时，球谐系数有三项：

{\displaystyle A\_{1,0}=G\iiint \limits \_{\text{earth}}z'\operatorname {d} \!M'=Gz\_{M}}
{\displaystyle A\_{1,1}=G\iiint \limits \_{\text{earth}}x'\operatorname {d} \!M'=Gx\_{M}}
{\displaystyle B\_{1,1}=G\iiint \limits \_{\text{earth}}y'\operatorname {d} \!M'=Gy\_{M}}
或
{\displaystyle C\_{1,0}={1 \over Ma}\iiint \limits \_{\text{earth}}z'\operatorname {d} \!M'={1 \over a}z\_{M}}
{\displaystyle C\_{1,1}={1 \over Ma}\iiint \limits \_{\text{earth}}x'\operatorname {d} \!M'={1 \over a}x\_{M}}
{\displaystyle S\_{1,1}={1 \over Ma}\iiint \limits \_{\text{earth}}y'\operatorname {d} \!M'={1 \over a}y\_{M}}
其中，{\displaystyle {\vec {r}}\_{M}=\left(x\_{M},y\_{M},z\_{M}\right)} 表示椭球的质心坐标。当坐标系的原点与椭球质心重合时，{\displaystyle x\_{M}=y\_{M}=z\_{M}=0}，所以球谐系数的一阶项一般都为零。

#### 二阶项

当 n=2 时，球谐系数有五项（仅以 {\displaystyle A\_{m,n}} 和 {\displaystyle B\_{m,n}} 为例）：

{\displaystyle A\_{2,0}=G\iiint \limits \_{\text{earth}}{1 \over 2}(2z'^{2}-x'^{2}-y'^{2})\operatorname {d} \!M'}
{\displaystyle A\_{2,1}=G\iiint \limits \_{\text{earth}}x'z'\operatorname {d} \!M'}
{\displaystyle B\_{2,1}=G\iiint \limits \_{\text{earth}}y'z'\operatorname {d} \!M'}
{\displaystyle A\_{2,2}=G\iiint \limits \_{\text{earth}}{1 \over 4}(x'^{2}-y'^{2})\operatorname {d} \!M'}
{\displaystyle B\_{2,2}=G\iiint \limits \_{\text{earth}}{1 \over 2}x'y'\operatorname {d} \!M'}
其中的五项积分，既可以四极矩张量 {\displaystyle \mathbf {Q} } 或惯性张量 {\displaystyle \mathbf {I} } 表达：

{\displaystyle \iiint \limits \_{\text{earth}}{1 \over 2}(2z'^{2}-x'^{2}-y'^{2})\operatorname {d} \!M'={1 \over 2}Q\_{zz}={1 \over 2}(I\_{xx}+I\_{yy})-I\_{zz}}
{\displaystyle \iiint \limits \_{\text{earth}}x'z'\operatorname {d} \!M'={1 \over 3}Q\_{xz}=-I\_{xz}}
{\displaystyle \iiint \limits \_{\text{earth}}y'z'\operatorname {d} \!M'={1 \over 3}Q\_{yz}=-I\_{yz}}
{\displaystyle \iiint \limits \_{\text{earth}}{1 \over 4}(x'^{2}-y'^{2})\operatorname {d} \!M'={1 \over 12}(Q\_{xx}-Q\_{yy})={1 \over 4}(I\_{yy}-I\_{xx})}
{\displaystyle \iiint \limits \_{\text{earth}}{1 \over 2}x'y'\operatorname {d} \!M'={1 \over 6}Q\_{xy}=-{1 \over 2}I\_{xy}}

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

参考文献
----

1. 1 2 3 4 5  宁津生. 管泽霖 , 编.. 测绘出版社. 1981.
2. 1 2 3 4 5 6 7 8 9 10 11 12 13 14  San Francisco W. H. Freeman and Company.. San Francisco: W. H. Freeman and Company. 1967 （英语）.  引文格式1维护：日期与年 (link)
3. 1 2  Jekeli, C.. . Elsevier. 2007: 11–42  [2020-04-15]. ISBN 978-0-444-52748-6. doi:10.1016/b978-044452748-6.00054-7. （原始内容存档于2018-07-01） （英语）.
4. 1 2 3 4 5 6  Sneeuw, Nico.  (PDF). Institute of Geodesy Universität Stuttgart. 2006  [2020-04-16]. （原始内容 (PDF)存档于2020-04-13）.

|  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |
| --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- |
| 

| * 
* 论
* 编
物理大地测量学 |
| --- |
|  |
| 
* **大地测量学：**应用大地测量学
* 椭球大地测量学
* **物理大地测量学**
* 空间大地测量学

 |
|  |
| 重要人物 | 
* 亚历克西斯·克劳德·克莱罗
* 皮埃尔-西蒙·拉普拉斯
* 約翰·彼得·古斯塔夫·勒熱納·狄利克雷
* 约翰·冯·诺伊曼
* 西莫恩·德尼·泊松
* 卡爾·弗里德里希·高斯
* 利斯廷
* 乔治·斯托克斯

* 弗里德里希·罗伯特·赫尔默特
* 海因里希·布伦斯
* 维宁·曼尼兹
* 米哈伊尔·莫洛坚斯基
* 阿恩·布耶哈马

 |
|  |
| 地球重力 | 

|  |  |
| --- | --- |
| 重力理论 | 
* 重力
	+ 正常重力
	+ 重力异常
	+ 重力扰动
* 重力位
	+ 正常重力位
	+ 扰动位
* 位基数

 |
|  |
| 地球形状 | 
* 参考椭球面
* 大地水准面
	+ 铅垂线
	+ 垂线偏差
	+ 大地水准面高
* 似地形表面
* 似大地水准面
* 海面地形
* 大地测量边值问题

 |
|  |
| 重力归算 | 
* 空间改正
* 布格改正
* 地形改正
* 地壳均衡改正

 |
|  |
| 相关定理 | 
* 斯托克斯定理
* 布隆斯公式
* 克莱罗定理
* 重力测量基本微分方程

 |

 |
|  |
| 测量技术 | 

|  |  |
| --- | --- |
| 测量方式 | 
* 水准测量
* 重力测量
* 卫星重力测量
	+ 卫星测高
	+ 卫星跟踪卫星
	+ 卫星重力梯度测量

 |
|  |
| 测量卫星 | 
* 小卫星挑战计划（CHAMP）
* 重力回溯及氣候實驗衛星（GRACE）
* 地球重力场和海洋环流探测卫星（GOCE）

 |

 |
|  |
| 高程系统 | 
* 力高
* 正高
* 正常高

 |
|  |
| 
* 分类 **大地测量学分类**
* 共享资源页面 **共享资源**
* 主题 **地球科学主题**

 |

 |