震级不只是数字:从里氏到矩震级,读懂地震能量与预警
发布时间:2026/9/9 8:17:49 作者:尧图编辑部 阅读量:1,286

我常年跟地震监测、灾害评估这类数据打交道身边朋友知道我搞这个逢年过节聊起热点新闻总爱问一句“刚那个地震到底几级啊”。说实话“magnitude”这个词在公开语料里几乎就是“地震震级”的代名词但它背后的门道远比一个数字复杂得多。震级不是“烈度”不是一个固定不变的测量结果更不是像温度计插进水里读个数那么简单。它背后牵涉到弹性波能量释放、仪器响应函数、震源机制反演甚至还有一套看似简单却极其容易踩坑的台网测报流程。这篇文章打算把“magnitude”彻底拆开来讲。我会先从震级的物理本质说起再梳理里氏、面波、体波、矩震级这套“家族谱系”然后带你过一遍一次地震从波形记录到最终定级的完整流程接着集中回答几个公众和初学者最爱问的“震级误区”最后从工程防灾视角聊聊为什么震级参数在地震预警和结构抗震设计里那么关键。无论你是做地球物理的学生、搞GIS的开发、做灾害保险的风控还是单纯想搞明白“为什么有时候官方说6.8级有时候说7.0级”这篇应该都能给到你想要的答案。1. 震级到底在测什么先搞清楚“能量”这件事很多人的第一反应震级是不是衡量“地面摇晃得有多厉害”不是。震级衡量的是地震在震源处释放了多少弹性应变能。换句话说它描述的是“这场地震有多大”而不是“你脚下这块地动得有多剧烈”。这就要回到震级的“鼻祖”——里氏震级Local MagnitudeML。1935年查尔斯·里克特Charles Richter当年在加州做地震目录时用的是伍德-安德森扭力地震仪定义了一个非常简单粗暴的公式ML log10(A) - log10(A0)其中A是地震波在距震中100公里处、由标准伍德-安德森地震仪记录到的最大振幅以微米计A0是某一个零级地震在该距离下的振幅参考值。由于对数直接参与计算震级每增加1意味着地面最大振幅放大10倍。但关键来了振幅放大10倍不等于能量放大10倍。能量的释放量跟振幅之间是 3/2 次幂的关系。为什么我尽量通俗地解释地震波携带的能量与质点运动速度的平方成正比而速度又与振幅和频率的乘积相关。通过一系列对地震波能量密度的积分推导最后得到一条工程上广泛使用的经验关系log10(E) 1.5 × ML 4.8其中E的单位是焦耳。算一下就知道震级每差1级log10(E) 就差了1.5意味着能量相差 10^1.5 ≈ 31.62 倍。差2级就是1000倍差3级就是31623倍。这也是为什么一次7级地震释放的能量远远不是6级地震“稍微大一点”那么简单而是31倍多的跨越。换句话说如果一次6级地震已经能让老旧房屋开裂那么7级地震的破坏潜力完全是指数级抬升。有了这个能量公式很多业务上遇到的现象就能讲通了。比如为什么6级以下的小地震遍布全球、几乎天天都有但7级以上就屈指可数因为每提升一个震级对地壳积累的弹性应变能的要求要放大30多倍断层得越长、破裂持续时间越久才能攒出这么大的能量。提示震级本身是对数标度没有上限。负震级也是合法存在的比如一些极微小的矿井微震或人工震源如爆破作业可以定到 M -1 甚至 M -2而理论上地壳能积累的最大震级通常认为在9.5级左右历史上观测到的最大纪录是1960年的9.5级大地震。2. 震级家族为什么不能只用一把尺子量所有地震任何用过地震目录的人都会发现震级后缀五花八门mb、Ms、ML、Mw、Me……不少人以为这纯粹是机构习惯不同其实背后是“不同频段的地震波对不同大小地震的敏感度不同”这一物理现实在起作用。标度全称测量依据适用场景主要局限MLLocal Magnitude里氏/近震震级短周期仪器记录的S波最大振幅震中距100km附近地方性中小地震对大地震会“饱和”大震测不准mbBody-wave Magnitude体波震级P波初至后数秒内最大振幅周期约1秒左右远震快速测定常用于地震监测速报易受震源深度和介质衰减影响大震同样饱和MsSurface-wave Magnitude面波震级周期20秒左右的面波最大振幅中等及较大地震的常规测定对深源地震不适用深震面波很弱MwMoment Magnitude矩震级地震矩M0反演得到大震、深震、海啸地震的最终定标需要较长时间和更多台站数据实时性差MeEnergy Magnitude能量震级地震辐射能量的直接估计与Mw配合用于描述大地震错动与能量特征依赖区域衰减模型不确定性较大从这张表就能看出没有任何单一震级标尺可以覆盖所有地震。小地震的高频信号很丰富但到了大地震1秒周期的短周期波会“封顶”——因为地震破裂面的物理尺寸限制了高频辐射的上限于是mb和ML在大震时会卡在一个值附近不再升高这就是常说的高频震级饱和现象。这也是为什么现代地震学遇到大地震最终都以 Mw矩震级来定案。Mw的物理基础是地震矩M0而M0的定义极其简洁M0 μ × A × D也就是断层介质的剪切模量 μ乘以破裂面面积 A再乘以平均位错量 D。它直接从震源机制出发描述的是断层整体破裂产生的“静力学等效作用力”几乎不受仪器频段限制。然后由 Hanks 和 Kanamori 在1979年左右确定的转换公式Mw (2/3) × log10(M0) - 6.06M0以N·m为单位时常数项通常用6.06若以dyn·cm为单位则为10.7。看到这个公式我很想多说一句2/3这个系数不是拍脑袋定的它就是为了让Mw和面波震级Ms在大震范围内尽量接近。正因为有这样的“向后兼容”Mw公布出来时公众不会觉得数字跨度诡异机构也能在历史目录中平滑衔接。在实际业务中一个大地震发生后的数分钟到数小时你会看到“速报震级”从6.6改到6.9最后又发布“最终矩震级7.0”。这并不代表之前测错了而是随着台网数据不断补充反演模型从区域模型过渡到全球模型震源机制从单一点源逐步展开成有限断层定级精度自然一步步收敛。业内资深测报员看到震级修订非但不会觉得难堪反而会认为这是数据质量正常的演进过程。3. 一次地震到底怎么被“定级”完整流程拆给你看从地震波形到最终发布一个震级大概要经过定位、测量、归算、复核四个环节。这里我以日常测报中心的视角展开讲你会发现这份工作很像做一次精密的测算而不是单纯“看仪器读数”。3.1 第一步先从波形图上找到地震波任何台站收到的原始记录都是连续时间序列地震到时检测靠的是P波和S波的到时差。P波纵波传播速度快一般6到8公里每秒S波横波速度慢大约3.5到4.6公里每秒。台站离震中越远P波和S波的到时差就越大。定震中位置最基本的方法就是三分法至少三个台站各自给出震中距然后画圆求交点。现代台网则直接通过震相走时表做网格搜索或非线性反演一般几秒内就能给出初步位置。但对震级估算来说位置精度比定位精度更影响振幅测量因为振幅随距离的衰减关系很敏感位置偏了算出来的震级偏差会非常大。3.2 第二步读振幅、测周期、套公式假定我们已经在一组垂直向短周期记录中识别出了清晰的P波震相接下来要量取P波初至后数秒内的最大振幅对mb而言以及对应的周期。然后把振幅、周期、震中距丢进标准的量规函数例如mb log10(A/T) Q(Δ, h)其中A是最大位移振幅T是周期Q(Δ, h)是依赖于震中距Δ和震源深度h的“量规函数”负责把振幅衰减和几何扩散效应扣除掉。如果用面波测Ms则取20秒左右周期的Rayleigh波最大振幅代入类似的公式。这里要注意台站的仪器响应必须事先做过标定数字波形还要严格转换到“真实地面运动位移”。实际工作中我见过不少新手把速度量velocity当成位移量displacement读结果算出来的震级整整偏大好几倍。仪器响应扣除这块真是一步错步步错。3.3 第三步多台平均剔除异常值单台震级往往会因为台站场地、方位性、路径衰减而波动所以业务上采取的原则是“多台平均”。比如某个地震触发了20个台站每个台都计算出一个单台震级会形成一个数值分布。统计上取中位数或截尾均值作为该机构的最终测定值。如果某几个台站的测定结果与其他台站系统性偏离超过0.3级就要检查是不是存在场地放大效应比如建在软土层上的台站容易记录到明显放大的振幅或者仪器限幅。这类异常台数据通常会被剔除并在测报日志上标注原因。3.4 第四步大震必做的矩震级反演对于5.5级以上的显著地震现代测报中心不会只报一个里氏震级或面波震级就交差。他们会在获得全球或区域地震波形数据后做点源/有限断层矩张量反演。矩张量反演通俗讲就是用理论格林函数合成出“假如这个地震按某种机制破裂各台站应该记录到什么波形”再与实际波形对比不断调整震源参数直到残差最小。这个过程输出地震矩M0进而算得Mw。由于反演结果基本摆脱了高频饱和的影响对于七八级的大地震特别可靠。这里顺便分享一个实操中的体会震级测定不是“越新算法就越好”。在很多地震速报系统里高频P波震级如 Pd、τc 等预警参数和传统量规函数震级是同时并行、交叉核验的。P波预警震级追求的是“快”——用前3秒到10秒的信号快速估算最终震级范围辅以几分钟内的面波/体波震级做校正。两条曲线相互印证才能让最终对外发布的震级既有速度又有精度。4. 最容易误解的5个“震级”问题一次说清楚网上关于震级的误区实在太多了我把高频出现的几个问题集中整理到这里有的适合普通读者有的对初学者和跨行业从业者尤其有用。问题一为什么有的报道说6.8级有的说7.0级极大可能是因为不同机构用了不同标度。比如同一场地震美国地质调查局可能先发布Mw 6.9国内台网发布Ms 6.8日本气象厅又发Mj 7.0日本自有震级标度。它们测的物理量有细微差别换算关系也不是严格恒定的。尤其是震级测定本身带 ±0.2 到 ±0.3 的不确定性所以数字有差异是常态。问题二震级和烈度是一回事吗绝对不是。震级是能量大小的客观量度一个地震只有一个震级。烈度Intensity描述的是地表破坏程度和人的感受同一个地震在不同地点烈度完全不同从“无感”到“建筑物毁灭”可以横跨好几个烈度级别。我国使用的是十二度烈度表如 I、II、III……XII日本用七度制欧洲用EMS-98互相不能直接划等号。问题三为什么有的6级地震破坏力比8级还大这种情况是存在的最典型的就是震源深度。6级浅源地震震源深度5公里距离地面人群极近能量几乎全部直接释放到地表而8级深源地震如果发生在600公里深处传到地表时能量已经大幅衰减地面破坏能力反而有限。“震级”只管能量总和不管能量传输路径上的耗散。问题四是不是只有震级大的地震才会引发海啸震级大是一个前提但不是充分条件。海啸取决于震源是否发生显著的垂直位错以及震源是否位于板块边界或海底地形突变处。有一次发生在深海沟的8.0级走滑型地震几乎没有触发海啸而有些7.5级左右的逆冲型海底地震却能掀出数米高的海啸。所以做海啸预警时不仅要看magnitude还要看震源机制。问题五为什么同一台地震仪两次记录到的振幅一样官方震级却不一样振幅一样不代表震级一样因为还要考虑震中距和震源深度。两台仪器记录的振幅即使完全相同离震中近的那个会被认为释放能量更小因为它的波没有走太远就衰减到同样的振幅。这其实是最基础的“振幅-距离补偿”思想搞清楚这一点震级公式就不再是一堆魔幻数字了。5. 工程与防灾视角震级不只是新闻里的一个数我这些年做抗震和灾害应急相关项目最深的一个体会是对普通人来说震级只是“新闻播报里的一个数”但对工程抗震、地震保险、应急决策系统来说震级是整套参数体系里最基座的变量。5.1 从震级到设计地震动参数工程结构抗震设计时不会直接拿“8级地震”去做结构验算而是用地震动参数如峰值地面加速度PGA、反应谱、地震影响系数曲线等来评估。做地震危险性分析时常用两种路径确定性地震危险性分析法DSHA和概率地震危险性分析法PSHA。PSHA的核心输入之一就是震级-频度关系。著名的古登堡-里克特Gutenberg-Richter关系log10(N) a - bM意味着震级大于等于M的地震数量呈指数衰减。b值通常接近1表示每降低1级地震数量约增为10倍。这个关系不仅用于长期地震活动性研究也用于评估某场地在多长时间内、可能遭遇到多大地震的概率。做重大工程核电站、特大桥、大坝时设计方甚至要通过“设定地震”方法反推一个最不利震中距和震级组合来生成人工地震动时程。注意工程抗震领域强调的是“设防地震动”而非“震级”因为同样的Mw可以因破裂方向性效应、场地效应而在地表产生截然不同的地震动。这也是在防灾项目中我为什么反复强调“不宜只看震级就给承灾体下结论”。5.2 地震预警如何用P波“抢跑”地震预警的原理其实一句话地震波传播需要时间而电磁波几乎瞬时到达。台网在震中附近捕捉到P波后依据P波前几秒的特征估算震级抢在S波到达更远城市之前发出警报。问题在于用P波前3秒到10秒的信号去估震级本质上是个“外推法”。P波刚开始破裂时你看到的只是整个断层破裂的一小段怎么预测整个事件最终有多大这是预警算法最大的挑战。实践中常用两类P波参数一是P波位移幅值Pd二是P波卓越周期τc。Pd对中等地震的震级估算比较准但对大地震容易低估τc则能更快地反映震源尺度但容易受信噪比影响。更前沿的做法是结合多台站P波信息做联合反演甚至用机器学习在历史地震P波波形上训练回归模型以期望把震级估算的误差再压一压。我在实际项目中踩过的坑是过度依赖单一P波特征参数结果在某些近场强震记录里由于P波和S波只差不到一两秒P波窗口还没读完S波就到了波形混叠严重参数提取直接失真。后来改成“P波初判 扩展窗口动态修正 触发后持续更新震级”的机制明显稳得多。5.3 震级参数未来还有哪些拓展空间深部学习尤其是图神经网络在震源参数反演上的应用正在变热。经典矩张量反演需要人工挑选波形、剔除坏道、设定滤波频带而当前一些研究机构已经在尝试让神经网络直接“端到端”从三维波形数据映射到Mw和震源机制全程不需要显式人工作业。这种方案一旦成熟对地震海啸预警和灾害评估系统的自动化提升会非常可观。另一个很值得关注的方向是与强震仪和结构健康监测数据融合。这些年不少建筑的屋顶、桥墩和土石坝顶部都装了强震加速度仪。当强震发生时这些结构台阵记录下来的不光是地面输入还包含结构自身的响应特征。如果能将这些数据与震级估算耦合起来就有可能在震后几分钟内直接给出某个特定工程场地遭受的实际地震动强度估计比单纯依赖“震级-烈度衰减关系”相对更贴近真实破坏状态。写在最后的几点经验震级这个概念看起来只是一个小数点带一位的数字但它从地震仪记录到的振幅一路走到灾害评估和抗震设计中间隔了量规函数、矩张量反演、概率危险性分析、场地效应校正等多层逻辑。我平时最怕看到的就是有人把震级当“唯一真理”直接套用——比如拿着某地历史上的最大震级直接断言“未来还会不会超过”或者“某个建筑是否安全”。正确的打开方式是把震级当成地震释放能量的标尺同时搭配震源深度、震源机制、破裂方向、场地条件和建筑结构特征才能做完整的风险判断。这些年处理数据的一个心得是永远要看震级的“后缀”区分ML、Ms、Mw再看测量台网和修正次数。只有理解了这套标度体系背后的物理和统计逻辑你在回答“这次地震到底有多大”时给出的才能不仅仅是一个数字而是一个真正经得起推敲的判断。