荐言献策

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

2020年北京顺义震群发震机理

宋程 马婷 徐小远 谭毅培

宋程,马婷,徐小远,等. 2020年北京顺义震群发震机理[J]. 华北地震科学,2022, 40(1):71-76. doi:10.3969/j.issn.1003−1375.2022.01.012
引用本文: 宋程,马婷,徐小远,等. 2020年北京顺义震群发震机理[J]. 华北地震科学,2022, 40(1):71-76. doi:10.3969/j.issn.1003−1375.2022.01.012
SONG Cheng,MA Ting,XU Xiaoyuan,et al. Mechanism Analysis of the 2020 Beijing Shunyi Earthquake Swarm[J]. North China Earthquake Sciences,2022, 40(1):71-76. doi:10.3969/j.issn.1003−1375.2022.01.012
Citation: SONG Cheng,MA Ting,XU Xiaoyuan,et al. Mechanism Analysis of the 2020 Beijing Shunyi Earthquake Swarm[J]. North China Earthquake Sciences,2022, 40(1):71-76. doi:10.3969/j.issn.1003−1375.2022.01.012

2020年北京顺义震群发震机理

doi: 10.3969/j.issn.1003-1375.2022.01.012
基金项目: 2021年度震情跟踪定向任务(2021010128);天津市地震局局内科研项目(yb202001);地震科技星火计划项目(XH21002Y)
详细信息
    作者简介:

    宋程(1992—),女,黑龙江牡丹江人,工程师,主要从事地震预报与地震活动性相关工作. E-mail: nibushiwode798@163.com

    通讯作者: 谭毅培(1982—),男,天津人,高级工程师, 主要从事地震预测预报和地震序列分析工作. E-mail: oivertan921@sina.cn
  • 中图分类号: P315.6

Mechanism Analysis of the 2020 Beijing Shunyi Earthquake Swarm

  • 摘要: 以2020年3—4月北京顺义小震群为研究对象,进行模板匹配遗漏地震检测、波形互相关震相检测、地震精定位、震源机制计算、破裂半径和应力降等震源参数计算,以探究其发震机理。计算结果显示:① 顺义震群发震构造为NW走向,震中呈现由NW向SE的迁移特征;② 除了震群中第一个ML2.4地震无法计算外,另外2次ML2.0以上地震的破裂面均处于之前地震产生的静态应力加载区内;③ 未检测到震群中存在重复地震活动。结合前人提出的震群发震机理模型,分析认为随机性模型可以较好地解释以上计算结果,顺义震群的发震机理可以解释为地震间的应力触发作用,此研究结果可为北京地区地震危险性分析提供科学参考。
  • 图  1  顺义震群震中及周边台站分布图

    注:红色实心圆为顺义震群震中位置;三角形为台站位置,蓝色三角表示参与模板匹配地震检测和精定位计算的台站,黑色实心三角表示仅参与地震精定位的台站,空心三角表示本文未使用的台站;灰色实心圆为2010年以来周边ML1.0以上地震震中位置

    图  2  6个ML≥1.0地震在东三旗台(DSQ)的记录波形比较图

    注:图中标号与表1中序号一致;以1号地震为模板,对6个地震P波和S波到时进行校正,灰色虚线表示对齐后震相到时

    图  3  顺义震群ML1.0以上地震精定位与ML2.5地震震源机制结果图

    注:红色圆点为精定位后的地震震中;黑色十字为台网给出的震中位置;放大图中标号与表1中序号一致,表示发震的先后顺序;左下角震源球节面图给出了求解过程中用到的台站以及初动方向

    图  4  断层面静态应力触发计算结果

    注:红色表示静态应力加载,蓝色表示卸载;黑色圆圈表示基于圆盘模型计算得到的地震的破裂面

    表  1  顺义震群ML1.0以上地震信息表

    序号发震时刻震中经度/
    (°E)
    震中纬度/
    (°N)
    深度/
    km
    震级
    ML
    12020-03-29 16:09:33.3116.62340.16632.4
    22020-03-29 17:36:21.5116.62540.16621.4
    32020-04-01 09:46:17.3116.62040.17862.5
    42020-04-09 08:08:17.3116.61640.17442.2
    52020-04-09 10:00:48.9116.60740.18631.0
    62020-04-09 11:24:21.2116.62440.17231.1
    下载: 导出CSV

    表  2  地震精定位所使用速度结构

    123456
    深度范围/km0~22~44~1515~3232~38≥38
    P波速度/(km/s)3.75.56.16.56.88.0
    下载: 导出CSV

    表  3  顺义震群ML1.0以上地震震源参数计算结果

    序号发震时刻震级/ML地震矩
    M0/(N·m)
    破裂半径
    R / m
    平均滑移量
    D / mm
    应力降
    Δσ/MPa
    12020-03-29 16:09:33.32.42.62×1012105.562.490.97
    22020-03-29 17:36:21.51.42.50×101153.840.920.70
    32020-04-01 09:46:17.32.53.31×1012112.912.761.01
    42020-04-09 08:08:17.32.21.64×101292.262.040.91
    52020-04-09 10:00:48.91.09.77×101041.130.610.62
    62020-04-09 11:24:21.21.11.24×101143.990.680.64
    下载: 导出CSV
  • [1] Vidale J E, Shearer P M. A survey of 71 earthquake bursts across southern California: exploring the role of pore fluid pressure fluctuations and aseismic slip as drivers[J]. Journal of Geophysical Research:Solid Earth, 2006, 111(B5): B05312.
    [2] Gomberg J. Unsettled earthquake nucleation[J]. Nature Geoscience, 2018, 11(7): 463-464. doi:  10.1038/s41561-018-0149-x
    [3] Kato A, Igarashi T, Obara K. Detection of a hidden Boso slow slip event immediately after the 2011 MW9.0 Tohoku-Oki earthquake, Japan[J]. Geophysical Research Letters, 2014, 41(16): 5868-5874. doi:  10.1002/2014GL061053
    [4] Shelly D R, Taira T, Prejean S G, et al. Fluid-faulting interactions: fracture-mesh and fault-valve behavior in the February 2014 Mammoth Mountain, California, earthquake swarm[J]. Geophysical Research Letters, 2015, 42(14): 5803-5812. doi:  10.1002/2015GL064325
    [5] Kato A, Fukuda J, Kumazawa T, et al. Accelerated nucleation of the 2014 Iquique, Chile MW 8.2 earthquake[J]. Scientific Reports, 2016, 6(1): 24792. doi:  10.1038/srep24792
    [6] Savage H M, Keranen K M, Schaff D P, et al. Possible precursory signals in damage zone foreshocks[J]. Geophysical Research Letters, 2017, 44(11): 5411-5417. doi:  10.1002/2017GL073226
    [7] Helmstetter A, Sornette D. Foreshocks explained by cascades of triggered seismicity[J]. Journal of Geophysical Research:Solid Earth, 2003, 108(B10): 2457.
    [8] Ruhl C J, Abercrombie R E, Smith K D, et al. Complex spatiotemporal evolution of the 2008 MW4.9 Mogul earthquake swarm (Reno, Nevada): interplay of fluid and faulting[J]. Journal of Geophysical Research:Solid Earth, 2016, 121(11): 8196-8216. doi:  10.1002/2016JB013399
    [9] Seif S, Zechar J D, Mignan A, et al. Foreshocks and their potential deviation from general seismicity[J]. Bulletin of the Seismological Society of America, 2019, 109(1): 1-18. doi:  10.1785/0120170188
    [10] Ellsworth W L, Bulut F. Nucleation of the 1999 Izmit earthquake by a triggered cascade of foreshocks[J]. Nature Geoscience, 2018, 11(7): 531-535. doi:  10.1038/s41561-018-0145-1
    [11] Yoon C E, Yoshimitsu N, Ellsworth W L, et al. Foreshocks and mainshock nucleation of the 1999 MW7.1 Hector Mine, California, earthquake[J]. Journal of Geophysical Research:Solid Earth, 2019, 124(2): 1569-1582. doi:  10.1029/2018JB016383
    [12] Wang J, Main I G, Musson R M W. Earthquake clustering in modern seismicity and its relationship with strong historical earthquakes around Beijing, China[J]. Geophysical Journal International, 2017, 211(2): 1005-1018. doi:  10.1093/gji/ggx326
    [13] Brodsky E E. The importance of studying small earthquakes[J]. Science, 2019, 364(6442): 736-737. doi:  10.1126/science.aax2490
    [14] Peng Z G, Zhao P. Migration of early aftershocks following the 2004 Parkfield earthquake[J]. Nature Geoscience, 2009, 2(12): 877-881. doi:  10.1038/ngeo697
    [15] Waldhauser F, Ellsworth W L. A Double-difference Earthquake location algorithm: method and application to the Northern Hayward Fault, California[J]. Bulletin of the Seismological Society of America, 2000, 90(6): 1353-1368. doi:  10.1785/0120000006
    [16] 祝治平, 张先康, 张建狮, 等. 北京—怀来—丰镇剖面地壳上地幔构造与速度结构研究[J]. 地震学报, 1997, 19(5): 499-505.
    [17] 王峻, 刘启元, 陈九辉, 等. 根据接收函数反演得到的首都圈地壳上地幔三维S波速度结构[J]. 地球物理学报, 2009, 52(10): 2472-2482. doi:  10.3969/j.issn.0001-5733.2009.10.006
    [18] Snoke J A, Munsey J W, Teague A C, et al. A program for focal mechanism determination by combined use of polarity and SV-P amplitude ratio data[J]. Earthquake Notes, 1984, 55(3): 15.
    [19] 赵翠萍, 陈章立, 华卫, 等. 中国大陆主要地震活动区中小地震震源参数研究[J]. 地球物理学报, 2011, 54(6): 1478-1489. doi:  10.3969/j.issn.0001-5733.2011.06.007
    [20] Andrews D J. A stochastic fault model: 1. Static case[J]. Journal of Geophysical Research:Solid Earth, 1980, 85(B7): 3867-3877. doi:  10.1029/JB085iB07p03867
    [21] 仲秋, 申文豪, 史保平. 近断层应力变化快速计算方法及对强余震空间分布的指示意义: 以MW 7.9汶川地震为例[J]. 地球物理学报, 2015, 58(1): 172-183. doi:  10.6038/cjg20150115
    [22] 王莹. 地震断层参数和投影断层参数对库仑破裂应力变化分布的影响[D]. 昆明: 云南大学, 2011.
    [23] Nadeau R M, Johnson L R. Seismological studies at Parkfield VI: moment release rates and estimates of source parameters for small repeating earthquakes[J]. Bulletin of the Seismological Society of America, 1998, 88(3): 790-814.
    [24] 谭毅培, 邓莉, 曹井泉, 等. 2015年河北滦县震群发震机理分析[J]. 地球物理学报, 2016, 59(11): 4113-4125. doi:  10.6038/cjg20161115
    [25] Horálek J, Fischer T. Role of crustal fluids in triggering the west Bohemia/Vogtland earthquake swarms: just what we know (a review)[J]. Studia Geophysica et Geodaetica, 2008, 52(4): 455-478. doi:  10.1007/s11200-008-0032-0
    [26] Borghi A, Aoudia A, Javed F, et al. Precursory slow-slip loaded the 2009 L’Aquila earthquake sequence[J]. Geophysical Journal International, 2016, 205(2): 776-784. doi:  10.1093/gji/ggw046
    [27] Jansen G, Ruhl C J, Miller S A. Fluid pressure-triggered foreshock sequence of the 2008 Mogul earthquake sequence: Insights from stress inversion and numerical modeling[J]. Journal of Geophysical Research:Solid Earth, 2019, 124(4): 3744-3765. doi:  10.1029/2018JB015897
    [28] Ross Z E, Cochran E S, Trugman D T, et al. 3D fault architecture controls the dynamism of earthquake swarms[J]. Science, 2020, 368(6497): 1357-1361. doi:  10.1126/science.abb0779
  • [1] 苏珊, 阎春恒, 向巍.  广西靖西MS5.2地震震源机制解及发震构造 . 华北地震科学, 2020, 38(S2): 13-18. doi: 10.3969/j.issn.1003-1375.2020.S2.003
    [2] 张素欣, 孙丽娜, 李永庆, 王晓山.  2019年12月5日河北唐山M4.5地震的发震背景 . 华北地震科学, 2020, 38(S1): 76-79. doi: 10.3969/j.issn.1003-1375.2020.S1.017
    [3] 聂冠军, 张忠利, 李海, 李蕾, 杨钦杰, 赵修敏.  2019年广西靖西MS5.2地震发震构造 . 华北地震科学, 2020, 38(S2): 1-6. doi: 10.3969/j.issn.1003-1375.2020.S2.001
    [4] 李蕾, 黄美丽, 阎春恒.  广西靖西地区2次强降雨触发的震群活动 . 华北地震科学, 2020, 38(S2): 61-67. doi: 10.3969/j.issn.1003-1375.2020.S2.010
    [5] 张志慧, 李亚军, 董腾超, 崔鑫.  山东长清M4.1地震序列重定位、发震构造及震级偏差 . 华北地震科学, 2020, 38(S1): 90-94. doi: 10.3969/j.issn.1003-1375.2020.S1.020
    [6] 张斌, 苏道磊, 申金超, 董晓娜, 程海泽.  山东济南4.1级地震序列重定位及其发震构造分析 . 华北地震科学, 2020, 38(2): 85-90. doi: 10.3969/j.issn.1003-1375.2020.02.015
    [7] 曾新福, 汤兰荣, 江春亮, 梁帅, 牛在坡.  九江-瑞昌5.7级地震地质灾害特征及发震构造 . 华北地震科学, 2018, 36(2): 8-17. doi: 10.3969/j.issn.1003-1375.2018.02.002
    [8] 李铂, 崔鑫, 苗庆杰, 穆娟.  CAP方法反演乳山震群3.0级以上地震震源机制解 . 华北地震科学, 2016, 34(3): 14-19. doi: 10.3969/j.issn.1003-1375.2016.03.003
    [9] 李守勇, 孙晴, 晁福林, 张新东.  利用静态库仑应力触发准则确定地震断层面 . 华北地震科学, 2016, 34(2): 67-72. doi: 10.3969/j.issn.1003-1375.2016.02.012
    [10] 郭斌, 张鹏, 王俊, 王斌, 刘琨.  2012年7月20日高邮M4.9级地震发震构造探讨 . 华北地震科学, 2015, 33(4): 24-29,36. doi: 10.3969/j.issn.1003-1375.2015.04.005
    [11] 高熹微, 万永革, 黄骥超, 李祥, 崔华伟.  内蒙古阿拉善左旗MS5.8地震的构造应力场和静态库伦应力变化分析 . 华北地震科学, 2015, 33(2): 48-54. doi: 10.3969/j.issn.1003-1375.2015.02.009
    [12] 张洪艳, 康建红, 卢燕红, 张宇, 陈聪.  2013年吉林前郭5.8级震群地震序列尾波Q值研究 . 华北地震科学, 2014, 32(3): 10-14. doi: 10.3969/j.issn.1003-1375.2014.03.003
    [13] 李冬梅, 周翠英, 董翔, 许萍, 赵金花.  山东地区震群活动与周围地区中强地震的关系 . 华北地震科学, 2011, 29(4): 21-26.
    [14] 吕子强, 肖兰喜, 许萍.  濮阳地震集中区地震精定位及其预测意义 . 华北地震科学, 2008, 26(2): 12-15.
    [15] 解用明, 韩和平, 鞠永, 袁国旭, 孙澎涛.  邢台7.2级地震的发震背景特征 . 华北地震科学, 2006, 24(2): 55-57.
    [16] 谷一山, 毛玉平.  1988年耿马7.2级地震断层特征及其发震构造研究 . 华北地震科学, 1992, 10(1): 1-11.
    [17] 高祥林.  邢台地震的发震构造 . 华北地震科学, 1991, 9(1): 51-58.
    [18] 郭文.  凹凸型发震构造与邢台地震 . 华北地震科学, 1987, 5(s1): 224-230.
    [19] 李绍柄.  唐山地震的发震构造与成因的探讨 . 华北地震科学, 1986, 4(2): 56-62.
    [20] 陶标生.  海原8.5级地震发震构造与应力场的讨论 . 华北地震科学, 1986, 4(1): 48-55.
  • 加载中
图(4) / 表 (3)
计量
  • 文章访问数:  379
  • HTML全文浏览量:  41
  • PDF下载量:  2
出版历程
  • 收稿日期:  2020-07-15
  • 网络出版日期:  2022-07-03
  • 刊出日期:  2022-01-15

2020年北京顺义震群发震机理

doi: 10.3969/j.issn.1003-1375.2022.01.012
    基金项目:  2021年度震情跟踪定向任务(2021010128);天津市地震局局内科研项目(yb202001);地震科技星火计划项目(XH21002Y)
    作者简介:

    宋程(1992—),女,黑龙江牡丹江人,工程师,主要从事地震预报与地震活动性相关工作. E-mail: nibushiwode798@163.com

    通讯作者: 谭毅培(1982—),男,天津人,高级工程师, 主要从事地震预测预报和地震序列分析工作. E-mail: oivertan921@sina.cn
  • 中图分类号: P315.6

摘要: 以2020年3—4月北京顺义小震群为研究对象,进行模板匹配遗漏地震检测、波形互相关震相检测、地震精定位、震源机制计算、破裂半径和应力降等震源参数计算,以探究其发震机理。计算结果显示:① 顺义震群发震构造为NW走向,震中呈现由NW向SE的迁移特征;② 除了震群中第一个ML2.4地震无法计算外,另外2次ML2.0以上地震的破裂面均处于之前地震产生的静态应力加载区内;③ 未检测到震群中存在重复地震活动。结合前人提出的震群发震机理模型,分析认为随机性模型可以较好地解释以上计算结果,顺义震群的发震机理可以解释为地震间的应力触发作用,此研究结果可为北京地区地震危险性分析提供科学参考。

English Abstract

宋程,马婷,徐小远,等. 2020年北京顺义震群发震机理[J]. 华北地震科学,2022, 40(1):71-76. doi:10.3969/j.issn.1003−1375.2022.01.012
引用本文: 宋程,马婷,徐小远,等. 2020年北京顺义震群发震机理[J]. 华北地震科学,2022, 40(1):71-76. doi:10.3969/j.issn.1003−1375.2022.01.012
SONG Cheng,MA Ting,XU Xiaoyuan,et al. Mechanism Analysis of the 2020 Beijing Shunyi Earthquake Swarm[J]. North China Earthquake Sciences,2022, 40(1):71-76. doi:10.3969/j.issn.1003−1375.2022.01.012
Citation: SONG Cheng,MA Ting,XU Xiaoyuan,et al. Mechanism Analysis of the 2020 Beijing Shunyi Earthquake Swarm[J]. North China Earthquake Sciences,2022, 40(1):71-76. doi:10.3969/j.issn.1003−1375.2022.01.012
    • 震群发震机理研究对区域地震活动和断层应力状态分析有重要作用。近十年来,华北地区5级以上地震频次较低,而小震群活动相对频繁。对群发小震物理机理的研究,成为了解华北现今地震活动状态的一个窗口。

      针对震群发震机理,前人根据不同震例研究提出了多种模型。Vidale等[1]基于对美国南加州震群活动的分析,将震群发震机理归纳为3种常见模型,即级联触发模型、断层慢滑动模型和流体侵入模型;Gomberg[2]进一步抽象概括为随机性模型(Stochastic model)和确定性模型(Deterministic model)两类。随机性模型主张地震序列的发生是由应力变化随机触发的,其机制与大地震改变区域应力状态从而触发余震活动类似;而确定性模型认为地震序列的发生和发展除了应力作用之外,还可能受流体渗流或无震滑动等其他物理过程的驱动。分析震群发震机理,可通过震中迁移特征分析[3-4]、重复地震识别[5-6]、序列统计特征分析[7-9]和静态应力触发作用计算[10-11]等方法和手段,判别其符合何种模型。

      近年来,对小震序列发震机理的研究逐渐向精细化和微观化发展。随着地震台站密度的不断增加,以及数字信号处理技术的不断发展,深度挖掘数字化地震波形所携带的震源信息,在更小的尺度上、更加精细地分析地震间的相互作用,逐渐成为发震机理研究的一个热点。Ellsworth等[10]运用了模板匹配地震检测,基于波形互相关的地震精定位、波形聚类和震源谱参数计算等方法,对1999年土耳其Izmit地震前震序列进行深入分析,认为其发震机理符合级联触发模型。Yoon等[11]对1999年美国Hector Mine地震的前震序列,使用模式识别地震检测、波形互相关地震精定位和谱比法震源参数计算,推测其发震机理以应力触发作用为主。此类研究的对象多为大地震的前震序列,而基于数字化地震波形对小震群地震间相互作用的精细化研究尚鲜有报道。

      2020年3—4月,北京市顺义区发生一次小震群活动(图1),震中位于张家口-渤海地震带。在100 km范围内历史上曾发生多次6级以上地震,包括1679年三河-平谷8级地震[12]。然而,自1996年12月16日北京顺义MS4.0地震后至2020年4月,北京地区未再次发生4级以上地震。鉴于小地震的数量要远多于中强地震,对小震活动的研究可在一定程度上弥补中强以上地震缺失的不足[13]。因而分析顺义小震群的发震机理,对了解本地区现今地震活动状态和地震危险性,有一定指示意义。本文以顺义震群为研究对象,基于数字地震波形进行地震精定位和静态应力触发作用计算,分析地震间的相互作用和此次震群的发震机理。

      图  1  顺义震群震中及周边台站分布图

    • 2020年3月29日至4月9日,北京顺义发生小震群活动。北京地震台网共记录到地震8个,其震中位置较为接近。其中ML0.0~0.9地震2个,ML1.0~1.9地震3个,ML2.0以上地震3个,最大地震为4月1日ML2.5地震(表1)。

      表 1  顺义震群ML1.0以上地震信息表

      序号发震时刻震中经度/
      (°E)
      震中纬度/
      (°N)
      深度/
      km
      震级
      ML
      12020-03-29 16:09:33.3116.62340.16632.4
      22020-03-29 17:36:21.5116.62540.16621.4
      32020-04-01 09:46:17.3116.62040.17862.5
      42020-04-09 08:08:17.3116.61640.17442.2
      52020-04-09 10:00:48.9116.60740.18631.0
      62020-04-09 11:24:21.2116.62440.17231.1

      采用模板匹配方法[14]对震群活动期间连续波形进行检测。选取ML2.0以上的3个地震作为模板,截取距离震中较近且记录波形信噪比较高的牛栏山台(NLS)、东三旗台(DSQ)和马道峪台(MDY)三分量直达S波到时前2 s到后2 s的波形,经过零相移4阶Butterworth滤波器2~8 Hz滤波后,在连续波形上进行互相关扫描。当满足3个台站互相关系数平均大于9倍绝对离差中位数或单台互相关系数大于0.7时,作为检测到的地震事件。扫描检测时间段(3月29日至4月9日)内,未检测到目录遗漏的地震事件。

      利用波形互相关震相检测技术,重新标定震群中ML1.0以上地震的直达P波和S波到时。鉴于区域台网提供的震相报告中存在一定人工量取震相到时误差,对精定位结果可靠性造成影响,因而加入波形互相关信息对震相到时进行校正。图2比较了6个ML1.0以上地震在东三旗台(DSQ)的记录波形,分别截取垂直向P波到时前0.5 s至后1.0 s波形,水平向S波到时前0.5 s至后1.5 s波形,经过2~8 Hz带通滤波。鉴于6个地震的波形相似度较高,可以利用波形互相关校正P波和S波到时。选取3月29日ML2.4地震(表1中1号地震)为模板,人工重新校核直达P波和S波到时(图1中参与精定位计算台站),再利用波形互相关标定了其他5个ML1.0以上地震的震相到时。

      图  2  6个ML≥1.0地震在东三旗台(DSQ)的记录波形比较图

      基于波形互相关重新标定过的震相到时数据,使用双差定位方法[15]对震群中ML1.0以上地震进行精定位。鉴于地震波形具有高度相似性,其震源位置接近,因而精定位初始位置设为震群最大地震4月1日ML2.5地震的震源位置。计算中使用奇异值分解(SVD)方法进行反演。参考基于深反射剖面[16]和接收函数[17]等顺义地区地下结构研究结果,精定位使用一维P波速度结构(表2),VP/VS波速比设为1.73。

      表 2  地震精定位所使用速度结构

      123456
      深度范围/km0~22~44~1515~3232~38≥38
      P波速度/(km/s)3.75.56.16.56.88.0

      地震精定位结果如图3所示。其中,地震的标号与表1中序号一致,表示地震发生的时间顺序。双差定位方法给出的平均水平误差约为11 m,平均深度误差约为39 m。鉴于反演使用的是SVD方法,其给出的反演误差是可靠的。对于ML1.0以上地震显示出2个特征:①地震震中呈现较为显著的NW向分布;②地震震中随时间向南东方向迁移。利用P波初动计算了震群中最大事件ML2.5地震的震源机制[18],得到了3组接近的震源机制解,取误差最小的一组为节面I走向132°、倾角56°、滑动角37°,节面II走向19°、倾角60°、滑动角140°。节面I走向与精定位结果走向一致,因而认为节面I可能为震群发震构造。

      图  3  顺义震群ML1.0以上地震精定位与ML2.5地震震源机制结果图

    • 对静态应力触发作用进行分析,首先要计算地震的破裂半径和应力降。由于震群中地震震级较小,波形记录信噪比较低,无法直接通过地震波形计算位移谱来拟合震源参数。故采用赵翠萍等[19]得到的中国大陆中小地震应力降与标量地震矩的拟合关系:

      $$ \mathrm{l}\mathrm{g}\Delta \sigma =0.14\times \mathrm{l}\mathrm{g}{M}_{0}-1.75 $$ (1)

      和区域震级ML与标量地震矩的拟合关系:

      $$ \mathrm{l}\mathrm{g}{M}_{0}=1.02\times {M}_{\mathrm{L}}+9.97 $$ (2)

      基于圆盘破裂断层模型计算标量地震矩计算破裂半径R

      $$ R={\left(\frac{7{M}_{0}}{16\Delta \sigma }\right)}^{1/3} $$ (3)

      再利用地震矩的定义计算平均滑移量D

      $$ D=\frac{{M}_{0}}{\mu \mathrm{{\text π} }{R}^{2}} $$ (4)

      式中:μ为剪切模量,设为30 GPa[11]表3给出震群中6个ML1.0以上地震震源参数计算结果。

      表 3  顺义震群ML1.0以上地震震源参数计算结果

      序号发震时刻震级/ML地震矩
      M0/(N·m)
      破裂半径
      R / m
      平均滑移量
      D / mm
      应力降
      Δσ/MPa
      12020-03-29 16:09:33.32.42.62×1012105.562.490.97
      22020-03-29 17:36:21.51.42.50×101153.840.920.70
      32020-04-01 09:46:17.32.53.31×1012112.912.761.01
      42020-04-09 08:08:17.32.21.64×101292.262.040.91
      52020-04-09 10:00:48.91.09.77×101041.130.610.62
      62020-04-09 11:24:21.21.11.24×101143.990.680.64

      ML1.0以上地震的破裂面投影到断层面上。地震破裂面上滑移量的分布采用Andrews[20]和Yoon等[11]使用的分布公式

      $$ \mathrm{d}\left(r\right)=\left\{\begin{array}{*{20}{l}}D{\left[1-{\left(\dfrac{r}{R}\right)}^{2}\right]}^{3/2}& r{\text{≤}} R\\ 0 & r {\text{>}} R\end{array}\right. {\text{。}} $$ (5)

      使用仲秋等[21]给出的计算静态剪切应力变化的解析式,计算每个地震在断层面上产生的静态应力变化。图4给出了4月1日ML2.5地震与4月9日ML2.2地震发生前断层面上静态应力的变化情况。由于3月29日ML2.4地震为震群的第一个地震,故无法计算其发生前断层面上的静态应力变化。图4a为4月1日ML2.5地震发生前断层面上静态应力的变化计算结果:背景染色表示之前地震(表1中地震1、2)产生的静态应力变化,红色表示应力加载区,即之前发生的地震使得红色区域的剪应力水平升高,反之蓝色表示应力卸载区,即之前发生的地震使得红色区域的剪应力水平降低;黑色圆圈表示4月1日ML2.5地震破裂面,圆点为震源位置。图4b为4月9日ML2.2地震发生前断层面上静态应力的变化计算结果,内容含义与图4a相同。计算结果显示,震群中ML1.0以上地震破裂面范围相互间重合面积较少,意味着震群中未发现重复地震活动;图4中2次ML2.0以上地震破裂面都处于之前地震产生的静态应力加载区内(红色),显示出之前发生的地震可能对这2次地震有应力加载作用。

      图  4  断层面静态应力触发计算结果

    • 根据地震精定位和静态应力触发作用计算结果,基于Gomberg[2]总结的随机性和因果性2类模型,对顺义震群可能的发震机理进行了分析。

      本文中计算的是投影在断层面上的库伦应力,震源机制解的结果对于库伦应力的计算结果具有一定影响。相比于倾角和滑动角,投影断层走向对于库伦应力分布起着较为重要的作用。断层走向是通过精定位结果和震源机制解结果共同约束的,且精定位是有误差范围的。而倾角和滑动角,主要影响静态应力在断层面的投影深度,对于库仑破裂应力变化分布形态上的整体影响较弱[22]

      静态应力触发作用计算结果可以用随机性模型解释。除了震群中第一个地震外,另外2次ML2.0以上地震的破裂面均处于之前地震产生的静态应力加载区内,意味着之前的地震对这2次地震的发生有加载作用,地震发生的物理机制可以用应力级联触发模型解释[10-11]。而震群中地震破裂面相邻而重合的面积很小,与地震应力触发作用在其破裂面边缘最强且随距离快速衰减的现象一致[21],也符合应力级联触发模型的特征。

      计算结果中并未出现支持因果性模型的典型现象。重复地震活动能够解释为稳滑区域包围的较强凹凸体的重复破裂[23],被认为是震群活动期间断层存在慢滑动过程的有力证据[324]。断层慢滑动触发震群活动时,震中随时间线性扩展,其扩展尖峰近似为直线[3]。流体侵入引发的震群活动震中扩展尖峰近似为二次型曲线,可用流体扩散方程解释[ 425]。本研究计算结果中,未发现重复地震活动、震中扩展速率呈直线或二次型曲线等现象,即未观测到支持震群发生过程中存在断层慢滑动或流体渗流作用的典型证据。震中单侧迁移现象可以用多种模型进行解释。在地震精定位计算中,观测到震中由北西向南东的震中单侧迁移现象。Ellsworth等[10]将地震由东向西单侧迁移现象作为支持应力触发模型的依据。但Gomberg[2]提出了不同的看法,认为断层慢滑动引发的地震成核区扩展也可能导致震中单侧迁移。对相同观测现象的解释有时具有非唯一性,可用两种不同物理模型予以解释。顺义震群震中单侧迁移无法作为支持的随机性模型或因果性模型的依据。

      级联触发模型认为震群中较早发生的地震改变了局部应力状态,是触发后续地震活动的决定性因素。而断层慢滑动和流体侵入模型并不否认应力的作用,而是强调在应力触发作用以外,还存在某种物理过程对震群活动起到推动作用。如动态应力变化引发断层慢滑动[3]、断层慢滑动引发库伦应力变化[26]、流体作用触发应力状态改变等[27-28]研究,都发现不同物理机制共同推动地震序列发展。如果地震序列可以用简单的应力触发机制描述,则没有必要引入其他物理过程予以解释[10]

      综合以上分析,本研究对顺义震群的计算结果可用应力触发模型解释,而未发现支持断层慢滑动和流体侵入模型的直接依据,其发震机理可以解释为地震间应力触发作用。

    • 针对2020年北京顺义震群,进行了地震精定位和静态应力触发作用计算,结果显示:① 发震构造为NW走向,震中呈现由北西向南东的迁移特征;② 除了震群中第一个地震无法计算外,另外两次ML2.0以上地震的破裂面均处于之前地震产生的静态应力加载区内;③ 未检测到震群中存在重复地震活动。

      基于随机性模型和确定性模型2类震群发震机理模型,认为随机性模型可以较好地解释以上计算结果,即应力级联触发作用在震群的发生和发展过程中起主要作用。顺义震群的发震机理可以解释为地震间应力触发作用。

参考文献 (28)

目录

    /

    返回文章
    返回