文章快速检索     高级检索
  大地测量与地球动力学  2020, Vol. 40 Issue (5): 534-539  DOI: 10.14075/j.jgg.2020.05.018

引用本文  

郁广程, 姚运生, 张丽芬, 等. 2013年吉林省松原地区MS≥5地震群成因研究[J]. 大地测量与地球动力学, 2020, 40(5): 534-539.
YU Guangcheng, YAO Yunsheng, ZHANG Lifen, et al. Study on Genesis of MS≥5 Songyuan Earthquake Swarm in Jilin Province in 2013[J]. Journal of Geodesy and Geodynamics, 2020, 40(5): 534-539.

项目来源

国家自然科学基金(41572354,41772840)。

Foundation support

National Natural Science Foundation of China, No. 41572354, 41772840.

通讯作者

姚运生,博士,研究员,主要研究方向为防灾减灾工程及防护工程,E-mail:yunsheng@eqhb.gov.cn

Corresponding author

YAO Yunsheng,PhD,researcher,majors in disaster prevention and mitigation engineering,E-mail:yunsheng@eqhb.gov.cn.

第一作者简介

郁广程,硕士生,主要研究方向为数字地震学,E-mail:guangchengyu@163.com

About the first author

YU Guangcheng,postgraduate,majors in broadband seismology,E-mail:guangchengyu@163.com.

文章历史

收稿日期:2019-05-24
2013年吉林省松原地区MS≥5地震群成因研究
郁广程1     姚运生1,2     张丽芬1     王秋良1     吴海波1     
1. 中国地震局地震研究所地震大地测量重点实验室,武汉市洪山侧路40号,430071;
2. 防灾科技学院,河北省三河市燕郊学院街465号,065201
摘要:2013-10-31~11-23吉林省松原地区连续发生5次MS≥5.0地震,其中最大地震为2013-11-23 MS5.8地震,该震群的突发性及密集性与区域历史构造地震特征不相符。通过计算2013年松原地区5次MS≥5地震零频极限值的相关系数发现,第1次和第5次地震的震源机制相似性较大,第2次、第3次和第4次地震的震源机制相似性较大。将松原地震同典型的构造地震和塌陷地震的地震参数进行比较发现,2013年松原地震的位移震源谱较典型构造地震呈现衰减早、衰减快的特征,波形复杂度和谱比值等时域、频域特征参数也与典型构造地震相差较大,更符合塌陷型地震的特征。
关键词构造地震塌陷地震地震成因频谱分析震源谱

2013-10-31~11-23吉林省松原市前郭尔罗斯蒙古族自治县连续发生了一系列地震,其中MS≥5.0地震就有5次,最大震级为MS5.8(表 1)。根据历史地震目录记载,松原地区地壳稳定,地震发生的频率远小于其他地区,强度也不大,该地震群发生之前震区附近仅发生过两次中强地震[1],震中位置见图 1。松原地区地震的成因,引起了研究者们的关注,高金哲等[2]推测可能与油气开采有关。对于2013年吉林松原地震群的发震机制,不少学者从地震精定位、震源机制解、发震构造等多个方面进行研究,得到不同的结论[3]

表 1 2013年松原地区MS≥5.0地震基本地震参数 Tab. 1 Basic seismic parameters of Songyuan MS≥5 earthquake swarm in 2013

图 1 松原地区MS≥5.0地震震中分布(截至2013-11-23) Fig. 1 Epicenter distribution of MS≥5.0 earthquakes in Songyuan

前人利用蕴含丰富地震信息的数字化地震波形资料,在不同类型地震的识别方面做了一定工作[4-5],本文在参考前人研究成果的基础上,从地震波形的时域、频域及震源参数等地震学特征入手,综合震区周围环境特征,探究2013年松原地区几次MS≥5地震的成因。

1 研究方法及数据 1.1 研究方法

由于震源体、传播路径、场地等的不同,台站接收到的地震波形信息十分丰富,通过对地震波进行分析,可对构造地震和因爆破、抽水、蓄水、采矿等人为因素引起的非天然地震加以区分。在台站j的地震仪记录的第i个地震在时间域里可表达为[6]

$ \begin{array}{l} U_{ij}^{obs}(t) = ({S_{ij}}(t) * P{}_{ij}{\rm{(}}t{\rm{)}} * \\ {L_j}{\rm{(}}t{\rm{)}} * {I_j}{\rm{(}}t{\rm{)}}) \cdot {K_j} \cdot Rij{\kern 1pt} \end{array} $ (1)

式中,$U_{ij}^{obs}(t)$为位移,Sij(t)为第i个地震的震源项,Pij(t)为第i个地震的震源与第j个台站之间的传播路径效应,Lj(t)为第j个台站的局部场地效应,Ij(t)为第j个台站的仪器响应,Kj为自由表面放大因子,Rij为地震波的辐射花样(取决于震源机制和台站相对于震源的位置),表示褶积。可通过傅里叶变换将式(1)在频率域中表达为:

$ \begin{array}{l} U_{ij}^{obs}(f) = {S_{ij}}(f) \cdot P{}_{ij}{\rm{(}}f{\rm{)}} \cdot {L_j}{\rm{(}}f{\rm{)}} \cdot \\ {I_j}{\rm{(}}f{\rm{)}} \cdot {K_j} \cdot \left| {{R_{ij}}} \right| \end{array} $ (2)

去除仪器响应后为:

$ U_{ij}^{obs}(f) = {S_{ij}}(f) \cdot P{}_{ij}{\rm{(}}f{\rm{)}} \cdot {L_j}{\rm{(}}f{\rm{)}} \cdot {K_j} \cdot \left| {Rij} \right| $ (3)

通过式(3)可以看出,若xy地震的震源足够接近,即两次地震震源的距离比震源距Rij小得多,则可认为台站j记录到的xy地震的观测位移谱UxjUyj只与各自的理论震源谱Sij(f)和辐射花样Rij中的震源机制部分有关,传播路径效应Pij(f)和台站的局部场地效应Lj(f)及自由表面放大系数Kj可以消除,也就是说,此时两次地震的观测位移谱可以反映各自的震源特征。

根据Brune模型[7],地震i的理论震源谱可表示为:

$ {S_i}{\rm{(}}f{\rm{)}} = \frac{{{\Omega _{{\rm{0}}i}}}}{{{\rm{1}} + {{{\rm{(}}f/{f_{Ci}}{\rm{)}}}^{\rm{2}}}}} $ (4)

式中,Ω0i为地震i的震源谱低频渐近线值,即零频极限值;fCi为低频渐近线和高频渐近线交点处的频率,即拐角频率。

对于位移震源谱的计算,本文首先截取地震波Sg到达至基本衰减结束为S波段,经带宽为0.2~22 Hz的4阶Butterworth带通滤波后,将S波段分段,对每小段扣除仪器响应,进行快速傅里叶变换,加汉宁窗以弥补能量损失,再按照邻段50%重叠的规则进行移窗。同样操作后得到每小段的傅里叶谱,叠加起来即得到S波段傅里叶谱。利用同样方法,对P波初动前256个采样点进行处理,得到噪声段的傅里叶谱,扣除噪声后将速度型数据除以2πf换算成位移型数据,即得到地震波的观测位移震源谱。根据Brune的ω-2震源模型[7],使拟合的理论位移震源谱和观测位移震源谱残差最小,即可得到零频极限值Ω0i和拐角频率fCi

结合地震的理论震源谱,通过计算2013年松原5次MS≥5地震的零频极限值Ω0i之间的相关系数分析这5次地震震源特性的相似性,若两次地震的震源特性相似,则相关系数rxy应较大,接近1;若震源特性不相似,则相关系数rxy较小[8]。地震xy的零频极限值相关系数rxy可表示为[9]

$ {r_x}y{\rm{ = }}\frac{{\sum\limits_{j = 1}^n {\sum\limits_{l = 1}^2 {({x_{jl}} - \bar x)({y_{jl}} - \bar y)} } }}{{\sqrt {\sum\limits_{j = 1}^n {\sum\limits_{l = 1}^2 {{{({x_{jl}} - \bar x)}^{\rm{2}}}} } } \sqrt {\sum\limits_{j = 1}^n {\sum\limits_{l = 1}^2 {{{({y_{jl}} - \bar y)}^{\rm{2}}}} } } }} $ (5)

式中,xjlyjl分别为台站j记录的地震xyl分量波形资料反演得到的零频极限值Ω0xjΩ0yjxy分别为xjlyjl的平均值,n为台站数目,l取1~2(分别为垂向和水平向S波)。此外,理论震源谱的拐角频率和高频衰减系数常用来区分不同成因的地震,本文计算了2013年5次松原地震和一些天然构造地震及塌陷地震的拐角频率和高频衰减系数,试图对2013年松原5次MS≥5地震的成因机制作出判断。

不同成因地震的震源特性不同,在地震波形时域分布特征及频率域的分布特征上都有所体现[10]。一般来说,在时域特征方面,塌陷地震较构造地震长周期面波发育、衰减较慢,波形比较简单;在频域上的差异体现在塌陷地震高频成分不发育,能量主要集中在低频段。在分析地震波形时域和频域的特征时,计算P波和S波的振幅比AP/AS及优势频率是一种常用的方法,但这些参数仅能反映某个时间点或某个频率点的幅值情况,具有较大的偶然性,分析结果不够准确。本文通过计算波形复杂度C及谱比值SR来反映地震波时域和频域的分布情况。波形复杂度C和谱比值SR的定义为[11]

$ {\rm{SR}} = \frac{{\int_{f = {L_1}}^{f = {H_1}} {\left| {X(f)} \right|} df}}{{\int_{t = {L_2}}^{t = {H_2}} {\left| {X(f)} \right|df} }} $ (6)
$ C = \frac{{\int_{t = L}^{t = H} {{x^2}} (t)dt}}{{\int_{t = {t_0}}^{t = L} {{x^2}} (t)dt}} $ (7)

式中,x(t)为地震信号t时刻的幅值,t0为P波到达台站的时刻,L为P波结束时刻,H为S波基本衰减完时刻,X(f)为地震信号在f频率处的幅值,L1H1L2H2分别为高频段和低频段的上、下限(本文采用的带限赫兹数为:L1=0.2 Hz,H1=1.2 Hz,L2=1.2 Hz,H2=6 Hz)。

地震的震源深度是地震学研究的核心问题之一,准确确定地震震源深度对于确定地震成因有非常重要的意义[12]。因此,本文根据震中距较小的流动台资料,利用Pg波、Sg波到时差的方法检验松原地震的震源深度。

1.2 数据选取

本文以2013年松原5次5级以上的地震事件作为研究对象,再选取两次构造地震事件及两次塌陷地震事件作为对比,具体信息见表 2。由于2013年松原地区5次5级以上地震震中位置接近,为尽量减少传播路径和场地对计算震源谱的影响,选用震中距较小的固定台站乾安台接收到的地震波数据。

表 2 研究序列的台站及仪器信息 Tab. 2 Information of station and seismometer
2 研究结果 2.1 震源谱特征分析

在位移震源谱中,拐角频率是对位移谱宽度的一种估计,反映低频与高频成分的分布特征,高频衰减系数则反映地震波中高频部分衰减的快慢程度,这两个参数基本上反映了震源体的特征[13]。震源谱参数受方位角的影响不大,因此选择震中距最小的固定台站波形数据进行震源参数的计算(表 3)。可以看出,2013年松原地区MS≥5.0地震的拐角频率和高频衰减系数较构造地震的小,与塌陷地震的较为接近。根据前人研究资料[5],构造地震拐角频率和高频衰减斜率比塌陷地震的大,因此从震源谱特征的角度来分析比较,松原地震的震源参数与典型构造地震的震源参数有一定差异,同典型构造地震相比,松原地震的震源谱具有衰减早、衰减快的特征,这种特征与塌陷地震更加类似。

表 3 研究序列的震源参数 Tab. 3 Source parameters of Songyuan swarm, tectonic earthquakes and collapsed earthquakes

此外,选择2013-10-31~11-23松原地区5次5级以上地震,根据乾安台(QAT)、长岭台(CLT)和松原台(SYT)等3个台站记录到的这5次地震的波形记录,反演各地震S波段垂向和水平向的零频极限值,并计算5次地震的相关系数,结果见表 4

表 4 2013年松原5次MS≥5地震零频极限值的相关系数 Tab. 4 Correlation coefficient of 2013 Songyuan MS≥5 earthquakes' spectral amplitude

表 4可知,在震源特性上,2013年松原5次MS≥5地震事件中,第1次和第5次地震事件比较相似,第2、3、4次地震事件比较相似。

2.2 时域特征分析

首先将松原地震波形与构造地震及塌陷地震的波形进行宏观上的比较(图 2),可以看出,松原地震的波形与构造地震的波形存在一定差异。一方面,松原地震衰减速度要比典型的构造地震衰减速度慢;另一方面,松原地震的波形周期较大,“毛刺”较少,主要是因为面波发育,频率成分比较单一,无高频波叠加,这些特征都与塌陷地震较为相似。通过计算波形复杂度(表 5)可以看出,塌陷地震的波形复杂度明显大于构造地震,而松原地震一些分量的波形复杂度又比塌陷地震大得多,说明虽然松原地震的波形在时域特征上与塌陷地震较为相似,但比典型的塌陷地震更为复杂。

图 2 松原地震波形与构造地震和塌陷地震波形的对比 Fig. 2 Waveform of Songyuan earthquake compared with tectonic earthquake and collapsed earthquake

表 5 研究序列的时域、频域参数 Tab. 5 Parameters of time domain and frequency domain
2.3 傅氏谱特征分析

一般塌陷地震较构造地震的震源深度浅,而地震波在松软的沉积层中传播时高频部分会被吸收,因此塌陷地震在频域中的特征也与构造地震有着明显的区别[4]。同样,将松原地震与构造地震及塌陷地震的傅氏谱进行宏观上的比较(图 3)可以发现,松原地震和塌陷地震的傅氏谱主要在低频段发育,优势频段集中在2 Hz以内;构造地震傅氏谱的优势频段较宽,大约为4 Hz以内,没有只集中在低频段的特点。对比地震的谱比值(表 5)发现,松原地震的谱比值明显小于构造地震的谱比值,与塌陷地震更为接近。由此可知,2013年松原地震与典型的构造地震明显不同,更加符合塌陷地震的特征。

2.4 震源深度检验

本文研究的5次松原地震中,只有2013-11-22发生的MS5.3地震有近震流动台站记录,因此选取此次地震事件的一个流动台数据对震源深度进行检验。当震中距较小时(小于1~2倍震源深度),可假设台站就位于震源的正上方,这样根据直达波法测得的震源距就可以看作震源深度的下限(真正的震源深度要浅于假定的震源深度)[14]。利用这种方法,选择震中距为2 km的流动台站L2202台记录到的地震波形,计算得到Pg波和Sg波的到时差Δt为2 s,结合当地速度模型[15],选择P波传播速度为vP=4 km/s,S波传播速度vS=1.32 km/s,根据式(8)计算得到震源深度的下限d为6.9 km,即震源深度不超过6.9 km,而不是8 km(表 1)。

$ {\rm{d = }}\frac{{{\rm{\Delta t}} \cdot {\rm{v_p}} \cdot {\rm{v_s}}}}{{{\rm{v_p - v_s}}}} $ (8)

地震勘探资料显示,松原市乾安地区的采油井深度多在3.8 km左右[2],松辽盆地中生代沉积底部、基岩顶层深度约为6 km[3]。此次地震的震源深度较浅,仍在油气圈闭的沉积盖层内。因此,在震源深度特征上,松原地区2013-11-22发生的MS5.3地震可能与该地区油气开采有关。由于近震台站数量较少,不能准确地对5次地震事件的震源深度进行精确检验,但综合地震的时域和频谱特性来看,5次地震呈现出一致性,且根据地震的零频极限频率的相关系数来看,第2、3、4次地震的震源特性具有相似性,大致推测松原地区地震群与该地区油气开采有关。若日后震区台站加密,对震源深度的计算可以更加准确,对地震动的时频特性和震源特性的认识会更加清晰[16]

图 3 傅氏谱对比 Fig. 3 Fourier spectrum comparison
3 结语

本文结合2013年松原地区MS≥5地震的台站记录信息,从不同方面对地震群的成因进行分析,得出以下结论:

1) 时空特征方面,震区附近历史上仅发生过两次5级以上地震,而2013年仅一个月内就发生5次5级以上地震,通过检验其中一次地震的震源深度发现,这种短暂突发、震群密集、震源浅的特点不符合该地区构造地震的特征。

2) 地震学参数方面,松原地震拐角频率为2.4 Hz左右,高频衰减斜率在-3.0左右,与构造地震相比呈现衰减早、衰减快的特征,而这种特征与塌陷地震类似。通过计算松原MS≥5地震零频极限频率的相关系数发现,5次地震中第1次和第5次地震的震源特性相似性较高,第2、3、4次地震的震源特性相似性较高。此外,松原地震的波形周期大、面波发育、频率成分单一,优势频率主要集中在2 Hz以内,通过对波形复杂度和谱比值的计算发现,松原地震明显不同于构造地震,而是更加符合塌陷型地震特征。

通过对2013年松原地区MS≥5地震的数字化波形时域、频域、震源参数进行分析,松原地震与典型构造地震有明显差别,其特征更加符合塌陷型地震。综合震区周边环境及震源深度推测,松原地震的成因与长期油气开采活动有关,但因近震台站数据及采油作业的具体资料有限,为进一步对地震成因机理进行分析,日后将加密震区观测,收集更多采油作业详细资料。

致谢: 感谢吉林省地震局和江西省地震局地震监测中心提供相关地震波形资料。

参考文献
[1]
李传友, 王一鹏, 张良怀, 等. 吉林省松原地区1119年6$\frac{3}{4}$级地震的发震构造条件[J]. 中国地震, 1999, 15(3): 237-246 (Li Chuanyou, Wang Yipeng, Zhang Lianghuai, et al. Causative Tectonic Conditions of the Historic Earthquake(M=6$\frac{3}{4}$) in 1119 in Songyuan Area[J]. Earthquake Research in China, 1999, 15(3): 237-246) (0)
[2]
高金哲, 李志伟, 包丰, 等. 2006年吉林乾安-前郭M5.0地震深度及其成因探讨[J]. 地球物理学进展, 2013, 28(5): 2328-2335 (Gao Jinzhe, Li Zhiwei, Bao Feng, et al. Resolving Focal Depth of the 31 March 2006 Jilin Qian'an-Qianguo Earthquake with Multiple Datasets and Its Implication for Seismogenesis[J]. Progress in Geophysics, 2013, 28(5): 2328-2335) (0)
[3]
刘双庆, 薛艳, 蔡宏雷, 等. 利用重测定的震源深度特征探讨2013年吉林前郭5.8级震群发震因素[J]. 地震研究, 2015, 38(2): 211-220 (Liu Shuangqing, Xue Yan, Cai Honglei, et al. Discussion on Cause Factor of Jilin Qianguo MS5.8 Earthquake Sequence in 2013 Using Focal Depth Characteristic of Relocation[J]. Journal of Seismological Research, 2015, 38(2): 211-220) (0)
[4]
张萍, 高艳玲, 肖健, 等. 辽宁台网记录爆破、矿震与地震的识别[J]. 地震地磁观测与研究, 2001, 22(5): 29-34 (Zhang Ping, Gao Yanling, Xiao Jian, et al. The Distinguish of Explosions Mine-Shocks and Earthquake Recording in Liaoning Seismic Network[J]. Seismological and Geomagnetic Observation and Research, 2001, 22(5): 29-34 DOI:10.3969/j.issn.1003-3246.2001.05.004) (0)
[5]
吴海波, 姚运生, 薛军蓉. 长江三峡地区地震事件频谱特征及分类问题探讨[J]. 大地测量与地球动力学, 2010, 30(6): 10-15 (Wu Haibo, Yao Yunsheng, Xue Junrong. Discuss of Frequency Spectrum Characteristics and Classified Problem of Seismic Event on Three Gorges[J]. Journal of Geodesy and Geodynamics, 2010, 30(6): 10-15) (0)
[6]
中国地震局监测预报司. 测震学原理与方法[M]. 北京: 地震出版社, 2017 (Department of Earthquake Monitoring and Prediction, CEA. Principles and Methods of Seismology[M]. Beijing: Seismological Press, 2017) (0)
[7]
Brune J N. Tectonic Stress and the Spectra of Seismic Shear Waves from Earthquakes[J]. Journal of Geophysical Research, 1970, 75(26): 4997-5009 DOI:10.1029/JB075i026p04997 (0)
[8]
张致伟, 乔慧珍, 吴朋, 等. 注水诱发地震的谱振幅相关系数及视应力研究[J]. 地震研究, 2015, 38(1): 42-50 (Zhang Zhiwei, Qiao Huizhen, Wu Peng, et al. Study on Correlation Coefficient of Spectral Amplitude and Apparent Stress of Water Injection Induced Earthquake[J]. Journal of Seismological Research, 2015, 38(1): 42-50 DOI:10.3969/j.issn.1000-0666.2015.01.006) (0)
[9]
Lund B, Bvarsson R. Correlation of Microearthquake Body-Wave Spectral Amplitudes[J]. Bulletin of the Seismological Society of America, 2002, 92(6): 2419-2433 DOI:10.1785/0119990156 (0)
[10]
林怀存, 王保平, 刘洪瑞, 等. 构造地震与塌陷地震对比研究[J]. 地震学报, 1990, 12(4): 448-455 (Lin Huaicun, Wang Baoping, Liu Hongrui, et al. Comparative Study of Tectonic and Collapse Earthquakes[J]. Acta Seismologica Sinica, 1990, 12(4): 448-455) (0)
[11]
刘代志, 李夕海. 核爆地震模式识别[M]. 北京: 国防工业出版社, 2010 (Liu Daizhi, Li Xihai. Pattern Recognition of Nuclear Explosion Earthquake[M]. Beijing: National Defense Industry Press, 2010) (0)
[12]
郑勇, 谢祖军. 地震震源深度定位研究的现状与展望[J]. 地震研究, 2017, 40(2): 167-175 (Zheng Yong, Xie Zujun. Present Status and Prospect of Earthquake Focal Depth Locating[J]. Journal of Seismological Research, 2017, 40(2): 167-175 DOI:10.3969/j.issn.1000-0666.2017.02.001) (0)
[13]
陈运泰, 吴忠良. 数字地震学[M]. 北京: 地震出版社, 2000 (Chen Yuntai, Wu Zhongliang. Broadband Seismology[M]. Beijing: Seismological Press, 2000) (0)
[14]
申学林, 王秋良, 魏贵春, 等. 基于近台数据的2017-06巴东M4.3地震序列震源深度研究[J]. 大地测量与地球动力学, 2018, 38(1): 14-17 (Shen Xuelin, Wang Qiuliang, Wei Guichun, et al. Study on Focal Depth of Badong M4.3 Earthquake Sequence in June 2017[J]. Journal of Geodesy and Geodynamics, 2018, 38(1): 14-17) (0)
[15]
吴微微, 杨建思, 苏金蓉, 等. 2013年吉林前郭-乾安震源区中强地震矩张量反演与区域孕震环境研究[J]. 地球物理学报, 2014, 57(8): 2541-2554 (Wu Weiwei, Yang Jiansi, Su Jinrong, et al. Moment Inversion of Moderate Earthquakes and Seismogenic Environment in Qianguo-Qian'an Source Region, 2013, Jilin Province[J]. Chinese Journal of Geophysics, 2014, 57(8): 2541-2554) (0)
[16]
孟庆筱, 吕健, 景鹏旭. 基于二维有限单元法的庐山M7[J]. 大地测量与地球动力学, 2019, 39(2): 131-136 (Meng Qingxiao, Lü Jian, Jing Pengxu. Simulation of Near Field Strong Ground Motion of Lushan M7.0 Earthquake Based on Two-Dimension Finite Element Method[J]. Journal of Geodesy and Geodynamics, 2019, 39(2): 131-136) (0)
Study on Genesis of MS≥5 Songyuan Earthquake Swarm in Jilin Province in 2013
YU Guangcheng1     YAO Yunsheng1,2     ZHANG Lifen1     WANG Qiuliang1     WU Haibo1     
1. Key Laboratory of Earthquake Geodesy, Institute of Seismology, CEA, 40 Hongshance Road, Wuhan 430071, China;
2. Institute of Disaster Prevention, 465 Xueyuan Street, Yanjiao High Tech Zone, Sanhe 065201, China
Abstract: A series of earthquakes occurred in Songyuan area, Jilin province between October 31 and November 23, 2013, including five MS≥5 earthquakes. The largest one, of MS5.8, occurred November 23, 2013. By calculating the correlation coefficients of the spectral amplitude of these five MS≥5 earthquakes in Songyuan area in 2013, it is found that the focal mechanisms of the first and fifth earthquakes are similar, while those of the second, third and fourth earthquakes are similar. Statistical analysis on these five earthquakes and compared with tectonic earthquakes and collapsed earthquakesshows that the displacement source spectrum of Songyuan swarm in 2013 presents the characteristics of early attenuation and fast attenuation compared with tectonic earthquakes. Characteristic parameters of time-domain and frequency-domain, such as waveform complexity and spectral ratio, are also quite different from those of tectonic earthquakes, but more in line with the characteristics of collapsed earthquakes.
Key words: tectonic earthquake; collapsed earthquake; seismogenesis; spectral analysis; source spectrum