地球内部密度_重力加速度与压强分布公式的拟合.doc

① 地球内部密度 、重力加速度与压强分布公式的拟合 李安生 (安阳师院 物理系 ,河南 安阳 455002) 摘要 :利用地震学家布伦由地震波速度分布得出的密度 、重力加速度 、压强的一组典型数据 ,采用拟合的方 法 ,得出了较简洁的密度 、重力加速度 、压强分布的近似公式 . 这些公式对研究地球内部的其他相关性质 ,具有 重要的意义 . 关键词 :地球的密度 ;地球内的压强 ;公式的拟合 中图分类号 : P 312 文献标识码 :A 文章编号 :0258 - 7971 (2007) 05 - 0480 - 05 1 5 ,由图 1 中的曲线可看出 , 密度 ρ 的分布是不 在地球科学的研究中 ,地球的许多特性都与地 球的密度分布和压强分布有关 ,因此地球内部的密 度公式和压强公式就显得十分重要. 历史上曾由求 解克雷若 ( A . C. Clairaut ) 微分方程得到几种不同 的密度公式1 , 但由于建立微分方程的假定和地 球的层状结构并不一致 ,而求解微分方程主要依据 的是边界条件和数学规律 ,因而所得公式都不令人 满意. 近年来 ,不断有学者探讨地球的密度及其与 地球的其它相关性质的关系2 ,3 ,也有学者尝试建 立地球密度结构的数学模型4 ,但是 ,由于其所建 数学模型严重偏离公认的地球密度分布 ,所得一系 列计算公式 ,基本都不成功 ( 见本文第 3 部分) . 为 了地球相关性质的研究方便 ,本文根据地震学家布 伦 ( K. E. Bullen) 对地球的地震纵波和横波的波速 分析所得出的数据 ,采用拟合的方法导出了密度 、 重力加速度 、压强随深度变化的近似公式 ,期望能 有利于对地球内部其他相关性质的分析和研究. 连续的 ,在上地幔 B 层与 C 层的分界面及附近有 小跃变 ,在核幔边界处 ( 深 2 898 km) 有一大跃变 , 密度随深度的变化关系应大致分 3 段来考虑. 由图 1 中 的 密 度 分 布 曲 线 , 对 照 已 知 函 数 的 图 像 类 型6 ,可知在深度 33～413 km 区间 ,密度函数的形 式为 ρ = ρ0 + b ( h - h0 ) , 由表 1 的数据很容易得 出 ρ = 3 . 32 + 8 . 42 ×10 - 7 ( h - 33 ×103 ) , ( 1) 表 1 地球的密度分布 、重力加速度分布和压强分布 Tab. 1 The densit y dist ributio n , t he gravit y acceleratio n and t he p ressure dist ributio n of t he eart h ρ/ ( g·cm - 3 ) g/ ( m·s - 2 ) h/ km p/ GPa 33 413 984 2 000 2 898 幔核界面 2 898 4 000 4 980 5 120 6 371 3 . 32 3 . 64 4 . 55 5 . 11 5 . 56 9 . 846 9 . 960 9 . 966 10 . 01 10 . 73 0 . 9 14 . 1 37 . 9 87 . 0 136 1 地球内部密度与重力加速度分布公式 表 1 是布伦提出的 A ″模式中 , 地球内部几个 特殊深度 h 处物质的密度ρ和重力加速度 g 与压 强 p 的数值[ 1 ] ,表中未列入地幔以上部分的数据 , 是因为地幔以上的部分 ,无论从横向还是纵向 ,都 有较大的变化 ,不易确定一个典型的代表值 ,故拟 合的公式从深度 33 km 做起 ,其上部分不做讨论. 根据表中的数据 ,ρ 和 g 随 深 度 的 变 化 曲 线 如 图 9 . 98 11 . 42 12 . 17 12 . 25 12 . 51 10 . 73 7 . 88 4 . 78 4 . 31 0 136 247 320 328 361  力加速度 g 的变化均可粗略地看成线性的 , 即 g = a + bh 的形式. 直接利用表 1 中的数据 , 和 EXC EL 程序中提供的拟合直线参数的方法 , 可方便地得出 地球核外 、核内 2 个区域中的重力加速度计算公式 g1 = 9 . 785 + 2 . 542 ×10 - 7 h , g2 = 20 . 00 - 3 . 094 ×10 - 6 h . ( 4) ( 5) 单位 m ·s- 2 . 使用公式 ( 4) 计算深度在 ( 0 ～ 2 898 km) 区间的重力加速度 g1 的相对误差小于 3 % , 使用公式 ( 5) 计算深度在 ( 2 898 ～ 6 371 km) 区间的 重 力 加 速 度 的