重力学：第七章 固体潮

n7.0 引言引言 n7.1 起潮力位及其调和展开起潮力位及其调和展开n7.2 刚体地球的理论固体潮刚体地球的理论固体潮n7.3 理论重力固体潮计算理论重力固体潮计算 在月亮和太阳的作用下，海水每天两次的周期性涨落称为海潮。海在月亮和太阳的作用下，海水每天两次的周期性涨落称为海潮。海潮现象非常明显，极易察觉。十九世纪末，英国天文学家潮现象非常明显，极易察觉。十九世纪末，英国天文学家 乔治乔治达尔文达尔文(George Howard Darwin)分析了当时积累的海潮观测资料，发现接近平分析了当时积累的海潮观测资料，发现接近平衡潮的月球半月潮的实际潮高比把地球看成刚体时的理论潮高小三分之衡潮的月球半月潮的实际潮高比把地球看成刚体时的理论潮高小三分之一。为了解释这种现象，他认为地球的固体表面也发生与海水类似的周一。为了解释这种现象，他认为地球的固体表面也发生与海水类似的周期性涨落，其涨落幅度约为海水涨落幅度的三分之一。后来，就把期性涨落，其涨落幅度约为海水涨落幅度的三分之一。后来，就把地球地球整体在月球和太阳引力的作用下的变形称为固体潮整体在月球和太阳引力的作用下的变形称为固体潮(Solid tide)。7.0 引言引言 达尔文于达尔文于1883年，从理论上对简化的地球模型在月亮和太阳作用下年，从理论上对简化的地球模型在月亮和太阳作用下的潮汐变形进行了计算，计算结果表明，若把地球看成均匀的、不可压的潮汐变形进行了计算，计算结果表明，若把地球看成均匀的、不可压缩的刚性球体，那么，缩的刚性球体，那么，要使地球表面上月亮和太阳的作用下的涨落幅度要使地球表面上月亮和太阳的作用下的涨落幅度为海水涨落幅度的三分之一，地球的切变模量约等于钢的切变模量为海水涨落幅度的三分之一，地球的切变模量约等于钢的切变模量。这。这是有关地球整体表现为固体的最早结论。是有关地球整体表现为固体的最早结论。地球固体潮伴随一系列地球物理地球固体潮伴随一系列地球物理现象现象：1相对地球表面的海潮；相对地球表面的海潮；2重力变化重力变化重力固体潮重力固体潮(gravity tide);3.地面倾斜地面倾斜地倾斜固体潮地倾斜固体潮;4.地面的变形地面的变形应变固体潮；应变固体潮；5经纬度变化经纬度变化经纬度固体潮；经纬度固体潮；6地球自转速度的变化等。地球自转速度的变化等。7.0 引言引言地球的固体潮地球的固体潮 固体潮（固体潮（Earth tide,Solid tide,Body tide）在太阳和月球引力的作在太阳和月球引力的作用下，固体地球产生的周期形变的现象。月球和太阳对地球的引力不但可用下，固体地球产生的周期形变的现象。月球和太阳对地球的引力不但可以引起地球表面流体的潮汐（如海潮、大气潮），还能引起地球固体部分以引起地球表面流体的潮汐（如海潮、大气潮），还能引起地球固体部分的周期性形变。的周期性形变。重力潮汐变化影响的最大幅度可达重力潮汐变化影响的最大幅度可达130 130 微伽，重力测量结果和精密微伽，重力测量结果和精密大地测量结果中应加入相应的修正。大地测量结果中应加入相应的修正。7.0 引言引言地球的固体潮地球的固体潮 地球固体潮是如何形成的？地球固体潮是如何形成的？起潮力。起潮力。以月球为例，由于月球是绕月以月球为例，由于月球是绕月-地共同质心地共同质心O 旋转，地球上各质点受到的离旋转，地球上各质点受到的离心力与月球引力的合力就是心力与月球引力的合力就是起潮力起潮力。月球和太阳相对地球位置不同时，地。月球和太阳相对地球位置不同时，地球上各质点受到的力的大小和方向都不同。球上各质点受到的力的大小和方向都不同。7.0 引言引言地球的固体潮地球的固体潮 起潮力主要来源？起潮力主要来源？太阳的质量虽然比月球的质量大，但月球同地球的距离比太阳同地球太阳的质量虽然比月球的质量大，但月球同地球的距离比太阳同地球的距离近，月球的引潮力比太阳的引潮力大。的距离近，月球的引潮力比太阳的引潮力大。日、月对地球的引力比较：日、月对地球的引力比较：日、月对地球（地面与球心）引力差比较：日、月对地球（地面与球心）引力差比较：gr=gAr gOr g=gA gO gmr/gsr2.18所以，所以，月球的引力是产生在固体潮主要原因。月球的引力是产生在固体潮主要原因。100/100.2/105.3/202222msmmmsFFNrGMmFNrGMmF7.0 引言引言地球的固体潮地球的固体潮 由于其他天体距地球甚远，对地球的引力甚微，在固体潮的研究中一般由于其他天体距地球甚远，对地球的引力甚微，在固体潮的研究中一般可略而不计。引潮力是作用在地球的单位质点上的太阳、月引力和地球绕可略而不计。引潮力是作用在地球的单位质点上的太阳、月引力和地球绕地月（和地日）公共质心旋转所产生的惯性离心力的合力。随着作用点的地月（和地日）公共质心旋转所产生的惯性离心力的合力。随着作用点的位置不同，引潮力的大小、方向也发生改变。位置不同，引潮力的大小、方向也发生改变。固体潮的重力响应固体潮的重力响应重力重力固体潮在一个固定点上的特征固体潮在一个固定点上的特征如图所示。图中可见。重力固如图所示。图中可见。重力固体潮变化由不同周期信叠加而体潮变化由不同周期信叠加而成。成。7.0 引言引言海潮海潮 地面上观测到的固体潮除了包括地球对起潮力的直接响应外，还包括地地面上观测到的固体潮除了包括地球对起潮力的直接响应外，还包括地球对海潮和大气潮的响应。球对海潮和大气潮的响应。地球对海潮和大气潮的响应称为地球的负荷潮地球对海潮和大气潮的响应称为地球的负荷潮，因为地球的固体潮、海潮和大气潮都来源于起潮力，因而地球的固体潮及因为地球的固体潮、海潮和大气潮都来源于起潮力，因而地球的固体潮及其负荷潮的谱线完全相同。其负荷潮的谱线完全相同。海水受引潮力作用而产生的海洋水体的长周期波动现象。它在铅直方向海水受引潮力作用而产生的海洋水体的长周期波动现象。它在铅直方向表现为潮位升降，在水平方向表现为潮流涨落。随着固体潮观测仪器精度表现为潮位升降，在水平方向表现为潮流涨落。随着固体潮观测仪器精度的提高以及固体潮观测资料的积累，证明了地球对海潮负荷产生影响，地的提高以及固体潮观测资料的积累，证明了地球对海潮负荷产生影响，地球对海潮负荷的这种响应后来就称为球对海潮负荷的这种响应后来就称为海潮负荷潮海潮负荷潮。在沿海地区海潮负荷潮。在沿海地区海潮负荷潮非常显著，负荷潮的幅度可能接近或者超过固体潮本身，在一般情况下在非常显著，负荷潮的幅度可能接近或者超过固体潮本身，在一般情况下在沿海地区重力负荷潮约占重力固体潮的沿海地区重力负荷潮约占重力固体潮的10%，应变负荷潮约占应变固体潮，应变负荷潮约占应变固体潮的的25%，地倾斜负荷潮约占地倾斜固体潮的，地倾斜负荷潮约占地倾斜固体潮的90%。7.0 引言引言大气潮大气潮 大气潮与海潮产生的原理相同，发生的规律也相似，在地球上某一点大气潮与海潮产生的原理相同，发生的规律也相似，在地球上某一点的海潮每天产生两次涨潮两次落潮的现象，大气潮同样会出现两次涨潮和的海潮每天产生两次涨潮两次落潮的现象，大气潮同样会出现两次涨潮和两次落潮，而且两次涨潮所经过的时间平均是两次落潮，而且两次涨潮所经过的时间平均是12小时小时25分，第二天涨潮的分，第二天涨潮的时间会比前一天平均推迟时间会比前一天平均推迟50分钟。在很多方面大气潮和海洋潮汐类似。大分钟。在很多方面大气潮和海洋潮汐类似。大气潮的激发机制包括：气潮的激发机制包括：大气辐射加热的日夜更替、月球的引力场影响、行大气辐射加热的日夜更替、月球的引力场影响、行星波和大气潮之间的非线性相互作用。简单来说月球和太阳对地球大气的星波和大气潮之间的非线性相互作用。简单来说月球和太阳对地球大气的摄引摄引（引力的波动）也会产生大气的（引力的波动）也会产生大气的“涨潮涨潮”和和“落潮落潮”，这就是大气潮。，这就是大气潮。大气潮是全球尺度的大气振荡。大气潮是全球尺度的大气振荡。根据万有引力定律，两个物体之间的引力和它们的质量成正比，因此根据万有引力定律，两个物体之间的引力和它们的质量成正比，因此对于潮汐来说，起潮力与被摄引物体的质量成正比。与海水相比大气质量对于潮汐来说，起潮力与被摄引物体的质量成正比。与海水相比大气质量（密度）小得多，所以大气潮引起的引力变化远不如海潮显著，只有用极（密度）小得多，所以大气潮引起的引力变化远不如海潮显著，只有用极精密的仪器测量才能发现。精密的仪器测量才能发现。7.0 引言引言7.0 引言引言 重力固体潮观测的主要目的：重力固体潮观测的主要目的：为重力测量和其他的精密测量（卫星定位和精密水准测量等）为重力测量和其他的精密测量（卫星定位和精密水准测量等）提供固体潮改正。提供固体潮改正。在消除模拟的海洋潮汐影响后，固体潮观测可用于研究全球地在消除模拟的海洋潮汐影响后，固体潮观测可用于研究全球地球模型。（研究固体地核平动振荡，液核共振效应，核幔边界黏性球模型。（研究固体地核平动振荡，液核共振效应，核幔边界黏性等）等）在扣除模拟的海洋潮汐影响和全球性固体潮影响后，可以用来在扣除模拟的海洋潮汐影响和全球性固体潮影响后，可以用来确定区域性的确定区域性的“地壳地壳-地幔模型地幔模型”，揭示岩石圈和上地幔中的横向不，揭示岩石圈和上地幔中的横向不均匀性。均匀性。检验全球和区域性的检验全球和区域性的“海洋潮汐模型海洋潮汐模型”。研究固体地球潮汐和海洋负荷潮汐与现代地球动力学过程（地研究固体地球潮汐和海洋负荷潮汐与现代地球动力学过程（地震、火山爆发）之间的联系。震、火山爆发）之间的联系。7.0 引言引言给定的给定的地球模型地球模型理论的理论的固体潮值固体潮值理论计算理论计算实际的实际的固体潮值固体潮值比较对比比较对比与实际接近的与实际接近的地球模型地球模型选择选择 n7.0 引言引言 n7.1 起潮力位及其调和展开起潮力位及其调和展开n7.2 刚体地球的理论固体潮刚体地球的理论固体潮n7.3 理论重力固体潮计算理论重力固体潮计算 7.1.1 起潮力起潮力 (1)地月系统的质心地月系统的质心 设地球质量为设地球质量为M，质心为，质心为Oe，月球的质量为，月球的质量为m，质心为，质心为Om，地月系统，地月系统质心为质心为O，地月质心之间距离为，地月质心之间距离为r，则根据合力矩定理，有，则根据合力矩定理，有 7.1 起潮力位及其调和展开起潮力位及其调和展开OeOOmrRMRRrmMmOOmrOOmMee73.0)3.60(15.811)(即 7.1.1 起潮力起潮力 (1)地月系统的质心地月系统的质心 根据万有引力定律，地球与月球的相互引力为根据万有引力定律，地球与月球的相互引力为 7.1 起潮力位及其调和展开起潮力位及其调和展开OeOOmrRM2rmMGFme 7.1.1 起潮力起潮力 (2)地月系统运动轨迹地月系统运动轨迹 相对地月质心运动轨迹（红色虚线）而言，地球质心的运动轨迹在基相对地月质心运动轨迹（红色虚线）而言，地球质心的运动轨迹在基本轨道上摆动。本轨道上摆动。7.1 起潮力位及其调和展开起潮力位及其调和展开OeOmO 7.1.1 起潮力起潮力 (3)地月系统围绕共地月系统围绕共同质心同质心O点旋转的惯性点旋转的惯性离心力离心力 以地月系统质心为参以地月系统质心为参考，看地、月的运动轨考，看地、月的运动轨迹。迹。7.1 起潮力位及其调和展开起潮力位及其调和展开 7.1.1 起潮力起潮力 若以地心为参考，若以地心为参考，看地、月运动轨迹。实看地、月运动轨迹。实际上，由于地球自转，际上，由于地球自转，地月系统质心地月系统质心O划过地划过地球内部球内部0.73R部位的点，部位的点，形成圆弧。形成圆弧。7.1 起潮力位及其调和展开起潮力位及其调和展开 7.1.1 起潮力起潮力 若以地月系统质心若以地月系统质心O为参考，在不考虑地球为参考，在不考虑地球自转的情况下（地球自转与潮汐没有直接关系），自转的情况下（地球自转与潮汐没有直接关系），考虑地球上考虑地球上任一点任一点A，它随地月系统绕，它随地月系统绕O点的运动点的运动轨迹为轨迹为AA圆弧，它与地心运动轨迹圆弧，它与地心运动轨迹OeOe 为相同为相同的圆弧。这些圆弧属圆心不同而半径相等的圆周，的圆弧。这些圆弧属圆心不同而半径相等的圆周，其半径等于其半径等于0.73R。7.1 起潮力位及其调和展开起潮力位及其调和展开 7.1.1 起潮力起潮力 月球的质心月球的质心Om沿以地月系统沿以地月系统的质心的质心O为焦点的椭圆旋转。为为焦点的椭圆旋转。为了保持月地系统的动态平衡，地了保持月地系统的动态平衡，地心心Oe绕绕O旋转，使得旋转，使得Oe、O、Om始终保持在一条直线上，始终保持在一条直线上，O在空在空间的位置不动。间的位置不动。7.1 起潮力位及其调和展开起潮力位及其调和展开 7.1.1 起潮力起潮力 此时，地球所做的运动就是平动。此时，地球所做的运动就是平动。理论力学中，平动的概念就是指物体内理论力学中，平动的概念就是指物体内部各点都具有相同的速度和加速度。在部各点都具有相同的速度和加速度。在这里，各点平动是圆周运动，即这里，各点平动是圆周运动，即A点与点与Oe点具有相同的向心加速度点具有相同的向心加速度。可以证明，。可以证明，地球内部任何一点都与地球内部任何一点都与A点一样，点一样，具有具有相同的向心加速度相同的向心加速度。显然，这里包括地球球心处。显然，这里包括地球球心处。7.1 起潮力位及其调和展开起潮力位及其调和展开OeOOmrO eO mA A 7.1.1 起潮力起潮力 这样，地球内任一点受到的这样，地球内任一点受到的惯性离心力相同，即等于地球惯性离心力相同，即等于地球质心处的离心力，且方向始终质心处的离心力，且方向始终平行于质心连线，即有平行于质心连线，即有 7.1 起潮力位及其调和展开起潮力位及其调和展开rrOOqemem2qm 7.1.1 起潮力起潮力 (3)地月系统围绕共同质心地月系统围绕共同质心O点旋转点旋转的惯性离心力的惯性离心力 可见，月球对地球提供的引力维可见，月球对地球提供的引力维持着地球绕持着地球绕O点作圆周运动，点作圆周运动，Oe点即点即受到月球的引力就等于该点的惯性离受到月球的引力就等于该点的惯性离心力，即心力，即 7.1 起潮力位及其调和展开起潮力位及其调和展开22rmGOOqemme 7.1.1 起潮力起潮力 但地球内各质点受但地球内各质点受到月球的引力在大小、到月球的引力在大小、方向上都不一样。方向上都不一样。7.1 起潮力位及其调和展开起潮力位及其调和展开 7.1.1 起潮力起潮力 (4)起潮力起潮力 起潮力是随着时间变化的，当起潮力是随着时间变化的，当A点处于不同位置时，起潮力大点处于不同位置时，起潮力大小不同，形成周期性变化。小不同，形成周期性变化。7.1 起潮力位及其调和展开起潮力位及其调和展开smsmmmmmmttttqfttqf总起潮力为地球上任一点所受到的。起潮力同理，可以导出太阳的为月球起潮力，即之矢量合与惯性离心力月球引力OeOOmrAr qmfmtmzm 7.1.2 起潮力位起潮力位 (1)月球起潮力位月球起潮力位 由于地球的自转，月球绕地球的旋转和地球绕太阳的旋转，因而在地由于地球的自转，月球绕地球的旋转和地球绕太阳的旋转，因而在地球内部任一点的起潮力是一个随时间变化的矢量场。球内部任一点的起潮力是一个随时间变化的矢量场。物体的引力是保守物体的引力是保守力，因而月球在力，因而月球在A点的引力点的引力fm(A)可用位函数表示可用位函数表示，即即Um(A)。地球内部的惯性离心力地球内部的惯性离心力qm(A)也是一个保守力，也可用地球在月球引力也是一个保守力，也可用地球在月球引力场内的平动加速运动在场内的平动加速运动在A点产生的惯性离心力位点产生的惯性离心力位Qm(A)的梯度表示，即的梯度表示，即 qm(A)gradQm(A)用用Tm(A)表示引力位表示引力位Um(A)和惯性离心力位和惯性离心力位Qm(A)的和，即的和，即 Tm(A)Um(A)+Qm(A)称为月球在地球内部任一点所产生的起潮力位。称为月球在地球内部任一点所产生的起潮力位。显然，有显然，有 tm(A)Tm(A)7.1 起潮力位及其调和展开起潮力位及其调和展开 7.1.2 起潮力位起潮力位 (1)月球起潮力位月球起潮力位 设地球、月球和太阳的质量分别为设地球、月球和太阳的质量分别为Me,Mm,Ms,，地月距、地日距分别为，地月距、地日距分别为rm,rs，r为地球内部为地球内部A点半径，则月球起潮力位可写成点半径，则月球起潮力位可写成7.1 起潮力位及其调和展开起潮力位及其调和展开220,2)(coscos1cos)()(cos11cos)(nmnnmmmmmmmmmnmnnmmAmAmmmmmmzPrrzrrrMGzrrMGATzPrrrrrMGzrrMGAT可得由OeOOmrmArm,Aqmfmrzm 7.1.2 起潮力位起潮力位 考虑到考虑到(r/rm)1/60.3，只取，只取(r/rm)3项，并考虑项，并考虑(GMm/rm)对于对于A点来说可近似为常数（仅万分之几的变化），可得点来说可近似为常数（仅万分之几的变化），可得7.1 起潮力位及其调和展开起潮力位及其调和展开)(cos)(cos)(3322mmmmmmmzPrrzPrrrMGAT22()cos1cos(cos)nmmmmmnmnmmmmMMrrTAGrzGzPzrrrr 7.1.2 起潮力位起潮力位 (2)太阳起潮力位太阳起潮力位 同理，并考虑到同理，并考虑到(r/rs)1/234000，且太阳与地球质量之，且太阳与地球质量之(Ms/Me)=333432，只取，只取(r/rm)2项，可得项，可得 由于由于zm,zs都是地球内部空间位置和时间的函数，所以，起潮位反映都是地球内部空间位置和时间的函数，所以，起潮位反映了起潮力在时间和空间上的分布规律了起潮力在时间和空间上的分布规律。7.1 起潮力位及其调和展开起潮力位及其调和展开)()()()(cos)(22ATATATzPrrrMGATsmsssss所以起潮位为 7.1.2 起潮力位起潮力位 日、月对地球的引力比较：日、月对地球的引力比较：日、月对地球（地面与球心）引力差比较：日、月对地球（地面与球心）引力差比较：gr=gAr gOr g=gA gO gmr/gsr2.18所以，所以，在太阳和月球的引力作用之中，月球的引力是产生在固体潮主要在太阳和月球的引力作用之中，月球的引力是产生在固体潮主要原因。原因。7.1 起潮力位及其调和展开起潮力位及其调和展开100/100.2/105.3/202222msmmmsFFNrGMmFNrGMmF 7.1.3 起潮力位的展开起潮力位的展开 (1)勒让德展开与杜森常数勒让德展开与杜森常数英国人杜森英国人杜森(Doodson)于于1922年引入一个常数年引入一个常数D，即，即7.1 起潮力位及其调和展开起潮力位及其调和展开20,)(cos)()(cos11nmnnmmmmnmnnmmAmzPrrrMGATzPrrrr可得由为地球半径。，为地心至月心平均距离为杜森常数，RcDcRMGDmmm3243OeOOmrmArm,Aqmfmrzm 7.1.3 起潮力位的展开起潮力位的展开 (1)勒让德展开与杜森常数勒让德展开与杜森常数 将杜森常数代入，若只取二次项或三次项，化简可得将杜森常数代入，若只取二次项或三次项，化简可得对于太阳而言，可得到类似的结果，只是杜森常数为对于太阳而言，可得到类似的结果，只是杜森常数为7.1 起潮力位及其调和展开起潮力位及其调和展开)cos2cos310()()32cos2()(33432232mmmmmmmmmmzzcRRrrcDATzRrrcDATmmmsssssDDMMccDRcDcRMGD45990.0433sm32且有为地球半径。距离，为地心至太阳质心平均为杜森常数，7.1.3 起潮力位的展开起潮力位的展开 (2)起潮位的拉普拉斯展开起潮位的拉普拉斯展开 首先讨论月球起潮位。选定天球坐标系，设首先讨论月球起潮位。选定天球坐标系，设P为北极，为北极，A为地面一点，为地面一点，M为月球在天球上的位置，为月球在天球上的位置，zm为月球对为月球对A点的地心天顶距，点的地心天顶距，为为A点纬度，点纬度，m为月球在天球上位置的赤纬，为月球在天球上位置的赤纬，Hm为月球时角。为月球时角。7.1 起潮力位及其调和展开起潮力位及其调和展开Hm 7.1.3 起潮力位的展开起潮力位的展开 (2)起潮位的拉普拉斯展开起潮位的拉普拉斯展开在球面三角形在球面三角形PAM中，有中，有7.1 起潮力位及其调和展开起潮力位及其调和展开mmmmHzcoscoscossinsincosHm 7.1.3 起潮力位的展开起潮力位的展开 (2)起潮位的拉普拉斯展开起潮位的拉普拉斯展开 根据球函数的加法定理，有根据球函数的加法定理，有若只取至二次项，化简可得若只取至二次项，化简可得7.1 起潮力位及其调和展开起潮力位及其调和展开mmknknnkmnnmnHPPknknPPzPcos)(sin)(sin)!()!(2)(sin)(sin)(cos1mmmmmmHPPHPPPPzP2cos)(sin)(sin121cos)(sin)(sin31)(sin)(sin)(cos22221212222 7.1.3 起潮力位的展开起潮力位的展开 (2)起潮位的拉普拉斯展开起潮位的拉普拉斯展开考虑到考虑到有有7.1 起潮力位及其调和展开起潮力位及其调和展开mmmmmmmPPPPPP22222212122222cos3)(sin,cos3)(sincossin3)(sin,cossin3)(sin21sin23)(sin,21sin23)(sinmmmmmmmmmmmmHHHHzP2coscoscos43cos2sin2sin431sin31sin3412coscoscos43coscossincossin321sin2321sin23)(cos222222222 7.1.3 起潮力位的展开起潮力位的展开 (2)起潮位的拉普拉斯展开起潮位的拉普拉斯展开 将杜森常数代入，化简可得将杜森常数代入，化简可得7.1 起潮力位及其调和展开起潮力位及其调和展开mmmmmmmmmmmmmmmmmHHRrrcDHHRrrcDT2coscoscoscos2sin2sinsin31sin313cos2sin2sin2coscoscos31sin31sin32222232222232 7.1.3 起潮力位的展开起潮力位的展开 (2)起潮位的拉普拉斯展开起潮位的拉普拉斯展开 同理，二阶太阳起潮位为同理，二阶太阳起潮位为7.1 起潮力位及其调和展开起潮力位及其调和展开ssssssssHHRrrcDT2coscoscoscos2sin2sinsin31sin3132222232 7.1.4 拉普拉斯三种潮拉普拉斯三种潮 (1)长周期潮长周期潮 考虑第一项考虑第一项 对于对于A点来说，点来说，sin2 是是常数常数，而，而sin2 m，sin2 s是周期波。显然，是周期波。显然，月球运行周期为一个月，其半周期约为月球运行周期为一个月，其半周期约为14天；而太阳（相对地球）天；而太阳（相对地球）运行周期为运行周期为1年，其半周期约为年，其半周期约为6个月。与此对应的潮汐现象以半个个月。与此对应的潮汐现象以半个月为周期的潮为半月潮，以半年为周期的潮为半年潮，两者统称为月为周期的潮为半月潮，以半年为周期的潮为半年潮，两者统称为长周期潮长周期潮。7.1 起潮力位及其调和展开起潮力位及其调和展开sm2222sin31sin31,sin31sin31 7.1.4 拉普拉斯三种潮拉普拉斯三种潮 (1)长周期潮长周期潮 显然，长周期潮在空间上的分布特征与纬显然，长周期潮在空间上的分布特征与纬度有关，度有关，=35 15 的纬度线将球面分成的纬度线将球面分成3个个带，即在赤道附近的这个带为正，在两极附带，即在赤道附近的这个带为正，在两极附近的带为负。近的带为负。可见，第一项为起潮位的带函数。可见，第一项为起潮位的带函数。7.1 起潮力位及其调和展开起潮力位及其调和展开35 15 7.1.4 拉普拉斯三种潮拉普拉斯三种潮 (1)长周期潮长周期潮 7.1 起潮力位及其调和展开起潮力位及其调和展开 7.1.4 拉普拉斯三种潮拉普拉斯三种潮 (2)日潮日潮 考虑第二项考虑第二项 由于地球自转，对于由于地球自转，对于A点来说，点来说，cosHm和和cosHs 是在是在sin2 sin2 m 和和sin2 sin2 s 基础上基础上,以日为以日为周期的波。周期的波。Hm和和Hs的周期为一个太阳日（的周期为一个太阳日（24小时小时50.47分），显然，分），显然，cosHm和和cosHs表现为以日为周期的表现为以日为周期的日潮日潮。7.1 起潮力位及其调和展开起潮力位及其调和展开ssmmHHcos2sin2sin,cos2sin2sin 7.1.4 拉普拉斯三种潮拉普拉斯三种潮 (2)日潮日潮 现在讨论日潮在北半球现在讨论日潮在北半球(sin2 0)某一瞬间的某一瞬间的空间分布。空间分布。当当 该项的起潮力位为零，时角等于该项的起潮力位为零，时角等于90的的2条经条经线将球面分为线将球面分为2个扇形区域，个扇形区域，若若 0，当，当cosH 0 时为涨潮区，反之（阴影区）为落潮区。为正时为涨潮区，反之（阴影区）为落潮区。为正。7.1 起潮力位及其调和展开起潮力位及其调和展开2,2smHH 7.1.4 拉普拉斯三种潮拉普拉斯三种潮 (2)日潮日潮 而当在南半球而当在南半球(sin2 0 时为涨潮区，反之（阴影区）为落潮时为涨潮区，反之（阴影区）为落潮区，且南北半球符号相同。区，且南北半球符号相同。可见，第三项为起潮位的扇函数可见，第三项为起潮位的扇函数。7.1 起潮力位及其调和展开起潮力位及其调和展开43,4smsmHHHH 7.1.5 起潮力位的杜森展开起潮力位的杜森展开潮汐波潮汐波 7.1 起潮力位及其调和展开起潮力位及其调和展开 7.1.5 起潮力位的杜森展开起潮力位的杜森展开潮汐波潮汐波 由于日、地、月三者之间的时空关系十分复杂，拉普拉斯展开给出由于日、地、月三者之间的时空关系十分复杂，拉普拉斯展开给出的三类潮汐波不能满足精细地描述固体潮特征的要求。的三类潮汐波不能满足精细地描述固体潮特征的要求。杜森利用杜森利用6个天文参数，将起潮力位展成振幅为常数的谐波的总合。个天文参数，将起潮力位展成振幅为常数的谐波的总合。在起潮力位的杜森展开式中，共有在起潮力位的杜森展开式中，共有386个谐波，每个谐波振幅都不大于个谐波，每个谐波振幅都不大于10-5D。7.1 起潮力位及其调和展开起潮力位及其调和展开 n7.0 引言引言 n7.1 起潮力位及其调和展开起潮力位及其调和展开n7.2 刚体地球的理论固体潮刚体地球的理论固体潮n7.3 理论重力固体潮计算理论重力固体潮计算 平衡潮的概念：平衡潮的概念：假设地球是一个刚体，则在起潮力的作用下，地球所产生的一切潮假设地球是一个刚体，则在起潮力的作用下，地球所产生的一切潮汐现象都称为平衡潮。汐现象都称为平衡潮。7.2 刚体地球的理论固体潮刚体地球的理论固体潮)32cos2()()32cos2()(22322232ssssmmmmzRrrcDATzRrrcDAT 7.2.1 刚体地球的重力固体潮刚体地球的重力固体潮 在地球是刚体的情况下，起潮力沿地球重力方向在地球是刚体的情况下，起潮力沿地球重力方向(即即垂直于大地水准面的方向垂直于大地水准面的方向)的分量称为刚体地球的重力的分量称为刚体地球的重力固体潮。固体潮。我们将我们将起起潮力位潮力位T(A)对对A点处向径方向求导，即点处向径方向求导，即7.2 刚体地球的理论固体潮刚体地球的理论固体潮1cos33432cos2)32cos2()(23223223zrcRDRzrcDzRrcDATRRR 7.2.1 刚体地球的重力固体潮刚体地球的重力固体潮 对于月球而言，有对于月球而言，有 对于太阳而言，有对于太阳而言，有7.2 刚体地球的理论固体潮刚体地球的理论固体潮1cos33423mmmmVmzrcRDg1cos33423ssssVszrcRDg 7.2.1 刚体地球的重力固体潮刚体地球的重力固体潮 若将上式按拉普拉斯展开，可得若将上式按拉普拉斯展开，可得 7.2 刚体地球的理论固体潮刚体地球的理论固体潮ssssssssVsHHrcRDg2coscoscoscos2sin2sinsin31sin313222223mmmmmmmmVmHHrcRDg2coscoscoscos2sin2sinsin31sin313222223ssssssssHHRrrcDT2coscoscoscos2sin2sinsin31sin3132222232 7.2.2 刚体地球的大地水准面潮汐刚体地球的大地水准面潮汐 在地球是刚体的情况下，大地水准面就是覆盖在在地球是刚体的情况下，大地水准面就是覆盖在刚体地球表面的静止海水面。当有起潮力作用时，此刚体地球表面的静止海水面。当有起潮力作用时，此大地水准面将因潮汐运动而发生形变。根据一般的原大地水准面将因潮汐运动而发生形变。根据一般的原理，形变后的大地水准而形状如图中的虚线椭圆。此理，形变后的大地水准而形状如图中的虚线椭圆。此时原来大地水准时原来大地水准面面上的上的A点上升至点上升至A 点，点，B点下降至点下降至B 点。点。AA 和和BB 称为大地水准面平衡潮高。这种现象称为大地水准面平衡潮高。这种现象称为刚体地球的称为刚体地球的大地水准面潮汐大地水准面潮汐。在在A 点的平衡潮高可表示为点的平衡潮高可表示为7.2 刚体地球的理论固体潮刚体地球的理论固体潮gATTTAWAW)()()(000斯公式可得可视为扰动位，由布容 7.2.3 刚体地球的地倾斜固体潮刚体地球的地倾斜固体潮 对于刚体地球，沿大地水准面水平方向的起潮对于刚体地球，沿大地水准面水平方向的起潮力分量力分量 gh 将引起垂线方向的变化，这就是说，地将引起垂线方向的变化，这就是说，地球上任一点的垂线方向应是重力球上任一点的垂线方向应是重力(减去起潮力垂直分减去起潮力垂直分量量)与引潮力水平分量的合力方向。由于起潮力水平与引潮力水平分量的合力方向。由于起潮力水平分量分量gh的作用使垂线方向的作用使垂线方向g偏离了偏离了 角角(一般很小），一般很小），因此可以写成：因此可以写成：=gh/g 由于垂线偏离，过由于垂线偏离，过A点的水平面发生同样的倾斜，点的水平面发生同样的倾斜，这就相当于地平面的倾斜，这种垂线方向的变化称这就相当于地平面的倾斜，这种垂线方向的变化称为为地倾斜固体潮地倾斜固体潮。7.2 刚体地球的理论固体潮刚体地球的理论固体潮 7.2.3 刚体地球的地倾斜固体潮刚体地球的地倾斜固体潮 对（月球或太阳）起潮力位水平方向求导，并对（月球或太阳）起潮力位水平方向求导，并令令=R，得，得7.2 刚体地球的理论固体潮刚体地球的理论固体潮zrcgRDggzzrcRDzTghh2sin2cossin433则有 7.2.3 刚体地球的地倾斜固体潮刚体地球的地倾斜固体潮 在实际应用和观测中，通常只能记录出地倾斜在实际应用和观测中，通常只能记录出地倾斜固体湖的南北和东西两个分量，即子午圈和卯酉固体湖的南北和东西两个分量，即子午圈和卯酉圈分量，即圈分量，即7.2 刚体地球的理论固体潮刚体地球的理论固体潮HHrcRDTgggggxyx2coscos2sincos2sin2cos231sinsin3,222300 7.2.3 刚体地球的地倾斜固体潮刚体地球的地倾斜固体潮 设设H=t0+，t0为格林威治子午圈的平月球时角为格林威治子午圈的平月球时角 为观测点的经度，则有起潮力位对为观测点的经度，则有起潮力位对卯酉圈方向求卯酉圈方向求导为导为7.2 刚体地球的理论固体潮刚体地球的理论固体潮HHrcRDTgy2sincoscossin2sinsin2cos123 n7.0 引言引言 n7.1 起潮力位及其调和展开起潮力位及其调和展开n7.2 刚体地球的理论固体潮刚体地球的理论固体潮n7.3 理论重力固体潮计算理论重力固体潮计算 7.3.0 刚体地球刚体地球理论重力固体潮计算理论重力固体潮计算 将将起起潮力位潮力位T(A)对对A点处向径方向求导可得到重力表达式，如点处向径方向求导可得到重力表达式，如 通过对通过对(c/r)及及 z 的不同方式展开，可得到地面的不同方式展开，可得到地面A点处月球、太阳的点处月球、太阳的起潮力与它们的在起潮力与它们的在A点的观测值与时间（世界时）的关系，从而可以点的观测值与时间（世界时）的关系，从而可以计算出地面各点不同时间的理论重力固体潮。计算出地面各点不同时间的理论重力固体潮。7.3 理论重力固体潮的计算理论重力固体潮的计算1cos334)(23zrcRDATR 7.3.1 刚体地球理论重力固体潮计算刚体地球理论重力固体潮计算天顶距计算方法天顶距计算方法 对于太阳，可用前面导出的公式对于太阳，可用前面导出的公式:对月球而言，为了提高精度，我们考虑在拉普拉斯展开中保留对月球而言，为了提高精度，我们考虑在拉普拉斯展开中保留n=3项，项，并对起潮力位进行求导，可得并对起潮力位进行求导，可得为了便于计算，将为了便于计算，将(c/r)用月球地平视差表示。用月球地平视差表示。7.3 理论重力固体潮的计算理论重力固体潮的计算)cos3cos5()1cos3(343423mmmmmmmmmmVmzzrccDzrcRDg1cos33423ssssVszrcRDg 7.3.1 刚体地球理论重力固体潮计算刚体地球理论重力固体潮计算天顶距计算方法天顶距计算方法 由图可知由图可知,设从月球中心刚好可以看到地球表面上设从月球中心刚好可以看到地球表面上A点，点，P被称为月球被称为月球的瞬间地平视差的瞬间地平视差(角角)，P0则为月球的平均地平视差。则为月球的平均地平视差。显然有显然有7.3 理论重力固体潮的计算理论重力固体潮的计算rRPcRPsinsin0，7.3.1 刚体地球理论重力固体潮计刚体地球理论重力固体潮计算算天顶距计算方法天顶距计算方法由由 可得可得对于太阳，有对于太阳，有 7.3 理论重力固体潮的计算理论重力固体潮的计算)31(cossinsin4)cos3cos5(sinsinsin2)31(cossinsin4sinsin23030402300ssssVsmmmmmmmmmmVmzPPRDgzzPPPRDzPPRDgPPrc 7.3.1 刚体地球理论重力固体潮计算刚体地球理论重力固体潮计算天顶距计算方法天顶距计算方法 7.3 理论重力固体潮的计算理论重力固体潮的计算算出理论固体潮值。轨道来确定，由此可计这些参数可以通过月球的函数，和固体潮理论值是按拉普拉斯展开，可知计算得到。平视差可从天文年历中而月球和太阳的瞬间地，HPPsm,79.852.207500 7.3.2 刚体地球理论重力固体潮计算刚体地球理论重力固体潮计算按杜森展开式计算按杜森展开式计算 根据杜森引入的根据杜森引入的6个天文参数：个天文参数：以角度表示的平月球地方时；以角度表示的平月球地方时；s 月球的平黄经；月球的平黄经；h 太阳的平黄经；太阳的平黄经；p 月球近地点的平黄经；月球近地点的平黄经；N N,N月球升交点的平黄经；月球升交点的平黄经；ps 地球近日点的平黄经。地球近日点的平黄经。这些参数都可表示为世界时这些参数都可表示为世界时t的函数。的函数。月球真黄经、纬度（月球真黄经、纬度（m,m）和和太阳的真黄经、纬度太阳的真黄经、纬度（s s,s s）都表示）都表示为上述参数的线性组合的幅角（见为上述参数的线性组合的幅角（见重力与固体潮教程重力与固体潮教程p233,11.40p233,11.40式）。式）。7.3 理论重力固体潮的计算理论重力固体潮的计算 7.3.2 刚体地球理论重力固体潮计算刚体地球理论重力固体潮计算按杜森展开式计算按杜森展开式计算 为了计算潮汐波方便，杜森利用天文参数和根据观测精度需要，将起为了计算潮汐波方便，杜森利用天文参数和根据观测精度需要，将起潮力位展开成潮力位展开成 386 个简单的谐波个简单的谐波（潮汐波）（潮汐波），写为，写为在月球和太阳的起潮力位中，在月球和太阳的起潮力位中，G0,1,2,3 只与观测点的纬度有关，称为大地系只与观测点的纬度有关，称为大地系数。而数。而Hi 只与月亮和太阳在天球上的位置、它们至地心的距离以及观测只与月亮和太阳在天球上的位置、它们至地心的距离以及观测点的经度有关。点的经度有关。7.3 理论重力固体潮的计算理论重力固体潮的计算3860,1,2,31cos(sin)iiiiTGH 7.3.2 刚体地球理论重力固体潮计算刚体地球理论重力固体潮计算潮汐波分类潮汐波分类 潮汐波按周期可分为潮汐波按周期可分为 长周期波长周期波，（全）日波，半日波和三分之一日，（全）日波，半日波和三分之一日波。波。7.3 理论重力固体潮的计算理论重力固体潮的计算2,equator1 3sin0,3515,pole（见（见重力学与固体潮重力学与固体潮）7.3.2 刚体地球理论重力固体潮计算刚体地球理论重力固体潮计算潮汐波分类潮汐波分类 潮汐波按周期可分为潮汐波按周期可分为 长周期波长周期波，（全）日波（全）日波，半日波和三分之一日，半日波和三分之一日波。波。7.3 理论重力固体潮的计算理论重力固体潮的计算0,sin2/4,/2,（见（见重力学与固体潮重力学与固体潮）7.3.2 刚体地球理论重力固体潮计算刚体地球理论重力固体潮计算潮汐波分类潮汐波分类 潮汐波按周期可分为潮汐波按周期可分为 长周期波，（全）日波，长周期波，（全）日波，半日波半日波和三分之一日和三分之一日波。波。7.3 理论重力固体潮的计算理论重力固体潮的计算20,cos/2,（见（见重力学与固体潮重力学与固体潮）7.3.2 刚体地球理论重力固体潮计算刚体地球理论重力固体潮计算潮汐波分类潮汐波分类 在起潮力位的调和展开中，可以得到无穷多个潮汐波。若从振幅大于在起潮力位的调和展开中，可以得到无穷多个潮汐波。若从振幅大于某个数值的潮汐波处截断，则起潮力位的调和展开由有限个潮汐波组成。某个数值的潮汐波处截断，则起潮力位的调和展开由有限个潮汐波组成。18831883年达尔文在他的调和展开中，用字母对潮汐进行了命名，如将主要年达尔文在他的调和展开中，用字母对潮汐进行了命名，如将主要长周期波长周期波取名为取名为Mm、Mf、Sa、Ssa，将，将主要日波主要日波取名为取名为K1、O1，将，将主要半主要半日波日波取名为取名为M2、S2等。这种命名方法简单并且便于记忆，但潮汐波很多，等。这种命名方法简单并且便于记忆，但潮汐波很多，字母有限，不可能用字母命名每个潮汐波。字母有限，不可能用字母命名每个潮汐波。杜森杜森19221922年在他的展开中只保留了振幅大于年在他的展开中只保留了振幅大于0.000100.00010D D 的潮汐波，每个的潮汐波，每个潮汐波的振幅计算至小数点后五位，共得潮汐波的振幅计算至小数点后五位，共得386386个潮汐波，其中长周期波个潮汐波，其中长周期波9999个，日波个，日波158158个，半日波个，半日波115115个，三分之一波个，三分之一波1414个。杜森用他提出的幅角数个。杜森用他提出的幅角数构成的编码对潮汐波进行了自然分类，每个潮汐波有一个编码，每一个编构成的编码对潮汐波进行了自然分类，每个潮汐波有一个编码，每一个编码代表一个潮汐波。码代表一个潮汐波。7.3 理论重力固体潮的计算理论重力固体潮的计算 7.3.2 刚体地球理论重力固体潮计算刚体地球理论重力固体潮计算潮汐波分类潮汐波分类 波频率相差越小的两个潮汐波，作调和分析时将它们分开所需的时间波频率相差越小的两个潮汐波，作调和分析时将它们分开所需的时间越长，将两个波分开所需的时间等于两个波相位相差越长，将两个波分开所需的时间等于两个波相位相差360360的时间间隔，的时间间隔，即即 波类码相差波类码相差1 1时，时，两个相邻波类两个相邻波类的潮汐波的角频率相差约的潮汐波的角频率相差约14.414.4/h/h，因而将不同波类潮汐被分开所需的时间为因而将不同波类潮汐被分开所需的时间为1 12 2天天；将；将两个相邻波群两个相邻波群（波群（波群码相差码相差1 1的波）分开所需的时间约等于的波）分开所需的时间约等于2929天天；将；将两个相邻分波两个相邻分波（分波码相（分波码相差差1 1的潮汐波）分开所需的时间约等于的潮汐波）分开所需的时间约等于1 1年年；而将杜森编码；而将杜森编码第六位编码相差第六位编码相差1 1的相邻潮汐波的相邻潮汐波彼此分开所需的时间约等于地球近日点的周期彼此分开所需的时间约等于地球近日点的周期2 2万年万年。7.3 理论重力固体潮的计算理论重力固体潮的计算360Tn重力固体潮在线计算网址重力固体潮在线计算网址nhttp:/www.megsystems.ca/webapps/tidecorr/tidecorr.aspxnMT80TW is a standard routine to compute tidal effects for gravity,tilt or strain.qhttp:/www.eas.slu.edu/GGP/mt80-wparicet/mt80w_exe/n地球固体潮计算程序地球固体潮计算程序nhttp:/ 地球的固体潮地球的固体潮 7.3.3 固体地球对起潮力的相应固体地球对起潮力的相应勒夫数勒夫数 如前所述，如前所述，固体潮的起潮力是唯一能够预先准确计算的作用于地固体潮的起潮力是唯一能够预先准确计算的作用于地球上的外力。球上的外力。在这一力的作用下，地球有什么样的响应？在这一力的作用下，地球有什么样的响应？实际地球并不是一个刚体，在起潮力的作用下，它会产生形变。实际地球并不是一个刚体，在起潮力的作用下，它会产生形变。因而若把地球视为一个弹性体，用球坐标表示其内部任意一点因而若把地球视为一个弹性体，用球坐标表示其内部任意一点A处的处的形变位移量形变位移量S向量为向量为7.3 理论重力固体潮的计算理论重力固体潮的计算eeeS)()()()(ASASASArr 7.3.3 固体地球对起潮力的相应固体地球对起潮力的相应勒夫数勒夫数 假设地球密度假设地球密度、拉梅常数、拉梅常数,只是半径只是半径 r 的函数，用大地水准面的函数，用大地水准面形变表示，即有形变表示，即有其中其中Hn(r),Ln(r)为为勒夫数。勒夫数。7.3 理论重力固体潮的计算理论重力固体潮的计算sin)()()()()()()()()(222gATrLASgATrLASgATrHASnnnnnnnnnr 7.3.3 固体地球对起潮力的响应固体地球对起潮力的响应勒夫数勒夫数 考虑有形变的地球模型考虑有形变的地球模型理论重力固体潮计算理论重力固体潮计算 当当r=R时，令时，令Hn(R)=hn,Ln(R)=ln，称，称hn,ln 为为n阶勒夫数。阶勒夫数。另一方面，由于地球内部发生形变，地球内部密度也将发生变化，由另一方面，由于地球内部发生形变，地球内部密度也将发生变化，由此将引起重力变化，即产生一个附加起潮力位，即此将引起重力变化，即产生一个附加起潮力位，即当当r=R时，令时，令Kn(R)=kn为为n阶勒夫数。阶勒夫数。显然，勒夫数显然，勒夫数h、k、l 与弹性地球内部密度和弹性模量有关，它们与弹性地球内部密度和弹性模量有关，它们的某种线性组合可以通过观测得到，从而可以计算出来，进而计算固的某种线性组合可以通过观测得到，从而可以计算出来，进而计算固体潮理论值。体潮理论值。7.3 理论重力固体潮的计算理论重力固体潮的计算)()()()(2ATrKATAVnnn基本概念基本概念地球的平动、起潮力、地月系统的惯性离心力、平衡潮、刚体地球的平动、起潮力、地月系统的惯性离心力、平衡潮、刚体地球的重力固体潮、大地水准面潮汐、地倾斜固体潮；地球的重力固体潮、大地水准面潮汐、地倾斜固体潮；基本关系（包括参数的意义）基本关系（包括参数的意义）起潮力位（包括太阳和月球）的（仅考虑二阶）基本表达式：起潮力位（包括太阳和月球）的（仅考虑二阶）基本表达式：现象与理解、基本技能现象与理解、基本技能三种拉普拉斯潮的名称及特征、固体潮表现形式，潮汐波分类。三种拉普拉斯潮的名称及特征、固体潮表现形式，潮汐波分类。第七章的主要知识点第七章的主要知识点)(cos)(cos)(2222ssssmmmmzPrRrMzPrRrMGAT