随机有限断层法合成地震动的研究与应用

目前很多研究者采用随机有限断层模型预测地震动并被证明是可行的.本文运用随机有限断层地震动叠加合成方法,模拟了1988年肃南5.7级地震,与实际记录进行对比分析,发现：用有限断层方法模拟得到的峰值加速度以及峰值加速度出现所对应的特征时刻与真实记录较为吻合;比较模拟的地震动反应谱和实际记录的加速度反应谱,总体上在工程感兴趣的频段内,有限断层方法拟合的结果基本能满足实际工作的需要.为进一步开拓该方法在工程应用中的前景,以兰州市柴家峡水电站为例,将此理论方法应用于缺乏强震记录地区的近场地震动估计中,模拟分析了马衔山北缘活动断裂发震时在坝址区产生的地震动特征,其近场合成结果与1125年兰州7.0级地震的烈度分布符合较好,可供工程抗震设计参考使用.     

基于强度和能量控制的随机地震动模型

考虑地震动的随机性和频率与强度非平稳性,通过理论分析,提出了一般随机地震动模型,并给出了确定模型参数的原则和方法。该模型以地震动强度、地震动能量以及地震动持时等宏观指标作为控制随机地震动模型参数的指标,而对其内在的频谱组成等指标只要求满足一般地震动的特征。该模型可以用于描述平稳随机过程、强度非平稳随机过程以及强度和频率完全非平稳随机过程。通过与常用功率谱模型的比较,验证了该模型的合理性。     

A ground motion selection algorithm based on the generalized conditional intensity measure approach

An algorithm is presented for the selection of ground motions for use in seismic response analysis. The algorithm is based on the use of random realizations from the conditional multivariate distribution of ground motion intensity measures, IM|IM_j, obtained from the generalized conditional intensity measure (GCIM) approach. The algorithm can be applied to the selection of both as-recorded amplitude-scaled and synthetic/simulated ground motions. A key feature is that the generality of the GCIM methodology allows for ground motion selection based on only explicit measures of the ground motions themselves, as represented by the various IM’s considered, rather than implicit causal parameters (e.g., source magnitude, source-to-site distance) which are presently used in other contemporary ground motion selection procedures. Several examples are used to illustrate the salient features of the algorithm, including: the effect of intensity measures considered; and the properties of ground motions selected for multiple exceedance probabilities. The flexibility of the proposed algorithm coupled with the GCIM methodology allows for objective and consistent ground motion selection as a natural extension of seismic hazard analysis.     

Seismic ground motion simulations from a class of spatial variability models

An improved version of the spectral representation method (References 12–14) for the generation of simulations of random processes and random fields by means of the Fast Fourier Transform algorithm is developed. The new formulation produces simulations that are ergodic in the mean without any restricting assumptions at the origin of the spectra and leads to a faster convergence rate of the stochastic characteristics of the simulations to those of the random field. As an example of the methodology, simulations in space and time based on a spatial variability model for the seismic ground motions are generated.     

Attenuation characteristics of ground motion intensity from earthquakes with intermediate depth

Data of the intermediate depth (the Geioyo and the Shizuoka) earthquakes in Japan recorded in a dense network is analysed and compared with various available attenuation relations. The approach of Midorikawa (Midorikawa S., 1993, Technophysics 218, 287–295) based on the empirical Greens function technique of Irikura (Irikura, K., 1986, Proceedings of the 7th Japan Earthquake Engineering Symposium, pp. 151–156.) has been used to model the rupture responsible for these earthquakes and peak ground acceleration are simulated at selected observation points. The method presented in this work includes the transmission effect in a multiple layer crustal model for a finite fault earthquake source model. Sharp attenuation rate is seen for such intermediate depth earthquakes which is difficult to explain through conventional attenuation relations. Detail study of the methodology and comparison of results shows that the transmission factor plays an important role for the sharp attenuation rate for intermediate-to deep-focus earthquakes.     

地震动强度对场地地震反应的影响

通过对场地地震反应的研究发现,地震动强度不同,同样的场地对地震动的影响也不同,得到的反应谱形状也不同。目前我国抗震设计规范中不同烈度下的设计反应谱形状完全一致,没有考虑地震动强弱对相同场地反应谱形状的影响。     

Theoretical research on intensity envelope function of strong ground motion

（金星，瘳振鹏）Ｔｈｅｏｒｅｔｉｃａｌｒｅｓｅａｒｃｈｏｎｉｎｔｅｎｓｉｔｙｅｎｖｅｌｏｐｅｆｕｎｃｔｉｏｎｏｆｓｔｒｏｎｇｇｒｏｕｎｄｍｏｔｉｏｎ￥ＸｉｎｇＪＩＮａｎｄＺｈｅｎ－ＰｅｎｇＬＩＡＯ（ＩｎｓｔｉｔｕｔｅｏｆＥｎｇｉｎｅｅｒｉｎｇＭｅｃｈ...     

Intelligent determination of peak ground motion parameters based on intensity assessed by artificial neural network

Artificial Neural Network (ANN) model of computation based on mathematical model of neural processes is applied to establish an intelligent computing network from seismic intensity to peak ground parameter instead of the conventional statistical relationship in this paper. For a give seismic intensity rating, the network formed with actual strong ground motion records directly produces the corresponding peak ground parameters and the effects of earthquake magnitude and epicentral distance are included. The computed results of the network trained with a number of strong motion records in the West America show that such networks have obtained good conversion relationship from seismic intensity to peak ground parameters. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 208–216, 1993.     

Rating damage potential of ground motion records

In this study, a total of 115,246 ground motions recorded during earthquakes of Moment magnitudes ranging from M_w 5.0 to M_w 9.0 are analyzed statistically. A total of 21 ground motion parameters characterising the recorded acceleration time histories are used in the analysis. Classification of these parameters through statistical correlation is reported and a parameter called "distance from zero-amplitude axis," or d_Z-A, is formulated in the principal component space. The ability for d_Z-A to rate the damage potentials of strong motion records is evaluated through correlation of d_Z-A with Japan Meteorological Agency(JMA) instrumental seismic intensities. This parameter can be used to rate damage potential of any strong motion record irrespective of the magnitude and location of the earthquake. It can also be used in selecting ground motion records of appropriate damage potential in seismic design and probabilistic analysis.     

The effect of glaciation on the intensity of seismic ground motion

Seismicity is known to contribute to landscape denudation through its role in earthquake‐triggered slope failure; but little is known about how the intensity of seismic ground motions, and therefore triggering of slope failures, may change through time. Topography influences the intensity of seismic shaking – generally steep slopes amplify shaking more than flatter slopes – and because glacial erosion typically steepens and enlarges slopes, glaciation may increase the intensity of seismic shaking of some landforms. However, the effect of this may be limited until after glaciers retreat because valley ice or ice‐caps may damp seismic ground motions. Two‐dimensional numerical models (FLAC 6.0) were used to explore how edifice shape, rock stiffness and various levels of ice inundation affect edifice shaking intensity. The modelling confirmed that earthquake shaking is enhanced with steeper topography and at ridge crests but it showed for the first time that total inundation by ice may reduce shaking intensity at hill crests to about 20–50% of that experienced when no ice is present. The effect is diminished to about 80–95% if glacier ice level reduces to half of the mountain slope height. In general, ice cover reduced shaking most for the steepest‐sided edifices, for wave frequencies higher than 3 Hz, and when ice was thickest and the rock had shear stiffness well in excess of the stiffness of ice. If rock stiffness is low and shear‐wave velocity is similar to that of ice, the presence of ice may amplify the shaking of rock protruding above the ice surface. The modelling supports the idea that topographic amplification of earthquake shaking increases as a result of glacial erosion and deglaciation. It is possible that the effect of this is sufficient to have influenced the distribution of post‐glacial slope failures in glaciated seismically active areas. Copyright © 2012 John Wiley & Sons, Ltd.     

从集集地震看建筑物震害与地震动参数的关系

本文根据9．21台湾集集大地震中南投县、台中县、台北县和苗栗县的建筑物震害调查资料确定出这4个主要受害地区的平均震害指数,然后义计算出这4个地区的地震动峰值加速度、峰值速度、峰值谱加速度、峰值谱速度、峰值谱位移的平均值。将这些均值分别与平均震害指数进行基于最小二乘法的线性回归,发现峰值速度和震害指数有非常好的相关性,相关系数接近1．0,而峰值加速度与震害指数相关性只有0．8,这说明地震动峰值速度能够体现地震动对脆性结构的破坏势,与以往将不同次地震的记录和资料混淆在一起研究而获得峰值加速度与震害有很好相关性形成了对比。     

工程爆破场地地震动强度预测研究

本文介绍了爆破地震动强度的预测方法和该领域研究的国内外现状,以岭澳核电站场地为例,利用该场地土石方工程爆破时实监测获得的大量实验数据进行统计回归分析,提出了一种估计(折算)爆破药量的方法,给出预测爆破地震动加速度、速度、的经验公式。按此地震动衰减规律并结合时实监测数据指导爆破,确保大亚湾核电站的安全运行。这种折算爆破药量的方法和预测爆破地震动强度的经验公式对指导工程爆破具有重要参考价值。     

Three-dimensional simulations of strong ground motion in the Shidian basin based upon the spectral-element method

The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical profiles within the basin are used to analyze the generation and propagation of surface waves and their relation to the basin structures in both the time and frequency domain.The amplification effects are analyzed by the distribution of peak ground velocity(PGV)and cumulative kinetic energy(Ek) in the basin.The results show that in the 3D basin,the excitation of the fundamental and higher surface wave modes are similar to that of the 2D model.Small bowls in the basin have great influence on the amplification and distribution of strong ground motion,due to their lateral resonances when the wavelengths of the lateral surface waves are comparable to the size of the bowls.Obvious basin edge effects can be seen at the basin edge closer to the source for constructive interference between direct body waves and the basin-induced surface waves.The Ek distribution maps show very large values in small bowls and some corners in the basin due to the interference of waves propagating in different directions.A high impedance contrast model can excite more surface wave modes,resulting in longer shaking durations as well as more complex seismograms and PGV and Ek distributions.     

Influence of ground motion intensity on the effectiveness of tuned mass dampers

The effectiveness of Tuned Mass Dampers (TMD) on buildings subjected to moderate and high-intensity motions is analysed. First, the response of a 22-storey four-bay reinforced concrete non-linear frame with a TMD is studied for motions with different intensities. Several values of the relevant parameters are assumed in the analyses. Then, equivalent single-degree-of-freedom systems with TMDs and without them are defined and analysed under the action of ground motions with intensities associated with different return intervals at the site where the structures are located. Vulnerability curves for the systems are obtained based on the probabilities of reaching two different performance limit states. The expected annual rate of exceedance of each limit state is calculated. The results show that the effectiveness of TMDs is higher for systems with small non-linearity produced by small and moderate earthquakes, than for systems with high non-linear behaviour, generally associated with high-intensity motions. Some recommendations about the applicability of TMD are given. Copyright © 1999 John Wiley & Sons Ltd.     

脉冲型地震动潜在破坏势参数的相关性分析

为了在众多参数中挑选其中最有代表性的参数,来解释和反映脉冲型地震动对结构的潜在破坏能力,以338条脉冲型地震动记录作为研究对象,分析地震动参数与中低层结构响应的相关性。选取了14个常用地震动参数,对各地震动参数之间的相关性进行分析,从中选出7个代表性地震动参数;并将脉冲型地震动输入中低层结构模型中计算结构响应,分析代表性地震动参数与结构响应的相关性,与基于非脉冲型地震动的相关性计算结果进行对比。选用了3层和7层2个RC框架结构作为中低层结构代表,其基本周期为0.62s和0.89s。结果表明:对于脉冲型地震动,对于3层结构时与结构响应相关性最好的为EPV,对于7层结构时与结构响应相关性最好的为PGV,因此可以用PGV和EPV作为表征脉冲型地震动对中低层结构潜在破坏能力的参数;而对于非脉冲型地震动,与结构响应相关性最好的参数为PGV,可以用PGV作为表征脉冲型地震动对中低层结构的潜在破坏能力的参数。因此,通过地震动参数来解释和表征脉冲型地震动对结构的破坏能力是可行的。     

基于假设检验的地震动强度(烈度)速报方法

地震发生后数分钟内,快速、可靠地判别出地震动强度(烈度)的空间分布,用以估计不同地区的受灾程度,可以为政府开展应急救援并合理分配救援力量提供决策依据,保证救援人员及时、准确地到达极震区并展开搜救,以减少生命财产损失。本文基于统计学中的假设检验方法,对历史震害资料进行统计,提出了一种利用强地震动参数判别地震动强度(烈度)的方法。比较结果表明,本方法所确定的地震动强度(烈度)与实际震害烈度对应较好,能较真实地反映实际震害情况。     

不同地震动参数与地震烈度相关性对比研究

本文以美国加州3次破坏性地震和2008年中国汶川Ms8.0地震中所获得的强震记录作为数据源,应用最小二乘法进行线性回归,得到了不同地震动参数与烈度间的相关性.统计发现,在地面加速度峰值、地面速度峰值等16种描述地震动强度的参数中,烈度识别正确率及相关性最好的两个参数分别是标准累积绝对速度和谱烈度.     

An index for measuring the effects of topography on seismic ground motion intensity

An index is proposed for the quantitative evaluation of topographical effects, as may be required in seismic microzoning and structural design applications. This topographical effects index (TEI) is defined as the relative Arias intensity for a given position. As an example, the TEI is calculated for a semicircular canyon excited by an earthquake signal that consists of plane elastic SH waves and is characterized by a theoretical spectrum of accelerations. The analytical solution by Trifunac¹⁵ provides the transfer functions for this problem. The behaviour and spacial distribution of the calculated TEI values are illustrated for different incidence angles and spectral parameters. The proposed index appears useful for generating quantitative criteria which, after proper confirmation by instrumental recordings, could be considered in earthquake code provisions.     

Modification of stochastic ground motion models for matching target intensity measures

Stochastic ground motion models produce synthetic time‐histories by modulating a white noise sequence through functions that address spectral and temporal properties of the excitation. The resultant ground motions can be then used in simulation‐based seismic risk assessment applications. This is established by relating the parameters of the aforementioned functions to earthquake and site characteristics through predictive relationships. An important concern related to the use of these models is the fact that through current approaches in selecting these predictive relationships, compatibility to the seismic hazard is not guaranteed. This work offers a computationally efficient framework for the modification of stochastic ground motion models to match target intensity measures (IMs) for a specific site and structure of interest. This is set as an optimization problem with a dual objective. The first objective minimizes the discrepancy between the target IMs and the predictions established through the stochastic ground motion model for a chosen earthquake scenario. The second objective constraints the deviation from the model characteristics suggested by existing predictive relationships, guaranteeing that the resultant ground motions not only match the target IMs but are also compatible with regional trends. A framework leveraging kriging surrogate modeling is formulated for performing the resultant multi‐objective optimization, and different computational aspects related to this optimization are discussed in detail. The illustrative implementation shows that the proposed framework can provide ground motions with high compatibility to target IMs with small only deviation from existing predictive relationships and discusses approaches for selecting a final compromise between these two competing objectives.     

有限断层地震波辐射能估算及其在合成强地面运动中的应用

对有限断层地震波能量辐射的估算通常采用断层面上子源能量的逐点求和方法。基于Brune圆盘模型,Anderson推导出有限断层地震波能量辐射S波的求解公式,即 ,其中 为断层面上地震矩, 为剪切模量, 和 分别为动态应力降和静态应力降,并指出在复合震源模型强地面运动预测应用中 以满足能量守恒。Rivera和Kanamori则从能量辐射表象定理出发,给出了有限断层中辐射能量的积分表达式,明确地指出了逐点求和所存在的问题。依据该积分表达式,本文推导出了复合源模型中新的辐射能完整的求解方法,指出Anderson方法实为断层面上点源辐射能量的简单叠加求和,后者则充分考虑了断层面上任一点在任一时刻能量传播过程中受到的断层面上所有位移破裂路径的交互影响。以1976年唐山7.6地震为例,应用上述方法分别计算了有限断层模型的辐射能量及近场强地面运动,如质点运动加速度,速度。结果表明如果模型参数满足 时,由本文给出的求解方法计算所得到的地震波辐射能已远远超出实际的辐射能量值,直接导致了对近场强地面运动参数如质点速度、加速度等的过高估算。因此,Zeng等和Anderson工作的局限性是非常明显的:地震矩守恒以及非物理的 无法准确地预测近场地面运动。未来工作中,对于有限断层模型的建立,在地震矩守恒这一约束条件的基础上,远场和近场能量解（或视应力）将可作为另一个重要的约束条件,为强地面运动的模拟提供一个更为恰当的求解方案。