Basin Extra Amplification Effects from Seismic Ground-Motion Residual Analysis: A Case Study of Kanto Basin, Japan
-
摘要: 我国地震动预测及地震危险性分析通常仅考虑局部场地浅层岩土层对地震动的放大效应,不能考虑较大范围的地质条件影响,如沉积盆地厚沉积层对地震动的附加放大效应通常被忽略,造成盆地内地震动及地震危险性预测结果普遍被低估。本文以地震动观测记录数据充足的日本关东盆地为例,采用地震动残差分析方法评估盆地附加放大效应,分析覆盖层厚度、盆地内空间位置、震级、震源距对地震动放大效应的影响,建立关东盆地附加放大效应经验评估模型。分析表明:关东盆地附加放大效应与反应谱周期相关,整体上从短周期的1.0逐渐增大至周期为5s时的1.5,附加放大效应与覆盖层厚度相关性较小,主要受盆地空间位置和震源距的影响;盆地北部边缘及西北部地区附加放大效应更强烈,盆地南部附加放大效应较小,这可能与盆地边缘效应密切相关。本文建立的关东盆地附加放大效应经验模型略高于BSSA14和ASK14模型的放大效应预测。相关研究结果可用于我国地震动预测、下一代地震动区划图修订等。Abstract: The ground motion prediction equations of China does not take the basin effects into account, resulting in the apparent underestimations for predicted ground motions and seismic hazards in the basin. The enhanced ground motions due to the thick sedimentary in basin amplifies the earthquake damages on engineering structures located in the basin. To evaluate the extra amplification effects in Kanto Basin, Japan, the residual analysis method of ground motion was applied to calculate the extra amplification effects in Kanto Basin. The factors, including potential influences of sediment thickness, spatial location in the basin, seismic magnitude and hypocentral distance were discussed. The result indicates that, the extra amplification effect is correlated with the spectral period, which gradually increases from 1.0 at short period to 1.5 at spectral period of 5s. The extra amplification effects are independent with the sediment thickness, while dependent on the hypocentral distance and spatial location. The north and northwest margins in the basin experienced much stronger amplification effects, being relatively smaller in the south Kanto Basin. Based on the extra amplifications at these stations, we developed the empirical model applied in the Kanto Basin. Compared with the BSSA14 and ASK14 models of NGA-West2, our model provided stronger predictions. This model could be used to revise Chinese ground motion prediction equations for the next generation of zoning map.
-
Key words:
- Ground motion /
- Residual analysis /
- Kanto basin /
- Basin extra amplification effect /
- Empirical mode
-
引言
中国大陆东部断裂系统是由一系列北东和北西走向断裂构成的共轭剪切破裂系统,其中北西向断裂处于破裂初期,活动性较强,多断续,集中成带分布并一定程度上控制着中强地震的发生(邵云惠,1980;李祖武,1983;冯希杰,1988;徐杰等,2003)。针对该组断裂的研究成果将对防震减灾产生重要意义。由于平原区NW向断裂都是隐伏断裂,所以对其具体位置尚需更深入的研究。通过测量断层气来探测隐伏断裂带已成为行之有效的方法之一,目前已在多条隐伏断裂探测工作中取得了良好的效果(汪成民等,1991;陈刚等,1995;苏鹤军等,2005;刘菁华等,2006;邵永新等,2007;姚道平等,2008;刘学领等,2011;冯军等,2011;张慧等,2013)。该方法将断裂视为天然地质体中连接深部与浅部的“通道”,大量深部气体如H2、Rn和Hg气等可经该通道向上逸出,导致以上气体浓度在断裂位置附近异常升高,利用该原理则可通过圈定断层气异常点位来探测断裂的空间位置。
研究区位于华北平原西部,区内第四系厚度达150—300m以上(河南省地质矿产局,1989)。新郑-太康断裂为一条北西向隐伏断裂(王志铄等,2017),由于缺乏系统的浅层地震勘探和钻孔资料,目前对其空间展布仍然存在争议。本次在利用石油勘探资料对新郑-太康断裂进行深部定位的基础上,在新郑地区布置2条与断层走向近垂直的地球化学测线,通过测量土壤中的H2浓度并分析H2分布特征,以确定该断裂的浅层位置。
1. 区域地震地质背景
太康隆起西接豫西隆起,东邻鲁西隆起,南北分别为周口坳陷和开封坳陷。该区是豫西隆起与鲁西隆起延伸、连接的过渡单元,其基底以元古代五指岭组花岗片麻岩、黑云母片麻岩为主,上覆震旦系坚硬的安山玢岩(河南石油勘探指挥部地质大队,1976),进入新生代以来该区整体下沉并伴随强烈的断裂活动(石油化工部物探局,1975)。南部周口坳陷是第三系及上古生界的沉积坳陷,基底埋深约3—6km,岩性以五指岭组片岩、嵩山组石英岩和登封杂岩为主(燃化部六四六厂研究所,1972)。其内部发育多条具有一定规模的北西向隐伏断裂(图 1),如登封-西华断裂(F4)和临颍-郸城断裂(F5),将周口坳陷切割为多个次级凹陷与凸起。
新郑-太康断裂(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)。
2. 探测实施与结果
2.1 探测仪器及测试方法
本次H2测试采用的仪器为德国德尔格公司生产的X-am 7000型多气体检测仪,传感器为电化学传感器,当气体扩散进入传感器后,在反应电极表面经过氧化或还原反应产生电流,并通过外电路流经2个电极,可通过外电路的负荷电阻测量气体浓度。该传感器可检测H2浓度上限为2000ppm,灵敏度为5ppm,测量误差±5%。该仪器稳定性和抗干扰能力均较强。
测线位置选择在新郑市解放路路侧和马庄村-前宫村的村间道路旁,测线位置见图 1(b),附近无污染源,地下水位稳定,埋深约10—15m。测量前,先对仪器进行校正,稳定后开始测量。测量时,先用锤将1.2m的钢钎打入土中成孔,拔出后立即将锥形麻花钻收集管拧入土中封孔,土壤气体通过锥形麻花钻收集管、胶皮管后进入测试仪器,由于H2比重小,逸出迅速。测量多次取峰值作为最终结果。
2.2 数据处理及断裂浅部定位
异常值的判定采用数值分析法,考虑到:① 测试区附近无污染源,空气中H2浓度背景值较低(2.113—5.998ppm);② 松散的第四系覆盖较厚,H2逸出条件复杂多变;③ 该断裂东段(太康段)向上延伸并分叉形成花状构造,分支断裂数目多且规模不等(中国地震局地球物理勘探中心,2015),同时断裂两侧裂隙发育,气体的逸出通道彼此相连或将导致产生于次级断裂深部的H2沿逸出条件更好的主断裂逸出,即气体发生偏移,我们认为异常值下限定得过高会导致部分次级断裂信息的丢失,因此将平均值与1倍标准差之和、平均值与2倍标准差之和综合考虑来划分异常。
断层的活动性质、几何特征(Ciotoli等,2007)以及破碎带和断层的发育程度(King等,1996;Annunziatellis等,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等,1996;Annunziatellis等,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 平面上,将各异常点位置与断裂深部资料相比对可以发现,北部断裂带位置与断裂深部位置吻合较好。该测线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 根据表 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浓度降低。
3. 结论与讨论
通过圈定和分析解放路与马庄村-前宫村测线上H2浓度的异常点,与深部石油勘探资料相对比,获得的主要结论如下:
(1)2条测线出现的高水平H2点位与石油勘探资料确定的断裂深部位置吻合度较高,准确地反映了断裂的浅层位置。
(2)2条测线同步解释出一条走滑断裂,以及其北部发育的一条宽约150m、产状变化的断裂带,产状变化与其走滑性质有关。
值得注意的是,就气体浓度沿断裂走向的变化规律而言,由北西向南东,北部断裂带H2浓度升高,而南部断裂的H2浓度降低。在气体逸出条件相近的前提下,南、北2条断裂的活动性呈“此消彼长”的关系。考虑到2条测线相距约20km,我们认为这2条测线可能位于一组左旋左阶断裂的阶区部位。但由于受到测线密度的限制,该结论有待进一步证实。
此外,关于该区氢气来源的问题亦值得讨论。华北地台具有由结晶基底和沉积盖层组成的“双重结构”(翟明国,2011;朱日祥等,2012;滕吉文等,2014):中条运动形成的结晶基底及上覆未经变质的中元古代以来的沉积盖层,盖层以碳酸盐岩和碎屑岩为主,而这些岩层中产生的H2微不足道(Ware等,1984)。目前已有大量实验表明,H2与CO2和CH4不同,其主要为无机过程的产物(Giardini等,1976;Ruff等,1980;Sugisaki等,1983)。构造应力作用下,破碎的矿物如石英(Hughes-Schrader等,1961;Wakita等,1980)、长石及硅酸盐矿物(Sugisaki,1980)表面的原子团与渗入裂隙中的水反应可以生成H2(Scharder机制,即Si+H2O→Si-OH+H2),含有黑云母的花岗岩可释放更多的H2(Sugisaki等,1983)。前已提及,太康隆起的沉积盖层下为震旦系火成岩、以花岗片麻岩及黑云母片麻岩为主的结晶基底,富含石英、黑云母等矿物,在动力破碎后通过反应释放出的大量H2,可能是本区土壤氢气的主要来源,由此可初步判断新郑-太康断裂至少为发育于上地壳的断裂,但由于缺乏更广泛、更精确的断层气及钻孔资料,该结论有待进一步研究和证实。
-
表 1 基础经验模型的拟合系数
Table 1. Fitness coefficients of empirical basic model
值 系数 C T1 T2 n1 n2 - 0.08750 1.6209 7.3731 0.7381 1.7743 表 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 -
刘启方, 于彦彦, 章旭斌, 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 -