Characterizing and quantifying earthquake-induced site nonlinearity

This study proposes a recording-based approach to characterize and quantify earthquake-induced site nonlinearity, exemplified as soil nonlinearity and/or liquefaction. Alternative to Fourier spectral analysis (FSA), the paper introduces time–frequency analysis of earthquake ground motion recordings with the aid of so-called Hilbert-Huang transform (HHT), and offers justification for the HHT in addressing the nonlinear features shown in the recordings. With the use of the 2001 Nisqually earthquake recordings, this study shows that the proposed approach is effective in characterizing site nonlinearity and quantifying the influences in seismic ground responses. Major results from this study are listed below.

1. HHT-based site amplification is defined as the ratio of marginal Hilbert amplitude spectra, similar to the Fourier-based one that is the ratio of Fourier amplitude spectra. The HHT-based site amplification can be used in effectively quantifying site nonlinearity in terms of frequency downshift in the low-frequency range and amplification-reduction factor in intermediate-frequency range in comparison with the Fourier-based one.

2. Instantaneous damping, and Hilbert and marginal Hilbert damping spectra are defined in ways similar to instantaneous frequency, and Hilbert and marginal Hilbert amplitude spectra, respectively. Consequently, the HHT-based site damping is found as the difference of marginal Hilbert damping spectra, which can be used as an alternative, complementary index to measure the influences of site nonlinearity in seismic ground responses.

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

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

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

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

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

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

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

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

A new approach to liquefaction hazard zonation: Application to Laoag City,Northern Philippines

To help mitigate liquefaction hazards in the Philippines, an inexpensive yet effective approach to liquefaction hazard zonation was developed in this study. The proposed approach is also useful in other areas especially where funds for more rigorous procedures may not be available. The approach utilizes the geomorphology-based criteria to identify liquefaction-prone deposits based on geology and grain characteristics, and generate a preliminary liquefaction susceptibility map. Then, microtremor recordings, popularly used in site effect estimation, are gathered to derive qualitative information on the density and thickness of these deposits and generate a site classification map. This latter map is also essentially a ground shaking hazard map in that it shows those areas where thick, soft deposits likely to amplify and prolong the duration of ground motion can be found. Therefore, it also identifies areas where seismic demand can be high that the possibility of liquefaction being triggered is likewise high. Combining the two maps, an integrated liquefaction hazard zonation map is produced which provides not only an improved characterization of the soils’ capacity to resist liquefaction but also integrates qualitative information on the seismic demand on these deposits as well. With information about the relative thickness of the deposits, the severity of potential damage can likewise be inferred from the map since thicker deposits relate to more serious damage. The proposed approach was applied to Laoag City, Northern Philippines, where it was shown to reliably identify areas that are vulnerable to the hazard.     

基于希尔伯特变换的非平稳地震动模拟方法的验证

基于希尔伯特变换的非平稳地震动模拟方法能够生成一系列地震动样本,具有与给定天然地震动相似的非平稳特征. 该方法在模拟天然地震动时域波形的形态和能量的时频分布方面,效果非常显著且已得到验证. 本文以理想单自由度体系弹塑性地震反应为基础,深入验证该方法在模拟天然地震动对结构影响效应方面的效果. 研究结果表明,该方法能够较好地模拟天然地震动的工程特性,其生成的地震动样本的峰值特性及其引起的结构弹塑性地震响应具有与天然地震动相似的特征.      

Ground roll attenuation using a time-frequency dependent polarization filter based on the S transform

The ground roll and body wave usually show significant differences in arrival time, frequency content, and polarization characteristics, and conventional polarization filters that operate in either the time or frequency domain cannot consider all these elements. Therefore, we have developed a time-frequency dependent polarization filter based on the S transform to attenuate the ground roll in seismic records. Our approach adopts the complex coefficients of the S transform of the multi-component seismic data to estimate the local polarization attributes and utilizes the estimated attributes to construct the filter function. In this study, we select the S transform to design this polarization filter because its scalable window length can ensure the same number of cycles of a Fourier sinusoid, thereby rendering more precise estimation of local polarization attributes. The results of applying our approach in synthetic and real data examples demonstrate that the proposed polarization filter can effectively attenuate the ground roll and successfully preserve the body wave.     

基于小波变换的层间隔震结构地震响应分析

利用小波多分辨率分析将地震动加速度分解为多频段小波分量,并运用复模态方法推导其计算层间隔震体系在地震作用下的动力响应公式,讨论各频段地震信号及结构响应的能量分配。同时利用小波时频工具分析地震动能量在时频域内的分布对层间隔震结构响应的影响,进而为考察地震动非平稳性对层间隔震结构非线性分析的影响提供方法。利用小波分析的以上优势,对一典型层间隔震结构分别进行弹性和弹塑性分析,结果表明弹性体系在地震作用下的响应可由该地震波各小波分量的响应叠加而得,地震动能量在时间上的集中会对层间隔震结构响应产生不利影响。     

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

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

Analysis of earthquake ground motions using an improved Hilbert–Huang transform

Some limitations of the Hilbert–Huang transform (HHT) for nonlinear and nonstationary signal processing are remarked. As an enhancement to the HHT, a time varying vector autoregressive moving average (VARMA) model based method is proposed to calculate the instantaneous frequencies of the intrinsic mode functions (IMFs) obtained from the empirical mode decomposition (EMD) of a signal. By representing the IMFs as time varying VARMA model and using the Kalman filter to estimate the time varying model parameters, the instantaneous frequencies are calculated according to the time varying parameters, then the instantaneous frequencies and the envelopes derived from the cubic spline interpolation of the maxima of IMFs are used to yield the Hilbert spectrum. The analysis of the length of day dataset and the ground motion record El Centro (1940, N–S) shows that the proposed method offers advantages in frequency resolution, and produces more physically meaningful and readable Hilbert spectrum than the original HHT method, short-time Fourier transform (STFT) and wavelet transform (WT). The analysis of the seismic response of a building during the 1994 Northridge earthquake shows that the proposed method is a powerful tool for structural damage detection, which is expected as the promising area for future research.     

人造地震动反应谱的非线性拟合

为了进一步提高人造地震动时程对目标反应谱的拟合精度和效率,本文提出了一种人造地震动反应谱非线性拟合方法.在人造地震动反应谱非线性拟合过程中,认为Fouricr幅值谱与反应谱调整项之间满足非线性(幂指数)关系.通过比例拉伸和折叠等数值技巧以促使模拟加速度向目标加速度逼近.算例分析结果表明,人造地震动反应谱非线性拟合方法具...     

Earthquake ground motion modeling of induced seismicity in the Groningen gas field

A physics‐based numerical approach is used to characterize earthquake ground motion due to induced seismicity in the Groningen gas field and to improve empirical ground motion models for seismic hazard and risk assessment. To this end, a large‐scale (20 km × 20 km) heterogeneous 3D seismic wave propagation model for the Groningen area is constructed, based on the significant bulk of available geological, geophysical, geotechnical, and seismological data. Results of physics‐based numerical simulations are validated against the ground motion recordings of the January 8, 2018, M_L 3.4 Zeerijp earthquake. Taking advantage of suitable models of slip time functions at the seismic source and of the detailed geophysical model, the numerical simulations are found to reproduce accurately the observed features of ground motions at epicentral distances less than 10 km, in a broad frequency range, up to about 8 Hz. A sensitivity analysis is also addressed to discuss the impact of 3D underground geological features, the stochastic variability of seismic velocities and the frequency dependence of the quality factor. Amongst others, results point out some key features related to 3D seismic wave propagation, such as the magnitude and distance dependence of site amplification functions, that may be relevant to the improvement of the empirical models for earthquake ground motion prediction.     

基于EMD,EEMD与CEEMD的信号 时频分析技术对比研究

时频分析方法在各个领域中广泛应用起来,特别是在工程领域中对结构和构件的变形以及内部损伤检测应用较多。HHT作为一种新兴的非线性,非平稳时频分析技术在探测和勘探领域中具有广阔的前景。文中重点对EMD,EEMD和CEEMD算法进行分析,通过对合成信号进行傅里叶变换和希尔伯特-黄变换,得到信号的频谱图,希尔伯特黄谱,边际谱等数据。通过对比分析,总结了三种方法的优缺点。结果表明,EMD,EEMD和CEEMD方法是进行非线性非平稳信号分析的有效工具。      

基于HHT的结构强震记录分析研究

本文利用HHT(Hilbert-HUang Transform)研究了结构强震记录的时频特性及结构动力特性。介绍了一座7层钢混框架结构及其强震观测台阵概况以及经历的地震情况,选择了在3次有代表性地震中的强震记录,利用一种新的非平稳信号处理方法HHT对记录进行了处理和分析,得到了该结构强震记录的时频幅值三维分布以及边际谱,并将边际谱与傅里叶谱进行了对比,识别了结构的自振频率。研究表明,对结构强震记录这种强非平稳信号,可以利用HHT分析得到能量集中分布的频段与时间范围。HHT边际谱与傅里叶谱相比,在低频部分幅值要大于傅里叶谱,而在高频部分,幅值要小于傅里叶谱。利用结构强震记录识别的自振频率比环境振动测试结果要小。     

稀疏时频分解方法的研究与运用（英文）

Gabor变换和S变换是常用的时频分析工具。根据测不准原理,它们的时频分解结果无法在时间域和频率域同时具有很高的分辨率。为了提高非平稳信号时频分解结果的分辨率,本文提出瞬时频率分布函数（IFDF）并利用它表达非平稳信号。当非平稳信号时频成分的分布满足测不准原理对信号可分辨的要求时,瞬时频率分布函数的支集和短时Fourier变换的小波脊支集是同一个集合。利用IFDF的该特征,本文提出一种迭代算法（Sparse-STFT）实现了信号的稀疏时频分解。该算法在每次迭代过程中利用残留信号的短时Fourier变换结果的脊支集更新信号的时频成分,每次迭代得到的时频成分的叠加结果即为最终的稀疏时频分解结果。文中的数值实验证明了Sparse-STFT可以有效地提高非平稳信号时频分解结果的分辨率。最后,本文将该方法应用于地震数据面波的压制中,取得了理想的处理结果。     

古田地震水口水电站重力坝强震观测记录HHT分析

利用HHT变换研究2008年3月6日福建古田ML4.6地震中,水口水电站重力坝强震反应台阵获取的强震反应观测资料的时频特性及重力坝结构动力特性,得到如下结论:(1)强震记录能量集中在0～15Hz频段和10～15s时段;(2)重力坝顺河向主频为3.7Hz。分析结果对认识水库地震近场地震动特性和重力坝地震反应有一定的实际意义。     

基于离散谐小波变换的地震波时变谱估计及非平稳地震波人工合成

采用在频域内具有极好紧支性和盒形特性的谐小波作为母小波,推导了基于离散谐小波变换的地震波时变谱,并给出了不同尺度下能量随时间变化的包络函数的近似解析表达式．文中以Northridge地震波为例分析了时频分布特征．研究发现,对于自振频率为f的结构,其地震反应与输入地震波的时变谱在f 处的时域最大值（即）以及地震波在f 频率点附近的信号分量在时域内的能量集中程度有很大的关系;与db4波基相比,利用谐小波作为母小波的小波变换,其频域具有较好的分辨率,但时域分辨率却较差．最后提出了两种基于离散谐小波逆变换的人工非平稳地震波仿真方法．      

地震勘探信号时频分析方法对比与应用分析

针对地震勘探信号,对非平稳信号时频分析几种适效方法:短时Fourier变换、小波变换、S变换、Wigner分布、平滑伪Wigner分布、锥形核时频分布、AOK(adaptive optimum kernel,自适应最优核函数)分布等进行对比与应用研究.在阐明各种方法基本原理的基础上,进行数值分析与应用研究.首先对非平稳地震勘探模拟信号进行试算及时频属性提取,结合各类方法的信号表示理论,在时频局部化的精度和交叉项抑制等方面对计算结果进行对比分析;进一步应用于实际二维地震数据,提取瞬时频率和瞬时带宽等时频属性,进行比较研究.研究表明:对于地震勘探信号,就本文涉及的几种时频分析方法而言,AOK分布是时频局部化精度最高、交叉项抑制最好、时频匹配最优的方法,值得在地震勘探信号分析和地震属性提取、频谱分解等应用中深入研究和应用.     

Hilbert-Huang变换与地震信号的时频分析

本文介绍了HHT时频分析方法及瞬时频率的概念，给出了已知信号的经验模态分解及其时频分布，并对实际地震波形信号进行了HHT时频处理与剖析。结果表明，HHT方法能准确描述地震波形信号的非线性时变特征，是地震信号时频分析的有效工具。