2011年新西兰Mw6.1地震震源过程及强地面运动特征初步分析

2011年2月21日新西兰南岛发生M_w6.1地震,震中位置43.58°S,172.70°E,震源深度约5.0 km,发震断层为新西兰第三大城市克赖斯特彻奇(Christchurch)南约9 km一条近东西走向的未知隐伏断层,为逆冲断层机制.地震已经造成163人遇难,Christchurch城内多处建筑物严重损毁,距震中约1 km的Heathcote Valley Primary School (HVSC)台站强地面运动峰值加速度(PGA)高达2.0g.针对新西兰M_w6.1地震近场强地面运动偏高这一现象,利用HVSC台站的强震观测记录,计算地震震源谱参数,应用Brune圆盘模型估算其近场强地面运动的理论值,并建立动态复合震源模型进行模拟计算.研究结果表明,新西兰M_w6.1地震近断层强地面运动偏高的主要原因为复杂震源破裂过程中有效应力降(动态应力降,Δσ_d)过高造成的.未来工作中,需要加强对可能发生的、距离城市较近的中小型地震的重视,防止地震对城市的加强型破坏.     

芦山地震强地面运动之阿里亚斯强度及Newmark位移特征研究

本文以芦山地震强地面运动记录为基础资料,研究了阿里亚斯强度和Newmark位移两个地震动参数的空间分布特征、衰减特征以及与其它地震动参数的相关性.研究结果表明:阿里亚斯强度的空间分布与地震断层空间展布和地震破裂方向具有相关性;阿里亚斯强度与峰值加速度(PGA)有较好的相关性,场地条件对二者的相关性具有显著影响,PGA相...     

地震动衰减模型在玛多Mw7.3级地震影响场快速评估中的应用

2021年5月22日青海省玛多县发生M_w7.3级地震。震后,根据初步估计的断层走向和破裂长度,基于YU15地震动衰减模型和三种NGA-West2（Next Generation Attenuation-West2）地震动衰减模型快速产出地震区震动图及理论烈度图。在获得强震记录和地表破裂长度信息后,对预测结果进行修正。通过比较理论烈度与调查烈度,并结合震动图分布形态以及衰减模型在2016年新疆呼图壁M_w6.0地震中的应用情况对四种地震动衰减模型的适用性进行了分析。结果表明:在台网稀疏地区,基于地震动衰减模型可在震后快速获得地震动分布,并产出具有应用价值的地震影响场;NGA-West2模型在断层破裂较长的大震中表现优于YU15模型,而在中强地震中后者适用性更强;近实时强震动记录可用来检验模型的适用性并对预测结果进行修正;断层破裂尺度、震源机制和破裂过程等信息的准确估计可有效提高地震影响场预测精度。     

New predictive equations for Arias intensity from crustal earthquakes in New Zealand

Arias Intensity (Arias, MIT Press, Cambridge MA, pp 438–483, 1970) is an important measure of the strength of a ground motion, as it is able to simultaneously reflect multiple characteristics of the motion in question. Recently, the effectiveness of Arias Intensity as a predictor of the likelihood of damage to short-period structures has been demonstrated, reinforcing the utility of Arias Intensity for use in both structural and geotechnical applications. In light of this utility, Arias Intensity has begun to be considered as a ground-motion measure suitable for use in probabilistic seismic hazard analysis (PSHA) and earthquake loss estimation. It is therefore timely to develop predictive equations for this ground-motion measure. In this study, a suite of four predictive equations, each using a different functional form, is derived for the prediction of Arias Intensity from crustal earthquakes in New Zealand. The provision of a suite of models is included to allow for epistemic uncertainty to be considered within a PSHA framework. Coefficients are presented for four different horizontal-component definitions for each of the four models. The ground-motion dataset for which the equations are derived include records from New Zealand crustal earthquakes as well as near-field records from worldwide crustal earthquakes. The predictive equations may be used to estimate Arias Intensity for moment magnitudes between 5.1 and 7.5 and for distances (both r_jb and r_rup) up to 300 km.     

Consistency of ground-motion predictions from the past four decades

Due to the limited observational datasets available for the derivation of ground-motion prediction equations (GMPEs) there is always epistemic uncertainty in the estimated median ground motion. Because of the increasing quality and quantity of strong-motion datasets it would be expected that the epistemic uncertainty in ground-motion prediction (related to lack of knowledge and data) is decreasing. In this study the predicted median ground motions from over 200 GMPEs for various scenarios are plotted against date of publication to examine whether the scatter in the predictions (a measure of epistemic uncertainty) is decreasing with time. It is found that there are still considerable differences in predicted ground motions from the various GMPEs and that the variation between estimates is not reducing although the ground motion estimated by averaging median predictions is roughly constant. For western North America predictions for moderate earthquakes have show a high level of consistency since the 1980s as do, but to a lesser extent, predictions for moderate earthquakes in Europe, the Mediterranean and the Middle East. A good match is observed between the predictions from GMPEs and the median ground motions based on observations from similar scenarios. Variations in median ground motion predictions for stable continental regions and subduction zones from different GMPEs are large, even for moderate earthquakes. The large scatter in predictions of the median ground motion shows that epistemic uncertainty in ground-motion prediction is still large and that it is vital that this is accounted for in seismic hazard assessments.     

南北地震带南段水平向地震动衰减特征

地震动衰减关系作为抗震救灾的重要依据,一直都是现代地震学研究的重点之一.本研究使用南北地震带南段区域2009—2016年共217个地震事件获得随震中距变化的水平向地震动峰值速度(PGV)和地震动峰值加速度(PGA)经验衰减关系,并计算场地响应.研究结果显示PGV衰减关系的拟合效果较PGA更好,两者的距离衰减系数会随事件矩震级增大呈线性减小;相较于大矩震级事件,中矩震级事件在近场可能产生较衰减关系理论值更大的PGV和PGA,同时衰减关系的拟合标准差会随事件矩震级的增大而减小.进行场地响应校正后的PGV和PGA更加符合经验衰减关系,PGA的场地响应影响较PGV更强但两者的趋势一致,并与该区域前人计算得到的地壳Qs值分布对应,表明地壳介质放大或压制地震波振幅和其传递地震波能量的能力是相关联的.本文结果一定程度上揭示了南北地震带南段的地震动强度衰减特征,为未来中国西南部的抗震减灾工作提供了重要的参考.     

地形效应影响下地震动参数与斜坡稳定性的相关性研究

地形对地震动的影响比较复杂, 考虑地形放大效应的地震滑坡稳定性分析需要选择合适的地震动参数. 本文使用自贡地形影响台阵记录到的2008年汶川M_S8.0地震主震加速度记录, 分析了地震动峰值加速度、 阿里亚斯烈度以及90%能量持时随地形高度的变化, 探讨了地形效应作用下峰值加速度和阿里亚斯烈度与地震动作用下斜坡稳定性的相关性. 结果表明: ① 地形场地对峰值加速度和阿里亚斯烈度均有显著的放大效应. 地形放大效应较为复杂, 其整体上随台站高度的增加而增大, 水平向的放大效应大于竖直向. 水平向峰值加速度的放大系数为1.1—1.8, 阿里亚斯烈度的放大系数为1.2—3.3; 竖直向相应放大系数分别为1.1—1.3和1.2—1.7. ② 地形对地震动持时也有一定的放大效应, 但不同高度、 不同分量的放大效应没有显著差异, 其放大系数均约为1.3. ③ 阿里亚斯烈度和峰值加速度均能很好地表征地形对地震动的影响, 与地震动对斜坡稳定性的影响具有很强的相关性. 与峰值加速度相比, 阿里亚斯烈度综合了地震动的多方面特征, 可以更好地表征地形对地震动的影响, 与地震动作用下斜坡稳定性的相关性更强.      

Simulation of strong earthquake characteristics of a scenario earthquake (MS7.5) based on the enlightenment of 2022 MS6.9 earthquake in Menyuan

The Menyuan area is an important transportation hub in the Hexi Corridor. The Menyuan M_S6.9 earthquake that occurred on January 8, 2022 had a major impact on the local infrastructure and transportation of this region. Due to the high possibility of similar strong earthquakes occurring in this area in the future, preliminary assessment of the seismic intensity characteristics of destructive earthquakes in this region is essential for effective disaster control. This paper uses the empirical Green′s function (EGF) method as a numerical simulation tool to predict the ground motion intensity of Datong Autonomous County under the action of the scenario earthquake (M_S7.5). Seismic records of aftershocks of the 2016 Menyuan M_S6.4 earthquake were used as Green’s functions for this simulation. The uncertainties associated with various source parameters were considered, and 36 possible earthquake scenarios were simulated to obtain 72 sets of horizontal ground motions in Datong County. The obtained peak ground acceleration (PGA) vs. time histories of the horizontal ground motion were screened using the attenuation relationships provided by the fifth-edition of China's Seismic Ground Motion Parameter Zoning Map and the NGA-West2 dataset. Ultimately, 32 possible acceleration-time histories were selected for further analysis. The screened PGA values ranged from 78.8 to 153 cm/s². The uncertainty associated with the initial rupture point was found to greatly affect the results of the earthquake simulation. The average acceleration spectrum of the selected acceleration-time history exceeded the expected spectrum of a intermediate earthquake, which means that buildings in Datong County might sustain some damage should the scenario earthquake occur. This research can provide reliable ground motion input for urban earthquake damage simulation and seismic design in Datong County. Growing the dataset of small earthquakes recorded in this region will facilitate the large-scale simulation of ground motions under different earthquake scenarios.     

2015年尼泊尔Gorkha地震强地面运动记录分析

2015年4月25日在尼泊尔Gorkha地区发生M_W7.8地震,距离发震断层约11 km的KATNP台站完整记录了主震的加速度时程.本文根据KATNP台站记录的加速度数据分析了Gorkha地震的地震动特征.结果表明Gorkha地震在KATNP台站处产生的水平向峰值加速度为0.17 g,竖直向峰值加速度为0.19 g,该数值小于科学家们对如此大规模地震产生的地震动的预期,初步推测这可能是由加德满都山谷产生的非线性响应造成的（Dixit et al.,2015）;地震在KATNP台站处产生了地表永久位移,其中竖向永久位移为131.9 cm,水平向永久位移的绝对值为159.2 cm,方向为南偏西19°（199°）,据此可简单推算出断层走向约为289°（109°）.地震产生了脉冲型地震动,影响因素有盆地效应、地震破裂的向前的方向性效应以及滑冲效应,其中盆地效应的周期约为5 s左右,方向性效应产生的速度脉冲的周期约为8 s左右.加速度反应谱显示在0.5 s和5.0 s左右各有一个峰值,前者是由地震破裂的脉冲式滑移产生的大量高频地震动造成的,后者是由于盆地效应和地震破裂的方向性效应造成的.基于阿里亚斯烈度计算的地震动持时约在36~46 s之间,小于与其规模相当的地震产生的地震动持时,并且不同方向上的地震动持时可能与地震破裂方向有关.阿里亚斯烈度随时间的变化比较简单,也反映了Gorkha地震是一次连续的、能量释放相对简单的地震事件.     

Rockburst induced ground motion—a comparative study

While rockbursts from underground copper mining in Western Poland normally produce surface peak ground accelerations (PGA) and velocities of 0.05–0.1 g and 1–3 cm/s, occasionally these peak motions may exceed 0.15 g and 10 cm/s, respectively. These larger motions are of considerable concern and an investigation has been undertaken to define the nature of these larger induced ground motions. This paper compares these rockburst motions with low intensity earthquakes. Various strong motion parameters such as PGA, peak ground velocity (PGV) and displacements as well as strong motion duration, Arias intensity, Fourier and response spectra are compared with those from earthquakes. It is concluded that although short duration is the most obvious parameter that differentiates rockbursts from earthquakes, in fact their high dominant frequencies, which result in high PGA/PGV ratios differentiate them the most. Two types of rockburst-induced ground motions are indicated in this paper: typical—with 3–6 months return period and characteristic, high frequency content—as well as rare events similar to shallow, low intensity earthquakes.     

APPLICATION OF TWO SIMPLIFIED NEWMARK MODELS TO THE ASSESSMENT OF LANDSLIDES TRIGGERED BY THE 2008 WENCHUAN EARTHQUAKE

At present, with the wide application of the Newmark method, various Newmark empirical formulas with different ground motion parameters have been fitted by many researchers based on global strong-motion records. However, the existing study about the Wenchuan earthquake does not quantitatively evaluate the applicability of different Newmark models based on the actual landslides distribution. The aim of this paper is to present a comparison between observed landslides from the 2008 Wenchuan earthquake and predicted landslides using Newmark displacement method based on different ground motion parameters. The factor-of-safety map and critical acceleration(a_c)map in the study area are obtained by using the terrain data and geological data. The distribution of Arias intensity(I_a)and PGA in the study area is obtained by using the attenuation formulas of Arias intensity(I_a)and PGA, which is regressed by Wenchuan ground motion records. Based on the distribution of Arias intensity(I_a)and PGA parameters, we obtained the predicted locations of landslide using Newmark regression equations which are generated using global strong-motion records. The results shows that the assessment results can better reflect the macroscopic distribution characteristics of co-seismic landslides, most predicted landslide cells are distributed on the two sides of the Beichuan-Yingxiu Fault, especially the Pengguan complex rock mass in the hanging wall. The abilities to predict landslide occurrence of the two Newmark simplified models are different. On the whole, the evaluated result of simplified model based on parameter I_a is better than that based on PGA parameter. The GFC values obtained by the Newmark model of I_a and PGA parameters are 65.7% and 34.9%respectively. The evaluated result based on I_a can better reflect the macro distribution of coseismic landslides. The Ls_Pred value based on the Newmark model of parameter I_a is 26.5%, and the Ls_Pred value based on the Newmark model of PGA parameter is 10.3%. However the total area of predicted landslides accounts for 2.4% of the study area, which indicates that the predicted landslide cells are greater than the observed landslide cells. This reminds us that depending on the current input of shear strength and ground-motion parameters, we can only conduct landslide hazard assessment in macro areas, the ability to predict landslide can be improved using more accurate topographic data and input parameters.     

Ground-motion prediction equations for the intermediate depth Vrancea (Romania) earthquakes

We present the regional ground-motion prediction equations for peak ground acceleration (PGA), peak ground velocity (PGV), pseudo-spectral acceleration (PSA), and seismic intensity (MSK scale) for the Vrancea intermediate depth earthquakes (SE-Carpathians) and territory of Romania. The prediction equations were constructed using the stochastic technique on the basis of the regional Fourier amplitude spectrum (FAS) source scaling and attenuation models and the generalised site amplification functions. Values of considered ground motion parameters are given as the functions of earthquake magnitude, depth and epicentral distance. The developed ground-motion models were tested and calibrated using the available data from the large Vrancea earthquakes. We suggest to use the presented equations for the rapid estimation of seismic effect after strong earthquakes (Shakemap generation) and seismic hazard assessment, both deterministic and probabilistic approaches.     

远场地震波辐射能计算以及对近断层地面运动的约束

从圆盘断层模型出发,根据地震波能量表象定理推导出了滑移弱化过程中远场S-波辐射能量表达式,并同已有的动力学模型作了比较.结果表明,得到的模型能量值或视应力的取值强烈地依赖于断层上的动态、静态应力降和破裂传播速度,而破裂速度则对应了断层带模型中断层破坏过程所耗散的能量.动摩擦应力上调和应力下调力学机制在能量求解中得到了充...     

Characterization of Stress Drops on Asperities Estimated from the Heterogeneous Kinematic Slip Model for Strong Motion Prediction for Inland Crustal Earthquakes in Japan

Dense strong motion observation networks provided us with valuable data for studying strong motion generation from large earthquakes. From kinematic waveform inversion of seismic data, the slip distribution on the fault surface of large earthquakes is known to be spatially heterogeneous. Because heterogeneities in the slip and stress drop distributions control the generation of near-source ground motion, it is important to characterize these heterogeneities for past earthquakes in constructing a source model for reliable prediction of strong ground motion. The stress changes during large earthquakes on the faults recently occurring in Japan are estimated from the detailed slip models obtained by the kinematic waveform inversion. The stress drops on and off asperities are summarized on the basis of the stress change distributions obtained here. In this paper, we define the asperity to be a rectangular area whose slip is 1.5 or more times larger than the average slip over the fault according to the previous study for inland crustal earthquakes. The average static stress drops on the asperities of the earthquakes studied here are in the range 6?C23?MPa, whereas those off the asperities are below 3?MPa. We compiled the stress drop on the asperities together with a data set from previous studies of other inland earthquakes in Japan and elsewhere. The static stress drop on the asperity depends on its depth, and we obtained an empirical relationship between the static stress drop and the asperity??s depth. Moreover, surface-breaking asperities seemed to have smaller stress drops than buried asperities. Simple ground motion simulations using the characterized asperity source models reveal that deep asperities generate larger ground motion than shallow asperities, because of the different stress drops of the asperities. These characteristics can be used for advanced source modeling in strong ground motion prediction for inland crustal earthquakes.     

Sensitivity of ground-motion scenarios to earthquake source parameters in the Tehran metropolitan area,Iran

The northern Tehran fault (NTF) is a principal active fault of the Alborz mountain belt in the northern Iran. The fault is located north of the highly populated Metropolitan Area of Tehran. Historical records and paleoseismological studies have shown that the NTF poses a high seismic risk for the Tehran region and the surrounding cities (e.g. Karaj). A series of ground-motion simulations are carried out using a hybrid kinematic-stochastic model to calculate broadband (0.1–20 Hz) ground-motion time histories for deterministic earthquake scenarios (M7.2) on the NTF. We will describe the source characteristics of the target event to develop a list of scenario earthquakes that are probably similar to a large earthquake on the NTF. The effect of varying different rupture parameters such as rupture velocity and rise time on the resulting broadband strong motions has been investigated to evaluate the range of uncertainty in seismic scenarios. The most significant parameters in terms of ground-shaking level are the rise time and the value of the rupture velocity. For the worst-case scenario, the maximum expected horizontal acceleration, and velocity at rock sites in Tehran range between 128 and 1315 cm/s/s and 11–191 cm/s, respectively. For the lowest scenario, the corresponding values range between 102 and 776 cm/s/s and 12 to 81 cm/s. Nonlinear soil effects may change these results but are not accounted for in this study. The largest variability of ground motion is observed in neighborhood of asperity and also in the direction of rupture propagation. The calculated standard deviation of all ground-motion scenarios is less than 30% of the mean. The capability of the simulation method to synthesize expected ground motions and the appropriateness of the key parameters used in the simulations are confirmed by comparing the synthetic peak ground motions (PGA, PGV and response spectra) with empirical ground-motion prediction equations.     

Empirical attenuation relationship for Arias Intensity

Arias Intensity is a ground motion parameter that captures the potential destructiveness of an earthquake as the integral of the square of the acceleration–time history. It correlates well with several commonly used demand measures of structural performance, liquefaction, and seismic slope stability. A new empirical relationship is developed to estimate Arias Intensity as a function of magnitude, distance, fault mechanism, and site category based on 1208 recorded ground motion data from 75 earthquakes in active plate‐margins. Its functional form is derived from the point‐source model, and the coefficients are determined through non‐linear regression analyses using a random‐effects model. The results show that for large magnitude earthquakes (M > 7) Arias Intensity was significantly overestimated by previous relationships while it was underestimated for smaller magnitude events (M ? 6). The average horizontal Arias Intensity is not significantly affected by forward rupture directivity in the near‐fault region. The aleatory variability associated with Arias Intensity is larger than that of most other ground motion parameters such as spectral acceleration. However, it may be useful in assessing the potential seismic performance of stiff engineering systems whose response is dominated by the short‐period characteristics of ground motions. Copyright © 2003 John Wiley & Sons, Ltd.     

地震动不确定性及其影响因素的初步分析

选取欧洲强震记录数据库中,同一台站记录到的两次地震以及两个不同台站记录到的同一地震的强震记录,分析了它们的PGA比值、卓越频率和反应谱的不确定性以及不确定性对结构反应的影响,得出：影响地震动不确定性的因素主要有震源机制、地震波传播途径、场地类别等。在此基础上,进一步分析了地震动不确定性对结构反应的影响。     

相对局部区域震源参数随机不确定性经验关系研究

在进行未来破坏性地震的强地面运动数值模拟时,震源参数选取的准确性对地震动预测的结果影响很大。震源参数的确定存在很多不确定性因素,既包含随机的不确定性因素,又包含认知的不确定性因素。本文在大量地震事件及文献调研的基础上,运用统计学方法对具备随机不确定性特征的震源参数进行统计研究,以震源参数经验公式的形态建立解释其随机性和不确定性的数学模型。为了研究局部地区震源参数的定标关系特征,获得更加适用于局部地震密集区域,尤其是包含中国大陆地区在内的局部区域的震源参数的经验关系,本文从GCMT地震目录中选取了1 700多个M_W≥5.5的地震事件,运用统计学方法研究地震密集地区的震源参数经验关系,包括震级、地震矩、破裂面积等,增加了相对较大的局部范围内凹凸体的地震样本数量,从统计学角度计算更加适合局部区域的震源参数的经验关系。统计结果表明:局部区域震例获得的震源参数的经验关系与不限区域震例获得的经验关系存在差异,尤其是涉及到断层破裂面积、凹凸体相关参数时差异较大,局部区域内震例获得的震源参数的经验关系将更具有代表性。应用本文获得的相对局部区域的经验公式计算未来破坏性地震的强地面运动所需的震源参数时,获得的地震动预测结果将更能体现目标区域真实的地震动特征,进而提高地震动预测结果的可靠性。      

俯冲带地震动预测模型SMA2020的评估检验实例——2021年2月13日日本福岛县东部海域MW7.1地震

地震动预测具有较大不确定性,利用强震动观测记录对已有地震动预测模型进行检验评估对于模型的合理使用具有重要参考意义。针对2021年2月13日发生在日本福岛县东部海域的M_W7.1地震,基于残差分析和对数似然函数法对俯冲带地震动预测模型SMA2020的路径项和场地项表征合理性、模型的预测精确性进行了评估。得到结果:（1） SMA2020模型对于地震动的路径衰减表征较为合理,在路径项中考虑了莫霍面反射效应,可以更真实地表征俯冲带深源地震的路径效应。（2） SMA2020模型对于PGV的场地效应表征较为合理,但对于PGA及PSA的场地效应表征有所偏差,还需要进一步检验验证。（3） SMA2020模型对于此次板内地震的短周期IM预测精确性要优于中长周期,对于T=1s的PSA预测表现最不佳。研究结果可为模型的潜在使用者提供有价值的辅助信息,也可以为我国地震动预测模型的开发研究提供理论参考。     

用于核工程地震安全性评价的中小地震水平向加速度反应谱衰减关系研究

利用美国西部的中小地震记录对基岩场地条件下的中小地震水平向地震动多阻尼加速度反应谱的衰减关系进行了初步研究,并分析了震级与距离对中小地震水平向地震动加速度反应谱衰减规律的影响.