• ISSN 1673-5722
  • CN 11-5429/P

地震动残差分析盆地附加放大效应:日本关东盆地为例

李春果 王宏伟 温瑞智 任叶飞

马兴全, 王志铄, 孙杰, 李涛, 李源, 蔡颖哲, 夏修军, 赵显刚. 通过断层气氢探测新郑-太康断裂的浅层位置[J]. 震灾防御技术, 2017, 12(2): 354-362. doi: 10.11899/zzfy20170212
引用本文: 李春果, 王宏伟, 温瑞智, 任叶飞. 地震动残差分析盆地附加放大效应:日本关东盆地为例[J]. 震灾防御技术, 2020, 15(4): 684-695. doi: 10.11899/zzfy20200403
Ma Xingquan, Wang Zhishuo, Sun Jie, Li Tao, Li Yuan, Cai Yingzhe, Xia Xiujun, Zhao Xian'gang. Exploration for the Shallow Location of Xinzheng-Taikang Fault via Fault Gas H2[J]. Technology for Earthquake Disaster Prevention, 2017, 12(2): 354-362. doi: 10.11899/zzfy20170212
Citation: Li Chunguo, Wang Hongwei, Wen Ruizhi, Ren Yefei. Basin Extra Amplification Effects from Seismic Ground-Motion Residual Analysis: A Case Study of Kanto Basin, Japan[J]. Technology for Earthquake Disaster Prevention, 2020, 15(4): 684-695. doi: 10.11899/zzfy20200403

地震动残差分析盆地附加放大效应:日本关东盆地为例

doi: 10.11899/zzfy20200403
基金项目: 

中国地震局工程力学研究所基本科研业务 2018B03

国家自然科学基金资助项目 51808514

国家自然科学基金资助项目 51878632

详细信息
    作者简介:

    李春果, 女, 生于1995年。博士研究生。主要从事地震动场地效应研究。E-mail: lcggzl007@163.com

    通讯作者:

    王宏伟, 男, 生于1990年。助理研究员。主要从事地震动特征及预测等方面研究。E-mail: whw1990413@163.com

Basin Extra Amplification Effects from Seismic Ground-Motion Residual Analysis: A Case Study of Kanto Basin, Japan

  • 摘要: 我国地震动预测及地震危险性分析通常仅考虑局部场地浅层岩土层对地震动的放大效应,不能考虑较大范围的地质条件影响,如沉积盆地厚沉积层对地震动的附加放大效应通常被忽略,造成盆地内地震动及地震危险性预测结果普遍被低估。本文以地震动观测记录数据充足的日本关东盆地为例,采用地震动残差分析方法评估盆地附加放大效应,分析覆盖层厚度、盆地内空间位置、震级、震源距对地震动放大效应的影响,建立关东盆地附加放大效应经验评估模型。分析表明:关东盆地附加放大效应与反应谱周期相关,整体上从短周期的1.0逐渐增大至周期为5s时的1.5,附加放大效应与覆盖层厚度相关性较小,主要受盆地空间位置和震源距的影响;盆地北部边缘及西北部地区附加放大效应更强烈,盆地南部附加放大效应较小,这可能与盆地边缘效应密切相关。本文建立的关东盆地附加放大效应经验模型略高于BSSA14和ASK14模型的放大效应预测。相关研究结果可用于我国地震动预测、下一代地震动区划图修订等。
  • 中国大陆东部断裂系统是由一系列北东和北西走向断裂构成的共轭剪切破裂系统,其中北西向断裂处于破裂初期,活动性较强,多断续,集中成带分布并一定程度上控制着中强地震的发生(邵云惠,1980李祖武,1983冯希杰,1988徐杰等,2003)。针对该组断裂的研究成果将对防震减灾产生重要意义。由于平原区NW向断裂都是隐伏断裂,所以对其具体位置尚需更深入的研究。通过测量断层气来探测隐伏断裂带已成为行之有效的方法之一,目前已在多条隐伏断裂探测工作中取得了良好的效果(汪成民等,1991陈刚等,1995苏鹤军等,2005刘菁华等,2006邵永新等,2007姚道平等,2008刘学领等,2011冯军等,2011张慧等,2013)。该方法将断裂视为天然地质体中连接深部与浅部的“通道”,大量深部气体如H2、Rn和Hg气等可经该通道向上逸出,导致以上气体浓度在断裂位置附近异常升高,利用该原理则可通过圈定断层气异常点位来探测断裂的空间位置。

    研究区位于华北平原西部,区内第四系厚度达150—300m以上(河南省地质矿产局,1989)。新郑-太康断裂为一条北西向隐伏断裂(王志铄等,2017),由于缺乏系统的浅层地震勘探和钻孔资料,目前对其空间展布仍然存在争议。本次在利用石油勘探资料对新郑-太康断裂进行深部定位的基础上,在新郑地区布置2条与断层走向近垂直的地球化学测线,通过测量土壤中的H2浓度并分析H2分布特征,以确定该断裂的浅层位置。

    太康隆起西接豫西隆起,东邻鲁西隆起,南北分别为周口坳陷和开封坳陷。该区是豫西隆起与鲁西隆起延伸、连接的过渡单元,其基底以元古代五指岭组花岗片麻岩、黑云母片麻岩为主,上覆震旦系坚硬的安山玢岩(河南石油勘探指挥部地质大队,1976),进入新生代以来该区整体下沉并伴随强烈的断裂活动(石油化工部物探局,1975)。南部周口坳陷是第三系及上古生界的沉积坳陷,基底埋深约3—6km,岩性以五指岭组片岩、嵩山组石英岩和登封杂岩为主(燃化部六四六厂研究所,1972)。其内部发育多条具有一定规模的北西向隐伏断裂(图 1),如登封-西华断裂(F4)和临颍-郸城断裂(F5),将周口坳陷切割为多个次级凹陷与凸起。

    图 1  区域地震、地震构造与测线分布图(图(b)为图(a)红框区的放大图)
    Figure 1.  Regional seismo-tectonic map and location of survey lines

    新郑-太康断裂(F1)起于新郑市北,向东经鄢陵,延伸至太康县以南,长147km,是太康隆起和周口坳陷的分界断裂(图 1)。最新的深部石油勘探资料表明该断裂整体倾向北,倾角较陡,在东段显示出正断特征,中段和西段表现为逆断层。新郑-太康断裂在浅部向上分叉,与北部次级断裂组合为负花状构造,具有明显的离散走滑特征,且新近纪以来活动明显,新近系底界断距可达150m(图 2,据河南油田勘探开发研究院内部资料)。

    图 2  新郑-太康断裂东段TK-90-311地震偏移剖面(剖面位置见图 1(a)
    Figure 2.  Seismic profile No.TK-90-311 across the eastern segment of Xinzheng-Taikang fault

    此外,周口坳陷内发生过多次5.0—6.0级地震,地震活动水平略高于太康隆起,区内最新的地震为2010年太康4.7级地震,震中位置见图 1(a),发震构造为新郑-太康断裂(孙杰等,2014)。

    本次H2测试采用的仪器为德国德尔格公司生产的X-am 7000型多气体检测仪,传感器为电化学传感器,当气体扩散进入传感器后,在反应电极表面经过氧化或还原反应产生电流,并通过外电路流经2个电极,可通过外电路的负荷电阻测量气体浓度。该传感器可检测H2浓度上限为2000ppm,灵敏度为5ppm,测量误差±5%。该仪器稳定性和抗干扰能力均较强。

    测线位置选择在新郑市解放路路侧和马庄村-前宫村的村间道路旁,测线位置见图 1(b),附近无污染源,地下水位稳定,埋深约10—15m。测量前,先对仪器进行校正,稳定后开始测量。测量时,先用锤将1.2m的钢钎打入土中成孔,拔出后立即将锥形麻花钻收集管拧入土中封孔,土壤气体通过锥形麻花钻收集管、胶皮管后进入测试仪器,由于H2比重小,逸出迅速。测量多次取峰值作为最终结果。

    异常值的判定采用数值分析法,考虑到:① 测试区附近无污染源,空气中H2浓度背景值较低(2.113—5.998ppm);② 松散的第四系覆盖较厚,H2逸出条件复杂多变;③ 该断裂东段(太康段)向上延伸并分叉形成花状构造,分支断裂数目多且规模不等(中国地震局地球物理勘探中心,2015),同时断裂两侧裂隙发育,气体的逸出通道彼此相连或将导致产生于次级断裂深部的H2沿逸出条件更好的主断裂逸出,即气体发生偏移,我们认为异常值下限定得过高会导致部分次级断裂信息的丢失,因此将平均值与1倍标准差之和、平均值与2倍标准差之和综合考虑来划分异常。

    断层的活动性质、几何特征(Ciotoli等,2007)以及破碎带和断层的发育程度(King等,1996Annunziatellis等,2008)在一定程度上决定了地下裂隙的产生、分布及闭合,进而通过改变地下通道的连通状况来影响气体由深至浅的迁移过程。断层气浓度曲线蕴含丰富的地质信息,以下将对2条测线分别进行讨论:

    (1)解放路测线

    该测线位于新郑市西北解放路路侧(图 1(b)),共39个测点,长2035m,平均测点间距50m,当H2浓度明显高时,加密到每25m一个测点。至测线两端,测点间距扩大至100m。

    各测点H2浓度的平均值μ=214.6ppm,标准差σ=341.4ppm,我们将浓度高于μ+2σ的异常点定为A类异常,表示异常幅度大,异常形态明显,对应规模较大的断裂;将高于μ+σ的异常点定为B类异常,该类异常幅度稍小,对应规模稍小的次级断裂。该测线异常点位置及类型见表 1,解释后的断裂见图 3。H2浓度在310m处异常高(1901ppm),为A类异常,且明显高于两侧,反映气体的逸出条件向两侧逐渐变差,破碎带集中在一个较窄的范围,故判断该处发育近直立的走滑断裂(Ciotoli等,2007);在760m与1035m处H2浓度接近(分别为789ppm和814.1ppm),为B类异常,在气体浓度曲线上表现出“双峰”态。由于新郑-太康断裂长期活动,带内充分发育断层泥等细粒物质导致断裂带渗透性降低,气体绕过渗透性低的断层带而从其两侧逸出,造成860—1010m间H2浓度低(340—360ppm)的现象,该处对应断裂带(King等,1996Annunziatellis等,2008),其宽约150m。H2浓度在1035m处向北迅速下降,判断该断裂带倾向北。根据异常类型判断,该断裂带的规模、活动性不及南部的走滑断裂(张大其,1991孟广魁等,1997焦德成等,2012)。

    表 1  解放路测线异常点位及类型
    Table 1.  Location and type of anomalies from the Jiefang road survey line
    距测线起始点位置/m H2浓度/ppm 异常类型 断层倾向
    310 1901 A 近直立
    760 789 B N
    1035 814.1 B N
    下载: 导出CSV 
    | 显示表格
    图 3  解放路测线H2测量结果与异常分析
    Figure 3.  The H2 concentration curve of the Jiefang road survey line
    (a)测线异常点位与石油勘探资料解译的深部断点位置对比结果;(b)测线H2浓度曲线及解释

    平面上,将各异常点位置与断裂深部资料相比对可以发现,北部断裂带位置与断裂深部位置吻合较好。该测线H2浓度的高异常准确地反映了断裂的浅层位置(图 3)。

    (2)马庄村-前宫村测线

    该测线位于解放路测线的东南(图 1(b)),测线长2125m,共35个测点,点间距在测线两端为100m,高值附近加密到每25m一个点。各测点间H2浓度平均值μ为256.8ppm,标准差σ为302.7ppm。测线异常点位及类型见表 2,解释后的断裂见图 4

    表 2  马庄村-前宫村测线异常点位及类型
    Table 2.  Location and types of anomalies from the Mazhuang-qiangong survey line
    距测线起始点距离/m H2浓度/ppm 异常类型 断层倾向
    900 600 B S
    1225 1250 A S
    1475 1210 A S
    下载: 导出CSV 
    | 显示表格
    图 4  马庄村-前宫村测线H2测量结果与异常分析
    Figure 4.  The H2 concentration curve of the Mazhuang-qiangong survey line
    (a)测线异常点位与深部石油勘探资料解译的断点位置对比结果;(b)测线H2浓度曲线及解释

    根据表 2,该测线共有2处A类异常:H2浓度在测线1225m和1475m处分别为1250ppm和1210ppm,在浓度曲线上表现出与解放路类似的“双峰”形态,1275—1425m之间H2浓度低(25—175ppm),对应断裂带宽约150m。自1175m处向南H2浓度逐渐下降,故判断该断裂带倾向南。900m附近存在一处B类异常,H2浓度为600ppm,对应活动性稍弱的次级断裂(图 4)。

    断裂在该段的深部位置由TK-94-237.5测线和TK-94-244.5测线控制,可以看到H2浓度异常点位置与石油勘探资料确定的断裂深部位置吻合较好,前者较准确地反映了断裂的浅层位置。

    对比解放路和马庄村-前宫村H2浓度曲线可知,2条H2浓度曲线解释出了相近的断裂组合,即北部存在宽约150m的断裂带,以及向南约325—450m处发育的单支断裂。不同之处在于断裂带的倾向:解放路H2浓度曲线显示断裂带北倾,而马庄村-前宫村H2浓度曲线显示断裂带倾向南,倾向的变化表明新郑-太康断裂具有明显的走滑特征(漆家福等,2003宋明春等,2015)。

    不仅如此,我们还发现如下现象:

    就单一测线而言,解放路H2浓度曲线表明该处南部的断裂规模更大、活动性更强,而马庄村-前宫村H2浓度曲线则表明该处北部断裂带的规模更大、活动性更强;就H2浓度沿断裂走向的变化规律而言,由北西向南东,北部断裂带的H2浓度升高,而南部断裂的H2浓度降低。

    通过圈定和分析解放路与马庄村-前宫村测线上H2浓度的异常点,与深部石油勘探资料相对比,获得的主要结论如下:

    (1)2条测线出现的高水平H2点位与石油勘探资料确定的断裂深部位置吻合度较高,准确地反映了断裂的浅层位置。

    (2)2条测线同步解释出一条走滑断裂,以及其北部发育的一条宽约150m、产状变化的断裂带,产状变化与其走滑性质有关。

    值得注意的是,就气体浓度沿断裂走向的变化规律而言,由北西向南东,北部断裂带H2浓度升高,而南部断裂的H2浓度降低。在气体逸出条件相近的前提下,南、北2条断裂的活动性呈“此消彼长”的关系。考虑到2条测线相距约20km,我们认为这2条测线可能位于一组左旋左阶断裂的阶区部位。但由于受到测线密度的限制,该结论有待进一步证实。

    此外,关于该区氢气来源的问题亦值得讨论。华北地台具有由结晶基底和沉积盖层组成的“双重结构”(翟明国,2011朱日祥等,2012滕吉文等,2014):中条运动形成的结晶基底及上覆未经变质的中元古代以来的沉积盖层,盖层以碳酸盐岩和碎屑岩为主,而这些岩层中产生的H2微不足道(Ware等,1984)。目前已有大量实验表明,H2与CO2和CH4不同,其主要为无机过程的产物(Giardini等,1976Ruff等,1980Sugisaki等,1983)。构造应力作用下,破碎的矿物如石英(Hughes-Schrader等,1961Wakita等,1980)、长石及硅酸盐矿物(Sugisaki,1980)表面的原子团与渗入裂隙中的水反应可以生成H2(Scharder机制,即Si+H2O→Si-OH+H2),含有黑云母的花岗岩可释放更多的H2Sugisaki等,1983)。前已提及,太康隆起的沉积盖层下为震旦系火成岩、以花岗片麻岩及黑云母片麻岩为主的结晶基底,富含石英、黑云母等矿物,在动力破碎后通过反应释放出的大量H2,可能是本区土壤氢气的主要来源,由此可初步判断新郑-太康断裂至少为发育于上地壳的断裂,但由于缺乏更广泛、更精确的断层气及钻孔资料,该结论有待进一步研究和证实。

  • 图  1  选取的地震事件及关东盆地内K-NET台站分布

    Figure  1.  Distribution of selected seismic events and K-NET stations in the Kanto basin

    图  2  震级-震源距分布及其0.3s、3s的水平向PSA

    Figure  2.  Magnitude vs. hypocentral distance and observed horizontal PSA at period of 0.3s and 3s for recordings under consideration

    图  3  覆盖层厚度对盆地附加放大效应的影响

    Figure  3.  The influence of sedimentary thickness on extra amplification effect

    图  4  盆地内K-NET台站的附加放大效应空间分布

    Figure  4.  Spatial distributions of extra amplification effects at various spectral periods in the Kanto basin

    图  5  震级对盆地附加放大效应的影响

    Figure  5.  The influence of magnitude on extra amplification effect

    图  6  震源距对盆地附加放大效应的影响

    Figure  6.  The influence of hypocentral distance on extra amplification effect

    图  7  基础经验模型拟合曲线及盆地附加放大效应经验模型

    Figure  7.  The fitness curve of basic empirical model and the empirical model of extra amplification effect in Kanto basin

    图  8  经验模型与BSSA14、ASK14模型盆地项对比

    Figure  8.  The comparisons of empirical model with model BSSA14 and model ASK14

    表  1  基础经验模型的拟合系数

    Table  1.   Fitness coefficients of empirical basic model

    系数
    C T1 T2 n1 n2
    - 0.08750 1.6209 7.3731 0.7381 1.7743
    下载: 导出CSV

    表  2  震源距调整基础模型的拟合系数

    Table  2.   Fitness coefficients of adjusted model scaling with hypocentral distance

    震源距/km 系数
    Di T3i n3i
    0<Rhyp≤100 0.3511 0.9702 -0.2742
    100<Rhyp≤200 -0.1264 1.5950 0.1258
    200<Rhyp≤300 -0.6790 1.3022 0.6643
    下载: 导出CSV
  • 刘启方, 于彦彦, 章旭斌, 2013. 施甸盆地三维地震动研究. 地震工程与工程振动, 33(4): 54-60. http://www.cnki.com.cn/Article/CJFDTotal-DGGC201304007.htm
    刘烁宇, 李英民, 2018. 盆地型长周期地震动的判别准则. 湖南大学学报(自然科学版), 45(5): 85-93. http://www.cnki.com.cn/Article/CJFDTotal-HNDX201805010.htm
    潘华, 张萌, 李金臣, 2017. 美国地震区划图的发展——地震危险性图与抗震设计图. 震灾防御技术, 12(3): 511-522. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20170307&journal_id=zzfyjs
    王海云, 2011. 渭河盆地中土层场地对地震动的放大作用. 地球物理学报, 54(1): 137-150. doi: 10.3969/j.issn.0001-5733.2011.01.015
    张萌, 潘华, 李金臣, 2017. 美国建筑抗震设计的法规体系与设计地震动的确定. 震灾防御技术, 12(2): 306-318. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20170207&journal_id=zzfyjs
    Abrahamson N. A., Silva W. J., Kamai R., 2014. Summary of the ASK14 ground motion relation for active crustal regions. Earthquake Spectra, 30(3): 1025-1055. doi: 10.1193/070913EQS198M
    Atik L. A., Abrahamson N., Bommer J. J., et al., 2010. The variability of ground-motion prediction models and its components. Seismological Research Letters, 81(5): 794-801. doi: 10.1785/gssrl.81.5.794
    Boore D. M., Stewart J. P., Seyhan E., et al., 2014. NGA-West2 equations for predicting PGA, PGV, and 5% damped PSA for shallow crustal earthquakes. Earthquake Spectra, 30(3): 1057-1085. doi: 10.1193/070113EQS184M
    Chandramohan R., Baker J. W., Deierlein G. G., 2016. Quantifying the influence of ground motion duration on structural collapse capacity using spectrally equivalent records. Earthquake Spectra, 32(2): 927-950. doi: 10.1193/122813eqs298mr2
    Denolle M. A., Miyake H., Nakagawa S., et al., 2014. Long-period seismic amplification in the Kanto Basin from the ambient seismic field. Geophysical Research Letters, 41(7): 2319-2325. doi: 10.1002/2014GL059425
    Eads L., Miranda E., Lignos D. G., 2015. Average spectral acceleration as an intensity measure for collapse risk assessment. Earthquake Engineering and Structural Dynamics, 44(12): 2057-2073. doi: 10.1002/eqe.2575
    Frankel A., Stephenson W., Carver D., 2009. Sedimentary basin effects in Seattle, Washington: ground-motion observations and 3D simulations. Bulletin of the Seismological Society of America, 99(3): 1579-1611. doi: 10.1785/0120080203
    Furumura T., Hayakawa T., 2007. Anomalous propagation of long-period ground motions recorded in Tokyo during the 23 October 2004 Mw 6.6 Niigata-ken Chuetsu, Japan, earthquake. Bulletin of the Seismological Society of America, 97(3): 863-880. doi: 10.1785/0120060166
    Hruby C. E., Beresnev I. A., 2003. Empirical corrections for basin effects in stochastic ground-motion prediction, based on the Los Angeles Basin analysis. Bulletin of the Seismological Society of America, 93(4): 1679-1690. doi: 10.1785/0120020121
    Koketsu K., Miyake H., Afnimar, et al., 2009. A proposal for a standard procedure of modeling 3-D velocity structures and its application to the Tokyo metropolitan area, Japan. Tectonophysics, 472(1-4): 290-300. doi: 10.1016/j.tecto.2008.05.037
    Marafi N. A., Eberhard M. O., Berman J. W., et al., 2017. Effects of deep basins on structural collapse during large subduction earthquakes. Earthquake Spectra, 33(3): 963-997. doi: 10.1193/071916eqs114m
    Meza-Fajardo K. C., Semblat J. F., Chaillat S., et al., 2016. Seismic-wave amplification in 3D alluvial basins: 3D/1D amplification ratios from fast multipole BEM simulations. Bulletin of the Seismological Society of America, 106(3): 1267-1281. doi: 10.1785/0120150159
    Moschetti M. P., Thompson E. M., Rekoske J., et al., 2020. Ground-motion amplification in Cook Inlet region, Alaska from intermediate-depth earthquakes, including the 2018 Mw 7.1 anchorage earthquake. Seismological Research Letters, 91(1): 142-152. doi: 10.1785/0220190179
    Skarlatoudis, A., Hosseini, M., Somerville, P. 2015. Basin Amplification Factors for Cascadia from the 2011 Tohoku Earthquake. 1-67.
    Wang H. Y., Wang S. Y., 2015. A new method for estimating VS(30) from a shallow shear-wave velocity profile (depth < 30m). Bulletin of the Seismological Society of America, 105(3): 1359-1370. doi: 10.1785/0120140103
    Wen R. Z., Xu P. B., Wang H. W., et al., 2018. Single-station standard deviation using strong-motion data from Sichuan Region, China. Bulletin of the Seismological Society of America, 108(4): 2237-2247. doi: 10.1785/0120170276
    Yuan, X. M., Liao, Z. P., 1995. Scattering of plane SH waves by a cylindrical alluvial valley of circular-acr cross-section. Earthquake Engineering and Structural Dynamics, 24(10): 1303-1313. doi: 10.1002/eqe.4290241002
    Zhao J. X., Zhang J., Asano A., et al., 2006. Attenuation relations of strong ground motion in Japan using site classification based on predominant period. Bulletin of the Seismological Society of America, 96(3): 898-913. doi: 10.1785/0120050122
  • 加载中
图(8) / 表(2)
计量
  • 文章访问数:  197
  • HTML全文浏览量:  64
  • PDF下载量:  19
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-21
  • 网络出版日期:  2021-04-07
  • 刊出日期:  2020-12-01

目录

/

返回文章
返回