荐言献策

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

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

2020年7月12日河北唐山MS5.1地震发震构造及唐山断裂带三维模型

王亚玲 蔡玲玲 李冬圣 王莉婵 赵英萍 王晓山

王亚玲,蔡玲玲,李冬圣,等. 2020年7月12日河北唐山MS5.1地震发震构造及唐山断裂带三维模型[J]. 华北地震科学,2021, 39(3):83-91. doi:10.3969/j.issn.1003−1375.2021.03.013
引用本文: 王亚玲,蔡玲玲,李冬圣,等. 2020年7月12日河北唐山MS5.1地震发震构造及唐山断裂带三维模型[J]. 华北地震科学,2021, 39(3):83-91. doi:10.3969/j.issn.1003−1375.2021.03.013
WANG Yaling,CAI Lingling,LI Dongsheng,et al. The Seismogenic Structure of the Hebei Tangshan Ms5.1 Earthquake on July 12, 2020 and Three-dimensional Modeling of Tangshan Fault Zone[J]. North China Earthquake Sciences,2021, 39(3):83-91. doi:10.3969/j.issn.1003−1375.2021.03.013
Citation: WANG Yaling,CAI Lingling,LI Dongsheng,et al. The Seismogenic Structure of the Hebei Tangshan Ms5.1 Earthquake on July 12, 2020 and Three-dimensional Modeling of Tangshan Fault Zone[J]. North China Earthquake Sciences,2021, 39(3):83-91. doi:10.3969/j.issn.1003−1375.2021.03.013

2020年7月12日河北唐山MS5.1地震发震构造及唐山断裂带三维模型

doi: 10.3969/j.issn.1003-1375.2021.03.013
基金项目: 河北省地震局地震科技星火重点项目(DZ20180319009)
详细信息
    作者简介:

    王亚玲(1992—),女,河北保定人,助理工程师,主要从事地震监测工作.E-mail:1625715932@qq.com

    通讯作者: 王晓山(1980—),男,河北唐山人,高级工程师,主要从事地震学研究. E-mail:409982349@qq.com
  • 中图分类号: P315.63

The Seismogenic Structure of the Hebei Tangshan Ms5.1 Earthquake on July 12, 2020 and Three-dimensional Modeling of Tangshan Fault Zone

  • 摘要: 基于2020年7月12日河北唐山MS5.1地震事件,以河北省测震台网测定的唐山5.1级地震序列目录中ML1.0以上地震事件为模板,采用匹配滤波技术对2020年07月12—16日的连续波形进行遗漏地震检测,并确定遗漏地震的发震时刻,进一步完善地震序列目录;计算了主震的震源机制解,得到节面I:走向 236°/倾角 77°/滑动角 −170°,节面Ⅱ:走向 143.7°/倾角 80.3°/滑动角 −13.2°;地震事件经震相到时校正后,采用双差定位方法进行精定位,重定位后的余震序列优势分布方向为NE向,结合发震构造背景,震源机制解节面Ⅰ与唐山-古冶断裂东北段的走向、倾角大体一致,由此推测唐山-古冶断裂为河北唐山MS5.1地震的发震断裂,节面I为同震破裂面;最后,基于 SKUA-GOCAD 软件,对发震断裂进行了三维模型的初步模拟和刻画,作为辅助工具为判定该地震的发震构造提供了更为直观的依据支撑。
  • 图  1  唐山地区构造及MS5.1地震序列分布图

    图  2  Eq712064434 模板滦县台的互相关扫描序列

    图  3  滦县台(LUX)震相标定过程示意图

    注:a图为地震事件在LUX台的垂直向记录波形,选取P波到时前0.5 s至到时后1.5 s的波形检测P波到时;b图地震事件在LUX台的东西向记录波形,选取S波到时前0.5 s至到时后1.5 s的波形检测S波到时;黑色短线代表台网震相报告给出的震相到时数据;红线表示基于波形互相关震相检测技术标定的震相到时

    图  4  序列精定位前后震中分布及震源机制解结果图

    图  5  唐山MS5.1地震参与计算台站的波形拟合结果及深度拟合残差分布图

    图  6  发震区域断层及余震三维空间展布

    图  7  研究区主要断裂带的三维几何模型

    表  1  研究区域周边断裂参数表[1924]

    断裂名称长度/km走向倾向倾角/(°)断层性质
    滦县-乐亭断裂(F190NNW(北段)
    NW(南段)
    NE(北段)
    SW(南段)
    35~50高角度逆断层(北段)
    正断层(南段)
    陡河断裂(F220NE25°~30°NW70正断层
    巍山-长山南坡断裂(F321NENW80挤压逆冲为主
    唐山-古冶断裂(F432NENW70~80逆冲性质
    榛子镇断裂(F580NESE60~80正断层
    宁河-昌黎断裂(F6170NESE35~50正断层
    蓟运河断裂(F760NW50ºSW70正断层
    下载: 导出CSV

    表  2  模板地震选取

    模板编号
    发震时刻(北京时间)ML
    纬度/(°N)
    经度/(°E)
    年-月-日时:分:秒
    Eq7120644342020-07-1206:44:34.81.339.758118.434
    Eq7120702252020-07-1207:02:25.82.939.772118.452
    Eq7120726502020-07-1207:26:50.12.639.781118.463
    Eq7121410342020-07-1214:10:34.31.239.762118.443
    Eq7122051032020-07-1220:51:03.92.439.775118.447
    Eq7122146072020-07-1221:46:07.51.039.750118.444
    Eq7130049282020-07-1300:49:28.41.239.747118.435
    Eq7130538002020-07-1305:38:00.01.339.752118.440
    Eq7131212482020-07-1312:12:48.71.139.761118.440
    Eq7141231082020-07-1412:31:08.21.439.757118.454
    Eq7152145272020-07-1521:45:27.71.839.671118.331
    Eq7160959252020-07-1609:59:25.11.339.749118.444
    下载: 导出CSV

    表  3  遗漏地震发震时刻和震级

    发震时刻震级ML
    纬度/(°N)
    经度/(°E)
    年-月-日时:分:秒
    2020-07-1206:45:01.10.539.72118.41
    2020-07-1206:45:29.70.439.77118.43
    2020-07-1206:46:29.80.539.76118.45
    2020-07-1207:00:18.90.439.82118.44
    2020-07-1208:34:47.60.039.69118.50
    2020-07-1210:16:47.60.339.79118.49
    2020-07-1306:38:34.40.539.68118.38
    2020-07-1521:52:51.51.439.66118.32
    下载: 导出CSV

    表  4  唐山地区速度模型

    深度/km05.010.015.020.025.035.0
    VP/(km.s−15.25.86.06.36.56.77.0
    下载: 导出CSV

    表  5  MS5.1主震震源机制解结果

    节面Ι/(°)节面Ⅱ/(°)P轴/(°)T轴/(°)NB)轴/(°)
    走向倾角滑动角走向倾角滑动角方位角俯仰角方位角俯仰角方位角俯仰角
    23677−170143.780.3−13.299.516.2190.22.3287.973.7
    下载: 导出CSV
  • [1] 邵霄怡, 王晓青, 窦爱霞, 等. 基于蒙特卡洛方法的地震目录模拟及相符性检验——以汾渭地震带为例[J]. 中国地震, 2018, 34(1): 35-47. doi:  10.3969/j.issn.1001-4683.2018.01.004
    [2] Woessner J, Wiemer S. Assessing the quality of earthquake catalogues: estimating the magnitude of completeness and its uncertainty[J]. Bulletin of the Seismological Society of America, 2005, 95(2): 684-698. doi:  10.1785/0120040007
    [3] 冯建刚, 蒋长胜, 韩立波, 等. 甘肃测震台网监测能力及地震目录完整性分析[J]. 地震学报, 2012, 34(5): 646-658. doi:  10.3969/j.issn.0253-3782.2012.05.006
    [4] Peng Z G, Vidale J E, Ishii M, et al. Seismicity rate immediately before and after main shock rupture from high-frequency waveforms in Japan[J]. Journal of Geophysical Research: Solid Earth, 2007, 112(B3): B03306.
    [5] 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
    [6] Schaff D. Improvements to detection capability by cross-correlating for similar events: a case study of the 1999 Xiuyan, China, sequence and synthetic sensitivity tests[J]. Geophysical Journal International, 2010, 180(2): 829-846. doi:  10.1111/j.1365-246X.2009.04446.x
    [7] Meng X F, Yu X, Peng Z G, et al. Detecting earthquakes around Salton Sea following the 2010 Mw7.2 El Mayor-Cucapah earthquake using GPU parallel computing[J]. Procedia Computer Science, 2012, 9: 937-946. doi:  10.1016/j.procs.2012.04.100
    [8] 谭毅培, 曹井泉, 卞真付, 等. 2013年8月河北蔚县小震群遗漏地震检测与发震构造分析[J]. 地震学报, 2014, 36(6): 1022-1031.
    [9] 谭毅培, 曹井泉, 刘文兵, 等. 2013年3月涿鹿微震群遗漏地震事件检测和发震构造分析[J]. 地球物理学报, 2014, 57(6): 1847-1856. doi:  10.6038/cjg20140616
    [10] 单连君, 张新东, 王晓山. 2016年1月9日河北怀来ML3.4震群发震构造分析[J]. 震灾防御技术, 2017, 12(3): 557-564. doi:  10.11899/zzfy20170312
    [11] 王宁, 王亚玲, 蔡玲玲, 等. 2017年9月临城震群遗漏地震检测及发震构造初探[J]. 地震地磁观测与研究, 2019, 40(6): 9-14. doi:  10.3969/j.issn.1003-3246.2019.06.002
    [12] 张素欣, 边庆凯, 张子广, 等. 唐山断裂北段地震分布特征及其构造意义[J]. 地震研究, 2020, 43(2): 270-277. doi:  10.3969/j.issn.1000-0666.2020.02.008
    [13] 韦士忠, 石汝斌, 李玉萍, 等. 用选频谱研究唐山余震序列较大地震前小震波谱和震源参数的某些特征[J]. 华北地震科学, 1989, 7(1): 43-51.
    [14] 吴培稚, 刘文龙, 陈宇卫. 用不等间距统计法对唐山余震序列的研究[J]. 地震, 1996, 16(3): 255-262.
    [15] 万永革, 沈正康, 曾跃华, 等. 唐山地震序列应力触发的粘弹性力学模型研究[J]. 地震学报, 2008, 30(6): 581-593. doi:  10.3321/j.issn:0253-3782.2008.06.004
    [16] 仲秋, 史保平. 1976年唐山MS7.8地震余震序列持续时间及对地震危险性分析的意义[J]. 地震学报, 2012, 34(4): 494-508. doi:  10.3969/j.issn.0253-3782.2012.04.007
    [17] 杨雅琼, 王晓山, 万永革, 等. 由震源机制解推断唐山地震序列发震断层的分段特征[J]. 地震学报, 2016, 38(4): 632-643. doi:  10.11939/jass.2016.04.009
    [18] 虢顺民, 李志义, 程绍平, 等. 唐山地震区域构造背景和发震模式的讨论[J]. 地质科学, 1977, 12(4): 305-321.
    [19] 刘芳晓. 利用横波地震反射技术实现滦县-乐亭断裂超浅层识别的方法研究[J]. 西部探矿工程, 2019, 31(5): 130, 133.
    [20] 王景明, 王文秋, 苏幼波. 唐山市区断裂活动与地质灾害[J]. 河北理工学院学报, 1998, 20(1): 75-83.
    [21] 王景明, 陈国顺, 郑文俊, 等. 唐山7.8级和7.1级地震地裂缝及地震成因探讨[J]. 长安大学学报: 地球科学版, 1981(2): 58-69.
    [22] 彭远黔, 孟立朋. 河北地震构造特征[M]. 石家庄: 河北人民出版社, 2017: 54-55.
    [23] 李赫, 董一兵, 王熠熙, 等. 唐山地区地震成因初探——以丰南2010年M4.1和2019年M4.5地震序列为例[J]. 地球物理学进展, 2020, 35(4): 1276-1283. doi:  10.6038/pg2020EE0139
    [24] 陈宏强, 专少鹏, 陈超, 等. 河北省唐山地区丰台-野鸡坨断裂第四纪活动性——来自14C和磁性地层年代学的证据[J]. 中国地质, 2021, 48(2): 605-617.
    [25] 易桂喜, 龙锋, 赵敏, 等. 2014年10月1日越西M5.0地震震源机制与发震构造分析[J]. 地震地质, 2016, 38(4): 1124-1136. doi:  10.3969/j.issn.0253-4967.2016.04.025
    [26] 易桂喜, 龙锋, 梁明剑, 等. 2016年9月23日四川理塘M4.9和M5.1地震震源机制与发震构造分析[J]. 国际地震动态, 2017(8): 50-51. doi:  10.3969/j.issn.0253-4975.2017.08.033
    [27] 张炜, 申春生, 邢玉安, 等. 地震短临异常新指标的探索——汞浓度探测[J]. 中国地震, 1989, 5(4): 13-19.
    [28] 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
    [29] 郭蕾, 宫猛, 王宁. 2012年唐山ML5.2地震前后的速度结构对比分析[J]. 中国地震, 2019, 35(4): 695-708. doi:  10.3969/j.issn.1001-4683.2019.04.010
    [30] 鲁人齐, 徐锡伟, 陈立春, 等. 2017年8月8日九寨沟MS7.0地震构造与震区三维断层初始模型[J]. 地震地质, 2018, 40(1): 1-11. doi:  10.3969/j.issn.0253-4967.2018.01.001
    [31] 王金艳, 鲁人齐, 张浩, 等. 郯庐断裂带江苏段新生界三维地质构造建模[J]. 地震学报, 2020, 42(2): 216-230. doi:  10.11939/jass.20190131
  • [1] 聂冠军, 张忠利, 李海, 李蕾, 杨钦杰, 赵修敏.  2019年广西靖西MS5.2地震发震构造 . 华北地震科学, 2020, 38(S2): 1-6. doi: 10.3969/j.issn.1003-1375.2020.S2.001
    [2] 肖本夫, 祁玉萍, 龙锋, 申源, 梁厚朗, 田优平.  2014年以来云南永善地区M≥4.7地震震源机制解及发震构造 . 华北地震科学, 2020, 38(3): 13-18. doi: 10.3969/j.issn.1003-1375.2020.03.003
    [3] 张志慧, 李亚军, 董腾超, 崔鑫.  山东长清M4.1地震序列重定位、发震构造及震级偏差 . 华北地震科学, 2020, 38(S1): 90-94. doi: 10.3969/j.issn.1003-1375.2020.S1.020
    [4] 张斌, 苏道磊, 申金超, 董晓娜, 程海泽.  山东济南4.1级地震序列重定位及其发震构造分析 . 华北地震科学, 2020, 38(2): 85-90. doi: 10.3969/j.issn.1003-1375.2020.02.015
    [5] 马婷, 谭毅培, 丁晶.  大震面波应力扰动对2014年四川康定地震的影响 . 华北地震科学, 2019, 37(4): 29-33. doi: 10.3969/j.issn.1003-1375.2019.04.006
    [6] 曾新福, 汤兰荣, 江春亮, 梁帅, 牛在坡.  九江-瑞昌5.7级地震地质灾害特征及发震构造 . 华北地震科学, 2018, 36(2): 8-17. doi: 10.3969/j.issn.1003-1375.2018.02.002
    [7] 王宁, 杨春利, 李雪英, 常亮, 赵英萍.  遗漏震相检测技术在2016年9月唐山震群中的应用 . 华北地震科学, 2017, 35(1): 38-43. doi: 10.3969/j.issn.1003-1375.2017.01.006
    [8] 郭斌, 张鹏, 王俊, 王斌, 刘琨.  2012年7月20日高邮M4.9级地震发震构造探讨 . 华北地震科学, 2015, 33(4): 24-29,36. doi: 10.3969/j.issn.1003-1375.2015.04.005
    [9] 高祥林.  邢台地震的发震构造 . 华北地震科学, 1991, 9(1): 51-58.
    [10] 李绍柄.  关于唐山地震的发震构造与成因的探讨一文的补充 . 华北地震科学, 1991, 9(3): 97-98.
    [11] 陈达生, 刘汉兴.  地震烈度椭圆衰减关系 . 华北地震科学, 1989, 7(3): 31-42.
    [12] 黄兴根, 张英礼, 王瑛.  唐山古地震遗迹的发现及其意义 . 华北地震科学, 1987, 5(s1): 347-354.
    [13] 业成之.  从跨断层测量资料探讨唐山地震孕育的断裂运动特征 . 华北地震科学, 1987, 5(s1): 292-299.
    [14] 卞兆银, 王进英.  华北地区地震活动特点及地震构造带的划分 . 华北地震科学, 1987, 5(s1): 239-248.
    [15] 郭文.  凹凸型发震构造与邢台地震 . 华北地震科学, 1987, 5(s1): 224-230.
    [16] 陶标生.  海原8.5级地震发震构造与应力场的讨论 . 华北地震科学, 1986, 4(1): 48-55.
    [17] 李绍柄.  唐山地震的发震构造与成因的探讨 . 华北地震科学, 1986, 4(2): 56-62.
    [18] 李介成, 陶建华.  四川西部温泉水温变化与地震前兆的研究 . 华北地震科学, 1985, 3(4): 77-85.
    [19] 于新昌, 靳雅敏, 孙迎.  晋冀鲁豫交界地区的地震活动性估计 . 华北地震科学, 1985, 3(4): 52-60.
    [20] 张四昌.  华北地区的地震分布图象与共轭孕震构造 . 华北地震科学, 1985, 3(2): 9-18.
  • 加载中
图(7) / 表 (5)
计量
  • 文章访问数:  1256
  • HTML全文浏览量:  1081
  • PDF下载量:  13
出版历程
  • 收稿日期:  2020-01-20
  • 网络出版日期:  2021-08-10
  • 刊出日期:  2021-07-20

2020年7月12日河北唐山MS5.1地震发震构造及唐山断裂带三维模型

doi: 10.3969/j.issn.1003-1375.2021.03.013
    基金项目:  河北省地震局地震科技星火重点项目(DZ20180319009)
    作者简介:

    王亚玲(1992—),女,河北保定人,助理工程师,主要从事地震监测工作.E-mail:1625715932@qq.com

    通讯作者: 王晓山(1980—),男,河北唐山人,高级工程师,主要从事地震学研究. E-mail:409982349@qq.com
  • 中图分类号: P315.63

摘要: 基于2020年7月12日河北唐山MS5.1地震事件,以河北省测震台网测定的唐山5.1级地震序列目录中ML1.0以上地震事件为模板,采用匹配滤波技术对2020年07月12—16日的连续波形进行遗漏地震检测,并确定遗漏地震的发震时刻,进一步完善地震序列目录;计算了主震的震源机制解,得到节面I:走向 236°/倾角 77°/滑动角 −170°,节面Ⅱ:走向 143.7°/倾角 80.3°/滑动角 −13.2°;地震事件经震相到时校正后,采用双差定位方法进行精定位,重定位后的余震序列优势分布方向为NE向,结合发震构造背景,震源机制解节面Ⅰ与唐山-古冶断裂东北段的走向、倾角大体一致,由此推测唐山-古冶断裂为河北唐山MS5.1地震的发震断裂,节面I为同震破裂面;最后,基于 SKUA-GOCAD 软件,对发震断裂进行了三维模型的初步模拟和刻画,作为辅助工具为判定该地震的发震构造提供了更为直观的依据支撑。

English Abstract

王亚玲,蔡玲玲,李冬圣,等. 2020年7月12日河北唐山MS5.1地震发震构造及唐山断裂带三维模型[J]. 华北地震科学,2021, 39(3):83-91. doi:10.3969/j.issn.1003−1375.2021.03.013
引用本文: 王亚玲,蔡玲玲,李冬圣,等. 2020年7月12日河北唐山MS5.1地震发震构造及唐山断裂带三维模型[J]. 华北地震科学,2021, 39(3):83-91. doi:10.3969/j.issn.1003−1375.2021.03.013
WANG Yaling,CAI Lingling,LI Dongsheng,et al. The Seismogenic Structure of the Hebei Tangshan Ms5.1 Earthquake on July 12, 2020 and Three-dimensional Modeling of Tangshan Fault Zone[J]. North China Earthquake Sciences,2021, 39(3):83-91. doi:10.3969/j.issn.1003−1375.2021.03.013
Citation: WANG Yaling,CAI Lingling,LI Dongsheng,et al. The Seismogenic Structure of the Hebei Tangshan Ms5.1 Earthquake on July 12, 2020 and Three-dimensional Modeling of Tangshan Fault Zone[J]. North China Earthquake Sciences,2021, 39(3):83-91. doi:10.3969/j.issn.1003−1375.2021.03.013
    • 强震发生后或者震群活动时,在震中附近或小区域范围内往往会有大量余震发生,其中震级较小的地震震相往往会淹没在背景噪声中难以识别,成为遗漏地震。地震目录是区域地震活动性研究、地震危险性分析的重要数据资料[1],其完整准确性直接影响区域地震活动、地震危险性的分析,并对后续科学研究的可信性和科学性产生影响[2-3],尤其对发震构造及震后趋势分析等研究带来直接影响。已有学者应用匹配滤波技术对中强地震余震序列完整性和余震分布特征[4-7]等方面进行研究,谭毅培等分别对2013 年3 月涿鹿微震群、2013 年8 月蔚县震群进行遗漏地震检测分析其发震构造[8-9];单连君等基于匹配滤波技术分析河北怀来ML3.4震群余震序列[10];王宁等将遗漏震相检测技术应用在2016 年9 月唐山震群中[11]。本文使用匹配滤波技术识别唐山MS5.1地震的遗漏地震进行其发震构造的确定。

      唐山断裂带处于近EW向阴山-燕山南缘活动构造带和NNE向冀鲁断块拗陷相汇的地方[12]。1976 年唐山7.8 级地震后,沿唐山断裂带余震十分活跃,并有大量关于该地震的发震构造、余震序列研究成果产出。韦士忠等通过计算选频谱,研究了唐山地震余震序列中1979 年5 月17 日5.1级地震事件,发现震源等效半径较小,应力降较高[13],地震序列的频谱规律为认识孕育较大地震的过程提供信息线索;吴培稚等采用不相等间距的统计方法,研究唐山的余震序列展布[14],从中得到一些预报的指标;万永革等通过唐山地震序列的应力变化研究,表明唐山地震序列中的前面大震会对后续小震的发生起到了调制作用[15];仲秋等通过研究1976年唐山7.8级大地震的余震序列,得到相应的余震序列持续时间公式[16],研究学者对唐山地震余震序列的研究领域及方法是多方面的,并且得到相应的数据指标,为认识主震与余震之间的相互关系,及对地震危险性估计有一定的意义;另外杨雅琼等基于255 个震源机制解对唐山地震序列的发震断层进行分段[17]。而余震展布特征不仅与断层走向有关,也可能与断层深部的几何展布有关,因此研究余震几何分布和构造特征具有重要的科学意义。

      本文基于2020 年7 月12 日6 时38 分唐山MS5.1地震,对7 月12—16 日的连续波形通过匹配滤波技术进行遗漏地震检测,估计遗漏地震的震级,完善地震序列目录;并对相似波形地震进行震相校正提高定位精度,结合震源机制解推测发震构造;再利用 SKUA-GOCAD 软件给出发震构造的初步三维模型,通过二维分析与三维模型相结合达到分析此次地震构造活动特征的目的,为进一步研究该地震及邻区地震活动性提供基础数据支持。

    • 唐山地区构造复杂如图1,受宁河-昌黎断裂(F6)、榛子镇断裂(F5)、滦县-乐亭断裂(F1)和蓟运河断裂(F7)的切割,唐山地区形成了一个 NE 向的菱形地块[18]。其中央是一条NE向的由陡河断裂(F2)、巍山-长山南坡断裂(F3)和唐山-古冶断裂(F4)组成的唐山断裂带,总体走向NE-NNE,倾向NW、SE,倾角70°~80°,为逆走滑断层,长约59 km。断裂详情见表1

      图  1  唐山地区构造及MS5.1地震序列分布图

      表 1  研究区域周边断裂参数表[1924]

      断裂名称长度/km走向倾向倾角/(°)断层性质
      滦县-乐亭断裂(F190NNW(北段)
      NW(南段)
      NE(北段)
      SW(南段)
      35~50高角度逆断层(北段)
      正断层(南段)
      陡河断裂(F220NE25°~30°NW70正断层
      巍山-长山南坡断裂(F321NENW80挤压逆冲为主
      唐山-古冶断裂(F432NENW70~80逆冲性质
      榛子镇断裂(F580NESE60~80正断层
      宁河-昌黎断裂(F6170NESE35~50正断层
      蓟运河断裂(F760NW50ºSW70正断层
    • 使用模版匹配滤波技术(matched filter technique,MFT)进行遗漏地震检测,可实现在低信噪比情况下提取地震信息。依据震相报告到时数据,挑选余震序列中震级较大的、信噪比好的地震事件,以S波到时为中心,从连续波形数据中截取S波到时前2 s至后2 s波形,筛选波形三分量信噪比的平均值大于3的地震事件作为模板。本文选取7 月12—16 日ML1.0以上震中距最近3 个台站震相清楚的余震(表2)为模板,台站选择遵循标准:①宽频带地震计的台站,该类型台站记录波形质量较好;②震中距较近的台站,信噪比较高。因此,选取了DOH、LUX、QIX 3 个台站参与互相关计算。

      表 2  模板地震选取

      模板编号
      发震时刻(北京时间)ML
      纬度/(°N)
      经度/(°E)
      年-月-日时:分:秒
      Eq7120644342020-07-1206:44:34.81.339.758118.434
      Eq7120702252020-07-1207:02:25.82.939.772118.452
      Eq7120726502020-07-1207:26:50.12.639.781118.463
      Eq7121410342020-07-1214:10:34.31.239.762118.443
      Eq7122051032020-07-1220:51:03.92.439.775118.447
      Eq7122146072020-07-1221:46:07.51.039.750118.444
      Eq7130049282020-07-1300:49:28.41.239.747118.435
      Eq7130538002020-07-1305:38:00.01.339.752118.440
      Eq7131212482020-07-1312:12:48.71.139.761118.440
      Eq7141231082020-07-1412:31:08.21.439.757118.454
      Eq7152145272020-07-1521:45:27.71.839.671118.331
      Eq7160959252020-07-1609:59:25.11.339.749118.444

      每个地震模板均对7 月12—16 日的连续波形进行扫描,计算互相关系数,将每个模板三分向数据的互相关结果进行组合,通过计算序列的绝对离差中位数(MAD)检测遗漏地震,其表达式为:

      $$ MAD = {\rm median}(\left| {{X_i}} \right| - \bar X)$$

      式中:Xi代表第i个互相关系数序列;$\overline X $为其平均值。本次运算中MAD的倍数阈值是9[5]。为了减少地震序列因仅有一两个台站距离震中较近而造成检测结果失误,还需设置单台互相关系数阈值,默认值为0.75[10]

      以Eq712064434地震模板的滦县台(LUX)扫描连续波形为例(图2),图中为扫描7 月12 日连续波形结果,扫描模板所在的连续波形互相关系数为1。连续波形通过阈值检测得到的结果为遗漏地震文件,对照原地震目录,序列未记录的事件作为疑似地震事件,最后拾取遗漏地震事件S波形的最大振幅,并利用其与模板地震的振幅比来估计遗漏事件震级。

      图  2  Eq712064434 模板滦县台的互相关扫描序列

    • 区域测震台网震相报告给出的震相到时存在一定的人工误差[10],这种误差对于地震序列时空演化的精细分析有较大的影响,因此本文利用少数几个震级较大的地震事件,通过波形互相关的方法对波形相似的地震事件进行震相校正,有效地减少误差,提高定位精度。

      以Eq712070225模板的滦县台(LUX)为例(图3),对相似波形地震事件进行震相校正。本程序得到的到时信息与台网人工读取震相无关,在一定程度上消除了地震事件在人工定位时产生的相对误差,从而提高了定位结果的精度。

      图  3  滦县台(LUX)震相标定过程示意图

    • 基于每个地震模板所识别出的地震事件的发震时刻,综合所有模板检测出来的地震事件,对比目录结果便可以找出遗漏地震。由于河北测震台网工作方式的改革提升,现采取截取实时波形的方式查看地震余震序列,一般肉眼可识别的地震已实时快报,并编写目录。针对此次地震来说余震序列并不发育,因此利用遗漏检测技术扫描7 月12—16 日连续波形共检测到的地震事件共8 个,遗漏地震震级、发震时刻、经纬度如表3所示。其中:ML0.0~0.5地震共有7 个,ML1.4的1 个。ML0.0~0.5地震编目条目由36增至42 条,表明ML0.0~0.5之间的地震目录的完善性得到明显的提升。虽然大部分检测出的遗漏地震较小,但是对分析此次MS5.1地震震中展布的几何特征具有一定意义,也为构建三维构造模型提供更为完善的数据支持。

      表 3  遗漏地震发震时刻和震级

      发震时刻震级ML
      纬度/(°N)
      经度/(°E)
      年-月-日时:分:秒
      2020-07-1206:45:01.10.539.72118.41
      2020-07-1206:45:29.70.439.77118.43
      2020-07-1206:46:29.80.539.76118.45
      2020-07-1207:00:18.90.439.82118.44
      2020-07-1208:34:47.60.039.69118.50
      2020-07-1210:16:47.60.339.79118.49
      2020-07-1306:38:34.40.539.68118.38
      2020-07-1521:52:51.51.439.66118.32
    • 精确定位后的主震及其余震序列空间展布和震源机制解,对于确定发震构造、判定未来地震趋势具有重要意义[25-26]。利用CAP[27]方法反演主震的震源机制解,该方法是将近震观测波形分解成Pnl和面波,分别进行拟合,对3 个断层面参数的各自取值区间进行格点的取值,将得到的理论地震图与实际观测波形进行拟合,波形拟合最好的一组参数即为求得的震源断层面解。反演本次唐山MS5.1地震的震源机制解采用的速度模型表4,计算结果如图4b。共选取10 个台站参与,其中体波和面波各部分波形拟合互相关系数大于80%的占70%(图5a),显示此次计算结果拟合较好,结果可靠,深度拟合残差显示(图5b),拟合深度为14 km时残差最小,此深度即为最佳震源深度。

      表 4  唐山地区速度模型

      深度/km05.010.015.020.025.035.0
      VP/(km.s−15.25.86.06.36.56.77.0

      图  4  序列精定位前后震中分布及震源机制解结果图

      图  5  唐山MS5.1地震参与计算台站的波形拟合结果及深度拟合残差分布图

      本文对余震序列采用双差定位方法[28]进行精定位,速度模型采用的是郭蕾等[29]使用的一维P波速度模型(表4)。为了便于对比精定位前后的震中分布趋势,删除陡河台、滦县台单台定位目录,图4a给出了精定位前地震目录(包括遗漏地震)的91(含主震)个地震事件的震中分布图。对经过波形互相关震相校正的地震事件进行重新定位,精定位结果包含84 个地震事件。对未能标定3 个以上台站震相到时的地震事件,其震中的位置可以定在与其互相系数最大的模板地震震中处,最终震中分布显示如图4b。由图中精定位结果可知,唐山余震序列分布趋势更为明显,主要集中在唐山-古冶断裂的东北部;震中NE向的优势展布也更为明显,与主震震源机制解的节面Ⅰ走向基本一致(表5)。推测节面Ⅰ为本次唐山MS5.1地震的发震断层面,该地震是发生在走向NE、倾向NW、高倾角的右旋兼走滑断层上,与唐山-古冶断裂(F4)的走向、倾向、倾角基本吻合,由此推测唐山-古冶断裂可能为其发震断裂。

      表 5  MS5.1主震震源机制解结果

      节面Ι/(°)节面Ⅱ/(°)P轴/(°)T轴/(°)NB)轴/(°)
      走向倾角滑动角走向倾角滑动角方位角俯仰角方位角俯仰角方位角俯仰角
      23677−170143.780.3−13.299.516.2190.22.3287.973.7
    • 由上文可知,唐山MS5.1地震发生后截止至2020 年7 月16 日,河北区域测震台网记录可定位余震82 条,经过遗漏检测之后增至90 条。结合图4b所呈现的断层地表几何特征与余震展布关系,对余震空间分布的最大优势分布进行分析和预处理[30-31],基于GOCAD软件初步建立断层在三维空间的地表展布模型,实现地震三维空间展布(图6)。

      图  6  发震区域断层及余震三维空间展布

      图6 清晰展示了余震序列在二维地表展布、三维空间分布特征,更加直观地显示余震空间上的优势展布方向。

    • 建立震区的三维构造模型,需要的数据资料包括以下几个方面:断层地表展露的几何特征、断层走向、倾角、倾向以及地下延展情况的详细描述、相关地震震中展布及震源机制解约束条件等。本文根据已有的断裂资料(图1),初步建立研究震区滦县-乐亭断裂(F1)、陡河断裂(F2)、巍山-长山南坡断裂(F3)、唐山-古冶断裂(F4)以及榛子镇断裂(F5)、宁河-昌黎断裂(F6)6条断裂的三维模型(图6)。

      基于已知的断裂数据资料,在图6的基础上通过GOCAD软件的“make suface”功能造面成图(图7)。

      图  7  研究区主要断裂带的三维几何模型

      图7所示:唐山-古冶断裂(F4)、陡河断裂(F2)、巍山-长山南坡断裂(F3)的空间走向较为准确地表达在三维空间上,实现对发震断裂的三维模拟与刻画。结合地震震源机制解结果及余震序列分布推测节面Ⅰ为同震断裂面,唐山-古冶断裂(F4)与节面Ⅰ的走向、倾角基本一致,结合地震的震中位置推测此次地震可能发生在唐山-古冶断裂(F4)上。在此基础上结合图7模拟的三维空间展示,可以看出余震空间展布趋势与唐山-古冶断裂(F4)NE至NEE走向延展具有一致性。据此,可以更直观、更有力地推测唐山-古冶断裂(F4)即为河北唐山MS5.1地震的发震断裂。

    • 本文通过匹配滤波技术,对2020 年7 月12—16 日唐山MS5.1 地震震中附近台站记录的连续波形进行遗漏地震检测、震相校正来完善地震序列目录,依据序列精定位、主震震源机制解分析发震构造,并基于GOCAD软件初步给出了发震区断层的三维构造模型。

      1)运用匹配滤波技术对余震进行遗漏检测,检测结果发现了遗漏的地震事件,因此测震台网人员可以通过此方法拾取震相,提高地震目录的完整性。

      2)通过波形互相关方法对震相校正,不仅可以估计遗漏地震的震级,还可以利用少数几个震级较大的地震事件对于波形相似的地震进行震相校正,有效地减少人工拾取震相误差,大大提高定位精度。

      3)精定位后的余震序列空间分布趋势更为明显,地震序列主要分布在唐山断裂带中NE向唐山-古冶断裂附近,序列优势分布方向为NE向,分析节面Ⅰ为同震断层面,与唐山-古冶断裂东北段走向、倾向、倾角大体一致,初步推测MS5.1地震发震构造为唐山-古冶断裂。

      4)基于GOCAD软件,结合已知的地质背景、震源机制解等资料,初步建立发震区的三维构造模型;在结论3的基础上,进行三维空间模拟更为直观展示了余震序列的空间分布特征;结合唐山-古冶断裂空间延展分布趋势,展现出一致性。作为辅助工具为此次地震的发震构造判定提供更为直观的论据支撑,为验证发震构造推断提供一种可视化的有效途径。

      基于GOCAD建立的三维构造模型为判断发震机理、震源断层的空间描述提供了一条可借鉴的有效途径。此外,由于唐山MS5.1地震余震序列不发育,样本数据量较小,无法有效、合理地建立三维断层面模型,这也表明初始的三维构造模型具有一定的局限性。未来成熟的三维模型,需要大量数据样本,不断提高各类数据质量,逐渐增加各种成熟的地震活动解译资料及探测资料,发展构造定量解析方法等多元数据条件约束的基础上进行丰富完善。

参考文献 (31)

目录

    /

    返回文章
    返回