近断层地震动的频谱周期参数和非平稳特征分析

分析地震动的频谱周期参数和频率非平稳特征,是地震危险性分析和结构抗震设计的重要内容.本文对近断层地震动记录按照不同的运动特征分组,进行了Hilbert-Huang变换分析.结合相关的定义得到了表征地震动频谱特性的7个周期参数,包括:Fourier幅值谱平均周期Tm、Hilbert边际谱平均周期Tmh、特征周期Tc、卓越周期等,并计算了地震动的Hilbert谱瞬时频率时程变异系数.计算结果表明,从整体上把握地震动频谱特性的三个周期参数Tm 、Tmh 和Tc适合于表征近断层地震动的频谱成份.近断层地震动的频率非平稳特征显著,Hilbert谱瞬时频率变异系数良好地反映了地震动的频率非平稳性质和程度.而且,考察了近断层地震动运动特征对地震动频谱周期参数和频率非平稳性的影响.     

Recording-based identification of site liquefaction

Reconnaissance reports and pertinent research on seismic hazards show that liquefaction is one of the key sources of damage to geotechnical and structural engineering systems. Therefore, identifying site liquefaction conditions plays an important role in seismic hazard mitigation. One of the widely used approaches for detecting liquefaction is based on the time-frequency analysis of ground motion recordings, in which short-time Fourier transform is typically used. It is known that recordings at a site with liquefaction are the result of nonlinear responses of seismic waves propagating in the liquefied layers underneath the site. Moreover, Fourier transform is not effective in characterizing such dynamic features as time-dependent frequency of the recordings rooted in nonlinear responses. Therefore, the aforementioned approach may not be intrinsically effective in detecting liquefaction. An alternative to the Fourier-based approach is presented in this study, which proposes time-frequency analysis of earthquake ground motion recordings with the aid of the Hilbert-Huang transform （HHT）, and offers justification for the HHT in addressing the liquefaction features shown in the recordings. The paper then defines the predominant instantaneous frequency （PIF） and introduces the PiF-related motion features to identify liquefaction conditions at a given site. Analysis of 29 recorded data sets at different site conditions shows that the proposed approach is effective in detecting site liquefaction in comparison with other methods.     

Long period microtremors, microseisms and earthquake damage: Northridge, CA, earthquake of 17 January 1994

Three studies of site amplification factors, based on the recorded aftershocks, and one study based on strong motion data, are compared one with another and with the observed distribution of damage from the Northridge, CA, earthquake of 17 January 1994 (M_L=6.4). In the epicentral area, when the peak ground velocities are larger than v_m≈15 cm/s, nonlinear response of soil begins to distort the amplification factors determined from small amplitude (linear) wave motion. Moving into the area of near-field and strong ground motion (v_m>30 cm/s), the site response becomes progressively more affected by the nonlinear soil response. Based on the published results, it is concluded that site amplification factors determined from small amplitude waves (aftershocks, small earthquakes, coda waves) and their transfer-function representation may be useful for small and distant earthquake motions, where soils and structures respond to earthquake waves in a linear manner. However in San Fernando Valley, during the Northridge earthquake, the observed distribution of damage did not correlate with site amplification determined from spectra of recorded weak motions. Mapping geographical distribution of site amplification using other than very strong motion data, therefore appears to be of little use for seismic hazard analyses.     

Site Response in Las Vegas Valley, Nevada from NTS Explosions and Earthquake Data

We report site response in Las Vegas Valley (LVV) from historical recordings of Nevada Test Site (NTS) nuclear explosions and earthquake recordings from permanent and temporary seismic stations. Our data set significantly improves the spatial coverage of LVV over previous studies, especially in the northern, deeper parts of the basin. Site response at stations in LVV was measured for frequencies in the range 0.2–5.0 Hz using Standard Spectral Ratios (SSR) and Horizontal-Vertical Spectral Ratios (HVR). For the SSR measurements we used a reference site (approximately NEHRP B ``rock' classification) located on Frenchman Mountain outside the basin. Site response at sedimentary sites is variable in LVV with average amplifications approaching a factor of 10 at some frequencies. We observed peaks in the site response curves at frequencies clustered near 0.6, 1.2 and 2.0 Hz, with some sites showing additional lower amplitude peaks at higher frequencies. The spatial pattern of site response is strongly correlated with the reported depth to basement for frequencies between 0.2 and 3.0 Hz, although the frequency of peak amplification does not show a similar correlation. For a few sites where we have geotechnical shear velocities, the amplification shows a correlation with the average upper 30-meter shear velocities, V₃₀. We performed two-dimensional finite difference simulations and reproduced the observed peak site amplifications at 0.6 and 1.2 Hz with a low velocity near-surface layer with shear velocities 600–750 m/s and a thickness of 100–200 m. These modeling results indicate that the amplitude and frequencies of site response peaks in LVV are strongly controlled by shallow velocity structure.     

应用HHT方法从竖向地震动记录识别场地液化的物理过程

利用HHT方法分析1987年美国Superstition Hills地震在Wildlife场地上获得的地震加速度记录,深入地论述了实际地震中土层液化对竖向地震动非平稳特性的影响,进而达到了应用HHT方法从地震动记录识别场地液化物理过程的目的。     

基于聚类经验模态分解（EEMD）的汶川MS8.0强震动记录时频特性分析

在2008年5月12日汶川MS8.0地震中,四川数字强震台网共获取了133组三分向加速度记录. 本文选取了一些不同断层距的台站所获取的强震动记录进行了处理和分析．在数据处理中,采用基于聚类经验模态分解（EEMD）提取信号时频特性的方法,有效获得了信号能量的时频分布,提取了中心频率、 Hilbert能量、最大振幅对应的时频等特性,并与傅里叶变换、小波变换进行了对比研究．研究结果表明, 对非线性的强震记录采用聚类经验模态分解（EEMD）能抑制经验模态分解（EMD）中存在的模态混叠问题;与傅里叶变换和小波变换相比发现, HHT边际谱在低频处幅值高于傅里叶谱;与小波变换受到所选取的母波强烈影响不同, HHT直接从强震记录中分离出固有模态函数（IMF）,更能反映出原始数据的固有特性, Hilbert谱反映出大部分能量都集中在一定的时间和频率范围内,而小波谱的能量却在频率范围内分布较为广泛．因此,基于EEMD的HHT在客观性和分辨率方面都具有明显的优越性,能提取到更多强震加速度记录的时频特性．      

International standard (ISO) on seismic actions for designing geotechnical works–An overview

The seismic performance of geotechnical works is significantly affected by ground displacement. In particular, soil–structure interaction and effects of liquefaction play major roles and pose difficult problems for engineers. An International Standard, ISO23469, is being developed for addressing these issues in a systematic manner within a consistent framework. The objective of this paper is to give an overview of this International Standard.In this International Standard, the seismic actions are determined through two stages. The first stage determines basic seismic action variables, including the earthquake ground motion at the site, the potential for earthquake-associated phenomena such as liquefaction and induced lateral ground displacement. These basic variables are used, in the second stage, for specifying the seismic actions for designing geotechnical works. In the second stage, the soil–structure interaction plays a major role. Types of analyses are classified based on a combination of static/dynamic analyses and the procedure for soil–structure interaction classified as follows:

– simplified: soil–structure interaction of a global system is modeled as an action on a substructure;

– detailed: soil–structure interaction of a global system is modeled as a coupled system.

Implications of the 2011 M9.0 Tohoku Japan earthquake for the treatment of site effects in large earthquakes

Site response in Japan is characterized using thousands of surface and borehole recordings from events of moment magnitude $(\mathbf{M}) > 5.5$ collected by the KiK-net network, including the 2011 M9.0 Tohoku earthquake. Site amplification is defined by the ratio of motions at the surface to those at depth (within the borehole), corrected for the depth effect due to destructive interference using a technique based on cross-spectral ratios between surface and down-hole motions. Site effects were particularly strong at high frequencies, despite the expectation that high-frequency response may be damped by nonlinear effects. In part, the large amplitudes at high frequencies are due to the prevalence of shallow soil conditions in Japan. We searched for typical symptoms for soil nonlinearity, such as a decrease in the predominant frequency and/or amplification, using spectral ratios of weak to strong ground motions. Localized nonlinearity occurred at some recording sites, but was not pervasive. We developed a general empirical model to express site amplification for the KiK-net sites as a function of common site variables, such as the average shear-wave velocity in the uppermost 30 m ( $\text{ V}_\mathrm{S30})$ and the horizontal-to-vertical (H/V) spectral ratio. We use the model to estimate site-corrected ground-motions for the Tohoku mainshock for a reference site condition; these motions are in reasonable agreement with the predictions of some of the published ground motion prediction equations for subduction zones.     

Hilbert-Huang变换在密频结构阻尼识别中的应用

Hilbert—Huang变换是一种新的数据处理方法，由经验模分解(Empirical Mode Decomposition)技术及Hilbert变换两部分组成。本文研究此方法对于密频结构阻尼识别的应用。首先对于两自由度系统模型，说明该方法用于阻尼识别的步骤。进而研究存在频率密集现象的高层建筑的阻尼识别问题。上述结果与理论值及由半功率带宽法的识别值进行了比较，对比显示Hilbert．Huang方法较传统方法具有良好的识别密频结构阻尼的性能，适用于大型结构的系统识别。     

线性SDOF体系共振动力响应的Hilbert谱分析

利用HHT方法分析了线性SDOF体系在共振状态下的强迫动力响应,系统地研究了线性SDOF体系在简谐波、线性调频波与正弦调频波输入下处于共振状态时结构动力响应的H ilbert谱与本征模态函数的特征,并且将H ilbert谱、Morlet小波谱及Fourier幅值谱所蕴含的物理意义及分辨率特性进行了比较。研究结果表明,与Fourier幅值谱相比,H ilbert谱在刻画非平稳时程的某些时变特征方面更加直观,其物理意义更加明确,而与Morlet小波谱相比H ilbert谱则具有较高的分辨率。     

Site effects by generalized inversion technique using strong motion recordings of the 2008 Wenchuan earthquake

The generalized inversion of S-wave amplitude spectra from the free-field strong motion recordings of the China National Strong Motion Observation Network System (NSMONS) are used to evaluate the site effects in the Wenchuan area. In this regard, a total of 602 recordings from 96 aftershocks of the Wenchuan earthquake with magnitudes of M3.7-M6.5 were selected as a dataset. These recordings were obtained from 28 stations at a hypocenter distance ranging from 30 km to 150 km. The inversion results have been verified as reliable by comparing the site response at station 62WUD using the Generalized Inversion Technique (GIT) and the Standard Spectral Ratio method (SSR). For all 28 stations, the site predominant frequency F _p and the average site amplification in different frequency bands of 1.0–5.0 Hz, 5.0–10.0 Hz and 1.0–10.0 Hz have been calculated based on the inversion results. Compared with the results from the horizontal-to-vertical spectral ratio (HVSR) method, it shows that the HVSR method can reasonably estimate the site predominant frequency but underestimates the site amplification. The linear fitting between the average site amplification for each frequency band and the V _s20 (the average uppermost-20 m shear wave velocity) shows good correlation. A distance measurement called the asperity distance D _Aspt is proposed to reasonably characterize the source-to-site distance for large earthquakes. Finally, the inversed site response is used to identify the soil nonlinearity in the main shock and aftershocks of Wenchuan earthquake. In ten of the 28 stations analyzed in the main shock, the soil behaved nonlinearly, where the ground motion level is apparently beyond a threshold of PGA > 300 cm/s² or PGV > 20 cm/s, and only one station coded 51SFB has evidence of soil nonlinear behavior in the aftershocks.     

HORIZONTAL-TO-VERTICAL SPECTRUM RATIO OF EARTHQUAKE GROUND MOTION FOR SITE CHARACTERIZATION

Nakamura's method, which uses a horizontal-to-vertical Fourier spectrum ratio of microtremor, has become popular to determine the predominant period and amplification of a site. In this study, this method is extended for earthquake ground motion recordings using new strong motion data recorded by JMA-87-type accelerometers. From the analysis of these accelerograms, horizontal-to-vertical Fourier spectrum ratios of a site for different earthquakes are also found to be stable irrespective of magnitude, distance and depth. To establish this fact, attenuation relations of velocity response spectra for horizontal and vertical components are derived for three damping ratios (0, 2 and 5 per cent) using the JMA data. Then the horizontal-to-vertical ratios of the velocity response spectra are obtained. The results show that the horizontal and vertical velocity response spectra are dependent on magnitude, distance and depth, but that their ratios are almost independent of magnitude, distance and depth. However, since the current data set consists of mostly intermediate to far field data, this observation should be limited to records of these distance ranges. Introducing station coefficients, representing site amplification, to this relation yields the value comparable to the horizontal-to-vertical Fourier spectrum ratio at a specific site. The stability of the spectrum ratio is explained by the transfer function between the ground surface and stiff-soil outcrop due to S-wave propagation. These results suggest that site amplification characteristics can be evaluated by one-point two-component surface recordings of earthquake ground motion, in a similar manner as proposed by Nakamura for microtremor. © 1997 by John Wiley & Sons, Ltd.     

Structural modal parameter identification and damage diagnosis based on Hilbert-Huang transform

Traditional modal parameter identifi cation methods have many disadvantages,especially when used for processing nonlinear and non-stationary signals.In addition,they are usually not able to accurately identify the damping ratio and damage.In this study,methods based on the Hilbert-Huang transform(HHT) are investigated for structural modal parameter identifi cation and damage diagnosis.First,mirror extension and prediction via a radial basis function(RBF) neural network are used to restrain the troublesome end-effect issue in empirical mode decomposition(EMD),which is a crucial part of HHT.Then,the approaches based on HHT combined with other techniques,such as the random decrement technique(RDT),natural excitation technique(NExT) and stochastic subspace identifi cation(SSI),are proposed to identify modal parameters of structures.Furthermore,a damage diagnosis method based on the HHT is also proposed.Time-varying instantaneous frequency and instantaneous energy are used to identify the damage evolution of the structure.The relative amplitude of the Hilbert marginal spectrum is used to identify the damage location of the structure.Finally,acceleration records at gauge points from shaking table testing of a 12-story reinforced concrete frame model are taken to validate the proposed approaches.The results show that the proposed approaches based on HHT for modal parameter identifi cation and damage diagnosis are reliable and practical.     

Side-effect of using response spectral amplification ratios for soft soil sites—Earthquake source-type dependent amplification ratios

Our previous studies show that site effects (amplification of rock motions), source and path effects are coupled when response spectra are used to characterize the amplification ratios for a soil site modelled as nonlinear or elastic. The coupling is referred to as a “side effect” of using response spectral amplification ratios. In the present study we use a suite of rock site records, well distributed with respect to magnitude and source distance, from crustal, subduction interface and slab earthquakes to evaluate the response spectral amplification ratio for soft soil sites. We compare these side-effects for ground motions generated by three types of earthquakes, and we find that, at periods much shorter or much longer than the natural period of a soil site modelled as elastic, the average amplification ratios with respect to rock site ground motions from three types of earthquakes are moderately different and are very similar for other spectral periods. These differences are not statistically significant because of the moderately large scatter of the amplification ratios. However, the extent of magnitude- and source-distance-dependence of amplification ratios differs significantly. After the effects of magnitude and source distance on the amplification ratios are accounted for, the differences in amplification ratios between crustal and subduction earthquake records are very large in some particular combinations of source distance and magnitude range. These findings may have potential impact in establishing design spectra for soft soil sites using strong motion attenuation models or numerical modelling.     

基于改进的希尔伯特-黄变换的地震强震动特征提取

2021年5月21日漾濞M_s6.4地震是云南近10 a以来继M_s6.5鲁甸地震和M_s6.6景谷地震后发生的又一次破坏性浅源地震,其震中位于滇西北地区维西-乔后断裂带附近。震源机制结果显示：此次地震属于走滑型破裂,符合区域构造特征。为进一步研究该地震强震动特征,特引入改进后的希尔伯特-黄变换(HHT),以不同角度客观分析强震动的尺度和频域特征。研究结果表明：希尔伯特-黄变换在对实际地震动记录进行特征提取后,得到的边际谱和时频谱可在一定程度上保留原始数据的主要信息,该方法可为类似破坏性浅源地震的强震动特征分析提供更有效的信息参数,保证了信号分解的稳定性,更适于处理非平稳信号,可为现代信号应用分析于强震动特征提供另一种新思路。     

2017年九寨沟MS7.0地震场地特征分析

场地条件对地震动具有较大影响, 研究不同场地条件下的地震动特征对地震动的校正具有重要意义。 本文以九寨沟地震为例, 收集了66个台站的198条三分量强震观测记录和SRTM(Shuttle Radar Topography Mission)公里格网的DEM(Digital Elevation Model)数据, 从多个角度对场地特征进行分析。 首先根据坡度法使用DEM数据对九寨沟特征进行了场地分类, 然后讨论了不同场地类型下的加速度时程、 加速度反应谱、 地震动衰减等地震动特征。 研究表明, 缺少实测资料的情况下, 地形坡度可以作为Vs30(地表以下30 m范围的平均剪切波速度)的一种替代指标, 利用坡度法可以较快速地对场地进行分类; 利用强震观测记录能够从多角度对台站的场地特征进行分析, 不同场地类型对地震动影响不同, 其中, 土层对地震动具有明显的放大效应。 该结论可以为地震动结果的校正提供依据。     

HHT方法在地震工程中的应用之初步探讨

详细介绍了一种新的非线性、非平稳信号的处理方法-HHT，讨论了其有效性和优越性，简单介绍了其在地震工程中的应用前景，并利用HHT处理了地震工程中常用El Centro地震记录，得到了该记录的振幅-频率-时间分布特性，并将得到的边际谱与傅立叶谱做了对比。     

Centrifuge modeling of seismic response of layered soft clay

Centrifuge modeling is a valuable tool used to study the response of geotechnical structures to infrequent or extreme events such as earthquakes. A series of centrifuge model tests was conducted at 80g using an electro-hydraulic earthquake simulator mounted on the C-CORE geotechnical centrifuge to study the dynamic response of soft soils and seismic soil–structure interaction (SSI). The acceleration records at different locations within the soil bed and at its surface along with the settlement records at the surface were used to analyze the soft soil seismic response. In addition, the records of acceleration at the surface of a foundation model partially embedded in the soil were used to investigate the seismic SSI. Centrifuge data was used to evaluate the variation of shear modulus and damping ratio with shear strain amplitude and confining pressure, and to assess their effects on site response. Site response analysis using the measured shear wave velocity, estimated modulus reduction and damping ratio as input parameters produced good agreement with the measured site response. A spectral analysis of the results showed that the stiffness of the soil deposits had a significant effect on the characteristics of the input motions and the overall behavior of the structure. The peak surface acceleration measured in the centrifuge was significantly amplified, especially for low amplitude base acceleration. The amplification of the earthquake shaking as well as the frequency of the response spectra decreased with increasing earthquake intensity. The results clearly demonstrate that the layering system has to be considered, and not just the average shear wave velocity, when evaluating the local site effects.     

Newmark design spectra considering earthquake magnitudes and site categories

Newmark design spectra have been implemented in many building codes, especially in building codes for critical structures. Previous studies show that Newmark design spectra exhibit lower amplitudes at high frequencies and larger amplitudes at low frequencies in comparison with spectra developed by statistical methods. To resolve this problem, this study considers three suites of ground motions recorded at three types of sites. Using these ground motions, influences of the shear-wave velocity, earthquake magnitudes, source-to-site distances on the ratios of ground motion parameters are studied, and spectrum amplification factors are statistically calculated. Spectral bounds for combinations of three site categories and two cases of earthquake magnitudes are estimated. Site design spectrum coefficients for the three site categories considering earthquake magnitudes are established. The problems of Newmark design spectra could be resolved by using the site design spectrum coefficients to modify the spectral values of Newmark design spectra in the acceleration sensitive, velocity sensitive, and displacement sensitive regions.     

Hilbert-Huang变换在提取地震信号动力特性中的应用

H ilbert-Huang变换(HHT)是一种处理非线性、非平稳信号的新方法。它通过经验模态分解将信号分解为有限个固有模态函数,并对每个固有模态函数进行H ilbert变换得到H ilbert谱。本文将这种方法应用于地震信号动力特性的提取,有效地获得了信号能量的时频分布,量化提取了中心频率、瞬时相位、瞬时能量、H ilbert能量、最大振幅对应的时频分布等动力特性,并与Fourier变换、小波变换等进行了比较,显示了HHT的优势以及对于进一步实现结构分析和控制的重要意义。