本页目录

深部 I · 用地震波看见地球内部

地震波不是一束直穿地球的探照灯。它会反射、折射、转换和衰减;但不同波型对介质的响应不同,正因为这种“挑剔”,到时差才可以把看不见的内部变成可检验的结构证据。

学习层:从一张地震记录反推一层介质

1. 具体情境:远处台站先收到哪一列波

一次地震发生后,台站记录到先到的短周期扰动和随后到达的较强横向运动。为了练习成像,我们先把震源到台站的路径压成一段长度,并给纵波、横波各一个平均速度;再比较到时和两波之间的间隔。

这是教学路径设定,不是某次地震的定位结果。真实路径可能穿过不同温压下的岩石,速度随深度和方向变化;实验故意把复杂性留在“模型边界”一节,而不是在默认数字里假装它已经被解决。

2. 揭示前预测:速度、介质和影区

打开实验前,先回答:

  1. 在同一路径上,若 \(v_P>v_S\),哪一个相位先到?两者间隔会怎样随距离变化?
  2. 横波需要剪切刚度。若路径进入液态外核,直达 S 波应继续传播,还是出现缺失?
  3. 仅凭一台站的一个到时差,能否唯一确定一个球壳界面?

点击揭示后,实验同时显示波的传播路径、到时条和模型判定。右侧的“不可见”表示该简化路径不产生直达 S 波,不表示台站完全没有任何后续波形。

3. 公式桥:从局部速度到旅行时间

均匀介质的最小模型是

\[ t_P=\frac{L}{v_P}, \qquad t_S=\frac{L}{v_S}, \qquad \Delta t=t_S-t_P. \]

若 \(v_S>0\),且速度保持常数,则

\[ \Delta t=L\left(\frac{1}{v_S}-\frac{1}{v_P}\right). \]

所以到时差对路径长度近似线性增长;但它只在相同起点、相同路径假设和平均速度下才有这个简单形式。一般旅行时间应写成

\[ t=\int_{\mathrm{path}}\frac{ds}{v(s)}, \]

真正的射线路径还要满足费马原理。平面界面的折射近似常用 Snell 关系

\[ p=\frac{\sin i}{v}, \]

其中 \(p\) 是射线参数。速度随深度变化会让直线路径弯曲,界面还会产生反射波和 P-S 转换波。

数量级上,地幔中的地震波速度是每秒几千米,地球半径约 \(6.4\times10^3\) km;跨越行星尺度的到时自然在分钟到小时量级。这里的数值只帮助建立刻度,不能替代具体区域的速度模型。

4. 互动实验:把 P/S 到时和内部剖面并排看

无 JavaScript 时的静态 fallback:教学设定为路径 \(L=4000\) km、\(v_P=8.0\) km/s、\(v_S=4.5\) km/s,并假定路径位于可传播横波的固体介质中。则 \(t_P=500\) s(8.33 min),\(t_S\approx888.9\) s(14.81 min),\(\Delta t\approx388.9\) s(6.48 min)。若把路径介质切换为液态外核,简化模型将 \(v_S\) 置为 0,直达 S 到时标为“不可用”。

量 默认读数 单位 / 解释
路径长度 \(L\) 4000 km,教学路径
\(v_P\) 8.0 km/s,平均纵波速度
\(v_S\) 4.5 km/s,平均横波速度
\(t_P\) 500.0 s,均匀路径近似
\(t_S\) 888.9 s,固体路径近似
\(\Delta t\) 388.9 s,S 减 P
液态路径 S 不可直达 剪切刚度代理为零

图示中的两条射线是概念路径,不是地震层析成像的真实射线追踪。

5. 边界:到时反演不是把地球剖开

  • 平均速度退化:更长的路径配合更快的速度,可能给出相似到时;单个到时通常不能唯一确定速度和界面。
  • 波型不等于单一物质标签:S 波缺失支持路径上存在低剪切刚度区域,但还需考虑几何、衰减、噪声、震源辐射和台站覆盖。
  • 路径会弯:把 \(L/v\) 当成完整模型只适用于教学或粗略估算。真实模型要处理速度梯度、各向异性、反射转换和三维结构。
  • 观测有基准:地震目录、台站位置、发震时刻、速度模型版本和拾取误差都要记录;动态地震事件数字必须标明事件日期、产品版本和核验日期。

因此,“液态外核”的结论来自大量台站、多个震相和影区几何的联合证据,不是一个按钮或一个单站读数单独证明的。

6. 迁移:从地震波到其他成像问题

把地震波换成医学超声、主动源地震或电磁波,迁移的核心都是:已知激励是什么,观测到传播时间或振幅的哪一部分,模型中的速度或阻抗怎样映射到空间结构?

迁移题:若两条不同路径给出相同到时,增加台站为什么有助于打破速度与界面深度的退化?若一个区域有明显各向异性,为什么只用一个标量 \(v\) 可能把织构误判成几何界面?下一章会把“内部运动”从波的路径转成板块边界的相对位移。

地球内部的分层与尺度

图 3.1球壳剖面把“层在哪里”与“波能否通过”放在同一画面:P 波可穿过液态外核,而 S 波在外核边界停止;各层半径与路径是概念示意,非真实比例。

波的路径如何留下影区

图 3.2沿“介质性质 → 波速 → 缺失相位 → 影区几何”读取推理链,理解为什么单台站不能独立切出地球内部;射线为教学路径,不是实际层析追踪。

1. P 波、S 波和弹性参数

在各向同性线性弹性介质中,纵波和横波速度分别受体积模量、剪切模量和密度控制。可用的结构性表达为

\[ v_P=\sqrt{\frac{K+4\mu/3}{\rho}}, \qquad v_S=\sqrt{\frac{\mu}{\rho}}. \]

液体近似没有静态剪切模量,即 \(\mu\approx0\),因此横波不能以直达体波形式传播;但液体仍可传播压缩型 P 波。高压下的相变、温度、成分和部分熔融会同时改变 \(K\)、\(\mu\) 和 \(\rho\)。

2. 影区是一组几何约束

地震学中的“影区”不是简单的黑洞,而是某类直接相位在某些角距范围内缺失或明显变弱。外核边界使 P 波折射、S 波无法穿过液态层;内核又提供了反射、折射和转换的附加相位。研究者用全球台站的相位到时和振幅共同拟合径向速度模型。

3. 层析成像与模型选择

地震层析把许多路径的旅行时间偏差写成速度扰动的积分问题。路径有限、分辨率不均和参数相关使反演需要正则化;一个平滑模型并不等于自然界真的平滑,而是数据与先验共同决定的解。

USGS Earthquake Hazards、IRIS/地震学数据中心以及地球物理教材是查阅波型、目录和速度模型的优先入口。实时事件列表会更新;若引用具体事件,必须同时标事件日期、数据产品或目录版本和核验日期。本页只使用教学输入。

4. 小结

地震成像的逻辑是“波型差异 → 到时与影区 → 多路径约束 → 速度模型”,每一步都有误差和边界。下一章将问:地幔的慢运动如何在地表表现为可测的板块相对速度?