从震源破裂过程角度分析中国大陆强震面波震级高于矩震级的原因.pdf

407 国家自然科学基金(42021003,42074058)资助 收稿日期:20220524;修回日期:20220623 北京大学学报(自然科学版)第 59 卷 第 3 期 2023 年 5 月 Acta Scientiarum Naturalium Universitatis Pekinensis,Vol.59,No.3(May 2023)doi:10.13209/j.0479-8023.2023.020 从震源破裂过程角度分析中国大陆强震 面波震级高于矩震级的原因 许月怡 张勇 北京大学地球与空间科学学院,北京 100871;通信作者,E-mail:zhang- 摘要 基于有限断层模型和数值模拟,研究震源参数对面波震级的影响。结果表明,在矩震级相同的情况下,地震破裂持续的时间越短,面波震级就越高,反之亦然。影响地震持续时间的震源参数包括子断层破裂持续时间、破裂速度和破裂长度。若地震矩不变,子断层破裂持续时间和破裂长度与面波震级负相关,破裂速度与面波震级正相关。中国大陆测定的面波震级多数情况下显著高于矩震级的原因可能有二:一是中国大陆发生的强震多为板内地震,由于应力降偏大,破裂持续时间和断层长度可能偏短;二是相对于美国国家地震信息中心测定的 Ms(20),中国地震台网中心测定的面波震级系统性地偏大。关键词 面波震级;矩震级;震源参数 Analysis of the Reasons why Surface Wave Magnitude is Higher than Moment Magnitude in the Mainland of China from the Perspective of Source Rupture Parameters XU Yueyi,ZHANG Yong School of Earth and Space Sciences,Peking University,Beijing 100871;Corresponding author,E-mail:zhang- Abstract The influence of the source parameters on Ms is studied through numerical simulations based on finite-fault models.The results show that the Ms is negatively correlated with earthquake rupture duration when moment magnitude remains unchanged.The source parameters that can impact the rupture duration contain the sub-fault duration,rupture velocity,and rupture length.Among them,the sub-fault duration and rupture length is negatively correlated with Ms,while the rupture velocity is positively correlated with Ms,if the moment magnitude remains unchanged.The fact that the measured Ms of strong earthquakes in the mainland of China are significantly larger than moment magnitudes may be due to two reasons.First,the strong earthquakes in the mainland of China are mostly intraplate earthquakes with relatively higher stress drop values,which may have short rupture duration and rupture length.Secondly,Ms determined by China Earthquake Networks Center(CENC)is systematically larger than Ms(20)determined by National Earthquake Information Center.Key words surface wave magnitude;moment magnitude;source parameters面波震级(Ms)和矩震级(Mw)是衡量强震(震级M 6.0)规模的两种常见标度,二者各有优缺点。面波震级的测定相对简单,但当 M 8.5 时,容易发生饱和;矩震级不会饱和,但与地震波能量的关系较弱,某些情况下无法较好地反映震感和灾害的程度。对一般地震而言,不同震级标度之间整体上存在线性经验关系13,但具体到某一个地震事件,中国地震台网中心(China Earthquake Networks Center,CENC)测定的面波震级(Ms)、美国国家地震信息中心(National Earthquake Information Center,NEIC)测北京大学学报(自然科学版)第 59 卷 第 3 期 2023 年 5 月 408 定的 Ms(20)以及全球矩心矩张量组(Global Centroid-Moment-Tensor,GCMT)测定的矩震级 Mw 之间存在较大的差异(图 1)。例如,2014 年云南景谷地震的Ms 比 Mw 大 0.8 级,2017 年西藏米林地震的 Ms 和 Ms(20)比 Mw大 0.4 级(表 1)。震级测定的一致性问题一直备受关注。矩震级完全取决于地震矩,与断层附近介质的剪切模量、错动面积和位错量相关,对断层错动的时间过程不敏感4;面波震级则由面波振幅与周期的比值决定,Ms 和 Ms(20)采用的面波周期范围也有差别4,所以受断层错动时间过程影响。一个较为简单且极端的例子是,假定两次矩震级相同的地震事件具有相同的位错量,但其中一次事件完成这一错动的时间为10 s,另一次事件完成这一错动的时间为 100 s,由于用于测定面波震级的远场地震波幅度与断层错动速率成正比,前者的面波幅度显著高于后者,导致其面波震级明显偏大。本文基于有限断层模型,采用不同的震源参数 Ms,Ms(20)和 Mw 分别来源于中国地震台网中心正式目录(http:/)、美国国家地震信息中心(https:/www.usgs.gov/)和全球矩心矩张量项目(https:/www.globalcmt.org/)图 1 2009 年 1 月至 2021 年 12 月中国大陆 Ms 6.5 的地震Ms与Ms(20)震级对比 Fig.1 Comparison between Ms and Ms(20)of earthquakes Ms 6.5 from January 2009 to December 2021 in the mainland of China 表 1 2009 年 1 月至 2021 年 12 月中国大陆Ms 6.5 地震目录 Table 1 Earthquake catalogue(Ms 6.5)from January 2009 to December 2021 in the mainland of China 发震时间(北京时)经度/()纬度/()震中 位置 Mw Ms(20)Ms 发震时间(北京时)经度/()纬度/()震中 位置 Mw Ms(20)Ms 20210522 02:04:13 98.37 34.61 青海 玛多 7.4 7.6 7.9 20150703 09:07:44 78.15 37.56 新疆 皮山 6.4 6.6 6.5 20210521 21:48:35 99.88 25.70 云南 漾濞 6.1 6.2 6.5 20141007 21:49:39 100.55 23.40 云南 景谷 6.1 6.3 6.9 20200723 04:07:18 86.81 33.19 西藏 尼玛 6.4 6.3 6.6 20140803 16:30:12 103.33 27.11 云南 鲁甸 6.2 6.3 6.6 20200626 05:05:17 82.34 35.69 新疆 于田 6.3 6.2 6.5 20140212 17:19:48 82.51 36.14 新疆 于田 6.9 7.2 7.3 20200119 21:27:56 77.21 39.83 新疆 伽师 6.0 6.1 6.5 20130722 07:45:56 104.21 34.54 甘肃 岷县 6.0 6.2 6.7 20171118 06:34:18 95.00 29.75 西藏 米林 6.5 6.9 6.9 20130420 08:02:47 102.99 30.30 四川 芦山 6.6 6.8 7.0 20170809 07:27:52 82.89 44.27 新疆 精河 6.3 6.4 6.6 20120630 05:07:31 84.74 43.42 新疆 新源 6.3 6.4 6.6 20170808 21:19:48 103.82 33.20 四川 九寨沟 6.5 6.8 7.0 20100414 07:49:36 96.59 33.22 青海 玉树 6.9 7.0 7.3 20161125 22:24:30 74.10 39.20 新疆 阿克陶 6.6 6.7 6.8 20090828 09:52:06 95.90 37.60 青海 海西 6.3 6.2 6.6 许月怡等 从震源破裂过程角度分析中国大陆强震面波震级高于矩震级的原因 409 模拟地震波,并基于模拟的地震波测定其面波震级,详细地探究不同的震源参数(包括子断层破裂(断层上质点振动)持续时间、破裂速度和破裂长度)如何造成面波震级与矩震级的差异。1 研究方法 1.1 基于有限断层模型的地震波模拟 根据震源表示定理,地震辐射的波场可以表示如下5:,(,)*(,;)d(),npqnp qutmtGt xx(1)其中,un(x,t)为观测点(台站)x 处 n 方向的位移,Gnp,q为格林函数,表示断层面 处单位强度的力偶在观测点 x 处激发的 n 方向位移。mpq 是断层面 处地震矩张量的分量,如果断层面上的震源机制不随时间改变,则可表示为 mpq(,t)=M0()s(,t)(ep vq+eq vp),e 和 v 分别为滑动方向和断层面法向方向的单位矢量,M0()为标量地震矩,s(,t)是震源时间函数。假定断层面由 K 个大小相等的子断层组成,台站总数为 L,根据式(1),台站 l 处的位移()lnu t可写成离散形式6:01,*(,;)()Klnkkknp qlkpqqputMstGt。xe ve v(2)在震源机制已知的情况下,本文基于 AK135 大 陆模型78,利用反透射系数法9计算 Gnp,q(x,t;)(ep vq+eq vp),所以只需给定断层面上的滑动分布和破裂模式,就可以根据式(2)计算台站处的观测波形。测定面波震级时,需要将观测资料转换为仿真资料,假设()lnu s为地面位移()lnu t的频谱(s=i2f,i=1),模拟地震仪的传递函数为 H(s),仿真记录的位移()lnd t可表示为 (1()e d2 illstnndtUsH ss。(3)中国基式(SK)中长周期地震仪的传递函数10为 32252.36()(0.44720.2693)(52.3625.),72sH sssss(4)美国世界标准台网长周期(WWSSN-LP)地震仪的传递函数11为 3220.97866()(0.80360.16877)(0.136570.00427)sH sssss。(5)两种仪器的幅频特性曲线(图 2)表明,SK 地震仪对 0.115 s 周期的地震波响应较好,其频带在0.211 s 是平坦的;WWSSN-LP 对 350 s 周期的地震波反映较好,频谱峰值出现在 15 s4。考虑到中国地震台网与境内地震的平均震中距小于 25,而Ms(20)需要震中距大于 20,我们分别设定震中距 图 2 中国基式(SK)中长周期地震仪和世界标准台网长周期(WWSSN-LP)地震仪的幅频特性曲线 Fig.2 Displacement response of typical Chinese Intermediate Period Broad Band SK Seismograph and World-Wide Standardized Seismograph Network Long-period Seismograph 图 3 测定Ms和Ms(20)所用台站分布 Fig.3 Distribution of stations for the determination of Ms and Ms(20)北京大学学报(自然科学版)第 59 卷 第 3 期 2023 年 5 月 410 在 525和 3090范围内均匀分布的各 100 个观测台站来测定 Ms 和 Ms(20),如图 3 所示。为了探究矩震级固定情况下震源参数对面波震级的影响,分别模拟矩震级为 6.8,7.2 和 7.5 的地震,包括单侧和双侧破裂模型(即表 2 中 6 种震例),震源机制固定,断层的走向、倾角和滑动角分别为 90,45和 45。图 4 和 5 以 Mw=7.2 为例,展示震源参数的定义。破裂长度为 L,矩震级为 6.8,7.2 和 7.5的破裂模型分别采用 1 km1 km,2 km2 km 和 3 km3 km 的子断层,震源时间函数均用一个简单的等腰三角形表示,破裂起始时间用子断层到震源的距离 l 除以破裂速度 Vr 来计算,子断层的破裂持续持续时间为 Tr(图 5(a),故整个地震的破裂持续时间为 Ta=lmax/Vr+Tr(lmax 与 L 正相关),双侧破裂(图5(c)持续时间少于单侧破裂(图 5(b)。以单侧破裂Ms(20)=Mw 时的震源参数为参照,6 种模拟震例都 表 2 破裂模型的震源参数 Table 2 Source parameters of rupture models Mw 破裂模式 子断层大小 破裂 宽度/km 参考模型 参考模型面波震级 震源 深度/km 破裂模型 Tr/s Vr/(kms1)L/kmMs(20)Ms Tr/s Vr/(kms1)L/km 6.8 单侧 1 km1 km 13 10 3 30 6.8 7.2 5.3 230 0.55.8 680 6.8 双侧 1 km1 km 13 10 3 30 7.0 7.4 5.3 230 0.55.8 680 7.2 单侧 2 km2 km 26 10 3 60 7.2 7.4 10.6 230 0.55.8 121607.2 双侧 2 km2 km 26 10 3 60 7.3 7.8 10.6 230 0.55.8 121607.5 单侧 3 km3 km 39 10 3 90 7.5 7.5 15.9 230 0.55.8 182407.5 双侧 3 km3 km 39 10 3 90 7.8 8.0 15.9 230 0.55.8 18240 图 4 参考模型的断层面滑动分布(Mw=7.2)Fig.4 Fault slip distribution of the reference model(Mw=7.2)(a)子断层震源时间函数示意图;(b)和(c)参考破裂模型的地震矩率函数(Mw=7.2,L=60 km,Vr=3 km/s,子断层破裂持续时间为 10 s)图 5 输入的震源时间函数 Fig.5 The input source time function 许月怡等 从震源破裂过程角度分析中国大陆强震面波震级高于矩震级的原因 411 在参考模型(表 2)的基础上,分别改变 L,Vr 和 Tr,根据式(2)(5)合成台站处的地震观测记录,模拟的观测数据采样率设定为 1 sps。1.2 震级公式 目前,CENC 测定浅源地震(深度小于 60 km)水平向面波震级时,需要将宽频带数字地震记录仿真为 SK 地震仪的模拟记录,然后根据下式10来计算:Ms=lg(AH/TH)max+1.66 lg()+3.5,2 130,3 s TH 25 s,(6)式中,(AH/TH)max 代表求取同一时刻或周期相差在1/8 周期之内的两水平分量面波位移的矢量和的模AH(单位:m)除以对应周期 TH(单位:s)的最大值;TH=(TN AN+TE AE)/(AN+AE),为东西和南北分量面波周期加权平均;是震中距(单位:)。NEIC 测定的 Ms(20)则根据 WWSSN-LP 地震仪的仿真记录,采用垂直向面波资料,然后根据下式11来计算:Ms(20)=lg(AZ/TZ)+1.66 lg()+3.3,20 160,18 s TZ 30 s)超出计算 Ms 和 Ms(20)所使用面波的周期范围,即此时利用式(6)和(7)计算的面波震级不能正确地反映面波的能量。2)随着 Vr 增加,面波震级增大,其变化速率从陡变缓,最后逐渐接近将整个断层面视为点源的计算结果。这是因为当 Vr 较小时,子断层的波场到达台站的时间相差较大,波形相互错开,导致叠加的面波幅度较小;当 Vr 较大时,所有子断层几乎同时破裂,加上震中距远大于破裂尺度,子断层波场叠加的结果与点源激发的波场相似。3)随着 L 增加,面波震级减小。原因在于破裂长度越大,断层面的平均滑动量越小,且相隔较远的子断层的波场到达同一个台站的时间差增加,导致面波振幅减小。4)震源参数相同时,双侧破裂的面波震级大于单侧破裂。这是因为可将双侧破裂视为两个破裂长度减半的单侧破裂的叠加,由于破裂长度减小,面波震级增加,因此双侧破裂的面波震级较高。总体而言,Mw 不变的情况下,如果地震的总破裂持续时间(Ta)减少(Vr 增大,或 L 和 Tr 减小),则面波震级增大;如果 Ta 增加(Vr减小,或 L 和 Tr 增大),则面波震级减小(图 7(a)和(b)。通过对比可以发现,Ms 比 Ms(20)平均大 0.23 级,存在系统性偏差(图 8)。Ms整体偏大的主要原因是CENC 与 NEIC 使用的量规函数不同,式(6)采用的数字 3.5 大于式(7)采用的数字 3.3。另一方面,图 7(c)显示 Ms Ms(20)差值与 lg(Ta)明显负相关,正交回归关系式为 Ms Ms(20)=1.13 lg(Ta)+1.86,相关系数为0.73。Ta 30 s 时,辐射的面波优势周期较小,此时 Ms 比 Ms(20)大 0.20.6 级(图 7(c),中国大陆强震 67 级居多(表 1),而 67 级地震的平均持续时间为 625 s13,这也可能是部分震例中 CENC测定的面波震级比 NEIC 显著偏大的原因之一。模拟结果表明,矩震级相同的地震,其面波震级差别可高达 1 级(图 6),这与不同震源参数下地震辐射能量不同有关。Vassiliou 等14给出地震辐射能E 的估算公式:2205511()d,1510EMTtt(10)其中,和 分别为 P 波和 S 波的波速,依据 AK135模型,=5.80 km/s,=3.46 km/s;T(t)为整个地震的归一化震源时间函数。我们通过(10)式来探究面波震级与辐射能 E 之间的关系,结果如图 9 所示,面 北京大学学报(自然科学版)第 59 卷 第 3 期 2023 年 5 月 412 从上往下,Mw分别为 6.8,7.2 和 7.5;实线和虚线分别表示单侧破裂和双侧破裂;黑线和灰线分别表示 Ms和 Ms(20)图 6 面波震级随震源参数的变化曲线 Fig.6 Surface-wave magnitude variation curves with source parameters 波震级与辐射能的对数 lgE 线性正相关。通过最小二乘拟合,得到两者之间的转换公式为 lgE=1.8678Ms(20)+0.9251,(11)lgE=1.7856Ms+1.1089,(12)相关系数分别为 0.96 和 0.88。Choy 等15计算了 NEIC 记录的 397 个实际震例的辐射能,推导出 E 与 Ms(20)的经验关系式为 lgE=1.5Ms(20)+4.4(图 9 中虚线)。式(11)估算的辐射能比 Choy 等15的结果小 0.51 个数量级(6.5 Ms(20)8),这是可以理解的。一方 面,本研究中波形的采样率为 1 sps,导致计算辐射能时的最高频率为 0.5 Hz;另一方面,本文假定子断层震源时间函数为一个简单的三角形,震源谱含较少的高频成分。实际上,地震波的高频部分(大于拐角频率)携带了超过 80%的能量16,所以式(11)和(12)中的 E 表征地震辐射的低频能量,与优势周期较长的“单色”面波震级17具有极强的相关性,而宽频带体波震级能够更好地反映地震辐射的高频能 量18。总而言之,中国大陆强震面波震级 Ms 高于国际主流机构测定的矩震级 Mw 和面波震级 Ms(20),这与震级计算公式和仿真仪器的差异,以及不同震 许月怡等 从震源破裂过程角度分析中国大陆强震面波震级高于矩震级的原因 413 图 7 Ms,Ms(20)以及两者震级差与地震总破裂持续时间(Ta)的关系 Fig.7 Relation among Ms,Ms(20),Ms Ms(20)and total rupture duration Ta 图 8 Ms与Ms(20)震级差数量分布 Fig.8 Histogram of the difference Ms Ms(20)级标度所示的物理意义不同有关19。为了与国际接轨,依据我国 2017 年发布的震级测定新标准地震震级的规定(GB177402017)10,应将 Mw 作为重点测定的震级,在及时测定地震矩的前提下,优先选择 Mw 对外发布19。另外,针对近震震级大于4.5 的浅源地震,应使用宽频带面波震级 Ms(BB)取代 Ms 作为对外发布的面波震级19。Ms(BB)是在垂直向的原始速度记录上测算的,量规函数与 Ms(20)相同,平均偏差为 0.054,并且 Ms(BB)的优势周期为 360 s,比 Ms 和 Ms(20)使用的频带更宽,能更好地衡量地震的强度19。需要注意的是,无论是面波震级还是矩震级,都存在“单色性”的局限,无法全面地反映地震激发的地震波能量大小。为此,未来首选的震级标度应当是快速且准确地确定地震波的能量震级,以便更好地对应地震造成的震感和破坏程度。3 结论 本文通过断层的数值模拟和理论格林函数计算 实线为本文模拟的面波震级与辐射能的正交最小二乘拟合线,虚线为 Choy 等15根据实际震例推导的经验关系 lgE=1.5Ms(20)+4.4 图 9 地震波辐射能与面波震级的关系 Fig.9 Relation between radiated seismic energy and surface-wave magnitude 北京大学学报(自然科学版)第 59 卷 第 3 期 2023 年 5 月 414 面波波场,研究矩震级不变情况下,子断层破裂持续时间(Tr)、破裂速度(Vr)和破裂长度(L)这 3 个震源参数对面波震级的影响,得到以下结论。1)总体而言,这 3 个参数都是通过改变地震的总破裂时长来影响面波震级,总破裂时间缩短(Tr 和L 减小,或 Vr 增大)时,面波震级增大,反之亦然。2)由于计算公式和仿真仪器的不同,Ms 总体上比 Ms(20)偏高 0.23 级,当地震持续时间较短时,Ms 显著大于 Ms(20),最大差值可达 0.6 级。3)矩震级对破裂模式和过程不敏感,面波震级强烈受震源谱影响,与地震辐射的低频能量有很强的相关性。不同的震级标度以各自的优势周期表征了地震波的基本特征,反映了地震波信号在该优势周期的强弱程度,都有其特定的意义和价值,彼此不宜强行换算和比较。我国测定的强震面波震级客观上存在比矩震级偏大的情况,这其中既有中国大陆强震特定震源特征的原因,也与我国特定的面波震级测定公式相关。一方面,中国大陆多数强震为板内地震,可能存在较高的应力降以及较短的破裂持续时间和破裂长度,导致面波震级相对于矩震级偏高;另一方面,我国测定面波震级公式的量规函数比 NEIC 大 0.2。上述两大因素共同导致我国测定的强震面波震级系统性地高于矩震级。参考文献 1 Scordilis E M.Empirical global relations converting Ms and mb to moment magnitude.Journal of Seis-mology,2006,10(2):225236 2 翟璐媛,张立文,任克新,等.中国台网数字化后与美国地震台网面波震级对比分析.地震地磁观测与研究,2015,36(2):1925 3 沙海军,吕悦军.中国地震台网面波震级与矩震级统计关系.地震地磁观测与研究,2018,39(6):3136 4 刘瑞丰,陈运泰.震级的测定.北京:地震出版社,2015 5 Aki K,Richards P G.Quantitative seismology.2nd edition.Sausalito:University Science Books,2002 6 Zhang Y,Feng W P,Chen Y T,et al.The 2009 LAquila Mw 6.3 earthquake:a new technique to lo-cate the hypocentre in the joint inversion of earth-quake rupture process.Geophysical Journal Interna-tional,2012,191(3):14171426 7 Kennett B L N,Engdahl E R,Buland R.Constraints on seismic velocities in the Earth from traveltimes.Geophysical Journal International,1995,122(1):108124 8 Montagner J P,Kennett B L N.How to reconcile body-wave and normal-mode reference Earth models.Geophysical Journal International,1996,125(1):229248 9 Wang R J.A simple orthonormalization method for stable and efficient computation of Greens functions.Bulletin of the Seismological Society of America,1999,89(3):733741 10 刘瑞丰,陈运泰,许绍燮,等.GB 177402017 地震震级的规定.北京:中国标准出版社,2017 11 IASPEI.Summary of magnitude working group reco-mmendations on standard procedures for determining earthquake magnitudes from digital data EB/OL.(20130327)20220622.http:/www.iaspei.org/commissions/commission-on-seismological-observation-and-interpretation/Summary_WG_recommendations_ 20130327.pdf 12 Kanamori H.The energy release in great earthquakes.Journal of Geophysical Research,1977,82(20):2981 2987 13 Bormann P,Wendt S,DiGiacomo D.Seismic sources and source parameters/New manual of seismolo-gical observatory practice 2(NMSOP2).Potsdam:Deutsches GeoForschungsZentrum,2013:1259 14 Vassiliou M S,Kanamori H.The energy release in earthquakes.Bulletin of the Seismological Society of America,1982,72(2):371387 15 Choy G L,Boatwright J L.Global patterns of radia-ted seismic energy and apparent stress.Journal of Geophysical Research:Solid Earth,1995,100(B9):1820518228 16 Ide S,Beroza G C.Does apparent stress vary with earthquake size?.Geophysical Research Letters,2001,28(17):33493352 17 陈运泰,刘瑞丰.地震的震级.地震地磁观测与研究,2004,25(6):112 18 宋潇潇,吴忠良,蒋长胜,等.矩震级及震源机制相似地震的辐射能量差别2014 年鲁甸地震与景谷地震宽频带体波震级的比较.中国地震,2016,32(4):653662 19 刘瑞丰,陈运泰,王丽艳,等.新的震级国家标 准的技术要点与主要特点.地震地磁观测与研究,2018,39(1):111