-
2020年10月22日11时03分37秒四川省绵阳市北川县发生了MS4.7地震，震中位于北川县坝底乡（东经104.18°，北纬31.83°），距离北川县城36 km，距绵阳市区62 km，距成都市区132 km，震中20 km内有坝底乡、马槽乡、禹里镇、开坪乡、白什羌族藏族乡、白什乡、桃龙藏族乡等乡镇. 绵阳地区震感明显，广元、德阳、成都、资阳、雅安等地有震感. 地震发生后，四川省地震局对震中坝底乡及周边地区展开现场烈度调查，震中附近30 km范围16个烈度调查点，仅有个别房屋出现裂缝、缩瓦等现象，无房屋倒塌，最大烈度为Ⅵ度.
北川县MS4.7地震发生在岷江断裂、虎牙断裂和龙门山断裂带的交汇区域，位于青藏高原东缘，属巴颜喀拉块体与扬子地块交界部位，历史上曾发生过多次中强震. 岷江断裂发生过1713年叠溪7级地震、1748年漳腊北6½级地震、1938年松潘6级地震、1960年漳腊6¾级地震和1933年叠溪7½级地震，虎牙断裂于1976年发生过两次7.2级地震，龙门山断裂带于2008年发生汶川8.0级地震. 根据四川地震台给出的震源机制解北川MS4.7地震为逆冲型地震，震源深度17 km. 2009年以来此区域发生了14次4级以上地震，北川MS4.7地震被认为是发生在龙门山断裂带汶川余震区的一组地震，为汶川余震区一次正常的起伏活动. 近年来，中国地震局开展了地震烈度速报台站等建设项目，地震活动性较高的区域强震动台网的密度得到了提升，如2020年7月2日唐山古冶MS5.1地震中共收集到363组三分量强震动数据，这些强震动记录对研究中小地震地震动的特征具有重要意义（温瑞智和王宏伟，2020）. 本文简要介绍了北川MS4.7地震的强震动记录，并初步进行了区域地震动衰减特征和场地特性等方面的分析.
-
“十五”期间，四川建成了由224个数字强震台组成的数字强震动观测台网，其中自由场固定数字强震动观测台站216个、地形影响台阵1个。观测台站分布于鲜水河—安宁河—则木河断裂带以及龙门山断裂带周边，台站以自由场地表观测为主，安装内置力平衡式加速度计的强震仪，采用阈值触发模式存储强震动记录，在汶川地震和芦山地震中获取到大量近场强震动记录，为我国开展工程地震学研究提供了重要基础. 地震烈度速报与预警工程在四川地区建成了1013个简易强震动台站（一般站、烈度速报台站），采用MEMS芯片传感器的烈度仪，选择乡镇政府内自由地表、通信铁塔站点的机房作为观测场地，部分台站烈度仪固定在不高于地面30 cm的墙体上，观测数据通过网络实时传输至数据中心.
北川MS4.7地震震中附近台站分布密集，震中距200 km内有369个台站. 截止2020年10月23日，北川MS4.7地震有547个台站触发，扣除个别台站垂直分量未获取有效记录和远距离台站信噪比过低的情况，共获取三分向强震动记录532组. 图1给出了北川MS4.7地震中有效获取记录的强震动台站分布情况，图2给出了强震动记录数量与震中距关系的统计情况，其中震中距50 km内获取三分向记录25组，200 km内获取三分向记录306组. 10月21日此区域还发生了一次MS4.6地震，获取到489组三分向强震动记录. 这两次地震观测到丰富的近场强震动观测记录，明显改善了此区域地震动记录匮乏的现状.
图 1 北川MS4.7地震中获取记录的强震动台站分布图
Figure 1. Strong-motion stations recorded the Beichuan MS4.7 earthquake
图 2 强震动记录数量与震中距统计图
Figure 2. Statistical chart of the number of strong motion records and epicenter distance
-
开展地震动特征分析前需对记录数据进行必要的处理，本文利用均值法对获取的强震动记录进行了基线调整（于海英等，2009），在仿真速度和频谱分析之前进行了巴特沃思带通滤波处理（金星等，2004），滤波范围0.1～20 Hz. 为更好地解读北川MS4.7地震的特点，本文选取了震中距50 km范围内的记录计算了相关参数，见表1. 本次地震的最大峰值地面加速度（PGA）记录由距离较近的51BCB台站获得，震中距仅2.9 km，东西、南北和垂直向PGA分别为−277.53 cm/s2、−250.72 cm/s2、103.51 cm/s2，峰值速度（PGV）分别为3.54 cm/s、4.12 cm/s、0.98 cm/s，计算仪器烈度7.1度. 该台址场地所处地貌单元为山区坡脚沟谷部位，地形起伏不大，但地层较复杂，地勘钻孔所揭露的地层主要为第四系全新统残坡积层和志留系中上统茂县群千枚岩，覆盖层厚度10.2 m，其中松散土以第四系全新统杂填土、第四系残坡积含碎石粉质粘土等构成，覆盖层等效剪切波速为215 m/s，属Ⅱ类场地. 51BCQ台站为距离震中最近的基岩场地台站，震中距26.7 km. 图3给出了51BCB和51BCQ记录的三分向加速度时程和5%阻尼比拟加速度反应谱（PSA）. 相比51BCQ的基岩强震动记录，51BCB的土层记录水平向地震动明显强于垂直向，记录相对持时较短，有明显尖锐脉冲现象，与近场地震动高频成分丰富的特点相一致，同时表明覆盖层对水平向地震动的放大作用. 根据建筑抗震设计规范（GB50011—2010），51BCB场址所在区域为第二组，设计地震基本加速度0.2g，特征周期值取0.40 s，我国中小城市和城镇主要建筑物的结构自振周期约为0.3～1 s（任叶飞等，2014），此范围内51BCB的PSA明显小于设计谱值，在实际震害调查中发现，台站周边房屋仅有少量裂缝，未造成严重破坏，也与反应谱特征基本吻合.
表 1 地震动记录及相关参数表
Table 1. Seismic records and related parameters
台站代码 纬度/(°) 经度/(°) 震中距/km 峰值加速度/（cm·s−2） 峰值速度/（cm·s−1） 仪器烈度 EW NS UD EW NS UD 51BCB 31.8 104.2 2.9 −277.53 −250.72 103.51 3.54 4.12 0.98 7.1 B2603 31.9 104.1 10.9 158.11 180.88 142.52 −3.85 −4.28 1.66 6.7 B0571 31.7 104.3 12.5 156.39 152.27 84.13 2.81 −2.27 −1.76 6.2 U2302 31.8 104.0 18.8 −84.89 58.35 −113.25 −1.05 0.46 −1.51 5.7 B2607 32.0 104.2 21.5 30.62 23.74 21.00 1.11 1.11 0.68 4.7 B0502 31.7 104.4 24.8 23.86 10.06 19.62 −1.62 0.49 0.67 4.5 B2604 32.0 104.0 26.5 40.26 36.04 −35.32 −0.87 −0.95 0.72 4.8 51BCQ 31.8 104.5 26.7 11.81 7.34 6.59 0.36 0.22 0.19 3.1 B2608 31.9 104.5 27.3 71.05 46.59 20.26 −1.73 0.89 −0.42 5.5 51BCY 31.7 104.4 29.1 64.28 37.23 −18.70 1.07 −0.69 −0.27 5.0 U2310 31.8 103.9 30.3 39.08 −52.45 −30.54 0.51 −1.07 0.27 4.8 F8303 31.6 104.1 30.3 −53.09 57.13 27.28 −1.22 1.46 −0.62 5.1 B0501 31.6 104.3 31.4 44.66 −17.62 43.60 0.50 −0.40 −0.51 4.4 B2606 31.8 104.6 36.4 −59.36 −35.03 16.99 −2.36 −0.96 −0.37 5.4 B0503 31.6 104.4 38.1 −25.13 −9.91 −37.91 −0.79 0.54 −1.32 4.9 B2601 31.9 104.6 39.4 28.32 −23.96 9.10 −2.01 0.76 0.36 4.8 B2602 32.0 104.5 39.8 15.67 13.63 −12.76 −0.14 −0.16 −0.17 3.0 F8301 31.5 104.2 41.1 −15.99 −15.52 −6.07 0.53 −0.55 0.29 3.7 U2301 31.8 103.7 42.1 21.87 30.39 −29.72 −0.19 0.40 −0.10 3.9 B8112 31.7 104.6 42.8 −19.51 −13.82 −5.24 0.95 0.59 −0.30 4.1 B2713 32.2 104.5 46.3 16.53 17.99 8.34 −0.23 −0.32 0.14 3.4 U2306 31.9 103.7 46.6 24.66 34.78 35.83 −0.30 −0.37 0.17 3.9 B2703 32.1 104.5 47.6 46.73 48.69 −28.52 −0.19 −0.42 0.13 4.3 51MZX 31.4 104.2 48.0 6.60 −4.60 −3.28 −0.45 −0.30 −0.15 3.0 B8106 31.7 104.7 49.5 −6.32 11.61 −2.67 0.48 0.28 −0.13 3.2 
图 3 51BCB台站和51BCQ台站加速度时程与5% 阻尼比拟加速度反应谱
Figure 3. Acceleration time histories and 5%-damped pseudospectral accelerations for the recorded at stations 51BCB and 51BCQ
PGA和PGV的空间分布特性是反映地震动特征的重要参数，能够为地震烈度评估和判断能量传播优势方向提供直观认识，见图4. 本文计算得到了532组记录的PGA和PGV参数，采用克里金方法对台站空区进行插值，图4给出了东西、南北和垂直3个方向的PGA和PGV等值线，为了便于更清晰显示等值线，将PGA和PGV取对数后再作插值. 从图中可以看出，地震动较强的记录主要集中在震中附近区域，PGA和PGV等值线长轴呈北西—南东方向，其中受场地影响较弱的垂直分量长轴方向展布更加明显，体现出地震动场在空间上的分布特征. 在震中外围区域PGA和PGV等值线均出现零散的异常增大区域，考虑是因布设简易强震动台的通信站点基本位于山坡和山顶位置，受到局部地形和场地的影响导致.
图 4 PGA和PGV空间分布图
Figure 4. Contours for PGA and PGV
图 5 PGA和PGV随震中距衰减分布
Figure 5. PGAs and PGVs attenuated with epicentral distance
地震波在传播过程中受震源和传播介质的影响，地震动随距离的增加可表现出不同的衰减特性（王文才等，2017）. 虽然近年来我国强震观测数据得到了快速积累，但全国大部分地区仍不足以用强震数据直接回归得到衰减关系（俞言祥等，2013），主要采用转换方法（胡聿贤和张敏政，1984）选择美国西部、日本和我国台湾地区为参考，建立该地区的地震动衰减关系，并利用我国强震观测数据加以控制，通过转换得到我国各分区地震动参数衰减关系.
本文统计了北川MS4.7地震水平分量 PGA、PGV 随震中距R的衰减特征，并与常用衰减关系预测值进行初步对比. 考虑震中位置和震级区间因素，选用第五代区划图中M≤6.5的青藏区长轴方向衰减关系作为对比预测值之一，另外选用了俞言祥和汪素云（2006）提出的中国西部地区水平向长轴基岩加速度衰减关系和雷建成等（2007）提出的西南地区衰减关系进行对比. 由于建立地震动参数衰减关系过程中所使用的强震记录多分布于中等震级和中等距离，俞言祥和汪素云（2006）以及雷建成等（2007）两种衰减关系主要适用于中强震级、震中距200 km以内的范围.
图5可看出此次地震PGA实测值随距离的衰减趋势与俞言祥和汪素云（2006）衰减关系的预测值一致性较好，PGA实测值相比3种衰减关系给出的预测值有一定程度偏大，主要是由于衰减关系给出的预测值为基岩地震动峰值，强震动台站基本布设在土层或结构物底层，受场地影响PGA会出现放大效应. 值得注意的是，大部分观测值较雷建成等（2007）和区划图衰减关系给出的预测值偏大，并且随着距离的增加，偏大趋势愈加明显. 这一现象表明，在此次地震中所体现的高频地震动衰减要慢于这两个衰减关系所采用的路径定标关系. 当然地震动具有较大随机性，受震源、路径和场地的综合作用影响，验证衰减关系预测结果合理性还需要更多次地震的观测记录. 图5还给出了实测PGV值与五代图青藏区PGV衰减关系（M≤6.5，长轴）预测值的对比情况，可以看出在80 km以内衰减关系预测效果较好，观测值基本分布在预测中位值加减一倍标准差范围内. 但在80～200 km范围内观测值分布较为离散，整体上要大于预测值，甚至出现了一个量级的偏差. 这是由于受场地放大效应影响，同时也与区域内非弹性衰减的横向不均匀性有关，具体原因还需要下一步工作深入分析.
-
本文简要介绍了北川MS4.7获取的532组地震强震动记录，从幅值、空间分布特征及衰减情况分析了地震动特征. 从地震动强度上看，表现出了震级小、仪器烈度高和震害轻的特点，震中附近2个台站仪器烈度相比震害调查烈度基本偏高1°左右. 51BCB的土层强震动记录和51BCQ的基岩记录在反应谱的对比结果体现了基岩场地下垂直向地震动记录放大效果有限，有助于基于垂直向基岩地震动不放大条件假设的谱比法（荣棉水等，2016）在分析地震动特征上的应用. 地震动幅值的空间分布长轴呈北西—南东方向，但在多个区域出现异常增大现象，目前分析主要受到局部场地影响，由此也建议台站场址勘选应避开局部不利地形条件.
北川MS4.7地震获取了大量强震动记录，其衰减特性与俞言祥和汪素云（2006）提出的中国西部地区水平向基岩加速度衰减关系一致性较好，相比第五代区划图建立的青藏区地震动衰减关系和雷建成等（2007）提出的西南地区衰减关系偏大，主要原因为衰减关系预测值为基岩地震动，实测值受场地影响导致，同时也为后期建立四川地区区域衰减关系提供了较好数据基础. 同时，北川MS4.7地震所获得的密集的、高质量的强震动记录，也可为以小震作为子震的经验格林函数方法（EGFM）再现大地震的强地震动场展布提供重要的数据支撑.
致谢
四川地震台苏金蓉研究员及其团队为本研究提供了北川MS4.7地震震源机制解和历史地震数据，四川省地震应急服务中心毛利工程师为本研究提供了调查资料，中国地震局工程力学研究所王宏伟助理研究员为本文部分图件的绘制提供了技术指导，在此一并感谢.
2020年10月22日四川北川MS4.7地震强震动特征分析
Strong-motion recordings in MS4.7 Beichuan, Sichuan earthquake on October 22, 2020
-
摘要: 2020年10月22日11时03分37秒四川省绵阳市北川县发生MS4.7地震，四川强震动台网与预警烈度速报台网在震区建成较密集的台站，获取了532组三分量加速度记录，有助于开展区域地震动衰减和地震动特征研究. 本文对强震记录进行常规处理后计算出强震动记录的相关参数，利用克里金插值方法得到地震动峰值加速度PGA和峰值速度PGV的空间分布图，长轴呈北西—南东方向. 分析强震动记录PGA、PGV 随距离的衰减规律，与常用衰减关系预测值进行对比，此次地震 PGA 的衰减特性与
俞言祥和汪素云（2006） 提出的中国西部地区水平向基岩加速度衰减关系有较好的一致性. 北川MS4.7地震获得的密集强震动记录为建立区域衰减关系，以及开展基于经验格林函数方法（EGFM）再现大震强地震动场展布等研究提供了重要的数据支撑.Abstract: On October 22, 2020, an earthquake with magnitude 4.7 occurred in Beichuan County, Sichuan Province. 532 sets of three-component acceleration records have been obtained by Sichuan strong motion network, which contribute to the study of regional ground motion attenuation and ground motion characteristics. After the conventional processing of strong motion records, the relevant parameters of strong motion records have been calculated. The spatial distribution maps of PGA and PGV are obtained by Kriging interpolation method, and the long axis is NE-SW direction. The attenuation law of PGA and PGV of strong motion records with distance is analyzed. Through comparison and analysis with three commonly used attenuation relationship prediction values, such as attenuation relationship of Eastern strong earthquake area established by the fifth generation zoning map, it can be seen that PGA and PGV of strong motion records are similar to those of this earthquake. The attenuation characteristics of the model are in good agreement with the attenuation relationship of horizontal bedrock acceleration proposed byYu and Wang (2006) . The dense strong motion records obtained from the Beichuan MS4.7 earthquake provide important data support for establishing regional attenuation relationships and carrying out studies based on the empirical Green function method (EGFM) to reproduce the distribution of strong ground motions of large earthquakes.-
Key words:
- Beichuan earthquake /
- strong-motion recording /
- attenuation characteristic
-
图 1 北川MS4.7地震中获取记录的强震动台站分布图
Figure 1. Strong-motion stations recorded the Beichuan MS4.7 earthquake
图 2 强震动记录数量与震中距统计图
Figure 2. Statistical chart of the number of strong motion records and epicenter distance
图 3 51BCB台站和51BCQ台站加速度时程与5% 阻尼比拟加速度反应谱
Figure 3. Acceleration time histories and 5%-damped pseudospectral accelerations for the recorded at stations 51BCB and 51BCQ
表 1 地震动记录及相关参数表
Table 1. Seismic records and related parameters
台站代码 纬度/(°) 经度/(°) 震中距/km 峰值加速度/（cm·s−2） 峰值速度/（cm·s−1） 仪器烈度 EW NS UD EW NS UD 51BCB 31.8 104.2 2.9 −277.53 −250.72 103.51 3.54 4.12 0.98 7.1 B2603 31.9 104.1 10.9 158.11 180.88 142.52 −3.85 −4.28 1.66 6.7 B0571 31.7 104.3 12.5 156.39 152.27 84.13 2.81 −2.27 −1.76 6.2 U2302 31.8 104.0 18.8 −84.89 58.35 −113.25 −1.05 0.46 −1.51 5.7 B2607 32.0 104.2 21.5 30.62 23.74 21.00 1.11 1.11 0.68 4.7 B0502 31.7 104.4 24.8 23.86 10.06 19.62 −1.62 0.49 0.67 4.5 B2604 32.0 104.0 26.5 40.26 36.04 −35.32 −0.87 −0.95 0.72 4.8 51BCQ 31.8 104.5 26.7 11.81 7.34 6.59 0.36 0.22 0.19 3.1 B2608 31.9 104.5 27.3 71.05 46.59 20.26 −1.73 0.89 −0.42 5.5 51BCY 31.7 104.4 29.1 64.28 37.23 −18.70 1.07 −0.69 −0.27 5.0 U2310 31.8 103.9 30.3 39.08 −52.45 −30.54 0.51 −1.07 0.27 4.8 F8303 31.6 104.1 30.3 −53.09 57.13 27.28 −1.22 1.46 −0.62 5.1 B0501 31.6 104.3 31.4 44.66 −17.62 43.60 0.50 −0.40 −0.51 4.4 B2606 31.8 104.6 36.4 −59.36 −35.03 16.99 −2.36 −0.96 −0.37 5.4 B0503 31.6 104.4 38.1 −25.13 −9.91 −37.91 −0.79 0.54 −1.32 4.9 B2601 31.9 104.6 39.4 28.32 −23.96 9.10 −2.01 0.76 0.36 4.8 B2602 32.0 104.5 39.8 15.67 13.63 −12.76 −0.14 −0.16 −0.17 3.0 F8301 31.5 104.2 41.1 −15.99 −15.52 −6.07 0.53 −0.55 0.29 3.7 U2301 31.8 103.7 42.1 21.87 30.39 −29.72 −0.19 0.40 −0.10 3.9 B8112 31.7 104.6 42.8 −19.51 −13.82 −5.24 0.95 0.59 −0.30 4.1 B2713 32.2 104.5 46.3 16.53 17.99 8.34 −0.23 −0.32 0.14 3.4 U2306 31.9 103.7 46.6 24.66 34.78 35.83 −0.30 −0.37 0.17 3.9 B2703 32.1 104.5 47.6 46.73 48.69 −28.52 −0.19 −0.42 0.13 4.3 51MZX 31.4 104.2 48.0 6.60 −4.60 −3.28 −0.45 −0.30 −0.15 3.0 B8106 31.7 104.7 49.5 −6.32 11.61 −2.67 0.48 0.28 −0.13 3.2 
下载: 导出CSV
-
[1] 胡聿贤,张敏政.缺乏强震观测资料地区地震动参数的估计方法[J] .地震工程与工程振动,1984,4(1):1-11. Hu Y X, Zhang M Z. A Method of Predicting Ground Motion Parameters for Regions with Poor Ground Motion Data[J]. Earthquake Engineering and Engineering Vibration, 1984,4(1):1-11(in Chinese). [2] 金星,马强,李山有.利用数字强震仪记录实时仿真地动速度[J] .地震工程与工程振动, 2004, 24(1):49- 54. Jin X, Ma Q, Li S Y. Real-time simulation of ground velocity using digital accelerograph record[J]. Earthquake Engineering and Engineering Vibration, 2004, 24(1):49- 54(in Chinese). [3] 雷建成, 高孟潭, 俞言祥. 四川及邻区地震动衰减关系[J]. 地震学报, 2007, 29(5): 500-511. doi: 10.3321/j.issn:0253-3782.2007.05.007 Lei J C, Gao M T, Yu Y X. Seismic motion attenuation relations in Sichuan and adjacent areas[J]. Acta Seismologica Sinica, 2007, 29(5): 500-511(in Chinese). doi: 10.3321/j.issn:0253-3782.2007.05.007 [4] 任叶飞, 温瑞智, 周宝峰, 等.2013年4月20日四川芦山地震强地面运动三要素特征分析[J].地球物理学报,2014,57(6):1836-1846. Ren Y F, Wen R Z, Zhou B F, et al.The characteristics of strong ground motion of Lushan earthquake on April 20,2013[J]. Chinese Journal of Geophysics, 2014, 57(6):1836-1846(in Chinese). [5] 荣棉水,李小军,王振明, 等.HVSR方法用于地震作用下场地效应分析的适用性研究[J].地球物理学报,2016, 59(8):2878-2891 Rong M S, Li X J, Wang Z M, et al. Applicability of HVSR in analysis of site-effects caused by earthquakes[J]. Chinese Journal of Geophysics, 2016,59(8):2878-2891(in Chinese). [6] 王文才, 李佐唐, 党红, 等.2017年8月8日四川九寨沟7.0级地震强震记录及特征分析[J].地震工程学报,2017,39(4):652- 656. Wang W C, Li Z T, Dang H, et al. Strong motion records and its characteristics in Jiuzhaigou MS7.0 earthquake on August 8,2017[J].China Earthquake Engineering Journal, 2017,39(4):652-656(in Chinese). [7] 温瑞智, 王宏伟. 2020年7月12日唐山古冶MS5.1地震的强震动特征[J]. 地震科学进展, 2020, 50(7): 9-13. Wen R Z, Wang H W. Strong-motion recordings in MS5.1 Guye, Tangshan earthquake on July 12, 2020[J]. Progress in Earthquake Sciences, 2020, 50(7): 9-13(in Chinese). [8] 于海英, 江汶乡, 解全才, 等.近场数字强震仪记录误差分析与零线校正方法[J]. 地震工程与工程振动, 2009, 29(6):1-12. Yu H Y, Jiang W X, Xie Q C, et al. Baseline correction of digital strong-motion records in near-field[J]. Earthquake Engineering and Engineering Vibration, 2009, 29(6):1-12(in Chinese). [9] 俞言祥, 汪素云. 中国东部和西部地区水平向基岩加速度反应谱衰减关系[J]. 震灾防御技术, 2006, 3(1): 206-217. Yu Y X, Wang S Y. Attenuation relations for horizontal peak ground acceleration and response spectrum in eastern and western China[J]. Technology for Earthquake Disaster Prevention, 2006, 3(1): 206-217(in Chinese). [10] 俞言祥, 李山有, 肖亮. 为新区划图编制所建立的地震动衰减关系[J]. 震灾防御技术, 2013, 8(1): 24-33. doi: 10.3969/j.issn.1673-5722.2013.01.003 Yu Y X, Li S Y, Xiao L. Development of ground motion attenuation relations for the new seismic hazard map of China [J]. Technology for Earthquake Disaster Prevention, 2013, 8(1): 24—33(in Chinese). doi: 10.3969/j.issn.1673-5722.2013.01.003 -
点击查看大图
图(5) / 表(1)
计量
- 文章访问数: 163
- HTML全文浏览量: 121
- PDF下载量: 14
- 被引次数: 0
