基于小田原试验场地数据的不同基岩输入模型对预测地震动特征的影响研究

随着强震台网的密布及观测记录的增加,为研究各类局部场地地震反应预测模型的合理性提供了有效的参考依据,也使利用强震记录及场地条件研究地震动特征成为可能。选取场地地质参数资料和地震记录数据齐全的日本小田原（Ashigara Valley）盲测试验场地,通过对比不同地震动输入方式及场地反应分析模型,研究地震动特征,分析现有模型的优劣。基于1990年8月5日M5.1强震事件的地表基岩记录和地下基岩地震记录,采用地下台强震记录直接输入、地表基岩台强震记录减半为基底地震动输入、地表基岩台强震记录反演为基底地震动输入作为3种基岩地震动输入。基于局部场地条件分别建立一维等效线性模型、二维黏弹性模型及二维时域等效线性化模型等工程中常用的场地数值分析模型,进行局部场地地震反应分析,预测该盲测场地的地表地震动特征,并与对应的实测强震记录结果进行对比,分析不同基岩地震动输入方式对预测地震动特征及地表土层反应谱特征的影响,重点分析地震动输入、土体非线性、场地横向不均匀性及几何与非线性特征共同作用等因素对地表地震动特征的影响,以期为地表地震动的合理预测提供参考。     

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.     

A centrifuge study of the influence of site response,relative stiffness,and kinematic constraints on the seismic performance of buried reservoir structures

The seismic performance of underground reservoir structures depends on the properties of the structure, soil, and ground motion as well as the kinematic constraints imposed on the structure. A series of four centrifuge experiments were performed to evaluate the influence of site response, structural stiffness, base fixity, and excitation frequency on the performance of relatively stiff reservoir structures buried in dry, medium-dense sand. The magnitude of seismic thrust increased and the distribution of seismic earth pressures changed from approximately triangular to parabolic with increasing structural stiffness. Heavier and stiffer structures also experienced increased rocking and reduced flexural deflection. Fixing the base of the structure amplified the magnitude of acceleration, seismic earth pressure, and bending strain compared to tests where the structure was free to translate laterally, settle, or rotate atop a soil layer. The frequency content of transient tilt, acceleration, dynamic thrust, and bending strain measured on the structure was strongly influenced by that of the base motion and site response, but was unaffected by the fundamental frequency of the buried structure (f_structure). None of the available simplified procedures could capture the distribution and magnitude of seismic earth pressures experienced by this class of underground structures. The insight from this experimental study is aimed to help validate analytical and numerical methods used in the seismic design of reservoir structures.     

Seismic response of deep Quaternary sediments in historical center of L’Aquila City (central Italy)

On 6 April 2009 a M_w=6.1 earthquake produced severe destruction and damage over the historic center of L’Aquila City (central Italy), in which the accelerometer stations AQK and AQU recorded a large amount of near-fault ground motion data. This paper analyzes the recorded ground motions and compares the observed peak accelerations and the horizontal to vertical response spectral ratios with those revealed from numerical simulations. The finite element method is considered herein to perform dynamic modeling on the soil profile underlying the seismic station AQU. The subsurface model, which is based on the reviewed surveys that were carried out in previous studies, consists of 200–400 m of Quaternary sediments overlying a Meso-Cenozoic carbonate bedrock. The Martin-Finn-Seed's pore-water pressure model is used in the simulations. The horizontal to vertical response spectral ratio that is observed during the weak seismic events shows three predominant frequencies at about 14 Hz, 3 Hz and 0.6 Hz, which may be related to the computed seismic motion amplification occurring at the shallow colluvium, at the top and base of the fluvial-lacustrine sequence, respectively. During the 2009 L’Aquila main shock the predominant frequency of 14 Hz shifts to lower values probably due to a peculiar wave-field incidence angle. The predominant frequency of 3 Hz shifts to lower values when the earthquake magnitude increases, which may be associated to the progressive softening of soil due to the excess pore-water pressure generation that reaches a maximum value of about 350 kPa in the top of fluvial-lacustrine sequence. The computed vertical peak acceleration underestimates the experimental value and the horizontal to vertical peak acceleration ratio that is observed at station AQU decreases when the earthquake magnitude increases, which reveals amplification of the vertical component of ground motion probably due to near-source effects.     

地铁车站的强地震反应分析及设计地震动参数研究

进行了地铁地下车站的地震反应分析,探讨了地铁车站地震反应的主要影响因素,介绍了地面与基岩间峰值相对位移的确定及其在地下结构抗震设计中的应用,初步研究了地铁车站埋深对结构地震反应的影响。分析结果表明,地震引起的地基变形是影响地下结构动力反应的决定性因素,结构峰值变形反应与自由场峰值变形反应之间近似存在简单的线性关系;相对于设计基本地震加速度,地面与基岩间峰值相对位移(PGRD)对于地下结构抗震分析及设计是一种更为合理的设计地震动参数。     

Non-linear site response simulations in Chang-Hwa region during the 1999 Chi–Chi earthquake, Taiwan

The destructive 1999 Chi–Chi earthquake (M_w 7.5) was the largest inland earthquake in Taiwan in the 20th century. Several observations witness the non-linear seismic soil response in sediments during the earthquake. In fact, large settlements as well as evidence of liquefaction attested by sand boils and unusual wet ground surface were observed at some sites. In this paper, we present a seismic response simulation performed with CyberQuake software on a site located within the Chang-Hwa Coastal Industrial Park during the 1999 Chi–Chi earthquake in Taiwan. A non-linear multi-kinematic dynamic constitutive model is implemented in the software. Computed NS, EW and UP ground accelerations obtained with this model under undrained and two-phase assumptions, are in good agreement with the corresponding accelerations recorded at seismic station TCU117, either for peak location, amplitudes or frequency content. In these simulations, liquefaction occurs between depths 1.3 and 11.3 m, which correspond to the observed range attested by in place penetration tests and other liquefaction analyses. Moreover, the computed shear wave velocity profile is very close to post-earthquake shear wave velocity profile derived from correlations with CPT and SPT data. Finally, it is shown that in non-linear computations, even though a 1D geometry is considered, it is necessary to take into account the three components of the input motion.     

Evidence for fault-related directionality and localized site effects from strong motion recordings of the 2003 Boumerdes (Algeria) earthquake: Consequences on damage distribution and the Algerian seismic code

The Algiers–Boumerdes region has been struck by a destructive magnitude 6.8 (M_w) earthquake on May 21, 2003. The study presented in this paper is based on main shock strong motions from 13 stations of the Algerian accelerograph network. A maximum 0.58g peak ground acceleration (PGA) has been recorded at 20 km from the epicenter, only about 150 m away from a PGA of 0.34g, with both a central frequency around 5 Hz, explained by a strong very localized site effect, confirmed by receiver function technique results showing peaks at 5 Hz with amplitudes changing by a factor of 2. Soil amplifications are also evidenced at stations located in the quaternary Mitidja basin, explaining the higher PGA values recorded at these stations than at stations located on firm soil at similar distances from the epicenter. A fault-related directionality effect observed on the strong motion records and confirmed by the study of the seismic movement anisotropy, in agreement with the N65 fault plan direction, explains the SW–NE orientation of the main damage zone. In the near field, strong motions present a high-frequency content starting at 3 Hz with a central frequency around 8 Hz, while in the far field their central frequency is around 3 Hz, explaining the high level of damage in the 3- to 4-story buildings in the epicentral zone. The design spectra overestimate the recorded mean response spectra, and its high corner frequency is less than the recorded one, leading to a re-examination of the seismic design code that should definitively integrate site-related coefficient, to account for the up to now neglected site amplification, as well as a re-modeling of the actual design spectra. Finally, both the proposed Algerian attenuation law and the worldwide laws usually used in Algeria underestimate the recorded accelerations of the 6.8 (M_w) Boumerdes earthquake, clearly showing that it is not possible to extrapolate the proposed Algerian law to major earthquakes.     

Vibration Effect of Near Earthquakes at Different Depths in a Shallow Medieval Mine

The shallow medieval Jeroným Mine is located at a distance of about 25 km southeast of the Nový Kostel focal zone where the most intensive seismic activity in West Bohemia (Czech Republic) has been documented. Permanent seismological monitoring has been carried out since 2004 in this mine. During the 2011 and 2014 seismic swarms, more than 1000 triggered records comprising almost 1500 earthquakes were recorded at the permanent station in the mine. Three short-term seismological experiments were accomplished during these swarms. Several temporary seismic stations were simultaneously placed in different parts of underground spaces which enabled comparison of vibration effect caused by near earthquakes in different parts of the mine. Although the depth of the lowest parts of mine is only about 60 m, a vibration effect generated by earthquakes from the Nový Kostel focal zone is not the same for the whole underground complex.     

地震动加速度反应谱场地条件影响调整系数研究

场地条件对地震动特性影响显著,在抗震设计反应谱的确定过程中,需根据场地条件对加速度反应谱予以相应的调整。已有场地条件影响调整方案研究成果,均基于数值模拟或局部地区强震动记录统计,多数仅给出了峰值加速度PGA场地条件影响调整系数,对非线性的考虑缺乏观测数据依据。为此在全球强震动记录统计获得的PGA归一化加速度反应谱和日本钻井台阵记录获得的加速度反应谱平台值非线性衰减指数的基础上,结合钻孔模型数值模拟和近期研究成果,建立了考虑场地条件影响非线性的地震动加速度反应谱场地条件影响调整系数方案。     

A comparative analysis of seismic response of shallow buried underground structure under incident P,SV and Rayleigh waves

In this study, A time-domain seismic response analysis method and a calculation model of the underground structure that can realize the input of seismic P, SV and Rayleigh waves are established, based on the viscoelastic artificial boundary elements and the boundary substructure method for seismic wave input. After verifying the calculation accuracy, a comparative study on seismic response of a shallow-buried, double-deck, double-span subway station structure under incident P, SV and Rayleigh waves is conducted. The research results show that there are certain differences in the cross-sectional internal force distribution characteristics of underground structures under different types of seismic waves. The research results show that there are certain differences in the internal force distribution characteristics of underground structures under different types of seismic waves. At the bottom of the side wall, the top and bottom of the center pillar of the underground structure, the section bending moments of the underground structure under the incidences of SV wave and Rayleigh wave are relatively close, and are significantly larger than the calculation result under the incidence of P wave. At the center of the side wall and the top floor of the structure, the peak value of the cross-sectional internal force under the incident Rayleigh wave is larger than the calculation result under SV wave. In addition, the floor of the underground structure under Rayleigh waves vibrates in both the horizontal and vertical directions, and the magnification effect in the vertical direction is more significant. Considering that the current seismic research of underground structures mainly considers the effect of body waves such as the shear waves, sufficient attention should be paid to the incidence of Rayleigh waves in the future seismic design of shallow underground structures.     

Spectral parameters of ground motion in different regions: comparison of empirical models

Fourier-amplitude spectrum is one of the most important parameters describing earthquake ground motion, and it is widely used for strong ground motion prediction and seismic hazard estimation. The relationships between Fourier-acceleration spectra, earthquake magnitude and distance were analysed for different seismic regions (the Caucasus and Taiwan island) on the basis of ground motion recordings of small to moderate (3.5≤M_L≤6.5) earthquakes. It has been found that the acceleration spectra of the most significant part of the records, starting from S-wave arrival, can be modelled accurately by the Brune's “ω-squared” point-source model. Parameters of the model are found to be region-dependent. Peak ground accelerations and response spectra for condition of rock sites were calculated using stochastic simulation technique and obtained models of source spectra. The modelled ground-motion parameters are compared with those predicted by recent empirical attenuation relationship for California.     

Deterministic Approach for the Seismic Microzonation of Bucharest

— The mapping of the seismic ground motion in Bucharest, due to the strong Vrancea earthquakes is carried out using a complex hybrid waveform modeling method which combines the modal summation technique, valid for laterally homogeneous anelastic media, with finite-differences technique, and optimizes the advantages of both methods. For recent earthquakes, it is possible to validate the modeling by comparing the synthetic seismograms with the records. We consider for our computations the frequency range from 0.05 to 1.0 Hz and control the synthetic signals against the accelerograms of the Magurele station, low-pass filtered with a cut-off frequency of 1.0 Hz of the 3 last major strong (M_w > 6) Vrancea earthquakes. Using the hybrid method with a double-couple seismic source approximation, scaled for the source dimensions and relatively simple regional (bedrock) and local structure models, we succeeded in reproducing the recorded ground motion in Bucharest at a satisfactory level for seismic engineering. Extending the modeling to the entire territory of the Bucharest area, we construct a new seismic microzonation map, where five different zones are identified by their characteristic response spectra.     

A brief theory and computing of seismic ground rotations for structural analyses

The seismic ground rotations are important with respect to spatial structural models, which are sensitive to the wave propagation. The rotational ground motion can lead to significant increasing of structural response, instability and unusual damages of buildings. Currently, the seismic analyses often take into account the rocking and torsion motions separately using artificial accelerograms. We present an exact analytical method, proposed by Nazarov [15] for computing of three rotational accelerograms simultaneously from given translational records. The method is based on spectral representation in the form of Fourier amplitude spectra of seismic waves, corresponding to the given three-component translational accelerogram. The composition, directions and properties of seismic waves are previously determined in the form of a generalized wave model of ground motion. It is supposed that seismic ground motion can be composed by superposition of P, SV, SH- and surface waves. As an example, the dynamic response analysis of 25-story building is presented. Here recorded (low-frequency) and artificial (high-frequency) accelerograms were used; each of them includes three translational and three rotational components. In this structural analysis, we have clarified primarily conditions under which rotational ground motion should be taken into account. Next, we have calculated three rotational components of seismic ground motion. Then they were taken as additional seismic loads components for further seismic analysis of the building. Note, soil–structure interaction (SSI) is not considered in this study. For computing, we use the special software for structural analyses and accelerogram processing (FEA Software STARK ES and Odyssey software, Eurosoft Co., Russia). It was developed and is used in engineering practice in the Central Research Institute of Building Constructions (TsNIISK, Moscow, Russia).     

Effects of near-fault and far-fault ground motions on nonlinear dynamic response and seismic damage of concrete gravity dams

As the forward directivity and fling effect characteristics of the near-fault ground motions, seismic response of structures in the near field of a rupturing fault can be significantly different from those observed in the far field. The unique characteristics of the near-fault ground motions can cause considerable damage during an earthquake. This paper presents results of a study aimed at evaluating the near-fault and far-fault ground motion effects on nonlinear dynamic response and seismic damage of concrete gravity dams including dam-reservoir-foundation interaction. For this purpose, 10 as-recorded earthquake records which display ground motions with an apparent velocity pulse are selected to represent the near-fault ground motion characteristics. The earthquake ground motions recorded at the same site from other events that the epicenter far away from the site are employed as the far-fault ground motions. The Koyna gravity dam, which is selected as a numerical application, is subjected to a set of as-recorded near-fault and far-fault strong ground motion records. The Concrete Damaged Plasticity (CDP) model including the strain hardening or softening behavior is employed in nonlinear analysis. Nonlinear dynamic response and seismic damage analyses of the selected concrete dam subjected to both near-fault and far-fault ground motions are performed. Both local and global damage indices are established as the response parameters. The results obtained from the analyses of the dam subjected to each fault effect are compared with each other. It is seen from the analysis results that the near-fault ground motions, which have significant influence on the dynamic response of dam–reservoir–foundation systems, have the potential to cause more severe damage to the dam body than far-fault ground motions.     

Regional and local seismic hazard assessment

Recent earthquakes such as the M_JMA 7.2 Hyogo-ken Nambu earthquake and the M 7.4 Kocaeli earthquake demonstrate once again the need to include detailed soil investigation into hazard evaluation, that is the need of microzonation. Seismic hazard assessment evaluated at a regional scale generally does not consider soil effects but only in a limited way using an attenuation law that can be ‘soft soil’ or ‘rock’. However, the relevant role of seismic hazard in the assessment of seismic coefficients for the definition of the actions in seismic codes must be properly considered. That is to say, the level of protection of buildings is proportional to a definite level of hazard (generally considered to be the ground motion with 10% probability of exceedence in 50 years). When a microzonation is performed, this criterion cannot be ignored, therefore, a clear linkage must be established between hazard (regional scale) and microzonation. The crucial point is represented by the reference motion (or input motion) to be used for site effects analysis, that must be compatible with the regional seismic hazard. In this paper, three different approaches for reference motion evaluation are analysed: probabilistic; stochastic; and deterministic. Through the case history of Fabriano microzonation the three approaches are compared. It is shown that each approach presents advantages and disadvantages with respect to the others. For example, the probabilistic approach (the reference motion is directly derived from the expected response spectra for a given return period) is linked with hazard, but produces an overestimation in short periods range, while the deterministic approach correctly simulates the wave propagation, but it ends with a kind of conditional probability. Until now, clear criteria to choose the right approach do not appear to exist and the expert experience is of fundamental importance.     

Critical evaluation of strong motion in Kocaeli and Düzce (Turkey) earthquakes

Turkey was struck by two major events on August 17th and November 12th, 1999. Named Kocaeli (M_w=7.4) and Düzce (M_w=7.2) earthquakes, respectively, the two earthquakes provided the most extensive strong ground motion data set ever recorded in Turkey. The strong motion stations operated by the General Directorate of Disaster Affairs, the Kandilli Observatory and Earthquake Research Institute of Bogazici University and Istanbul Technical University have produced at least 27 strong motion records for the Kocaeli earthquake within 200 km of the fault. Kocaeli earthquake has generated six motions within 20 km of the fault adding significantly to the near-field database of ground motions for M_w>=7.0 strike–slip earthquakes. The paper discusses available strong motion data, studies their attenuation characteristics, analyses time domain, as well as spectral properties such as spectral accelerations with special emphasis on fault normal and fault parallel components and the elastic attenuation parameter, kappa. A simulation of the Kocaeli earthquake using code FINSIM is also presented.     

Site response study in Abruzzo (Central Italy): underground array versus surface stations

In this site response study we examined local earthquakes recorded at surface stations of a local seismic network and at a temporary underground seismic array installed in a tunnel underneath the Gran Sasso Massif in Abruzzo (central Italy). This allowed us to compare the seismic site response beneath the mountain and on the surface in similar geological environment (soft rock sites). We applied spectral ratios method on different segments of the seismograms and used different reference spectra in the 1–20 Hz frequency band. We found little or no amplification effects at most of the surface stations whereas site transfer functions evaluated with respect to underground sites show an amplification factor up to 6 in the 1–8 Hz frequency range. Coda spectral ratios estimated at soft rock sites are confirmed as good estimates of shear wave transfer function.     

Assessment of liquefaction potential based on peak ground motion parameters

Conventionally, evaluation of liquefaction potential of loose saturated cohesionless deposits as specified in Japanese design codes employs peak ground acceleration (PGA). However, recent large-scale earthquakes in Japan revealed that liquefaction at some sites did not occur even though large PGAs were recorded at or near these sites. As an alternative approach, an evaluation procedure based on peak ground motion parameters, i.e. incorporating both PGA and the peak ground velocity (PGV), is proposed. By performing parametric studies using one-dimensional seismic response analysis and formulating regression models, seismic-induced shear stresses within the deposit are expressed in terms of peak ground motion parameters at the surface, and these are used to calculate the factor of safety against liquefaction. Application to case histories in Japan indicates that the proposed two-parameter equation can adequately account for the occurrence and non-occurrence of liquefaction at various sites as compared to the conventional PGA-based approach. Moreover, analyses of several strong motion records at various sites show that liquefaction may occur when PGA≥150 gal and PGV≥20 kine, indicating that these values can serve as thresholds in assessing the possible occurrence of liquefaction.     

汶川地震远场地震动特征及其对长周期结构影响的分析

汶川8．0级大地震中,中国数字强震动台网获得了大量的数字强震动记录,这些记录特别是远场记录具有丰富的长周期地震动分量。本文根据东南强震动中心获取的区域数字强震动记录,分析汶川8．0级地震的远场地震动特征,研究了基岩场地及深厚软弱场地的长周期地震动及其差异;根据长周期结构的特点,选取建于深厚软弱场地上的江苏A050强震台的超长地震记录进行结构地震反应分析,研究汶川8．0级地震对远场长周期结构的影响,并结合实际震害特点,提出了长周期结构抗震及地震安全对策中一些值得注意和思考的问题。     

用断裂动力学制定地震烈度区划初探

根据第一作多年来从事用断裂动力学的解析解来表示强震地面运动的振幅、速度、加速度,所得的结果与两次名强震地面运动的记录进行了比较,二符合较好,故提出应用断裂动力学的研究成果来制定确定性的地震烈度区划方法的初步探讨意见。