随机有限断层法及其工程应用中的问题分析

随机有限断层法作为半经验半理论的地震动合成方法,综合考虑了震源、传播路径以及场地条件对地震动的影响,可在工程关注的频率范围内模拟地震动时程,为实现较大区域地震动模拟提供了一种有效的方法。本文系统阐述了随机有限断层法的基本原理、静力学拐角频率模型及动拐角频率模型的发展;分析了主要模型参数（断层尺度、破裂速度、剪切波速及应力降）、Kappa因子和介质品质因子的取值原则,并结合中国大陆的实际情况给出了相应的取值范围;最后讨论了现有模拟方法在考虑参数取值、子断层划分和子断层之间相互作用等方面的不足以及能量处理方式上的缺陷,明确了随机有限断层法发展的方向。     

汶川地震水平向地震动高频衰减参数的统计分析

高频衰减参数kappa是近场强地震动模拟中的一个重要参数,本文利用汶川地震的强震记录计算了主余震的kappa值,并分别给出了3种场地类别的kappa值随震中距和震级变化的统计参数。统计结果表明,水平方向地震动的高频衰减较明显;当震中距或震级较小时,部分强震记录的衰减缓慢;3种场地类别的kappa值统计结果均与震级和震中距有关,并随着场地类别的增大,相关系数相应增加。计算得到的kappa值离散性很大,本文在线性拟合的基础上给出了kappa的取值范围,这一范围基本控制在拟合直线两倍标准差之间。本文的研究结果对kappa值的估计有重要参考价值。     

考虑椭圆模型MW≥5.5地震震源参数时域计算方法

为提高时效性,地震预警系统对震源的处理一般采用“点源”假定,不考虑震源尺度、破裂方向和震源区应力降,对震源参数简化处理使得预测地震动参数分布过于粗略,特别是对破坏性大震.本文应用日本地震(M_W≥5.5)强震动数据,考虑椭圆震源模型,提出了一种适用于中、大地震矩震级、拐角频率和应力降等震源参数时域计算方法,并分析了场地效应、地震动传播衰减、震源效应和滤波频带等因素对时域方法的影响.结果表明,该方法能利用P波信息快速获取震源能量释放过程,量化随破裂传播震源谱拐角频率、应力降和破裂面积的演化过程,为地震动预测提供更多重要的震源参数.     

汶川地震强地面运动模拟

基于确定性震源模型的方法主要用于计算低频(11 Hz)地震动常用经验格林函数法或随机方法,对低频地震动模拟不够准确．本文在确定性震源模型方法基础上,尝试采用分解给定的震源模型的方法来模拟宽频带(0.1——10 Hz)强地面运动,即采用分级离散断层面和分解断层面破裂单元上升时间的方法,增加震源时间函数中的高频信号,从而避免了对地震记录丰富程度和准确性的依赖．文中模拟计算了汶川MS8.0地震在8个地震观测台的地震动,将模拟结果和观测记录进行了加速度时程曲线和傅里叶振幅谱的对比分析．对比结果显示,模拟估计的地震动峰值加速度和持续时间与观测记录的数据基本在plusmn;50%的精度范围内相同,傅里叶振幅谱显示模拟结果有得到10 Hz左右的高频成分. 四川盆地中的台站模拟结果高频衰减比观测记录要快,原因是模拟过程没有考虑场地效应.对强地震动模拟还是要综合考虑震源、传播路径和场地的影响．研究结果表明,此改变震源输入的确定性方法可应用于模拟近断层宽频强地面运动．      

地震动随机合成中与震源谱相关的动力学拐角频率

介绍了基于有限断层震源的地震动随机合成中静力学拐角频率和动力学拐角频率,讨论了两者的缺点;进而提出了一种改进的震源谱模型,它能够表达破裂面上地震波频率辐射不均匀性、大震时拐角频率随破裂面积增加而有所下降的趋势．据此,可以形成与该震源谱模型相关联的动力学拐角频率,应用于地震动的随机合成,避免结果对子源尺寸的依赖．通过对美国北岭Mw6.7地震近场6个基岩台站合成的地震动与实际记录的比较,验证了本文方法的有效性.      

2013年乌鲁木齐MS5.6和MS5.1地震强地震动模拟研究

利用基于动力学拐角频率的随机有限断层法,针对2013年3月29日和8月30日乌鲁木齐发生的M_S5.6和M_S5.1两次地震,选择不同机构的震源参数对23个台站的45条强震记录进行强地震动模拟,对比分析加速度时程和反应谱,并计算模型偏差.通过与实际场点记录进行对比的结果显示:选取的不同机构震源参数模型得到的模拟加速度时程在持时和形状上与实际记录有一定差距,对于近震源的台站,模拟结果的峰值加速度(PGA)比观测结果小,对于稍远的台站,结果基本能够保持一致;而加速度反应谱,模拟结果与观测结果具有较好一致性; 随机有限断层模型的误差在±0.5以内,拟合效果较好,而且高频段的拟合效果要好于低频段;采用不同机构震源参数模型得到的模拟结果与观测结果得到的PGA分布特征较为一致,但模拟结果的峰值加速度要低于观测结果.      

地震动随机模拟方法主要影响参数分析

简要介绍了地震动随机模拟方法的原理,总结了影响地震动随机模拟结果的输入参数,重点分析了应力降、品质因子、高频衰减参数对PGA和反应谱Sa的影响。研究结果表明:应力降对地震动模拟结果影响最大;随着应力降的增加,PGA增大,各周期处的反应谱均有较大变化;矩震级增大时,PGA和反应谱的受应力降的影响也相应增大;品质因子Q0和指数n越大,PGA也越大,短周期处的反应谱受两者影响大,长周期处受其影响小;高频衰减参数k增大时,PGA和反应谱相应减小,k对PGA和短周期处的反应谱影响较大。     

考虑区域特征的地震动模拟 ——以2020年伽师MS6.4地震为例

为研究伽师地区地震动的区域特征,提高地震动预测的可靠性,文中利用伽师地区近10a收集的强震动观测记录,采用非参数谱反演方法分离了地震动的震源、传播路径及场地影响,建立了反映伽师地区地震动区域特征的路径衰减、路径持时、场地线性反应等经验模型,并给出了伽师地区46次地震的震源参数(矩震级、拐角频率和应力降).文中确定的伽师...     

设计地震加速度图的合成

本文提出了一种模拟发震断层附近地震动时程的半经验模型。如给定在该断层上可能发生的地震震级和场地位置,则可由此模型用计算机提供一组设计地震加速度图。模型中提供的经验格林函数计算公式是用所研究地区的强震记录按统计方法在频域内建立,而不采用在时域内修正个别小震记录的方法。模型中包括的其他重要参数,如断层破裂长度和震源谱,可以直接从该地区的强震地面运动的经验衰减规律得到。用本文方法综合的地震加速度时程是与该地区地球介质内地震波的频散和衰减的统计特征相一致的非平稳过程,并能考虑场地相对于震源的方位和震源破裂方式的影响。由于未来的破坏性地震在发震断层上确切的破裂区段一般不能预测,就重大工程场地的抗震设计而言,可以按本文方法提供一组具有不同方位角的设计地震加速度时程。     

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方向,受局部场地条件影响沿洱海西侧出现高烈度异常区.     

天津地区场地响应和震源参数联合反演

根据天津数字地震台网记录的发生在天津地区的中小地震波形资料,采用Moya方法,反演天津地区39个地震的震源谱低频水平、拐角频率和28个台站的场地响应,计算了参与反演的地震矩、拐角频率、应力降和震源半径等参数,并讨论了它们之间的关系.结果表明: 27个井下台的场地响应曲线形态基本一致,低频有不同程度的放大,高频衰减较快;蓟县地表基岩台的场地响应低频基本在1附近变化,高频也有明显衰减.地震矩与震级及震源谱拐角频率的关系密切;拐角频率随震级的增大,有明显的下降;对于3级以上的地震,应力降、震源半径随地震震级的增大有明显的增大.     

强震动记录选取中目标谱的研究及应用

随着强震动观测技术的发展,强震动观测数据的累积推动了结构抗震设计理论的发展。目前,在许多国家的抗震规范中都将结构动态时程分析方法列为静态设计方法的必要补充,但对于高层建筑、大跨结构等特殊结构形式及重要建设工程,要求强制使用动态时程分析方法。我国结构抗震设计规范及相关标准关于强震动记录选取的规定过于笼统,结构设计人员在面对大量强震动记录数据时缺乏明确的选取依据。本文搜集整理了国内外现有强震动记录数据,通过对地震参数的统计,提出基于场地类型的改进反应谱谱值相关系数经验公式,并将其应用于条件均值反应谱(conditions mean spectrum,CMS)的建立。同时,本文基于地震动衰减关系建立地震信息与抗震设计参数的对应关系,将其作为匹配规范标准谱选波方法的筛选依据,进而明确了抗震规范中强震动记录选取的初选条件。研究中发现,规范标准谱不宜作为谱形匹配选取强震动记录的目标谱,因此,提出两种建立合理目标谱的方法:一是基于地震动衰减关系建立PGA目标谱;另一是基于地震参数建立EPA均值目标谱。最后探讨了在地震安全性评价工作中实际强震动记录的选取及应用。本文主要研究内容及结论如下:(1)通过对20 000余条强震动记录的地震动参数进行分析,发现记录的加速度、速度及位移峰值和反应谱谱值的相关系数的分界周期点随着场地条件而变化,并提出基于场地类别改进反应谱谱值相关系数经验公式。改进后的反应谱谱值相关系数经验公式改善了目标谱的场地特征,使其更符合场地条件的要求。(2)针对抗震设计规范中关于强震动记录选取的相关条款,提出强震动记录选取的基本原则和具体方法,为完善相关规定提供参考依据。针对一般建设工程,依据地震动加速度峰值(PGA)衰减关系,提出同等设防PGA目标值对应的地震参数范围,建立断层距和地震设计分组的对应关系,进而确定强震动记录选取的初选条件。讨论了结构响应和强震动记录反应谱谱值的关系,以典型框架结构为例确定匹配周期范围。结构位移反应与强震记录反应谱相关性分析表明,两者在结构自振周期前后相关关系最为显著,距离自振周期越远的反应谱值与结构反应的相关系数越小。基于以上分析,建立了选取强震记录的基本原则:以距离和场地条件作为强震动记录选取第一限制条件。根据结构自振周期确定匹配规范标准目标谱周期范围,建议的匹配周期范围为0.4T1~1.5T1,最后由实际记录与规范标准谱在匹配周期范围内的均方差由小到大的顺序选取记录,记录的放缩系数限定在0.2~4.0之间。(3)本文基于地震动衰减关系,建立PGA目标谱并匹配目标分布的强震动记录的选取方法。对于未进行地震安全性评价,无法提出地震参数的一般建筑工程拟建场地,先利用PGA衰减关系,确定同等PGA目标值下的震级和距离,再代入加速度反应谱衰减关系,确定与PGA目标值对应的反应谱即为PGA目标谱。PGA目标谱不仅是确定的均值谱,而且具有概率含义。PGA目标谱的分布及概率特性反应了地震动的不确定性,进而为结构反应的离散性提供了一个合理的参考范围。(4)本文按场地条件和断层距将现有强震动记录数据分组,建立EPA均值目标谱。通过将EPA均值谱、PGA目标谱和规范标准谱进行比较,EPA均值谱与PGA目标谱谱形相似均具有良好谱形,匹配效果优于规范标准谱。建议规范标准谱作为强震记录选取的目标谱时,需要:(1)增加特征周期;(2)减小地震设计分组1、2组,场地类型为I类场地和II类场地的中长周期段反应谱值。(5)探讨了地震安全性评价工作中强震动记录的选取,结合不同等级的工程场地地震安全性评价,提出与之相适应的强震记录选取的基本方法及原则,建议如下:(1)I级安评的工程场地,采取基于设定地震的强震动记录选取方法。首先根据设定地震确定地震参数的初选条件,然后以场地设计谱作为目标谱,匹配场地设计谱选取记录。(2)II级和III级安评的工程场地,采取建立条件均值目标谱的方法。即以场地设计谱上结构特征周期点对应的加速度反应谱值为条件,利用设定地震的震级和距离等基于加速度反应谱衰减关系建立条件均值目标谱,并采取匹配条件均值目标谱分布的方法选取强震动记录。(3)Ⅳ级安评的工程场地,根据目标PGA值,采取PGA目标谱的方法选取强震动记录。     

2010年4月14日玉树Ms7.1地震加速度场预测

基于有限断层震源、且使用动力学拐角频率的地震动随机模拟方法预测玉树地震近断层的加速度场.首先,基于有限断层震源建模方法建立该次地震的震源模型;然后,基于上述地震动模拟方法预测玉树地震近断层191个节点的加速度时程.在此基础上,取每个结点的加速度峰值绘制该次地震的近断层加速度场.结果表明:(1)近断层加速度场主要受震源破裂过程和断层面上滑动分布的影响.断层面上凹凸体投影到地表的区域附近,加速度峰值最大,也是震害最严重的区域;(2)对于走滑地震,断层沿线附近的场地并非均会发生破裂方向性效应;发生破裂方向性效应的场地与凹凸体在断层面上的位置有关.     

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

核工程是对地震安全性要求极高的工程,在其工程场地地震危险性评价工作中,弥散地震地震动估计往往对厂址设计基准地震动具有关键性影响。然而,弥散地震不是由特定构造引起的,因此,评定其地震影响时会遇到较大困难。常规方法是采用适当的衰减关系来计算工程场地震级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大致上是保守的,可以作为随机有限断层法模拟结果的推荐反应谱。     

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

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

滇西南地区品质因子、场地反应与震源参数研究

为深入了解滇西南地区地震动影响因素特征,更好地服务于目标区的抗震设防工作与地震预测等相关研究,利用两步非参数化广义反演方法对滇西南地区2007-2019年242组三分量强震记录进行分析,通过对地表观测地震动在频域内的分离求解,给出了该地区0.5～20 Hz的品质因子经验关系,获得了26个强震台站在研究频段内的场地反应,并与HVSR法揭示的场地放大结果进行对比分析。基于网格搜索法确定了35次地震事件的地震矩M0、拐角频率fc及震源模型参数γ,拟合得到了部分震源参数之间的定标关系。研究结果表明,滇西南地区品质因子Q值与频率的关系为Q=94.23f 0.43;相较于HVSR法,非参数化广义反演方法确定的场地反应整体偏高,但2种方法得到的场地反应谱形较一致;拐角频率、地震矩及应力降等震源参数互相具有一定依赖性。     

云南漾濞6.4级地震强地面运动的模拟和空间分布特征分析

北京时间 2021 年05月21日 21 时48分云南大理州漾濞县(北纬 25.67°,东经 99.87°)发生了 6.4 级地震.本文利用震源破裂过程控制的 NNSIM随机有限断层方法,模拟对比了近场 6 个强震台的加速度、速度、反应谱记录,据此确定了开展近场地震动模拟所需要的参数的大小,进而建立了近断层 200 km范围内的地面运动.分析了漾濞地震近场地面运动的空间分布特征,结果显示此次地震地面运动的峰值速度PGV、峰值加速度PGA以及反应谱分布均表现为圆形分布,无明显的上下盘特征和走向特征,震中极震区峰值加速度超过了 400 cm·s－2 ,对应国家标准烈度Ⅸ度.     

基于随机有限断层方法的唐山地震强地面运动模拟

基于动力学拐角频率的随机有限断层模型,以1976年唐山大地震的极震区为研究目标,通过梳理前人的研究成果,建立了目标区的计算模型并确定了震源模型参数,依靠目标区的局部场地资料和场地动力学参数,利用四分之一波长法,获取局部场地放大效应参数;利用随机有限断层模型地震动模拟方法,给出了唐山大地震近场强地震动场模拟结果,并与历史地震宏观烈度进行对比,结果表明:在Ⅷ、Ⅸ度影响烈度的范围与宏观资料吻合较好。     

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

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

基于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^-...