Some Technical Notes on the Data Processing of the Spectral Ratio Based on the Strong-motion Records
-
摘要: 强震动记录H/V谱比法自提出以来，已广泛应用于地震工程各研究领域中。但对数据处理的2个关键环节——Taper预处理及傅氏谱平滑未有详细研究。因此，本文针对强震动记录H/V谱比法计算涉及的重要数据处理过程，对四川地区19个强震动台站在汶川地震余震中获取的642组强震动记录开展研究。研究实例表明：S波H/V谱比振幅在周期 < 1s时高于全时程，论证计算H/V谱比时截取S波窗口的必要性；S波窗口的截断会引起傅里叶振幅谱的边瓣效应，若不加以处理，将显著影响低频部分H/V谱比结果，Taper预处理对于消除这种截断误差具有良好效果；在兼顾平滑效果及卓越周期识别准确度的基础上，通过不同带宽的Parzen窗试算，认为0.5Hz带宽宜用于傅里叶H/V谱比曲线平滑。Abstract: The H/V spectral ratio method based on strong motion record has been widely used in the field of earthquake engineering since it was developed; however, there are few research papers about the key steps of tapering the truncated S-wave time series and smoothing the Fourier amplitude spectrum in the strong-motion records. Based on 642 strong-motion records collected at 19 stations during Wenchuan aftershocks in Sichuan area, we illustrated some technical notes about data processing. The results show that the amplitude of the S-wave spectral ratio less than 1.0s period is higher than those calculated using the whole time series, and extracting S-wave window in calculating H/V spectral ratio is needed. The truncation of S-wave window could lead to the side lobe effect of Fourier amplitude spectrum, and the H/V spectral ratio at the low frequency part will be significantly contaminated if without tapering, which could effectively eliminate truncation error. On the basis of taking into account both the smoothing effect and the excellent accuracy of period identification, through the trial calculation of Parzen windows with different bandwidths, we consider that the 0.5Hz bandwidth is suitable for smoothing the Fourier H/V spectral ratio curve..
-
Key words:
- Strong-motion record /
- Spectral ratio method /
- Data processing /
- S-wave window
-
图 2 震级-震中距和震级-PGA分布
Figure 2. Distribution of magnitude and epicentral distance and peak ground acceleration
图 3 强震动记录S波窗口拾取示例
Figure 3. Illustration of selecting S wave window in strong motion recordings
图 4 强震动记录全时程、S波H/V谱比对比
Figure 4. H/V spectral ratio of whole-record and S wave window and the ratio of FASR
图 6 Taper处理前后的时程对比
Figure 6. The difference of S wave time recording between tapering(b) and without tapering(a)
图 7 Taper处理前后的傅氏谱对比
Figure 7. The difference of S wave window Fourier spectrum between tapering and without tapering
图 8 Taper处理对H/V谱比计算的影响
Figure 8. The influence of calculating H/V spectral ratio with tapering S wave window
图 9 不同带宽Parzen窗平滑后的S波H/V谱比曲线对比
Figure 9. The distinctness of H/V spectral ratio caves of S wave window smoothed with different band Parzen window
图 10 S波H/V平均谱比曲线在卓越周期处的标准差及所采用的强震动记录数量
Figure 10. The standard deviation of average curve of S wave H/V spectral ratio in predominant period and recording number in different stations
图 11 S波H/V平均谱比曲线在0.05—3s周期段内标准差的平均值
Figure 11. The mean value of standard deviation of S wave H/V spectral ratio in 0.05s-3s period
-
程佩青, 2001.数字信号处理教程.2版.北京:清华大学出版社. 大崎顺彦, 2008.地震动的谱分析入门.2版.北京:地震出版社. 何先龙, 佘天莉, 高峰, 2016.一种地震动P波和S波初至时间自动拾取的新方法.地球物理学报, 59(7):2519-2527. http://www.cnki.com.cn/Article/CJFDTotal-DQWX201607017.htm 冀昆, 2014.基于H/V谱比法的场地特征研究.哈尔滨: 中国地震局工程力学研究所. 罗桂纯, 2015.基于强震动观测的场地效应非线性特性研究.北京: 中国地震局地球物理研究所. 任叶飞, 温瑞智, 周宝峰等, 2014.2013年4月20日四川芦山地震强地面运动三要素特征分析.地球物理学报, 57(6):1836-1846. http://www.cqvip.com/QK/94718X/20146/661852468.html 荣棉水, 李小军, 王振明等, 2016.HVSR方法用于地震作用下场地效应分析的适用性研究.地球物理学报, 59(8):2878-2891. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb201608014 温瑞智, 冀昆, 任叶飞等, 2015a.基于谱比法的我国强震台站场地分类.岩石力学与工程学报, 34(6):1236-1241. http://www.cnki.com.cn/Article/CJFDTotal-YSLX201506018.htm 温瑞智, 王宏伟, 任叶飞等, 2015b.芦山余震震源参数及震源区品质因子反演.哈尔滨工业大学学报, 47(4):58-63. http://d.old.wanfangdata.com.cn/Periodical/hebgydxxb201504011 王伟君, 刘澜波, 陈棋福等, 2009.应用微动H/V谱比法和台阵技术探测场地响应和浅层速度结构.地球物理学报, 52(6):1515-1525. doi: 10.3969/j.issn.0001-5733.2009.06.013 王勤彩, 刘杰, 郑斯华等, 2005.云南地区与频率有关的P波、S波衰减研究.地震学报, 27(6):588-597, 699. doi: 10.3321/j.issn:0253-3782.2005.06.002 Akaike H., 1992. Information theory and an extension of the maximum likelihood principle. Breakthroughs in Statistics, 1:610-624. Baer M., Kradolfer U., 1987. An automatic phase picker for local and teleseismic events. Bulletin of the Seismological Society of America, 77(4):1437-1445. http://www.researchgate.net/publication/215755718_An_automatic_phase_picker_for_local_and_teleseismic_events Bonilla L. F., Steidl J. H., Gariel J. C., et al., 2002. Borehole response studies at the garner valley downhole array, Southern California. Bulletin of the Seismological Society of America. 92(8):3165-3179. doi: 10.1785/0120010235 Dimitriu P. P., Papaioannou C. A., Theodulidis N. P., 1998. EURO-SEISTEST strong-motion array near Thessaloniki, Northern Greece:A study of site effects. Bulletin of the Seismological Society of America. 88(3):862-873. Hassani B. Zafarani H., Farjoodi J., Ansari A., 2011. Estimation of site amplification, attenuation and source spectra of S-waves in the East-Central Iran. Soil Dynamics and Earthquake Engineering, 31(10):1397-1413. doi: 10.1016/j.soildyn.2011.05.017 Husid R., 1967. Gravity effects on the earthquake response of yielding structures. Pasadena, California: California Institute of Technology. Lermo J., Francisco J. Chávez-García F. J., 1993. Site effect evaluation using spectral ratios with only one station.Bulletin of the Seismological Society of America, 83(5):1574-1594. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=fbe8db440669882caac22472a98d6f26 Lee C. T., Cheng C. T., Liao C. W., et al., 2001. Site classification of Taiwan free-field strong-motion stations. Bulletin of the Seismological Society of America, 91(5):1283-1297. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=484aa1c629414dfc6c696fecfc77ee6d Mccann M. W. J., Shah H. C., 1979. Determining strong-motion duration of earthquakes. Bulletin of the Seismological Society of America, 69(4):1253-1265. http://www.researchgate.net/publication/245903749_Determining_Strong-Motion_Duration_of_Earthquakes?ev=auth_pub Ren Y. F., Wen R. Z., Yamanaka H., et al., 2013. Site effects by generalized inversion technique using strong motion recordings of the 2008 Wenchuan earthquake. Earthquake Engineering and Engineering Vibration, 12(2):165-184. doi: 10.1007/s11803-013-0160-6 Satoh T., Kawase H., Matsushima S., 2001.Differences between site characteristics obtained from microtremors, S-waves, P-waves, and codas. Bulletin of the Seismological Society of America, 91(2):313-334. doi: 10.1785/0119990149 Saragiotis C. D., Hadjileontiadis L. J., Panas S. M., 2002. PAI-S/K:A robust automatic seismic P phase arrival identification scheme. IEEE Transactions on Geoscience and Remote Sensing, 40(6):1395-1404. doi: 10.1109/TGRS.2002.800438 Wang H. W., Ren Y. F. Wen R. Z., 2018. Source parameters, path attenuation and site effects from strong-motion recordings of the Wenchuan aftershocks (2008-2013) using a non-parametric generalized inversion technique. Geophysical Journal International, 212(2):872-890. doi: 10.1093/gji/ggx447 Yadab P. D., Shin A., Takashi K., et al., 2017. Assessment of nonlinear site response at ocean bottom seismograph sites based on S-wave horizontal-to-vertical spectral ratios:A study at the Sagami Bay area K-NET sites in Japan. Earth Planets and Space, 69(29):1-7. doi: 10.1186/s40623-017-0615-5 Yamazaki F., Ansary M. A., 1997. Horizontal-to-vertical spectrum ratio of earthquake ground motion for site characterization. Earthquake Engineering and Structural Dynamics, 26(7):671-689. doi: 10.1002/(SICI)1096-9845(199707)26:7<671::AID-EQE669>3.0.CO;2-S Zare M., Bard P. Y., Ghafory-Ashtiany M., 1999. Site characterizations for the Iranian strong motion network. Soil Dynamics and Earthquake Engineering, 18(2):101-123. -
下载:
点击查看大图
计量
- 文章访问数: 106
- HTML全文浏览量: 51
- PDF下载量: 20
- 被引次数: 0
