地球化学：第四讲 地球内部压力分布特征.ppt

四 地球内部压力分布特征 1 压力的单位 压力是指作用于单位面积上的力，国际单位是Pa，1 Pa=1 N/m2，即1N的力作用于1 m2面积上所产生的力。 在高压实验技术中常用以下单位： 标准大气压，atm，1 atm=1.01325×105 Pa； 公斤力·厘米-2，kgf·cm-2，1 kgf·cm-2=9.80665×104 Pa； 巴，bar，1 bar=1.0×105 Pa； 兆帕，MPa，1 MPa=10 bar； 千巴，Kbar，1 Kbar=1000 bar； 吉帕，GPa，1 GPa=10 Kbar= 10000 bar； 2 弹性力学基础 2.1 体应变和剪应变 应力是物体中各部分之间相互作用的内力。设想在物体中存在一截面?S，设被此面元分开的两部分物质之间的作用力和反作用力分别为?f和-?f，则在此截面上的应力定义为如下矢量： ?=df/ds 在MKS单位制中应力的单位为N/m2，即Pa。 在短暂的作用下，地球介质可以看作是弹性体，根据地震学测量直接获得的弹性波的传播速度，就可以进一步得到其它的力学常数。 在固体中一个截面上的应力一般不与此截面垂直，可以将其分解为法向分量?┴和切向分量?//，即正应力和剪切应力。 与此相对应，固体的应变也存在两种基本形式，即与正应力相对应的体应变和纯剪切应力相对应的剪切应变。 体应变定义为体积的相对变化： ε体应变=ΔV/V 在弹性限度内，正应力与体应变成正比：?┴=K·ε体应变=K·ΔV/V 其中K为体变模量，单位与压力相同，其倒数k=1/K称为压缩系数。 剪切应变定义为剪变角： ε剪应变=ε 在弹性限度内，剪切应力与剪切应变成正比： ?//=G·ε剪应变 其中G为剪切模量，单位与压力相同。 2.2 弹性波 地震波在岩石中传播时，即有纵波即P波，也有横波即S波。 对于纵波：VP2=（K+4G/3）/ρ 对于横波：VS2=（G）/ρ 由于K和G都是正的，这表明纵波波速总是大于横波波速。 对于液体和气体，其G=0，VS=0。 根据测得的弹性波速，就可以进一步得出一系列地球深部的物理性质参数。如： 地震参数Ф：Ф = K/ρ = VP2 - 4/3·VS2 根据该参数可以计算得到密度； 体变模量K：K=ρ[ VP2 - 4/3·VS2] 从K的数量级可以区分气体和流体，一般105Pa为气体，109Pa或更大的为液体或固体； 剪切模量G：G =ρVS2 一般流体的G数值很小，对于理想流体G=0；金属和岩石的G值一般为1013-1015Pa。根据地球外核的G接近于零，由此推断外核应当为液态。 此外，还可以进一步计算得到扬氏模量和泊松比。 图1 地球各层圈的深度、密度、波速和温度分布（Codie,1982） 3 地球深部的弹性波特征及地球内部的层圈结构 表1 地球分层简表 4 地球深部的密度分布 [1] Williamson-Adama方程 该方程的假设： a 假定地球为球对称型； b 地球内部由均匀层圈组成（各层圈内部是均匀的）； c 地球内部的物质基本处于流体静压状态； d 只考虑随深度增加的静压力对密度的影响，忽略温度变化对密度的影响。 dρ/dr = -Gmρ/r2Ф 其中G为重力常数，m为半径r内的地球质量，ρ为半径r处的地球密度，Ф为地震参数。 [2]PREM地球模型 PREM地球模型的全称为Preliminary Reference Earth Model，是Dziewonski和Anderson（1981）在前人工作的基础上完成的，也是目前被广泛应用的一种地球模型。该模型在形式上更为简化，对于地球最上部1000 km范围采用了平均的方法，如洋壳占地球表面的2/3，莫霍面深度为11 km；陆壳占地球表面的1/3，莫霍面深度为35 km，因此取莫霍面平均深度为19 km。 图2 PREM地球模型的图解（Poirier,1991） 5 地球内部重力和压力分配模式 根据得到的密度分布，就可以根据以下公式分别计算出地球内部重力和压力在半径方向的分布： dg/dr = 4πG(ρ-2ρm/3) 其中，g为重力，G为重力常数，ρm为半径r范围内地球质量为m时的平均密度，ρ为随r变化的密度。 dP/dr = -G·mρ/r2 其中，P为压力，G为重力常数，m为半径r范围内地球的质量，ρ为随r变化的密度。 图3 地球内部的压力-深度曲线（Mason and Moore,1982） 压力对物质性质和结构的影响