龙门山地区的kappa(k_0)模型及汶川M_s8.0地震的强地震动模拟

在地震动模拟中,高频衰减参数(κ_0)是一个重要的参数,它控制了傅里叶谱高频部分的衰减特性.本文利用汶川M_S8.0地震和芦山M_S7.0地震主余震的1671组强震动观测记录,计算了龙门山地区50个断层距小于150 km的强震动台站的κ_0,基于随机有限断层法模拟了汶川地震中这些台站的加速度时程、傅里叶振幅谱和反应谱,并与前人的研究进行比较.结果表明,合理计算的κ_0可以有效地改善加速度时程振幅和高频谱(1 Hz)的模拟结果.另一方面,本文基于κ_0与地形高程的相关性,建立了龙门山地区的κ_0经验模型.分别采用该经验模型和κ_0=0.04 s模拟了汶川地震的峰值加速度分布,并与观测记录进行比较.结果表明,采用本文提出的κ_0模型可以更好地重现汶川地震的峰值加速度分布,特别是在断层破裂方向的反方向区域和山区.综上,汶川地震中山区的峰值加速度明显大于盆地地区的现象,不仅与断层破裂产状有关,还与山区和盆地地区的κ_0之间显著的差异有关.     

基于GP14.3运动学混合震源模型和SPECFEM 3D谱元法的宽频地震动模拟

基于确定性物理模型的全过程地震动模拟是现代地震工程的重要发展方向.然而受限于合理震源模型和计算资源需求,目前模拟的有效频率还多处于低频范围,难以满足工程结构敏感频带(5～10 Hz或更高)需求.本文即借助运动学混合震源模型能激发宽频地震波和谱元法空间高精度及计算收敛快的优势,首先将确定性的凹凸体震源模型与GP14.3随机震源模型结合得到有限断层运动学混合震源模型,进而将上述混合震源模型开发到SPECFEM 3D谱元法开源代码中,实现了基于谱元法和运动学混合震源模型的全过程宽频带地震动模拟.将方法首先应用于一维波速结构模型0～10 Hz地震动模拟,通过与频率波数域(FK)方法结果进行比较,验证了方法的精度;进而应用于2021年5月21日云南漾濞6.4级地震0.1～5 Hz地震动模拟,通过与4个台站的时程记录和相应反应谱的比较,以及与NGA-West2地震动衰减方程在频率0.1～5 Hz的反应谱的比较,检验了方法的适用性;最后给出了漾濞地区的地震动峰值加速度(PGA)和峰值速度(PGV)云图,分析了漾濞地震下近场强地面运动的空间分布特征.结果显示,震中PGA接近400 cm·s^-...     

龙门山地区的kappa（κ0）模型及汶川MS8.0地震的强地震动模拟

在地震动模拟中,高频衰减参数（κ₀）是一个重要的参数,它控制了傅里叶谱高频部分的衰减特性.本文利用汶川M_S8.0地震和芦山M_S7.0地震主余震的1671组强震动观测记录,计算了龙门山地区50个断层距小于150 km的强震动台站的κ₀,基于随机有限断层法模拟了汶川地震中这些台站的加速度时程、傅里叶振幅谱和反应谱,并与前人的研究进行比较.结果表明,合理计算的κ₀可以有效地改善加速度时程振幅和高频谱（>1 Hz）的模拟结果.另一方面,本文基于κ₀与地形高程的相关性,建立了龙门山地区的κ₀经验模型.分别采用该经验模型和κ₀= 0.04 s模拟了汶川地震的峰值加速度分布,并与观测记录进行比较.结果表明,采用本文提出的κ₀模型可以更好地重现汶川地震的峰值加速度分布,特别是在断层破裂方向的反方向区域和山区.综上,汶川地震中山区的峰值加速度明显大于盆地地区的现象,不仅与断层破裂产状有关,还与山区和盆地地区的κ₀之间显著的差异有关.     

2021年5月21日云南漾濞MS6.4地震随机有限断层三维地震动模拟

2021年5月21日云南省大理州漾濞县发生MS6.4 地震,震中附近遭受了强烈地震破坏.为预测此次地震的地震动影响场,利用震源运动学破裂随机模型,基于随机有限断层三维地震动模拟方法,给出了此次地震中 28个触发强震动台站的三分量加速度时程模拟记录,并结合强震动观测记录,估计了此次地震的地震应力降及震源破裂过程,进一步模拟给出了此次地震中 2823 个虚拟观测点的三分量加速度时程.结果表明,模拟记录的峰值地面加速度(PGA)、峰值地面速度(PGV)与观测值接近,并体现了地震动峰值的衰减规律、近场饱和效应和破裂方向性效应;模拟与观测记录的 5%阻尼比拟加速度反应谱(PSA)的幅值接近、谱形相似,在0.05～10 s周期段,模拟记录可以很好地预测地震动.基于三分量模拟记录给出了漾濞MS 6.4 地震的仪器测定地震烈度图,与云南省地震局发布的烈度图接近,极震区烈度最高可达Ⅷ度,震源破裂方向性导致震中 SE方向的烈度普遍高于 NW方向,受局部场地条件影响沿洱海西侧出现高烈度异常区.     

2022年9月5日泸定MS6.8地震宽频带地震动谱元法模拟

2022年9月5日12时52分四川甘孜州泸定县发生M_S6.8地震,此次地震造成泸定县及其周边地域的严重破坏和人员的重大伤亡。为重现此次地震的地震动影响场,分析近场强地面运动的空间分布特征,将确定性的凹凸体震源模型与随机震源模型结合得到有限断层运动学混合震源模型,进而将上述混合震源模型开发到SPECFEM 3D谱元法开源代码中,实现了基于谱元法和运动学混合震源模型的泸定M_S6.8地震的全过程宽频带(0.1—5 Hz)地震动模拟,通过与六个台站的时程记录、对应的反应谱以及NGA-West2地震动衰减曲线进行比较检验了方法的精度和适用性;进而给出了此次地震的三分量速度波场快照图,再现了地震波传播时近场地震动的方向性效应和局部场地效应;最后给出了震中100 km×100 km范围内的地震动峰值加速度(PGA)和峰值速度(PGV)云图,分析了泸定地震近场强地面运动的空间分布特征,并给出了基于模拟结果的地震烈度分布图。结果显示,震中PGA接近600 cm/s²,PGV接近50 cm/s,烈度达到Ⅸ度,且由于泸定地区内高山峡谷地形对地...     

一种随机有限断层的三维地震动模拟方法——鲁甸地震为例

传统随机有限断层方法只能模拟水平单向地震动,不能完全满足工程抗震设计对三维(两个水平向及竖向)地震动需求。文中在传统方法的基础上发展了一种可模拟三维地震动的随机有限断层方法,利用水平向和/或垂直向P、SV、SH波在频域上的经验表达并结合随机相位合成三维地震动加速度时程。为验证方法的可靠性,模拟了2014年云南鲁甸地震的三维地震动。结果表明,改进方法合成的加速度时程波形与观测记录较相似,且模拟记录PGA、0.05～10 s周期段PSA幅值和谱形与观测记录也具有良好一致性,能合理反映地震动随距离的衰减规律;鲁甸地震三维地震动的可靠模拟也验证了路径衰减、路径持时经验模型、高频衰减参数、震源谱衰减率等模拟输入参数的可靠性,可进一步应用于该区域设定地震地震动预测。     

新疆北天山地区宽频带地震动模拟关键技术研究

近年来,我国强震动观测记录数据库有了大幅扩充,一定程度上满足了重大工程抗震设防对宽频带地震动输入的迫切需求。然而,对于地震频发、强震动观测记录少的地区,宽频带地震动的数值模拟仍是确定合理地震动输入的重要途径。北天山地区是新疆经济社会发展的核心区域,现有强震动观测记录难以满足该区重大建设工程的抗震设防需求。研究构建考虑北天山地区震源机制、地质构造以及场地特征的宽频带地震动模拟方法对于高地震风险的北天山地区具有重要工程应用价值。宽频带地震动模拟方法由具有强区域特征的高频地震动随机模拟方法及基于求解无限域波动问题的低频地震动模拟方法组成。本文围绕北天山宽频带地震动模拟方法构建中的关键问题展开研究:针对高频地震动模拟,主要研究中小震震源模型的构建及非平稳地震动合成所需的群速度模型;针对低频地震动模拟,主要研究无限域波动数值模拟中高精度人工边界的构建,以尽可能缩小模拟工作区和提高计算效率。取得的研究成果如下:(1)北天山地区震源和路径传播模型构建。利用该地区2009—2015年843个M_L≥3.0地震记录,采用多台多地震联合反演方法建立了该地区的路径衰减、台站场地响应以及震源谱模型。结果表明,北天山地区中小震震源谱具有分段特性:4级以下地震的震源谱与引入震源谱高频衰减因子(κ=0.003)的Brune模型在统计意义上基本一致,4级以上地震的震源谱逐渐偏离Brune模型,呈现较为明显的双拐角频率特性;另外,震源谱存在分区差异,库拜盆地、塔里木盆地内地震的震源谱在高频段衰减更为显著。(2)北天山地区等效群速度模型构建。在已有基于等效群速度模型合成地震动的工作基础上,引入震相微分计算方法改进了等效群速度的计算方法;提出应用测震记录P波到时构建最快波群速度模型,分析了最快波群速度选取对等效群速度模型构建的影响;基于北天山地区的强震动观测记录,统计得到该地区的等效群速度模型。等效群速度模型引入到高频随机模拟方法的算例结果表明,基于等效群速度模型考虑地震动非平稳特征可给出更为合理的地震动模拟结果,对于面波发育的深厚覆盖层场地影响尤为突出。(3)高精度人工边界条件完美匹配层构建。完美匹配层内场方程和界面条件通常分别采用复坐标延伸技术变换强形式无限域内波动方程和界面条件得到,其构建过程相互独立,可能出现匹配不合理而引发的数值失稳、计算精度低下等问题。为此提出采用复坐标延伸技术变换弱形式无限域波动方程以构建完美匹配层的方法。弱形式波动方程耦合了波动方程及界面条件,规避了变换后所得场方程与界面条件的匹配不合理问题。新方法可直接给出弱形式匹配层,结合勒让德谱元建立了弹性介质近场波动谱元模拟方案。利用算例验证了新方案的精度及数值稳定性。(4)新源—和静6.6级地震地震动场的数值模拟。2012年6月30日新源—和静地震是北天山地区近50年来发生的最大地震,测震及强震动观测记录均较为丰富。本文应用反演所得滑动分布和随机滑动分布分别拟合了强震动台站的地震动时程,对比分析了拟合反应谱与观测记录反应谱的差异性;模拟了该地震的烈度分布图,分析其与实际调查烈度图和基于新疆常用地震烈度衰减关系计算的烈度图之间的差异,探讨了破坏性地震发生后应用高频随机有限断层方法生成地震烈度分布图的可行性。     

兰州观象台存放台阵汶川M_S8.0强震动记录与分析

兰州观象台存放台阵在汶川特大地震中获取了10组三分量数字强震动加速度记录.本文依据这些实际观测资料研究了黄土场地在地震发生后的自由场地面运动加速度过程及其特征.结果表明在震源、场地和震中距相同情况下,场地上不同测点记录的加速度峰值存在一定程度的差异,较明显地反映了场地地形对EW分向地面自由场加速度的影响.10个测点的峰值加速度(PGA)均显示EW向大于NS向,而垂直向最小;各测点的傅里叶谱值,尤其是在小于1 Hz频段范围内,差异极小.     

傅里叶谱和反应谱之间的关系及其对经验地震动预测方程(GMPE)调整的影响研究(Ⅰ)

经验反应谱地震动预测方程(GMPEs)的函数形式通常由地震动傅里叶谱模拟的概念导出,随后这些GMPEs由经验观测数据校准,所以,对特定地震场景的地震动预测不会构成重大问题。然而,当调整反应谱GMPEs来计算未被原始经验数据集覆盖的条件时,傅里叶谱模拟导出反应谱的假设可能会产生意想不到的结果。因此,几个问题出现了。例如,地震动傅里叶谱和反应谱的区别和相似之处是什么?如果它们是不同的,什么机制可以解释这样的差异?对傅里叶振幅谱(FAS)所做的调整在反应谱中怎么表示?本文利用随机振动理论(RVT)探讨了地震动傅里叶谱和反应谱之间的关系。借助简单的Brune震源模型~([1-2]),在固定震级和距离的情况下,用RVT生成了加速度谱。RVT分析表明,反应谱低频谱值标定可以视为与相应傅里叶谱标度值等同。然而,反应谱高频谱值受到了带宽很大的傅里叶谱的控制。实际上,地震动峰值加速度不能反映高频地震动特征,它受到了整个地震动傅里叶谱的控制。此外,本文说明了对FAS做的调整如何相似或区别于对反应谱坐标所做的相同的调整。为此,我们研究了调整应力参数(Δσ)(震源项)和调整反应场地响应的属性(V_S-κ_0)这两种情况。     

随机有限断层法模拟中强地震近场地震动的研究及应用

核工程是对地震安全性要求极高的工程,在其工程场地地震危险性评价工作中,弥散地震地震动估计往往对厂址设计基准地震动具有关键性影响。然而,弥散地震不是由特定构造引起的,因此,评定其地震影响时会遇到较大困难。常规方法是采用适当的衰减关系来计算工程场地震级5.0和5.5的弥散地震在距离厂址5 km处的地震动参数。本文针对核电厂地震危险性评价需求,探讨弥散地震评价的新方法,即通过随机有限断层法来模拟中强地震近场地震动参数。通过参数敏感性分析技术筛选了对拟合结果有决定性影响的关键模型参数;研究并改善了有限断层模型拟合中强地震模型应力降参数的方法,并基于多个国家和地区中强地震近场加速度记录估算了应力降参数值;此外,通过建立中强地震类属模型来模拟特定震级-距离的地震动,从而对经验统计方法进行补充和验证。(1)详细研究了模拟地震动时程的随机有限断层法以及近年来对该方法的改进,改进后的随机有限断层法适合模拟中强地震。(2)对随机有限断层法模拟中强地震近场强震动的参数影响情况进行了分析。将随机有限断层法应用于中强地震时,由于震源信息的准确性较差从而使模型参数具有较大的不确定性。为此,本文首先比较了不同场地方位角中强地震近场地震动的时程和平均拟加速度反应谱(PSA),并定量分析了模型主要参数对中强地震近场地震动模拟结果的参数敏感性。结果表明:不同场地方位角的中强地震近场PSA在短周期部分差别较大;应力降是模型最重要的参数,它对反应谱短周期部分影响最大;几何扩散系数对PSA的整体影响也较为明显。将随机有限断层方法应用到工程安全性评价工作时,应当重点关注对反应谱短周期部分影响较大的应力降和该区域的几何扩散系数,同时要调查该区域优势场地方位角的分布,从而更合理地控制中强地震近场强震动的模拟。(3)开展了中强地震模型应力降参数值的估算。基于美国小头骨山(Little Skull Mountain)MW5.6地震两个近场台站记录的地震动模拟,详细研究并改善了采用随机有限断层法拟合中强地震地震动拟加速度反应谱(PSA)来确定应力降参数值的方法。结果表明:采用不同频段反应谱残差和计算得到的应力降值差别较大,确定中强地震应力降较为合适的反应谱频段是中高频,采用该频段确定的应力降参数值模拟的反应谱和峰值加速度与实际记录较为符合。基于该方法并根据美国西部、日本、意大利以及汶川等国家和地区的浅源中强地震近场加速度时程记录,估算了相应的模型子断层应力降参数值。结果表明:各个国家和地区计算得到的中强地震随机有限断层模型应力降参数具有地区差异性;模拟的峰值加速度大致上正确反映了实际记录PGA的衰减规律,并且模拟结果在近场随距离的衰减较快,这可能与计算所使用的中强地震高频成分丰富有关。(4)针对中强地震震源及路径参数难以准确确定的情况,提出了通过建立随机有限断层类属模型参数来考虑随机不确定性,从而进行中强地震近场地震动模拟的方法,即在考虑区域特定参数、选择具有代表性的震源参数的同时,考虑关键参数随机不确定性,开展大量模拟计算。完成了对核电厂工程弥散地震地震动参数评价的工程应用分析。对模拟PGA的统计结果表明:考虑核电厂弥散地震时,MS5.0地震震中距5 km处0.15g的峰值加速度是合适的;MS5.5地震震中距5 km处0.2g的峰值加速度是合适的;MS6.0地震震中距5 km处0.3g的峰值加速度是合适的。该结果为现有弥散地震地震动评价提供了验证信息。模拟PSA统计结果与衰减关系PSA的比较表明,模拟的平均PSA大致上是保守的,可以作为随机有限断层法模拟结果的推荐反应谱。     

2015年4月25日尼泊尔地震波场传播及烈度初步模拟分析

2015年4月25日尼泊尔发生Mw7.9地震.本文采用三维曲线坐标网格有限差分方法,依据USGS给出的震源运动学反演结果和震源区域地形数据,模拟了尼泊尔地震波场传播过程,得到震源区域烈度分布模拟结果,与实际观测的地表峰值速度(PGV)大体吻合.结果表明:地震烈度的空间分布整体上受控于震源的单边破裂特征,高烈度区域主要分布在破裂传播方向上,即震源东半部.震源南侧到东南侧近场,由于受到震源和地形的双重影响,形成最大烈度分布区域,最大烈度约为IX.南侧平原受低速沉积层影响形成高烈度区域.震源西侧及盆地内烈度相对较低.     

Testing the Validity of Simulated Strong Ground Motion from the Dynamic Rupture of a Finite Fault, by Using Empirical Equations

This paper is concerned with testing the validity of the ground motions estimated by combining a boundary integral equation method to simulate dynamic rupture along finite faults with a finite difference method to compute the subsequent wave propagation. The validation exercise is conducted by comparing the calculated ground motions at about 100 hypothetical stations surrounding the pure strike-slip and pure reverse faults with those estimated by recent ground motion estimation equations derived by regression analysis of observed strong-motion data. The validity of the ground motions with respect to their amplitude, frequency content and duration is examined. It is found that the numerical simulation method adopted leads to ground motions that are mainly compatible with the magnitude and distance dependence modelled by empirical equations but that the choice of a low stress drop leads to ground motions that are smaller than generally observed. In addition, the scatter in the simulated ground motions, for which a laterally homogeneous crust and standard rock site were used, is of the same order as the scatter in observed motions therefore, close to the fault, variations in source propagation likely contribute a significant proportion of the scatter in observed motions in comparison with travel-path and site effects.     

近场强地震动数值模拟的简化计算方法

近场强地震动除受场地条件的影响外,还受到震源破裂面上子源的空间分布特点、子源破裂先后顺序的强烈控制,基于数值格林函数法的近场强地震动数值模拟方法可以综合考虑震源、传播途径及局部场地条件的影响,对计算过程进行合理简化,分2步完成地震动模拟:第1步,在介质均匀区采用矩张量的解析解计算所有子源在盖层底面的位移,形成下一步有限元计算的输入场;第2步,在盖层介质不均匀区,结合局部人工透射边界技术,采用时、空解耦的波动显式有限元方法计算地表强地震动。在有限断层模型中,采用具有9个力偶的等效地震矩张量表达断层产状、滑动方向等的影响,采用Brune模型定义各子源的滑动时间函数,描述滑动的时、空不均匀分布特征,从而细化震源模型。通过对Northridge地震中4个基岩台站地表地震动的模拟结果和强震记录,验证了此简化计算方法的可行性     

Spatial distribution of near-fault ground motion

Introduction In order to explain severe damage of structures and ground surface near earthquake fault, many theoretical and numerical methods have been developed to study strong ground motion near the fault (Somerville, 1998) via the kinematic and dynamic source models. Moreover, a great number of valuable strong ground motion records, especially near the fault, have been acquired in recent years. These valuable records provide useful data for improving simulation approaches. Based on those r…     

断层附近地面地震动空间分布

运用震源位错模型, 分析矩形垂直断层及倾斜断层走向滑动和倾向滑动的近场地震动场，以地表地震动的傅立叶振幅谱比为参量考察断层附近地震动空间分布的特点. 结果表明, 断层附近的地震动强度主要受近旁子断层的控制，高强度的地震动分布在紧靠断层两侧有限的带状区域内，长周期分量受断层破裂传播方向性的影响. 走向滑动的方向性影响主要表现在垂直于断层走向的分量，倾向滑动则表现在平行于断层走向的分量，且深震在地面上引起的地震动强度分布比浅震要平缓，影响范围宽. 倾斜断层产生的地震动有明显的上盘效应，空间分布不对称，与观测结果相符. 最后给出了断层附近近场地震动强度分布拟合函数的表达式，并与美国的1997统一建筑规范规定的近场因子作了比较.      

特征地震强地面运动参数衰减公式拟合

计算机仿真模拟设定地震断层动态破裂传播和近断层强地表运动响应的结果表明， 对于特征地震而言，近断层附近的地表运动特征与断层破裂传播的方向性有着强烈的依赖关系. 当场地（观测点）至断层的距离给定时，正向于破裂传播方向的场地（场地A）的地表质点运动（位移、速度、加速度），远远大于震中附近（场地B）和反向于破裂传播方向的场地（场地C）的地表质点运动，而且沿断层垂直分量所辐射的SH波的传播起到了主导作用. 对应于场地A，B和C，统计分析结果表明，峰值加速度的几何平均值之比为2.15:1.5:1， 而且各自的均方差分别为0.12, 0.11和0.13. 如果将所得的研究结果应用于概率地震危险性分析中，对于较低的年超越频度，近断层附近的地表峰值加速度的估算值可下降15%~30%. 因此，考虑到断层破裂传播方向性对地表运动的影响，区域衰减曲线的回归分析模型应该给予恰当的修正.      

近断层强地震动场预测

以1997年4月11日新疆伽师地震（M_w6.1）为例,详细介绍了近断层强地震动场的预测方法.首先,用有限断层震源建模方法建立了该次地震的震源模型;然后,基于动力学拐角频率的地震动随机模拟方法,模拟了该次地震仅有主震加速度记录、且位于巨厚土层上的三个台站的加速度时程,并用实际地震记录进行了验证.在此基础上,基于预测的近断层77个节点的加速度时程的峰值绘制了该次地震的加速度场.结果表明,上述方法模拟的加速度时程在0.5 Hz以上的高频段是可行的、实用的;预测的近断层加速度场具有非常明显的上盘效应.地表最大加速度的范围与断层面上最大凹凸体位置相对应,说明与断层面上凹凸体相对应的地面上的建（构）筑物将会遭受到较为严重的震害.     

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

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

Kinematic source model for simulation of near-fault ground motion field using explicit finite element method

1 Introduction The prediction of broadband near-fault ground motion is an important aspect of seismic risk analysis, and is essential for the evaluation of urban seismic safety. Many methodologies are used to model the near-fault ground motions, such as the Green function method, ?nite element (FE) / ?nite difference (FD) method and a hybrid method (Aagaard, 1999; Oprsal and Zahradnik, 2002; Oprsal et al., 2003 a and b). In the hybrid method, the FE / FD method is used to model the low fr…     

Rupture directivity and hanging wall effect in near field strong ground motion simulation

A random synthesis procedure based on finite fault model is adopted for near field strong ground motion simulation in this paper.The fault plane of the source is divided into a number of sub-sources,the whole moment magnitude is also divided into more sub-events.The Fourier spectrum of ground motion caused by a sub-event in given sub-source,then can be derived by means of taking the point source spectrum,attenuation with distance,energy dissipation,and near surface effect,into account.A time history is synthesized from this amplitude spectrum and a random phase spectrum,and being combined with an envelope function.The ground motion is worked out by superposition of all time histories from each sub-event in each sub-source,with time lags determining by the differences between the triggering times of sub-events and ddistances of the sub-sources.From the example of simulations at 21 near field points in a scenario earthquake with 4 dip angles of the fault plane,it is illustrated that the procedure can describe the rupture directivity and hanging wall effect very well.To validate the procedure,the response spectra and time histories recorded at three near fault stations MCN,LV3 and PCD during the Northridge earthquake in 1994,are compared with the simulated ones.