根据渐进谱的统计规律生成地震加速度时程

基于美国西部80条基岩上的近场强震记录, 采用Nakayama方法生成记录的渐进谱, 并参照Kameda方式,用统计方法建立了根据震级、距离等地震参数预测渐进谱的统计模型. 提出一种以渐进谱为目标谱的生成幅值和频率非平稳地震加速度时程的迭代方法. 由于考虑了渐进谱幅值和相位的相互影响,所生成的时程的相位也是时频非平稳的,并在相位调整中识别了相位谱增量符号以加速迭代收敛进程. 最后根据统计回归的目标渐进谱模型和本文提出的拟合目标渐进谱的方法,可生成不同震级、距离条件下的幅值和频率均非平稳的地震加速度时程.      

Generation of spectrum compatible accelerograms

In this paper a new method is presented for generating earthquake accelerograms which have pre-established response spectra. The non-stationary random nature and other salient features of the accelerograms can be taken care of by the procedure developed. The method leads to a sample spectrum which lies above the given spectrum. The generation of records to suit several spectra simultaneously can also be handled by this approach. The method is detailed first. This is followed by several numerical examples.     

Simulation of three-dimensional strong ground motions along principal axes,San Fernando earthquake

Power spectral density which describes frequency content is considered one of the most significant properties to be taken into account when generating ground motions through the use of stochastic processes. Using a smoothed and normalized Fourier amplitude spectrum, frequency content for components of motion along a set of principal axes is estimated. Fourier amplitude spectra obtained by this moving-window technique are presented which show the time dependency of frequency content for motions produced by the San Fernando earthquake of 9 February 1971. A mathematical model to simulate ground motion processes is proposed for which both the intensity and frequency content are non-stationary. Using this mathematical model with parameter characteristics along principal axes similar to those of the motions recorded during the San Fernando earthquake, three-dimensional ground motions are synthetically generated. The properties of the simulated motions show general characteristics similar to the characteristics observed in real accelerograms. The suggested model is considered adequate for engineering purposes.     

地震动随机函数模型研究(Ⅰ)——模型建立

基于随机过程的随机函数解释,发展了一种随机地震动模型.模型包括两个部分,分别是基于理论物理关系的Fourier幅值谱模型和基于经验关系的Fourier相位谱模型.基于拟加速度谱概念建立的幅值谱模型能够反映强震记录的幅值特性,而Fourier相位谱模型则利用累积相位谱的单值连续特性,以线性多项式与三角级数组合的形式给出.模拟记录与实测记录对比结果显示两者不仅在波形结构上相似性良好,而且反应谱层次亦吻合良好.     

Generating multiple spectrum compatible accelerograms using stochastic neural networks

A new neural‐network‐based methodology for generating artificial earthquake spectrum compatible accelerograms from response spectra was proposed in 1997, in which, the learning capabilities of neural networks were used to develop the knowledge of the inverse mapping from the response spectra to earthquake accelerograms. Recently, this methodology has been further extended and enhanced. This paper presents a new stochastic neural network that is capable of generating multiple earthquake accelerograms from a single‐response spectrum. A new stochastic feature to the neural network has been combined with a new scheme for data compression using the replicator neural networks developed in the original method. A benefit of this extended methodology is gaining efficiency in compressing the earthquake accelerograms and extracting their characteristics. The proposed method produces a stochastic ensemble of earthquake accelerograms from any response spectra or design spectra. An example is presented that used 100 recorded accelerograms to train the neural network and several design spectra and response spectra to test this improved methodology. Copyright © 2001 John Wiley & Sons, Ltd.     

地震动功率谱与反应谱的转换关系

本文详细评述了现今常用的几种加速度反应谱与功率谱的转换关系。对于小阻尼单质点体系而言，考虑到输入地震动是一非平稳的随机过程，又由于其传递函数的窄频带滤波特性，它的加速度反应将是一窄频带的非平稳随机过程。对于峰值系数水平的超越不是独立的，而是成群超越。据此本文考虑非平稳效应和对峰值系数水平的成群效应，对前人的转换关系进行了修正，并基于随机振动理论，给出了对工程上常用的频率平稳、强度非平稳的地震;动模型的功率谱和反应谱的转换关系。此转换关系对于长、短持时的地震动记录和反应谱长、短周期部分以及不同阻尼比的反应谱都能给出精度较高的结果。     

A modified method for simulating non-stationary multi-point earthquake ground motion

A spectral-representation-based algorithm is proposed to simulate non-stationary and stochastic processes with evolutionary power,according to a prescribed non-stationary cross-spectral density matrix. Non-stationary multi-point seismic ground motions at different locations on the ground surface are generated for use in engineering applications. First,a modified iterative procedure is used to generate uniformly modulated non-stationary ground motion time histories which are compatible with the prescribed power spectrum. Then,ground motion time histories are modeled as a non-stationary stochastic process with amplitude and frequency modulation. The characteristic frequency and damping ratio of the Clough-Penzien acceleration spectrum are considered as a function of time in order to study the frequency time variation. Finally,two numerical examples are presented to validate the efficiency of the proposed method,and the results show that this method can be effectively applied to the dynamic seismic analysis of long and large scale structures.     

地震波动强度非平稳特征提取数学建模分析

地震对人类生活安全有很大的威胁,对建筑结构有明显的影响,研究地震波动强度非平稳特征可为地震的预防提供参考。研究地震波动强度非平稳特征提取模型,分别设定确定性参数和随机参数的取值。对地震波动强度非平稳特性进行傅里叶分析,根据分析结果选取数学模型,并对模型参数进行设定;将记录的地震波动强度所对应的数振幅谱进行分解,拆分成一组不同尺度的分量,并对其重组创建连续分量,提取地震波振幅谱;利用地震波的振幅谱和相位谱之间的关系对地震波动强度非平稳特征进行提取。经过仿真实验证明,本方法提取地震波动强度非平稳特征的准确度更高,对分析高分辨率地震资料具有重要意义。     

Generation of artificial strong motion accelerograms

A method for generating synthetic strong motion accelerograms for use in engineering design is presented. This method utilizes the model proposed by Trifunac in 1971²⁷ in conjunction with the recent empirical scaling functions for characterization of amplitudes and duration of strong shaking in terms of (i) earthquake magnitude, M, and epicentral distance, R, or (ii) Modified Mercalli Intensity (MMI) at the recording station. The method also enables one to consider the desired levels of confidence that the synthetic motion will not be exceeded, direction of ground motion (horizontal or vertical) and the dispersive properties of geologic environment beneath and surrounding the station. The principal features of this approach are that the resulting accelerograms have non-stationary frequency and amplitude characteristics which are in full agreement with known principles of wave propagation through a stratified medium, and that the Fourier amplitudes and the frequency-dependent duration are scaled in accordance with known trends as in recorded accelerograms.     

Synthesis of accelerograms compatible with the Chinese GB 50011-2001 design spectrum via harmonic wavelets： artificial and historic records

A versatile approach is employed to generate artificial accelerograms which satisfy the compatibility criteria prescribed by the Chinese aseismic code provisions GB 50011-2001. In particular, a frequency dependent peak factor derived by means of appropriate Monte Carlo analyses is introduced to relate the GB 50011 -2001 design spectrum to a parametrically defined evolutionary power spectrum (EPS). Special attention is given to the definition of the frequency content of the EPS in order to accommodate the mathematical form of the aforementioned design spectrum. Further, a one-to-one relationship is established between the parameter controlling the time-varying intensity of the EPS and the effective strong ground motion duration. Subsequently, an efficient auto-regressive moving-average (ARMA) filtering technique is utilized to generate ensembles of non-stationary artificial accelerograms whose average response spectrum is in a close agreement with the considered design spectrum. Furthermore, a harmonic wavelet based iterative scheme is adopted to modify these artificial signals so that a close matching of the signals' response spectra with the GB 50011-2001 design spectrum is achieved on an individual basis. This is also done for field recorded accelerograms pertaining to the May, 2008 Wenchuan seismic event. In the process, zero-phase high-pass filtering is performed to accomplish proper baseline correction of the acquired spectrum compatible artificial and field accelerograms. Numerical results are given in a tabulated format to expedite their use in practice.     

Wavelet-based response spectrum compatible synthesis of accelerograms—Eurocode application (EC8)

An integrated approach for addressing the problem of synthesizing artificial seismic accelerograms compatible with a given displacement design/target spectrum is presented in conjunction with aseismic design applications. Initially, a stochastic dynamics solution is used to obtain a family of simulated non-stationary earthquake records whose response spectrum is on the average in good agreement with the target spectrum. The degree of the agreement depends significantly on the adoption of an appropriate parametric evolutionary power spectral form, which is related to the target spectrum in an approximate manner. The performance of two commonly used spectral forms along with a newly proposed one is assessed with respect to the elastic displacement design spectrum defined by the European code regulations (EC8). Subsequently, the computational versatility of the family of harmonic wavelets is employed to modify iteratively the simulated records to satisfy the compatibility criteria for artificial accelerograms prescribed by EC8. In the process, baseline correction steps, ordinarily taken to ensure that the obtained accelerograms are characterized by physically meaningful velocity and displacement traces, are elucidated. Obviously, the presented approach can be used not only in the case of the EC8, for which extensive numerical results/examples are included, but also for any code provisions mandated by regulatory agencies. In any case, the presented numerical results can be quite useful in any aseismic design process dominated by the EC8 specifications.     

Wavelet-based generation of spatially correlated accelerograms

For the seismic analysis of complex or nonlinear extended structures, it is useful to generate a set of properly correlated earthquake accelerograms that are consistent with a specified seismic hazard. A new simulation approach is presented in this paper for the generation of ensembles of spatially correlated accelerograms such that the simulated motions are consistent with (i) a parent accelerogram in the sense of temporal variations in frequency content, (ii) a design spectrum in the mean sense, and (iii) with a given instantaneous coherency structure. The formulation is based on the extension of stochastic decomposition technique to wavelet domain via the method of spectral factorization. A complex variant of the modified Littlewood-Paley wavelet function is proposed for the wavelet-based representation of earthquake accelerograms, such that this explicitly brings out the phase information of the signal, besides being able to decompose it into component time-histories having energy in non-overlapping frequency bands. The proposed approach is illustrated by generating ensembles of accelerograms at four stations.     

Application of arma models to estimate earthquake ground motion and structural response

This paper deals with the use of ARMA models in earthquake engineering. Tools and methods applied to strong ground motion are discussed emphasizing simulation of probabilistic earthquake response spectra. The ARMA models are applied to Icelandic earthquake data and a tentative model for Icelandic earthquakes is presented. This model, which is derived using 54 accelerograms, is based on a low-order, time-invariant ARMA process excited by Gaussian white noise and amplitude modulated using a simple envelope function to account for the non-stationary characteristics. This simple model gives a reasonable fit to the observed ground motion. Further, this model produces accurate earthquake response spectra, which, combined with accompanying attenuation and duration formulae, might be useful in earthquake hazard and risk assessment.     

非平稳地震动过程模拟方法(Ⅱ)

在总结回顾具有同一统计特征非平稳地震动过程模拟的研究基础上, 提出了一个模拟产生与某个给定地震动记录具有同一统计特征[CD2]时变功率谱的非平稳地震动过程的方法. 方法具有如下一些特点: 模拟产生的样本过程在幅值上和频率上均是非平稳的, 幅值和频率的非平稳性均取决于目标时变功率谱; 任意一个样本过程的时变功率谱未必符合目标谱, 但在统计意义上却严格符合目标谱. 最后, 论文以Landers地震的一个加速度记录为例对方法进行了验证.      

Simulation of artificial earthquakes

A procedure is developed for the simulation of artificial earthquake accelerograms, The time variation of amplitude and frequency content is preserved in the simulation procedure. Sixteen artificial earthquake accelerograms are simulated and compared with a target accelerogram. The time variation of amplitude and frequency content for 26 historical earthquake accelerograms is characterized.     

Simulation of non-stationary ground motion processes (II)

Introduction When doing Monte Carlo analysis of stochastic earthquake response of nonlinear structures, it is necessary to have a great number of ground motion processes with an identical statistical feature; In theoretical analysis of earthquake resistance of structures, to study the effect of one parameter, it is required to keep other parameters fixed, and it is often needed one set of ground motion proc-esses having an identical statistical feature; Similarly, while doing structural model…     

基于新《电力设施抗震设计规范》的地震动随机模型参数研究

根据新《电力设施抗震设计规范》(GB50260-2013)对随机地震动功率谱参数的取值进行分析。介绍新《电力设施抗震设计规范》中设计反应谱,推导设计谱到功率谱的转换过程。选取常用的Clough-Penzien修正过滤白噪声模型作为拟合函数,通过功率谱参数拟合,得到不同场地类别、不同抗震设防烈度下的地震动模型参数。该研究成果可为电力工程抗震设计分析提供参考。     

New method of generating spectrum compatible accelerograms using neural networks

A new method is proposed for generating artificial earthquake accelerograms from response spectra. This method uses the learning capabilities of neural networks to developed the knowledge of the inverse mapping from the response spectra to earthquake accelerogram. In the proposed method the neural networks learn the inverse mapping directly from the actual recorded earthquake accelerograms and their response spectra. A two-stage approach is used. In the first stage, a replicator neural network is used as a data compression tool. The replicator neural network compresses the vector of the discrete Fourier spectra of the accelerograms to vectors of much smaller dimension. In the second stage, a multi-layer feed-forward neural network learns to relate the response spectrum to the compressed Fourier spectrum. A simple example is presented, in which only 30 accelerograms are used to train the two-stage neural networks. This example demonstrates how the method works and shows its potential. © 1998 John Wiley & Sons, Ltd.     

时变功率谱模型参数衰减规律研究

基于美国西部基岩强震加速度记录资料,采用多重滤波技术计算水平向和竖向地震加速度记录时程的时变功率谱值,并利用非线性最小二乘法拟合了对应于各条加速度记录时程的时变功率谱模型参数值.采用经验统计方法,分析了水平向和竖向加速度时变功率谱模型参数随震级、震中距及频率的变化规律,建立其衰减模型,并由多随机变量一致加权最小二乘法拟合得到了各衰减模型的系数值.通过与现有时变功率谱衰减规律的对比分析,证明了本文提出的衰减模型的合理性.     

A non-parametric wave type based model for real-time prediction of strong ground motion accelerogram

A wave type based method for real-time prediction of strong ground motion (SGM) accelerogram is developed. Real-time prediction of SGM is requested in predictive building control systems to trigger and control actuator systems achieving the goal of reduction of the structural deformations during an on-going earthquake. It is well known that SGM is a classic example of non-stationary stochastic process with temporal variation of both amplitude and frequency content. The developed non-parametric model considers the non-homogeneity of the seismic process which contains different wave types with the individual frequency contents and time-dependency amplitude distribution pattern. Therefore, an important part of the method is to detect dominant seismic wave phases. Prediction of seismic signal is undertaken by applying frequency adaptive windowing approach, which leads to predict the on-coming signal in time window t+Δt based on the measured data in the time window t. Besides use of the frequency adaptive windowing, constant windowing and semi-adaptive windowing approaches are deployed. The results show that use of the adaptive time windows relevant to dominant frequency of the signal will enable the model to catch and predict the most dominant frequencies. Performance of the proposed model is verified by the use of 97 free-field accelerograms, which were applied to train and validate the prediction model. The selected accelerograms were measured above the soil type C and D according Eurocode 8 and their Moment magnitude are ranging between 6.2 and 7.7. The learning capability of the radial basis function Artificial Neural Network is used to reconstruct the SGM accelerogram. The most significant advantage of the proposed model is the concept of wave type based modeling which has the advantage of a conceptual physical modeling of the seismic process. Comparison of the real-time predicted and the observed accelerograms shows a high correlation when the frequency adaptive approach is applied. This paper lays a foundation for more effective use of real-time predictive control systems and potential for future extension in active structural control as well as in real-time seismology.