用改进的经验格林函数方法模拟唐山地震动

建筑物的抗震设防需要尽可能地掌握未来大地震强震动记录信息,但大地震强震动记录的匮乏阻碍了抗震设防实践的发展。经验格林函数方法作为模拟地震动的主要方法,可以提供可靠的大地震强震动记录,但也存在着许多问题,如缺乏对大地震断层滑动分布不均匀的描述、用经验确定小震数目、模拟方法受到大小地震相似条件的限制等。文中对上述经验格林函数方法存在的问题进行了研究,改进的经验格林函数方法,有效地解决了上述问题。并用其对唐山大地震进行了模拟,并把模拟的地震动时程和反应谱与实际记录相比较,发现用改进方法模拟的地震动加速度反应谱比用未改进方法模拟结果更接近实际的地震动记录加速度反应谱。由此说明改进的经验格林函数可更准确的模拟地震动。     

HORIZONTAL-TO-VERTICAL SPECTRUM RATIO OF EARTHQUAKE GROUND MOTION FOR SITE CHARACTERIZATION

Nakamura's method, which uses a horizontal-to-vertical Fourier spectrum ratio of microtremor, has become popular to determine the predominant period and amplification of a site. In this study, this method is extended for earthquake ground motion recordings using new strong motion data recorded by JMA-87-type accelerometers. From the analysis of these accelerograms, horizontal-to-vertical Fourier spectrum ratios of a site for different earthquakes are also found to be stable irrespective of magnitude, distance and depth. To establish this fact, attenuation relations of velocity response spectra for horizontal and vertical components are derived for three damping ratios (0, 2 and 5 per cent) using the JMA data. Then the horizontal-to-vertical ratios of the velocity response spectra are obtained. The results show that the horizontal and vertical velocity response spectra are dependent on magnitude, distance and depth, but that their ratios are almost independent of magnitude, distance and depth. However, since the current data set consists of mostly intermediate to far field data, this observation should be limited to records of these distance ranges. Introducing station coefficients, representing site amplification, to this relation yields the value comparable to the horizontal-to-vertical Fourier spectrum ratio at a specific site. The stability of the spectrum ratio is explained by the transfer function between the ground surface and stiff-soil outcrop due to S-wave propagation. These results suggest that site amplification characteristics can be evaluated by one-point two-component surface recordings of earthquake ground motion, in a similar manner as proposed by Nakamura for microtremor. © 1997 by John Wiley & Sons, Ltd.     

Analysis of nonlinear site response using the LSST downhole accelerometer array data

Both linear and nonlinear behaviors of soil deposits were evaluated by strong and weak motion data observed on the surface and at depths of 6, 11, 17, 47 m at the Large Scale Seismic Test (LSST) array in Lotung, Taiwan. The soil properties measured by well logging and by the shear wave velocity profile measured by uphole and cross-hole methods are available. Both one-dimensional equivalent-linear method and nonlinear method are used for the evaluation have been used. The synthetic records at various depths are obtained by using the records at the bottom as input motion. These synthetic records are then compared with actual records at corresponding depths. Records of 13 earthquakes are used. We find that the synthetic records obtained from a linear model match well with actual records for small input motions, but the results obtained from a nonlinear model match poorly. On the other hand, the synthetic records using both the nonlinear model and equivalent-linear model are in good agreement with the observed records for large input motions. In these cases, the predicted response spectra using the linear model consistently overestimate the observed records. The threshold distinguishing the large and small input motions is 0.04 g at depth of 47 m for the LSST data. Thus, the nonlinearity started at 0.04 g and occurred unequivocally at 0.075 g. Furthermore, the dominant frequencies shift toward lower values when input motions become large. Clearly, the observed records at the LSST site manifest nonlinearity of soil response. The hysteresis loops evaluated by the nonlinear method show a permanent strain of about 0.01% in soil layers at higher ground motion input levels in this case.     

Applicability of the H/V spectral ratio of microtremors in assessing site effects on seismic motion

The authors examine the reliability of site response estimations obtained by the horizontal to vertical (H/V) spectral ratios of microtremors by means of cross‐validation with the ratio of the horizontal spectra of earthquake motion with respect to reference sites. The data comprise microtremor and ground motion records recorded at 150 sites of Yokohama strong motion array. The use of non‐supervised pattern recognition techniques aims to group the sites with more objectivity. Attributes defining the overall shape of the amplification spectra serve as input in the computation of Euclidean distance similarity coefficients amongst sites. The implementation of the Ward clustering scheme leads to the attainment of a meaningful tree diagram. Its analysis shows the possibility of summarizing the results into six general patterns. A good coincidence of site effects estimates at 80 per cent of the sites becomes apparent. However, this coincidence appears poor for sites characterized by H/V amplification ratios around 2 or smaller and predominant periods longer than 0.5 s. In such cases, the presence of stiff, sandy sediments in the soil profile proves common. To proscribe H/V estimations, relying solely on the small spectral ratios criterion seems inadequate. Copyright © 2001 John Wiley & Sons, Ltd.     

基于NGA模型的芦山余震地震动衰减关系

四川省芦山县于2013年4月20日发生了7.0级地震,基于国家强震动台网中心发布的此次地震的2 106组余震强震动记录,利用基于NGA衰减关系的简化模型,使用随机效应方法拟合得到其PGA以及0.01~8 s周期下的伪加速度反应谱（阻尼比5%）的地震动衰减关系。与其他学者研究得出的四川其他地区衰减关系进行对比,本文模型能够较好地预测芦山余震强震动衰减。不论基岩场地和软土场地,本文在衰减关系在短周期均与卢大伟模型结果一致。但是在长周期情况下,本文模型在软土场地大于卢大伟结果,而在基岩场地上低于俞言祥结果。因此,在使用本文衰减关系时,要特别注意场地和周期对其的影响。     

Simulation of accelerograms from simplified deterministic approach for the 23rd October 2004 Niigata-ken Chuetsu,Japan earthquake

A simple hybrid approach for the simulation of strong ground motion is presented in this paper. This approach is based on the deterministic modelling of rupture plane initially started by Midorikawa, Tectonophysics 218:287–295, (1993) and further modified by Joshi, Pure Appl Geophys (PAGEOPH) 8:161, (2004). In this technique, the finite rupture plane of the target event is divided into several subfaults, which satisfy scaling relationship. In this paper, simulation of strong ground motion due to a rupture buried in a earth medium consisting of several layers of different velocities and thicknesses is made by considering (1) transmission of energy at each layer; (2) frequency filtering properties of medium and earthquake source; (3) correction factor for slip of large and small magnitude earthquakes and (4) site amplification ratio at various stations. To test the efficacy of the developed technique, strong motion records were simulated at different stations that have recorded the 2004 Niigata-ken Chuetsu, Japan earthquake (M _s 7.0). Comparison is made between the simulated and observed velocity and acceleration records and their response spectra. Distribution of peak ground acceleration, velocity and displacement surrounding the rupture plane is prepared from simulated and observed records and are compared with each other. The comparison of synthetic with the observed records over wide range of frequencies shows that the present technique is effective to predict various strong motion parameters from simple deterministic model which is based on simple regression relations and modelling parameters.     

2012年唐山4.8级地震强震动记录特征及场地相关特性分析

2012年5月28日河北省唐山市古冶区与滦县交界发生4.8级地震, 国家强震动台网中心在河北、 天津和北京的94个强震动台站记录到了本次地震的加速度。 本文给出了获取记录的强震动台站分布及强震动记录结果, 统计了强震动记录数量随震中距的变化, 给出了3个较小震中距台站记录到的加速度时程; 绘制了空间地震动峰值加速度等值线图及周期0.2 s、 2.0 s加速度反应谱值的等值线, 发现峰值加速度等值线与长周期加速度反应谱等值线极值分布具有明显地域差异, 分析认为是由于厚沉积层对长周期地震动具有放大作用造成的。 通过强震动记录与适用于本区的三个衰减关系对比, 分析了此次地震的峰值加速度衰减特征, 同时研究了周期0.2 s、 2.0 s加速度反应谱值的衰减特征, 周期2.0 s反应谱值随震中距的衰减与衰减关系能较好地对应, 然而在震中距100～130 km沉积层较厚的集中地区, 表现出了实际记录较衰减关系值偏大的现象, 认为同样是由于厚沉积层对地震动加速度反应谱长周期的放大作用导致的。 研究了震中距差别不大的情况下, 场地类型与沉积层厚度对反应谱特征周期的影响, 对比基岩台站与软弱地基土层台站的强震动记录反应谱, 发现软弱土层台站的土层对地震动有一定的放大作用, 导致中长周期地震动被放大, 对比位于沉积层较薄的隆起区台站与位于沉积层较厚的凹陷区台站强震动记录反应谱, 发现厚的沉积层不仅对反应谱长周期有放大的作用, 同时也会使得反应谱特征周期值变大。     

汶川地震不同场地反应谱平台值统计分析

加速度反应谱平台值表征地震动的强度特性,场地条件是影响反应谱平台值的一个重要因素.本文选取6组同一地区断层距相近而场地条件不同的强震记录,对其标准化的加速度反应谱形状及平台值进行分析;并以汶川地震中173个有详细地勘场地上的强震记录为基础,统计分析了不同场地类别和断层距区间内的加速度反应谱平台值.本文研究结果显示,场地条件对加速度反应谱平台值有较大影响,随着场地变软,加速度反应谱平台值增大.本文定义了场地影响系数,计算并给出了不同地面峰值加速度对应的场地影响系数.     

Evaluation of seismic site factors in the Mississippi Embayment. II. Probabilistic seismic hazard analysis with nonlinear site effects

An integrated probabilistic seismic hazard analysis procedure that incorporates nonlinear site effects, PSHA-NL, is developed and used to characterize the influence of thick deposits of the upper Mississippi Embayment (ME) on seismic site coefficients. PSHA-NL follows the methodology of the 2002 USGS hazard maps and generates a compatible set of ground motion records. The motions are propagated using nonlinear and equivalent linear site response analyses and ME properties developed in a companion paper and used to derive surface uniform hazard response spectra. A set of generic site coefficients are derived and summarized in a format similar to NEHRP site coefficients, with an added dimension of ME deposits thickness to the Paleozoic rock, a physically meaningful impedance boundary. These coefficients compare well with NEHRP site coefficients for 30 m profiles. For thicker soil profiles, developed site coefficients are lower at short periods and higher at long periods than NEHRP site coefficients.     

Using satellite altimetry data to augment flow estimation techniques on the Mekong River

Satellite altimetry is routinely used to provide levels for oceans or large inland water bodies from space. By utilizing retracking schemes specially designed for inland waters, meaningful river stages can also be recovered when standard techniques fail. Utilizing retracked waveforms from ERS‐2 and ENVISAT along the Mekong, comparisons against observed stage measurements show that the altimetric measurements have a root mean square error (RMSE) of 0·44–0·65 m for ENVISAT and 0·46–0·76 m for ERS‐2. For many applications, however, stage is insufficient because discharge is the primary requirement. Investigations were therefore undertaken to estimate discharges at a downstream site (Nakhon Phanom (NP)) assuming that in situ data are available at a site 400 km upstream (Vientiane). Two hypothetical, but realistic scenarios were considered. Firstly, that NP was the site of a de‐commissioned gauge and secondly, that the site has never been gauged. Using both scenarios, predictions were made for the daily discharge using methods with and without altimetric stage data. In the first scenario using a linear regression approach the altimetry data improved the Nash‐Sutcliffe r² value from 0·884 to 0·935. The second scenario used known river cross‐sections while lateral inflows were inferred from a hydrological model: this scenario gave an increase in the r² value from 0·823 to 0·893. The use of altimetric stage data is shown to improve estimated discharges and further applications are discussed. Copyright © 2010 John Wiley & Sons, Ltd.     

离散Fréchet距离在强震台站场地分类中的应用

当强震台站场地资料不完整时,所收集到的强震数据因缺乏准确的场地类别信息而难以有效利用。为解决这一问题,本文提出一种基于离散Fréchet距离的强震台站场地分类方法。将获取到的664个KiK-net台站场地按照《建筑抗震设计规范（GB 50011—2010）》进行分类,并构建2个数据集。利用数据集1得到Ⅰ、Ⅱ、Ⅲ类场地标准谱比曲线,并结合离散Fréchet距离对数据集2中的台站进行场地分类,统计分类成功率与误判率。统计结果表明,本文方法能较准确地对场地进行分类,且误判率在可接受范围内。将本文方法分类统计结果与斯皮尔曼秩相关系数法分类的成功率与误判率进行对比,结合本文方法分类后得到的平均谱比曲线,均可表明应用本文方法进行强震台站场地分类具有合理性。     

A note on the numerical evaluation of Duhamel's integral

General expressions are derived for the numerical evaluation of Duhamel's integral and its derivative. The work comprises an extension (to unequal time steps) and an application (to a piecewise linear forcing function) of the numerical integration approach adopted by Cronin.¹ The application is particularly relevant to the digital computation of response spectra from strong motion earthquake records.     

汶川地震余震震源参数及地震动衰减与场地影响参数反演分析

基于强震动观测台钻孔数据, 统计得到了四川、 甘肃地区20 m深的浅硬土层场地平均剪切波速模型． 通过与美国加州地区的波速模型对比, 结合四川、 甘肃地区地震预报的地壳模型, 延拓剪切波速模型至40 km． 应用四分之一波长法计算了浅硬土层场地的平均场地放大系数． 利用场地放大系数, 消除13次余震中浅硬场地的台站场地响应, 反演了龙门山断层上、 下盘的介质品质因子以及13次余震的有效应力降． 通过与普通土层场地的强震动记录对比, 对浅硬土层场地的放大系数进行调整, 得到了深厚土层场地的平均放大系数．      

2008年汶川MS8.0地震反应谱平台值衰减特性研究

2008年汶川M_S8.0地震获得了大量的强震动记录, 为研究反应谱特征参数衰减特性提供了重要的基础资料. 本文对已知场地条件的174个台站的记录进行研究, 计算其加速度反应谱并按照最小二乘分段拟合方法进行标定, 进而拟合了反应谱平台值的衰减关系;对比分析了水平方向与竖直方向反应谱平台值的衰减特性, 提出用“平台值的场地衰减影响系数”来定量研究不同场地类型对反应谱平台值的影响, 用“衰减曲线下降速率”来定量分析不同衰减曲线的衰减速率. 通过计算得出水平向Ⅰ, Ⅱ, Ⅲ类场地的平台值场地衰减影响系数平均值分别为0.5358, 1和1.579, 且Ⅲ类场地的加速度平台值衰减曲线的衰减速率最小.      

Earthquake Induced Site Effect in the Algiers–Boumerdes Region: Relation Between Spectral Ratios Higher Peaks and Observed Damage During the May 21st M w 6.8 Boumerdes Earthquake (Algeria)

The Algiers–Boumerdes regions were hit by an M _w 6.8 destructive earthquake on May 21, 2003. The accelerometric and seismometric networks successfully recorded the main shock and many aftershocks at some locations where the damage was most extensive. A microtremor measurement was performed in the same locations; some of them are localized on the Mitidja basin. In this paper, we propose to analyze earthquake-induced site effect derived from horizontal to vertical spectral ratio from ambient noise (noise horizontal to vertical—NHV), or from very weak, weak, moderate and strong ( peak ground acceleration—PGA >10 % g) seismic motions (EHV), and transfer function evaluation from soil velocity profile data at four sites. H/V spectral ratios are computed by using both Fourier and response spectra. Compared to the transfer function, the obtained results show that in the case of soft soils, NHV as well as EHV give a good estimation of the soil's fundamental frequency, whereas the NHV underestimate the H/V amplitude and the EHV amplitude increase with the seismic motion intensity. In the case of firm soils, whereas the NHV gives flat curves synonymous for a rock site or a bump, the EHV is more appropriate as seen by identifying clear peaks with non negligible amplitude. In the case of soft sites as well as in the case of firm sites, strong peaks at frequencies higher than the fundamental one are found from EHV curves. Those peaks would not be found when looking at NHV ratios alone, are evidenced by the computed transfer function as well as by an analytical formulation, and are in agreement with the observed distribution of damage during the M _w 6.8 2003 Boumerdes earthquake. Finally, the same analysis is performed by using response spectra rather than Fourier spectra and leads to the same conclusions. Moreover, the calculation of the response spectra is more easily compared with the smoothing operation of the Fourier spectra.     

Near-source ground motion attenuation relationship for PGA and PSA of the vertical and horizontal components

In this paper, empirical ground-motion models for the vertical and average horizontal components of peak ground-motion and acceleration response spectra from shallow crustal earthquakes are derived using near-source database. These attenuation relationships were derived using a worldwide dataset consisted of corrected and processed accelerograms of 678 strong-motion records recorded with 60 km of the rupture plane of earthquakes between M_w 5.2 and 7.9. Ground motion models are functions of earthquake mechanism, distance from source to site, local average shear wave velocity, nonlinear soil response, sediment depth, depth-to-top of the rupture, hanging wall effects and faulting mechanism.     

烟台海岸软土场地特征及对地震动参数的影响

本文利用剪切波速和钻孔资料，根据中国GB50011-2001规范、美国FEMA-NEHRP规范、欧洲EUROCODE 8规范，对烟台滨海软土场地类别进行评定。并通过场地地震动效应分析，研究滨海场地土特征对地震动的影响。结果表明，简单的场地类别划分难以反映软土场地的实际情况；滨海软土场地对地震动加速度的幅值和频谱的影响非常明显，其对峰值加速度的放大倍数较高，但受土的非线性变形影响，随地震动输入的增大其放大倍数呈非线性降低；软土场地使得加速度反应谱谱形变宽，特征周期变大。     

Statistical analysis of peaks and directivity of earthquake ground motion

Conversion factors are useful for attenuation and damage estimation relationships. These factors among different definitions of peaks (i.e. larger, average and resultant) for peak ground motion indices and acceleration response spectrum were investigated. A large number of horizontal acceleration records recorded at 76 free-field sites of the Japan Meteorological Agency were used in this study. Two orthogonal horizontal components were combined in the time domain to get the maximum resultant peak ground motion indices and acceleration response spectrum in the horizontal plane. From the analysis, the means of the larger/resultant ratio were found to be 0·934 for acceleration, 0·926 for velocity, and 0·913 for displacement. A similar decreasing trend was observed for the means of the average/resultant ratio of the ground motion indices and acceleration response spectrum. The directivity of peak ground motion indices was also examined. It was found that the peak ground motion is more likely to occur in the transverse direction than in other directions. This trend is more prominent in the long-period contents of ground motion.     

施甸5.9级地震强震观测资料的时-频反应谱分析

通过对施甸5．9级地震在太平等观测点的强震记录进行时一频反应谱计算,给出了三个方向的时一频反应谱图,根据时一频反应谱的特点,结合场地条件,对这些场点的震害情况进行了分析,指出了在结构抗震方面应注意的环节,并对时一频反应谱的应用进行了探讨。     

Short‐term monitoring of foredune formation on the east coast of South Africa

Although dunes fronted by sandy beaches constitute approximately 80 per cent of South Africa's coastline, few studies have addressed the formation and life cycle of coastal foredunes, the small, ephemeral shore‐parallel dune ridges typically less than 5 m high and 20 m wide, which form seaward of the storm line. This study used regular, detailed topographic surveys of embryo and foredunes at Tugela mouth, an aggrading stretch of shoreline on the subtropical east coast of South Africa, over a 32‐month period, to gain insight into the formation and motion of these highly mobile landforms over the short term. Average wind drift potential at Tugela mouth during the study period, at 2·35 m s^?1, was an order of magnitude lower than that typical of most parts of the eastern South African coast. The dominant sand‐moving wind for the region was from the southwest to west‐southwest at 10·7 to 13·8 m s^?1, with a secondary vector from north to north‐northeast at 10·8 to 13·8 m s^?1. Signi?cant shoreline retreat, a result of the low sediment yield of the Tugela River during the study period, was one of the main results. This provided the context for redistribution of sand from the inland to the seaward side of the study area, a consequence of the dominant wind direction, and for frequent creation and destruction of short‐lived embryo dunes. Those foredunes which survived the whole study period tended to increase in height, but there was no consistent directional trend in foredune crest movement throughout the 32 months. The study results generally supported Psuty's model of foredune development, but could not con?rm his contention of landward retreat of dune forms under conditions of shoreline erosion. This may be due to the relatively short duration of the study, or possibly to low wind drift potential at the site. Copyright © 2003 John Wiley & Sons, Ltd.