**扰动位**（英語：），也称**异常位**（英語：），指地球的真实重力位与正常重力位之间的差异。

扰动位是建立地球重力场模型过程中的关键变量，与大地水准面高和高程异常有着紧密的关系。在求解地球形状和地球重力位的问题的过程中，可以先定义一个简单的、能够直接计算得到的正常重力位和正常椭球体，再通过求解扰动位得到大地水准面或似大地水准面与正常椭球体之间的差距（如大地水准面高和垂线偏差），从而得到地球的近似形状和真实重力位。在选取正常椭球体时，通常定义其与大地水准面密合，扰动位的量级很小（仅占真实重力位的百万分之五），对真实重力位起到改正项的作用，通常可以用线性近似和球面近似的方法进行求解。

数学表达
----

从数学上的定义来看，扰动位 T 通常表达成真实重力位 W 与正常重力位 U 之间的差距：

{\displaystyle T=W-U}
其中，两个重力位都由引力位部分和离心力位的部分组成，且两者的离心力位部分可以视作是相同的，因此扰动位表现的是两者的引力位差，具有引力位满足拉普拉斯方程的性质：

{\displaystyle \nabla T\equiv 0}
### 球谐函数

因此，扰动位 T 在边界面的外部（{\displaystyle r>R}）展开为球谐函数：

{\displaystyle T(r,\theta ,\lambda )=\sum \_{n=0}^{\infty }\left({R \over r}\right)^{n+1}T\_{n}(\theta ,\lambda )}
其中球坐标 {\displaystyle (r,\theta ,\lambda )} 表示外部空间某一点的径向距离、极角和经度。{\displaystyle T\_{n}(\theta ,\lambda )} 则表示 {\displaystyle n} 阶面谐函数，且有完全形式：

{\displaystyle T\_{n}(\theta ,\lambda )={GM \over R}\sum \_{m=0}^{n}\left(\Delta C\_{nm}\cos {m\lambda }+\Delta S\_{nm}\sin {m\lambda }\right)P\_{nm}(\cos {\theta })}
上式中各项的含义如下： 

* {\displaystyle GM} 地心引力常数
* {\displaystyle n}为面谐函数的阶数

* R 为用以近似的地球半径
* {\displaystyle P\_{nm}} 是 **解析失败 (SVG（MathML可通过浏览器插件启用）：从服务器"http://localhost:6011/zh.wikipedia.org/v1/"返回无效的响应（"Math extension cannot connect to Restbase."）：): n\** 
阶 m\ 次的缔合勒让德多项式
* {\displaystyle \Delta C\_{nm}} 和 {\displaystyle \Delta S\_{nm}} 是真实重力位与正常重力位所采用的球谐系数的差值

又根据正常椭球体的定义，其产生的正常重力位应当与真实重力位在球谐展开式的最大项上相同，且具有一阶系数为零的性质，所以扰动位也常写为：

{\displaystyle T(r,\theta ,\lambda )=\sum \_{n=2}^{\infty }\left({R \over r}\right)^{n+1}T\_{n}(\theta ,\lambda )}

径向导数
----

由扰动位的球谐表达式，可以求出其一阶和二阶径向导数的相应表达式：

{\displaystyle {\partial T \over \partial r}=-{1 \over r}\sum \_{n=2}^{\infty }(n+1)\left({R \over r}\right)^{n+1}T\_{n}(\theta ,\lambda )}
{\displaystyle {\partial ^{2}T \over \partial r^{2}}={1 \over r^{2}}\sum \_{n=2}^{\infty }(n+1)(n+2)\left({R \over r}\right)^{n+1}T\_{n}(\theta ,\lambda )}
利用缩写 {\displaystyle T'\_{n}(\theta ,\lambda )=(n+1)(n+2)T\_{n}(\theta ,\lambda )}，可以得到调和函数 {\displaystyle r^{2}{\partial ^{2}T \over \partial r^{2}}}的球谐展开式： 

{\displaystyle r^{2}{\partial ^{2}T \over \partial r^{2}}=\sum \_{n=2}^{\infty }\left({R \over r}\right)^{n+1}T'\_{n}(\theta ,\lambda )}

球面近似
----

在球近似的条件下，以下三个偏微分被视作相同：

{\displaystyle {\partial \over \partial n}={\partial \over \partial h}={\partial \over \partial r}}
{\displaystyle n}、{\displaystyle h} 和 {\displaystyle r} 分别表示重力矢量的方向、高程方向和地心方向。 

### 大地水准面上的扰动位及其径向导数

以球面对大地水准面进行近似（即假设 {\displaystyle r=R}），不考虑球谐函数是否收敛的问题，则大地水准面上的扰动位可以表达为：

{\displaystyle T=T(R,\theta ,\lambda )=\sum \_{n=2}^{\infty }T\_{n}(\theta ,\lambda )}

与其他物理量的关系
---------

### 重力扰动

条目：重力扰动
重力扰动是指大地水准面上的一点 \mathbf {P}  处，真实重力 {\displaystyle {\vec {\text{g}}}\_{P}} 与同一位置上的正常重力 {\displaystyle {\vec {\gamma }}\_{P}} 的差异，即

{\displaystyle \delta {\vec {\text{g}}}={\vec {\text{g}}}\_{P}-{\vec {\gamma }}\_{P}}
利用球面近似，重力扰动可以通过扰动位的一阶径向导数来表述：

{\displaystyle \delta {\text{g}}=-({\partial W \over \partial n}-{\partial U \over \partial n'})\longrightarrow \delta {\text{g}}=-({\partial W \over \partial n}-{\partial U \over \partial n})=-{\partial T \over \partial n}\longrightarrow \delta {\text{g}}=-{\partial T \over \partial r}}
其中 {\displaystyle n} 和 {\displaystyle n'} 分别表示真实重力与正常重力方向，即铅垂线方向和椭球的法线方向。

将其展开作级数，得：

{\displaystyle \delta {\text{g}}={1 \over r}\sum \_{n=2}^{\infty }(n+1)\left({R \over r}\right)^{n+1}T\_{n}(\theta ,\lambda )}
### 重力异常

条目：重力异常
重力异常与重力扰动的区别在于，重力异常比较的重力是点 \mathbf {P}  处的真实重力 {\displaystyle {\vec {\text{g}}}\_{P}} 和其在椭球面上的投影 {\mathbf {Q}} 处的正常重力 {\displaystyle {\vec {\gamma }}\_{Q}} ，即

{\displaystyle \Delta {\vec {\text{g}}}={\vec {\text{g}}}\_{P}-{\vec {\gamma }}\_{Q}}
其与重力扰动和扰动位的一阶径向导数的关系由重力测量基本微分方程给出，这一方程的球面近似形式为：

{\displaystyle {\partial T \over \partial r}+{2 \over r}T+\Delta {\text{g}}=0\longrightarrow \Delta {\text{g}}=-{\partial T \over \partial r}-{2 \over r}T}
将其展开作级数，得：

{\displaystyle \Delta {\text{g}}={1 \over r}\sum \_{n=2}^{\infty }(n-1)\left({R \over r}\right)^{n+1}T\_{n}(\theta ,\lambda )}
### 斯托克斯公式

条目：斯托克斯公式（大地测量学）
通过泊松积分公式，可以在表达出大地水准面上的重力异常 {\displaystyle \Delta {\text{g}}} 和扰动位 {\displaystyle T} 之间的关系，即：

{\displaystyle T={R \over 4\pi }\iint \limits \_{\sigma }\Delta g\,S(\psi )\operatorname {d} \!\sigma }
这一公式由爱尔兰数学物理学家斯托克斯在1849年给出，又称为斯托克斯公式。其中的 {\displaystyle S(\psi )} 被称为斯托克斯函数，该项由单位球面上的被计算点与重力异常观测值所在的角元素之间的夹角 {\displaystyle \psi } 决定：

{\displaystyle S(\psi )={1 \over \sin(\psi /2)}-6\sin {\psi \over 2}+1-5\cos {\psi }-3\cos {\psi }\ln {(\sin {\psi \over 2}+\sin ^{2}{\psi \over 2})}}

相关条目
----

* 重力测量

* 重力归算
* 大地测量边值问题

参考文献
----

1. 1 2 3 4 5  San Francisco W. H. Freeman and Company.. San Francisco: W. H. Freeman and Company. 1967 （英语）.  引文格式1维护：日期与年 (link)
2. ↑  Janák, J.; Slovík, S.; Fašková, Z.; Mikula, K.. Slovak Journal of Civil Engineering. 2010-09-01, **18** (3): 27–32  [2020-04-14]. ISSN 1210-3896. doi:10.2478/v10189-010-0014-4. （原始内容存档于2018-06-23）.
3. 1 2 3 4 5 6 7 8 9 10 11 12 13 14  Torge, Wolfgang.. Walter de Gruyter GmbH & Co KG. 2001. ISBN 978-3-11-017072-6 （英语）.
4. ↑  孔祥元; 郭际明; 刘宗泉.. 武汉大学出版社. 2001. ISBN 978-7-30-707562-7.
5. ↑  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） （英语）.
6. ↑  宁津生. 管泽霖 , 编.. 测绘出版社. 1981.
7. ↑  Survey, U. S. Coast and Geodetic; Lambert, Walter Davis; Darling, Frederic Warren.. U.S. Government Printing Office. 1936  [2020-04-15]. （原始内容存档于2020-08-21） （英语）.

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

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

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

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

 |
|  |
| 地球重力 | 

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

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

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

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

 |

 |
|  |
| 测量技术 | 

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

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

 |

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

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

 |

 |