查看“︁扰动位”︁的源代码

'''扰动位'''（{{Lang-en|Disturbing potential}}），也称'''异常位'''（{{Lang-en|Anomalous potential}}），指地球的真实[[重力位]]与[[正常重力位]]之间的差异。<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|82}}

扰动位是建立[[地球]][[重力场模型]]过程中的关键变量，与[[大地水准面高]]和[[高程异常]]有着紧密的关系。<ref>{{Cite journal|title=Disturbing potential and its geometrical properties|url=http://content.sciendo.com/view/journals/sjce/18/3/article-p27.xml|last=Janák|first=J.|last2=Slovík|first2=S.|date=2010-09-01|journal=Slovak Journal of Civil Engineering|issue=3|doi=10.2478/v10189-010-0014-4|volume=18|pages=27–32|issn=1210-3896|last3=Fašková|first3=Z.|last4=Mikula|first4=K.|access-date=2020-04-14|archive-date=2018-06-23|archive-url=https://web.archive.org/web/20180623010514/https://content.sciendo.com/view/journals/sjce/18/3/article-p27.xml|dead-url=no}}</ref><ref name=":12">{{Cite book|chapter=|url=https://books.google.co.jp/books/about/Geodesy.html?id=pFO6VB_czRYC&redir_esc=y|publisher=Walter de Gruyter GmbH & Co KG|date=|isbn=978-3-11-017072-6|language=en|first=Wolfgang|last=Torge|title=Geodesy|year=2001|location=|pages=}}</ref>{{Rp|214}}在求解地球形状和地球重力位的问题的过程中，可以先定义一个简单的、能够直接计算得到的[[正常重力位]]和[[正常椭球体]]，再通过求解扰动位得到[[大地水准面]]或[[似大地水准面]]与正常椭球体之间的差距（如大地水准面高和[[垂线偏差]]），从而得到地球的近似形状和真实重力位。'''<ref name="whu2">{{cite book|author1=孔祥元|author2=郭际明|author3=刘宗泉|title=大地测量学基础|publisher=武汉大学出版社|ISBN=978-7-30-707562-7|pages=|last=|first=|year=2001|isbn=|location=}}</ref>{{Rp|20}}'''在选取正常椭球体时，通常定义其与大地水准面密合，扰动位的量级很小（仅占真实重力位的百万分之五<ref>{{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=":0">{{Cite book|title=地球形状及外部重力场|author=宁津生|publisher=测绘出版社|year=1981|isbn=|location=|authorlink=宁津生|editor=管泽霖|first=}}</ref>{{Rp|243}}，通常可以用[[线性近似]]和球面近似的方法进行求解<ref name=":12" />{{Rp|64}}。

== 数学表达 ==
从数学上的定义来看，扰动位 <math>T</math> 通常表达成真实重力位 <math>W</math> 与[[正常重力位]] <math>U</math> 之间的差距：<ref name=":2" />{{Rp|82}}

: <math>
T = W - U
</math>

其中，两个重力位都由[[引力位]]部分和离心力位的部分组成，且两者的离心力位部分可以视作是相同的<ref name=":12" />{{Rp|214}}，因此扰动位表现的是两者的引力位差，具有引力位满足[[拉普拉斯方程]]的性质：<ref name=":2" />{{Rp|86}}

: <math>
\nabla T \equiv 0
</math>

=== 球谐函数 ===
因此，扰动位 <math>T</math> 在边界面的外部（<math>
r>R
</math>）展开为[[球谐函数]]：<ref name=":2" />{{Rp|88}} 

:<math>
T(r,\theta,\lambda) = \sum_{n=0}^{\infin}\left({R \over r}\right)^{n+1}T_n(\theta,\lambda)
</math>

其中[[球坐标]] <math>
(r,\theta,\lambda)
</math> 表示外部空间某一点的径向距离、极角和经度。<math>
T_n(\theta,\lambda)
</math> 则表示 <math>
n
</math> 阶[[面谐函数]]，且有完全形式：<ref name=":12" />{{Rp|215}} 

:<math>
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})
</math>

上式中各项的含义如下： 

* <math>
GM
</math> [[標準重力參數|地心引力常数]]
* <math>
n
</math>为面谐函数的阶数

* <math>R</math> 为用以近似的地球半径
* <math>P_{nm}</math> 是 <math>n\ </math>阶 <math>m\ </math>次的[[缔合勒让德多项式]]
* <math>\Delta C_{nm}</math> 和 <math>\Delta S_{nm}</math> 是真实重力位与正常重力位所采用的球谐系数的差值

又根据[[正常椭球体]]的定义，其产生的[[正常重力位]]应当与真实重力位在球谐展开式的最大项上相同，且具有一阶系数为零的性质，所以扰动位也常写为：<ref name=":12" />{{Rp|215}}

:<math>
T(r,\theta,\lambda) = \sum_{n=2}^{\infin}\left({R \over r}\right)^{n+1}T_n(\theta,\lambda)
</math>

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

:<math>
{\partial T\over\partial r} = -{1 \over r}\sum_{n=2}^{\infin}(n+1)\left({R \over r}\right)^{n+1}T_n(\theta,\lambda)
</math>
:<math>
{\partial^2 T\over\partial r^2} 
= {1 \over r^2}\sum_{n=2}^{\infin}(n+1)(n+2)\left({R \over r}\right)^{n+1}T_n(\theta,\lambda)
</math>

利用缩写 <math>
T'_n(\theta,\lambda) = (n+1)(n+2)T_n(\theta,\lambda)
</math>，可以得到[[调和函数]] <math>
r^2{\partial^2 T\over\partial r^2}
</math>的球谐展开式： 

:<math>
r^2{\partial^2 T\over\partial r^2} = \sum_{n=2}^{\infin}\left({R \over r}\right)^{n+1}T'_n(\theta,\lambda)
</math>

== 球面近似 ==
在球近似的条件下，以下三个[[偏微分]]被视作相同：<ref name=":2" />{{Rp|87}}

:<math>
{\partial\over\partial n} = {\partial\over\partial h} = {\partial\over\partial r}
</math>

<math>
n
</math>、<math>
h
</math> 和 <math>
r
</math> 分别表示重力[[矢量]]的方向、高程方向和地心方向。 

=== 大地水准面上的扰动位及其径向导数 ===
以[[球面]]对[[大地水准面]]进行近似（即假设 <math>
r = R
</math>），不考虑球谐函数是否收敛的问题，则大地水准面上的扰动位可以表达为：<ref name=":12" />{{Rp|215}}

:<math>
T = T(R,\theta,\lambda) = \sum_{n=2}^{\infin}T_n(\theta,\lambda)
</math>

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

=== 重力扰动 ===
{{Main|重力扰动}}
重力扰动是指大地水准面上的一点 <math>\mathbf{P}</math> 处，真实重力 <math>\vec{\text{g}}_P</math> 与同一位置上的正常重力 <math>\vec{\gamma}_P</math> 的差异<ref name=":12" />{{Rp|84}}，即

:<math>\delta \vec{\text{g}} = \vec{\text{g}}_P - \vec{\gamma}_P  </math>

利用球面近似，重力扰动可以通过扰动位的一阶径向导数来表述：<ref name=":12" />{{Rp|85}}

:<math>\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}  </math>

其中 <math>
n
</math> 和 <math>
n'
</math> 分别表示真实重力与正常重力方向，即[[铅垂线]]方向和椭球的[[法线]]方向。

将其展开作[[级数]]，得：<ref name=":12" />{{Rp|88}}

:<math>\delta \text{g} = {1 \over r}\sum_{n=2}^{\infin}(n+1)\left({R \over r}\right)^{n+1}T_n(\theta,\lambda)  </math>

=== 重力异常 ===
{{Main|重力异常}}
重力异常与重力扰动的区别在于，重力异常比较的重力是点 <math>\mathbf{P}</math> 处的真实重力 <math>\vec{\text{g}}_P</math> 和其在椭球面上的投影 <math>\mathbf{Q}</math> 处的正常重力 <math>\vec{\gamma}_Q</math> <ref name=":12" />{{Rp|83}}，即

:<math>\Delta \vec{\text{g}} = \vec{\text{g}}_P - \vec{\gamma}_Q  </math>

其与重力扰动和扰动位的一阶径向导数的关系由[[重力测量基本微分方程]]给出<ref name=":12" />{{Rp|86}}，这一方程的球面近似形式为：

:<math>{\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  </math>

将其展开作[[级数]]，得：<ref name=":12" />{{Rp|89}}

:<math>\Delta\text{g} = {1 \over r}\sum_{n=2}^{\infin}(n-1)\left({R \over r}\right)^{n+1}T_n(\theta,\lambda)  </math>

=== 斯托克斯公式 ===
{{Main|斯托克斯公式（大地测量学）}}
通过[[泊松积分]]公式，可以在表达出大地水准面上的重力异常 <math>\Delta\text{g}  </math> 和扰动位 <math>T  </math> 之间的关系，即：<ref name=":12" />{{Rp|93-94}}

:<math>
T = {R\over4\pi}\iint\limits_{\sigma}\Delta g \,S(\psi)\operatorname{d}\!\sigma 
</math>

这一公式由爱尔兰数学物理学家[[斯托克斯]]在1849年给出，又称为斯托克斯公式。其中的 <math>
S(\psi)
</math> 被称为[[斯托克斯函数]]<ref>{{Cite book|chapter=|url=https://books.google.co.jp/books?id=yG4iAQAAIAAJ&printsec=frontcover&hl=zh-CN#v=onepage&q&f=false|publisher=U.S. Government Printing Office|date=1936|language=en|first=U. S. Coast and Geodetic|last=Survey|first2=Walter Davis|last2=Lambert|first3=Frederic Warren|last3=Darling|title=Tables for Determining the Form of the Geoid and Its Indirect Effect on Gravity|year=|isbn=|location=|pages=|access-date=2020-04-15|archive-date=2020-08-21|archive-url=https://web.archive.org/web/20200821034906/https://books.google.co.jp/books?id=yG4iAQAAIAAJ&printsec=frontcover&hl=zh-CN#v=onepage&q&f=false|dead-url=no}}</ref>，该项由单位球面上的被计算点与重力异常观测值所在的角元素之间的夹角 <math>
\psi
</math> 决定：<ref name=":12" />{{Rp|94}}

: <math>
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})}
</math>

== 相关条目 ==

* [[重力测量]]

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

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