Energy magnitude determination of moderate-strong earthquakes
-
摘要: 本文根据地震波衰减特性，开展了利用宽频带地震波形数据测定地震波能量ES和能量震级Me的方法研究。利用震中距处于20°—98°范围内的宽频带远震P波波形数据，测定了4次国外和4次国内中强震的能量震级Me，并对其面波震级MS、矩震级MW及能量震级Me进行了分析对比。结果表明：面波震级MS表示的是地震在某一固定频率所辐射的地震波能量大小；矩震级MW与地震所产生的断层长度、断层宽度、震源破裂的平均位错量等静态构造效应密切相关；而能量震级Me反映的是震源动态特征，与地震震源的动力学特性密切相关。由于地震是以地震波形式辐射，能量主要集中在震源谱的拐角频率附近，因此能量震级Me更适合描述地震的破坏性。由此可见，联合测定面波震级MS，矩震级MW和能量震级Me对于地震定量研究以及地震灾害与风险评估具有重要作用。
-
关键词:
- 地震波能量 /
- 能量震级 /
- 宽频带远震P波波形数据 /
- 震源动态特征 /
- 震源谱
Abstract: This paper, based on the seismic wave attenuation characteristics, explores a method for utilizing broadband seismic waveform data to measure seismic energy ES and energy magnitude Me. The energy magnitudes Me of four overseas and four domestic moderate-strong earthquakes are measured by using broadband teleseismic P waveform data with epicentral distance in the range of 20°−98°. The comparison of surface wave magnitude MS, moment magnitude MW and energy magnitude Me of several earthquakes shows that: surface wave magnitude MS represents the magnitude of the energy of seismic wave radiating at a fixed frequency; the moment magnitude MW is closely related to the static structural effects such as the length of fault, the width of fault and the average dislocation amount of the source rupture produced by the earthquake; the energy magnitude Me reflects source dynamic characteristics and is closely related to the dynamic characteristics of the seismic source. The energy radiated by the earthquake in the form of seismic waves is mainly concentrated near the corner frequency of the source spectrum, so the energy magnitude Me is more suitable to describe the destructiveness of the earthquake. Therefore, joint determination of surface wave magnitude MS, moment magnitude MW and energy magnitude Me are of great significance for quantifying and evaluating the disasters caused by earthquakes or tsunamis. -
引言
为了尽可能地减少地震灾害造成的损失，多年来地震学家一直专注于地震前兆信息提取和地震预测研究（屈春燕等，2006；马晓静等，2009；张建国等， 2012， 2013；Zhang et al，2017）。Gorny等（1988）分析得出，1984年加兹利MS7.1地震震前曾出现显著的卫星热红外亮温异常，此后，与地震相关的卫星热红外亮温异常研究发展迅速，逐步产生了多种卫星热红外数据分析方法。例如卫星技术稳健估计方法（robust satellite technique，简写为RST）（Tronin et al，2002；Filizzola et al，2004；Ouzounov，Freund，2004；邓志辉等， 2005， 2014；单新建等，2005）、内外温差法（陈顺云等， 2006， 2013；Pulinets et al，2006；陈梅花等，2007；Saraf et al，2008）和相对功率谱方法（张元生等，2010，2011；Blackett et al，2011；Xie et al，2013；解滔等，2015a，b；张治广等，2017a）等。Filizzola等（2004）认为在地震孕育过程中，震中的地表温度高于其它区域，因此采用一种基于统计的RST方法，该方法将研究区的长期背景温度值从研究时段内研究区的亮温数据中减去，从而得到了研究时段的显著温度异常变化，并于1999年雅典MS5.9地震前提取到了显著的异常。为了进一步研究地震孕育过程中断裂带内外温度差异，陈梅花等（2007）以断裂带为界，分析了断裂带内外的温度差异及其与地震的关系，并分析得出2000年姚安MS6.5地震前，红河断裂内外温差呈稳定上升趋势，震后迅速回落至正常水平；张元生等（2010，2011）利用相对功率谱方法提取到了2011年日本MS9.0和2008年汶川MS8.0地震前的热红外亮温异常；Xie等（2013）沿用该相对功率谱方法，采用具有可变时频窗的连续小波变换计算了亮温数据的能谱，不同频段热红外相对小波能谱时空演化特征表明，该方法能够同时兼顾能谱分析的频率和时间分辨率，得到的相对小波能谱变化时频域更为平滑，并能够很好的抑制边界效应。
据中国地震台网中心测定，2018年5月28日1时50分，吉林省松原市宁江区发生MS5.7地震，震中位置为（124.71°E，45.27°N），震源深度为13 km。震中位于松原油气开发区，且处于第二松花江断裂前郭部分与扶余—肇东断裂的交会区（图 1）。地震孕育过程中，断层活动必然会导致CO2，CH4等“温室气体”释放增强，而此类气体极易吸收太阳的直接辐射和地表再发射的红外辐射，故使得地表温度逐渐升高，因此若能够检测到红外异常，一方面可以支持前人关于震前热红外亮温高值异常的研究结果（张元生等，2010；解滔等，2015a，b；张治广等，2017a），另一方面可以更有力地说明红外异常的温室效应机制。此外，由于此次地震震中位于两条断裂带的交会区，不同断裂对本次地震的贡献可通过相应地区热红外演化特征来进行揭示。
![]()
图 1 2018年5月28日松原MS5.7地震震中位置及周边的断层分布Figure 1. Location of the epicenter of Songyuan MS5.7 earthquake on May 28，2018 and faults distribution in the vicinity of the source region基于此，本文拟利用相对小波能谱分析方法，对松原MS5.7地震前（119°E—134°E，40°N—50°N）范围内的静止气象卫星FY-2G，FY-2E在2016年5月29日至2018年5月28日的连续观测亮温数据予以分析，以获取与此次地震相关的地震热红外异常演化特征，以期为进一步探索卫星热红外亮温异常与地震的关系积累研究震例。此外，由于本次地震的余震空间分布并未呈现出与扶余—肇东断裂与第二松花江断裂前郭部分相近的优势分布，因此此次地震的发震构造一直无法确定。本文研究分析得出，热红外相对小波能谱沿着两条断裂衰减，其中沿扶余—肇东断裂衰减较慢，因此可初步认为，扶余—肇东断裂对本次地震的贡献较为显著，可能为进一步揭示发震构造提供佐证。
1. 数据选取和处理方法
1.1 数据选取
本文选取中国静止气象卫星FY-2E和FY-2G （FY-2G为FY-2E的替代卫星）观测的热红外遥感亮温产品数据进行分析（中国气象局国家卫星气象中心，2008）。FY-2E和FY-2G两颗卫星分别于2008年6月15日和2014年12月31日发射，定点于105°E赤道上空，每1小时或30分钟对全球约1/3面积观测1次，红外波段分别为10.3—11.3 μm和11.5—12.5 μm （张璇等，2016；张治广等，2017a，b ）。为了避免白天太阳照射导致的地表大幅升温造成的干扰，选取地方时夜间即北京时间00:00—04:00 （GMT 16:00—20:00）时段的5次亮温观测资料进行分析。
1.2 数据处理
小波变换作为一种线性时频分析方法，在频率域和时间域均具有很好的局部化性质，能够较好地解决频率分辨率与时间分辨率的矛盾，因此被广泛应用于地球物理勘探、地震信号处理及地震前兆分析等各项研究中（解滔等，2015a，b；王黎宏，2016）。对于一组有限的时间序列f （t），定义其小波变换为
$ {W_\varPsi }f(a{\text{，}}b) {\text{＝}} \int\nolimits_{ {\text{－}} \infty }^{ {\text{＋}} \infty } {f(t)\varPsi ^*_{a{\text{，}}b}} (t)\;{\rm{d}}t{\text{，}} $
(1) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
$ \varPsi {\rm{(}}\omega {\rm{)}} {\text{＝}}\frac {1} {\sqrt[{4}]{{\text{π}}}}{\rm{exp}}{\rm{[}}{ {\text{－}} {\frac {{{(\omega - {\omega _0})}^2}} {2}}}{\rm{]}} {\text{，}} $
(2) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
$ {R_W}(a{\text{，}}b) {\text{＝}}{\frac {|{W_\varPsi }f(a{\text{，}}b){|^2}} {{\overline W ^2}(a{\text{，}}b)}}{\text{，}} $
(3) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
${\overline W ^2}(a{\text{，}}b) {\text{＝}} \displaystyle\frac{1}{N}\sum\limits_{l {\text{＝}} 0}^{N {\text{－}} 1} {|{W_l}f(a{\text{，}}b){|^2}} {\text{，}} $
(4) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
Torrence和Compo （1998）对连续小波变换的边界效应的研究结果表明：小波能谱在接近信号边界处呈明显衰减状态，且衰减时间随着频率降低而变长；而在边界处的实际能谱值相较于利用小波变换获取的相对能谱值高很多，即如果信号在边界表现出能谱高值异常，实际的能谱异常幅值将大于计算结果，因此利用连续小波变换进行数据处理时并不会由于边界效应导致虚假异常的产生（解滔等，2015a，b）。连续小波变换为非正交变换，不仅在时域和频域有混叠，频域内相邻频段信息中也含有重叠的成分，因此本文取频率间隔对尺度因子a＝0.5进行离散化处理（Torrence，Compo，1998；解滔等，2015b）；然后对所选观测资料的每个像元进行计算，数据时间长度取N＝730 （闰年时N＝731）；之后应用Xie等（2013）改进的连续小波变换对FY-2E和FY-2G卫星亮温观测数据进行处理，计算得到信号相对小波能谱时频信息后，对研究区域（119°E—134°E，40°N—50°N）进行时空扫描，分析异常区域的时空演化特征及其与本次松原MS5.7地震的关系。
当研究区域的像元上空有云层覆盖时，静止卫星FY-2E和FY-2G观测到的数据为云顶的温度，而云层下地表温度要远高于云顶的温度，为了去除云层影响，采用补窗法进行处理（张元生等，2011；解滔等，2015b），即以补窗后整点值的平均值作为像元的日均值。对所有像元进行相同处理，计算得到像元的日均值序列，采用低通滤波器提取温度日均值数据中周期T≥1年的趋势成分，计算二者的均方差，并以1.5倍均方差作为数据下限阈值，去除云层的影响（Xie et al，2013；解滔等，2015b）。本文按照上述数据处理方法，计算了松原MS 5.7地震前两年间所有像元的相对小波能谱值，并对计算结果进行了时频扫描，得到地震前后热红外亮温小波相对能谱的时空演化结果。
2. 结果分析
2.1 2018年5月28日松原MS5.7地震热红外异常分析
为了尽可能地排除离散噪声异常像元的影响，分析可能与地震等构造活动有关的热红外亮温高值异常现象，本文按照以下3条准则进行异常信息判别（解滔等，2015a，b）：① 异常像元须聚集成一定范围的异常区域，而非零星离散的分布；② 显著异常持续时间须在15天以上；③ 异常区域须沿构造断裂带分布，特别是活动断裂带。
对所选数据的小波相对能谱的时空变化特征（图 2）进行分析可以看出：2018年5月28日松原MS5.7地震前，震中附近的断裂带出现了明显的热红外亮温异常现象。从空间上看，2018年2月初在震中附近NW向的第二松花江断裂前郭部分出现异常；2018年3月，异常区开始沿NE向的扶余—肇东断裂展布，3月中旬异常空间分布范围达到峰值，随后异常区域逐渐缩小；至2018年4月仅NE向的扶余—肇东断裂区域存在条带状分布的异常区域，直到2018年4月中旬异常消失。异常空间分布呈现出NW向与NE向的条带交会状，即“V”型分布特征，未来震中位于两条带状热红外异常区域在迁移过程中形成的交会带附近。从相对小波能谱异常幅度上看，2018年2月初在震中附近NW向的第二松花江断裂前郭部分出现异常，随后异常幅度逐渐增大，3月中旬异常幅度达到最大，能谱相对变化达到均值的7倍以上，之后异常幅值逐渐衰减，直至2018年4月中旬异常消失，异常持续时间达3个月。本次地震热异常呈现出与以往震例类似的演化过程，即由初始小范围异常发展至大区域高能量异常，随后开始急速减弱的过程（张璇等，2016）。因此，2018年2月至4月于松原地区出现的热红外亮温异常可能与此次松原MS5.7地震有关。
![]()
图 2 2018年5月28日松原MS5.7地震前亮温小波相对能谱的异常时空演化Figure 2. Spatio-temporal evolution of relative energy spectrum of bright temperature for Songyuan MS5.7 earthquake on May 28，20182.2 研究区历史震例分析
尽管地震前震中及附近区域表现出的地表增温异常现象已被大量的观测结果所证实，但由于地震孕育机理的复杂性和地震热红外异常的多样性，地震热红外成因机理的研究仍然是一个难点（解滔等，2015a）。郭增建和秦保燕（1979）认为，大震震源区锁住的断层会发生预滑移，沿断层带的温度会升高，因而产生震前热红外亮温异常现象。Cicerone等（2009）讨论了地震前卫星热红外亮温异常现象，结果显示，在富含地热和油气的区域地震前出现热红外亮温异常现象的可能性更大，因此解滔等（2015b）分析得出，地震孕育过程中应力逐渐增大，当应力接近岩石临界破裂状态时，会出现大量贯通地表的裂隙，地下逸出气体也随之增多，导致地下岩层与地表的热对流加强。大量裂隙的出现及热对流的增强导致了地震前热红外异常的产生，且异常区沿震中附近的断裂带集中分布（Tronin，1996，2000；马瑾等，2006；强祖基等，2009，2010；解滔等，2015b）。
需要特别说明，热红外异常只是表明在异常时段该区域地表温度较正常背景出现了异常变化，而温度变化并非仅由地震所致，因此，地震热异常分析过程中不可避免地会出现虚假地震异常现象。本文将2013年以来研究区域内热红外亮温异常与地震对应关系列于表 1。在2013年1月至2018年5月期间，依据本文的异常判别准则，该研究区内共出现了7组热红外亮温相对能谱异常，其中3组异常出现后无地震对应，另外4组分别对应于2013年4月22日内蒙古与辽宁交界处的MS5.3地震、2013年6月20日内蒙古与黑龙江交界处MS5.0地震、2013年10月31至11月23日吉林松原市前郭MS5.8震群和本次2018年5月28日松原MS5.7地震，这4组地震前均在震中附近出现了亮温热红外相对能谱高值异常，且地震均发生在异常区边缘。在2013年1月至2018年5月期间，研究区域内共发生MS≥5.0地震5组，除上述4组地震外，还有2016年1月2日发生在黑龙江省牡丹江市林口县的MS6.4地震，该次地震的震源深度为580 km，为深源地震，但风云卫星没有记录到该地震产生的亮温异常。为便于对比分析，选取异常幅值最大的日期作为异常极值日，有震异常和无震异常的极值空间分布如图 3所示。可以看出，有震异常空间分布较为聚集，异常幅值较大，而无震异常空间分布较有震异常更零散，异常幅值相对较小。
表 1 2013—2018年研究区热红外亮温异常与地震对应情况Table 1. The thermal infrared anomaly and corresponding earthquakes in the studied area异常出现时间 异常空间范围 对应地震 2013年1—3月 121°E—123°E，42°N—44°N 2013年4月22日内蒙古与辽宁交界MS5.3地震 2013年2—3月 124°E—126°E，46°N—50°N 2013年6月20日内蒙古与黑龙江交界MS5.0地震 2013年10—12月 121°E—124°E，42°N—46°N 2013年10月31日—11月23日吉林省松原市前郭MS5.8震群 2014年4—7月 129°E—131°E，44°N—46°N 无 2015年3—5月 121°E—122°E，47°N—49°N 无 2016年6—7月 121°E—122°E，40°N—41°N 无 2018年2—4月 125°E—127°E，44°N—46°N 2018年5月28日吉林松原MS5.7地震 ![]()
图 3 历史地震及热红外相对能谱极值的空间分布（a） 有震异常相对能谱极值空间分布；（b） 无震异常相对能谱极值空间分布Figure 3. Spatial distribution of thermal infrared extreme value anomaly and history earthquakes（a） Spatial distribution of thermal infrared extreme value anomaly and corresponding earthquakes；（b） Spatial distribution of thermal infrared extreme value anomaly without corresponding earthquakes距离松原MS5.7地震最近的为2013年10月31日至11月23日松原MS5.8震群，两次地震震中相距仅为88 km。本文以该震群为例，对比分析了两组地震前后热红外亮温异常特征及其差异，如图 4所示。从异常及发震时间上看，松原MS5.7地震发生于热红外亮温异常消失后约1个月，而2013年10月31日至11月23日松原MS5.8震群发生在异常出现至异常幅值逐渐增大的过程中，地震后异常幅值逐渐减小。从空间上看，两组地震均发生在异常区边缘，且邻近异常衰减最慢区域。不同的是，2018年5月28日松原MS5.7地震前出现的热红外亮温异常呈现出显著的沿发震构造分布的特点，即异常主要沿第二松花江断裂前郭地区和扶余—肇东断裂展布，而2013年松原震群前异常首先出现在震中距较远的嫩江断裂，其后才出现在北东向的扶余—肇东断裂处。这种大范围远距离的热红外亮温异常可能与地下逸出气体增多以及地下岩层与地表的热对流增强有关（徐秀登等，1991；强祖基等，1992；Tronin，1996，2000）。
![]()
图 4 2013年松原前郭MS5.8震群前后亮温小波相对能谱的异常时空演化Figure 4. Temporal-spatial evolution of relative energy spectrum of brightness temperature for Songyuan Qianguo MS5.8 earthquake in 20133. 讨论与结论
本文针对2018年5月28日吉林松原MS5.7地震，利用连续小波相对能谱方法对发震前两年期间中国静止气象卫星FY-2E （FY-2G）连续观测的亮温资料进行分析，分析结果显示，在震中附近的第二松花江断裂前郭段及扶余—肇东断裂区域均在震前出现显著的相对小波能谱高值异常。震前两个月左右异常幅值和规模开始衰减并于4月中旬消失。地震震中位于异常区边缘的活动断层上，发震时间为异常峰值后的两个月内，与前人的研究结果基本一致。因而认为本文讨论的2018年2月至4月松原地区热红外异常可能与松原MS5.7地震有关。通过对2013年以来研究区历史地震及其与热红外亮温相对小波能谱高值异常的统计可初步发现，有震异常的相对小波能谱高值异常在空间上通常较为集中，呈聚集分布，且异常幅值较高；而无震异常的相对小波能谱高值异常空间上分布较为零星，异常幅值也较低。但由于统计样本个数的限制，目前还无法定量给出两者之间的差别，更多异常特征和异同需要基于更大样本的统计。因此，下一步拟基于中国大陆地区历史地震开展有震异常和无震异常的定量分析。
-
![]()
图 1 计算九寨沟地震能量震级所用台站分布图
Figure 1. Distribution of the stations used for calculation of energy magnitude of Jiuzhaigou earthquake
![]()
图 2 九寨沟地震能量震级Me残差分布
Figure 2. Distribution of the residuals for energy magnitude Me of the Jiuzhaigou earthquake
![]()
图 3 计算芦山地震能量震级所用台站分布图
Figure 3. Distribution of the stations used for calculation of energy magnitude of Lushan earthquake
![]()
图 4 芦山地震能量震级Me残差分布
Figure 4. Distribution of the residuals for energy magnitude Me of the Lushan earthquake
![]()
图 5 计算哥斯达黎加地震（a）和尼加拉瓜地震（b）能量震级所用台站分布图
Figure 5. Distribution of the stations used for calculation of energy magnitude of Costa Rica earthquake （a） and Nicaragua earthquake （b）
![]()
图 6 牡丹江地震台记录到的哥斯达黎加地震（a）和尼加拉瓜地震（b）的波形记录图
Figure 6. Comparison of P-waveforms from Costa Rica earthquake （a） and Nicaragua earthquake （b） recorded by the station MDJ
表 1 计算九寨沟地震能量震级Me的单台数据
Table 1 Single station data for calculating energy magnitude Me of Jiuzhaigou earthquake
台站名 震中距/° Me 台站名 震中距/° Me 台站名 震中距/° Me IU.MAKZ 21.46 7.02 IU.MBWA 56.42 6.22 IU.PAB 81.55 6.18 BK.CMB 97.75 6.66 IU.WAKE 57.19 6.47 IU.FUNA 82.63 6.21 II.AAK 24.86 6.70 IU.KBS 57.87 6.24 IU.KIP 84.85 6.31 II.KURK 25.44 6.45 IU.ADK 58.32 6.46 II.TAU 85.86 6.34 IU.MAJO 28.23 6.36 IU.PMG 59.25 6.16 IU.LSZ 86.81 6.12 II.SIMI 28.57 6.50 II.MSEY 59.38 6.17 II.MSVF 87.12 6.66 II.BRVK 31.06 6.19 II.WRAB 60.58 6.06 IU.POHA 87.71 6.53 II.ERM 32.05 6.08 II.KWAJ 63.29 6.58 II.FFC 89.31 6.47 IU.YSS 32.37 6.12 IU.KONO 63.83 6.15 IU.COR 91.20 6.34 IU.DAV 32.98 6.76 II.ALE 64.10 6.08 IU.AFI 92.97 6.75 IU.YAK 33.20 6.30 IU.FURI 64.42 6.24 GT.LBTB 94.43 6.58 II.ARU 38.59 6.09 IU.MIDW 66.27 6.45 II.CMLA 94.85 6.48 IU.MA2 40.40 6.09 IU.CTAO 66.90 6.05 IU.XMAS 96.14 6.28 II.KAPI 41.07 6.31 IU.NWAO 67.31 6.56 IU.KOWA 96.41 6.95 IU.GUMO 42.06 6.53 IU.COLA 67.41 6.12 IU.RAO 96.80 6.22 II.UOSS 42.10 6.25 IU.HNR 68.24 6.57 IU.TSUM 97.39 7.07 IU.PET 43.32 5.89 II.KDAK 68.99 6.36 碾子山 20.30 6.63 IU.GNI 47.05 6.22 II.BFO 69.18 5.94 五常 21.50 6.55 II.KIV 47.89 6.29 IU.KMBO 71.25 5.92 讷河 21.77 6.01 IU.BILL 49.33 6.47 IU.TARA 71.83 6.85 宾县 21.97 6.35 II.DGAR 50.38 6.52 II.ESK 71.91 6.02 宝清 21.98 6.47 II.OBN 50.66 6.18 II.BORG 73.18 6.28 靖宇 21.06 6.14 II.RAYN 51.60 6.23 II.ABPO 75.12 6.18 柳河 20.95 6.21 II.LVZ 51.76 6.03 II.MBAR 76.27 6.63 丰满 21.55 6.42 IU.KEV 54.58 5.97 IU.SFJD 77.95 6.41 抚松 20.65 6.31 IU.ANTO 55.99 6.19 IU.JOHN 78.17 6.01 牡丹江 22.91 6.55 平均能量震级 6.3 注：表中10个中国台站的地震引自国家测震台网数据备份中心（2007）和郑秀芬等（2009）。 
下载: 导出CSV
表 2 不同机构测定的4次国外和4次国内地震的各类震级
Table 2 Various magnitudes measured by different institutes for four overseas and four domestic earthquakes
发震时间 发震地点 震源深度
/km震源破裂
时间/s震源机制
类型MS
（CENC）MW
（GCMT）Me
（IRIS）Me
（本文）1991-04-22 哥斯达黎加 10.0 81 T 8.1 7.6 7.4 7.4 1992-09-02 尼加拉瓜 15.0 172 T 7.7 7.6 6.7 6.6 2006-11-15 昆士兰 39.0 133 T 8.0 8.3 7.7 7.7 2007-01-13 昆士兰 10.0 101 N 7.9 8.1 8.6 8.1 2008-05-12 四川汶川 12.8 126 TS 8.0 7.9 8.1 8.0 2013-04-20 四川芦山 12.2 39 T 7.0 6.6 6.7 6.6 2014-08-03 云南鲁甸 10.0 46 S 6.5 6.2 6.4 6.4 2017-08-08 四川九寨沟 9.0 43 S 7.0 6.5 6.4 6.3 注：T为逆断层，N为正断层，S为走滑断层，TS为包含走滑成分的逆断层。
下载: 导出CSV
-
陈培善. 1990. 地震定量的国际现状[J]. 地震地磁观测与研究, 11（3）: 33-39. Chen P S. 1990. The international status of earthquake quantification[J]. Seismological and Geomagnetic Observation and Research, 11（3）: 33-39 （in Chinese）.
国家测震台网数据备份中心. 2007. 国家测震台网地震波形数据[DB/OL]. 中国地震局地球物理研究所. http://www. seisdmc.ac.cn. doi: 10.11998/SeisDmc/SN. Data Management Centre of China National Seismic Network. 2007. Waveform data of China national seismic network[DB/OL]. Institute of Geophysics, China Earthquake Administration. http://www.seisdmc.ac.cn. doi: 10.11998/SeisDmc/SN （in Chinese）.
李志强, 侯建盛, 李洋, 苏桂武, 姜立新, 孙柏涛, 袁一凡, 林均岐. 2013. 2013年4月20日四川芦山MS7.0地震灾害特点分析[J]. 地震地质, 35（2）: 398-410. doi: 10.3969/j.issn.0253-4967.2013.02.018 Li Z Q, Hou J S, Li Y, Su G W, Jiang L X, Sun B T, Yuan Y F, Lin J Q. 2013. Analysis on the Characteristics of the MS7.0 Lushan, Sichuan province, earthquake hazard on April 20, 2013[J]. Seismology and Geology, 35（2）: 398-410 （in Chinese）.
赵仲和. 2014. 2014年8月3日云南鲁甸M6.5地震的地震波能量[J]. 国际地震动态, （9）: 24-28, 34. doi: 10.3969/j.issn.0235-4975.2014.09.007 Zhao Z H. 2014. Seismic energy of Ludian M6.5 earthquake, August 3, 2014[J]. Recent Developments in World Seismology, （9）: 24-28, 34 （in Chinese）.
郑秀芬, 欧阳飚, 张东宁, 姚志祥, 梁建宏, 郑洁．2009. " 国家测震台网数据备份中心”技术系统建设及其对汶川大地震研究的数据支撑[J]．地球物理学报, 52（5）: 1412—1417． doi: 10.3969/j.issn.0001-5733.2009.05.031 Zheng X F, Ouyang B, Zhang D N, Yao Z X, Liang J H, Zheng J. 2009. Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake. Chinese Journal of Geophysics, 52（5）:1412—1417（in Chinese）.
中国地震局. 2013. 中国地震局发布四川省芦山" 4·20”7.0级强烈地震烈度图[EB/OL]. [2018−11−15]. https://www.cea.gov.cn/cea/xwzx/fzjzyw/5195229/index.html. China Earthquake Administration. 2013. China Earthquake Administration release the seismic intensity map of Lushan Sichuan M7.0 strong earthquake on April 20, 2017[EB/OL]. [2018-11-15]. https://www.cea.gov.cn/cea/xwzx/fzjzyw/5195229/index.html （in Chinese）.
中国地震局. 2017. 中国地震局发布四川九寨沟7.0级地震烈度图[EB/OL]. [2018−11−15]. https://www.cea.gov.cn/cea/xwzx/fzjzyw/5206273/index.html. China Earthquake Administration. 2017. China Earthquake Administration release the seismic intensity map of Jiuzhaigou Sichuan M7.0 earthquake[EB/OL]. [2018-11-15]. https://www.cea.gov.cn/cea/xwzx/fzjzyw/5206273/index.html （in Chinese）.
Aki K, Richards P G. 1980. Quantitative Seismology: Theory and Methods[M]. San Francisco: W. H. Freeman: 932.
Ammon C J, Kanamori H, Lay T. 2008. A great earthquake doublet and seismic stress transfer cycle in the central Kuril islands[J]. Nature, 451（7178）: 561-565. doi: 10.1038/nature06521
Boatwright J, Choy G L. 1986. Teleseismic estimates of the energy radiated by shallow earthquakes[J]. J Geophys Res, 91（B2）: 2095-2112. doi: 10.1029/JB091iB02p02095
Bormann P, Liu R F, Ren X, Gutdeutsch R, Kaiser D, Castellaro S. 2007. Chinese National Network magnitudes, their relation to NEIC Magnitudes, and recommendations for new IASPEI magnitude standards[J]. Bull Seismol Soc Am, 97（1B）: 114-127. doi: 10.1785/0120060078
Bormann P, Saul J. 2008. The new IASPEI standard broadband magnitude mB[J]. Seismol Res Lett, 79（5）: 698-705. doi: 10.1785/gssrl.79.5.698
Bormann P, Saul J. 2009a. Earthquake magnitude[G]//Encyclopedia of Complexity and Systems Science. Heidelberg: Springer: 2473−2496.
Bormann P, Saul J. 2009b. A fast, non-saturating magnitude estimator for great earthquakes[J]. Seism Res Lett, 80（5）: 808-816. doi: 10.1785/gssrl.80.5.808
Choy G L, Boatwright J L. 1995. Global patterns of radiated seismic energy and apparent stress[J]. J Geophys Res, 100（B9）: 18205-18228. doi: 10.1029/95JB01969
Choy G L, Boatwright J. 2007. The energy radiated by the 26 December 2004 Sumatra-Andaman earthquake estimated from 10-minute P-wave windows[J]. Bull Seismol Soc Am, 97（1A）: S18-S24. doi: 10.1785/0120050623
Convers J A, Newman A V. 2011. Global evaluation of large earthquake energy from 1997 through mid-2010[J]. J Geophys Res, 116（B8）: B08304.
Di Bona M, Rovelli A. 1988. Effects of bandwidth limitation on stress drop estimated from integrals of the ground motion[J]. Bull Seismol Soc Am, 78: 1818-1825.
Di Giacomo D, Grosser H, Parolai S, Bormann P, Wang R J. 2008. Rapid determination of Me for strong to great shallow earthquakes[J]. Geophys Res Lett, 35（10）: L10308.
Di Giacomo D. 2010. Determination of the Energy Magnitude ME: Application to Rapid Response Purposes and Insights to Regional/Local Variabilities[D]. Potsdam: Universität Potsdam: 49−57.
Duda S J, Yanovskaya T B. 1993. Spectral amplitude-distance curves for P-waves: effects of velocity and Q-distribution[J]. Tectonophysics, 217（3/4）: 255-265.
Goes S D B, Velasco A A, Schwartz S Y, Lay T. 1993. The April 22, 1991, Valle de la Estrella, Costa Rica （MW=7.7） earthquake and its tectonic implications: A broadband seismic study[J]. J Geophys Res, 98（B5）: 8127-8142. doi: 10.1029/93JB00019
Gutenberg B, Richter C F. 1954. Seismicity of the Earth and Associated Phenomena[M]. 2nd ed. Princeton: Princeton University Press: 1−15.
Haskell N A. 1964. Total energy spectra density of elastic waves for propagating faults[J]. Bull Seismol Soc Am, 54（6）: 1811-1841.
Ide S, Beroza G C. 2001. Does apparent stress vary with earthquake size?[J]. Geophys Res Lett, 28（17）: 3349-3352. doi: 10.1029/2001GL013106
Kanamori H, Mori J, Hauksson E, Heaton T H, Hutton L K, Jones L M. 1993. Determination of earthquake energy-release and ML using TERRAscope[J]. Bull Seismol Soc Am, 83（2）: 330-346.
Kanamori H, Kikuchi M. 1993. The 1992 Nicaragua earthquake: a slow tsunami earthquake associated with subducted sediments[J]. Nature, 361（6414）: 714-716. doi: 10.1038/361714a0
Kennett B L N, Engdahl E R, Buland R. 1995. Constraints on seismic velocities in the Earth from traveltimes[J]. Geophys J Int, 122（1）: 108-124. doi: 10.1111/gji.1995.122.issue-1
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
Montagner J P, Kennett B L N. 1996. How to reconcile body-wave and normal-mode reference Earth models[J]. Geophys J Int, 125（1）: 229-248. doi: 10.1111/gji.1996.125.issue-1
NEIC. 1991. Today in earthquake history[EB/OL]. [2018−11−15]. https://earthquake.usgs.gov/learn/today/index.php? month=4&day=22.
NEIC. 1992. M7.7−near the coast of Nicaragua[EB/OL]. [2018−11−15]. https://earthquake.usgs.gov/earthquakes/eventpage/ usp0005ddn/executive.
NEIC. 2006. M8.3−Kuril Islands[EB/OL]. [2018−11−15]. https://earthquake.usgs.gov/earthquakes/eventpage/usp000exfn/executive.
NEIC. 2007. M8.1−east of the Kuril Islands[EB/OL]. [2018−11−15]. https://earthquake.usgs.gov/earthquakes/eventpage/ usp000f2ab/executive.
Newman A V, Okal E A. 1998. Teleseismic estimates of radiated seismic energy: the E/M0 discriminant for tsunami earthquakes[J]. J Geophys Res, 103（B11）: 26885-26898. doi: 10.1029/98JB02236
Parolai S. 2009. Denoising of seismograms using the S transform[J]. Bull Seismol Soc Am, 99（1）: 226-234. doi: 10.1785/0120080001
Richter C F. 1935. An instrumental earthquake magnitude scale[J]. Bull Seismol Soc Am, 25（1）: 1-32.
Rudnicki J W, Freund L B. 1981. On energy radiation from seismic sources[J]. Bull Seismol Soc Am, 71（3）: 583-595.
Singh S K, Ordaz M. 1994. Seismic energy release in Mexican subduction zone earthquakes[J]. Bull Seismol Soc Am, 84（5）: 1533-1550.
Stockwell R G, Mansinha L, Lowe R P. 1996. Localization of the complex spectrum: the S transform[J]. IEEE Trans Signal Process, 44（4）: 998-1001. doi: 10.1109/78.492555
Venkataraman A, Kanamori H. 2004a. Effect of directivity on estimates of radiated seismic energy[J]. J Geophys Res, 109（B4）: B04301.
Wang R J. 1999. A simple orthonormalization method for stable and efficient computation of Green’s functions[J]. Bull Seismol Soc Am, 89（3）: 733-741.
-
期刊类型引用(3)
1. 李继业，胡澜缤，李营，马龙辰，王强，张思萌，李冬妮. 松原M_S5.1地震前断层土壤气H_2、Hg地球化学特征与热红外异常响应研究. 地震工程学报. 2023(04): 933-945 . 
百度学术
2. 李利波，张翔，李智蓉. 2021年云南省3次中强地震前后热红外异常特征分析. 地震地磁观测与研究. 2022(05): 16-22 . 百度学术
3. 米日阿依·买土地，玉素甫江·如素力，喀迪阿依·阿力木，艾尔肯·图尔孙. 基于STL和GESD算法的于田两次Ms7.3地震热异常时空分布特征识别研究. 河南科学. 2021(03): 421-428 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 1700
- HTML全文浏览量: 1056
- PDF下载量: 97
- 被引次数: 5
