斯涅尔定律(地震波)
地震波在每层边界处发生折射:sin(i) / v₁ = sin(r) / v₂
S波在液态外核中消失,只有P波可以穿过。
实时动画演示P波和S波在地球内部各层(地壳、地幔、外核、内核)中的传播、反射和折射。可视化地震阴影区和核幔边界效应。
地震波在穿过不同地球圈层(如地幔到外核)的界面时,其传播路径会发生偏折,这遵循斯涅尔折射定律。该定律描述了入射角与速度变化的关系。
$$ \frac{\sin i}{v_1}= \frac{\sin r}{v_2}= p $$其中,$i$ 是入射角,$r$ 是折射角,$v_1$ 和 $v_2$ 分别是波在上下两层介质中的传播速度。$p$ 称为射线参数,对于同一条射线是常数。当波从高速层进入低速层(如地幔进入外核)时,折射角会小于入射角,导致射线向法线方向弯曲。
P波和S波的传播速度由介质的弹性性质决定。对于均匀、各向同性的固体,其速度公式如下:
$$ v_P = \sqrt{\frac{K + \frac{4}{3}G}{\rho}}, \quad v_S = \sqrt{\frac{G}{\rho}} $$$v_P$ 是P波速度,$v_S$ 是S波速度。$K$ 是体积模量(反映抗压缩能力),$G$ 是剪切模量(反映抗剪切变形能力),$\rho$ 是密度。在液态外核中,$G = 0$,因此 $v_S = 0$,S波无法传播。
地震预警系统:利用P波传播速度快、破坏力小的特点。当地震仪首先检测到P波时,系统能抢在破坏性更大的S波和面波到达前几秒到几十秒,向远方城市发出预警,为高铁减速、电梯停靠、工厂停机等争取关键时间。
地球内部结构探测:这是地震波最经典的应用。通过分析全球地震台网记录到的P波和S波走时、振幅以及阴影区模式,科学家绘制出了地球内部的地壳、地幔、液态外核和固态内核的分层结构图,这被称为“给地球做CT”。
石油与矿产勘探:工程现场常见的是人工地震勘探。通过在地表制造小型震动(如爆破),并接收从地下不同岩层界面反射回来的地震波,可以推断地下岩层的形状和性质,从而定位油气储层或矿藏位置。
核试验监测:《全面禁止核试验条约》的国际监测系统包含一个全球地震台网。通过精确分析地下核爆炸产生的地震波信号(其波形特征与天然地震不同),可以监测和鉴别是否发生了违规的地下核试验。
开始使用本模拟器时,有几个需要注意的要点。首先,请理解“波速在各层中固定”这一前提。实际地球中,即使在同一地幔内,速度也会随深度对应的压力、温度、密度变化而连续增加。模拟器将其简化为“分层”,因此边界处的急剧折射会被凸显。实际数据中的路径可能会更加平滑。
其次,请勿混淆“震中距”与“至观测点的距离”。模拟器中的“震中距”是指从震源到地表某点的角距离(例如103度)。但实际震源深度较深时,波传播至地表的路径会发生变化。例如对于深度600公里的地震,阴影区的范围会与模拟器(假设震源位于地表)的情况不同。请务必养成时刻关注“震源深度”的习惯。
最后,请勿过度解读振幅(震级)的变化。虽然可以通过滑块调整振幅,但这仅是简化显示。实际地震波的振幅不仅受传播距离影响,还会因地基放大效应、散射和衰减而发生显著变化。即使在模拟器中增大S波振幅,其“无法在液态外核中传播”的根本原理也不会改变。关注这种“能否传播”的本质差异,才是学习的要诀。
本模拟器背后的波动传播原理,是CAE工程师必须了解的多领域基础。首当其冲的是“无损检测”。例如,在利用超声波探伤检测结构内部缺陷时,会区分使用类似P波的纵波和类似S波的横波,以确定缺陷类型和位置。波在金属中的传播,以及在异质材料边界处的折射、反射现象,与地震波遵循完全相同的物理规律。
另一个领域是“岩土工程”与“土木结构抗震设计”。这里直接应用了模拟器中看到的“地层分层结构”。需要计算在地表附近软弱层(相当于地壳)传播的波,如何在坚硬岩层(相当于地幔)边界处折射,以及如何在地表被放大。这构成了评估建筑物或桥梁固有周期是否与之发生共振的地基响应分析的起点。
进一步拓宽视野,其原理也相通于“声学工程”与“水下声呐”。在水中S波无法传播(因为是液体),只有P波(声波)可以传播。潜艇声呐探测海中物体时,需要考虑由水温、盐度差异导致的声速分层结构(深海声道),这同样可以用斯涅尔定律描述。地球物理学中积累的波动传播技术,正是多物理场模拟的通用语言。
熟悉本模拟器后,下一步可以尝试挑战“反问题”的思维方式。模拟器是在已知地球结构(速度模型)的前提下计算波的路径(正问题)。但在实际工作中,需要根据观测点获得的P波、S波到达时间数据,反推震源位置和地下结构。解决此反问题的基本方法便是分析“走时曲线”。通过绘制震中距与到达时间关系的曲线形态,可以推断速度结构。在模拟器中尝试改变各种参数,并亲手绘制简化的走时曲线草图,会极大地加深理解。
如果想稍微深入探究数学背景,可以接触描述波动的偏微分方程。地震波的传播可以用由弹性体运动方程推导出的波动方程来表示。一维简单情况下的形式如下: $$ \frac{\partial^2 u}{\partial t^2} = v^2 \frac{\partial^2 u}{\partial x^2} $$ 其中 $u$ 是位移,$v$ 是波速。下一步可以研究该方程的解如何呈现行波形式,以及波速 $v$ 如何与前述 $v_P$ 和 $v_S$ 的公式相关联。理解这一点后,就应该能开始阅读基于FEM(有限元法)的地震波模拟论文了。
最后,推荐一个具体的进阶主题:“面波(瑞利波、勒夫波)”。模拟器处理的是在地球内部传播的“体波”(P波、S波),但实际地震中造成强烈摇晃和长周期振动的,是在地表附近传播的面波。通过分析这种波的频散特性(速度随周期变化的性质),可以详细推断地壳厚度和浅部结构。将视角从内部结构扩展到地表现象,你会更切实地感受到其对地震工程的应用价值。