固体地球物理学导论讲稿.ppt
第一页，讲稿共四十一页哦 地球重力由两部分组成，地球上任何一个物体，都同时受到地球的引力地球重力由两部分组成，地球上任何一个物体，都同时受到地球的引力F和因随地球自转而产生的惯性离心力和因随地球自转而产生的惯性离心力C的作用。的作用。由牛顿万有引力定律，有物体由牛顿万有引力定律，有物体m所受万有引力为所受万有引力为物体所受的惯性离心力为物体所受的惯性离心力为两者的矢量合为重力，即两者的矢量合为重力，即(G=6.67210-11m3/(kgs2)2mMGFrmC2CFg第二页，讲稿共四十一页哦 将物体质量去除其所受重力，可得单位质量所受到的重力将物体质量去除其所受重力，可得单位质量所受到的重力重力场，其中重力场，其中 在地球物理学中所称的重力就是指重力场强度，重力场强度实际上就是重力在地球物理学中所称的重力就是指重力场强度，重力场强度实际上就是重力加速度。加速度。衡量重力大小的单位有两个系统，一个是高斯制（衡量重力大小的单位有两个系统，一个是高斯制（CGSM)CGSM)，另一个是国际制（，另一个是国际制（SISI）。）。CGSM:Gal(CGSM:Gal(伽），伽），mGal mGal（毫伽），（毫伽），GalGal（微伽）（微伽）SI:m/sSI:m/s2 2 g.u.(10 g.u.(10-6-6 m/s m/s2 2)目前，最好的重力仪测量精度可达到微伽级。目前，最好的重力仪测量精度可达到微伽级。rC2,2MGF 第三页，讲稿共四十一页哦3.1.1 3.1.1 地球重力的积分表达地球重力的积分表达 设笛卡儿坐标系统的原点位于地设笛卡儿坐标系统的原点位于地球球心，球球心，Z Z轴为地轴，点质量元轴为地轴，点质量元A坐坐标为标为(,)，计算点，计算点p坐标为坐标为(x,y,z)，两点之间的距离为，两点之间的距离为，如图所示，如图所示。p(x,y,z)dm=dvXYZOA(,)RpRA第四页，讲稿共四十一页哦惯性离心力可写成惯性离心力可写成即重力为即重力为地球引力可写成万有引力积分形式，即地球引力可写成万有引力积分形式，即 dddzyxzGkdddzyxyGjdddzyxxGidvGdmGFeeeeeVVVVM23222232222322233)()()()()()()()()()()()(zyxzyxFkyFjxFigkgjgig)()(22jyixC22第五页，讲稿共四十一页哦3.1.2 3.1.2 地球外部重力的球地球外部重力的球谐函数表达谐函数表达 设球坐标系统的原点位设球坐标系统的原点位于地球球心，于地球球心，Z Z轴为地轴，轴为地轴，地球半径地球半径R,为为计算点计算点p坐标为坐标为(r,)，两点之间的距离为，两点之间的距离为，如图所示。，如图所示。第六页，讲稿共四十一页哦 地球重力（位）在地球外部空间的重力满足拉普拉斯方程，因此，可求解球坐地球重力（位）在地球外部空间的重力满足拉普拉斯方程，因此，可求解球坐标中的拉普拉斯方程，得到地球外部空间重力表达式。标中的拉普拉斯方程，得到地球外部空间重力表达式。0sin1sinsin112222222VrVrrVrrr222001sin21)(cos)sin()cos(),(rPmDmCrRrGMrWmnnnmmnmnn220022sin)(cos)sin()cos()1(),(rPmDmCrRnrGMrgmnnnmmnmnnr第七页，讲稿共四十一页哦 2.1.3 2.1.3 地球重力变化及主要原因地球重力变化及主要原因 显然，由于惯性离心力的作用，地球形状、内部密度等原因，地球重力是变化的，且显然，由于惯性离心力的作用，地球形状、内部密度等原因，地球重力是变化的，且总体有随纬度变化的特征。两极处最大，赤道处最小，重力并不总指向地心。总体有随纬度变化的特征。两极处最大，赤道处最小，重力并不总指向地心。引起地球表面重力变化的主要原因：引起地球表面重力变化的主要原因：地球的形状地球的形状扁椭球体引力，最大变化达扁椭球体引力，最大变化达1800mGal1800mGal；地球自转地球自转惯性离心力，最大变化达惯性离心力，最大变化达3400mGal3400mGal；地球内部物质密度分布不均匀；地球内部物质密度分布不均匀；地球表面起伏不平，最大变化达地球表面起伏不平，最大变化达1000mGal1000mGal以上；以上；太阳与月球的引力，最大变化达太阳与月球的引力，最大变化达0.3mGal0.3mGal。第八页，讲稿共四十一页哦3.3.1 3.3.1 地球重力位地球重力位 重力位是一个标量函数，可由重力各分量沿着力的方向积分得到，即重力位是一个标量函数，可由重力各分量沿着力的方向积分得到，即由重力位函数可导出重力在各个方向上的分量，即由重力位函数可导出重力在各个方向上的分量，即n为为W(x,y,z)W(x,y,z)函数内法线方向上的单位矢量。重力为沿重力位函数内法线方向上的函数内法线方向上的单位矢量。重力为沿重力位函数内法线方向上的偏导数就等于该处的重力值。偏导数就等于该处的重力值。2221rdmGWeVgnWgzWgyWgxWzyx,第九页，讲稿共四十一页哦3.3.2 3.3.2 重力等位面重力等位面 在地球重力场空间内，我们总能找到一个常数在地球重力场空间内，我们总能找到一个常数C C，使得，使得 可以证明，上面的方程为一个可以证明，上面的方程为一个空间曲面方程空间曲面方程。所以，它所代表的曲面就是一。所以，它所代表的曲面就是一个等位面。由于个等位面。由于C C的任意性，因此，地球重力位有无穷多个等位面，且彼此不的任意性，因此，地球重力位有无穷多个等位面，且彼此不相交。相交。等位面上所有点的重力位相等，但重力值不一定相等。等位面上任意一点的内法线（等位面上所有点的重力位相等，但重力值不一定相等。等位面上任意一点的内法线（场增大法线方向）为该处重力方向。场增大法线方向）为该处重力方向。重力等位面在生活中的体现：任意一个平静的水面重力等位面在生活中的体现：任意一个平静的水面水准面，即为重力等水准面，即为重力等位面。位面。CzyxW),(第十页，讲稿共四十一页哦 3.3.3 3.3.3 参考椭球面参考椭球面地球基本形状地球基本形状 由于地球是一个两极压扁的椭球体，斯托克斯在理论上证明了如果地球由于地球是一个两极压扁的椭球体，斯托克斯在理论上证明了如果地球表面重力已知，可以推导出地球表面理论公式，即与地球表面最接近的重力表面重力已知，可以推导出地球表面理论公式，即与地球表面最接近的重力等位面方程等位面方程参考椭球面。参考椭球面。由克雷罗给出的二级近似公式由克雷罗给出的二级近似公式其中其中f 为为J2 为与地球形状和转动惯量有关参数，为与地球形状和转动惯量有关参数，m赤道离心力与引力之比值。赤道离心力与引力之比值。根据参考椭球面，可以建立经纬度系统，以致地球上任何一点的位置可以用根据参考椭球面，可以建立经纬度系统，以致地球上任何一点的位置可以用经纬度来描述。经纬度来描述。经度线：过地轴的的平面与参考椭球面之交线。经度线：过地轴的的平面与参考椭球面之交线。纬度线：垂直地轴的平面与参考椭球面之交线。纬度线：垂直地轴的平面与参考椭球面之交线。)sin1(2fRrmJf21232)/(22RGMRm第十一页，讲稿共四十一页哦 纬度的不同定义：纬度的不同定义：地理纬度：地理纬度：地面任意一点上参考椭球面法线与赤道面之夹角。地面任意一点上参考椭球面法线与赤道面之夹角。地心纬度：地心纬度：地面任意一点与地心连线与赤道面之夹角。地面任意一点与地心连线与赤道面之夹角。天文纬度：天文纬度：地面任意一点上铅垂线（大地水准面法线）与赤道面之夹角。地面任意一点上铅垂线（大地水准面法线）与赤道面之夹角。第十二页，讲稿共四十一页哦 3.3.4 3.3.4 大地水准面及其形状大地水准面及其形状 大地水准面是指与大地水准面是指与“平均平均”海平面重合的水准面或重力等位面，其延伸到陆地之下海平面重合的水准面或重力等位面，其延伸到陆地之下所形成的一个封闭曲面。所形成的一个封闭曲面。确定大地水准面的形状可分两步进行：第一步是确定地球的基本形状，第二步是确定大地水准面的形状可分两步进行：第一步是确定地球的基本形状，第二步是确定大地水准面与基本形状或参考椭球面的偏差，即大地水准面的高度确定大地水准面与基本形状或参考椭球面的偏差，即大地水准面的高度N N高程异高程异常。斯托克斯首先证明了常。斯托克斯首先证明了N N可以由重力的分布计算出来。其基本思想如下：可以由重力的分布计算出来。其基本思想如下：假设实测重力位与参考面上重力位之差为假设实测重力位与参考面上重力位之差为这里这里T T称为重力干扰位。由布容斯公式可计算出大地水准面的高度称为重力干扰位。由布容斯公式可计算出大地水准面的高度N N，即，即其中其中g g0 0为参考椭球面上的（正常）重力值。为参考椭球面上的（正常）重力值。0WWT0/gTN 第十三页，讲稿共四十一页哦第十四页，讲稿共四十一页哦 3.3.5 3.3.5 垂线偏差与高程异常垂线偏差与高程异常 大地水准面与参考椭球面的差异，大地水准面与参考椭球面的差异，反映在法线方向上的差异称为垂线偏差反映在法线方向上的差异称为垂线偏差，反映在垂向距离的差异称为高程异常，反映在垂向距离的差异称为高程异常。大地水准面与参考椭球面差异的分布大地水准面与参考椭球面差异的分布是不均匀的，最大的差异可达是不均匀的，最大的差异可达 117 m 117 m，它，它与地球表面地形以及地下物质分布有关。与地球表面地形以及地下物质分布有关。由于大多数地区大地水准面与参考椭球由于大多数地区大地水准面与参考椭球面差异不大，因此在很多情况下，可将两者面差异不大，因此在很多情况下，可将两者视为相同，这时，天文纬度近似等于地理纬视为相同，这时，天文纬度近似等于地理纬度。度。第十五页，讲稿共四十一页哦第十六页，讲稿共四十一页哦第十七页，讲稿共四十一页哦3.4.1 3.4.1 正常重力正常重力 由于地球内部物质存在不均匀，地球表面也不光滑，准确地计算地球的引力由于地球内部物质存在不均匀，地球表面也不光滑，准确地计算地球的引力是不可能的。若把地球内部物质分布和表面形状理想化，即假设是不可能的。若把地球内部物质分布和表面形状理想化，即假设 地球是一个两极压扁的旋转椭球体且表面光滑；地球是一个两极压扁的旋转椭球体且表面光滑；地球内部物质密度呈层状均匀（层面共焦点，层内均匀）；地球内部物质密度呈层状均匀（层面共焦点，层内均匀）；地球是一个刚性球体，内部各质点位置不变；地球是一个刚性球体，内部各质点位置不变；地球的质量、自转角速度不变。地球的质量、自转角速度不变。利用实际观测结果，可以导出一个近似公式，称为参考椭球面（大地水准利用实际观测结果，可以导出一个近似公式，称为参考椭球面（大地水准面）上正常重力公式，即面）上正常重力公式，即)2sinsin1()(2120gge第十八页，讲稿共四十一页哦式中式中g g0 0为正常重力值，其随纬度为正常重力值，其随纬度 变化；变化；g ge e，g gp p 分别称为赤道处和两极处重分别称为赤道处和两极处重力平均值；力平均值；称为地球重力扁度称为地球重力扁度(g gp pg ge e)/g)/ge e；为地球几何扁度。为地球几何扁度。418121第十九页，讲稿共四十一页哦赫尔默特赫尔默特19191919公式：公式：卡西尼卡西尼19301930公式：公式：19791979年国际地球物理和大地测量联合会颁布的公式：年国际地球物理和大地测量联合会颁布的公式：2220m/s )2sin000007.0sin005285.0(1780520.9g2220m/s )2sin0000059.0sin0052884.0(1780490.9g2220m/s )2sin000005.0sin0053024.0(1780327.9g第二十页，讲稿共四十一页哦 3.4.2 3.4.2 重力异常重力异常 在这里，我们所指的重力值不包括重力随时间变化的那部分。一般意义上重在这里，我们所指的重力值不包括重力随时间变化的那部分。一般意义上重力异常是指在地面上观测到的重力值与正常重力值相比较后的差值。力异常是指在地面上观测到的重力值与正常重力值相比较后的差值。这样的重力异常往往意义不明确，因为造成重力异常的原因有两个方面这样的重力异常往往意义不明确，因为造成重力异常的原因有两个方面，一是地形起伏，二是地下物质分布不均。为了使重力异常具有明确的物理，一是地形起伏，二是地下物质分布不均。为了使重力异常具有明确的物理意义，根据不同的目的，可对异常进行不同的校正，可得出不同物理意义的意义，根据不同的目的，可对异常进行不同的校正，可得出不同物理意义的重力异常。重力异常。0ggg第二十一页，讲稿共四十一页哦 自由空间重力异常自由空间重力异常 自由空间重力异常是对一般意义下的重力异常进行一项自由空间重力异常是对一般意义下的重力异常进行一项“自由空间自由空间”校正，也称校正，也称为高度校正，即为高度校正，即 g gh h=0.3086h (h=0.3086h (h为测点高程，单位为测点高程，单位 m)mGal m)mGal 自由空间校正系数是正常重力在大地水准面附近垂向梯度的近似值。自由空间自由空间校正系数是正常重力在大地水准面附近垂向梯度的近似值。自由空间(Free Air)(Free Air)异常为异常为 g f=g测测 g0+gh自由空间校正只是将高程的影响去掉，其物理意义在于：反映地表物质质量和地下物质自由空间校正只是将高程的影响去掉，其物理意义在于：反映地表物质质量和地下物质密度差异的引力效应。自由空间异常是研究大地水准面密度差异的引力效应。自由空间异常是研究大地水准面地球形状的重要资料。地球形状的重要资料。第二十二页，讲稿共四十一页哦 布格重力异常布格重力异常 如果在自由空间校正的基础上，把地形引起的引力效应也去掉，得到单纯反映地下如果在自由空间校正的基础上，把地形引起的引力效应也去掉，得到单纯反映地下物质密度分布的重力异常，这个异常叫布格重力异常。物质密度分布的重力异常，这个异常叫布格重力异常。为得到布格异常，必须再进行消除地形影响的两项校正。为得到布格异常，必须再进行消除地形影响的两项校正。布格校正：布格校正：g gB B=-2=-2 G G h h =-0.0419h mGal =-0.0419h mGal (h(h为海拔高程，为海拔高程，单位单位m,为地表物质平均密度，单位为地表物质平均密度，单位 g/cm3)地形校正（地形校正（TC）：计算出测点周围地形相对平板层的起伏物质）：计算出测点周围地形相对平板层的起伏物质所引起的引力效应所引起的引力效应.布格重力异常：布格重力异常：g gB B=g=g测测 g g0 0+g gh h+g gB B+g gTCTC 第二十三页，讲稿共四十一页哦 3.4.3 3.4.3 布格重力异常与地球内部构造布格重力异常与地球内部构造 布格重力异常在地质构造上的反映布格重力异常在地质构造上的反映第二十四页，讲稿共四十一页哦 布格重力异常与地形的关系布格重力异常与地形的关系 第二十五页，讲稿共四十一页哦第二十六页，讲稿共四十一页哦第二十七页，讲稿共四十一页哦 3.5.1 3.5.1 地壳均衡概念的由来地壳均衡概念的由来 1854 1854年，英国人普拉特（年，英国人普拉特（J.H.Pratt)J.H.Pratt)在尼泊尔喜玛拉雅山附近进行大地测量在尼泊尔喜玛拉雅山附近进行大地测量时发现，喜玛拉雅山的物质不能抵消巨大的垂向偏差，由此推测喜玛拉雅山底下有时发现，喜玛拉雅山的物质不能抵消巨大的垂向偏差，由此推测喜玛拉雅山底下有巨大的低密度物质存在。巨大的低密度物质存在。第二十八页，讲稿共四十一页哦 3.5.2 3.5.2 艾里假说艾里假说第二十九页，讲稿共四十一页哦 2.5.3 2.5.3 普拉特假说普拉特假说第三十页，讲稿共四十一页哦 艾里假说是根据阿基米德原理，假设单位面积地壳的质量应等于艾里假说是根据阿基米德原理，假设单位面积地壳的质量应等于 同等面积排开岩浆同等面积排开岩浆的质量，即有等式的质量，即有等式 （T T0 0+h+t+h+t）0 0 =（T T0 0+t+t）式中式中T T0 0、t t 分别为地壳正常厚度和插入岩浆的厚度，分别为地壳正常厚度和插入岩浆的厚度，0 0、分别为地壳正常密度分别为地壳正常密度和地幔岩浆密度，和地幔岩浆密度，h h为海拔高程。为海拔高程。第三十一页，讲稿共四十一页哦因此有因此有由于由于T T0 0总是不变的，故有总是不变的，故有 对于这个问题，艾里解释说：在陆地，任意单位面积内对于这个问题，艾里解释说：在陆地，任意单位面积内“伸入伸入”到岩浆里的到岩浆里的物质所要补偿的质量等于海平面以上物质的质量。物质所要补偿的质量等于海平面以上物质的质量。若地壳平均密度为若地壳平均密度为2.672.67，岩浆平均密度为，岩浆平均密度为3.273.27，则有，则有 t=4.45 h t=4.45 h海洋地区，同理有海洋地区，同理有 h h（0 01.031.03）=t=t（0 0），则），则 t t =2.73 h=2.73 h 这里这里t t和和 h h 分别为正常地壳厚度减小量和海水深度。在艾里假说中，分别为正常地壳厚度减小量和海水深度。在艾里假说中，t t 或或 t t 被称为补偿深度。被称为补偿深度。hTt000ht00第三十二页，讲稿共四十一页哦hTT000 普拉特普拉特海福特假说是基于地壳密度变化，单海福特假说是基于地壳密度变化，单位面积地壳物质密度为位面积地壳物质密度为海洋地区有海洋地区有 无论是普拉特假说还是艾里假说，都认为有一个等压无论是普拉特假说还是艾里假说，都认为有一个等压力面存在，即地壳存在着均衡现象，尽管力面存在，即地壳存在着均衡现象，尽管“地壳均衡地壳均衡”这这个概念不是由他们提出的，但是仍然被誉为地壳均衡现象个概念不是由他们提出的，但是仍然被誉为地壳均衡现象的首次解释者。正是由于地壳均衡问题的提出，引发了人的首次解释者。正是由于地壳均衡问题的提出，引发了人们对地球内部构造的研究。们对地球内部构造的研究。hTh.T000031第三十三页，讲稿共四十一页哦3.5.4 3.5.4 均衡异常均衡异常 若把地壳视为均衡的，按照均衡理论，根据高程变化，将均衡地壳物质密度或厚度若把地壳视为均衡的，按照均衡理论，根据高程变化，将均衡地壳物质密度或厚度补偿正常地壳补偿正常地壳计算补偿物质产生的引力效应，作为均衡校正，把它从布格重计算补偿物质产生的引力效应，作为均衡校正，把它从布格重力异常中减去，即可以得到均衡异常。力异常中减去，即可以得到均衡异常。均衡正异常均衡正异常 地壳物质盈余。地壳物质盈余。均衡负异常均衡负异常 地壳物质亏缺。地壳物质亏缺。均衡异常是研究地震和地壳运动的重要资料均衡异常是研究地震和地壳运动的重要资料。第三十四页，讲稿共四十一页哦3.6 3.6 潮汐作用与固体潮潮汐作用与固体潮3.6.1 3.6.1 潮汐现象及其实质潮汐现象及其实质潮汐类型潮汐类型Half-day tide two periods of rise-fall in a lunar day(太阴日太阴日2424小时小时5454分分)Whole-day tide one period of rise-fall in a lunar day(太阴日太阴日24小时小时54分分)Mixture tide in a month,the two types of tides alternatively occur.起潮力起潮力100/100.2/105.3/202222msemmemessesFFNrMGMFNrMGMF第三十五页，讲稿共四十一页哦 起潮力起潮力We can decompose the gravity from the moon into two components.One is along the radius r and another is take the direction perpendicular to r.For the position A on the surface they can be written in gAr=GMmcos/R2 gA=-GMmsin/R2For the center of the earth O,we similarly havegOr=GMmcos/R2 gO=-GMmsin/R2100/100.2/105.3/202222msemmemessesFFNrMGMFNrMGMF第三十六页，讲稿共四十一页哦 The differences of gravities between them can be obtained fromgr=gAr-gOr=Gm(cos/R2cos/R2)Gmr(3cos2-1)/R3g=gA-gO=Gm(sin/R2sin/R2)3Gmr(sin(2)/(2R3)The gravity differences from the sun are similar to the lunar.The ratio of gravity changes at direction of r is gmr/gsr2.18Therefore,the effect form the moon is primary and that from the sun is secondary.第三十七页，讲稿共四十一页哦潮汐波潮汐波第三十八页，讲稿共四十一页哦2.6.2 2.6.2 固体潮固体潮第三十九页，讲稿共四十一页哦2.6.3 2.6.3 潮汐对重力的影响潮汐对重力的影响2.6.4 2.6.4 潮汐摩擦与地月系统的演变潮汐摩擦与地月系统的演变第四十页，讲稿共四十一页哦待续第四十一页，讲稿共四十一页哦
