近断层地震动模拟现状

地震动是由3个物理过程（震源破裂过程、波传播过程、场地反应）组成的一种复杂系统的产物,地震动模拟均是围绕这3个物理过程的建模开展的。地震动模拟目前仍然是一门相对较新的科学,强震观测中不断发现的新情况、新问题及其深入研究进一步推动近断层地震学理论和实践的发展。减少建模中的不确定性,用基于观测物理学的统计特征逐渐取代基于现象的假设描述,以改善地震动模拟的精度。基于大量地震动模拟的研究文献和资料,归纳、评述了近断层地震动模拟方法的现状、3个物理过程的建模方法及其发展趋势。     

Analysis of nonstationary properties of ground motion recorded during the Montenegro earthquake of 15 April 1979

The accelerogram of a strong ground motion represents the detailed characteristics of an earthquake and contains a large number of data both for earthquake engineers and seismologists. It shows not only the intensity, the frequency content and the duration of the strong ground motion, but also data referring to source origin, wave propagation mechanism, geological and site conditions at the recording site and its closer surrounding. Therefore, recently, accelerograms analysis has become a widely accepted method for fundamental investigations referring to understanding the background of the seismic waves generation and propagation mechanism and investigation of the strong ground motion nature and the parameters defining the ground motion for the needs of seismic resistant design.Presented in this paper are the results from the nonstationary analysis of the horizontal components of ground motion records obtained during the Montenegro earthquake of 15 April 1979. For the purpose of the analysis, the physical spectra concept was applied, since it enables simultaneous investigation of the time varying nature of the amplitude and the frequency content of the recorded ground motion. The application of the physical spectra for the design of a realistic synthesized accelerogram is also presented.     

典型场地在不同地震动输入水平下的反应

强震动作用下场地土的非线性对地震动反应特征影响明显。场地土的非线性效应首先是由岩土工程学家在室内试验中发现的,继而地震学家在强震动记录分析中亦发现了场地土的非线性效应。一般来讲,随着地震动输入强度增加,场地土非线性效应越明显。但对于不同场地类型,相同地震动水平作用下,场地土非线性效应是不同的。针对这一问题,选择了几个典型场地,采用等效线性化方法分析了不同地震动强度作用下各类场地土的非线性效应,并提出了一种以场地等效剪切波速随地震动输入水平变化为依据的场地土非线性效应确定方法,且定义当等效剪切波速变化率不足20%时,场地地震反应分析可不计入场地土的非线性效应。     

A technical note on seismic microzonation in the central United States

Microzonation is an effort to evaluate and map potential hazards found in an area, urban area in particular, that could be induced by strong ground shaking during an earthquake. These hazards include: ground motion amplification, liquefaction, and slope failure. The microzonation maps, depicting ground-motion amplification, liquefaction, and landslide potentials, can be produced if the ground motion on bedrock (input) and the site conditions are known. These maps, in combination with ground-motion hazard maps (on bedrock), can be used to develop a variety of hazard mitigation strategies such as seismic risk assessment, emergency response and preparedness, and land-use planning. However, these maps have certain limitations that result from the nature of regional mapping, data limitations, generalization, and computer modeling. These microzonations show that when strong ground shaking occurs, damage is more likely to occur, or be more severe, in the higher hazard areas. The zones shown on the hazard maps should not serve as a substitute for site-specific evaluations.     

近源地震动峰值加速度衰减关系影响因素分析

本文收集了丰富的强震资料, 以峰值加速度为例, 采用简单且体现近场峰值加速度PGA震级饱和和距离饱和特性的衰减模型, 研究了表征震级与距离饱和效应的R0(M)的性质。R0与震级相关, 同时与震源性质、地震波频谱有关。在单个地震的R0(M)的求取中, 由于R0与系数d几乎呈线性关系, 所以要求单个地震R0的值, 必须先根据理论约束确定d的大小。在检验衰减方程的预测效果时, 不仅要判断衰减曲线是否反映了实测资料的平均变化趋势, 而且要判断实测资料是否绝大多数落在84%及16%概率水平的预测曲线之内(之间).     

The SISMA prototype system: integrating Geophysical Modeling and Earth Observation for time-dependent seismic hazard assessment

An innovative approach to seismic hazard assessment is illustrated that, based on the available knowledge of the physical properties of the Earth structure and of seismic sources, on geodetic observations, as well as on the geophysical forward modeling, allows for a time-dependent definition of the seismic input. According to the proposed approach, a fully formalized system integrating Earth Observation data and new advanced methods in seismological and geophysical data analysis is currently under development in the framework of the Pilot Project SISMA, funded by the Italian Space Agency. The synergic use of geodetic Earth Observation data (EO) and Geophysical Forward Modeling deformation maps at the national scale complements the space- and time-dependent information provided by real-time monitoring of seismic flow (performed by means of the earthquake prediction algorithms CN and M8S) and permits the identification and routine updating of alerted areas. At the local spatial scale (tens of km) of the seismogenic nodes identified by pattern-recognition analysis, both GNSS (Global Navigation Satellite System) and SAR (Synthetic Aperture Radar) techniques, coupled with expressly developed models for interseismic phase, allow us to retrieve the deformation style and stress evolution within the seismogenic areas. The displacement fields obtained from EO data provide the input for the geophysical modeling, which eventually permits to indicate whether a specific fault is in a “critical state.” The scenarios of expected ground motion (shakemaps) associated with the alerted areas are then defined by means of full waveforms modeling, based on the possibility to compute synthetic seismograms by the modal summation technique (neo-deterministic hazard assessment). In this way, a set of deterministic scenarios of ground motion, which refer to the time interval when a strong event is likely to occur within the alerted area, can be defined both at national and at local scale. The considered integrated approach opens new routes in understanding the dynamics of fault zones as well as in modeling the expected ground motion. The SISMA system, in fact, provides tools for establishing warning criteria based on deterministic and rigorous forward geophysical models and hence allows for a well-controlled real-time prospective testing and validation of the proposed methodology over the Italian territory. The proposed approach complements the traditional probabilistic approach for seismic hazard estimates, since it supplies routinely updated information useful in assigning priorities for timely mitigation actions and hence it is particularly relevant to Civil Defense purposes.     

Early Estimation of Seismic Hazard for Strong Earthquakes in Taiwan

A shakemap system providing rapid estimates of strong ground shaking could be useful for emergency response providers in a damaging earthquake. A hybrid procedure, which combines site-dependent ground motion prediction models and the limited observations of the Real-Time Digital stream output system (RTD system operated by Central Weather Bureau, CWB), was set up to provide a high-resolution shakemap in a near-real-time manner after damaging earthquakes in Taiwan. One of the main factors that affect the result of ground motion prediction analysis is the existence of site effects. The purpose of this paper is to investigate the local site effects and their influence in the ground shaking and then establish an early estimation procedure of potential hazard for damaging earthquakes. Based on the attenuation law, the site effects of each TSMIP station are discussed in terms of a bias function that is site and intensity-level dependent function. The standard deviation of the site-dependent ground motion prediction model can be significantly reduced. The nonlinear behavior of ground soil is automatically taken into account in the intensity-level dependent bias function. Both the PGA and the spectral acceleration are studied in this study. Based on the RTD data, event correctors are calculated and applied to precisely estimate the shakemap of damaging earthquakes for emergency response.     

基于袋装砂井排水固结法处理软基的沉降过程流固耦合模拟

在软土分布地区,高速公路的“桥头跳车”、路面开裂等一直是困扰工程建设者的最主要问题,而对软土路基的沉降过程进行模拟预测则是其中的难点。利用改进的剑桥模型,采用FLAC3D有限差分软件,以珠江三角洲某高速公路试验段的地层数据进行袋装砂井排水固结法处理软土地基沉降过程的流固耦合模拟,通过模拟结果来对比实际检测的应力、应变情况。结果表明,表面沉降和分层沉降的模拟值与实测值能很好地吻合,模型能较为真实地反映路基沉降变形的实际情况。     

Estimation of seismic spectral acceleration in Peninsular India

Peninsular India (PI), which lies south of 24°N latitude, has experienced several devastating earthquakes in the past. However, very few strong motion records are available for developing attenuation relations for ground acceleration, required by engineers to arrive at rational design response spectra for construction sites and cities in PI. Based on a well-known seismological model, the present paper statistically simulates ground motion in PI to arrive at an empirical relation for estimating 5% damped response spectra, as a function of magnitude and source to site distance, covering bedrock and soil conditions. The standard error in the proposed relationship is reported as a function of the frequency, for further use of the results in probabilistic seismic hazard analysis.     

Modeling of strong motion generation areas of the 2011 Tohoku,Japan earthquake using modified semi-empirical technique

Modification in the semi-empirical technique for the simulation of strong ground motion has been introduced to incorporate the strong motion generation areas (SMGA) in the modeled rupture plane. Strong motion generation areas identified within the rupture plane of the Tohoku earthquake of March 11, 2011 (M _w = 9.0), have been modeled using this modified technique. Two different source models having four and five SMGAs, respectively, are considered for modeling purpose. Strong motion records using modified semi-empirical technique have been simulated at two near-field stations located at epicentral distance of 137 and 140 km, respectively, using two different source models. Comparison of the observed and simulated acceleration waveforms is made in terms of root mean square error (RMSE) at both stations. Minimum root mean square error of the waveform comparison has been obtained at both the stations for source model having five SMGAs. Simulations from same rupture model have been made at other four stations lying at epicentral distance between 154 and 249 km. Comparison of observed and simulated records has been made in terms of RMSE in acceleration records, velocity records and response spectra at each six station. Simulations have been made at six other stations to obtain distribution of peak ground acceleration and peak ground velocity with hypocentral distance. Peak ground acceleration and velocity from simulated and observed records are compared at twelve stations surrounding the source of Tohoku earthquake. Comparison of waveforms and parameters extracted from observed and simulated strong motion records confirms the efficacy of the developed modified technique to model earthquake characterized by SMGAs.     

Near-source strong motion database catalog for Iran

This paper discusses a newly developed high-quality integrated dataset of shallow earthquake ground motions that occurred in Iran, from 1976 to 2013. A total of 860 three-component strong motion records are processed from 183 earthquake events, moment magnitudes 5.0?≤?M _w ?≤?7.4, and rupture distances of R _RUP  ≤ 120 km. Strong motion data from Iran having special tectonic features and shallow earthquakes with depths less than 35 km are included. This paper presents a thorough procedure used to collect and to generate a database following the Next-Generation Attenuation-West research projects. This database can be used in the development and ranking of ground motion models and for seismological and engineering hazard and risk analyses. Unprocessed strong motion records are obtained from the Iranian Strong Motion Network (ISMN). The time series collected were thoroughly examined through several rounds of quality reviews. The newly generated database includes the peak ground acceleration, peak ground velocity, and pseudo-spectral acceleration for the 5% damped with periods ranging from 0.01 to 10 s. The database also includes ground motion information and source characterization and parameters. This study is the near-source compiled ground motion database that can be used for Iran, and it is consistent with standard worldwide databases.     

A review of strong motion studies in Gujarat State of western India

The present work reviews the strong motion studies done in Gujarat State of western India. Prior to the 2001 Bhuj earthquake, no strong motion instrument was in operation in Gujarat. After the earthquake, number of research institutes/universities from India and abroad deployed strong motion instruments to study aftershock activity, source dynamics, path and site effects. The strong motion recordings have enhanced the general understanding of the physics of earthquakes in the region. An attempt has been made to develop attenuation relationship for the Gujarat region from the actual ground motions recorded by the strong motion networks. The Government of Gujarat with the help from Asian Development Bank, World Bank (WB), Ministry of Science and Technology and Ministry of Earth Sciences, Government of India, has established a permanent dense network of strong motion accelerograph (SMA) all over Gujarat. In addition, the Institute of Seismological Research has been established in Gandhinagar, Gujarat, with the help of WB for carrying out seismological research. Recently, many important studies have been carried out using actual acceleration data obtained from a dense network of 54 SMA, as well as synthetic data generated using region-specific ground motion parameters. The recorded data are used to obtain region-specific ground motion parameters and ground motion prediction equation. A deterministic hazard analysis for the entire state of Gujarat has been carried out using site-specific ground motion parameters. The estimated peak ground acceleration and modified Mercalli intensity values have been used to estimate the vulnerability of the different types of buildings in 31 cities of Gujarat. As Gujarat has three distinct regions having varied geological conditions, the recorded strong motion data gave an opportunity to study the effect of geological and local-site conditions on the response spectra. This study for an intra-plate region like Gujarat is a pioneer work. Still, lots of research work need to be carried out as more and more data are available, such as development of more robust ground motion prediction equations and a 3D-velocity structure of Gujarat. Generation of shake maps in real time and a credible early earthquake warning system is need of the hour for disaster mitigation and management.     

Simulation of non-stationary conditional ground motion fields in the time domain

This paper addresses the topic of simulating spatially variable ground motion fields conditioned on a known accelerogram. The conditional ground motion fields can be used in design or verification studies where seismic analysis has to be performed for a couple of natural accelerograms that have been preselected by seismologists or other experts. The methodology is based on conditional densities. In contrast to most authors, the conditional densities method is not applied to the Fourier coefficients, but it is used for the construction of a conditional Gaussian process model in the time domain. This has the advantage that fully non-stationary conditional time histories can be simulated directly in the time domain. The cross-correlation functions needed for this approach are evaluated from commonly used ground motion models expressed as evolutionary power spectral densities. An application to the El Centro earthquake record is presented. The properties of the simulated ground motion fields are analysed and compared to the data and the theoretical model.     

Strong ground motion attenuation in Aswan area, Egypt

Recent and paleo seismicity indicate that moderate seismic activity is relatively large for Aswan area. This is a warning on the possibility of occurrence of earthquakes in the future too. No strong motion records are available in Aswan area for engineers to rely upon. Consequently, the seismological modeling is an alternative approach till sufficient instrumental records around Aswan become available. In the present study, we have developed new ground motion attenuation relationship for events spanning 4.0?≤? M _w?≤?7.0 and distance to the surface projection of the fault up to 100 km for Aswan based on a statistically simulated seismological model. We generated suites of ground motion time histories using stochastic technique. The ground motion attenuation relation describes the dependence of the strength of the ground motions on the earthquake magnitude and distance from the earthquake. The proposed equation for peak ground acceleration (PGA) for the bed rock is in the form of: $ {\mathbf{log}}{\text{ }}\left( {{\mathbf{PGA}}/{\mathbf{gal}}} \right){\text{ }} = {\mathbf{1}}.{\mathbf{24}} + {\mathbf{0}}.{\mathbf{358}}{M_{\mathbf{w}}} - {\text{ }}{\mathbf{log}}\left( {\mathbf{R}} \right){\text{ }}-{\text{ }}{\mathbf{0}}.{\mathbf{008}}{\text{ }}{\mathbf{R}}{\text{ }} + {\text{ }}{\mathbf{0}}.{\mathbf{22}}{\text{ }}{\mathbf{P}} $ . Where PGA is the peak ground acceleration in gal (cm/s²); M_w, its moment magnitude; R is the closest distance between the rupture projection and the site of interest; and the factor P is a dummy variable. It is observed that attenuation of strong motion in Aswan is correlated with those used before in Egypt.     

强地震动作用下层状岩体破坏的物理模拟研究

岩体地震动力破坏现象是常见但研究较少的复杂课题。利用物理模拟试验,研究了层状岩体边坡在强地震动作用下的变形破坏问题。按结构面走向平行和垂直地震动方向,对不同岩体结构特征的模型施加不同的幅值和频率的水平向地震动荷载,进行岩体的振动破坏试验。结果表明,结构面空间展布方向与地震动荷载的相互关系不同、结构面软弱程度不同,岩体变形破坏的形式和分布特征也不相同。岩体的动力破坏主要是岩体结构的破坏,即岩体结构的宏观破坏和微观损伤,岩体结构特征是控制其动力变形破坏的主要因素。极大的不均匀性是岩体动力破坏的显著特点,这与大量极震区岩体破坏的现象类似。但岩体动力破坏模式的建立和定量的力学分析还需更多研究。     

Estimation of the spatial distribution of ground motion parameters for two recent earthquakes in Japan

A recent development in strong motion instrumentation in Japan provides an opportunity to collect valuable data sets, especially after moderate and large magnitude events. Gathering and modeling these data is a necessity for better understanding of regional ground motion characteristics. Estimations of the spatial distribution of earthquake ground motion plays an important role in early-stage damage assessments for both rescue operations by disaster management agencies as well as damage studies of urban structures. Subsurface geology layers and local soil conditions lead to soil amplification that contributes to the estimated ground motion parameters of the surface. We present a case study of the applicability of the nationally proposed GIS-based soil amplification ratios [J. Soil Dyn. Earthqu. Eng. 19 (2000) 41–53] to the October 6, 2000 Tottori-ken Seibu (western Tottori Prefecture) and the March 24, 2001 Geiyo earthquakes in Japan. First, ground motion values were converted to those at a hypothetical ground base-rock level (outcrop) using an amplification ratio for each 1×1 km area, based on geomorphological and subsurface geology information. Then a Kriging method, assuming an attenuation relationship at the base-rock as a trend component, is applied. Finally, the spatial distribution of ground motion at ground surface is obtained by applying GIS-based amplification factors for the entire region. The correlation between the observed and estimated ground motion values is reasonable for both earthquakes. Thus, the proposed method is applicable in near real-time early-damage assessments and seismic hazard studies in Japan.     

近断层速度脉冲地震动引起的边坡永久位移预测模型

基于小波分析方法,从NGA数据库的3 551条地震记录中选取189条速度脉冲地震动,地震动均转换成发生最强脉冲的方向。基于Newmark方法,分析了近断层速度脉冲地震动作用引起的边坡永久位移值。结果表明：近断层速度脉冲地震动对边坡产生特殊的破坏作用,表现在滑动位移值大、滑动体破坏力强等方面;边坡永久位移值与速度脉冲地震动的峰值速度具有高度相关性,位移值较大时尤为明显。建立了基于单变量形式的峰值速度及双变量形式的峰值速度、峰值加速度两种边坡永久位移预测模型,模型简单实用,与回归数据具有很好的相关性,前者更适用于预测对实际工程影响较大的永久位移值,且离散性较小。提出的预测模型为考虑近断层地震动速度脉冲特性影响的边坡永久位移值的概率地震灾害分析提供了基础。     

A Multi-Parameter Correlation for Predicting the Seismic Displacement of an Earth Dam or Embankment

Predictive displacement-based methods provide a useful index of the seismic performance of earth dams and embankments and can be used in preliminary assessments of these structures. In practice, simplified Newmark-type sliding block methods are commonly used for this purpose. Using a database of 122 previously published case histories of permanent deformations of earth dams and embankments, the performance of six simplified sliding block models was examined. The results show that all six simplified methods underpredict seismic displacement for many of the embankment and earth dam cases that were examined, sometimes by a significant amount. An empirical correlation was developed by performing linear multiple regression analysis utilizing multiple slope and ground motion input parameters. This approach is believed to more properly reflect strong ground motion characteristics than the use of a single ground motion parameter such as the peak ground acceleration, the approach that has been previously employed in other correlations of this type. After exploring numerous functional forms, the final resulting seismic displacement correlation that was proposed was determined to be a function of the critical acceleration, the critical acceleration ratio, the slope height, the peak ground acceleration, the peak ground velocity, the spectral acceleration, and the predominant period of earthquake shaking. The proposed empirical equation yields better correlation with the case history database than does other existing empirical correlations or simplified sliding block models.     

Generation of stochastic earthquake ground motion in western Saudi Arabia as a first step in development of regional ground motion prediction model

Earthquake ground motion model is an essential part of seismic hazard assessment. The model consists in several empirical ground motion prediction equations (GMPEs) that are considered to be applicable to the given region. When the recorded ground motion data are scarce, numerical modeling of ground motion based on available seismological information is widely used. We describe results of stochastic simulation of ground motion acceleration records for western Saudi Arabia. The simulation was performed using the finite fault model and considering peak ground acceleration and amplitudes of spectral acceleration at natural frequencies 0.2 and 1.0 s. Based on the parameters of the input seismological model that were accepted in similar previous studies, we analyze influence of variations in the source factor (stress drop) and in the local attenuation and amplification factors (kappa value, crustal amplification). These characteristics of the model are considered as the major contributors to the ground motion variability. The results of our work show that distribution of simulated ground motion parameters versus magnitude and distance reveals an agreement with the GMPEs recently used in seismic hazard assessment for the region. Collection of credible information about seismic source, propagation path, and site attenuation parameters using the regional ground motion database would allow constraining the seismological model and developing regional GMPEs. The stochastic simulation based on regional seismological model may be applied for generation of ground motion time histories used for development of analytical fragility curves for typical constructions in the region.     

黄土高填方场地工后沉降预测模型性能评估方法

工后沉降预测结果是黄土高填方场地变形稳定性评价和建筑物规划布局的重要参考依据。为遴选适合黄土高填方场地的工后沉降预测模型,基于某典型黄土高填方工程的实测沉降数据,分析了工后沉降曲线的变化规律和发展趋势,建立了17种回归参数模型,提出了模型预测效果的评价指标和方法。结果表明：（1）该工程填方区工后沉降历时曲线呈“缓变型”变化,土方填筑完工初期无陡增段,随时间增加沉降速率逐步降低,尚未出现沉降趋于稳定的水平段;（2）将外推预测误差、内拟合误差和后验误差比最小化作为综合控制目标,可遴选出理想的回归参数模型;（3）MMF模型（Ⅱ型）和双曲线模型具有较高的预测精度、较好的稳定性和较强的适应性,在17种模型中的预测效果最佳;（4）沉降数据的变化越平稳,模型预测效果越好;（5）增大建模数据的时间跨度,会提升预测精度,但增大至一定值后,预测精度提升效果不再显著。