传统与现代震级标度_刘同辉.pdf

第 卷第期 年月地震工程学报 ，收稿日期：第一作者简介：刘同辉（），男，硕士研究生，主要从事震级和地震地质研究方面的工作。：。通信作者：王墩，男，教授，博士生导师，主要从事观测地震学、地震地质及大震实时减灾方面研究。：。刘同辉，王墩 传统与现代震级标度 地震工程学报，（）：，（）：传统与现代震级标度刘同辉，王墩，（中国地质大学（武汉）地球科学学院，湖北 武汉 ；中国地质大学（武汉）湖北巴东地质灾害国家野外科学观测研究站，湖北 武汉 ）摘要：首先阐述并比较传统震级（、及）所用数据、计算方法及适用范围，系统梳理两种现代震级标度和的理论基础及计算方法。然后根据使用数据及适用场景不同，介绍、及 等矩震级测定方法；同时，随着越来越丰富的 观测数据，基于高频 观测也越来越多地被用于地震震级快速测定工作。最后对这些震级测定方法的时效性进行对比分析，讨论不同震级标度的适用场景及其稳定性。关键词：震级测定；传统震级；现代震级；矩震级中图分类号：文献标志码：文章编号：（）：，（.，（），；.，（），）：，（，），（）（，），：；引言地震作为自然界中最严重的自然灾害之一，其震级大小一直都是人们关注的热点。震级是衡量震源产生地震波大小或强度的标度，由震级可以进一步求得其他震源参数，如地震能量、破裂持续时间、地震烈度分布等，所以震级是应急救援、地震危险性评估的重要依据。里克特（）在 年提出了地方性震级，也称里氏震级。该震级是里克特根据古登堡（）与和达清夫（）的建议得出的，这个震级是由伍德安德森（）地震仪上记录的水平方向上的最大振幅来决定的，受限于震中距范围和仪器精度，无法在全球范围内统一。为了测定远距离浅源地震的震级，年古登堡通过测量周期约为 面波的水平位移基础上，发展了面波震级，同时有地震学家提出了与地震能量更加相关的质点最大运动速度（）替换振幅。为了确定震源深度超过 的地震震级，古登堡和里克特通过计算波、波和波的振幅，得到体波震级。和都是根据相对短周期地震波的振幅计算出来的，通常和的周期小于，而的周期在 左右（也称为（）。随着数字地震仪的发展，有限频带震级和（）可以由宽频带震级（）和（）来补充。后两种震级都是通过直接量取垂向上地震波最大速度 而得到的，以此方便统一标准，进行地震速报。、和等传统震级标度本质上是经验性的，与地震的震源物理过程没有直接关系。同时，在不同震级范围内都存在震级饱和现象，难以准确测定大震震级。由于地震矩与地震能量密切相关，和 从地震矩角度出发，在 年提出了矩震级，将震级与地震破裂、断层作用等联系起来，它是与位错量的大小和破裂面积相关的力学量度，可以准确测定任意大小的地震 。矩震级根据使用波形及适用情形等不同，又可以分为、等多种标度。震相（）是介于波和波之间的长周期震相，其反演具有不易限幅、结果相对稳定、反演速度快等特点，在大震震后应急响应中得到广泛应用。主要通过对波的矩张量反演确定海上大地震的海啸潜力，后来扩展到一般的远震。是通过引入全球宽频带地震台记录的直达波最大位移与台阵技术得到的震源持续时间相结合，提出的一种新的震级标度，用以快速准确测定大震震级。能量震级则直接对应地震能量，能更好反映震源的动态特征，它与地震波速度、拐角频率等重要震源物理参数相关，对大地震的破坏性可以进行更好的描述。随着现代全球定位系统（，）观测技术的进步，利用高频 记录的峰值地动位 移（，）来测定地震震级的方法又为地震学提供新的思路，它对低频信号更敏感，且不会限幅，能够更真实记录地震的位移记录。图不同震级标度及所用数据 传统震级传统震级标度从整体上来看主要有两个共性：一是在同一震源引起的地震中，震级越大，其对应地震波的振幅越明显；二是地震波从震源扩散至观测位置的过程中衰减是已知的，都可用以下公式表示：（）（，）（）式中：表示震级；为最大振幅（最大地动位移）；为周期；（，）是用于对振幅随震中距和震源深度的变化作校正的因子；是台站校正因子，是震源校正因子，它们都是用来消除地壳结构变化、岩性差异等因素的影响。随着现代科学技术的进步，地震学观测的精度和要求越来越高，数字地震仪以其频带范围宽、测定时效高的特点为研究提供了更强有力的支持。年，国际地震学和地球内部物理学协会（，）利用宽频带地震记录提出了新的震级标度，并在我国首次成功应用，中国地震台网中心（，）也于 年在参照 新震级标度的基础上确定了国家震级标准（），因此以下使用最新国家震级标准公式。地方性震级（）测定地方性震级应使用仿真的 短周期地震仪，记录周期为 的波或 的最大振幅。公式如下：（）（）（）（）式中：为最大振幅；、分别为南北向和东西向的最大振幅；为震中距；（）为量规函数。地震工程学报 年主要用于测定震中距小于 的地震，其作为首个提出的地震震级，只考虑了最大位移的变化，未考虑到周期的影响，在量规函数（）中也没有考虑与震源深度关系，再加上震中距的限制，所以主要适用于区域性的浅源地震。面波震级面波震级测定面波震级应使用仿真的基式（）中长期地震仪，记录周期为 水平向面波质点运动的最大位移。公式如下：（）.（）式中：为水平方向上地动位移的矢量和；为周期；为震中距。主要用于测定 范围内的地震，由于水平向上是瑞利波和勒夫波叠加的混 合波，可 能导 致 震级 测定 结 果 不稳定，所以 选用仿真 （长周期）地震仪测定 左右的垂向面波。震级的震中距范围较大，但可观测到的震源深度较浅（），适用于测定浅源远震的震级。宽频带面波震级（）测定宽频带体波震级（）应使用宽频带地震仪，在垂向上记录周期为 面波的最大运动速度。公式如下：（）（）（）式中：为测得垂向面波质点运动的最大速度；为相应周期；为震中距。（）主要用于测定 ，且 范 围 内 的 地 震。相 较 于 颁布的面波震级（），宽频带震级（）有着更大的周期范围，对于区域性的浅源小地震、大地震都同样适用。体波震级短周期体波震级测定短周期体波震级应使用仿真的 短周期地震仪，记录周期为以内的 波在垂向上的最大运动位移。公式如下：（）（，）（）式中：为波的最大运动位移；为周期；为震中距；（，）为量规函数。主要用于测定 范围内的地震，由于只记录短周期体波，无法对大震震级进行准确测定，比较容易出现震级饱和的情况。同时，波震相基本不受震源深度的影响（），可以被清楚地识别，再通过量规函数对震中距和深度进行校正，可以快速测定中深源小震的震级。宽频带体波震级（）测定宽频带体波震级（）应使用宽频带地震仪，在垂向上记录周期为 体波的最大运动速度。公式如下：表传统震级的计算公式及其适用范围 震级标度公式仪器地震波类型记录方向周期震中距深度（）（）：仿真伍德安德森地震仪 ：（）（）仿真 短周期地震仪 ：（）（）仿真 （长周期）地震仪 ：（）仿真基式（）中长期地震仪（）（）宽频带面波记录仪 （）（，）仿真 短周期地震仪仿真 （短周期）地震仪（）（）（，）宽频带体波记录仪波或 波水平 波垂向 波和 波水平 波垂向 波垂向波垂向 震源距 第 卷 第期刘同辉，等：传统与现代震级标度（）（）（，）（）式中：为测得垂向体波质点运动的最大速度；为相应周期；为震中距；为深度。（）主要用于测定 范围内的地震，一般只需要 波的初始地震波信号，测定较大震级的时间较短，因此在地震应急、海啸预警中应用广泛。日本气象厅震级 日本是世界上地震、海啸灾害最频发的国家之一，居民的生命财产安全遭受着严重威胁。巨大海啸的产生很大程度上取决于地震的震级，海啸预警需要在短时间内准确测定震级，为应急减灾提供指导。日本气象厅（）依托于本土高密度的地震台网，使用周期小于 的体波精确计算出日本海岸 范围内局部地震事件的震源位置和震级，通过该方法确定的震级被称为 ，公式如下：（）（，）（）当地震事件满足震级大于，深度小于 ，并发生在海岸附近（）以内时，使用地震预警（）震级，公式为：（）（）式中：是通过积分加速度记录得出的最大位移；为震中距；为深度；为常数。现代震级地震能量和地震矩是目前对于震源物理特征描述最清楚的参量。地震波能量是用于量化震源的特征量，地震以地震波形式辐射的能量主要集中在震源频谱的拐角频率附近。对于地震危险性和风险评估，地震能量的大小更受关注，特别是对于其中的高频成分。地震矩是仅次于地震波能量的第二个量化震源的特征量，是描述地震大小的绝对力学量，它是源区不可恢复的非弹性变形的量度，这表明地震矩是断层滑动引起地震强度的直接量度。地震矩和地震波能量是两个重要的物理量，由地震矩可以得到矩震级，由 地 震 波 能 量可 以 得 到 能 量 震级。地震发生时绝大部分能量转化为机械能（断层位移和岩石破裂）并以热能形式散失，很小一部分以地震波能量（）形式辐射出来，被台站记录到。尽管地震波能量（）占比很小，但仍是定量化判别地震大小的重要手段，也是地震学研究的重要参数之一 。古登堡和里克特 根据震级与地震能量的关系得到经验公式：（）（）值得注意的是，利用震级与地震能量之间的经验关系来估计地震能量具有一定的局限性，仅仅是对地震能量的粗略估计。震级是针对单一频段地震波的测定结果，例如测定短周期体波震级一般使用周期是左右的体波信号，面波震级一般是使用 左右周期的面波信号，而地震能量的测定则应考虑所有频段的地震波能量。的求解地震矩最早是 年由日本学者安芸敬一（）首次提出，它既可以通过波长远大于震源尺度的地震波远场位移谱测定，也可以用近场地震波、地质与大地测量等资料测定。我们假设地震过程中应变能等于地震波能量，根据 等 的研究，地震矩和辐射能量之间存在一定的比例关系，可以从地震的长周期频段中近似得到。等 和 等 利用这种比例关系提出了一个简单的震源模型，得到了辐射能量和地震矩之间的关系。在地壳与地幔中，应力降 约为，可得到地震能量与地震矩关系式：（）（）对和都取平均值进行计算，可近似得到，针对大地震，我们可以测得断层面积、断层长度以及地震的平均位错量等静态震源参数，进而测得地震矩。公式如下：（）式中：为地壳与上地幔间存在介质的剪切模量，约为（）；为断层面的平均位错量；为断层面面积。地震矩是根据地震波形计算出来的。为计算地震矩，需要进行傅里叶变换，将位移波形从时间域转换到频率域。从图中可以看出，近场位移脉冲的频谱在低频时表现为水平的，等于阴影部分的面积。随着频率增大，位移频谱持续下降，角频率与脉冲宽度成反比。地震矩与频谱水平成正比，计算公式如下：地震工程学报 年 （）式中：是震源区岩石的平均密度；是平均波速（通常是波或波）；是震源距；为辐射模式项（地震波和波的振幅与方位角和离源角的关系。图近场位移与远场位移和速度之间的关系 等 于 年提出了一个完全由地震矩决定的、新的震级标度，利用式（）中与的关系代入到式（）中，用替换可得到公式：（）（）式中：单位为，由于矩震级是由地震矩计算出来的，所以不会出现震级饱和现象。地震矩是一个静态物理量，它是由地震波振幅低频分量的大小决定的，反映震源处的破裂强度。断层面积越大，激发的长周期地震波的能量越大。所以地震矩与断层长度、断层宽度、震源断裂平均位错等静态构造效应密切相关。矩震级根据不同的地震波类型和适用范围又划分为多种，包括、和 等标度（表）。表不同的震级分类 震级标度震级范围震中距地震波备注 相（）地震的结果很稳定，作为 首选发布震级 主要是长周期体波（）、中周期面波（）在没有稳定的情况下，发布此震级 体波（波 波）（）将其发布优先顺序排在和的情况之后 区域范围内的矩张量反演（）利用区域地震台网，对大于地震的结果可靠。基于高密度宽频带地震台网，可以计算小于地震的 波使用宽带地震仪垂直分量上记录的直达波，主要用来快速估计俯冲大震震级 波综合 台站垂直向记录的直达波最大位移和区域密集台阵叠加得到的破裂持续时间，得到大震地震矩。主要适用于级及以上大震、作为美国地质调查局（，简写为 ）三种主要的矩震级对外发布，为全球范围内大震的实时减灾和科学研究提供了重要参考。是一种通过波震相快速测定震级的标度，最初是由 等 开发的，目的是为了快速估计矩震级，从而准确评估日本近海大型地震的海啸潜力。后将其扩展至远震，现已成为太平 洋 海 啸 预 警 中 心（，）的主要测定震级。从此，就成为 快速估计的优选震级。其公式如下：（）（）（）（）是在日本建筑研究所 引入全球宽频带地震台（）记录直达 波最大位移的工作基础上，王墩等将密集台阵技术和反投影方法得到的震源持续时间（）相结合，提出一种新的震级标度，来快速精确测定大震震级。其公式如下：（）（）式中：是事件的震级；为事件在第个台站记录到的远震波的最大位移，震中距为；是参与计算的台站数量；为震源持续时间；、均为参数项。将该方法应用于 年间发生在日本周边的地震，得到的地震震级与美国地质调查局最终确定的矩震级（，）基本相当（标准偏差小于）。在时间效率上，该方法速度快于基于全球地震台网数据的 反演震第 卷 第期刘同辉，等：传统与现代震级标度级方法，与基于区域地震台网数据的 反演震级方法相当。由于远场直达波最大位移正比于最大地震矩变化率（），因此在假定地震震源持续时间函数为规则三角形或长方形的情况下，远场 波最大位移与破裂持续时间的乘积给出了地震的总地震矩（）（图）。通过地震矩，就得到了矩震级。该方法操作简单，不涉及反演，能实现自动、快速（）地获得稳定的震级估计。因此将新震级标度 与 反演结果相互验证、综合利用，将能更好地实现快速准确地测定大震震级。图 震级原理示意图 此方法主要利用两部分数据：（）全球台站的数据计算波在垂向位移的最大幅值；（）区域密集台网的数据，用反投影方法来计算地震的震源持续时间；（）将两者通过理论或经验公式结合，得出震级。的求解能量震级与一个明确的震源物理参数有关，即辐射的地震能量。地震辐射的能量以地震波的形式集中在震源频谱的角频率附近，这使得比更适合用来描述地震的破坏潜力 。是与震源谱的低频渐近线有关，描述了地震源的整体构造效应，而则是在震源频谱的更大频率范围内计算的，与地震工程研究中感兴趣的频率更有关系。弹性波释放的能量与地面运动速度的平方成正比，假设地震的震源为点源，震源周围为均匀球面，就可以用远震波的垂直向记录测定地震能量。其公式为：?（）（）式中：、和分别是波速度、波速度和震源处介质密度；、分别代表积分截取频率的最小值和最大值，代表地震矩率的导数（地震矩的二次导数）；为频率。经过乔伊和博特莱特等的研究，确定了新的的经验关系为 ：.（）将替换为，建立了能量震级公式：（.）（）随着数字地震学的发展，地震资料的观测精度不断提高，方法不断更新，一方面可以利用宽频带地震记录直接计算出地震能量，另一方面通过震源机制的校正和传播效应的修正可以建立新的理论公式。美国国家地震信息中心（，简写为 ）利用 和 在 年提出的方法，测定全球范围内中强地震的辐射能量，进而测得能量震级。但其程序需要对震源机制进行校正，因而测定能量震级的时间较长。德国国家地球科学研究中心（，简写为 ）也在积极探索能量震级的测定方法，等 提出一个新的程序，建立了新的辐射能量计算公式，使用发震后短时间内震中距 范围内的远震波来计算能量震级。考虑地震波从震源到接收器传播过程中所经历的频率衰减，通过应用不同频率的频谱振幅衰减函数，来实现对传播路径效应的修正。该程序避免了时间窗口饱和的问题，并通过研究证明辐射图样对辐射能量的测定结果影响较小，不需要根据震源机制解对结果进行修正，因此在时效性上有了明显提高。地震以地震波形式辐射的能量主要集中在震源谱的拐角频率附近，因此能量震级更适合描述地震的潜在破坏性。对于地震灾害与风险评估，人们更关注的是地震辐射能量的大小，特别是其高频成分的地震辐射能量。从断裂动力学、运动学以及地震所造成灾害的角度来说，地震能量与应力降、震源破裂速度密切相关，因此联合测定和对于量化和评估地震或海啸造成的灾害评估具有十分重要的意义。基于高频 观测的地震震级随着应急地震学的发展，开发可靠的地震预警地震工程学报 年（，）方法，量化地面运动和地震源参数之间的关系显得尤为关键。在时间域中，典型的地面运动指标包括地面加速度峰值、有效地面加速度峰值、地面速度峰值和地面位移峰值（分别为、和 ）。频率域中的指标包括主要周期 以及频谱加速度、速度和位移。不同震级和距离范围的地震引起的地面运动记录被合成为一系列参数，而后这些参数被应用于地震学的研究中。在传统的震级测定中，地震仪的传感器是惯性框架下产生的直接观测值，因此对微弱地震动信号敏感。但当大地震发生时，基于惯性原理设计的地震仪难以记录地面强震动，测得的地震动值常出现限幅现象。而全球定位系统（）的位置和位移均是在地球参考框架下计算得出的，不会受到基线偏移的影响，可以直接测量地表强震动位移。尽管 测量的噪声基底限制了其对于地震事件的敏感度，但在一定范围内地面位移峰值与震级还是有着良好的相关性，并且 观测记录不会出现震级饱和现象，可以用来测定大震震级。实时 观测比使用波初动的地震预警方法耗时长，但得到的震级更 为 准确，可以及 时评 估震 害 和 海啸 预警。等 提出了如下的峰值地动位移统计关系，描述了不同震级在近、中、远场的地震动衰减关系。其公式为：（）（）式中：、和是回归系数；是震源距离。（峰值地动位移的简称）是未经滤波的 地震图上的动态位移峰值，单位为。利用高频 记录的峰值地动位移（）来测定地震震级的方法可以为地震学提供新的思路，它对低频信号更敏感，不会限幅，能够更真实的记录地震的位移。讨论和结论震级饱和利用观测到的地震波振幅确定震级时，由于不同震级标度测定的是特定频段的地震波振幅，当震级大到一定的级别时，测得的最大振幅不再增加，致使测得的地震震级不再随地震的增大而增大的现象，称为震级饱和（）。传统的震级标度，如地方性震级、体波震级、长周期体波震级及面波震级等，不一定能反映出大地震的真实规模，因为它们来自对大地震的研究，由相对较短的周期振幅（约）得出的，超过一定限度就会饱和（：，：，：；见图）。对于大地震，单一频带的振幅不能代表整个断层的特征，地震破裂过程涉及一个广泛的周期（或频率）范围，至少从到，振幅范围很宽，大约从 到 。较大的地震在较长的周期（较低的频率）范围内会释放出更多的能量。图不同震级标度的震级饱和现象 一般来讲，在较低的频率范围，地震越大，地震矩随之增大。对级以下地震，不同震级对应的振幅值反映了地震矩。随着地震震级的增加，越来越难测得真实的地震矩。如图所示，只有当使用频率低至 的地震信号时，才可能正确测定级地震的真实大小。根据源模型预测的波频谱，假设恒定的应力降为 ，低频的振幅等于地震矩，圆圈表示角频率图预测的波频谱 震级求解的时效性近些 年 来，无 论 是 年 汶 川 大 地 震 还 是 年日本东北大地震等破坏性大震，都体现出震第 卷 第期刘同辉，等：传统与现代震级标度级和烈度评估不准而导致公民的生命和财产遭受了巨大损失。因此，加强建筑物抗震强度，重视应急地震学研究迫在眉睫。应急地震学的重点研究内容之一就是利用波形信息快速确定震级等震源参数。基于准确地震震级，结合震灾区的强地面震动记录，就可以快速确定地震的破裂过程，科学评估各地区可能的受灾程度 ，为抗震减灾提供精准指导。所以震级测定的时效性，显得尤为重要。传统的地方性震级、面波震级、体波震级等在震后数分钟内由中国地震台网中心测定得到，根据不同情况进行速报：对于的浅源地震，一般选用地方性震级进行发布；对于的浅源地震，一般选择面波震级进行发布；对于中深源地震，一般选择短周期体波震级进行发布。借助于数字地震仪的进步，宽频带面波震级（）和宽频带体波震级（）可以在原始的速度宽频带记录仪上直接测定，这给我国地震速报工作提供了很大帮助。而矩震级和能量震级由于需要计算地震矩和地震能量，测定方法较为复杂，在震后快速测定过程中往往存在因观测数据不足而存在不稳定性。因此，针对大震，采用和 等多种震级标度相互验证，可更精准有效地判定大震震级，为大震应急救援与海啸预警提供基础科学数据。参考文献（）刘瑞丰震级测定技术的发展与展望城市与减灾，（）：，（）：陈运泰，刘瑞丰地震的震级地震地磁观测与研究，（）：，（）：，（）：，（）：，（）：，（）：，：，：，：，（）：，（）：，（）：陈运泰，刘瑞丰矩震级及其计算地震地磁观测与研究，（）：，（）：，：，（）：，（）：，（）：袁乃荣，刘瑞丰，李赞，等能量震级及其测定地震地磁观测与研究，（）：，（）：陈锋，杨建思，王伟平，等 基于高频 峰值地动位移的震级标度探讨中国地震，（）：，（）：，（）：国家质量监督检验检疫总局，中国国家标准化管理委员会 地震震级的规定：北京：中国标准出版社，：，地震工程学报 年 ，：，：，：；，（）：，（）：，（）：，：，（）：，（）：，（）：，：，：，（）：，（）：，（）：，（）：，：，（）：，（），（）：，（），（）：，（），（）：，（）：，（）：李赞能量震级的测定及其应用研究北京：中国地震局地球物理研究所，：，：，：，：，（）：，（）：，（）：，：，倪四道 应急地震学的研究进展 中国科学院院刊，（）：，（）：，（）：，（），（）：（本文编辑：张向红）第 卷 第期刘同辉，等：传统与现代震级标度