Research on the Repositioning and Seismic Source Characters of the Changdao Coastal Earthquake Sequences, Shandong Province
-
摘要: 基于山东数字化地震台站的震相资料和胶东半岛地区的速度结构模型,利用双差定位方法对2017年2月14日开始的山东长岛地区双震群进行重新定位,结果显示主震群余震序列分布沿一级断裂分布,震群位置比较集中。并利用CAP方法(Cut and Paste)对长岛海域2个震群中的几个较大地震进行震源机制解的反演分析,其中,2017年3月3日的主震震源机制解的节面Ⅰ走向为320°,倾角57°,滑动角为12°;节面Ⅱ走向为223.4°,倾角80°,滑动角为146.4°;主震的最佳拟合深度在9.7km,6次较大地震的节面Ⅰ走向也基本一致。Abstract: Based on the data collected by the Shandong Digital seismic Network and the speed-structure model, we use the accurate double-difference location method to relocate two sequences starting from the February 14, 2017, the result show that the distribution of the aftershocks of the main sequence is along the class-Ⅰ fault in earthquake swarm, In this paper also applied the CAP method to analyze inversely the focal mechanism solutions of several earthquake in the Changdao costal area. The focal mechanism solution of the main earthquake on March 3, 2017 includes the flowing parameters:nodal plane Ⅰ (trend of 320°, dip angle 57° and slide angle 12°); plane Ⅱ (trend of 223.4°, dip angle 80° and slide angle 146.4), The optimum fitting depth is 9.7km The solutions of nodal plane of these six large earthquakes agree with each other well.
-
Key words:
- Changdao /
- Accurate double-difference location /
- CAP /
- Focal mechanism solution
-
引言
长岛是山东省唯一的海岛县,包括南北长山岛在内的10个居民岛屿,地处构成渤海和黄海分界线的庙岛群岛的南端。长岛位于华北地台区,与胶辽地块其它地区具有类似的地质构造演化历史,由于受到中国东部新生代扩展断陷活动以及第四纪海侵共同影响而成为了独立的岛屿。长岛位于NW向张家口-渤海-威海断裂带的东段,即蓬莱-威海断裂带内,带内多条晚更新世活动断裂均从长岛附近海域通过。同时在其西侧还发育了晚更新世以来强烈活动的郯庐断裂带,与活动构造的分布特点相吻合,长岛附近地震活动水平较高,现代中小震密集分布。值得注意的是,在附近海域曾发生过1548年渤海海峡7级地震和1948年威海6级地震。因此,有必要对长岛及其附近海域震群进行重定位和震源特征研究。
长岛县及附近海域发育有NNE、NE、NW等多个方向的断裂。长岛区内发育多组断裂,并被断裂围限成大小不等的块体。NW断裂也是以蓬莱-张家口断裂为主,长岛-芝罘岛和蓬莱-张家口断裂等多条NW、近EW向断裂组合成1条倾角较大的断裂带,不同方向的断裂交汇造成该区域地震活动水平较高。
2017年2月14日山东台网记录到长岛海域第一个地震以来,又陆续记录到余震3000余次。由于此次震群序列属于海域地震,所以在密集台阵的布设上有一定困难,只有3个海岛台,此次震群附近台站的分布和附近断裂情况,如图 1所示。
1. 定位方法及资料选取
本文采用的双差定位(Double-Difference)方法由Waldhauser等(2000)提出,该方法是地震定位学中的1种相对定位方法,是主事件方法的应用推广,但在震群相对位置的反演过程中,无需把震群的某个地震作为主事件,而是震群的某一地震和震群矩心的相对位置进行反演,它在程序参数的选择上有一定的要求。该方法在国际和国内已经被广泛使用,Waldhauser等(2000)利用双差定位算法在北加州北海沃德断层上的2个震群上进行了应用,利用震群矩心200km内的台站和10km范围内的地震,计算了地震对之间的走时残差,使得震群在垂直的条带分布更加集中;周龙泉等(2003)利用双差定位方法对云南大姚双主震序列进行了重新定位研究;黄媛(2008)对汶川地震及其余震进行重定位;张勇(2014)、张广伟等(2014)对云南鲁甸地区进行重新定位并通过震源机制分析其地震活动构造意义;李铂(2017)对山东乳山地区进行重新精定位并分析其断层活动性;罗佳宏(2017)利用双差定位方法对三峡水库的地震活动进行了研究。
双差定位方法的优点主要是利用谱域中的波形互相关技术提高定位质量和定位精度,相对于主事件方法,数据量大大增加,在一定的范围内进行相近地震事件对的配对,要求配对的地震事件在一定空间范围至少可以搜索到至少1次及以上的地震事件,故台网的定位质量相比会有一定提高。双差定位方法的基本原理为:相对于同1个地震台站k,i和j分别代表不同的地震震源,2个震源的地震波观测到时和理论到时之间的残差为rki和rkj,Tki是第1个震源i到地震台站k的体波观测到时,ti是发震时间,u是波的慢度,ds代表地震波沿射线路径的线元(Helmberger, 1980, 2000)。
对2个地震事件组i和j的残差做差,即双重残差:
$$ {r_{kj}}\;{r_{ki}} = \left({\sum\limits_{i = 1}^3 {\frac{{\partial {T_{kl}}}}{{\partial {x_j}}} + \Delta {t_j} + \int_l^k {\delta u{\rm{d}}s} } } \right)\left({\sum\limits_{i = 1}^3 {\frac{{\partial {T_{ki}}}}{{\partial {x_i}}}} + \Delta {t_j} + \int_i^k {\delta u{\rm{d}}s} } \right) $$ (1) 公式(1)即为双差定位方法的联合反演公式。
对于距离相近的事件,可以忽略2个事件震中间的速度(慢度)变化,即$\int_l^k {\delta u{\rm{d}}s} $和$\int_i^k {\delta u{\rm{d}}s} $可以相互抵消,得到:
$$ {r_{kj}}\;{r_{ki}} = \left({\sum\limits_{i = 1}^3 {\frac{{\partial {T_{kl}}}}{{\partial {x_j}}} + \Delta {t_j}} } \right)\left({\sum\limits_{i = 1}^3 {\frac{{\partial {T_{ki}}}}{{\partial {x_i}}}} + \Delta {t_j}} \right) $$ (2) 利用式(2)进行震中反演。对于式(2),由于忽略了$\int_l^k {\delta u{\rm{d}}s} $和$\int_i^k {\delta u{\rm{d}}s} $,即认为事件对之间不存在速度变化,这就要求事件对之间的距离相隔不能太远,因此在实际的地震定位工作中,选择合适的事件对之间的最大距离尤为重要。该距离上限太小,则组成的不同事件对数据太少,反之,则会引入较大的结构速度误差。
山东地震台网自从2006年开始进行数字化改造,又陆续对台网进行网络升级、台站改造,目前已经有133个数字化测震台站。在长岛海域震群附近50km范围内有6个台站,由于该震群是比较少见的海域震群,而在海域台站的布设方面有一定困难,造成了部分区域的观测空区。
山东省烟台市长岛县海域自2017年2月14日开始发生地震以来,3月3日又发生了4.5级地震,之后又发生了一系列的地震震群活动,9月2日,在该震群南方10km左右又出现了1个新的震群活动。截至2018年1月1日,山东数字化台网已经记录到近2700次余震,其中1级以下余震2175次,1—2级余震435次,2.0—2.9级余震73次,3.0—3.9级余震4次,4级以上余震1次(图 2)。
在2017年9月出现的主震群南侧的小震群也出现了200余次余震,其中最大地震为2017年3月3日的4.4级地震,这也是山东2017年记录到境内及周边海域最大的地震。
截至2018年1月,此次长岛地震的震群活动一直持续,从中选择了2711次有3个及以上台站记录的地震,最终参与双差定位的地震包括4个台站记录到的地震1002个,5个台站记录到的地震775个,6个台站记录到的地震681个,其中2017年3月3日长岛海域4级地震有48个地震台站参与定位。在重定位时,由于残差原因,最终选择了参与定位地震数量为2581次,包括11232次P波和S波的到时记录定位残差也从0.172下降到了0.133。为了验证数据的时效性,绘制了几个震相的走时曲线,如图 3所示,其中Pg、Sg、Pn和Sn的走时曲线离线性均较低,可以说明地震数据质量比较稳定。
根据折合走时曲线图(图 4)可以发现,Pg和Pn的数据点与理论折合走时曲线较为符合,在一定程度上反映了当前速度模型在本研究区域的正确性,可以将当前速度模型作为研究的基础数据资料。
采用目前山东地震数字化台网正在使用的一维速度模型(李铂等,2012),如表 1所示。
表 1 地壳速度模型Table 1. Crustal velocity model深度/km vp/km·s-1 vs/km·s-1 22 6.13 3.54 33 6.88 3.98 70 7.93 4.58 图 5显示了长岛海域震群的地震频次和定位深度的分布,图 5(a)为双差定位之前的台网记录结果,可以看出在6—18km地震频次比较集中,但整体分布比较分散;图 5(b)为定位之后的深度结果,深度优势分布在8—18km,其中10—16km范围更为明显。
图 6为利用双差定位方法对长岛2个震群重新定位后的结果分布,由图可以看出,蓝色点为台网原始记录,地震整体分布比较分散,没有明显的条带状分布,经过重新定位后,2个震群的外围发散地震向内收敛明显,震中分布更加集中,优势分布也更加明显。3月3日发生的北侧震群(以下称主震群)走向与一级断裂——大竹山岛-威海北断裂走向基本一致,断裂全长70km,且靠近断裂的西北端口处,余震呈NWW向分布,A’B’的主破裂长度约在17km,其共轭方向的C’D’破裂长度约6km,几次较大的地震都发生在该破裂轴上,4级以上地震4次(图中红白球标识),4级以下地震2次(图中蓝白球标识);9月2日发生的主震群南侧13km的震群(以下称副震群),地震个数较少,但也基本看出走向与主震群走向基本一致,在2个二级断裂之间,不排除存在隐伏断层的可能性,震源深度分布略低于主断层,优势深度集中在6—10km范围内,这也与副震群靠近陆地有一定关系,在一定程度上表明该深度的应力介质强度多局限于这一范围。
2. 震源机制解
采用CAP(Cut and Paste)方法反演了2017年3月3日长岛海域M 4.4地震的震源机制解由于震群属于海域型震群,在对其进行反演的过程中使用的台站较少,而且相对于震中,台站方位角的覆盖性也较差,所以利用CAP方法对其进行震源机制解反演,其优势之一是在台站利用率较低的情况下可以得到较好的震源机制解结果。
在反演长岛海域震群的记录时,一般对波形记录进行权重设置,对比前期对胶东半岛其它震群设置的权重比例(一般设置为2:1),对首波和长周期的面波分别选取0.07—0.15Hz和0.04—0.08Hz频段。在反演过程中对观测波形进行滤波,扣除仪器响应和波形旋转,反演后的最佳拟合波形如图 7(a)所示,黑色波形曲线为观测值,红色波形曲线为理论值,左侧为台站代码、震中距以及该台理论P波初至与观测P波初至的差值,拟合波形的下方为相对应的时间偏移量,5个数值代表台站的3分量5个分向(体波垂向(PV)、体波径向(PR)、面波垂向(Surf.V)、面波径向(Surf.R)、面波切向(SH))。由图可见,所选台站的拟合系数均较高,基本在85%以上,这也客观地证明了结论的正确性。图 7(b)给出了此次地震最佳的拟合深度,从图中看出对应深度拟合曲线中的斜率最大点处的深度为最佳深度,该深度约9.7km,误差为0.2km,也与双差定位测定的优势深度相一致。
本文反演了长岛海域震群6次较大的震源机制解(表 2,图 7)。由表 2可以看出,长岛几次较大地震都发生在主震群范围内,表中定位结果采用双差定位,几次较大余震的破裂位置均在3km范围内,几次地震的主压应力轴方向基本一致,均以走滑型为主,节面Ⅰ走向也较一致,几次较大地震均发生在主震群所在区域,节面Ⅰ走向也与一级断裂竹山岛-威海北断裂走向基本一致,其中2017年3月3日的主震震源机制解的节面Ⅰ走向为320°,倾角57°,滑动角为12°;节面Ⅱ走向为223.4°,倾角80°,滑动角为146.4°;几次较大地震的最佳拟合深度为8.7—13.7km。
表 2 长岛海域震群序列6次较大地震震源机制解Table 2. Focal mechanisms of six larger earthquake in the Changdao earthquake sequences发震时间 东经/° 北纬/° 震级 节面Ⅰ/° 节面Ⅱ/° 最佳拟合
震源深度/km走向 倾角 滑动角 走向 倾角 滑动角 2017-03-03
02:48120.8056 38.1051 4.4 320 57 12 223.4 80 146.4 9.7 2017-03-21
00:54120.7839 38.1014 3.9 321 53 -10 57.1 82 -142.6 10.1 2017-03-27
22:06120.7721 38.1050 4.1 321 50 161 63.5 75.6 41.6 8.7 2017-04-08
02:08120.8051 38.1032 4.0 323 63 0 233 90 153 9.8 2017-04-08
07:44120.7897 38.1012 3.7 351 60 -57 118.6 43.4 -133.3 13.7 2017-06-18
16:33120.7899 38.0955 4.0 331 55 -12 61 90 -143.2 10.7 3. 结论与讨论
采用双差定位方法对山东长岛附近海域2次可定位震群的2711次地震事件进行了重定位,由于海域震群运用定位方法的实例较少,所以在附近台站数量和方位角的分布上存在一定局限,通过研究发现重定位后的2个震群地震向内收敛明显,震中分布更加集中,优势分布也更加明显,更清晰地明确了条带分布的基本特性。余震呈NWW向展布,破裂长度约为17km,其共轭方向的破裂长度约为6km;主震群的震源位置与一级断裂大竹山岛—威海北断裂走向具有较高的一致性,副震群震例较少,但也能看出走向与主断层的走向基本一致。由于副震群周围没有较大的地质断裂结构作为破裂依据,也不排除有较小隐伏断层的存在,主副震群在定位后都有一定的收敛,没有明显的震群交集;2个震群的优势深度为6—15km,也与山东内陆及近海海域的地震发震深度范围相一致。
利用CAP方法反演了2017年3月3日长岛海域M 4.4地震的震源机制解,在震源深度9.7km附近对应的双力偶解为最佳的震源机制解。其余几次较大地震震中位置与震群主震的震中位置一致,震源机制解中节面Ⅰ的走向也具有一致性。
综合长岛海域2次震群的定位结果、几次较大地震震源机制解反演结果和区域地震资料,可以看出,长岛震群几次显著地震的震源机制解较为一致,均表现为走滑机制,其中NW向的节面与NWW向的张家口-蓬莱断裂带的走向一致。同时,从双差定位的结果也可看出,长岛震群表现出出现NW向的优势展布特征,进一步印证了发震断层可能是NW向的假设。因此,初步判定长岛震群的发震构造可能是NW向的张家口-蓬莱断裂的分支断裂。
山东长岛海域2次震群都处于山东近海,2个震群的活动性也在逐步减弱,但目前仍有小的地震持续发生,在后续的工作中,还需要进一步收集地震和地质断裂构造数据,以判断此次震群的后期持续性和断层的破裂特征。
致谢: 感谢Waldhauser F.提供的双差定位程序、崔仁胜提供的CAP程序以及评审专家给出的宝贵建议。 -
表 1 地壳速度模型
Table 1. Crustal velocity model
深度/km vp/km·s-1 vs/km·s-1 22 6.13 3.54 33 6.88 3.98 70 7.93 4.58 表 2 长岛海域震群序列6次较大地震震源机制解
Table 2. Focal mechanisms of six larger earthquake in the Changdao earthquake sequences
发震时间 东经/° 北纬/° 震级 节面Ⅰ/° 节面Ⅱ/° 最佳拟合
震源深度/km走向 倾角 滑动角 走向 倾角 滑动角 2017-03-03
02:48120.8056 38.1051 4.4 320 57 12 223.4 80 146.4 9.7 2017-03-21
00:54120.7839 38.1014 3.9 321 53 -10 57.1 82 -142.6 10.1 2017-03-27
22:06120.7721 38.1050 4.1 321 50 161 63.5 75.6 41.6 8.7 2017-04-08
02:08120.8051 38.1032 4.0 323 63 0 233 90 153 9.8 2017-04-08
07:44120.7897 38.1012 3.7 351 60 -57 118.6 43.4 -133.3 13.7 2017-06-18
16:33120.7899 38.0955 4.0 331 55 -12 61 90 -143.2 10.7 -
黄媛, 吴建平, 张天中等, 2008.汶川8.0级大地震及其余震序列重定位研究.中国科学D辑:地球科学, 38(10):1242-1249. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgkx-cd200810007 李铂, 崔鑫, 叶庆东等, 2012.山东地区地震波一维速度模型研究.华北地震科学, 30(4):1-6. http://d.old.wanfangdata.com.cn/Periodical/hbdzkx201204001 李铂, 崔鑫, 娄世平等, 2017.山东乳山地区震群特征及发震背景再研究.震灾防御技术, 12(3):491-500. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20170305&journal_id=zzfyjs 罗佳宏, 马文涛, 李春政, 2017.结合波形互相关的双差定位方法在三峡水库地震中的应用.震灾防御技术, 12(1):56-67. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20170106&journal_id=zzfyjs 张广伟, 雷建设, 梁姗姗等, 2014.2014年8月3日云南鲁甸MS 6.5级地震序列重定位与震源机制研究.地球物理学报, 57(9):3018-3027. http://www.cnki.com.cn/Article/CJFDTotal-DQWX201409027.htm 张勇, 许力生, 陈运泰等, 2014.2014年8月3日云南鲁甸MW 6.1(MS 6.5)地震破裂过程.地球物理学报, 57(9):3052-3059. http://www.cnki.com.cn/Article/CJFDTotal-DQWX201409031.htm 周龙泉, 刘杰, 马宏生等, 2009.2003年大姚6.2级、6.1级地震序列震源位置及震源区速度结构的联合反演.地震, 29(2):12-24. doi: 10.3969/j.issn.1000-3274.2009.02.002 Helmberger D., Ni S. D., Wen L. X., et al., 2000. Seismic evidence for ultralow-velocity zones beneath Africa and Eastern Atlantic. Journal of Geophysical Research:Atmospheres, 105 (B10):23865-23878. doi: 10.1029/2000JB900143 HelmbergerD. V., Engen G. R., 1980. Modeling the long-period body waves from shallow earthquakes at regional ranges. Bulletin of the Seismological Society of America, 70 (5):1699-1714. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=Y7gZYeE4k5gj7wltB7IsRx6NC8fESoJSDY09ucXIj+E= Waldhauser F., Ellsworth W. L., 2000. A double-difference earthquake location algorithm:method and application to the northern Hayward Fault, California. Bulletin of the Seismological Society of America, 90 (6):1353-1368. doi: 10.1785/0120000006 期刊类型引用(2)
1. 张正帅,周晨,郑建常,刘承雨. 长岛震群震源谱参数反演和震源参数特征. 地震. 2022(01): 85-98 . 百度学术
2. 李霞,陈时军,张正帅,戴宗辉,李小晗,卢仲斌. 山东庙岛群岛地区P波三维速度结构反演与2017年震群的发震构造分析. 地震地质. 2020(05): 1188-1204 . 百度学术
其他类型引用(0)
-