汶川地震远场地震动场地相关性与分析方法评价

为考查远场地震动的场地相关性并评价一些场地特性分析方法的适用性,采用不同方法对汶川地震山东省12个远场台站的强震记录进行了分析.选取台站分别位于按建筑抗震设计规范(CBC)场地划分中的Ⅰ—Ⅲ类场地上.地震动记录的分析方法包括傅里叶幅值谱法,地震反应谱法,水平与竖向谱比率法,参考点谱比率法,以及尾波分析等.结果表明,按傅里叶幅值谱法,地震反应谱法,水平与竖向谱比法计算得到的卓越周期均远大于台站场地的卓越周期,不同方法得到的结果之间也有较大差别,且主要反映长周期地震动的卓越频率;参考点谱比率法的结果未反映地震动的卓越周期,也与场地的卓越周期差别较大;对完整记录尾波分析所得的结果比较接近场地的卓越周期.希望本文能为考虑远场地震作用时设计谱的建立,以及场地特性估计时地震动分析方法的选取提供参考依据.     

昭通巡龙强夯土非天然地震记录特征及识别研究

通过分析昭通巡龙2022年11月—2023年5月多次重型地面平整机施工和小震级天然地震事件的记录特征及差异后认为：昭通巡龙强夯土事件的P波初动向下、持续时间短、衰减快、P波与S波最大振幅小于天然地震,强夯土事件波形记录的频谱特征与天然地震差异较大,天然地震P波与S波震相明显,强夯土事件震相不明显。可以利用Pg/Sg谱比、交叉Pg/Sg谱比实现强夯土非天然地震和天然地震事件的分类识别。分类结果表明：当频率大于5 Hz时,Pg/Sg谱比判据基本实现天然地震和强夯土事件分类;交叉Pg/Sg谱比判据对应的分类识别正确率为天然地震事件91.3%,强夯土事件90.7%。交叉Pg/Sg谱比判据比单一频带谱比判据能更好地反映天然地震和强夯土事件的特征差异。     

安徽省国家强震动台站场地响应及背景噪声分析

利用安徽省部分强震动台站记录的地脉动数据和历史地震数据,采用H/V谱比法和噪声功率谱密度函数法,对比分析台站的场地响应及背景噪声,结果表明:全椒、马鞍山地震台的场地卓越周期小且背景噪声低;沈巷地震台卓越周期大且背景噪声高;三山地震台无明显卓越周期;地震台站背景噪声与场地响应曲线之间不存在一一对应关系。     

特殊长周期地震动的参数特征研究

近断层脉冲型地震动和远场软土层场地类谐和地震动是两类特殊的长周期地震动,当前的规范均很少对这两类地震作用进行具体的规定。研究了近断层脉冲型和远场类谐和地震动的幅值、幅值比（V／A,D／V）、傅里叶幅值谱和反应谱的差别,分析了相位角和作用循环周期数对简单脉冲的影响,并用于解释两类特殊地震动的工程特征。以集集地震动为数据基础,分析了两类长周期地震动的傅里叶谱和反应谱特征;将平均加速度和位移规准反应谱分别与规范设计谱进行了比较。建议设计谱在长周期段考虑近断层作用和软土场地面波效应的影响。     

利用宁洱地震强震记录分析场地响应

采用Nakamura方法,利用2007年6月3日宁洱地震强震动记录S波资料和脉动观测资料,计算了宁洱震区一些强震动台的场地响应.结果表明,当台站场地为基岩或覆盖层较薄时,由脉动数据和S波数据分别计算的谱比曲线一致性很好.当覆盖层大于35 m的时候,两者的谱比值有差异,即用S波计算的场地放大响应要明显大于用脉动计算的场地响应.     

汶川地震强地面运动模拟

基于确定性震源模型的方法主要用于计算低频(11 Hz)地震动常用经验格林函数法或随机方法,对低频地震动模拟不够准确．本文在确定性震源模型方法基础上,尝试采用分解给定的震源模型的方法来模拟宽频带(0.1——10 Hz)强地面运动,即采用分级离散断层面和分解断层面破裂单元上升时间的方法,增加震源时间函数中的高频信号,从而避免了对地震记录丰富程度和准确性的依赖．文中模拟计算了汶川MS8.0地震在8个地震观测台的地震动,将模拟结果和观测记录进行了加速度时程曲线和傅里叶振幅谱的对比分析．对比结果显示,模拟估计的地震动峰值加速度和持续时间与观测记录的数据基本在plusmn;50%的精度范围内相同,傅里叶振幅谱显示模拟结果有得到10 Hz左右的高频成分. 四川盆地中的台站模拟结果高频衰减比观测记录要快,原因是模拟过程没有考虑场地效应.对强地震动模拟还是要综合考虑震源、传播路径和场地的影响．研究结果表明,此改变震源输入的确定性方法可应用于模拟近断层宽频强地面运动．      

汶川地震强余震场地非线性反应

采用汶川地震强余震26个强震动台站记录,基于H/V谱比法,计算台站场地的强、弱震作用下场地卓越频率之比R_(fp)和体现强、弱震作用下场地H/V谱比曲线差异程度的DNL,进而识别并分析场地非线性反应特征。结果显示,R_(fp)、DNL与PGA之间存在显著相关性;其中10个台站出现明显的场地非线性反应特征;大部分台站发生场地非线性反应的PGA阈值为100 cm/s~2,部分为50 cm/s~2;实例证明,因某些台站不易通过H/V谱比法识别场地卓越频率,因此采用R_(fp)识别场地非线性反应具有一定局限性。     

2019年10月28日甘肃夏河5.7级地震强地面运动特征分析

2019年10月28日甘肃省夏河县发生5.7级地震,中国数字强震动台网的18个专业台站在此次地震中触发。本文处理捕获的54条三分向加速度记录,给出近场台站的地震动参数,绘制了震中附近区域峰值加速度等值线图,其长轴呈WN-ES方向展布。将实际观测数据与几种常用地震动衰减关系对比,发现俞言祥^[1]短轴衰减预测模型能更好地反映此次地震的影响场。将振幅较大的62LBL、62BLX台的反应谱与我国抗震设计反应谱比较,采用最小二乘法拟合出不同震中距5个台站各周期谱加速度衰减特性,总结了此次地震的反应谱基本特征。运用H/V谱比法对4个典型台站进行场地地震反应分析,研究了局部场地条件对峰值加速度和峰值速度的影响过程。     

震级、震源距对地震信号H/V影响分析

本文采用喜马拉雅科学探测台阵四川省境内10个台站,计算了所有地震信号0.5—15 Hz的H/V谱比,分析了震级、震源距对H/V谱比的影响,并与各台站场地地脉动H/V谱比进行了比较. 结果表明: 震级对各台站地震信号的H/V谱比振幅值以及峰值频率几乎没有影响;震源距对谱峰明显的H/V谱比曲线的谱形和峰值频率影响不大,但对谱峰不明显的H/V谱比曲线的峰值频率有一定影响;H/V谱峰明显的台站,地震信号与地脉动信号的H/V谱比曲线的峰值频率具有很高的一致性,H/V谱峰不明显的台站,二者的峰值频率有一定差异.      

基于ITD的强震动加速度记录数据处理方法

基于本征时间尺度分解算法(ITD)对比了二次拟合法、EMD算法对人工合成地震波的处理结果,证明了ITD算法对于强震动加速度记录数据的基线校正有较好的校正效果;进一步基于ITD算法对台湾集集地震和某强震动台阵未校正加速度记录进行处理,并采用带通滤波方法进行进一步滤波处理和变换得到校正后的加速度、傅里叶频谱、反应谱;最后通过梯形数值积分计算得到速度和位移.对比分析了ITD算法、EMD算法采用上述计算流程和强震分析软件的计算结果,验证了文中提出的强震动数据处理方法是稳定、可靠的,可以作为实时、快速的处理强震动加速度记录的一种有效方法.     

全球地震灾害发展趋势综述

基于EM-DAT（OFDA/CRED）国际灾害数据库及美国地质勘探局（USGS）给出的数据,提取全球1970—2018年干旱、地震、极端温度、极端气候、洪水、滑坡、火山活动、火灾和块体运动9种自然灾害的相关数据,并进行综合分析,着重分析死亡人数和经济损失。分析结果表明,21世纪以来地震灾害造成的损失最大。此外,2000—2017年全球重大地震尤其是9次特大地震的分析结果显示,可有效减轻地震灾害损失的措施至少包括建立地震预警机制、建设韧性城乡、建立完善的风险排查机制。     

Influence of parameters selection in Chebyshev filters on the strong motion data processing

This contribution discusses the application of Chebyshev Type I filter for processing real earthquake records. Consideration is given to the effects of filtering parameters (passband amplitude ripple and order of the filter) on the time series, strong-motion parameters, Fourier Amplitude Spectrum of acceleration, and elastic displacement response spectra. Time histories of five earthquakes with different moment magnitudes have been examined (from stations located close to the epicenters). Data processing is based on application of bandpass Chebyshev filtering over frequency range with substantial signal to noise ratio (level of 3 or approximately 3 dB). Applying different filters, we have monitored several important strong-motion parameters: peak values of acceleration, velocity, and displacement; Arias intensity, acceleration/velocity spectrum intensity, significant duration, etc. Some new results and conclusions concerning the influence of Chebyshev filter in data processing of records have been summarized. The graphical and numerical outcomes obtained, as well as the comparison with a Butterworth causal filter, are included in the work. The results could be potentially useful to engineering seismologists who need to evaluate and better understand the merits of this type of filtering for strong-motion data processing.     

Prediction of earthquake energy input from smoothed fourier amplitude spectrum

This paper attempts to show analytically that the energy-input spectra of damped SDOF systems and undamped MDOF systems excited by an earthquake motion can be predicted by smoothing the Fourier amplitude spectrum of the base acceleration. The spectral window for smoothing in the frequency domain for a damped SDOF system is identical with the probability density function of the time-variant or instantaneous vibration frequency resulting from non-linear hysteresis. The spectral window for an undamped MDOF system is identical with the set of squared participation factors associated with vibration modes. It was found that the increase in damping factor and the increase in participation of higher modes provide wider spectral windows, resulting in more flattened or unaltered energy-input spectra due to enhanced smoothing effects.     

Features of different types of seismic events in China’s Capital Region

Seismic records produced by different seismic sources vary. In this study, we compared the waveform records and time-frequency characteristics of tectonic earthquakes, artificial explosions, and mine collapses in China’s Capital Region. The results show that tectonic earthquakes are characterized by stronger S-wave energy than P-wave energy, obvious high-frequency components, and wide frequency bands of P and S waves. Artificial explosions are characterized by greater P-wave amplitude than S-wave amplitude and near-station surface wave development. Mine collapses are characterized by lower overall frequency, more obvious surface waves, and longer duration. We extracted quantitative discriminants based on the analysis of different event records, with 31 feature values in 7 categories (P/S maximum amplitude ratio, high/low frequency energy ratio, P/S spectral ratio, corner frequency, duration, the second-order moment of spectrum, and energy strongest point). A comparison of the ability of these feature values to recognize distinct events showed that the 6–17 Hz P/S spectral ratio was able to completely distinguish artificial explosions from the other two types of events. The S-wave corner frequency performed relatively well in identifying all three types of events, with an accuracy of over 90%. Additionally, a support vector machine was used to comprehensively distinguish multiple features, with an accuracy for all three types of events reaching up to 100%.     

Empirical evaluation of microtremor H/V spectral ratio

The objective of this work is to perform a purely empirical assessment of the actual capabilities of the horizontal-to-vertical (H/V) spectral ratio technique to provide reliable and relevant information concerning site conditions and/or site amplification. This objective has been tackled through the homogeneous (re)processing of a large volume of earthquakes and ambient noise data recorded by different research teams in more than 200 sites located mainly in Europe, but also in the Caribbean and in Tehran. The original recordings were first gathered in a specific database with information on both the sites and recorded events. Then, for all sites close to an instrumented reference, average site-to-reference spectral ratios (“spectral ratio method” (SSR)) were derived in a homogeneous way (window selection, smoothing, signal-to-noise ratio threshold, averaging), as well as H/V ratios (“HVSRE–RF”) on earthquake recordings. H/V ratios were also obtained from noise recordings at each site (either specific measurements, or extracted from pre- or post-event noise windows). The spectral curves resulting from these three techniques were estimated reliable for a subset of 104 sites, and were thus compared in terms of fundamental frequency, amplitude and amplification bandwidth, exhibiting agreements and disagreements, for which interpretations are looked for in relation with characteristics of site conditions. The first important result consists in the very good agreement between fundamental frequencies obtained with either technique, observed for 81% of the analyzed sites. A significant part of the disagreements correspond to thick, low frequency, continental sites where natural noise level is often very low and H/V noise ratios do not exhibit any clear peak. The second important result is the absence of correlation between H/V peak amplitude and the actual site amplification measured on site-to-reference spectral ratios. There are, however, two statistically significant results about the amplitude of the H/V curve: the peak amplitude may be considered as a lower bound estimate of the actual amplification indicated by SSR (it is smaller for 79% of the 104 investigated sites), and, from another point of view, the difference in amplitude exhibits a questioning correlation with the geometrical characteristics of the sediment/basement interface: large SSR/HV differences might thus help to detect the existence of significant 2D or 3D effects.     

A note on autoregressive spectral estimates for frequency-wavenumber analysis of strong-motion array data

Frequency-wavenumber (f-k) spectra of seismic strong-motion array data are useful in estimating back-azimuth and apparent propagation velocity of seismic waves arriving at the array. Such estimates are required to model wave passage effects while studying spatial variability of strong ground motion. Although periodogram-based spectral estimates are commonly used, practical applications based on them encounter limitations, such as, lack of objective criteria for selecting a proper smoothing window and its associated bandwidth, and relatively large variance of the estimated spectral quantities. We present an alternative spectral estimate based on parametric time series modelling approach. The well-known autoregressive (AR) time series model is used in a system-based approach to estimate the spectral matrix of auto- and cross-spectral densities. Such spectral estimates are found to be smoother than the windowed periodogram estimates, and can directly be used in f-k spectral analysis. We present an example application of the proposed technique using strong-motion data recorded by the SMART-1 array in Taiwan during the January 29 1981 $M_{L}$ 6.3 earthquake. Our results, in terms of back azimuth and apparent propagation velocity, are found to be in excellent agreement with those reported in the literature.     

Determination of horizontal and vertical design spectra based on ground motion records at Lali tunnel,Iran

Most acceleration diagrams show high levels of unpredictability, as a result, it is the best to avoid using diagrams of earthquake acceleration spectra, even if the diagrams recorded at the site in question. In order to design earthquake resistant structures, we, instead, suggest constructing a design spectrum using a set of spectra that have common characteristics to the recorded acceleration diagrams at a particular site and smoothing the associated data. In this study, we conducted a time history analysis and determined a design spectrum for the region near the Lali tunnel in Southwestern Iran. We selected 13 specific ground motion records from the rock site to construct the design spectrum. To process the data, we first applied a base-line correction and then calculated the signal-to-noise ratio (R_SN) for each record. Next, we calculated the Fourier amplitude spectra of the acceleration pertaining to the signal window (1), and the Fourier amplitude spectra of the associated noise (2). After dividing each spectra by the square root of the selected window interval, they were divided by each other (1 divided by 2), in order to obtain the R_SN ratio (filtering was also applied). In addition, all data were normalized to the peak ground acceleration (PGA). Next, the normalized vertical and horizontal responses and mean response spectrum (50%) and the mean plus-one standard deviation (84%) were calculated for all the selected ground motion records at 5% damping. Finally, the mean design spectrum and the mean plus-one standard deviation were plotted for the spectrums. The equation of the mean and the above-mean design spectrum at the Lali tunnel site are also provided, along with our observed conclusions.