zhenbo

ISSN 2096-7780 CN 10-1665/P

留言板

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

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

利用人工震源实验揭示日本樱岛火山之下结构演变

TomokiTsutsui MasatoIguchi TakeshiTameguri HaruhisaNakamichi 王林

TomokiTsutsui, MasatoIguchi, TakeshiTameguri, HaruhisaNakamichi, 王林. 利用人工震源实验揭示日本樱岛火山之下结构演变[J]. 地震科学进展, 2018, (10): 20-36. doi: 10.3969/j.issn.0253-4975.2018.10.003
引用本文: TomokiTsutsui, MasatoIguchi, TakeshiTameguri, HaruhisaNakamichi, 王林. 利用人工震源实验揭示日本樱岛火山之下结构演变[J]. 地震科学进展, 2018, (10): 20-36. doi: 10.3969/j.issn.0253-4975.2018.10.003

利用人工震源实验揭示日本樱岛火山之下结构演变

doi: 10.3969/j.issn.0253-4975.2018.10.003
详细信息
    通讯作者:

    Tomoki Tsutsui,e-mail: tomoki@gipc.akita-u.ac.jp

  • 中图分类号: P315.2

  • 摘要: 通过人工震源实验,研究地震反射率变化用于监测活火山之下结构演变,认为其与岩浆卸载有关。樱岛火山是日本最活跃的火山之一,也是目前研究的主要目标。每年在火山区域开展人工震源实验,共计完成6期实验。通过对每期实验的数据进行分析,获取了两个地震反射剖面。实验清晰地反映出樱岛东北部深度6.2 km处的地震反射率年变化情况。2009年12月首次岩浆入侵时出现最大反射率,随后逐渐下降,直到2013年12月,与樱岛火山膨胀压缩周期相一致。反射率变化发生在地壳深部嵌入的明显反射区内。位于中间层的夹层结构被作为反射区模型。低速岩浆嵌入中间层及其后续波速的增加可用于解释反射率的变化范围。这也解释了地壳深部岩浆喷发与温度下降的关系。该研究描述了一种用于即时获取岩浆特性以及监测活火山的新方法。

     

  • 图  1  樱岛火山。左上角箭头表示樱岛火山的位置。五角星表示主要的爆破点,圆点表示台站位置。灰色虚线表示岩浆活动通道[16]。KT、MD和SY分别表示北岳峰顶、南岳和昭和火山口。利用日本地理空间信息局(GSI)提供的50 m分辨率的数字栅格地图提取了100 m间距等高线图

    图  2  自2007年1月起伸缩仪记录值和火山落灰记录结果。曲线表示图1中HART的伸缩仪记录月值(由日本京都大学提供)。柱状表示鹿儿岛县每月火山落灰结果(由日本气象厅提供)。箭头表示每次地震实验。罗马数字表示活跃周期

    图  3  每个爆破点的炮孔分布。“+”代表炮孔和2009年原始炮孔。2013年未实施爆破实验的8个爆破点,分别为KURE、KRKE、JGKE、UTOE、UTKE、KMME、FKRE和WIZE

    图  4  2009年实验获取的震动图,除了图5中的爆破点SHRE和KURE。图中已标注正极

    图  5  每年爆破记录。上图为EW测线上SHRE爆破点记录,下图为NS测线上KURE爆破点记录。爆破点KURE未获取2013年12月震动图。箭头表示到达波A。图形排列方式与图4相同

    图  6  年震动图。(a)位于爆破点SHRE附近3.5 km处的站点X087B获取的震动图。图中,A表示清晰的到达波,与来自深度6.2 km的反射相一致。较大振幅出现在到达波后4 s,表明到达波穿过浅部岩层。(b)位于爆破点SHRE附近0.05 km处的站点X040B获取近场震动图

    图  7  初至波与到达波A的振幅均方根比值年变化。图中显示95%置信区间,值非常小,位于符号后

    图  8  数据处理流程。表1表2中罗列了校正参数。表4显示在地震反射分析中所使用的参数

    图  9  基于不同的参考值得到爆破点SHRE年相关性。粗线表示2009年至2011年连续3年实验结果与叠加震动图的最大互相关中值。细线表示2009年至2011年连续3年实验结果与2009年震动图的最大互相关中值。置信区间为95%,值很小,位于标记后

    图  10  参考震动图。上面6幅图是沿着EW测线的爆破记录,下面9幅图是沿着NS测线的爆破记录

    图  11  地震剖面位置。中间点沿着EW测线为0—188,沿着NS测线为0—250。高程如图1所描述

    图  12  深部地壳剖面。上行为EW测线剖面。下行为NS测线剖面。垂直虚线表示剖面的交叉点。参考剖面位于最左侧。参考图的右侧为每年特定剖面。参考剖面中,红色表示正反射,蓝色表示逆反射。在特定剖面中,粉色表示正偏差,绿色表示负偏差。颜色深浅表示相对差异程度(原图及电子版为彩色图)

    图  13  特定震动图。上部为沿着EW测线爆破点SHRE的爆破记录。下部为沿着NS测线爆破点KURE的爆破点记录。箭头表示到达波A

    图  14  特定时间剖面。上行为EW测线。下行为NS测线。箭头表示到达波A

    图  15  偏差剖面方法图解。左上角图表示偏差,右下角图表示偏差剖面图例。偏差剖面中,色调表示相对尺度上沿着色环,由方向确定的极性和由距离确定的偏差度

    图  16  偏差剖面。偏差剖面体现出与参考值的差异。上部表示EW测线的偏差剖面,下部表示NS测线的偏差剖面。相应的参考剖面显示在各行的最左侧。所应用的色调变化值与图15中相同

    图  17  反射区α分布。反射区α显示在有界区域内。实线圆表示有意义的中间点,黑点表示无意义的中间点,这些中间点都是地震剖面中确切的点。图中显示了海岸线和海平面以上500 m等高线

    图  18  2009年和2010年沿着EW剖面包含爆破点URAE的剖面图。(a)特定时间段。(b)偏差剖面。上行的箭头表示到达波A,下行的箭头表示反射区α。每个剖面范围都与之前相应剖面相同

    图  19  α的结构模型。(a)几何结构;(b)不同厚度在6 Hz时波速与表面反射率的关系

    表  1  历年实验与2009年实验振幅比的关系

    爆破点 2009年12月 2010年12月 2011年12月 2012年12月 2013年12月 2014年12月
    KURE 1.000 0.803 0.669 0.674 0.959
    KRKE 1.000 0.990 1.379 1.481 1.244
    JGKE 1.000 1.181 1.326 1.200 1.238
    UTOE 1.000 1.069 1.017 0.922 0.881
    UTKE 1.000 1.097 0.941 1.097 1.069
    KMME 1.000 0.993 1.024 1.142 1.010
    KOME 1.000 1.634 1.367 1.368 1.106 1.696
    FKRE 1.000 0.883 1.678 0.986 1.298
    WIZE 1.000 0.701 0.505 1.390 0.978
    URAE 1.000 1.188
    URNE 1.000 1.061 1.023 1.228 1.138 1.241
    SHRE 1.000 0.994 1.238 1.333 1.163 1.176
    FTME 1.000 0.807 0.538 0.917 0.578 0.536
    KBNE 1.000 1.225 1.122 1.255 1.488 1.453
    下载: 导出CSV

    表  2  历年实验与2009年实验的爆破延时关系(s)

    爆破点 2009年12月 2010年12月 2011年12月 2012年12月 2013年12月 2014年12月
    KURE 0.000 –0.004 –0.007 0.011 0.004
    KRKE 0.000 0.004 0.022 0.034 –0.040
    JGKE 0.000 0.006 0.019 0.017 0.012
    UTOE 0.000 0.008 0.018 0.010 0.006
    UTKE 0.000 0.008 0.020 0.020 0.010
    KMME 0.000 0.003 0.009 0.008 –0.003
    KOME 0.000 –0.002 0.006 0.008 0.008 0.004
    FKRE 0.000 –0.008 0.020 0.026 0.020
    WIZE 0.000 –0.004 –0.010 0.020 0.014
    URAE 0.000 0.001
    URNE 0.000 0.009 0.010 0.016 0.002 0.016
    SHRE 0.000 0.000 0.024 0.032 0.024 0.020
    FTME 0.000 0.016 –0.008 0.020 –0.004 0.000
    KBNE 0.000 0.030 0.024 0.038 0.040 0.040
    下载: 导出CSV

    表  3  历年实验与2009年实验结果的相关关系

    爆破点 2009年12月 2010年12月 2011年12月 2012年12月 2013年12月 2014年12月
    KURE 1.000 0.789 0.486 0.548 0.655
    KRKE 1.000 0.776 0.597 0.631 0.407
    JGKE 1.000 0.742 0.605 0.666 0.578
    UTOE 1.000 0.790 0.690 0.703 0.668
    UTKE 1.000 0.856 0.744 0.796 0.699
    KMME 1.000 0.837 0.830 0.812 0.794
    KOME 1.000 0.809 0.796 0.803 0.719 0.763
    FKRE 1.000 0.751 0.251 0.329 0.287
    WIZE 1.000 0.560 0.400 0.257 0.358
    URAE 1.000 0.894
    URNE 1.000 0.845 0.760 0.726 0.735 0.640
    SHRE 1.000 0.703 0.737 0.670 0.690 0.619
    FTME 1.000 0.848 0.462 0.673 0.446 0.462
    KBNE 1.000 0.668 0.576 0.543 0.429 0.449
    下载: 导出CSV

    表  4  地震反射分析

    序号 参数
    1. 轨迹混合 邻近的3个轨迹
    2. 中间点种类 间距:EW测线为35.5 m,NS测线为23.2 m,EW测线和NS测线的最长距离分别为6 690 m和5 457 m
    3. 带通滤波 2—8 Hz
    4. 增益恢复 自动增益控制时间:0.6 s
    5. 预测去卷积 预测距离:0.11—0.15 s,因子长度:0.16 s
    6. 弱音器 初至波弱音器(所有),地滚波弱音器(特定)
    7. 静态校正 折射和高层,数据:海平面之上372 m
    8. NMO校正 速度函数;表5
    9. CMP叠加 最大叠加数:EW为14次,NS为10次
    10. 偏移 叠后2D基尔霍夫深度偏移,波速模型;表5
    下载: 导出CSV

    表  5  波速函数(引自Tsutsui等[20]的研究成果)

    双向时间
    /s
    叠加波速
    /(m·s–1
    深度
    /m
    地层波速
    /(m·s–1
    0.95 1 459 693 1 459
    1.43 1 880 1 344 2 506
    2.02 2 385 2 409 3 313
    3.41 3 410 5 814 4 495
    4.8 4 283 10 280 5 901
    下载: 导出CSV
  • [1] McNutt S R. Eruption Response and Mitigatio[M].// Sigurdsson H et al.(Eds.). Encyclopedia of Volcanoes. Second ed., Academic Press, New York, 2015: 1070-1071
    [2] Iguchi M,Ishihara K,Surono M, et al. Learn from 2010 eruptions at Merapi and Sinabung volcanoes in Indonesia[J]. Ann. Disaster Prev. Res. Inst.,Kyoto Univ.,2011,54B:185-194
    [3] Iizuka S, Hasegawa I, Ito K, et al. Precise measurements of changes in seismic wave velocities by means of explosion-seismic method[R]. Reports, Geological Survey of Japan, Kawasaki, 1975, 254: 46-74(in Japanese with English abstract)
    [4] Nishimura T,Tanaka S,Yamawaki T, et al. Temporal changes in seismic velocity of the crust around Iwate volcano,Japan,as inferred from analyses of repeated active seismic experiment data from 1998 to 2003[J]. Earth Planets Space,2005,57:491-505 doi: 10.1186/BF03352583
    [5] Poupinet G,Elisworth W L,Frechet J. Monitoring velocity variations in the crust using earthquake doublets:an application to the Calaveras Fault,California[J]. J. Geophys. Res.,1984,89:5719-5731 doi: 10.1029/JB089iB07p05719
    [6] Anggono T,Nishimura T,Sato H, et al. Spatio-temporal changes in seismic velocity associated with the 2000 activity of Miyakejima volcano as inferred from cross-correlation analyses of ambient noise[J]. J. Volcanol. Geotherm. Res.,2012:247-248,93-107
    [7] Matsushima J,Yokota T,Okubo Y, et al. Repeated seismic reflection measurements in the Kakkonda geothermal field[J]. J. Volcanol. Geotherm. Res.,2004,129:343-356 doi: 10.1016/S0377-0273(03)00286-5
    [8] Aramaki S. Formation of the Aira Caldera,Southern Kyushu,—22,000 years ago[J]. J. Geophys. Res.,1984,89:8485-8501 doi: 10.1029/JB089iB10p08485
    [9] Kobayashi T, Miki D, Sasaki H, et al. Geological Map of Sakurajima Volcano[M]. Second ed., Geologcal Survey of Japan, AIST., 2013(in Japanese with English abstract)
    [10] Japan Meteorological Agency. National Catalogue of the Active Volcanoes in Japan[M]. Fourth ed., The Japan Meteorological Agency and the Volcanological Society of Japan, 2013
    [11] Omori F. The Sakurajima eruptions and earthquakes I-VI[J]. Bull. Imp. Earthq. Inv. Comm.,1922,8:525
    [12] Ishihara K. A Quantitative Relation Between the Ground Deformation and the Volcanic Materials Ejected[C]. Abstract 98 IAVCEI Symposium--Arc Volcanism, 1981: 143
    [13] Eto T,Takayama T,Yamamoto K, et al. Reupheaval of the ground surface at the Aira caldera — December 1991 October 1996[J]. Ann. Disaster Prev. Res. Inst.,Kyoto Univ.,1997,40B-1:49-60
    [14] Yamamoto K,Sonoda T,Takayama T, et al. Vertical ground deformation associated with the volcanic activity of Sakurajima Volcano,Japan during 1996-2010 as revealed by repeated precise leveling surveys[J]. Bull. Volcanol. Soc. Japan,2013,58:137-151
    [15] Yamamoto K,Matsushima T,Yoshikawa S, et al. Vertical ground deformation in Sakurajima volcano and around Aira Caldera measured by precise leveling survey conducted in October and November 2013[J]. Ann. Disaster Prev. Res. Inst.,Kyoto Univ.,2013,57B:116-124
    [16] Hidayati Sri,Ishihara K,Iguchi M. Volcano-tectonic earthquakes during the stage of magma accumulation at the Aira Caldera,Southern Kyushu,Japan[J]. Bull. Volcanol. Soc. Japan,2007,52(6):289-309
    [17] Tsuya H,Minakami T. Minor activity of volcano Sakura-zima in October,1939[J]. Bull. Earthq. Res. Inst.,Univ. Tokyo,1940,18:318-339
    [18] Iguchi M,Yokoo A,Tameguri T. Intensity of volcanic explosions at Showa Crater of Sakurajima Volcano[J]. Ann. Disaster Prev. Res. Inst.,Kyoto Univ.,2010,53B:233-250
    [19] Iguchi M,Tameguri T,Ohta Y, et al. Characteristics of volcanic activity at Sakurajima Volcano’s Showa crater during the period 2006 to 2011[J]. Bull. Volcanol. Soc. Japan,2013,58:115-136
    [20] Tsutsui T,Yagi N,Iguchi M, et al. Structure of Northeastern Sakurajima,South Kyushu,Japan,revealed by seismic reflection survey[J]. Bull. Volcanol. Soc. Japan,2013,58:239-250
    [21] Kurashimo E,Hirata N,Morita Y, et al. A Dense Seismic Observation System With Small High-performance Off-line Data Loggers[J]. Zishin,ser. 2,2006,59:107-116
    [22] Yilmaz O. Seismic Data Analysis. Processing, Inversion, and Interpretation of Seismic Data Volume I[M]. Society of Exploration Geophysicists, Tulsa, USA, 2001(1000pp.)
    [23] Aki K, Richards PG. Quantitative Seismology: Theory and Methods, Volume II[M]. Freeman and company, San Francisco, 1980: 661-675
    [24] Clearbout J F. Synthesis of a layered medium from its acoustic transmission response[J]. Geophysics,1968,33:264-269 doi: 10.1190/1.1439927
    [25] Gudmundsson A. Magma chambers:formation,local stresses,excess pressures,and compartments[J]. J. Volcanol. Geotherm. Res.,2012:237-238,19-41
    [26] Miller C F,Furbish D J,Walker B A, et al. Growth of plutons by incremental emplacement of sheets in crystal-rich host:evidence from Miocene intrusions of the Colorado River region,Nevada,USA[J]. Techtonophysics,2011,500(1):65-77 doi: 10.1016/j.tecto.2009.07.011
    [27] Thomson K. Volcanic features of the North Rockall Trough:application of visualization techniques on 3D seismic reflection data[J]. Bull. Volcanol.,2005,67(2):116-128 doi: 10.1007/s00445-004-0363-9
    [28] Okada Y. Internal deformation due to shear and tensile faults in a half-space[J]. Bull. Seismol. Soc. Amer.,1992,82:1018-1040
    [29] Spera F J. Physical properties of magmas[M]// Sigurdsson H(Ed.). Encyclopedia of Volcanoes, Academic Press, New York, 2000: 171-190
    [30] Williams H, McBirney A R.Volcanology[M]. Freeman, Cooper and Co., 1979
    [31] Murase T,McBirney A R. Properties of some common igneous rocks and their melts at high temperatures[J]. Geol. Soc. Am. Bull.,2.0.CO;2">1973,84:3563-3592 doi: 10.1130/0016-7606(1973)84<3563:POSCIR>2.0.CO;2
    [32] Iguchi M. Eruptive Activity of Sakurajima Volcano During the Period From July 2011 to June 2012[M]// Report of Research Program for Preparation Process of Volcanic Eruptions Through the Multiple Discipline Observations 2012. Volcanological Research Center, Kyoto Univ., 2013: 1-8(in Japanese with English abstract)
  • 加载中
图(19) / 表(5)
计量
  • 文章访问数:  538
  • HTML全文浏览量:  196
  • PDF下载量:  10
出版历程
  • 收稿日期:  2016-09-06
  • 修回日期:  2016-09-21
  • 刊出日期:  2018-10-01

目录

    /

    返回文章
    返回
    本系统由北京仁和汇智信息技术有限公司设计开发 百度统计