土层平均剪切波速确定及其可靠性验证——基于钻井台阵强震动记录

震害资料显示,场地条件对地震动特性以及工程结构破坏程度影响显著。为减少因场地效应而造成的经济损失和社会影响,在进行场地地震反应分析时,需最大限度地减小因场地土层模型参数的不确定性引起的地震动评估偏差,为工程结构地震反应分析选取并生成适当的地震动输入。随着强震动观测技术的逐渐发展,大量可靠的钻井台阵记录为地震过程中场地观测点的动力反应提供了直接数据。以美国加州地区La Cienega钻井台阵强震动观测数据为基础,利用互相关函数,对不同强度地震作用下场地土层的平均剪切波速进行分析,并在此基础上,以Cyclic 1D为模拟平台,建立一维自由场地地震反应有限元分析模型。分析结果表明：通过钻井台阵地震动观测数据识别,得到场地平均剪切波速,能够反映该场地的动力特性,数值模拟计算结果和台阵地震动记录基本吻合,可为数值模型参数选取提供依据。     

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

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

Comparison of 1D linear,equivalent-linear,and nonlinear site response models at six KiK-net validation sites

Vertical seismometer arrays represent a unique interaction between observed and predicted ground motions, and they are especially helpful for validating and comparing site response models. In this study, we perform comprehensive linear, equivalent-linear, and nonlinear site response analyses of 191 ground motions recorded at six validation sites in the Kiban–Kyoshin network (KiK-net) of vertical seismometer arrays in Japan. These sites, which span a range of geologic conditions, are selected because they meet the basic assumptions of one-dimensional (1D) wave propagation, and are therefore ideal for validating and calibrating 1D nonlinear soil models. We employ the equivalent-linear site response program SHAKE, the nonlinear site response program DEEPSOIL, and a nonlinear site response overlay model within the general finite element program Abaqus/Explicit. Using the results from this broad range of ground motions, we quantify the uncertainties of the alternative site response models, measure the strain levels at which the models break down, and provide general recommendations for performing site response analyses. Specifically, we find that at peak shear strains from 0.01% to 0.1%, linear site response models fail to accurately predict short-period ground motions; equivalent-linear and nonlinear models offer a significant improvement at strains beyond this level, with nonlinear models exhibiting a slight improvement over equivalent-linear models at strains greater than approximately 0.05%.     

Response of the Tarzana strong motion site during the 1994 Northridge earthquake

Exceptionally high ground motions (horizontal peak ground acceleration (PGA) of 1.82g) were recorded at the Tarzana Station during the main shock of the 1994 Northridge earthquake (moment magnitude 6.7 at an epicentral distance of 6 km). At the time of the main shock, the instrument was located near the edge of a 21 m-high ridge with side slopes ranging from 3H:1V to 15H:1V. The ridge is underlain by shallow fill and soft rocks of Medelo Formation.

The objectives of this study were to (1) identify the relative contributions of various factors such as local geology, topography, source mechanism, and travel path on the large ground motions recorded at Tarzana Station and (2) develop an analytical model that could adequately predict observed ground motions at the Tarzana site during the Northridge earthquake and at similar sites during future earthquakes. This study is an integral part of a series of inter-related studies referred to as the ROSRINE research (Resolution of Site Response Issues during Northridge Earthquake) project.

The PGA at the surface of competent bedrock (1 km/s shear wave velocity found about 100 m below ground surface) is estimated by Silva [ROSRINE Study (2000)] at 0.46 gravity (g). To identify the source of ground motion amplification, one-dimensional ( ), two-dimensional (TELDYN and SASSI), and three-dimensional (SASSI) analyses were conducted using both recorded aftershock data and an estimated ground acceleration time histories at a 100 m depth.

The results of the analyses indicate that (1) local geology and topography could only partially account for the observed ground motion amplification, and (2) the PGA and response spectra at a point near the edge of the ridge (the location of the instrument at the time of the main shock) is in good agreement with recorded values when the angle of incident of shear waves (SV waves) at 100 m depth is assumed at 30° from vertical. Considering the local geology and variation of shear wave velocity with depth, the 30° incident angle at 100 m depth corresponds to an 8° incident angle of shear waves at the ground surface. This observation is, in general, consistent with the incident angles of shear waves reported from study of the recorded aftershock data.     

Site dependence of far-source ground motions during the Wenchuan earthquake

This paper aimed to examine the site dependence and evaluate the methods for site analysis of far-source ground motions. This was achieved through the examination of frequency content estimated by different methods based on strong ground motions recorded at twelve far-source stations in Shandong province during the Wenchuan earthquake. The stations were located in sites with soil profiles ranging from code classes Ⅰ to Ⅲ. Approaches used included the Fourier amplitude spectrum (FAS), the earthquake response spectrum (ERS), the spectral ratio between the horizontal and the vertical components (H/V), the spectral ratio between the spectra at the site and at a reference site (SRRS), and coda wave analysis (CWA). Results showed that major periods of these ground motions obtained by FAS, ERS and H/V ratio methods were all evidently larger than site dominant periods; the periods were also different from each other and mainly reflected the frequency content of long period components. Prominent periods obtained by the SRRS approach neither illuminated the long period aspect nor efficiently determined site features of the motions. The CWA resulted in a period close to site period for stations with good quality recordings. The results obtained in this study will be useful for the evaluation of far-source effect in constructing seismic design spectra and in selecting methods for ground motion site analysis.     

Anisotropic seismic ground responses identified from the Hualien vertical array

This paper presents the results of ground response analyses using the free-field motions of six earthquakes recorded at the site of the Hualien Large Scale Seismic Test. From the phase spectra of transfer functions between the ground motions recorded at the surface and various depths, the predominant frequencies of soil layers are identified and used to estimate the shear wave velocities of soil layers. Results show that the soils at the Hualien LSST site are anisotropic, such that the earthquake responses exhibit different characteristics in two horizontal principal directions.     

Site response in the Qionghai Basin in the Wenchuan earthquake

Amplification effects of soil site response can significantly impact ground motions, and must be considered in the seismic fortification of buildings/structures to prevent or mitigate this potential seismic hazard. Utilizing acceleration time histories from the main shock of the Wenchuan earthquake recorded at four stations (i.e., one on bedrock and three on soil) in the Qionghai Basin, the site responses from three soil sites are studied by using the traditional spectral ratio method. The bedrock site is selected as a reference site. This study found that peak ground accelerations (PGAs) on the soil sites are much larger than on bedrock, with EW, NS and UD components of 3.96–6.58, 6.27–10.98, and 3.17–6.66 times those of the bedrock site, respectively. The amplification effects of the soil sites on ground motions in the frequency range of 0.1 Hz to 10 Hz are significant, depending on the thickness of the soil layer and the frequency content of the site. A significant amplification occurs with high frequency components of ground motion at shallow soil sites, and low and high frequency components of ground motion at intermediate soil sites.     

Topographic proxies from 2-D numerical analyses

We use 2D numerical analyses at ground motion stations from the NGA-West2 dataset to develop parameters to predict the effect of surface topography on response spectra. The simplistic numerical analyses use sinusoidal input motions, uniform soil profiles, elastic soil properties, and absorbing boundary conditions. We define several topographic parameters for stations using the natural logarithm of 2D amplifications in peak ground acceleration of a sinusoidal input motion in different orientations. The natural log of 2D amplifications when averaged over multiple orientations is found to have the most predictive power among the studied parameters. We also explore the relationship between the topographic parameters developed in this study, and the topographic parameters developed at the same sites in an earlier study (Rai et al. in Earthq Spectra, 2016b).     

一种深厚覆盖层场地地震动参数的确定方法

基于欧美规范确定了坐落在深厚覆盖层上KH抽水蓄能电站上、下库场地基本运行和最大设计地震动峰值加速度、反应谱和时程等动参数。首先依据场地区域地震烈度区划图、特征周期区划图和依据场地地质地震条件选取的5条种子实测地震动确定场地基岩输入加速度时程、峰值加速度和设计反应谱,进而基于各土层地质参数和一维弹性波传播模拟程序确定覆盖层表面的平均峰值加速度、平均反应谱和5条地震动时程,对所得到的平均反应谱和峰值加速度进行光滑处理后确定可用于各建筑物结构抗震设计的地震动参数,包括覆盖层表面水平向动力响应加速度时程、峰值加速度和设计反应谱。该方法可较好地保留输入地震动的真实动力特性,如持时、相位和频率等,为我国规范中建议的确定场地地震动参数的方法提供有益的补充。     

Validation of the Deep Velocity Structure of the Tokachi Basin Based on 3-D Simulation of Long-Period Ground Motions

The construction of 3-D basin velocity structures is ongoing in many regions of Japan. The structure models are constructed mainly for the prediction of long-period ground motions from future large earthquakes. In this paper, we validate the 3-D velocity structure model of the Tokachi basin, a deep sedimentary basin located in eastern Hokkaido, Japan, based on 3-D simulation of long-period (2–20 s) ground motions from three nearby intermediate-depth earthquakes; this model was constructed by the National Research Institute for Earth Science and Disaster Prevention (NIED). We make comparisons between the observed and synthetic long-period ground motions for the basin-induced surface waves as well as the direct S-wave. We also try to revise the 3-D velocity structure in the western part of the Tokachi basin based on 1-D velocity structures estimated using long-period S-wave modeling and the microtremor survey method. We then perform the 3-D simulation again to validate the revised model. Based on quantitative comparisons of the long-period ground motions from these simulations with those observed, we conclude that the NIED and revised velocity structure models are generally good at the central basin sites, but that both models require modification at the basin edges to explain the details of the observed basin-induced surface waves.     

覆盖层厚度对反应谱峰值的影响

在研究多个工程场地钻孔资料的基础上,选取和构造了若干典型场地剖面,利用目前工程上广泛应用的场地地震反应分析的一维等效线性化波动方法,计算了各剖面在不同基底地震动输入的反应谱值及地表加速度峰值,分析了覆盖层厚度对反应谱峰值及峰值周期、地表加速度峰值和放大倍数的影响,得出了一些有意义的结论。     

Analysis of ground motions at Treasure Island site during the 1989 Loma Prieta earthquake

Site response analyses and coherence studies were conducted at the Treasure Island site where surface motions were recorded during the Loma Prieta earthquake. The analyses were conducted using a nonlinear dynamic effective stress method which took into account the effects of the liquefaction that occured at the site. The rock motions recorded at nearby Yerba Buena Island were used as input motions. Computed and recorded ground motions transverse to the direction of wave propagation and associated response spectra were in good agreement. Agreement was also good in the radial direction, except in certain frequency bands higher than 1·25 Hz. Coherence studies showed that some of these discrepancies may be due to low coherence between the Treasure Island and Yerba Buena motions in these same frequency bands.     

Modeling of spatially correlated,site-reflected,and nonstationary ground motions compatible with response spectrum

Seismic risk analysis and mitigation of spatially extended structures require the synthesis of spatially varying ground motions in the response history analysis of these structures. These synthetic motions are usually desired to be spatially correlated, site reflected, nonstationary, and compatible with target design response spectra. In this paper, a method is presented for simulating spatially varying ground motions considering the nonstationarity, local site effects, and compatibility of response spectra. The scheme for generating spatially varying and response spectra compatible ground motions is first established for spatial locations on the ground surface with varying site conditions. The design response spectrum is introduced as the “power” spectrum at the base rock. The site amplification approach is then derived based on the deterministic wave propagation theory, by assuming that the base rock motions consist of out-of-plane SH wave or in-plane combined P and SV waves propagating into the site with assumed incident angles, from which tri-directional spatial ground motions can be generated. The phase difference spectrum is employed to model ground motions exhibiting nonstationarity in both frequency and time domains with different site conditions. The proposed scheme is demonstrated with numerical examples.     

Evaluating nonlinear effective stress site response analyses using records from the Canterbury earthquake sequence

A widely used one-dimensional nonlinear effective stress site response analysis program is used to model the response of potentially liquefiable soils during strong shaking. Ground motion records from six events of the 2010–2011 Canterbury earthquake sequence and the extensive site investigation data that have been obtained for the Christchurch area provide the basis for the analyses. The results of the analyses depend significantly on the input motions and soil profile characterization, so these important aspects are examined. Deconvolved Riccarton Gravel input motions were generated, because recorded rock or firm layer motions were not available. Nonlinear effective stress seismic site response analyses are shown to capture key aspects of the observed soil response through the comparison of acceleration response spectra of calculated surface motions to those of recorded surface motions; however, equivalent-linear and total stress nonlinear analyses capture these aspects as well. Biases in the computed motions compared to recorded motions were realized for some cases but they can be attributed primarily to the uncertainty in the development of the input motions used in the analyses.     

Near-fault ground motion bi-normalized pseudo-velocity spectra and its applications

Ground motions with forward-directivity effect in the near-fault region are obviously different from ordinary far-field ground motions. Design spectral models for this kind of motions have been proposed by correlating sim-ple pulses with parameters attenuation relationships in a previous study of the authors. To further test the applica-bility of the established design spectral model, we analyze ground motion pseudo-velocity response spectra (PVS), normalized pseudo-velocity spectra (NPVS) and bi-normalized pseudo-velocity spectra (BNPVS) of 53 typical near-fault forward-directivity ground motions. It is found that BNPVS not only has more salient features to reflect the difference between soil and rock sites, but also has less scattering to reveal the nature of forward-directivity motions. And then, BNPVS is used for prediction of design spectra accounting for the influence of site conditions, and the constructed design spectra are compared with those spectra established previously. It is concluded that site condition can heavily affect ground motions, buildings on rock can be even more dangerous than those on soil sites, in particular for ordinary buildings with short to middle vibration periods. Finally, pulse models are also suggested for structural analyses in the near-fault region.     

渭河盆地中土层场地对地震动的放大作用

土层场地对地震动的影响较大,建(构)筑物的选址及其抗震设防必须考虑土层场地的放大作用,以避免或减轻其震害.汶川地震中,布设在渭河盆地中的数字强震动台网共有27个台站(包括2个基岩台站和25个土层台站)获得良好的主震加速度时程.利用这些加速度时程,选择汤峪台做为参考场地,基于考虑几何衰减的传统谱比法分析研究了25个土层场...     

Equivalent stationary motion and average response spectra

This study not only attempts to justify an equivalent stationary motion model, but also endeavours to illustrate its potentials as an attractive analysis tool. Using 38 ground accelerations recorded on rock sites as data base, this study first demonstrates that the equivalent stationary motions do give average response spectra compatible with real records. This study also shows that through the equivalent stationary model, a set of motions may be defined by specifying the joint distributions of the model parameters. Finally, by converting a relationship between earthquake magnitude and liquefaction potentials of ground motions, this study has derived normalized magnitude-dependent response spectra for rock sites.     

三种土层结构反应谱特征周期的统计分析

本文人工合成了若干条能够反映不同地震动特征的基岩加速度时程,并将其作为土层地震反应分析的地震动输入。选取了若干个软弱夹层分别在底部、中部和顶部的场地剖面,利用工程上广泛应用的场地地震反应分析的一维等效线性化波动方法,计算了不同场地类别,3种土层结构在不同地震动输入下的地表加速度反应谱的特征周期。给出了不同土层结构在统计意义上的反应谱特征周期的平均值,并通过与正常剖面反应谱特征周期的比较给出了不同场地类别,3种土层结构的反应谱特征周期的影响系数。本文的研究结果对进一步考虑场地分类具有一定的意义。     

Simulation of synthetic ground motions for specified earthquake and site characteristics

A method for generating a suite of synthetic ground motion time‐histories for specified earthquake and site characteristics defining a design scenario is presented. The method employs a parameterized stochastic model that is based on a modulated, filtered white‐noise process. The model parameters characterize the evolving intensity, predominant frequency, and bandwidth of the acceleration time‐history, and can be identified by matching the statistics of the model to the statistics of a target‐recorded accelerogram. Sample ‘observations’ of the parameters are obtained by fitting the model to a subset of the NGA database for far‐field strong ground motion records on firm ground. Using this sample, predictive equations are developed for the model parameters in terms of the faulting mechanism, earthquake magnitude, source‐to‐site distance, and the site shear‐wave velocity. For any specified set of these earthquake and site characteristics, sets of the model parameters are generated, which are in turn used in the stochastic model to generate the ensemble of synthetic ground motions. The resulting synthetic acceleration as well as corresponding velocity and displacement time‐histories capture the main features of real earthquake ground motions, including the intensity, duration, spectral content, and peak values. Furthermore, the statistics of their resulting elastic response spectra closely agree with both the median and the variability of response spectra of recorded ground motions, as reflected in the existing prediction equations based on the NGA database. The proposed method can be used in seismic design and analysis in conjunction with or instead of recorded ground motions. Copyright © 2010 John Wiley & Sons, Ltd.     

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.