2023年甘肃积石山MS6.2地震随机有限断层三分向地震动模拟

摘要：2023年12月18日23时甘肃省临夏州积石山县发生MS6.2地震，震中附近造成大量人员伤亡。基于不同数据和方法反演得到的震源破裂过程还有争议，为模拟此次地震的地面运动，分别采用SW倾向与NE倾向的震源破裂模型，利用随机有限断层法模拟三分向地震动。选取震中距约100 km内的57个触发强震台站，模拟其三分向加速度时程、速度时程及5%阻尼比拟加速度反应谱（PSA）。结果表明，模拟的峰值地面加速度（PGA）、峰值地面速度（PGV）能够体现地震动峰值的衰减规律，NE倾向的震源破裂模型会导致更大的地震动。模拟与观测记录的幅值接近、谱形相似，其中SW倾模型的PSA残差更小。基于三分向模拟地震动给出积石山MS6.2地震的仪器烈度分布，SW倾向与NE倾向模型极震区烈度均达Ⅷ度，其中SW倾向模型与应急管理部发布的烈度图更加接近。积石山MS6.2地震的模拟说明随机有限断层法可用于计算三分向地震动，并验证模拟输入震源、路径与场地参数的可靠性。

关键词：积石山MS6.2地震; 随机有限断层法; 地震动模拟; 地震烈度; 强震记录

中图分类号： P319.9文献标志码：A文章编号： 1000-0844（2024）04-0778-06

DOI：10.20000/j.1000-0844.20240311002Three-component ground motion simulation for the 2023

MS6.2 earthquake in Jishishan， Gansu， based on the

stochastic finite-fault methodMAO Lan SHI Yucheng LU Yuxia LIU Kun

（1. Lanzhou Institute of Seismology， CEA， Lanzhou 730000， Gansu， China;

2. Gansu Earthquake Agency， Lanzhou 730000， Gansu， China;

3. Key Laboratory of Loess Earthquake Engineering of CEA amp; Gansu Province， Lanzhou 730000， Gansu， China;

4. Gansu Provincial Research Center for Conservation of Dunhuang Cultural Relics， Dunhuang 736200， Gansu， China）Abstract：

At 23：00 on December 18， 2023， an MS6.2 earthquake occurred in Jishishan County， Linxia Prefecture， Gansu Province. A large number of casualties were found near the epicenter. The source rupture process due to inversion based on different data and methods remains controversial. The source rupture models dipping southwest and northeast were selected to simulate the ground motion of the earthquake， and the stochastic finite-fault method was employed to simulate three-component ground motions. The time-histories of three-component acceleration and velocity and 5% damped ratio pseudo-spectral acceleration （PSA） were simulated by selecting 57 triggering strong motion stations within 100 km of the epicenter. Results show that the simulated peak ground acceleration and peak ground velocity can demonstrate the attenuation law of peak ground motion， and the source rupture model dipping northeast leads to large ground motion. The amplitudes of simulated and observed PSAs are close， the spectral shapes are similar， and the PSA residual of the model dipping southwest is small. The instrument intensity distribution of the Jishishan MS6.2 earthquake is obtained based on the simulated three-component records. The intensities of the two models in the meizoseismal area both reach Ⅷ degree， and the southwest-dipping model is closer to the intensity map released by the Ministry of Emergency Management. Simulation results of the Jishishan MS6.2 earthquake reveal that the stochastic finite-fault method can be employed in the simulation of the three-component ground motion， and the reliability of the simulated source， path， and site parameters can be verified.

Keywords：Jishishan MS6.2 earthquake; stochastic finite-fault method; ground motion simulation; seismic intensity; strong motion records

0引言

据中国地震台网测定，2023年12月18日23时59分，甘肃省临夏州积石山县发生MS6.2地震，震中位于积石山县柳沟乡（35.70°N，102.79°E），震源深度约10 km。此次地震造成老旧农村建筑破损相对严重，大量房屋损毁或倒塌。截至12月22日，此次地震共造成151人死亡，983人受伤。

对地震造成的地面运动进行模拟，是指导抗震设防、制定应急预案、减少地震灾害的迫切需求。地震动模拟的随机方法是模拟设定地震地震动的主要方法之一。随机有限断层法提出以来，经过对拐角频率、高频标定因子等参数计算方法的不断改进［1-4］，目前已广泛应用于设定地震的地震动模拟，如2008年汶川地震［5］、2017年九寨沟地震［6］、1992年开罗地震［7］、2011年东日本大地震［8］。传统随机有限断层法仅能实现S波的水平单向模拟。曹泽林［9］探讨了使用随机法合成三分向地震动的可行性，认为辐射因子是表达方向性的关键参数。王宏伟等［10］通过对P、SV、SH波在频域上的经验表达式并结合高频地震动随机振动特性，实现了使用随机有限断层法模拟三分向地震动，并在2021年漾濞地震［11］和2021年玛多地震［12］得到应用。

不同反演方法和数据得到的地震P波初动节面表明，积石山MS6.2地震可能发生在一条SW倾向的逆冲断层上，也有可能发生在一条NE倾向的逆冲断层上［13-15］。在相同应力降、路径与场地参数的条件下，本文基于不同倾向的震源模型，采用随机有限断层法模拟此次地震的三分向地震动，得到目标场地的加速度时程、速度时程及5%阻尼比拟加速度反应谱（Pseudo-Spectral Acceleration，PSA），并与强震记录对比，验证模拟方法与输入参数的可靠性；同时基于模拟结果给出此次地震的模拟仪器烈度分布，并与台站仪器烈度及正式发布的烈度图进行比较。

1地震概况

应急管理部发布了甘肃积石山MS6.2地震烈度图，如图1所示。由图1可知最大烈度为Ⅷ度，等震线长轴呈NNW走向，长轴124 km，短轴85 km。截至12月24日，此次地震共记录到3.0级以上余震12次，其中3.0～3.9级余震9次，4.0～4.9级余震3次，最大余震为4.1级。震源深度集中在8～10 km，分布呈NNW走向。根据中国地震台网中心发布的震源机制［16］，可知此次地震为逆冲型地震，是青藏高原东北缘拉脊山断裂带上有现代地震记录以来发生的最大地震。

拉脊山北缘断裂位于青藏高原东北缘拉脊山脉的北侧，位于青藏地块、阿拉善地块、鄂尔多斯地块和华南地块的交接处［17］。积石山MS6.2地震发生在青藏高原东缘向东挤出的构造体系内，是近年来我国陆内强震活动的主要控震构造［18］。拉脊山地区地震活动表现出小震频繁而大震稀少的特点，此次地震震中周围50 km范围内自1900年以来没有5级以上中强地震发生［19］。

根据甘肃地震台提供的数据，本文统一对强震记录进行去除仪器响应、基线校正及巴特沃斯0.08～30 Hz滤波处理后，共获得距离此次地震震中100 km范围内57组完整的三分向强震记录，记录时间间隔为0.01 s。其中，记录峰值最大的是距离震中约18 km的JSNGJ台站东西向记录，其EW、NS、UD向的地面峰值加速度分别为-389.4 cm/s2、-329.0 cm/s2、-194.7 cm/s2 。

2模拟方法

3积石山MS6.2地震的模拟参数

随机有限断层法模拟三分向地震动，需要确定震源、路径和场地等关键模拟参数。

3.1震源参数

Zhang等［22］采取中国地震台网中心给定的震源位置（102.79°E，35.70°N;-10 km），GCMT提供的矩张量解，矩震级（MW）为6.0，利用近震资料和远震体波资料反演确定积石山MS6.2地震的震源运动学破裂模型，如图2所示，其中，SW倾向模型［图2（a）］破裂持续时间约为9 s，NE倾向模型［图2（b）］破裂持续时间约为7 s。震源破裂模型显示此次地震为单侧破裂，与余震分布方向一致，主要滑动量位于震中NNW向的破裂面上。

震源至场点的地震波传播距离与传播速度，由速度结构采用试射法计算得出，同时计算得到地震波在震源处的离源角与场点处的入射角。P波与S波的水平成层速度结构根据Crust 1.0 全球地壳模型［28］，在研究区内将不同深度的平均波速作为模拟输入，如图3所示。地震动路径持时采用活动地壳区［29-30］的地震动路径持时经验模型。

3.3场地参数

场地放大与场地地下30 m平均剪切波速（vS30）相关，分为地壳介质和地表岩土层放大两个方面。地壳介质放大模型采用美国国家地震减灾计划中vS30=760 m/s的经验放大效应［31］;地表岩土层放大模型以vS30=760 m/s作为基岩面，根据基岩PGA（Peak Ground Acceleration，PGA）设计场地线性和非线性放大效应［32］。vS30根据基于地形数据的全球估计值确定［33］。

场地衰减表现为当地震波的频率高于一定频率时，加速度谱值会急剧减小，可用kappa滤波器表示。稂子平等［34］研究发现kappa的中位数随vS30增大而减小的趋势明显，采用对数线性函数拟合了kappa与vS30的相关性模型，本文采用此模型计算场点处的kappa值。

综合上述震源、路径、场地信息，设定积石山MS6.2地震模拟的相关参数如表1所列。

4模拟结果及讨论

4.1地震动峰值

地震动峰值是衡量地震强烈程度的重要指标，PGA和PGV反映了地震对地面的冲击力和振动能量的大小。基于不同倾向震源破裂模型模拟的台站三分向PGA、PGV与记录峰值对比如图4所示。在三个方向上，基于SW倾模型和NE倾模型的模拟地震动峰值整体上与记录吻合程度较好，说明几何扩散、品质因子等输入参数能够较好地表达地震动峰值随距离的衰减规律。但记录峰值较强的离散性在模拟结果中未很好地体现出来，这是由于随机有限断层法对模型的简化以及其输入参数难以表达复杂的场地效应和地形效应。部分台站的模拟结果高于记录值，可能是仅由坡度拟合vS30值造成的误差，如基岩台站YDGQS和JYGST，输入vS30值约350 m/s时导致了基岩场地不合理的地表岩土层放大。竖直方向模拟峰值的高估现象更加明显，特别是远场台站的竖直方向。水平成层速度结构计算得到的地震波射线无法体现地层复杂性，模拟输入的射线入射角精度随震中距加大而减小，可能导致竖直向能量分配的误差。

三个方向上NE倾模型整体大于SW倾模型。较SW倾模型，NE倾模型的断层倾角略大，会造成NE倾模型上盘效应不足，达到地面的地震波能量略小；NE倾模型的破裂面积更大，但破裂时间更少，如图5所示，NE倾模型的能量释放更加集中，由此造成更大的脉冲效应。

4.2PSA残差

地震动加速度反应谱反映了地震动的频谱特性，可用来评估和预测结构的地震响应。基于不同倾向震源破裂模型的模拟结果的PSA残差如图6所示。图6反映了PSA谱值在0.05～10 s周期的残差分布，其中PSA残差定义为观测PSA与模拟PSA之比的对数。绝大多数台站两个模型的模拟与记录的PSA残差分布在±1之间。以PSA残差均值是否更接近0为判断模拟效果好坏的标准，三个方向上SW倾模型的整体模拟效果更好。竖直向PSA残差的波动幅度更加明显，可能是由于vS30计算的地表岩土层效应对竖直向放大效果并不理想的缘故。应力降主要影响地震动的高频成分。在相同应力降条件下，与SW倾模型相比，NE倾模型的短周期比中长周期的增量更大，说明使用不同震源破裂模型需要采用不同应力降才能获得符合记录的效果。

4.3地震动时程与反应谱

BX045、JSNGJ、NL001、BX003是距离震中最近的强震仪台站，震中距在30 km以内;BM043位于震中北偏西约39 km;NG001位于震中东偏南约75 km。其中，BX045、JSNGJ、BX003、BM043为基岩台站，NL001、NG001为土层台站。实测与基于不同倾向震源破裂模型模拟的三分向加速度时程、速度时程、PSA如图7所示。模拟与实测波形与谱形相似，整体上谱值相近。加速度时程呈现由弱渐强再渐弱的趋势，速度时程中仍然有明显的脉冲，说明这次地震产生了更丰富的长周期地震动。近场台站的记录波形呈现多峰值、幅值大的特征，且至少有2个明显的脉冲，模拟结果也呈现多峰值的特征，但第一个脉冲并不显著。近场台站波形受震源破裂过程的影响很大，说明实际发震断层破裂的滑动量分布可能存在2个或以上区域集中。NG001作为远场台站，P波与S波的到时差异较大，从地震动时程中也可以明显看出，峰值的到达时间较晚，由此可以区分P波与S波讨论随机有限断层法模拟三分向地震动的优劣。当射线入射角合适时，地震能量在三个方向上的分配合理，模拟结果与记录波形基本吻合。

4.4烈度分布

地震烈度能够直接反映地震对地面和建筑物的破坏程度，地震烈度的分布情况可以为应急救援、灾后重建提供重要依据。基于不同倾向震源破裂模型模拟的积石山MS6.2地震仪器烈度分布如图8所示。仪器烈度分布均呈现出以震中偏北为中心的近圆形分布，极震区具有沿断层走向分布的特点。不同倾向震源破裂模型的烈度分布均呈现显著的上盘效应，符合逆冲型地震的特征。震中附近烈度等值线较为密集，向外随着烈度逐渐减小，变化幅度也较为缓和。破裂面上的滑动分布控制着高烈度区的位置，极震区烈度为Ⅷ度，NE东倾模型的Ⅷ度区面积更大。模拟烈度中出现的较多零散烈度异常区，是随机法模拟地震动和场地效应导致的。在不同倾向模型的模拟中，积石山县的烈度均达到Ⅷ度，民和县的烈度在Ⅵ度到Ⅶ度之间。

远场台站的模拟仪器烈度和记录仪器烈度差异较小，近场台站的模拟仪器烈度较大。近场烈度主要由PGV决定，而图4显示近场台站的模拟PGV整体高于记录值。局部场地放大效应也可能是模拟烈度较高的原因。比较模拟地震烈度与发布烈度图的分布范围，其中更接近发布烈度图的是SW倾模型。

5结论

本文基于积石山MS6.2地震SW倾和NE倾的震源破裂模型，利用随机有限断层法模拟三分向地震动，考虑路径与场地效应，对比模拟与记录强震动台站的三分向PGA、PGV、加速度时程与PSA，给出震中附近100 km范围的仪器烈度分布，并与台站仪器烈度、正式发布的烈度图比较，结果表明：

（1） 三分向模拟PGA和PGV能够反映地震动随距离的衰减规律，PSA在0.05～10 s周期内整体模拟效果较好，地震动波形与谱形一致程度较高，可说明模拟输入震源、路径和场地参数的合理性。

（2） 模拟结果显示发震断层SW倾和NE倾模型极震区均为Ⅷ度，其中发震断层NE倾会造成比SW倾更大的地震动。从PSA残差均值、模拟烈度与发布烈度图对比结果看，发震断层SW倾更符合本次地震的实际情况。

（3） 选取适合的输入参数，利用随机有限断层法可以得到合理的三分向地震动模拟结果。本文的场地放大与衰减均依赖于vS30的取值，仅依据坡度计算的vS30值造成部分台站模拟的误差;三分向地震动模拟需要地震射线的出射角分配水平向和竖直向能量，本文仅采用平均水平地层速度结构会忽视复杂地层与地形的作用;地震应力降的不确定性也会影响模拟结果的可靠性。

致谢：感谢甘肃地震台的数据支持，感谢审稿专家为本文的修正提供宝贵意见。参考文献（References）

［1］BERESNEV I A，ATKINSON G M.Modeling finite-fault radiation from the ωn spectrum［J］.Bulletin of the Seismological Society of America，1997，87（1）：67-84.

［2］BERESNEV I A，ATKINSON G M.FINSIM-a Fortran program for simulating stochastic acceleration time histories from finite faults［J］.Seismological Research Letters，1998，69（1）：27-32.

［3］BOORE D M.Simulation of ground motion using the stochastic method［J］.Pure and Applied Geophysics，2003，160（3）：635-676.

［4］MOTAZEDIAN D，ATKINSON G M.Stochastic finite-fault modeling based on a dynamic corner frequency［J］.Bulletin of the Seismological Society of America，2005，95（3）：995-1010.

［5］傅磊，李小军.龙门山地区的kappa （k0）模型及汶川MS8.0地震的强地震动模拟［J］.地球物理学报，2017，60（8）：2935-2947.FU Lei，LI Xiaojun.The kappa （k0） model of the Longmenshan region and its application to simulation of strong ground-motion by the Wenchuan MS8.0 earthquake［J］.Chinese Journal of Geophysics，2017，60（8）：2935-2947.

［6］周红.基于NNSIM随机有限断层法的7.0级九寨沟地震强地面运动场重建［J］.地球物理学报，2018，61（5）：2111-2121.ZHOU Hong.Reconstruction of strong ground motion of Jiuzhaigou M7.0 earthquake based on NNSIM stochastic finite fault method［J］.Chinese Journal of Geophysics，2018，61（5）：2111-2121.

［7］OMAR K，ATTIA M，FERGANY E S，et al.Modeling of strong ground motion during the 1992 Cairo earthquake in the urban area northern Greater of Cairo，Egypt［J］.NRIAG Journal of Astronomy and Geophysics，2013，2（1）：166-174.

［8］GHOFRANI H，ATKINSON G M，GODA K，et al.Stochastic finite-fault simulations of the 2011 Tohoku，Japan，earthquake［J］.The Bulletin of the Seismological Society of America，2013，103（2B）：1307-1320.

［9］曹泽林.随机合成三个分量高频地震动的探讨［D］.哈尔滨：哈尔滨工业大学，2014.CAO Zelin.Preliminary study on random synthesis of three components high frequency ground motion［D］.Harbin：Harbin Institute of Technology，2014.

［10］王宏伟，任叶飞，温瑞智.一种随机有限断层的三维地震动模拟方法：鲁甸地震为例［J］.地震工程与工程振动，2021，41（2）：181-191.WANG Hongwei，REN Yefei，WEN Ruizhi.A modified stochastic finite-fault method for simulating ground motions in three dimensions：a case study of Ludian earthquake［J］.Earthquake Engineering and Engineering Dynamics，2021，41（2）：181-191.

［11］强生银，王宏伟，温瑞智，等.2021年5月21日云南漾濞MS6.4地震随机有限断层三维地震动模拟［J］.地球物理学报，2021，64（12）：4538-4547.QIANG Shengyin，WANG Hongwei，WEN Ruizhi，et al.Three-dimensional ground motion simulations by the stochastic finite-fault method for the Yangbi，Yunnan MS6.4 earthquake on May 21，2021［J］.Chinese Journal of Geophysics，2021，64（12）：4538-4547.

［12］李春果，王宏伟，温瑞智，等.2021年青海玛多MS7.4地震随机有限断层三维地震动模拟［J］.地震地质，2021，43（5）：1085-1100.LI Chunguo，WANG Hongwei，WEN Ruizhi，et al.Three-component ground motion simulations based on the stochastic finite-fault method for the 2021 Maduo MS7.4 earthquake，Qinghai Province［J］.Seismology and Geology，2021，43（5）：1085-1100.

［13］郭祥云，韩立波，张旭，等.2023年12月18日甘肃积石山6.2级地震震源参数和破裂特征［J］.地震科学进展，2024，54（1）：75-85.GUO Xiangyun，HAN Libo，ZHANG Xu，et al.Source parameters and rupture characteristics of the M6.2 Jishishan earthquake in Gansu Province on December 18，2023［J］.Progress in Earthquake Sciences，2024，54（1）：75-85.

［14］李振洪，韩炳权，刘振江，等.InSAR数据约束下2016年和2022年青海门源地震震源参数及其滑动分布［J］.武汉大学学报（信息科学版），2022，47（6）：887-897.LI Zhenhong，HAN Bingquan，LIU Zhenjiang，et al.Source parameters and slip distributions of the 2016 and 2022 Menyuan，Qinghai earthquakes constrained by InSAR observations［J］.Geomatics and Information Science of Wuhan University，2022，47（6）：887-897.

［15］杨九元，温扬茂，许才军.InSAR观测揭示的2023年甘肃积石山MS6.2地震发震构造［J/OL］.武汉大学学报（信息科学版），（2024-01-23）［2024-05-23］.https：//doi.org/10.13203/j.whugis20230501.YANG Jiuyuan，WEN Yangmao，XU Caijun.Seismogenic fault structure of the 2023 Jishishan （Gansu） MS6.2 earthquake revealed by InSAR observations［J/OL］.Geomatics and Information Science of Wuhan University，（2024-01-23）［2024-05-23］.https：//doi.org/10.13203/j.whugis20230501.

［16］杨志高，徐泰然，梁建宏等.中国地震台网中心大震震源机制CMT 产品［EB/OL］.（2024-01-01）［2024-03-03］.https：//https//data.erathquake.cn/datashare/report.shtml？PAGEID=earthquake_dzzyjzamp;cmtype=review.YANG Zhigao，XU Tairan，LIANG Jianhong et al.China Earthquake Networks Center （CENC） moment tensor （CMT） product［EB/OL］.（2024-01-01）［2024-03-03］.https：//https//data.erathquake.cn/datashare/report.shtml？PAGEID = earthquake_dzzyjzamp;cmtype = review.

［17］张培震，邓起东，张国民，等.中国大陆的强震活动与活动地块［J］.中国科学（D辑），2003，33（增刊1）：12-20.ZHANG Peizhen，DENG Qidong，ZHANG Guomin，et al.Strong earthquake activity and active blocks in Chinese mainland［J］.Scientia Sinica （Terrae），2003，33（Suppl01）：12-20.

［18］陆诗铭，吴中海，李智超.2023年12月18日甘肃省积石山6.2级地震的控震构造及特征［J］.地震科学进展，2024，54（1）：86-93.LU Shiming，WU Zhonghai，LI Zhichao.Seismic structure characteristics of the 18 December 2023 M6.2 Jishishan earthquake，Gansu Province［J］.Progress in Earthquake Sciences，2024，54（1）：86-93.

［19］左可桢，赵翠萍.2023年甘肃积石山6.2级地震序列精定位［J］.地震，2024，44（1）：204-208.ZUO Kezhen，ZHAO Cuiping.Relocation of the 2023 MS6.2 Jishishan earthquake sequence in Gansu Province［J］.Earthquake，2024，44（1）：204-208.

［20］ONISHI Y，HORIKE M.The extended stochastic simulation method for close-fault earthquake motion prediction and comments for its application to the hybrid method［J］.Journal of Structural and Construction Engineering，2004（586）：37-44.

［21］JIANG J J，BAIRD G，BLAIR D.Polarization and amplitude attributes of reflected plane and spherical waves［J］.Geophysical Journal International，1998，132（3）：577-583.

［22］ZHANG Y，WANG R J，ZSCHAU J，et al.Automatic imaging of earthquake rupture processes by iterative deconvolution and stacking of high-rate GPS and strong motion seismograms[J].Journal of Geophysical Research （Solid Earth），2014，119（7）：5633-5650.

［23］赵翠萍，陈章立，华卫，等.中国大陆主要地震活动区中小地震震源参数研究［J］.地球物理学报，2011，54（6）：1478-1489.ZHAO Cuiping，CHEN Zhangli，HUA Wei，et al.Study on source parameters of small to moderate earthquakes in the main seismic active regions，China mainland［J］.Chinese Journal of Geophysics，2011，54（6）：1478-1489.

［24］ATKINSON G M，BOORE D M.Ground-motion relations for eastern North America［J］.The Bulletin of the Seismological Society of America，1995，85（1）：17-30.

［25］WEN L X，YU S.The China seismological reference model project［J］.Earth and Planetary Physics，2023，7（5）：521-532.

［26］蒲举，尹欣欣.甘肃地区S波非弹性衰减Q值再研究［J］.四川地震，2019（1）：12-15.PU Ju，YIN Xinxin.Further study on Q-value of S-wave inelastic attenuation in Gansu area［J］.Earthquake Research in Sichuan，2019（1）：12-15.

［27］SATO H，FEHLER M，WU R S.13 scattering and attenuation of seismic waves in the lithosphere［J］.International Geophysics，2002，81：195-208.

［28］LASKE G，MASTERS G，MA Z，et al.Update on CRUST1.0-A 1-degree global model of Earth's crust［C］//Geophysical Research Abstracts.Vienna，Austria：EGU General Assembly 2013，2013，15（15）：2658.

［29］BOORE D M，THOMPSON E M.Path durations for use in the stochastic-method simulation of ground motions［J］.Bulletin of the Seismological Society of America，2014，104（5）：2541-2552.

［30］BOORE D M，THOMPSON E M.Revisions to some parameters used in stochastic-method simulations of ground motion［J］.Bulletin of the Seismological Society of America，2015，105（2A）：1029-1041.

［31］ATKINSON G M，BOORE D M.Earthquake ground-motion prediction equations for eastern North America［J］.The Bulletin of the Seismological Society of America，2007，97（3）：1032.

［32］SEYHAN E，STEWART J P.Semi-empirical nonlinear site amplification from NGA-West2 data and simulations［J］.Earthquake Spectra，2014，30（3）：1241-1256.

［33］HEATH D C，WALD D J，WORDEN C B，et al.A global hybrid vS30 map with a topographic slope-based default and regional map insets［J］.Earthquake Spectra，2020，36（3）：1570-1584.

［34］稂子平，俞瑞芳，肖亮，等.局部场地地震动高频衰减系数估计模型［J］.地震学报，2023，45（5）：919-928.LANG Ziping，YU Ruifang，XIAO Liang，et al.An estimation model of high frequency attenuation coefficient of ground motion for local site［J］.Acta Seismologica Sinica，2023，45（5）：919-928.

［35］国家市场监督管理总局，国家标准化管理委员会.中国地震烈度表：GB/T 17742—2020［S］.北京：中国标准出版社，2020.State Administration for Market Regulation，Standardization Administration of the People's Republic of China.The Chinese seismic intensity scale：GB/T 17742—2020［S］.Beijing：Standards Press of China，2020.

（本文编辑：任栋）