荐言献策

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

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

邢台老震区库仑应力演化及地震危险性分析

陈婷 章阳 王晓山 孙丽娜

陈婷,章阳,王晓山,等. 邢台老震区库仑应力演化及地震危险性分析[J]. 华北地震科学,2021, 39(4):50-56. doi:10.3969/j.issn.1003−1375.2021.04.008
引用本文: 陈婷,章阳,王晓山,等. 邢台老震区库仑应力演化及地震危险性分析[J]. 华北地震科学,2021, 39(4):50-56. doi:10.3969/j.issn.1003−1375.2021.04.008
CHEN Ting,ZHANG Yang,WANG Xiaoshan,et al. Coulomb Stress Evolution and Seismic Hazard Analysis on the Xingtai Old Earthquake Region[J]. North China Earthquake Sciences,2021, 39(4):50-56. doi:10.3969/j.issn.1003−1375.2021.04.008
Citation: CHEN Ting,ZHANG Yang,WANG Xiaoshan,et al. Coulomb Stress Evolution and Seismic Hazard Analysis on the Xingtai Old Earthquake Region[J]. North China Earthquake Sciences,2021, 39(4):50-56. doi:10.3969/j.issn.1003−1375.2021.04.008

邢台老震区库仑应力演化及地震危险性分析

doi: 10.3969/j.issn.1003-1375.2021.04.008
基金项目: 河北省地震局强震发震构造与机理研究创新团队(DZ20180319009);2019年震情跟踪定向工作任务(2019010122);河北省地震科技星火计划青年项目(DZ20190419028)
详细信息
    作者简介:

    陈婷(1990—),女,山东聊城人,工程师,主要从事数字地震学及中长期预测等研究. E-mail:chenting_ouc@163.com

  • 中图分类号: P315.7

Coulomb Stress Evolution and Seismic Hazard Analysis on the Xingtai Old Earthquake Region

  • 摘要: 基于Burgers粘弹性层状模型,综合考虑同震位错和震后粘滞松弛效应,利用PSGRN/PSCMP程序计算地震破裂面上产生的库仑破裂应力变化,分析邢台地震序列中每次6级以上地震前的所有6级以上地震对本次地震的综合库仑破裂应力变化,探讨强震间的应力触发关系。结果表明:除3月26日6.2级地震外,其余地震均发生在库仑应力增加的应力触发区,与之前的6级以上地震存在很好的应力触发关系;根据邢台地震序列的库仑应力演化与5级以上地震活动关系讨论邢台老震区的地震危险性,认为5级以上地震发生在应力触发区概率达到71.9%,邢台老震区中强地震活动与库仑应力变化有很好的相关性,可将库仑应力动态演化特征应用到本区地震的中长期预测中,为后续中强震发生地点提供参考。
  • 图  1  邢台老震区地震地质背景

    图  2  邢台地震序列强震间的库仑应力触发作用

    注:a MS6.8地震对MS6.7地震的触发情况;b MS6.8地震和MS6.7地震对MS7.2地震的触发情况;c MS6.8地震、MS6.7地震和MS7.2地震对MS6.2地震的触发情况;d MS6.8地震、MS6.7地震、MS7.2地震及MS6.2地震对MS6.0地震的触发情况;“☆”表示触发地震的震中位置;“○”表示被触发地震的震中位置;1 bar=0.1 MPa

    图  3  各个时间点邢台地震序列库仑应力变化

    注:a 3月29日MS6.0地震同震响应库仑应力变化;b 震后100年库仑应力变化;c 震后300年库仑应力变化;d 震后500年库仑应力变化;1 bar=0.1 MPa。

    图  4  各个时间段库仑应力变化与5级地震的触发关系

    注:a 3月8日MS6.8地震库仑应力变化对3月8—22日之间5级地震的触发关系;b 3月8—22日3个MS6.0以上地震综合库仑应力变化对3月22—26日之间5级地震的触发关系;c 3月8—26日4个MS6.0以上地震综合库仑应力变化对3月26—29日之间5级地震的触发关系;d 3月8—29日5个MS6.0以上地震综合库仑应力变化对3月29日至今5级地震的触发关系;图中“○”表示5到6级地震的震中位置;1 bar=0.1 MPa。

    表  1  地壳介质模型

    层名厚度/kmVP/(km∙s−1VS/(km∙s−1ρ/(kg∙m−3η1/(Pa∙s)η2/(Pa∙s)
    上地壳3.53.201.552 600完全弹性
    4.56.003.452 670
    5.56.203.452 670
    中地壳8.06.653.852 800
    下地壳5.56.803.953 1001.0×10193.0×1020
    5.06.253.253 1001.0×10193.0×1020
    5.07.304.203 1001.0×10193.0×1020
    上地幔8.104.653 36002.0×1021
     注:VP为P波速度;VS为S波速度;ρ为壳幔密度;η1为瞬时粘滞系数;η2为稳态粘滞系数。
    下载: 导出CSV

    表  2  邢台地震序列(MS≥6.0)的震源参数及同震位错模型

    发震时间震中位置震级
    MS
    震源深度/
    km
    断层破裂长度/
    km
    破裂宽度/
    km
    走向滑动量/
    m
    倾向滑动量/
    m
    摩擦系数
    纬度/(˚N)经度/(˚E)走向/(˚)倾角/(˚)滑动角/(˚)
    1966-03-0837.35114.926.81027801704016−0.925−0.1630.4
    1966-03-2237.50115.086.794901753514−0.805−0.070
    1966-03-2237.53115.057.2919901806625−1.5570
    1966-03-2637.68115.276.21529070159198−0.427−0.164
    1966-03-2937.35115.036.0254772231560.334−0.142
    下载: 导出CSV
  • [1] King G C P, Stein R S, Lin J. Static stress changes and the triggering of earthquakes[J]. Bulletin of the Seismological Society of America, 1994, 84(3): 935-953.
    [2] Harris R A. Introduction to special section: stress triggers, stress shadows, and implications for seismic hazard[J]. Journal of Geophysical Research: Solid Earth, 1998, 103(B10): 24347-24358. doi:  10.1029/98JB01576
    [3] 万永革, 沈正康, 曾跃华, 等. 青藏高原东北部的库仑应力积累演化对大地震发生的影响[J]. 地震学报, 2007, 29(2): 115-129. doi:  10.3321/j.issn:0253-3782.2007.02.001
    [4] 邵志刚, 周龙泉, 蒋长胜, 等. 2008年汶川MS8.0地震对周边断层地震活动的影响[J]. 地球物理学报, 2010, 53(8): 1784-1795. doi:  10.3969/j.issn.0001-5733.2010.08.004
    [5] 单斌, 熊熊, 金笔凯, 等. 松潘-甘孜块体东北端强震间相互作用及地震危险性研究[J]. 地球物理学报, 2012, 55(7): 2329-2340. doi:  10.6038/j.issn.0001-5733.2012.07.018
    [6] 程佳, 徐锡伟, 刘杰. 2014年鲁甸6.5级地震成因、破裂特征和余震分布特征的库仑应力作用[J]. 地球物理学报, 2016, 59(2): 556-567. doi:  10.6038/cjg20160214
    [7] 徐晶, 季灵运, 姬存伟, 等. 川滇菱形块体东边界断裂带内库仑应力演化及危险性[J]. 地震地质, 2017, 39(3): 451-470. doi:  10.3969/j.issn.0253-4967.2017.03.002
    [8] 万永革, 吴忠良, 周公威, 等. 几次复杂地震中不同破裂事件之间的“应力触发”问题[J]. 地震学报, 2000, 22(6): 568-576. doi:  10.3321/j.issn:0253-3782.2000.06.002
    [9] 陈连旺, 陆远忠, 刘杰, 等. 1966年邢台地震引起的华北地区应力场动态演化过程的三维粘弹性模拟[J]. 地震学报, 2001, 23(5): 480-491. doi:  10.3321/j.issn:0253-3782.2001.05.004
    [10] 李守勇, 邹立晔, 孙晴, 等. 1966年邢台地震静态库仑应力触发分析[J]. 大地测量与地球动力学, 2016, 36(10): 918-921, 926.
    [11] 万永革, 沈正康, 曾跃华, 等. 唐山地震序列应力触发的粘弹性力学模型研究[J]. 地震学报, 2008, 30(6): 581-593. doi:  10.3321/j.issn:0253-3782.2008.06.004
    [12] Wang R J, Lorenzo-Martín F, Roth F. PSGRN/PSCMP-a new code for calculating co-and post-seismic deformation, geoid and gravity changes based on the viscoelastic-gravitational dislocation theory[J]. Computers & Geosciences, 2006, 32(4): 527-541.
    [13] 邵志刚, 傅容珊, 薛霆虓, 等. 以Burgers体模型模拟震后粘弹性松弛效应[J]. 大地测量与地球动力学, 2007, 27(5): 31-37.
    [14] 徐锡伟, 于贵华, 王峰, 等. 1966年邢台地震群的发震构造模型——新生断层形成? 先存活断层摩擦粘滑?[J]. 中国地震, 2000, 16(4): 364-378. doi:  10.3969/j.issn.1001-4683.2000.04.008
    [15] 王椿镛, 段永红, 吴庆举, 等. 华北强烈地震深部构造环境的探测与研究[J]. 地震学报, 2016, 38(4): 511-549.
    [16] 邵学钟, 张家茹, 章思亚, 等. 邢台地震区深部构造背景的地震转换波探测和研究[J]. 地球物理学报, 1993, 36(5): 609-620. doi:  10.3321/j.issn:0001-5733.1993.05.008
    [17] 侯遵泽, 杨文采, 于常青. 华北克拉通地壳三维密度结构与地质含义[J]. 地球物理学报, 2014, 57(7): 2334-2343. doi:  10.6038/cjg20140727
    [18] 姜文亮, 张景发. 首都圈地区精细地壳结构——基于重力场的反演[J]. 地球物理学报, 2012, 55(5): 1646-1661. doi:  10.6038/j.issn.0001-5733.2012.05.022
    [19] 柳畅, 石耀霖, 郑亮, 等. 三维黏弹性数值模拟华北盆地地震空间分布与构造应力积累关系[J]. 地球物理学报, 2012, 55(12): 3942-3957. doi:  10.6038/j.issn.0001-5733.2012.12.007
    [20] Ryder I, Bürgmann R, Pollitz F. Lower crustal relaxation beneath the Tibetan Plateau and Qaidam Basin following the 2001 Kokoxili earthquake[J]. Geophysical Journal International, 2001, 187(2): 613-630.
    [21] 尹凤玲, 蒋长胜, 韩立波, 等. 红河断裂带库仑应力演化及未来地震危险性估计[J]. 地球物理学报, 2018, 61(1): 183-198. doi:  10.6038/cjg2018L0369
    [22] 河北省地震局. 一九六六年邢台地震[M]. 北京: 地震出版社, 1986.
    [23] 龙锋, 闻学泽, 徐锡伟. 华北地区地震活断层的震级-破裂长度、破裂面积的经验关系[J]. 地震地质, 2006, 28(4): 511-535. doi:  10.3969/j.issn.0253-4967.2006.04.001
    [24] 朱航, 闻学泽. 小江-则木河断裂带大地震序列的静应力触发作用[J]. 地球科学, 2012, 37(1): 199-206.
    [25] 沈正康, 万永革, 甘卫军, 等. 华北地区700年来地壳应力场演化与地震的关系研究[J]. 中国地震, 2004, 20(3): 211-228. doi:  10.3969/j.issn.1001-4683.2004.03.001
    [26] Deng J S, Sykes L R. Evolution of the stress field in southern California and triggering of moderate-size earthquakes: A 200-year perspective[J]. Journal of Geophysical Research: Solid Earth, 1997, 102(B5): 9859-9886. doi:  10.1029/96JB03897
  • [1] 靳雅敏, 于新昌.  唐山地震余震区地震活动性综合评判 . 华北地震科学, 1987, 5(s1): 56-62.
    [2] 彭克中, 李沙白.  邢台地震与我国强震观测工作 . 华北地震科学, 1987, 5(s1): 34-48.
    [3] 郭文.  凹凸型发震构造与邢台地震 . 华北地震科学, 1987, 5(s1): 224-230.
    [4] 易学发, 许国昌.  西安地区活动断裂带及地震危险性探讨 . 华北地震科学, 1987, 5(4): 85-92.
    [5] 黄雪范.  华北地震危险性的预测 . 华北地震科学, 1987, 5(s1): 112-116.
    [6] 靳雅敏, 于新昌, 扈志刚.  由小震应力场的动态变化估计河北及邻近地区的地震危险性 . 华北地震科学, 1987, 5(1): 83-89.
    [7] 林邦惠, 董大南.  正交和斜交破裂过程的研究及其在邢台地震的应用 . 华北地震科学, 1986, 4(2): 16-25.
    [8] 侯立臣.  邢台地震是我国地震事业发展的里程碑——纪念邢台地震二十周年 . 华北地震科学, 1986, 4(2): 5-10.
    [9] 靳雅敏, 于新昌.  邢台地震序列趋势特征 . 华北地震科学, 1986, 4(4): 120-126.
    [10] 常克贵.  应用极值理论对鄂尔多斯北缘地震危险性的估计 . 华北地震科学, 1985, 3(2): 30-40.
  • 加载中
图(4) / 表 (2)
计量
  • 文章访问数:  2140
  • HTML全文浏览量:  1295
  • PDF下载量:  20
出版历程
  • 收稿日期:  2021-01-08
  • 网络出版日期:  2021-09-27
  • 刊出日期:  2021-10-20

邢台老震区库仑应力演化及地震危险性分析

doi: 10.3969/j.issn.1003-1375.2021.04.008
    基金项目:  河北省地震局强震发震构造与机理研究创新团队(DZ20180319009);2019年震情跟踪定向工作任务(2019010122);河北省地震科技星火计划青年项目(DZ20190419028)
    作者简介:

    陈婷(1990—),女,山东聊城人,工程师,主要从事数字地震学及中长期预测等研究. E-mail:chenting_ouc@163.com

  • 中图分类号: P315.7

摘要: 基于Burgers粘弹性层状模型,综合考虑同震位错和震后粘滞松弛效应,利用PSGRN/PSCMP程序计算地震破裂面上产生的库仑破裂应力变化,分析邢台地震序列中每次6级以上地震前的所有6级以上地震对本次地震的综合库仑破裂应力变化,探讨强震间的应力触发关系。结果表明:除3月26日6.2级地震外,其余地震均发生在库仑应力增加的应力触发区,与之前的6级以上地震存在很好的应力触发关系;根据邢台地震序列的库仑应力演化与5级以上地震活动关系讨论邢台老震区的地震危险性,认为5级以上地震发生在应力触发区概率达到71.9%,邢台老震区中强地震活动与库仑应力变化有很好的相关性,可将库仑应力动态演化特征应用到本区地震的中长期预测中,为后续中强震发生地点提供参考。

English Abstract

陈婷,章阳,王晓山,等. 邢台老震区库仑应力演化及地震危险性分析[J]. 华北地震科学,2021, 39(4):50-56. doi:10.3969/j.issn.1003−1375.2021.04.008
引用本文: 陈婷,章阳,王晓山,等. 邢台老震区库仑应力演化及地震危险性分析[J]. 华北地震科学,2021, 39(4):50-56. doi:10.3969/j.issn.1003−1375.2021.04.008
CHEN Ting,ZHANG Yang,WANG Xiaoshan,et al. Coulomb Stress Evolution and Seismic Hazard Analysis on the Xingtai Old Earthquake Region[J]. North China Earthquake Sciences,2021, 39(4):50-56. doi:10.3969/j.issn.1003−1375.2021.04.008
Citation: CHEN Ting,ZHANG Yang,WANG Xiaoshan,et al. Coulomb Stress Evolution and Seismic Hazard Analysis on the Xingtai Old Earthquake Region[J]. North China Earthquake Sciences,2021, 39(4):50-56. doi:10.3969/j.issn.1003−1375.2021.04.008
    • 静态库仑应力变化最早应用在分析大地震对余震分布的影响中,研究认为当静态库仑应力改变量超过0.01 MPa(0.1 bar)就可以影响余震的发生位置[1]。强震或大地震引起的库仑应力变化不仅对余震的发生有影响,更重要的是对周边断层上的强震同样能产生触发作用,使处于库仑应力增加区域的地震提前发生[1-2]。近年来,国内众多学者综合分析断裂带上库仑应力动态演化与强震活动的关系,进而分析区域地震危险性[3-7]

      1966年3月邢台震群共发生5次6.0级以上地震,前人对邢台地震的库仑应力变化展开大量研究。万永革等[8]研究3月8日MS6.8地震对3月22日MS7.2地震的触发作用;陈连旺等[9]研究邢台地震长时间尺度的库仑应力变化;李守勇等[10]分析邢台地震序列2个地震间的相互触发关系,由于使用弹性模型计算产生较大误差,7.2级地震以后的触发关系不好;万永革等[11]研究认为粘弹性模型计算效果优于完全弹性;Wang R等[12]编写了粘弹性应力应变程序PSGRN/PSCMP程序,该程序基于Burgers粘弹性模型[13],考虑地球重力对形变场的影响,可计算同震和震后长时间尺度的应力变化。

      前人对邢台地震库仑应力的计算基于弹性模型假设和粘弹性有限元模型,研究内容也仅限于2个地震间的触发分析。为了更加接近岩石圈的真实地壳结构,本文基于Burgers粘弹性模型,综合考虑同震位错和震后粘滞松弛效应,采用PSGRN/PSCMP程序计算地震破裂产生的库仑应力变化,分析邢台地震序列中6级以上地震受该地震之前所有6级以上地震的综合库仑应力的叠加作用。结合该区5级以上地震活动,分析库仑应力演化与中强地震活动的关系,为区域地震危险性分析乃至中长期地震预报提供一种有效途径。

    • 1966年3月邢台发生震群型地震,地震主要发生在华北平原西部的束鹿断陷盆地,其中6.0级以上地震5次,分别为3月8日MS6.8地震、3月22日MS6.7和MS7.2地震、3月26日MS6.2地震和3月29日MS6.9地震(图1)。震群总体呈N35˚E的椭圆状分布,长轴约110 km,短轴约60 km。7.2级地震位于震区中部,其余地震等间距地对称分布在南北两侧。邢台震群是单个深断裂上端点向上撕裂性破裂扩展、斜列状深断裂间相互作用导致局部应力迁移和对相邻断裂的加载而诱发的三维破裂过程[14]。震源断层均以高角度走滑分量为主,最大主应力为NEE-SWW向,最小主应力为NNW-SSE向[15]

      图  1  邢台老震区地震地质背景

    • 根据库仑应力破裂准则,断层面上的库仑应力变化近似可以表示为

      $$ {{\Delta }}CFS{\rm{ = \Delta }}\tau {\rm{ + }}\mu {\rm{'\Delta }}{\sigma _{\rm{n}}} $$ (1)

      式中:△CFS为库仑应力变化量;△$\tau $为断层面上的剪应力变化量;△${\sigma _n}$为断层面上的正应力变化量;$\mu '$为视摩擦系数,依靠介质的应力变化率,取值不同只会影响应力变化的大小,但对计算库仑应力变化的空间分布影响不大。

      本文采用汪荣江研究员开发的PSGRN/PSCMP程序计算强震的震后库仑应力变化,模拟同震响应和震后松弛效应。

    • 地壳介质模型采用分层半无限空间粘弹性模型,分为完全弹性的上地壳和中地壳、Burgers粘弹性的下地壳和上地幔。壳幔分层结构、P波速度和S波速度参考邵学钟等[16]的研究成果,地层介质密度参考侯遵泽等[17]的研究结果。由于他们未给出上地幔层的密度,此层密度参考姜文亮等[18]的密度和P波速度的经验公式求得。根据柳畅等[19]研究的邢台震区岩石介质粘滞系数随深度的变化结果,邢台地区下地壳的稳态粘滞系数取为3.0×1020 Pa∙s,上地幔的粘滞系数取为2.0×1021 Pa∙s,瞬时与稳态粘滞系数一般相差1~1.5个数量级[20-21],因此下地壳的瞬时粘滞系数取为1.0×1019 Pa∙s。计算库仑应力变化时的具体地壳介质模型见表1

      表 1  地壳介质模型

      层名厚度/kmVP/(km∙s−1VS/(km∙s−1ρ/(kg∙m−3η1/(Pa∙s)η2/(Pa∙s)
      上地壳3.53.201.552 600完全弹性
      4.56.003.452 670
      5.56.203.452 670
      中地壳8.06.653.852 800
      下地壳5.56.803.953 1001.0×10193.0×1020
      5.06.253.253 1001.0×10193.0×1020
      5.07.304.203 1001.0×10193.0×1020
      上地幔8.104.653 36002.0×1021
       注:VP为P波速度;VS为S波速度;ρ为壳幔密度;η1为瞬时粘滞系数;η2为稳态粘滞系数。
    • 邢台地震序列的5次6级以上地震断层参数选取河北省地震局[22]给出的震源机制解。同震位错模型参数中,破裂长度和破裂宽度参考龙锋等[23]经验公式。走向滑动量参考朱航等[24]走滑地震地表破裂长度和平均同震位错的经验公式,估算出走向滑动量后,根据断层面的滑动角估算出倾向滑动量。改变摩擦系数只会对库仑应力值的大小有影响,对库仑应力变化的空间图式影响不大[25],Deng和Sykes等[26]研究发现摩擦系数µ´取值在0~0.6之间时,地震产生的库仑应力对后续地震有较好的应力触发。本文计算中摩擦系数µ´参考沈正康等[25]计算华北地区库仑应力时的取值取为0.4。具体参数见表2

      表 2  邢台地震序列(MS≥6.0)的震源参数及同震位错模型

      发震时间震中位置震级
      MS
      震源深度/
      km
      断层破裂长度/
      km
      破裂宽度/
      km
      走向滑动量/
      m
      倾向滑动量/
      m
      摩擦系数
      纬度/(˚N)经度/(˚E)走向/(˚)倾角/(˚)滑动角/(˚)
      1966-03-0837.35114.926.81027801704016−0.925−0.1630.4
      1966-03-2237.50115.086.794901753514−0.805−0.070
      1966-03-2237.53115.057.2919901806625−1.5570
      1966-03-2637.68115.276.21529070159198−0.427−0.164
      1966-03-2937.35115.036.0254772231560.334−0.142
    • 根据表1的地壳介质模型和表2的震源参数及同震位错模型,采用基于粘弹性模型的PSGRN/PSCMP程序计算库仑应力变化,研究邢台地震序列中6级以上地震的综合库仑应力触发情况(图2)。

      图  2  邢台地震序列强震间的库仑应力触发作用

      研究强震间应力触发关系时,以被触发地震的震源深度作为计算库仑应力变化的深度,被触发地震的断层参数作为接收断层。以MS6.7地震的断层参数作为主接收断层,计算9 km深度处MS6.8地震库仑破裂应力变化(图2a),MS6.7地震位于地震库仑破裂应力增加的应力触发区,震中位置的库仑破裂应力为0.4 bar,即MS6.8地震有利于触发MS6.7地震;以MS7.2地震的断层参数作为主接收断层,计算9 km深度处MS6.8地震和MS6.7地震综合库仑破裂应力变化(图2b),MS7.2地震发生在地震综合库仑破裂应力增加的应力触发区,震中位置的库仑破裂应力为5.9 bar,即MS6.8地震和MS6.7地震对MS7.2地震有触发作用;以MS6.2地震的断层参数作为主接收断层,计算15 km深度处MS6.8地震、MS6.7地震及MS7.2地震综合库仑破裂应力变化(图2c),MS6.2地震处于地震综合库仑破裂应力减少的应力影区,震中位置的库仑破裂应力为-3 bar,即MS6.8地震、MS6.7地震和MS7.2地震推迟了MS6.2地震的发生;以MS6.0地震的断层参数作为主接收断层,计算25 km深度处MS6.8地震、MS6.7地震、MS7.2地震及MS6.2地震综合库仑破裂应力变化(图2d),MS6.0地震在地震综合库仑破裂应力增加的应力触发区发生,震中位置的库仑破裂应力为2.5 bar,即MS6.8地震、MS6.7地震、MS7.2地震及MS6.2地震对MS6.0地震有触发作用。

      3月22日MS6.7地震之后,邢台地震序列4个6级以上地震中有3个受综合库仑应力触发作用,触发率为75%。李守勇等[10]基于弹性模型计算邢台地震序列中4个6级以上地震受前个地震的库仑应力触发情况,研究认为有2个地震被触发,触发率为50%。本文采用粘弹性模型模拟邢台老震区地下介质结构,与地下真实情况更加接近,计算6级以上地震库仑应力时综合考虑之前发生的所有6级以上地震,计算得到的库仑应力对地震的触发率高于基于弹性模型计算6级以上地震受之前单个6级以上地震的库仑应力触发率。

    • 计算邢台老震区库仑应力动态演化,选取区域内主断层参数为接收断层[22]:走向24.6°,倾角81.2°,滑动角141.4°。计算从1966年3月29日MS6.0地震后10 km深度处的同震、震后100年、震后300年以及震后500年的库仑应力变化(图3)。

      图  3  各个时间点邢台地震序列库仑应力变化

      图3分别为3月29日MS6.0地震后同震响应、震后100年、震后300年以及震后500年的库仑应力变化图。地震发生时在其附近区域释放应力产生应力影区,在两端造成应力加载产生应力触发区,呈现“蝴蝶状”分布特征。从图中可以看出,随着时间的演化,应力影区和应力触发区都在向外扩散,库仑应力长时间尺度变化主要与地层的粘滞系数有关,粘滞系数越大,应力传递速度越慢,达到稳态时间越长,应力变化幅度越小;反之,粘滞系数越小,传递越快,达到稳态时间越短,应力变化幅度越大[21]。长期构造加载通常会使地震断层接近其破裂应力水平,整个邢台老震区内北东部位、西南部位、西侧和东侧均位于库仑应力增加的应力触发区,后面应重点关注这些区域有发生强震的可能。

    • 通常认为5级以上地震更具有破坏性,为了更好地分析邢台老震区的地震危险性,探讨库仑应力演化与5级以上地震活动关系,统计邢台地震后发生的5级以上地震受库仑应力变化的触发情况。1966年3月8日MS6.8地震发生至今,总共记录到28个5~6级地震,其中8个发生在3月8日MS6.8地震与3月22日MS6.7地震间,11个发生在3月22日MS7.2地震与3月26日MS6.2地震间,1个发生在3月26日MS6.2地震与3月29日MS6.0地震间,8个发生在3月29日MS6.0地震至今。计算各个时间段的库仑应力变化,统计5到6级地震发生在库仑破裂应力触发区(≥0.1 bar)的频次,进而分析5级以上地震活动与库仑应力演化的关系(图4)。

      图  4  各个时间段库仑应力变化与5级地震的触发关系

      计算库仑应力变化时选取深度为10 km。图4a计算3月8日MS6.8地震库仑应力变化对3月8—22日之间5级地震的触发关系,8个地震中5个震中位于应力触发区,3个震中位于应力影区,受库仑破裂应力触发率为62.5%;图4b计算3月8日MS6.8地震、3月22日MS6.7地震和3月22日MS7.2地震综合库仑应力变化对3月22—26日之间5级地震的触发关系,11个地震中5个震中位于应力触发区,6个震中位于应力影区,受库仑破裂应力触发率为45.5%;图4c计算3月8日MS6.8地震、3月22日MS6.7地震、3月22日MS7.2地震和3月26日MS6.2地震综合库仑应力变化对3月26—29日之间5级地震的触发关系,该时间段内的这个地震位于应力触发区,受库仑破裂应力触发率为100%;图4d计算3月8日MS6.8地震、3月22日MS6.7地震、3月22日MS7.2地震、3月26日MS6.2地震和3月29日MS6.0地震综合库仑应力变化对3月29日至今5级地震的触发关系,8个地震中6个震中位于应力触发区,2个震中位于应力影区,受库仑破裂应力触发率为75%。

      综合以上结果,28个5.0$ \leqslant $MS<6.0地震有17个受库仑破裂应力触发,综合触发率为60.7%。结合前文研究的6级以上地震的触发关系,3月8日MS6.8地震发生后,32个5级以上地震中有20个发生在应力触发区,12个发生在应力影区,地震发生在应力触发区的概率为62.5%。12个未被触发的地震中有3个发生在应力触发区的边缘,原因可能有以下3点:①历史地震破裂参数(震中位置、走向、倾角、滑动角)估计产生误差;②受限于观测资料和接收断层参数的选取产生误差;③破裂带区域次级破裂的结果。如果考虑以库仑破裂应力>0[25]作为应力触发区,则这些地震将被触发。排除这3个库仑破裂应力触发边缘的地震,综合触发率可达到71.9%,可以为后续破坏性地震的发生提供地点参考。

    • 邢台地震序列6级以上强震间的综合触发关系较好,3月8日MS6.8地震之后4个6级以上地震中仅有1个地震(3月26日MS6.2地震)落入库仑应力影区,其余3个地震均受库仑应力触发作用,且5级以上地震发生在应力触发区概率达到71.9%,说明邢台老震区库仑应力变化对中强地震的发生触发作用明显,可将库仑应力动态演化应用到本区地震中长期危险性预测中。由于邢台老震区粘滞系数大,1966年3月29日最后一个6级地震发生后,小尺度时间范围内库仑应力动态变化不明显,库仑应力变化呈现“蝴蝶状”分布特征,地震发生后在其附近区域释放应力产生应力影区,在两端造成应力加载产生应力触发区。邢台老震区的北东部位、西南部位、西侧和东侧均位于库仑应力增加的应力触发区,应重点关注这些区域发生中强震可能性。

      在计算库仑应力变化时,存在着很多误差。由于没有MS6.8地震之前的库仑破裂应力变化数据,假定MS6.8地震之前的累积库仑破裂应力为0;建立地壳层状模型时,由于缺少本地区精细的三维速度模型及粘滞系数,模型参数均是根据前人研究成果选取的;对于历史地震破裂参数、初始应力场和长期加载场认识有限,对于库仑应力估计比较粗糙。以上这些都会给计算带来难度,虽不能准确地预测破裂时间,但根据库仑应力动态演化仍然能够为地震中长期预测提供重要信息。

参考文献 (26)

目录

    /

    返回文章
    返回