由计算机软件自动测定地震参数并完成发布已经成为全球地震台网的普遍做法,中国地震台网也发展了自己的地震参数自动测定软件系统,实现了对中国及全球地震时空强参数的自动测定,并自2013年4月1日起正式向公众提供准实时地震信息服务.相对于传统的人工测定地震参数流程,由计算机软件自动测定地震参数时效优势显著,其中对全球地震通常在震后10 min左右完成地震时空强基本参数测定,并通过多种手段实时向政府、媒体和社会公众发布.我国作为地球科学研究和防震减灾大国,对于全球地震信息的准实时服务,有助于提高我国地震科学、地球科学在国际地学界的地位和声望,扩大我国在国际防震减灾领域的影响,为国家争取更多的外交主动.
随着地震监测技术的发展和国家对防震减灾技术要求的日益提高,矩震级、震源机制(胡幸平等,2008;吕坚等,2011;张广伟等,2014)和地震破裂过程(杜海林等,2009;张勇等, 2009, 2013;王为民等,2013;刘成利等,2014;刘志鹏和盖增喜,2015)等新参数不断地丰富了传统地震时空强参数的测量.在地震参数自动测定中,自动测定震级是关键技术之一.地震灾害与震级的大小密切相关,震级的大小直接关系地震灾害评估和应急处置.近年来,国内外学者提出了一系列快速测定震级的方法(Tsuboi et al,1995, 1999;Wu, Teng,2002;Allen, Kanamori,2003;Wu et al,2006, 2007;Kanjo et al,2006;Zollo et al,2006;Hara,2007;Lomax et al,2007;Bormann, Saul,2007;Yamada, Ide,2008;Yamada, Mori,2009;Lomax, Michelini,2009;Delouis et al,2009;Allen, Ziv,2011;Earthwormcentral,2011;Vallée et al,2011;Katsumata et al,2013;Seiscomp3,2014),这些方法在地震预警等实时自动测定地震参数的应用中发挥着重要作用.其中MWP震级测定方法(Tsuboi et al,1995, 1999)以其快速稳定获得矩震级初步估计的优势,在海啸预警、地震参数准实时测定等工作中得到了广泛应用.
通过对国际国内震级测定方法的研究,以及吸取2013年芦山MS 7.0地震自动测定震级偏差较大经验的基础上(梁建宏等,2015),中国地震台网中心研发的地震参数自动测定软件系统(以下简称自动系统)对震级测定方法进行了改进,改进之一是使用MWP震级测定方法测定中强地震的震级.本文应用MWP震级测定方法对全球中强地震震级进行自动测定,结果显示,使用该方法可在震后快速稳定地获得地震事件矩震级MW的初步估计,改善自动测定中强地震震级的精准度,提升中国地震台网对全球地震信息的准实时服务能力.
1 MWP震级测定方法据Tsuboi (1995, 1999)的研究,宽频带垂直向P波位移可视为近似的震源时间函数,地震矩可以表示为P波位移的积分,由P波位移的积分可以获得每个台站的标量地震矩.标量地震矩表示为
(1) |
式中M0为标量地震矩,p1和p2分别为位移积分后的地震图上第一个峰值和第二个峰值的大小,ρ和α分别为P波传播路径上介质的平均密度和P波平均速度,r为震中距,FP为双偶极子震源机制确定的辐射图型.使用矩震级计算公式(Kanamori,1977)
(2) |
可得矩震级MW.根据Tsuboi等(1995)研究,对辐射图型进行校正,等价于平均矩震级加上0.2.为了快速得到矩震级,在计算单台矩震级时不作辐射图型校正,即忽略FP项,计算出每个台站的MW后,将所有台站MW的平均值加上0.2,对辐射图型进行校正.此震级被称为MWP (Tsuboi et al,1995),即宽频带P波矩震级.
由于MWP震级测定方法只使用P波,因此该方法可使用S波或面波限幅记录,靠近震中台站资料可供使用.另外,测定震级选取的时间窗位于S波之前,测定速度更快.因此,MWP方法得到了广泛应用,现已经成为太平洋海啸预警中心(PTWC)、美国国家地震信息中心(NEIC)等地震自动处理系统估计震级的方法.
本文在计算单台MWP时,按照Tsuboi (1995, 1999)的做法,不去除仪器响应,直接对宽频带垂直向记录进行两次积分得到位移积分记录,采用与Tezel和Yanik (2013)类似的做法,将定位时自动拾取的初至P波到时作为TP,将初至S波理论到时作为TS,在时间窗[TP,TS-3]内,寻找第一个峰值和第二个峰值p1和p2,传播路径上介质的平均密度ρ取为3400 kg/m3,使用与震中距有关的P波视速度α=0.16Δ+7.9 km/s (Kanjo et al,2006),根据式(1)和式(2)以及MWP_corrected=(MWP_initial -1.03)/0.843(Whitmore et al,2002)得到校正过的每个台站的MWP.例如2015年9月17日智利MW 8.2地震,震中距为42.1°的全球地震台网RCBR台原始宽频带记录、所取的P波窗口(图 1a中蓝色窗内的波形)、所取的P波窗口内原始宽频带记录、位移记录、位移积分绝对值记录及MWP随时间的变化见图 1b,在图 1b中红色竖线位置处取得位移积分绝对值的最大值和相应的MWP最大值,可得该台MWP=8.2.
由于准实时自动测定地震参数对时效性要求较高,故发布地震参数时,可能有些台站的实时波形还未满足时间窗[TP,TS-3],此时则用该时间窗内位移积分记录的最大值表示p1,计算MWP.
2 MWP震级测定方法的应用基于中国地震台网和全球地震台网共享的实时波形数据,自动系统实现了对中国及全球地震时空强参数的自动测定,数据流程见图 2,其中calmag模块根据locator模块自动定位结果和波形数据缓冲区中的实时波形数据测定震级,其使用的其中一种震级测定方法就是MWP震级测定方法.对于全球中强地震,自动系统使用全球地震台网共享的415个台站和中国地震台网47个台站的实时波形数据进行自动定位和测定震级,如图 3所示.一共462个台站绝大多数装配了宽频带地震仪(见表 1),多数仪器的速度平坦响应为120 s.
自2014年9月22日至2015年12月22日,自动系统共测定全球MW_GCMT≥5.5(据全球矩心矩张量(GCMT)的MW)地震254个,如图 3所示,这些地震主要分布在环太平洋地震带上,最大地震为2015年9月17日发生的智利近海MW 8.2地震.对于254个地震,共计算了6402个MWP,平均每个地震事件有25.2个台站参加MWP计算,每个台站的MWP与MW_GCMT比较见图 4,MWP-MW_GCMT的平均值为-0.08,标准差为0.39.总体而言,在越靠近MW_GCMT的地方,MWP越集中,表明上述自动测定MWP的方法可用.另外,由图 4也可见单台确定地震事件震级的准度.
由于地震产生的地震波辐射具有方向性,处于不同方位、震中距的地震台站测得的震级也会有较大差别,因此与其他震级测定方法一样,我们将多台MWP的平均值作为地震事件的震级.为此,我们计算全部台站MWP的算术平均值和标准差,误差超过标准差的值将被丢弃,重新计算剩余MWP的算术平均值MWP_AVG,将其作为事件最终MWP.MWP_AVG与MW_GCMT的对比见图 5a和图 5b,自动测定地震参数的用时(完成参数测定的时刻与发震时刻的差)见图 5c.
由图 5a和图 5b可见,MWP_AVG与MW_GCMT的一致性较好,无明显系统偏差,对于中深源地震,亦无明显系统偏差,对于MW≥7.5的地震,MWP_AVG偏小的居多.与MW_GCMT相比,MWP_AVG的平均偏差为0.16,标准差为0.20,绝大多数地震的MWP_AVG偏差位于±0.3的偏差范围内,只有个别地震的MWP_AVG偏差绝对值大于0.3,偏差绝对值小于等于0.3的地震有219个,占86.2%,在±0.3的震级偏差范围外,偏差小于-0.3的地震个数略多于偏差大于0.3的地震个数,这可能是由于对于部分地震,在计算MWP时,部分台站实时波形数据还没有满足所需的P波时间窗,致使这些台站的MWP偏小,从而导致MWP_AVG偏小.自动测定地震参数的平均用时为531 s,标准差为165 s (见图 5c中的两条虚线),用时小于600 s的地震有187个,占73.6%.
由以上分析可知,如此计算MWP可能会带来震级偏差,但与人工震后30~60分钟测定地震参数相比,时效性优势显著.将MWP震级测定方法应用于自动测定全球中强地震震级,所得结果作为准实时发布的初步震级估计,其偏差是可以接受的.
3 讨论与结论 3.1本文使用MWP震级测定方法,对地震参数自动测定系统进行了改进,实现了对全球中强地震更快更精确的自动测定.对2014年9月22日至2015年12月22日自动系统测定的254个全球MW≥6.0地震的MWP的分析结果表明,MWP与MW_GCMT的一致性较好,无明显系统偏差,对于中深源地震,亦无明显系统偏差,86.2%的MWP与MW_GCMT的偏差绝对值小于等于0.3,自动测定地震参数的平均用时为531 s,显示MWP震级测定方法在自动测定全球中强地震震级中的稳定性和可用性.
3.2准实时自动测定与发布地震参数是当今全球各个地震台网的主流,快速稳定的震级测定方法是必然的选择.相比传统的ML、MS等震级测定方法,MWP震级测定方法有2个显著优点:①由于只使用P波列,因此该方法可使用S波限幅记录,可用波形资料较多;②选取的时间窗位于S波之前,测定速度更快.由于我们强调较高的时效,MWP是自动测定的,定位的偏差、较少的台站和实时波形未满足所取的P波时间窗都可能导致计算的MWP有偏差,这也是自动测定地震参数中“快”和“准”不可兼得的普遍规律.
3.3对于MW_GCMT≥7.5得到强震,MWP偏小的可能性大,这是因为强震破裂过程复杂,可能包含多个子事件,MWP只是其中一个子事件的震级,故造成震级低估.另外,仪器频带影响不能忽略,本文中我们使用的台站大多数装配的是速度平坦响应达到120 s的仪器(见表 1),而非Tsuboi (1995, 1999)研究所使用的速度平坦响应达到360 s的仪器.即使使用速度平坦响应达到360 s的宽频带仪器,记录的P波的频带有限性也可能造成MWP饱和.虽然我们使用了MWP_corrected=(MWP_initial-1.03) / 0.843(Whitmore et al,2002)对震级进行校正,一定程度上克服了使用MWP震级测定方法测定强震时的震级偏小问题,但本文所用台站装配的仪器频带和震级范围等与Whitmore (2002)所用资料有所不同,MWP_corrected=(MWP_initial-1.03) / 0.843的适用问题有待进一步研究.由于我们在计算MWP时对速度平坦型记录进行了两次积分,从而得到地震矩轨迹,对时间序列进行积分等同于低通滤波,低频噪声不易去除,会造成震级高估,因此,无论是在同一区域还是不同区域发生的较大前震,其持续时间很长的长周期地震波会导致测定的MWP失准,使MWP偏大,影响的大小视相继发生的两个地震的时空距离以及相对震级大小.对于这些可能导致MWP失准的情况,我们需要针对性的研究作进一步改进.
3.4MWP震级测定方法是将宽频带P波位移近似为远场震源时间函数,从而确定标量地震矩和矩震级,其本身就是一个近似测定矩震级的方法.利用计算机软件自动测定MWP,只使用了P波到S波之间数十秒的波形资料,选取的时间窗口会因定位误差和理论走时误差而有偏差,另外,P波传播路径上介质的平均密度ρ和P波平均速度α也并非一个常数或一个线性函数,这些因素也会造成计算的MWP存在偏差,因此MWP与MW_GCMT存在少许偏差是正常的.而我们要精确地测定矩震级需要更确定的方法,选取更多更长、信噪比较高的波形资料,使用更接近于真实情况的速度模型,而这需要花费较长的时间,适合于人工分析与复核.
3.5地震灾害是人类共同面对的自然灾害,重大自然灾害事件应急与救援的国际化已成为一个国际合作的趋势.中国作为一个饱受地震灾害影响的国家,对全球中强地震参数的自动测定和准实时地震信息服务,有助于提高中国地震台网面对全球地震的应急处置水平,提升中国地震台网在全球的影响力,推动防震减灾科技的国际合作.
致谢 本文作图使用了GMT (Wessel et al,2013)和matplotlib网站(http://matplotlib.org/)提供的软件包,作者在此一并表示感谢.[] | Allen R M, Kanamori H .2003. The potential for earthquake early warning in southern California[J]. Science, 300 (5620) : 786–789. DOI:10.1126/science.1080912 |
[] | Allen R M, Ziv A .2011. Application of real-time GPS to earthquake early warning[J]. Geophys Res Lett, 38 : L16310. |
[] | Delouis B, Charléty J, Vallée M .2009. A method for rapid determination of moment magnitude MW for moderate to large earthquakes from the near-field spectra of strong-motion records (MWSYNTH)[J]. Bull Seismol Soc Am, 99 (3) : 1827–1840. DOI:10.1785/0120080234 |
[] | Earthworm. 2013. Earthworm[EB/OL].[2013-05-15]. http://www.earthwormcentral.org/documentation2/ovr/localmag_ovr.html. |
[] | Hara T .2007. Magnitude determination using duration of high frequency energy radiation and displacement amplitude:Application to tsunami earthquakes[J]. Earth Planets Space, 59 (6) : 561–565. DOI:10.1186/BF03352718 |
[] | Kanamori H .1977. The energy release in great earthquakes[J]. J Geophys Res, 82 (20) : 2981–2987. DOI:10.1029/JB082i020p02981 |
[] | Kanjo K, Furudate T, Tsuboi S .2006. Application of MWP to the great December 26, 2004 Sumatra earthquake[J]. Earth Planets Space, 58 (2) : 121–126. DOI:10.1186/BF03353368 |
[] | Katsumata A, Ueno H, Aoki S, et al .2013. Rapid magnitude determination from peak amplitudes at local stations[J]. Earth Planets Space, 65 (8) : 843–853. DOI:10.5047/eps.2013.03.006 |
[] | Lomax A, Michelini A, Piatanesi A .2007. An energy-duration procedure for rapid determination of earthquake magnitude and tsunamigenic potential[J]. Geophys J Int, 170 (3) : 1195–1209. DOI:10.1111/gji.2007.170.issue-3 |
[] | Lomax A, Michelini A .2009. MWPd:A duration-amplitude procedure for rapid determination of earthquake magnitude and tsunamigenic potential from P waveforms[J]. Geophys J Int, 176 (1) : 200–214. DOI:10.1111/gji.2008.176.issue-1 |
[] | Saul J, Bormann P, . 2007. Rapid estimation of earthquake size using the broadband P-wave magnitude mB[EB/OL].[2013-05-15]. ftp://ftp.gfz-potsdam.de/pub/home/st/saul/SaulBormannPosterAGU2007.pdf. |
[] | SeisComp3. 2013. SeisComp3[EB/OL].[2013-05-15]. http://www.seiscomp3.org/wiki/doc/applications/scmag. |
[] | Tezel T, Yanik K .2013. Improvement in MWP magnitude determinations and applications to earthquakes in Turkey[J]. Seismol Res Lett, 84 (6) : 991–996. DOI:10.1785/0220120160 |
[] | Tsuboi S, Abe K, Takano K, et al .1995. Rapid determination of MW from broadband P waveforms[J]. Bull Seismol Soc Am, 85 (2) : 606–613. |
[] | Tsuboi S, Whitmore P M, Sokolowski T J .1999. Application of MWP to deep and teleseismic earthquakes[J]. Bull Seismol Soc Am, 89 (5) : 1345–1351. |
[] | Vallée M, Charléty J, Ferreira A M G, et al .2011. SCARDEC:A new technique for the rapid determination of seismic moment magnitude, focal mechanism and source time functions for large earthquakes using body-wave deconvolution[J]. Geophys J Int, 184 (1) : 338–358. DOI:10.1111/gji.2010.184.issue-1 |
[] | Wessel P, Smith W H F, Scharroo R, et al .2013. Generic mapping tools:Improved version released[J]. EOS Trans AGU, 94 (45) : 409–410. |
[] | Whitmore P M, Tsuboi S, Hirshorn B, et al .2002. Magnitude dependent correction for MWP[J]. Sci Tsunami Hazards, 20 (4) : 187–192. |
[] | Wu Y M, Kanamori H, Allen R M, et al .2007. Determination of earthquake early warning parameters, τc and Pd from southern California[J]. Geophys J Int, 170 (2) : 711–717. DOI:10.1111/gji.2007.170.issue-2 |
[] | Wu Y M, Teng T L .2002. A virtual subnetwork approach to earthquake early warning[J]. Bull Seismol Soc Am, 92 (5) : 2008–2018. DOI:10.1785/0120010217 |
[] | Wu Y M, Yen H Y, Zhao L, et al .2006. Magnitude determination using initial P waves:A single-station approach[J]. Geophys Res Lett, 33 (5) : L05306. |
[] | Yamada M, Mori J .2009. Using τc to estimate magnitude for earthquake early warning and effects of near-field terms[J]. J Geophys Res, 114 (B5) : B05301. |
[] | Yamada T, Ide S .2008. Limitation of the predominant-period estimator for earthquake early warning and the initial rupture of earthquakes[J]. Bull Seismol Soc Am, 98 (6) : 2739–2745. DOI:10.1785/0120080144 |
[] | Zollo A, Lancieri M, Nielsen S .2006. Earthquake magnitude estimation from peak amplitudes of very early seismic signals on strong motion records[J]. Geophys Res Lett, 33 (23) : L23312. DOI:10.1029/2006GL027795 |
[] | 杜海林, 许力生, 陈运泰.2009. 利用阿拉斯加台阵资料分析2008年汶川大地震的破裂过程[J]. 地球物理学报, 52 (02) : 372–378. |
[] | 胡幸平, 俞春泉, 陶开, 等.2008. 利用P波初动资料求解汶川地震及其强余震震源机制解[J]. 地球物理学报, 51 (06) : 1711–1718. |
[] | 梁建宏, 孙丽, 徐志国, 等.2015. 2013年四川芦山MS 7.0地震自动速报震级偏差分析及方法改进[J]. 地震学报, 37 (6) : 983–996. |
[] | 刘成利, 郑勇, 熊熊, 等.2014. 利用区域宽频带数据反演鲁甸MS 6.5级地震震源破裂过程[J]. 地球物理学报, 2014 (09) : 3028–3037. DOI:10.6038/cjg20140927 |
[] | 刘志鹏, 盖增喜.2015. 利用多台阵压缩传感方法反演尼泊尔MW 7.9地震破裂过程[J]. 地球物理学报, 2015 (06) : 1891–1899. DOI:10.6038/cjg20150605 |
[] | 吕坚, 郑勇, 马玉虎, 等.2011. 2010年4月14日青海玉树MS 4.7、MS 7.1、MS 6.3地震震源机制解与发震构造研究[J]. 地球物理学进展, 26 (05) : 1600–1606. DOI:10.3969/j.issn.1004-2903.2011.05.012 |
[] | 王为民, 郝金来, 姚振兴.2013. 2013年4月20日四川芦山地震震源破裂过程反演初步结果[J]. 地球物理学报, 56 (04) : 1412–1417. DOI:10.6038/cjg20130436 |
[] | 张广伟, 雷建设, 梁姗姗, 等.2014. 2014年8月3日云南鲁甸MS 6.5级地震序列重定位与震源机制研究[J]. 地球物理学报, 2014 (09) : 3018–3027. DOI:10.6038/cjg20140926 |
[] | 张勇, 许力生, 陈运泰.2009. 2008年汶川大地震震源机制的时空变化[J]. 地球物理学报, 52 (02) : 379–389. |
[] | 张勇, 许力生, 陈运泰.2013. 芦山4.20地震破裂过程及其致灾特征初步分析[J]. 地球物理学报, 56 (04) : 1408–1411. DOI:10.6038/cjg20130435 |