地震波到达时刻模拟器 返回
地球科学·地震工程

地震波到达时刻模拟器

通过P波·S波速度差推估震源距离。直观理解大森公式和紧急地震速报的工作原理。

参数设置

地震波传播动画(震源 → 观测点)
P波到达 t_P [s]
S波到达 t_S [s]
S-P 时间差 Δt [s]
推估震源距离 d [km]

快速P波波前(蓝)与较慢S波(红)从震源同心扩展,观测点先记录P波,后记录S波。下方走时图为距离对时间的两条直线。距离用大森关系 d=(t_S−t_P)/(1/V_S−1/V_P) 推估。
验证:V_P=8, V_S=4.6 km/s, d=100 km → t_P=12.5 s, t_S=21.74 s, Δt=9.24 s, 复原 d=100.0 km。

计算结果
P波到达时刻
25.0 s
S波到达时刻
42.9 s
PS时间差 Δt
17.9 s
大森公式 推估距离
133 km
走时曲线
地震波形
Vp/Vs比分析
走时

蓝色:P波走时曲线,红色:S波走时曲线。竖点线表示所选震源距离。

理论·主要公式

$t_P = D / V_P$, $\quad t_S = D / V_S$

初期微动继续时间(PS时差)
$\Delta t = t_S - t_P = D \left(\dfrac{1}{V_S} - \dfrac{1}{V_P}\right)$

大森公式(近似)
$D \approx 7.42 \cdot \Delta t$ [km]

💬 讨论地震波

🙋
地震发生时,最初有小摇晃,然后是大摇晃,这是为什么?
🎓
这就是P波和S波的区别。P波(Primary wave)是纵波,像声波一样通过压缩和膨胀传播岩石,速度约6km/s。S波(Secondary wave)是横波,通过左右上下摇晃岩石产生大摇晃,速度约3.5km/s。快速的P波先到达导致"初期微动",慢速的S波后到达导致"主要动"。如果震源在东京地下60km,P波约10秒后到达,S波约17秒后到达。
🙋
紧急地震速报如何在S波到达前发出警报?
🎓
因为P波比S波快,震源附近的地震计捕获P波时,S波还在远处。气象厅全国约4000个地震计网络首先捕获P波,数秒内计算震源位置、深度、震级,然后在S波到达前通过手机和紧急广播发送警报。距离震源越远,预留时间越长。最危险的是"直下型"地震,震源直下几乎没有预留时间。
🙋
根据初期微动继续时间推估距离的"大森公式",这个7.42数字从何而来?
🎓
大森房吉(1900年前后)统计日本各地地震记录开发的经验公式。PS时差 $\Delta t = D(1/V_S - 1/V_P)$ 关于D的解为 $D = \Delta t / (1/V_S - 1/V_P)$。代入日本地壳平均值Vp≈7km/s、Vs≈4km/s,得到 $1/4 - 1/7 = 0.107$,倒数约9.3。7.42对应Vp≈6、Vs≈3.5,从实测数据拟合确定。现代紧急地震速报使用更精确的地壳结构模型,但大森公式至今在教育和概算中使用。
🙋
地震波知识在CAE结构分析中有用吗?
🎓
直接相关。耐震设计中"地震动输入频谱"需与建筑固有周期对应,这取决于地盘增幅特性,由S波速度决定(VS30:从地表向下30m的平均横波速度)。Abaqus、Ansys等进行地震应答分析时,从地盘S波速度剖面输入时刻历地震动。地震探查(石油天然气勘探)使用P波反射法将地下结构"CAE网格化"可视化——整个地球变成有限元模型。
🙋
S波为什么不能在液体中传播?
🎓
S波利用"剪切变形"。固体相邻粒子间有横向拉力(剪切刚度G),所以横波能传播。但液体和气体的G为零——变形后不能恢复。所以S波无法通过液体。这是发现地球内部结构的重大发现。1906年地震学家莱曼根据S波的"阴影区"发现"地球有液体外核"。外核是熔融的铁镍液体。

常见问题

可以。通过滑块或输入框调整P波速度(Vp)和S波速度(Vs)。改变速度时,到达时刻和大森公式系数实时变化,直观理解地震波传播的物理关系。
本模拟器假设理想均质地盘的理论值。实际地震因地质结构和震源深度影响而产生误差,但对于大森公式原理理解和紧急地震速报机制学习的教育目的,精度已足够。
可以。根据设置的震源距离和速度差自动计算P波检出到S波到达的预留时间。修改警报发令阈值(例如距震源○km以内),可体验实际速报系统的判断基准。
没有直接的数据链接功能,但可将输出的震源距离或到达时差手动转用为结构分析软件的输入参数(如地震动位相特性或输入延迟时间)。用于地震应答分析的事前检讨。
P波和S波的速度比(Vp/Vs)表示什么?

Vp/Vs与泊松比直接相关。$\nu = (Vp^2 - 2Vs^2) / (2(Vp^2 - Vs^2))$。花岗岩(Vp/Vs≈1.73、ν≈0.25)、饱和砂岩(Vp/Vs≈2以上、ν≈0.3~0.4)等,岩石种类和含水状态不同,Vp/Vs比不同。地震波层析成像利用此比探测火山岩浆库和流体存在。

如何推估地震深度(震源深度)?

从多个地震计的到达时差反向逆算三维位置(震源确定法)。震源距D是震央距(水平)和震源深度h的合成 $D = \sqrt{r^2 + h^2}$。深源地震(300km以上)发生在沉降板块内,日本特有深源地震面(和达-贝尼奥夫带)有观测记录。

从"初期微动继续时间"能推估规模(M)吗?

仅从时间差无法推估规模。规模推估需要用P波振幅(最大变位)或P波波形陡峭度(频率内容)。紧急地震速报在P波到达3秒内从波形振幅进行"震级推估"和"震源距推估"结合,预测摇晃强度。

什么是地震波层析成像?

多个地震·地震计对的走时数据可三维成像地球内部速度结构。医学CT扫描的地球版。沉降板块(俯冲板)因冷温度、密度高,P波速度高呈明亮;热点(地幔柱)因温度高速度低呈暗色。地球规模的CFD/结构分析。

本计算工具的精度如何?

本工具假设均一速度结构(单层模型)。实际地球是地壳、上部地幔、过渡带、下部地幔、外核、内核的多层结构,深部处地震波会屈折反射。精确走时计算使用Jeffreys-Bullen表(JB表)或iaspei91速度模型(IASP91)。

地震波到达时刻模拟器简介

地震波到达时刻模拟器基于从震源辐射的P波(纵波)和S波(横波)的传播速度差异,计算观测点各波的到达时刻。设P波速度为$V_p$,S波速度为$V_s$,震源距为$D$,则P波到达时刻$t_p$和S波到达时刻$t_s$分别表示为$t_p = D / V_p$和$t_s = D / V_s$。此时,两者的到达时差$\Delta t = t_s - t_p$被称为大森公式,整理为$\Delta t = D (1/V_s - 1/V_p)$。本模拟器利用这一关系,可视化从观测点测得的$\Delta t$反推震源距$D$的过程。另外,紧急地震速报原理也基于P波的早期检出,在S波到达前发警报,这依赖于$V_p > V_s$的速度差。通过本模拟器,可直观理解地震波传播特性和防灾技术的物理基础。

实际应用

业界实际使用范例
建筑业,大成建设和大林组等将本模拟器用于耐震设计的高度化。例如超高层建筑"阿倍野哈尔卡斯"基础设计中,利用P波·S波到达时间差,根据假定震源距离等,优化免震·制振装置的配置。铁道领域,JR东日本将此应用于新干线早期停止系统,模拟地震检知到制动作动的时间延迟。与紧急地震速报联动的列车控制可靠性得到提高。

研究·教育应用
东京大学地震研究所将大森公式可视化教材采纳。学生直观理解震源距与P波·S波速度差(约1.7倍)的关系,成为地震学入门实习的必用工具。防灾科学技术研究所对比实时波形数据和模拟结果,研究地盘结构不均质性对到达时刻的影响。教育现场,中学理科课使用本工具进行"紧急地震速报工作原理"的体验学习。

与CAE分析的联动及实务定位
本模拟器与结构分析软件"ANSYS"、"Abaqus"的输入条件生成相联动。例如桥梁的地震应答分析,用本工具推估的震源距计算出的地震动波形作为边界条件设置。相比传统的统一地震波输入,更符合现象的非线性动态分析成为可能。在实务中,建筑许可申请时的性能评价书作成中,将本工具作为假定地震动合理性的视觉证据,促进设计者与审查机关的共识形成。

常见误解与注意事项

"P波和S波到达时差越大,震源越远"似乎是这样的思路,但实际上由于地盘结构的变化会改变波的传播速度,所以并非简单的正比例关系,需要注意。另外,"P波到达后到S波到达的时间直接对应震源距"的想法也很普遍,但大森公式仅是假设均质地盘的近似公式。实际地震中因地壳不均质性及反射屈折等影响,会产生误差。此外,"紧急地震速报在P波检知后立即确定震源"的理解也不准确。实际上需要数秒的处理时间来分析多个观测点数据,在震源附近的地区,S波到达后速报才发出的可能性是存在的。需要理解这一点。

使用指南

  1. 通过震源距滑块(distSlider)设置0~500km范围,确定P波·S波传播距离
  2. 调整P波速度滑块(vpSlider)在6~8km/s,花岗岩层设为6.5km/s,玄武岩层设为7.8km/s以符合实测值
  3. 调整S波速度滑块(vsSlider)在3.5~4.5km/s,确认其约为P波速度的0.6倍,运行模拟

具体计算示例

震源距离100km情况下,P波速度6.8km/s、S波速度4.0km/s时计算,P波到达时间为14.7秒,S波到达时间为25.0秒,初期微动持续时间Δt为10.3秒。大森公式 D≈7.42×Δt 可估算 D≈7.42×10.3≈76km,Δt与震源距离成正比(距离加倍则Δt约加倍)。

实务中的注意事项