地球物理学报  2021, Vol. 64 Issue (10): 3619-3631   PDF    
2021年5月21日漾濞MS6.4地震近场地震动特征和方向性效应分析
张斌1, 李小军1,2, 林国良3, 荣棉水1, 俞言祥2, 王玉石1, 傅磊2     
1. 北京工业大学城市建设学部, 北京 100124;
2. 中国地震局地球物理研究所, 北京 100081;
3. 云南省地震局, 昆明 650224
摘要:利用中国国家强震动观测台网(NSMONS)和云南地震预警台网获得的158组强震动记录,分析了漾濞MS6.4地震观测的近场地震动幅值、反应谱值特征以及不同周期反应谱值的空间分布.将近场台站反应谱值与NGA-West2预测方程进行对比,结合现场震害调查结果验证此次地震近场地震动的方向性效应.研究结果表明,近场30 km内地震动峰值加速度(PGA)表现出南北向(NS)高于东西向(EW)的方向性差异,PGA最大相差达1.86倍;在水平向加速度反应谱值(SA)卓越周期附近的谱值普遍高于基本设计谱,漾濞台NS向卓越周期对应的SA是罕遇设计谱平台值的3.87倍;空间分布显示,断层破裂前方(震中东南方向)观测的水平向SA显著大于破裂后方,在短周期更为明显,随着周期的增大差异有所减小;破裂前方台站观测的短周期SA大于NGA-West2预测结果,而破裂后方台站观测的SA在整个周期范围内均小于NGA-West2预测结果,表现出显著的方向性效应;震害调查结果也表明近场区域的建筑物破坏呈现NS向破坏较EW向严重的方向性差异.
关键词: 漾濞地震      近场地震动特征      空间分布      方向性效应     
Analysis of strong ground motion characteristics and directivity effect in the near-field for the May 21, 2021 MS6.4 Yangbi earthquake
ZHANG Bin1, LI XiaoJun1,2, LIN GuoLiang3, RONG MianShui1, YU YanXiang2, WANG YuShi1, FU Lei2     
1. Faculty of Urban Construction, Beijing University of Technology, Beijing 100124, China;
2. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
3. Yunnan Earthquake Agenoy, Kunming 650224, China
Abstract: Using the 158 sets of strong ground motion recordings obtained from the National Strong Motion Observation Network System (NSMONS) of China and the Local Earthquake Early Waring Network of Yunnan Province, we analyse the characteristics of observed near-field strong ground motion amplitude, acceleration response spectra and the spatial distribution of spectral accelerations at different periods for the MS6.4 Yangbi earthquake. The spectral accelerations of the near-field stations were compared with the NGA-West2 prediction models, and the on-site earthquake damage investigation results were combined to verify the directivity effect of the near-field strong ground motion of the MS6.4 Yangbi earthquake. The research results show that the peak ground acceleration (PGA) within 30 km in the near-field shows a difference in directionality from north-south (NS) to east-west (EW), and the maximum difference in PGA is 1.86 times; The spectral accelerations near the predominant period of the horizontal acceleration response spectra (SA) are generally higher than the basic design spectra. The SA corresponding to the predominant period in the NS direction of 53YBX station is 3.87 times the value of the rare design spectra platform; The spatial distribution contours show that the horizontal SAs observed in the rupture forward (south-east of epicenter) are significantly larger than that in the rupture backward, which is more obvious in a short period, and the difference gradually decreases with the increase of the period; The short-period SAs observed in the rupture forward is greater than the median predictions of four NGA-West2 prediction models, while the SAs observed in the rupture backward is smaller than the median predictions of four NGA-West2 prediction models throughout the period, showing a significant directivity effect; The earthquake damage investigation results also show that the damage of buildings in the near-field shows a more severe directivity difference in NS direction than EW direction.
Keywords: Yangbi earthquake    Near-field strong ground motion characteristics    Spatial distribution    Directivity effect    
0 引言

据中国地震台网测定,北京时间2021年5月21日21时45分,云南省大理州漾濞县发生MS6.4地震,震中位于25.67°N,99.87°E,震源深度8 km.漾濞MS6.4地震是近十年来我国云南地区继MS6.5鲁甸地震、MS6.6景谷地震后发生的又一次6级以上破坏性浅源地震.据央视新闻报道,记者从云南省大理州委获悉,截止2021年5月22日6时,此次地震共造成30人伤亡,房屋倒塌192间,一般损坏13090间.经查询中国地震台网公布数据,截止2021年5月28日21时,共发生5级及以上余震2次,4级及以上余震11次,3级及以上余震35次.此次地震能量主要集中在前8 s内释放,矩震级为MW6.1,以向东南方向单侧走滑破裂为主,滑动极值区位于震中东南方向,主要破裂长度约为15 km(张旭等,个人通讯).云南省地震局现场工作队对灾区218个调查点进行震害调查,参照震区地震构造背景、余震分布、震源机制、仪器烈度等科技支撑成果,结合地表强震动观测记录,发布了《云南漾濞6.4级地震烈度图》(www.yndzj.gov.cn/yndzj/_300559/_300651/629959/index.html).此次地震最高烈度为Ⅷ度,Ⅵ度区及以上面积约6600 km2(含洱海),共涉及大理州6县市.

在地震过程中,地震动特征的不同往往会引起震中区建筑结构破坏表现出不同的特点,深入分析地震动观测特征有利于加深对震害的认识.本文采用漾濞地震中NSMONS和云南地震预警台网一般站获得的158组强震动记录,分析了近场台站观测记录的地震动幅值、反应谱值特征以及断层距200 km内不同周期反应谱值的空间分布;基于近场台站反应谱值与NGA-West2预测方程的对比,结合现场震害调查进一步揭示此次地震近场地震动的方向性效应特征.

1 强震动观测记录及其数据处理

NSMONS是从国家第十个“五年计划”开始建设的,包括1154个永久自由场地台站、12个专用观测阵列,于2008年3月正式运行(Li et al., 2008a).NSMONS记录了汶川MS8.0大地震和芦山MS7.0地震的大量强震动记录,极大地丰富了中国强震动数据库,尤其是增加了大量的大震近场记录(Li et al., 2008b; Xie et al., 2014).国家地震烈度速报与预警工程于2015年8月正式启动,计划建成1800个配置宽频带地震仪和加速度仪的基准站、3600个配置加速度仪的基本站、超过10000个配置烈度计的台站,将为重灾区位置和灾情分布的快速判定,社会公众的避震逃生,重大工程的制动、关闭的紧急处置等提供服务,减轻直接地震灾害及次生灾害损失(蒋长胜和刘瑞丰,2016).漾濞地震发生后,中国地震局工程力学研究所强震动观测组公布了NSMONS获得的29组自由场地的强震动记录,云南地震预警台网一般站获得了146组自由场地记录.本文以这两个台网获得的158组断层距200 km内数据为基础开展研究.经查询,本文选用的NSMONS的12组强震动记录均由中国地震局工程力学研究所生产的SLJ-100型3分向力平衡加速度计记录,选用的云南地震预警台网一般站的146组加速度记录均由北京港震科技股份有限公司生产的GL-P2C烈度计记录,GL-P2C烈度计采用的是高灵敏度3分向微电机Mems加速度计,两种加速度计的主要技术指标如表 1所示.由表 1可知,除了采样率有所差异,两种加速度计的主要技术指标基本相同,具有较好的一致性.将烈度计的记录增采至200 SPS后与强震动数据同时应用于本研究.地震的震中及其触发的断层距200 km内台站分布如图 1所示.

表 1 加速度计主要技术指标 Table 1 Main technical indicators for two types of accelerometer
图 1 漾濞MS6.4地震震中及触发的强震动台站和云南地震预警台网一般站分布,黑色实线圈为云南省地震局发布的宏观烈度分布 Fig. 1 Locations of epicenter and recording stations of the MS6.4 Yangbi earthquake. Triangles indicate the strong motion stations in the NSMONS; Circles indicate the general stations in the Local Earthquake Early Waring Network of Yunnan Province; Star indicates the epicenter; Red box indicates the fault plane; Black solid coils are the macroscopic intensity distribution issued by Yunnan Earthquake Agenoy

本研究中的近场强震动记录采用近断层强震动记录的基线校正改进方法(张斌等,2020)进行处理,避免了滤波对位移时程和永久位移信息的影响.剩余记录均遵循PEER NGA数据处理流程(Ancheta et al., 2014),应用4阶Butterworth非因果带通滤波进行处理,低通截止频率为30 Hz.为了保证确定的滤波高通截止频率的有效性和尽量减小滤波对峰值位移(PGD)的影响,首先根据震源谱模型的理论频率、加速度的傅里叶振幅谱在低频段的衰减正比于f2,噪声的出现会造成傅里叶振幅谱在低频段的翘起(Joyner and Boore, 1988Atkinson and Silva, 2000)以及记录的信噪比不小于3(Boore and Bommer, 2005)等判定标准确定高通截止频率的范围.然后在该范围内以0.01 Hz的间隔逐一选取高通截止频率进行处理,最后根据处理后的速度和位移时程的末尾是否在零线确定最终的高通滤波截止频率;再进行后处理(Boore et al., 2012).因此,处理后的地震动获得的PGA、峰值速度(PGV)、PGD和SA是可靠的.

2 漾濞地震观测的强震动特征 2.1 幅值特征

漾濞地震有158组强震动记录位于断层距200 km以内,其中30 km内的10组.断层距(台站到断层破裂面最近的距离)、断层投影距(台站到断层破裂面在地表投影的水平距离)的计算采用了张旭等(个人通讯)有限断层模型,由于其初始破裂点与中国地震台网公布的震中位置有所差别,本文中采用的是中国地震台网公布的震中位置,并将其有限断层模型投影到该震中位置,如图 1中断层破裂面所示.漾濞地震的有限断层模型参数如表 2所示.记录到最大PGA的台站是断层距4.67 km的53YBX台,其EW、NS和UD向的PGA分别为374.40 Gal、-695.65 Gal和-407.42 Gal,三个方向的PGV分别为-27.18 cm·s-1、-27.27 cm·s-1和-7.90 cm·s-1.表 3给出了近场30 km内10个台站记录及三个方向的相关地震动参数.图 2绘制了53YBX台、L2203台和L2205台处理后的EW、NS和UD方向的加速度和速度时程.由图 2表 3可知,近场30 km内大部分台站NS向的PGA高于EW向,最大相差1.86倍,呈现显著的方向性差异.

表 2 漾濞MS6.4地震破裂尺度 Table 2 The rupture dimensions of the MS6.4 Yangbi earthquake
表 3 漾濞MS6.4地震近场30 km内台站记录及其相关参数 Table 3 Ground motion parameters of strong-motion recording of stations at distance less than 30 km observed in the MS6.4 Yangbi earthquake
图 2 漾濞MS6.4地震近场台站(a)53YBX、(b)L2203、(c)L2205记录的EW、NS和UD向加速度和速度时程 Fig. 2 Acceleration and velocity time histories in the EW, NS, and UD directions of (a) 53YBX, (b) L2203, (c) L2205 stations
2.2 反应谱特征

经查询《建筑抗震设计规范》(GB50011-2010),53YBX台、L2203台和L2205台所在地区抗震设防烈度均为Ⅷ度,设计地震动分组为第三组(中华人民共和国住房和城乡建设部和中华人民共和国国家质量监督检验检疫总局,2010).根据《中国地震动参数区划图》(GB18306-2015)规定,53YBX台所在地区Ⅱ类场地基本地震动PGA和SA特征周期分别为0.2g和0.45 s,Ⅱ类场地罕遇地震动PGA和SA特征周期分别为0.38g和0.5 s;L2203台和L2205台的分别为0.15g和0.45 s,0.285g和0.5 s(中华人民共和国国家质量监督检验检疫总局和中国国家标准化管理委员会,2015).分别计算出三个台站强震动记录的5%阻尼比EW和NS向加速度反应谱值(SA(T=0.01~6.0 s)),并将其与该地区的设计反应谱进行对比,如图 3所示.

图 3 近场台站EW和NS向5%阻尼比加速度反应谱与设计反应谱的对比 Fig. 3 Comparison of horizontal 5% damping acceleration response spectra from near-field stations with design spectra

53YBX台EW、NS向加速度反应谱在0.14 s前明显高于罕遇地震动设计谱,分别在0.14~1.60 s、0.14~1.10 s高于基本地震动设计谱而低于罕遇地震动设计谱,在其他周期范围均低于基本地震动设计谱.反应谱卓越周期在0.1 s,0.1 s后反应谱值迅速下降,EW、NS向卓越周期对应的反应谱值分别是罕遇地震动设计谱平台值的1.54倍和3.87倍.L2203台EW、NS向反应谱值分别在0.12 s和0.16~0.32 s、0.22~0.28 s高于罕遇地震动设计谱,卓越周期分别为0.24 s和0.28 s,卓越周期对应的反应谱值分别是罕遇地震动设计谱平台值的1.21倍和1.59倍.L2205台NS向反应谱值在0.18~0.32 s高于基本地震动设计谱而低于罕遇地震动设计谱,NS向在其他周期范围和EW向在几乎所有周期范围内均低于基本地震动设计谱;NS向反应谱卓越周期在0.24 s,卓越周期对应的反应谱值是基本地震动设计谱平台值的1.66倍.近场台站的加速度反应谱值在卓越周期附近呈现NS向显著大于EW向的方向性差异.

2.3 观测的水平向地震动的空间分布

上述研究表明近场地震动幅值和反应谱值表现出EW向和NS向的方向性差异,为了更直观地观察两个方向地震动强度的差异,图 4给出了EW和NS向PGA、SA(T=0.1 s、0.5 s、1.0 s、3.0 s、5.0 s)的空间分布对比.从图 4可知,极震区NS向的PGA、SA(T=0.1 s)显著大于EW向,这与前两节中幅值特征和反应谱特征相符;极震区NS向的SA(T=0.5 s、1.0 s)显著小于EW向,而NS向的SA(T=3.0 s、5.0 s)略大于EW向.此次地震以向东南方向单侧走滑破裂为主(张旭等,个人通讯),无论是EW向还是NS向,不同周期的地震动在断层破裂的前方(震中东南方向)和后方(震中北西方向)表现出明显的差异.在相同的距离处,断层破裂前方观测的加速度反应谱值明显大于破裂后方.在短周期(PGA~0.1 s)时,随着周期的增大,断层破裂前方和破裂后方的观测加速度反应谱值的差异越来越大,破裂的方向性更为突出;当周期大于0.1 s时,随着周期的增大,断层破裂前方与破裂后方的观测加速度反应谱值的差异逐渐减小,破裂的方向性没有像短周期时那么突出.为了更好地研究断层破裂前方和后方的水平向地震动强度的方向性差异,图 5给出了EW向和NS向地震动PGA和SA(T=0.1 s、0.5 s、1.0 s、3.0 s、5.0 s)的几何均值的空间分布.图 5也显示出与图 4相同的断层破裂前方观测的加速度反应谱值明显大于破裂后方的方向性差异,在短周期(PGA~0.1 s)时,随着周期的增大,断层破裂前方和破裂后方的观测加速度反应谱值的差异越来越大,破裂的方向性更为突出;当周期大于0.1 s时,随着周期的增大,断层破裂前方与破裂后方的观测加速度反应谱值的差异逐渐减小,破裂的方向性没有像短周期时那么突出.

图 4 EW和NS向PGA、SA(T=0.1 s、0.5 s、1.0 s、3.0 s、5.0 s)的空间分布对比 Fig. 4 Comparison of the spatial distribution of PGA and SA (T=0.1 s, 0.5 s, 1.0 s, 3.0 s, 5.0 s) between EW and NS direction. Star indicates the Yangbi MS6.4 earthquake epicenter; Red box indicates the fault plane
图 5 EW和NS向PGA、SA(T=0.1 s、0.5 s、1.0 s、3.0 s、5.0 s)几何均值的空间分布 Fig. 5 Contours of geometric mean of PGA and SA (T=0.1 s, 0.5 s, 1.0 s, 3.0 s, 5.0 s) in the direction of EW and NS. Star indicates the Yangbi MS6.4 earthquake epicenter; Red box indicates the fault plane

震中东南方向40多公里处的长周期地震动(SA(T=3.0 s、5.0 s))明显大于震中及附近区域的长周期加速度反应谱值,表现出明显的长周期放大特征.经查询,该区域预警台网一般站布设在洱海周边,洱海南、西、北三面受河流冲积形成洪积平原,十八条溪溪口横向连接成洪积裙,深厚的沉积层放大了漾濞地震的长周期地震动.震后地震科考发现位于洱海边上的云南信创大理产业园大楼(25.619°N,100.287°E)的22~26层的楼梯间和电梯间之间的填充墙上出现剪切裂缝,且随着楼层增加,剪切裂缝越明显,初步推测与深厚沉积层对长周期地震动放大以及建筑物结构的地震动放大有关.

3 近场地震动的方向性效应

上述近场台站获取的地震动幅值、反应谱值表现出明显的方向性差异,不同周期的水平向地震动加速度反应谱值在断层破裂前方(震中东南方向)和后方(震中北西方向)也存在显著差异.初步认为是受此次走滑地震震源破裂的方向性效应影响.为了进一步验证这一推测,将断层破裂前方、后方近场台站记录的水平向加速度反应谱值与NGA-West2预测方程进行对比分析.

NGA-West2预测方程是在NGA-West1的基础上,利用NGA-West2数据库建立的考虑了震源特性(震级、断层类型、上/下盘效应)、距离(几何衰减、非弹性衰减)、场地(土层线性/非线性效应、盆地响应)和区域的差异性等影响因素的水平向地震动预测模型(Abrahamson et al., 2014Boore et al., 2014Campbell and Bozorgnia, 2014Chiou and Youngs, 2014Idriss, 2014).这些考虑多个影响因素的地震动预测方程很大程度上减小了地震动估计中的不确定性.强震动台站的场地VS30来源于谢俊举等(个人通讯)收集的四川、云南地区强震动台站的钻孔波速资料给出的值;由于云南地震预警台网一般站布设时没有钻孔,根据Allen和Wald(2009)提出的地形坡度与VS30的修正关系获得场地VS30,强震动台站和云南预警台网一般站的场地VS30列于表 3和附录中.图 6给出了断层破裂前方台站(53YBX、L2205、L0101)和后方台站(L2201、L2204、L3003)记录的水平向加速度反应谱值与四个NGA-West2预测方程的对比,张旭等(个人通讯)有限断层模型参数列于表 2中,6个台站的断层距和场地VS30列于表 3中.从图 6a可以看出,断层破裂前方台站观测的水平向加速度反应谱值在短周期处高于四个NGA-West2预测方程的结果,而图 6b显示断层破裂后方台站观测的水平向加速度反应谱值在整个周期范围内均低于四个NGA-West2预测方程的结果.由水平向地震动加速度反应谱值的空间分布以及与NGA-West2预测方程的对比结果可知,此次漾濞地震的地震动表现出显著的方向性差异,断层破裂前方观测的地震动高频成分明显高于预测结果,断层破裂后方的则总体上明显低于预测结果.

图 6 断层破裂前方(a)和破裂后方台站(b)记录的水平向加速度反应谱值与4个NGA-West2预测方程的对比 Fig. 6 Comparison of observed horizontal spectral accelerations in the forward direction (a) and the backward direction (b) with the median predictions from four NGA-West2 ground motion prediction equations
4 震害的方向性特征

根据现场震害调查结果,震中区破坏最严重的是1层老旧穿斗木构架结构,部分土墙倒塌、木梁未破坏;而2~3层砖混、框架结构破坏很轻或几乎没有破坏.53YBX台附近的下关完小(25.665°N,99.965°E)一栋EW走向的三层砖混结构宿舍楼(坐北朝南),其墙体多处脱落,如图 7(ac),而另一栋NS走向的两层砖混结构教学楼(坐西朝东)基本完好无损,如图 7b.宏观震中淮安村的房屋震害也明显体现了震害的方向性,如图 8中两个紧邻的房屋(25.691°N,99.942°E),其中坐北朝南的房屋1梭瓦现象严重,墙皮大面积脱落,土坯墙体呈贯通竖向裂缝,而东西朝向的房屋2仅见小面积墙皮脱落和填充墙开裂,梭瓦不明显,其震害远轻于房屋1.相似的震害现象在微观震中秀岭村和距离震中8 km左右的漾濞县城苍山西镇(震中东南方位)多处可见.地震近场区的建筑物破坏呈现NS向破坏较EW向严重的方向性差异,很好地验证了前述近场地震动特征以及断层破裂前方和后方台站观测到地震动强度的方向性差异.

图 7 下关完小EW向宿舍楼(a、c)和NS向教学楼(b)震害情况(经纬度25.665°N,99.965°E,照片a镜向正北,照片b镜向正西) Fig. 7 The earthquake damage of EW to the dormitory building (a, c) and NS to the teaching building (b) of Xiaguan Primary School
图 8 宏观震中淮安村调查点震害情况(经纬度25.691°N,99.942°E,照片1镜向正西,照片2镜向北西) Fig. 8 The earthquake damage at the survey point of Huai′an village in the macro-epicenter
5 结论与讨论

本文采用漾濞地震中NSMONS和云南地震预警台网一般站获得的158组强震动记录,分析了近场台站观测记录的地震动幅值、反应谱值特征以及断层距200 km内不同周期反应谱值的空间分布;基于近场台站记录的反应谱值与NGA-West2预测方程的对比研究,并结合现场震害调查结果进一步揭示此次地震近场地震动的特点.得到了如下主要结论:(1)近场30 km内大部分台站记录NS向的PGA高于EW向的方向性差异,最大可达1.86倍;(2)近场台站记录的EW、NS向加速度反应谱在卓越周期附近高于基本地震动设计谱,在NS向更为显著,53YBX台NS向卓越周期对应的反应谱值是罕遇地震动设计谱平台值的3.87倍;(3)EW向、NS向的加速度反应谱值和两者几何均值的空间分布图表现出在相同的距离处,断层破裂前方观测的不同周期处加速度反应谱值均明显大于破裂后方的方向性差异.在短周期(PGA~0.1 s)时,随着周期的增大,两者的差异越来越大,破裂的方向性更为突出;当周期大于0.1 s时,随着周期的增大,两者的差异逐渐减小,破裂的方向性没有像短周期时那么突出;(4)断层破裂前方台站观测的水平向加速度反应谱值在短周期处高于NGA-West2预测方程的结果,而断层破裂后方台站观测的水平向加速度反应谱值在整个周期范围内均低于NGA-West2预测方程的结果,表现出显著的方向性差异;(5)近场区域的1~3层建筑结构破坏也呈现NS向破坏较EW向严重的方向性差异.很好地佐证了此次走滑地震近场地震动的方向性效应.

此次地震引起的极震区地震动的空间分布,在PGA、0.1 s周期处NS方向地震动明显高于EW方向地震动,在周期0.5 s、1.0 s处EW方向地震动反而高于NS方向地震动,在长周期3.0 s、5.0 s处NS方向地震动略高于EW方向但二者差别不大.这一现象说明由于震源、传播路径、场地的复杂性,不同周期处地震动的空间分布规律并不是在所有频段都一致,而是呈现出分频段不一致的特征.由于地震动等值线的高点以及中心主要位于漾濞县城,从场地条件来看整个县城几乎都分布在南北长约5 km,东西长约2 km的山间盆地内,EW与NS方向沉积土层厚度和水平方向分布(即盆地效应不同方向)的差异对周期0.5 s和1.0 s处地震动放大程度有差别,这可能是造成其它周期处NS向地震动高于EW向,而0.5 s和1.0 s处EW方向地震动反而高于NS方向地震动的主要原因.

近场地震动的方向性效应主要是由于震源破裂的方向性引起的,往往会放大断层破裂前方的地震动(谢俊举等,2017安昭等,2019).张旭等(个人通讯)利用远场体波数据反演得到的漾濞地震震源破裂过程显示此次地震以向东南方向单侧走滑破裂为主,滑动极值区位于震中东南方向.张克亮等(2021)通过弹性半空间断裂位错模型结合GNSS同震变形场约束反演得到的此次地震最大滑动量也位于主震震源以南区域.这与本文得到的断层破裂前方(震中东南方向)观测的水平向加速度反应谱值明显大于破裂后方(震中北西方向)的方向性效应具有很好的一致性.

云南省地震局发布的《云南漾濞6.4级地震烈度图》显示,等震线长轴呈北北西走向,此次地震最高烈度为Ⅷ度,主要集中在漾濞县苍山西镇、漾江镇、太平乡3个乡镇.此次地震科考发现下关完小、宏观震中淮安村和微观震中秀岭村的建筑物破坏呈现NS向破坏较EW向严重的方向性差异,这些地区均位于断层破裂前方的苍山西镇.断层破裂前方台站53YBX、L2205、L0101、L2203也位于苍山西镇和太平乡.本文显示的近场台站幅值和反应谱值呈NS向显著大于EW向的方向性差异、断层破裂前方(震中东南方向)观测的水平向SA显著大于破裂后方(震中北西向)、近场断层破裂前方台站观测的水平向加速度反应谱值在短周期处高于NGA-West2预测方程结果以及近场区域建筑结构破坏的方向性差异等特征与该宏观烈度分布图等震线长轴方向和最高烈度分布区基本一致.近场地震动的方向性效应往往会放大断层破裂前方的地震动,造成建筑结构更为严重的破坏.因此,在建筑结构的抗震设计和地震动衰减关系的研究中需考虑近场地震动方向性效应的影响.

附录
附表 1 强震动台站和云南地震预警台网一般站的场地VS30 Appendix Table 1 Site VS30 of the strong motion stations in the NSMONS and the general stations in the Local Earthquake Early Waring Network of Yunnan Province
致谢  感谢中国地震局工程力学研究所为本研究提供数据支持.审稿专家对本文的完善提出了非常宝贵的意见,在此一并表示感谢.
References
Abrahamson N A, Silva W J, Kamai R. 2014. Summary of the ASK14 ground motion relation for active crustal regions. Earthquake Spectra, 30(3): 1025-1055. DOI:10.1193/070913EQS198M
Allen T I, Wald D J. 2009. On the use of high-resolution topographic data as a proxy for seismic site conditions (VS30). Bulletin of the Seismological Society of America, 99(2A): 935-943. DOI:10.1785/0120080255
An Z, Xie J J, Li X J, et al. 2019. Directivity effects on strong ground motion from the February 6, 2018 MW6.4 Hualian earthquake. Chinese Journal of Geophysics (in Chinese), 62(12): 4658-4672. DOI:10.6038/cjg2019M0579
Ancheta T D, Darragh R B, Stewart J P, et al. 2014. NGA-West2 database. Earthquake Spectra, 30(3): 989-1005. DOI:10.1193/070913EQS197M
Atkinson G M, Silva W J. 2000. Stochastic modeling of California ground motions. Bulletin of the Seismological Society of America, 90(2): 255-274. DOI:10.1785/0119990064
Boore D M, Bommer J J. 2005. Processing of strong-motion accelerograms: needs, options and consequences. Soil Dynamics and Earthquake Engineering, 25(2): 93-115. DOI:10.1016/j.soildyn.2004.10.007
Boore D M, Sisi A A, Akkar S. 2012. Using pad-stripped acausally filtered Strong-Motion data. Bulletin of the Seismological Society of America, 102(2): 751-760. DOI:10.1785/0120110222
Boore D M, Stewart J P, Seyhan E, et al. 2014. NGA-West2 equations for predicting PGA, PGV, and 5% damped PSA for shallow crustal earthquakes. Earthquake Spectra, 30(3): 1057-1085. DOI:10.1193/070113EQS184M
Campbell K W, Bozorgnia Y. 2014. NGA-West2 ground motion model for the average horizontal components of PGA, PGV, and 5% damped linear acceleration response spectra. Earthquake Spectra, 30(3): 1087-1115. DOI:10.1193/062913EQS175M
Chiou B S J, Youngs R R. 2014. Update of the Chiou and Youngs NGA model for the average horizontal component of peak ground motion and response spectra. Earthquake Spectra, 30(3): 1117-1153. DOI:10.1193/072813EQS219M
General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, China National Standardization Administration Committee. 2015. GB 18306-2015 Seismic ground motion parameters zonation map of China (in Chinese). Beijing: China Standards Press.
Idriss I M. 2014. An NGA-West2 empirical model for estimating the horizontal spectral values generated by shallow crustal earthquakes. Earthquake Spectra, 30(3): 1155-1177. DOI:10.1193/070613EQS195M
Jiang C S, Liu R F. 2016. National seismic intensity rapid reporting and early warning project: opportunity and challenge of seismic network. Journal of Engineering Studies (in Chinese), 8(3): 250-257.
Joyner W B, Boore D M. 1988. Measurement, characterization, and prediction of strong ground motion. //Earthquake Engineering and Soil Dynamics Ⅱ-Recent Advances in Ground-Motion Evaluation: Proceedings of the Specialty Conference, Park City, UT, USA: Geotechnical Special Publication, 43-102.
Li X J, Zhou Z H, Huang M, et al. 2008a. Preliminary analysis of strong-motion recordings from the magnitude 8.0 Wenchuan, China, earthquake of 12 May 2008. Seismological Research Letters, 79(6): 844-854. DOI:10.1785/gssrl.79.6.844
Li X J, Zhou Z H, Yu H Y, et al. 2008b. Strong motion observations and recordings from the great Wenchuan Earthquake. Earthquake Engineering and Engineering Vibration, 7(3): 235-246. DOI:10.1007/s11803-008-0892-x
Ministry of Housing and Urban Rural Development of the People's Republic of China, General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China. 2010. GB 50011-2010 Code for seismic design of buildings (in Chinese). Beijing: China Architecture & Building Press.
Xie J J, Li X J, Wen Z P, et al. 2014. Near-source vertical and horizontal strong ground motion from the 20 April 2013 MW6.8 Lushan earthquake in China. Seismological Research Letters, 85(1): 23-33. DOI:10.1785/0220130121
Xie J J, Li X J, Wen Z P, et al. 2017. Long period strong ground motion from the April 15, 2016 Kumamoto MW7.0 earthquakein Japan. Chinese Journal of Geophysics (in Chinese), 60(11): 4431-4446. DOI:10.6038/cjg20171129
Zhang B, Yu Y X, Xiao L. 2020. An improved method for near-fault strong ground motion records' baseline correction. Journal of Vibration and Shock (in Chinese), 39(5): 137-142, 163.
Zhang K L, Gan W J, Liang S M, et al. 2021. Coseismic displacement and slip distribution of the 2021 May 21, MS6.4, Yangbi earthquake derived from GNSS observation. Chinese Journal of Geophysics (in Chinese), 67(7): 2253-2266. DOI:10.6038/cjg2021O0524
安昭, 谢俊举, 李小军, 等. 2019. 2018年2月6日花莲MW6.4地震近场地震动方向性效应. 地球物理学报, 62(12): 4658-4672. DOI:10.6038/cjg2019M0579
蒋长胜, 刘瑞丰. 2016. 国家地震烈度速报与预警工程——测震台网的机遇与挑战. 工程研究-跨学科视野中的工程, 8(3): 250-257.
谢俊举, 李小军, 温增平. 2017. 日本熊本MW7.0地震的长周期地震动. 地球物理学报, 60(11): 4431-4446. DOI:10.6038/cjg20171129
张斌, 俞言祥, 肖亮. 2020. 近断层强震记录基线校正的改进方法. 振动与冲击, 39(5): 137-142, 163.
张克亮, 甘卫军, 梁诗明, 等. 2021. 2021年5月21日MS6.4漾濞地震GNSS同震变形场及其约束反演的破裂滑动分布. 地球物理学报, 67(7): 2253-2266. DOI:10.6038/cjg2021O0524
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会. 2015. GB 18306-2015中国地震动参数区划图. 北京: 中国标准出版社.
中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局. 2010. GB 50011-2010建筑抗震设计规范. 北京: 中国建筑工业出版社.