logo

Journals Books
Advanced Search
Sign in Register

中文

SCIENTIA SINICA Terrae
Download PDF
SCIENTIA SINICA Terrae, Volume 47, Issue 11: 1329-1338(2017) https://doi.org/10.1360/N072017-00268

2017年M7.0级九寨沟地震同震库仑应力变化及其与2008年汶川地震的关系

单斌1,*, 郑勇1, 刘成利1, 谢祖军1, 孔军2
More info
  • ReceivedAug 15, 2017
  • AcceptedOct 13, 2017
  • PublishedOct 25, 2017
Previous Table of Contents Next
Share
Rating

Abstract

2017年8月8日, 四川九寨沟地区发生M7.0级强震, 造成了重大人员伤亡和财产损失. 判断余震可能的发生趋势、圈定震区周边未来可能的强震危险区, 对于指导灾民安置、灾后重建等具有重要的参考价值. 本文利用弹性位错理论和分层岩石圈模型, 计算了汶川地震和九寨沟地震引起的应力场变化, 讨论了九寨沟地震与汶川地震之间的关系、九寨沟地震对余震分布以及周边主要活动断层上应力积累的影响. 结果表明: (1) 汶川地震导致的同震及震后应力变化显著增强了九寨沟地震震源处的应力积累, 有效地促进了九寨沟地震的发生; (2) 九寨沟地震所导致的同震库仑应力场变化较好地解释了目前所记录到的余震的空间分布; (3) 东昆仑断裂玛曲-荷叶段、虎牙断裂北端受九寨沟地震影响, 应力积累状态提升, 地震危险性有所增强.


Funded by

国家自然科学基金项目(41674106)


Acknowledgment

感谢中国地震局监测预报处对本工作的大力支持. 感谢中国地震局地球物理研究所“国家数字测震台网数据备份中心”和四川省地震局提供了近震地震波形资料, IRIS为本研究提供了波形数据. 在本文完成过程中, 主编、编委和审稿人提出建设性的修改意见, 特此感谢. 本文图件均用GMT软件绘制.


References

[1] 邓启东, 张培震, 冉勇康, 杨晓平, 闵伟, 楚全芝. 2002. 中国活动构造基本特征. 中国科学D辑: 地球科学, 32: 1020−1030. Google Scholar

[2] 罗艳, 曾祥方, 倪四道. 2013. 震源深度测定方法研究进展. 地球物理学进展, 28: 2309−2321. Google Scholar

[3] 单斌, 熊熊, 郑勇, 刁法启. 2009. 2008年5月12日Mw7.9汶川地震导致的周边断层应力变化. 中国科学D辑: 地球科学, 39: 537−545. Google Scholar

[4] 单斌, 熊熊, 郑勇, 金笔凯, 刘成利, 谢祖军, 许厚泽. 2013. 2013年芦山地震导致的周边断层应力变化及其与2008年汶川地震的关系. 中国科学: 地球科学, 43: 1002−1009. Google Scholar

[5] 邵志刚, 傅容珊, 薛霆虓, 黄建华. 2008. 昆仑山Ms8.1级地震震后变形场数值模拟与成因机理探讨. 地球物理学报, 51: 805−816. Google Scholar

[6] 王卫民, 赵连锋, 李娟, 姚振兴. 2008. 四川汶川8.0级地震震源过程. 地球物理学报, 51: 1403−1410. Google Scholar

[7] 谢祖军, 郑勇, 姚华建, 房立华, 张勇, 刘成利, 王毛毛, 单斌, 张会平, 任俊杰, 季灵运, 宋美琴. 2017. 2017年Ms7.0级九寨沟地震震源性质及发震构造初步分析. 中国科学: 地球科学, doi: 10.1360/N072017-00279. Google Scholar

[8] 张培震, 邓启东, 张国民, 马瑾, 甘卫军, 闵伟, 毛凤英, 王琪. 2003. 中国大陆的强震活动与活动地块. 中国科学D辑: 地球科学, 33(增刊): 12−20. Google Scholar

[9] Bai D, Unsworth M J, Meju M A, Ma X, Teng J, Kong X, Sun Y, Sun J, Wang L, Jiang C, Zhao C, Xiao P, Liu M. Crustal deformation of the eastern Tibetan plateau revealed by magnetotelluric imaging. Nat Geosci, 2010, 3: 358-362 CrossRef ADS Google Scholar

[10] Bassin C, Laske G, Masters G. 2000. The current limits of resolution for surface wave tomography in North America. EOS Trans AGU, 81: 897. Google Scholar

[11] Deng J, Sykes L R. Triggering of 1812 Santa Barbara Earthquake by a Great San Andreas Shock: Implications for future seismic hazards in southern California. Geophys Res Lett, 1996, 23: 1155-1158 CrossRef ADS Google Scholar

[12] Freed A M. Earthquake triggering by static, dynamic, and postseismic stress transfer. Annu Rev Earth Planet Sci, 2005, 33: 335-367 CrossRef ADS Google Scholar

[13] Harris R A, Simpson R W, Reasenberg P A. Influence of static stress changes on earthquake locations in southern California. Nature, 1995, 375: 221-224 CrossRef ADS Google Scholar

[14] Harris R A, Simpson R W. Suppression of large earthquakes by stress shadows: A comparison of Coulomb and rate-and-state failure. J Geophys Res, 1998, 103: 24439-24451 CrossRef ADS Google Scholar

[15] Heidbach O, Ben-Avraham Z. Stress evolution and seismic hazard of the Dead Sea Fault System. Earth Planet Sci Lett, 2007, 257: 299-312 CrossRef ADS Google Scholar

[16] Ji C, Wald D J, Helmberger D V. 2002. Source description of the 1999 Hector Mine, California earthquake, Part I: Wavelet domain inversion theory and resolution analysis. Bull Seismol Soc Amer, 92: 1192−1207. Google Scholar

[17] King G C P, Stein R S, Lin J. 1994. Static stress changes and the triggering of earthquakes. Bull Seismol Soc Amer, 84: 935−953. Google Scholar

[18] Nalbant S S, Steacy S, Sieh K, Natawidjaja D, McCloskey J. Seismology: Earthquake risk on the Sunda trench. Nature, 2005, 435: 756-757 CrossRef PubMed ADS Google Scholar

[19] Okada Y. 1992. Internal deformation due to shear and tensile faults in a half-space. Bull Seismol Soc Amer, 82: 1018−1040. Google Scholar

[20] Pollitz F, Vergnolle M, Calais E. 2003. Fault interaction and stress triggering of twentieth century earthquakes in Mongolia. J Geophys Res, 108: 2503. Google Scholar

[21] Ren J, Xu X, Yeats R S, Zhang S. Millennial slip rates of the Tazang fault, the eastern termination of Kunlun fault: Implications for strain partitioning in eastern Tibet. Tectonophysics, 2013, 608: 1180-1200 CrossRef ADS Google Scholar

[22] Shan B, Xiong X, Wang R, Zheng Y, Yang S. Coulomb stress evolution along Xianshuihe-Xiaojiang Fault System since 1713 and its interaction with Wenchuan earthquake, May 12, 2008. Earth Planet Sci Lett, 2013, 377-378: 199-210 CrossRef ADS Google Scholar

[23] Shen W, Ritzwoller M, Kang D, Kim Y, Lin F, Ning J, Wang W, Zhou L. 2016. A seismic reference model for the crust and uppermost mantle beneath China from surface wave dispersion. Geophys J Int, 206: 954−979. Google Scholar

[24] Stein R S. Earthquake conversations. Sci Am, 2003, 288: 72-79 CrossRef ADS Google Scholar

[25] Toda S, Stein R S, Reasenberg P A, Dieterich J H, Yoshida A. Stress transferred by the 1995 Mw = 6.9 Kobe, Japan, shock: Effect on aftershocks and future earthquake probabilities. J Geophys Res, 1998, 103: 24543-24565 CrossRef ADS Google Scholar

[26] Wang R. 1999. A simple orthonormalization method for the stable and efficient computation of Green’s functions. Bull Seismol Soc Amer, 89: 733−741. Google Scholar

[27] Wang R, 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. Comp Geosciences, 2006, 32: 527-541 CrossRef ADS Google Scholar

[28] Xiong X, Shan B, Zheng Y, Wang R. Stress transfer and its implication for earthquake hazard on the Kunlun Fault, Tibet. Tectonophysics, 2010, 482: 216-225 CrossRef ADS Google Scholar

[29] Zhang J. 2017. On fault evidence for a large earthquake in the late fifteenth century, Eastern Kunlun fault, China. J Seismol, 1−9. Google Scholar

[30] Zhao G, Unsworth M J, Zhan Y, Wang L, Chen X, Jones A G, Tang J, Xiao Q, Wang J, Cai J, Li T, Wang Y, Zhang J. Crustal structure and rheology of the Longmenshan and Wenchuan Mw 7.9 earthquake epicentral area from magnetotelluric data. Geology, 2012, 40: 1139-1142 CrossRef ADS Google Scholar

[31] Zhao L S, Helmberger D V. 1994. Source estimation from broadband regional seismograms. Bull Seismol Soc Amer, 84: 91−104. Google Scholar

[32] Zheng Y, Yang Y, Ritzwoller M H, Zheng X, Xiong X, Li Z. Crustal structure of the northeastern Tibetan plateau, the Ordos block and the Sichuan basin from ambient noise tomography. Earthq Sci, 2010, 23: 465-476 CrossRef ADS Google Scholar

[33] Zhu L P, Helmberger D V. 1996. Advancement in source estimation techniques using broadband regional seismograms. Bull Seismol Soc Amer, 86: 1634−1641. Google Scholar

  • Click to view Download high-quality image

    图1

    九寨沟地震及周边地区构造背景

    红色沙滩球为2017年九寨沟M7.0级地震震源机制解, 灰色沙滩球为2008年Mw7.9级汶川地震和2013年Ms7.0级芦山地震震源机制解, 绿色圆圈为该地区内主要城镇, 黑色实线为汶川地震地表破裂面. 黑色细线为该区域内主要活动断层, 参考新构造的研究结果(邓启东等, 2002; 张培震等, 2003).

  • Click to view Download high-quality image

    图2

    九寨沟地震动态破裂过程反演结果

    (a) 图中白色箭头为破裂滑移方向和大小, 黑色等高线为破裂传播时间(单位为秒), 水平坐标为332°走向, 纵坐标为沿断层破裂方向距离地表的距离; (b) 地震矩释放的时间函数

  • Click to view Download high-quality image

    图3

    远震波形拟合情况

    黑线为观测波形, 红线为合成波形; 左侧字符为地震波波形和台站名缩写; 左侧上方数据为方位角, 下方数据为震中距; 每个波形右上方的数字为最大振幅

  • Click to view Download high-quality image

    图4

    汶川地震导致的库仑应力变化

    (a) 同震库仑应力变化; (b) 同震+震后粘弹松弛库仑应力变化. 接收断层参数取strike=332°, dip=84°, rake=5°; 计算深度为6km; 有效摩擦系数取0.4

  • Click to view Download high-quality image

    图5

    九寨沟地震同震静态库仑应力变化场与余震分布的关系

    库仑应力计算时, 有效摩擦系数取0.4, 白色圆点是九寨沟地震余震重定位的结果(谢祖军等, 2017). (a) 深度范围5~20km库仑应力变化的最大值, 绿色线段为九寨沟地震破裂面在地表的投影, 红色虚线为地震剖面的水平投影; (b)、(c) 为沿深度剖面AA′和BB′的库仑应力变化情况及与余震的关系

  • Click to view Download high-quality image

    图6

    九寨沟地震同震静态库仑应力对周边主要活动断层的影响

  • 表1   岩石圈分层结构模型及参数a)

    序号

    深度范围(km)

    P波波速(km s‒1)

    S波波速(km s‒1)

    密度(kg m‒3)

    粘滞系数(Pa s)

    1

    0~23.3

    6.1

    3.55

    2740

    ∞

    2

    23.3~44

    6.3

    3.65

    2780

    5×1017

    3

    44~53

    7.0

    3.99

    2950

    5×1017

    4

    >53

    8.05

    4.47

    3320

    1×1020

    a) 地壳各层厚度、地震波速、密度参考Crust1.0结构模型, 各层粘滞系数的大小参考邵志刚等(2008)数值模拟的研究结果

  • 表2   汶川地震在九寨沟地震震中位置所导致的应力变化(MPa)a)
    深度(km) 时刻

    μ=0.0

    μ=0.2

    μ=0.4

    μ=0.6

    μ=0.8

    6

    同震

    0.018

    0.016

    0.015

    0.013

    0.011

    同震+震后

    0.052

    0.045

    0.038

    0.031

    0.024

    20

    同震

    0.017

    0.015

    0.014

    0.012

    0.011

    同震+震后

    0.047

    0.041

    0.035

    0.029

    0.023

    a)μ为计算中选取的有效摩擦系数

  • 表3   不同粘滞系数对九寨沟地震震中位置库仑应力变化(MPa)的影响a)

    ρ c

    ρm=1×1019

    ρm=1×1020

    1×1017

    0.059

    0.058

    5×1017

    0.042

    0.038

    1×1018

    0.037

    0.034

    1×1019

    0.024

    0.022

    a) 表中ρc、ρm分别为中下地壳和地幔的粘滞系数(Pa s)

  • 表4   相应所研究的主要断层及断层参数a)

    断层名称

    走向(°)

    倾角(°)

    滑动角(°)

    滑动速率(mm/yr)

    有效摩擦系数μ

    岷江

    180

    45~56

    45

    < 1

    0.6~0.8

    虎牙

    150

    75

    45

    1.4

    0.6~0.8

    玛曲-荷叶*

    100~145

    89

    0

    1.4~3.2

    0.6~0.8

    哈南-青山湾

    69

    70

    ‒45

    ‒

    0.1~0.8

    迭部-白龙江

    290

    68~85

    45

    ‒

    0.1~0.8

    龙日坝

    205~229

    60

    135

    5.1 ± 1.2

    0.4~0.6

    青川

    70

    75

    ‒170

    ‒

    0.6~0.8

    礼县-罗家堡

    65

    65

    45

    0.95

    0.6~0.8

    a)*东昆仑断裂玛曲-荷叶段(塔藏断裂)参考Ren等 (2013)的研究结果. ‒ 表示没有可靠的滑移速率用于确定该断层的有效摩擦系数

Copyright 2019 Science China Press Co., Ltd. 《中国科学》杂志社有限责任公司 版权所有

京ICP备18024590号-1