Research on strong ground motion simulation and ShakeMap of 2014 Yutian MS7.3 earthquake
-
摘要: 基于地震动的时空衰减规律和传播特征,采用邻近地震监测台站地震动时程对1 km×1 km尺度的网格点进行近实时插值计算,同时结合场地效应对震区地震动参数进行修正,并以2014年2月12日新疆于田MS7.3地震为例,计算震区格网内各点的地震时程,同时以8 s为时间间隔绘制出地震动峰值等值线图并将其连续播放,得到了于田MS7.3地震峰值地震动(PGV,PGA)的空间分布.结果表明,于田县东部至民丰县北部地区受场地条件影响,震区震害在软弱地基土层及浅地下水位等因素作用下对震区地震动具有明显的放大效应,预测的地震动特征与现场宏观调查结果是一致的.在当前强震台网分布不均匀的情形下,本文方法能较好地描述震区地震动特征,较客观地反映灾区的强地面运动特征.Abstract: Based on attenuation law and space-time propagation characteristics of seismic ground motion, we use recorded time histories from adjacent monitoring stations to make near real-time interpolation calculation in a 1 km×1 km grid, and correct ground motion parameters according to the site effects. Taking the 12th February, 2014 Yutian MS7. 3 earthquake as an example, we calculated seismic time history in each grid near the epicenter, and then plot peak ground motion contour maps with 8 s intervals continuously. Finally we showed spatial distribution of peak ground motions (PGV and PGA) of this earthquake. The results show that, influenced by site conditions as well as soft ground soil and shallow groundwater level, ground motions in the region from eastern Yutian County to northern Minfeng County were amplified obviously, which were well consistent with field survey results. On the conditions of uneven distribution of strong ground motion observation stations, the method can be used to describe the characteristics of earthquake ground motions, and reflect the characteristics of strong ground motion in disaster areas objectively.
-
Keywords:
- Yutian MS7. 3 earthquake /
- seismic ground motion /
- monitoring station /
- ShakeMap
-
引言
场地卓越周期是描述场地特性的重要指标. 地表土层对不同周期的地震波有选择放大作用,致使在地震记录图上显示某些周期的波形特别多且好,即显得“卓越”,故称为地震卓越周期(《岩土工程手册》编写委员会,1994). 确定场地卓越周期的方法有以下3种(童广才,刘康和,2000): ① 当场地内有强震动记录时,通过频谱分析确定; ② 由常时微动测试分析确定; ③ 根据场地分层剪切波速测试结果按其子层周期求和公式计算. 前两种方法均是通过确定场地的地震动卓越周期以进一步确定场地的卓越周期. 高广运等(2000)和陈鹏等(2009)将卓越周期分为记录卓越周期Tr、 测试卓越周期Tm和波速卓越周期Tv. Tr为真实反应地震动的卓越周期(利用地震动卓越周期间接反映场地卓越周期,该周期也在不同程度上受到地震震源谱特性和区域地壳介质对地震波传播的影响),当工程场地范围内有适宜的强震动记录时,抗震设计应首选Tr; Tm为接近场地固有周期的卓越周期; Tv为与场地固有周期相比有一定误差的卓越周期,除非地基土层基本满足均匀平行的条件,否则应尽量避免选用.
场地卓越周期在数值上与卓越频率互为倒数,因此卓越周期可由卓越频率求倒得出. 目前日本强震动观测台网KiK-net记录中既有场地地下记录(NS1,EW1,UD1),也有地表记录(NS2,EW2,UD2),因此对场地内强震记录通过频谱分析确定其卓越频率时可以采用以下3种方法: ① 地表记录的傅里叶谱分析法; ② 地表水平/垂直傅里叶谱比法(Wen et al,2006; 任叶飞等,2013); ③ 地表/地下傅里叶谱比法.
本文拟选取日本强震动观测台网中两个基岩台站和两个Ⅲ类场地台站的数百条记录,分别采用上述3种方法对两个基岩台站和两个Ⅲ类场地台站的卓越频率进行分析,并对比分析各种方法的优缺点及其适用情况.
1. 强震动数据
本文从日本强震动观测台网中选取了FKOH06和KGSH03两个基岩台站以及TCGH16和IBRH17两个Ⅲ类场地台站,这些台站分别有地表和相应的地下观测点,各台站分别获得了大量的地震事件记录. 在选择台站地震动记录时,对地震事件随机选取,但是需注意所选的地震动记录峰值加速度的分布应尽量能覆盖各个等级. 最终用于本文研究的强震动记录为: FKOH06台站29次地震事件的地表和地下记录,共174条; KGSH03台站25次地震事件的地表和地下条记录,共150条; TCGH16台站45次地震事件的地表和地下记录,共270条; IBRH17台站33次地震事件的地表和地下记录,共198条.
2. 数据分析
首先使用ViewWave软件读取地震记录,该软件默认对记录进行滤波和基线校正; 其次进行傅里叶变换,分别绘制出各条记录的傅里叶谱以及地表/地下和地表水平分量/垂直分量的傅里叶谱比曲线,读取傅里叶谱最大幅值以及傅里叶谱比曲线中谱比最大值所对应的频率值,即为卓越频率; 然后分别绘制出各场地经上述3种方法分析得到的卓越频率的散点分布图; 最后对比卓越频率的分布情况,分析在确定场地卓越频率时各方法的适用性.
图 1—6给出了FKOH06和KGSH03两个基岩台站所在场地的卓越频率分布. 由图 1和图 4可以看出,基岩台站地下和地表记录的3个分量的傅里叶谱卓越频率分布均很分散. 由图 2和图 5可以看出,台站所在场地的地表水平/垂直傅里叶谱比卓越频率分布也都很分散. 由图 3和图 6可以看出,场地的水平分量和垂直分量的地表/地下傅里叶谱比卓越频率分布不同: 对于FKOH06台站所在场地,水平分量南北向和东西向卓越频率分布比较分散,而垂直分量卓越频率分布相对集中; 对于KGSH03台站所在场地则得到相反的结论. 综上,对于基岩场地卓越频率,上述3种方法所得结果均很分散,无法给出一个准确的卓越频率值,这表明基岩地震动的卓越周期主要反映的是地震震源谱特性和区域地壳介质对地震波传播的影响,而局部基岩场地对地震动卓越周期的影响则不明显.
图 7—12为Ⅲ类场地台站的卓越频率分布图. 由图 7和图 10可以看出: Ⅲ类场地台站TCGH16和IBRH17地下记录的3个分量的傅里叶谱卓越频率分布比较分散; TCGH16台站地表水平分量的傅里叶谱卓越频率分布比较集中,约为4.47 Hz,而其垂直分量的傅里叶谱卓越频率分布比较分散; IBRH17台站地表水平分量和垂直分量的傅里叶谱卓越频率分布均比较集中,分别约为8.96 Hz和10.26 Hz. 由图 8和图 11可以看出,TCGH16台站所在场地的地表水平/垂直傅里叶谱比卓越频率分布比较集中,约为4.53 Hz,而IBRH17台站所在场地则比较分散. 由图 9和图 12可以看出,TCGH16和IBRH17台站所在场地的地表/地下傅里叶谱比卓越频率分布均比较集中,其所在场地的水平分量的地表/地下傅里叶谱比卓越频率分别约为4.54 Hz和9.35 Hz,相应垂直分量约为11.29 Hz和11.39 Hz. 由上述结果可知: 对于TCGH16台站所在场地,3种方法均能给出一个相对确定的卓越频率值,分别为4.47 Hz,4.53 Hz和4.54 Hz; 对于IBRH17台站所在场地,地表水平/垂直傅里叶谱比卓越频率分布比较分散,而由地表记录傅里叶变换所得到的卓越频率为8.96 Hz,由地表/地下傅里叶谱比法得到的卓越频率为9.35 Hz.
3. 讨论与结论
对于基岩台站FKOH06和KGSH03所在场地,由3种方法所得的卓越频率分布均比较分散,无法给出一个确定的数值,初步认为是由于场地岩层坚硬、 局部场地条件对地震动影响较小所致,因此场地地表和地下记录所体现的主要还是地震震源谱特性和区域地壳介质对地震波传播的影响,而场地上不同地震的震源及地震波传播路径存在较大差异,从而导致所得的场地卓越频率分布比较分散.
对于Ⅲ类场地,TCGH16和IBRH17台站所在场地的土层条件不同,上述3种方法在具体场地的适用性也不同. 对于TCGH16台站所在场地,3种方法均能给出一个确定的卓越频率数值,且3种方法所得数值很接近. 对于IBRH17台站的所在场地,地表水平/垂直傅里叶谱比卓越频率分布比较分散,无法给出确定数值,而另外两种方法均可给出确定数值,但数值相差较大. 结合两个台站所在场地的具体土层条件分析,TCGH16和IBRH17台站所在场地土层均为砂层、 砂砾夹有黏土层,但等效剪切波速(前者172 m/s左右,后者272 m/s)及覆盖土层厚度(前者约75 m,而后者深达235 m)差别较大,初步认为是地表土层剪切波速及土层厚度的影响所导致的.
综合上述3种方法的分析结果,可以得到如下结论:
1)地表记录傅里叶谱分析法所得的卓越频率所包含信息较多,既有地震震源和区域地壳介质的影响,也有局部场地特性影响. 对于基岩场地,由于场地条件对地震动影响较小,而地震震源和区域地壳介质的影响比较突出,因此同一场地上由不同记录得出的结果偏差较大,如上述两个基岩台站所在场地. 但对于地表有土层的场地,场地条件对地震动影响明显,可以得出较为准确的结果.
2)地表水平/垂直傅里叶谱比法的前提是认为场地条件对地震动垂直分量的影响相对于水平分量要小得多,但当场地地表为较厚土层时,场地条件对地震动垂直分量的影响也较大,因此该方法在某些场地适用(如TCGH16台站所在场地),但在某些场地所得结果却很分散(如IBRH17台站所在场地).
3)地表/地下傅里叶谱比法理论上最能反映场地条件对地震动的影响,本文认为该方法所得场地卓越频率是比较准确的. 作者还选取了20次不同地震事件在TCGH16台站的地震动记录与相邻的TCGH12台站的地震动记录进行了对比分析,结果显示两个台站地下记录的傅里叶谱差异很小(图 13),但在地表差异却较大(图 14),这说明地表/地下傅里叶谱比法确实能很好地体现场地条件对地震动的影响,从而验证了该方法的准确性.
对于基岩场地,由于局部场地条件的影响相对较弱,地震动频谱中包含的场地特性影 响信息很少,所以3种方法所得结果均很分散; 对于地表存在土层的场地,地表水平/垂直傅里叶谱比法所得结果存在不确定性,而地表/地下傅里叶谱比法则能给出一个相对较为准确的场地卓越频率值.
-
图 1 插值计算方法(a)及插值原理(b)示意图(张红才,2008)
Figure 1. Schematic diagram of interpolation calcutation method (a) and interpolation principle (b) (after Zhang, 2008)
表 1 场地放大系数
Table 1 Site amplification factor
场地条件 Ⅱ类场地PGA/(cm·s-2) 50 100 150 200 300 400 Ⅰ0 0.72 0.74 0.75 0.76 0.85 0.90 Ⅰ1 0.80 0.82 0.83 0.85 0.95 1.00 Ⅱ 1.00 1.00 1.00 1.00 1.00 1.00 Ⅲ 1.30 1.25 1.15 1.00 1.00 1.00 Ⅳ 1.25 1.20 1.10 1.00 0.95 0.90 表 2 本文所用台站及其主要参数
Table 2 Stations and main parameters used in this paper
台站名称 震中距/km 岩性 主要仪器类型 于田 61.1 灰绿岩 CTS-1E, EDAS-24IP 应急B 115.2 变质岩 CMG-40T, EDAS-24IP 和田台阵0 214.6 砂岩 CTS-1E, EDAS-24IP 叶城 432.3 石英砾岩 CMG-3ESPC, EDAS-24IP 且末 267.4 灰白色花岗岩 CTS-1E, EDAS-24IP 若羌 423.4 变质泥岩 CTS-1E, EDAS-24IP 狮泉河 419.6 变质泥岩 CTS-1E, EDAS-24IP 改则 432.8 变质泥岩 CTS-1E, EDAS-24IP 阿拉尔 483.4 混合花岗岩 FSS-3BDH, EDAS-24IP 巴楚 493.4 变质岩 CTS-1E, EDAS-24L6 西克尔 499.2 灰岩 BBVS-60, EDAS-24IP 乌什 513.6 变质岩 CTS-1E, EDAS-24IP 阿克苏 502.8 变质岩 CTS-1E, EDAS-24IP 神木园 519.4 变质岩 CTS-1E, EDAS-24IP -
陈鲲, 俞言祥, 高孟潭. 2010.考虑场地效应的ShakeMap系统研究[J].中国地震, 26 (1): 92-102. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD201001010.htm Chen K, Yu Y X, Gao M T. 2010. Research on ShakeMap system in terms of the site effect[J]. Earthquake Research in China, 26 (1): 92-102 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD201001010.htm
霍俊荣, 胡聿贤. 1992.地震动峰值参数衰减规律的研究[J].地震工程与工程振动, 12 (2): 1-11. http://www.cnki.com.cn/Article/CJFDTOTAL-DGGC199202000.htm Huo J R, Hu Y X. 1992. Study on attenuation laws of ground motion parameters[J]. Earthquake Engineering and Engineering Vibration, 12 (2): 1-11 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DGGC199202000.htm
李铂, 董翔, 陈亚红, 何琳, 于澄. 2012.线源地震影响场计算模型[J].地震地磁观测与研究, 33 (3): 14-19. http://www.cnki.com.cn/Article/CJFDTOTAL-DZGJ2012Z1006.htm Li B, Dong X, Chen Y H, He L, Yu C. 2012. Realization of calculation model of line-source effect fields[J]. Seismic Geomagnetic Observation and Research, 33 (3): 14-19 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DZGJ2012Z1006.htm
李海兵, 孙知明, 潘家伟, 刘栋梁, 张佳佳, 李成龙, 刘亢, 云锟, 龚正. 2014. 2014年于田MS7.3地震野外调查:特殊的地表破裂带[J].地球学报, 35 (3): 391-394. http://www.oalib.com/paper/4868084 Li H B, Sun Z M, Pan J W, Liu D L, Zhang J J, Li C L, Liu K, Yun K, Gong Z. 2014. Field study of the 12 February 2014 Yutian MS7.3 earthquake: A special surface rupture zone[J]. Acta Geoscientica Sinica, 35 (3): 391-394 (in Chinese). http://www.oalib.com/paper/4868084
刘军, 宋立军, 胡伟华, 李帅. 2014.新疆地区场地条件分区及其矢量化[J].内陆地震, 28 (1): 70-74. http://www.cnki.com.cn/Article/CJFDTOTAL-LLDZ201401010.htm Liu J, Song L J, Hu W H, Li S. 2014. Zoning and vectorize of site condition in Xinjiang[J]. Inland Earthquake, 28 (1): 70-74 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-LLDZ201401010.htm
金星, 张红才, 韦永祥, 李军. 2010.基于地震监测台网资料近实时插值计算震动图的初步研究[J].防灾减灾学报, 26 (1): 1-11. http://www.cnki.com.cn/Article/CJFDTOTAL-DDYJ201001002.htm Jin X, Zhang H C, Wei Y X, Li J. 2010. Preliminary study of near real-time interpolation calculation of ShakeMap based on seismic monitoring network's information[J]. Journal of Disaster Prevention and Reduction, 26 (1): 1-11 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DDYJ201001002.htm
孟令媛, 周龙泉, 刘杰. 2014. 2014年新疆于田MS7.3地震近断层强地面运动模拟及烈度分布估计[J].地震学报, 36 (3): 362-371. http://www.dzxb.org/Magazine/Show?id=28966 Meng L Y, Zhou L Q, Liu J. 2014. Estimation of near-fault strong ground motion and intensity distribution of the 2014 Yutian, Xinjiang MS7.3 earthquake[J]. Acta Seismologica Sinica, 36 (3): 362-371 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28966
姚远, 伊力亚尔·阿不力孜, 吴传勇, 吴国栋, 陈建波. 2014. 2014年2月12日新疆于田MS7.3地震发震构造初步研究[J].内陆地震, 28 (2): 121-126. Yao Y, Yiliyar·Abulizi, Wu C Y, Wu G D, Chen J B. 2014. Preliminary study on seismogenic structure of Yutian MS7.3 earthquake on Feb.12, 2014[J]. Inland Earthquake,28 (2):121-126 (in Chinese).
张红才. 2008.基于地震监测台网资料的震动图及震动烈度研究[D].哈尔滨:中国地震局工程力学研究所: 17-23. Zhang H C. 2008. Research on ShakeMaps and Shaking Intensity Based on Records of Monitoring Network[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration: 17-23 (in Chinese).
张勇, 许力生, 陈运泰, 汪荣江. 2014. 2014年2月12日于田MW6.9地震破裂过程初步反演:兼论震源机制对地震破裂过程反演的影响[J].地震学报, 36 (2): 159-164. http://www.dzxb.org/Magazine/Show?id=28961 Zhang Y, Xu L S, Chen Y T, Wang R J. 2014. Fast inversion for the rupture process of the 12 February 2014 Yutian MW6.9 earthquake: Discussion on the impacts of focal mechanisms on rupture process inversions[J]. Acta Seismologica Sinica, 36 (2): 159-164 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28961
张振斌, 唐丽华. 2010.新疆南天山地区土层场地地震动峰值加速度衰减关系[J].中国地震, 26 (4): 401-406. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD201004007.htm Zhang Z B, Tang L H. 2010. The attenuation relationship of peak ground motion acceleration on soil site at the south Tianshan region in Xinjiang[J]. Earthquake Research in China, 26 (4): 401-406 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD201004007.htm
张振国, 张伟, 孙耀充, 朱耿尚, 文健, 陈晓非. 2014. 2014年2月12日新疆于田地震强地面运动初步模拟及烈度预测[J].地球物理学报, 57 (2): 685-689. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201402032.htm Zhang Z G, Zhang W, Sun Y C, Zhu G S, Wen J, Chen X F. 2014. Preliminary simulation pf strong ground motion for Yutian, Xinjiang earthquake of 12 February 2014, and hazard implication [J]. Chinese Journal of Geophysics, 57 (2): 685-689 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201402032.htm
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会. 2008.中国地震烈度表(GB/T17742-2008)[S].北京:中国标准出版社:5-7. General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, Standardization Administration of the People's Republic of China. 2008. The Chinese Seismic Intensity Scale (GB/T 17742-2008)[S].Beijing: Standards Press of China: 5-7.
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会. 2015.中国地震动参数区划图(GB/T18306-2015)[S].北京:中国标准出版社: 6-10. General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, Standardization Administration of the People's Republic of China. 2015. China Seismic Parameter Zoning Map (GB/T 18306-2015)[S]. Beijing: Standards Press of China: 6-10.
Borcherdt R D. 1994. Estimates of site-dependent response spectra for design (methodology and justification)[J]. Earthquake Spectra, 10 (4): 617-653. doi: 10.1193/1.1585791
-
期刊类型引用(10)
1. 陈钰鑫,李平,高志寅,田兆阳,薄景山,李孝波. 基于地脉动数据的北川老县城河谷场地效应研究. 振动与冲击. 2025(06): 254-262 . 百度学术
2. 吴雄伟,杨江,夏界宁,范涛. 红沿河核电厂地震仪表系统震后数据分析. 大地测量与地球动力学. 2024(09): 985-990 . 百度学术
3. 李小军,钱玉,荣棉水,孔小山. 基于地震动反应谱与傅里叶谱确定场地卓越频率的差异性. 震灾防御技术. 2024(03): 421-435 . 百度学术
4. 李孝波,宋霖君,宣雨童,吴义文,欧阳刚垒. 基于HVSR和VRSR法的场地地震响应分析. 振动与冲击. 2023(22): 303-311 . 百度学术
5. 兰景岩,宋锡俊,王婷. 上覆海水层对自由场基本周期影响的离心模型试验研究. 岩土工程学报. 2021(04): 768-775 . 百度学术
6. 潘章容,李同林,崔建文,田秀丰,石文兵,张卫东,袁洁. 2021年5月21日云南漾濞M6.4及相关地震强震动记录特征分析. 地震工程学报. 2021(04): 791-798 . 百度学术
7. 黄俊阁,罗永红,王运生,朱兴貌,欧剑锋,张跃跃,南凯. 基于环境噪声测试表征斜坡地震动响应:以自贡西山斜坡为例. 成都理工大学学报(自然科学版). 2020(01): 16-27 . 百度学术
8. 华永超,齐文浩,薄景山,郑桐. 场地周期的确定方法综述. 地震工程与工程振动. 2020(02): 239-251 . 百度学术
9. 戴苗,王秋良,刘坚,黄俊,魏贵春,申学林. 2017-06三峡库区M4.3、M4.1双地震强震记录分析. 大地测量与地球动力学. 2018(01): 18-23 . 百度学术
10. 杜修力,袁雪纯,黄景琦,许紫刚. 典型土层场地随机地震反应规律分析. 震灾防御技术. 2017(03): 574-588 . 百度学术
其他类型引用(15)