Study on the rupture directivity of the 2021 Yangbi earthquake sequence
-
摘要:
基于2021年云南漾濞MS6.4地震序列的强震动记录，建立了地震动参数预测方程，采用破裂方向性效应拟合方法估计了强震动记录丰富且空间分布均匀的九次地震的震源破裂方向性特征。结果表明：其中的四次地震（1号、4号、8号、9号）表现出破裂方向性效应，且均为不均匀双侧破裂，但优势破裂方向不同（4号和8号为东南向，1号西北向，9号西南向），说明2021年漾濞地震序列的地震破裂较为复杂；由于主震（4号地震）的破裂速度较慢（约为2.2 km/s），其破裂方向性效应较弱且主要对峰值速度有影响，其它三次地震破裂方向性效应十分显著，破裂速度大于主震；此外，四次地震的破裂方向性效应还存在一定的周期相关性。
Abstract:This study used the strong-motion recordings of the MS6.4 Yangbi earthquake sequence in Yunnan Province to establish the prediction equations of ground motion parameters, and then estimated the focus rupture directivities characteristics for nine Yangbi earthquakes with abundant recordings in a good spatial station coverage based on the rupture directivity effect fitting method. The results indicated that the rupture directivity effects are observed in four of these earthquakes （i.e., the 1st, 4th, 8th, and 9th earthquakes）. The four earthquakes are all characterized by the bilateral ruptures with various predominant rupture directions （i.e., southeast for the 4th and 8th events, northwest for the 1st event, and southwest for the 9th event）, which illustrates the rupture complexities of the Yangbi earthquake sequence. Since the estimated rupture velocity （about 2.2 km） of the mainshock （the 4th event） is very slow, its rupture directivity effects are not strong and mainly affect the peak ground velocity. However, the rupture directivity effects for the other three earthquakes are very strong and the rupture velocities are faster than that of the mainshock. The dependency of the rupture directivity effects on the period was also observed in the four earthquakes.
-
引言
高精度重力梯度测量作为一种新兴的地球物理探测技术，不仅在海洋重力辅助导航、 海底潜艇探测和海底资源勘探等领域具有重要意义(徐遵义等，2007； DiFrancesco et al，2009； Araya et al，2011)，而且在陆地资源勘探、 大型水库监测和水底坑洞探测等领域也是一种有效的地球物理探测技术(王谦身，2003; Woolrych et al，2015)． 此外，航空和卫星重力梯度测量在地球重力场等基础科学研究领域也具有重要的作用(宁津生等，1996，2002)．
自1971年美国空军首次提出精度为1 E的移动级重力梯度仪以来，高精度重力梯度仪的研究得到了迅速发展，其主要设计原理为差分加速度计和扭矩测量． 基于差分加速度计的航空重力梯度仪，由于其自身的高稳定性和高精度得到了迅速发展和应用(李红军等，2002； 刘凤鸣等，2009)，但大部分还处于实验样机阶段，目前唯一定型的商用重力梯度仪是美国贝尔实验室研制的旋转加速度计重力梯度仪(DiFrancesco et al，2009)． 基于扭矩测量的重力梯度仪，由于体积和稳定性问题，其发展受到了限制． 彭益武等(2006)设计的基于扭矩测量的二维簧片重力梯度仪在体积和稳定性等方面取得了突破． 另外，超导重力梯度仪由于其测量精度、 体积和整体质量等方面的优势，越来越多地被应用于航空重力测量中(吴琼等，2013)，具有非常好的发展前景(Moody et al，2002)．
上述各类梯度仪的设计原理均属于相对测量，需要对参考坐标及变化的物理量进行标定，不适用于发射到宇宙外星球进行长期的内部质量分布探测研究，因此，在美国宇航局20世纪末的火星探测计划中，Micro-G公司开展了基于激光干涉原理的绝对重力梯度仪的研制，其测量精度的设计指标为1 E，其目的是探测火星的深部结构(Brown et al，1999; Robertson，2001)．
受此启发，激光绝对重力梯度仪作为研究地球深部质量迁移的一种新型有效工具，更能满足地震的监测需求． 目前，我国差分加速度计的最高精度为10-6 m/s2，这与梯度仪要求的测量精度还有很大差距； 而现有的绝对重力测量技术的测量精度可达10-8 m/s2，完全能够满足重力梯度仪的精度要求(胡华等，2012； 滕云田等，2013)，故基于创新设计及软件开发，利用激光干涉法完成地表重力垂直梯度测量的方案是完全可行的． 虽然目前该仪器的体积较大，但从固定台站观测的角度看，该缺点可以忽略； 而且随着激光干涉重力梯度仪研究的进一步开展，将来可以进行小型化设计，使之满足移动重力垂直梯度仪的测量需求．
本文拟基于激光干涉法对新型地表重力垂直梯度测量系统进行研究并初步构建原理样机. 该测量系统采用激光干涉原理，通过测定两个落体同步自由下落时相对于刚性连接的两个参考点的位置信息确定各自的重力加速度，然后通过差分计算得到测点的重力垂直梯度． 由于两个参考点为刚性连接，其感受到地面的振动信息完全相同，且两个落体利用同一套伺服控制系统，由振动引入的测量误差将在差分时消除，这样即可得到某测点的高精度重力垂直梯度值． 与其它重力梯度仪相比，激光干涉重力梯度仪具有以下优点： ① 有较大的测量范围，无需标定，无漂移； ② 地面振动属于共模误差，无需设计特殊的隔振系统，对动态测量非常有利； ③ 地面倾斜属于共模误差； ④ 对地球浅层异常质量体敏感．
1. 设计原理
1.1 基本原理
基于激光干涉法重力垂直梯度测量的基本原理是，当两个落体在一个真空舱中同步自由下落时，运动过程中所受到的干扰因素可以等效为一个合力，并且该合力对每个落体的作用均相同，唯一不同的是所受重力场重力垂直梯度的影响． 因此，通过解算垂直向间隔为h的两个落体(较高位置的为上落体P1，较低位置的为下落体P2)在真空中作自由下落运动时测量得到的干涉条纹，分析重建各自的自由下落轨迹(吴琼等，2012)，并获取上落体P1在其有效测量高度href_U处的第一绝对重力加速度值gU和下落体P2在其有效测量高度href_D处的第二绝对重力加速度值gD，再根据下式得到单次下落的重力垂直梯度γ，即
(1) 图 1给出了重力梯度测量算法原理，主要是为了确定式(1)中的两个未知量href_U和href_D． 落体P1的初始位置为A00，其对应的初始速度和位移均为0，但在实际工作中一般将初始速度和位移无限接近于0的A0点作为其初始位置，在计算中近似认为A00与A0重合，因此hinst_U对应的是P1的初始位置． 在heff_U位置处，由于上落体P1在自由下落过程中所受到的重力梯度影响的综合效应为0，则测量得到的重力加速度值所对应的高度为
(2) 对于下落体P2，同样有
(3) ![]()
图 1 重力梯度测量算法原理Figure 1. Principle for the algorithm of the gravity vertical gradient measurement在计算重力垂直梯度时，首先计算heff_U和heff_D，然后将式(2)和(3)代入式(1)求得最终测点的垂直向重力梯度值． 对于heff_U和heff_D，二者采用的推导过程相同，这里仅以heff_U为例进行说明：
1)设定初始条件． 考虑到重力垂直梯度的测量精度要求，在heff_U的计算过程中，其计算精度须精确到mm； 对于上落体P1，从初始位置A00到位置A0的下落运动过程为非自由下落运动，距离一般小于5 μm，时间小于30 ms． 这段非自由下落运动过程对heff_U值的影响在误差允许范围内，因此在计算中可被忽略． 也就是说，假设上落体在t00时刻从初始位置A00开始进入自由下落运动，则t0=0，v0=0，x0=0； 初始位置A00的绝对重力加速度g0为已知设定值． 此外，(t1，x1)和(tn，xn)已根据所获取的干涉信号解算得到．
2) 落体的自由下落运动方程为
(4) 根据假设的初始条件，解上述方程，可得
(5) 由于gU是t1—tn时段内重力场平均作用于上落体P1的结果，则其表达式为
(6) 将式(5)代入式(6)，得
(7) 又由于gU为某个位置的绝对重力加速度值，将与其对应的位置记为C，则初始位置A00到位置C的距离记为heff_U，则
(8) 对比式(7)与式(8)，可得到
(9) 由式(9)可知，在利用同一套落体伺服控制系统对上、 下落体进行同步控制的过程中，heff_U与heff_D在对应的开始时间t1和结束时间tn的取值完全相同，即heff_U=heff_D，故将式(2)和式(3)代入式(1)，得
(10) 即利用式(1)进行重力垂直梯度计算时与g0的选取无关，仅取决于上、 下落体自由下落时重力加速度的计算和结构设计时所确定的两个落体光心之间的距离． 结构设计过程中，设上、 下落体光心间的水平距离为95 mm，若水平梯度按300—500 E(朱英，李卢玲，1982)计算，梯度测量误差为3—5 E，而本套梯度测量系统精度为100 E，因此水平梯度引入的误差可以暂时忽略，故式(1)中上、 下落体的光心间距即为上、 下落体光心间的垂直距离．
1.2 光路设计
激光干涉法重力垂直梯度测量系统的光路设计如图 2所示． 可以看到，激光器发出激光束入射至分光镜BS1的上表面，分别得到BS1上表面透射和反射两束同源激光，透射光入射至分光镜BS2的上表面，反射光经反射镜M反射后入射至分光镜BS3上表面． 分光镜BS2、 下落体P2和参考镜MS2等构成下落体P2的干涉测量光路，所形成的干涉条纹利用雪崩光电二极管(advanced photo diodes，简写为APD，下文分别用APD1和APD2表示上落体P1和下落体P2自由下落时形成的干涉条纹的接收器)接收并转换为电压信号被后端的高速数字化仪接收； 分光镜BS3、 上落体P1和参考镜MS1等构成P1的干涉测量光路，形成的干涉条纹被APD1接收并转换为电压信号被后端的高速数字化仪接收(吴琼等，2011)．
![]()
图 2 干涉法重力梯度测量原理Figure 2. Measuring principle of the gravity vertical gradient based on interferometry这套装置中，激光器、 分光镜(BS1，BS2，BS3)、 反射镜(M)、 参考镜(M1，M2)和光电接收器(APD1，APD2)均为刚性连接，真空舱中的P1和P2在同一套伺服电机控制系统作用下作自由下落运动，因此APD1和APD2接收到的干涉条纹所受到的地震动、 真空度、 气压、 极移和潮汐等环境干扰均保持相同． 源自这些因素的误差在测点重力梯度的差分计算时将被自动消除，不会影响最终测量结果的精度．
2. 测量系统设计
2.1 双落体伺服控制系统结构设计
基于激光干涉法的重力垂直梯度测量系统是利用一套电机伺服控制系统控制垂直方向间距为50 cm的两个落体作同步自由下落运动，故需要构建一个如图 3所示的真空系统，为落体的自由下落提供高真空环境，并为双落体伺服控制系统提供机械支撑．
真空系统安装在由3个支撑腿支撑的角架上，其中心旋转轴通过角架的质心，以确保角架完成水平调节后，真空系统的中心旋转轴为竖直状态． 该系统由上盖、 上二通底座、 上真空舱、 下真空舱和下二通底座等构成． 由于真空系统腔体较大，在构建真空环境时将分子泵直接安装在下二通底座上，并将其预抽至2×10-4 Pa以上时打开离子泵阀门，维持真空腔体内真空度优于5×10-4 Pa，满足落体自由下落对真空度的要求．
图 4为图 3所示真空系统内部安装的落体伺服控制系统． 图 4a显示的是落体伺服控制系统的整体结构，电机的输出力矩通过大气侧齿轮组和过真空传动机构驱动真空侧的齿条带动双落体作上下往复运动； 齿条经齿条运动导座约束后，可确保其沿竖直方向上下运动． 上托架内放置上落体P1，齿条通过万向联轴节驱动上托架，上托架沿上导轨作上下往复运动，对于下落体P2的驱动也采用同样的控制方式． 图 4b所示结构是为保证上、 下落体驱动的同步性而进行的特殊设计，由连接板、 连接块、 连接柱以及万向联轴节等构成． 连接板通过4个螺钉与上托架底端固定，末端与连接块固定； 连接块通过螺母锁定连接柱，三者构成“7”型结构； 连接柱通过万向联轴节与下托架固定． 该设计不仅可以保证上、 下托架在电机驱动下同步运动，还可抵消安装误差和制造误差所造成的径向力作用，使整套控制机构运行平稳可靠．
![]()
图 4 落体伺服控制系统设计(a)安装在真空舱内的落体伺服控制系统设计；(b)上、 下落体间刚性连接机构设计Figure 4. Design of the free-fall body servo control mechanism(a)Design of the free-fall body servo control mechanism in vacuum chamber;(b)Design of the rigid linkage between the up and down free-fall bodies2.2 双光路干涉测量系统设计
双光路干涉测量系统首先通过两个干涉光路测量单元分别采集上落体和下落体在自由下落过程中形成的激光干涉信号，再通过对其解算和差分运算得到地表重力垂直梯度．
依据图 2所示的测量原理设计双光路干涉测量系统的机械结构，如图 5所示． 激光器采用碘稳频激光系统，波长稳定度可达10-12，为整套测量系统提供长度基准． 高速数字化仪的时间基准由外部铷原子时钟提供，铷原子时钟输出频率的稳定度可达10-11，为整套仪器提供时间基准．
![]()
图 5 双光路激光干涉测量系统机构设计Figure 5. Design of the measuring structure of dual optical path measuring system based on laser interferometry干涉条纹的处理解算方面，由于上、 下落体自由下落的时间和位移信息可追溯到时间和长度基准，所以最后通过差分获得的重力垂直梯度值就具有了大动态范围、 无需标定和无零漂等特点； 构成双光路干涉测量系统所需的各个部件均刚性安装在双光路干涉测量系统的底板上，因此振动产生的误差属于共模误差，不会对最终的梯度测量结果产生影响．
光束垂直性调节方面，与激光干涉绝对重力仪不同，双光路干涉测量系统需要对两条测量光束进行垂直性调节． 首先，利用酒精液面提供的水平面和光学自准直仪，通过调整双光路干涉测量系统底板的支撑腿，调整入射上落体P1的测量光束的垂直性； 然后，同样利用酒精液面提供的水平面和光学自准直仪，通过调整反射镜M上的微调螺母，调整入射下落体P2的测量光束的垂直性． 这种调整方案的误差与光学自准直仪的误差相同，可以保证入射落体P1和P2的测量光束的垂直性误差小于0.2″．
2.3 测量系统构建与初步试验
目前，根据2.1和2.2节的设计已构建完成了一套可以正常工作的激光干涉法地表重力垂直梯度测量系统，如图 6a所示． 由于这套测量系统机构比较庞大，目前仅在实验室获取了初步的试验数据，如图 6b所示． 试验数据共有35组，每组进行16次测量.
![]()
图 6 激光干涉法地表重力垂直梯度测量系统(a)及其初步试验数据(b)Figure 6. The apparatus of the gravity vertical gradient measurement based on laser interferometry(a) and its preliminary test data(35 groups and 16 measurements per group)(b)在落体自由下落过程中，选取0.045—0.130 s时段的干涉信号进行计算． 对应落体自由下落距离为1—8 cm，生成的干涉信号是调频正弦波，频率为1.4—4.0 MHz(吴琼等，2012)． 考虑到采样频率引入的时间测量误差，选取采样频率为100 MHz，使得由采样引入的时间误差导致最终的测量误差小于1×10-8 m/s2．
图 7给出了单次测量中上落体自由下落时生成的干涉信号及其展开，可以看出，获取的干涉信号无明显干扰信号，信噪比高，可以满足高精度重力加速度值计算的要求．
![]()
图 7 实测上落体自由下落时得到的干涉条纹(a)及其展开(b)Figure 7. The interference fringes of the upper free-fall body during freely falling(a) and the expansion of the part of these fringes(b)3. 讨论与结论
在重力垂直梯度的测量过程中，潮汐、 极移、 气压、 真空度、 温度梯度、 涡流以及参考棱镜的振动等相较于明显的误差源，在绝对重力测量时不会影响重力垂直梯度的测量结果，但双落体自由下落的同步性控制、 参考点以及各个光学部件间的刚性和双测量光路的垂直性等则是影响其测量精度的重要因素．
2012年至今，在原有设计研发激光干涉重力仪的基础上，初步完成了激光干涉法地表重力垂直梯度测量系统的原理研究、 结构设计和系统构建等，并获取了本套测量系统的第一批试验数据． 从图 7所示的干涉信号可以看出，本套测量系统获取的上、 下落体在自由下落过程中所产生的干涉信号的信噪比高，干涉条纹清晰完整、 无畸变，完全满足后期的计算要求． 图 6b显示的是本次试验完成的35组测量结果的组平均值，试验中每组设定的测量次数为16次，共完成560次测量． 对测量获得的35组平均值再进行求平均计算，确定本次试验样机测量的测点的垂直向重力梯度为2886 E，精度为99 E．
2012年11月曾利用两台拉科斯特相对重力仪对本次试验所在测点的重力垂直梯度进行测量，得到的梯度值为2402 E(精度100 E)，对比分析初步确定本套梯度测量系统的测量准确度约为400 E．
本套梯度测量系统的精度和准确度不仅需要更多试验数据的测定，还需要制定更完备的测试精度和准确度的试验方案． 另外，针对伺服控制系统运行过程中的降噪问题、 数据处理算法和整套仪器的小型化设计等方面尚需进一步改进．
-
![]()
图 1 （a） 用来建立地震动预测方程的记录所对应的震级-距离分布图；（b）可用记录数随周期的变化
Figure 1. （a） Magnitude-distance distribution of the recordings for developing ground motion prediction equations； （b） The number of the usable recording varied with the period
![]()
图 2 PGA和PGV观测值随距离衰减及其预测中位值
Figure 2. The distance attenuation of the observed PGA and PGV values with distance and the predicted medians
![]()
图 3 本文选用的9次地震及其记录的台站空间分布
Figure 3. The nine earthquakes analyzed in this study and the spatial distribution of stations for each event
![]()
图 4 PGA距离校正事件内残差随方位角的变化以及Cd最佳拟合结果
Figure 4. The path-corrected intra-event residuals varied with azimuth for PGA in each event and the azimuthal dependence of the best-fitted Cd
![]()
图 5 PGV距离校正事件内残差随方位角的变化以及Cd最佳拟合结果
Figure 5. The path-corrected intra-event residuals varied with azimuth for PGV in each event and the azimuthal dependence of the best-fitted Cd
![]()
图 6 基于不同地震动强度指标给出的$\max C_{\rm{d}}^{0.5}/\min C_{\rm{d}}^{0.5} $，纵坐标数值表示不同周期的PSA
Figure 6. The $\max C_{\rm{d}}^{0.5}/\min C_{\rm{d}}^{0.5} $ values based onvarious ground motion intensity measures，numbers on the ordinate indicate the period of PSA
![]()
图 7 基于不同周期PSA估计的破裂方向及其标准差范围
Figure 7. Rupture directions and standard deviations based on PSAs at various periods
表 1 地震动参数预测方程回归系数及标准差
Table 1 Regression coefficients and standard deviations of the parameter prediction equations of ground motion
地震动强度指标 回归系数 a1 a2 a3 a4 a5 a6 a7 φ τ σ PGA 2.907 4 −0.089 5 0.028 5 −1.719 0 0.142 1 −0.009 6 −0.085 5 0.340 8 0.192 4 0.391 4 PGV 2.226 0 −0.106 2 0.046 6 −1.481 1 0.099 0 −0.008 6 −0.521 2 0.358 9 0.244 6 0.434 3 PSA T＝0.10 s 2.581 5 0.066 6 0.026 1 −1.456 6 0.043 9 −0.005 4 −0.023 1 0.357 2 0.171 4 0.396 2 T＝0.15 s 2.650 7 0.082 0 0.025 9 −1.081 5 0.041 3 −0.008 2 −0.241 3 0.368 8 0.205 0 0.421 9 T＝0.20 s 2.818 1 0.090 8 0.029 2 −1.117 1 0.027 4 −0.006 6 −0.333 8 0.419 9 0.230 2 0.478 9 T＝0.26 s 2.902 3 0.030 4 0.042 1 −0.963 7 −0.000 5 −0.006 7 −0.405 1 0.456 1 0.237 0 0.514 0 T＝0.30 s 2.764 9 0.046 3 0.044 1 −0.900 2 −0.013 7 −0.006 5 −0.426 5 0.448 0 0.230 3 0.503 7 T＝0.36 s 3.021 2 0.068 7 0.043 0 −1.003 7 −0.008 1 −0.005 4 −0.556 9 0.424 4 0.224 9 0.480 3 T＝0.40 s 3.165 0 0.103 1 0.039 3 −1.1328 0.005 1 −0.005 3 −0.624 6 0.420 0 0.223 8 0.475 9 T＝0.46 s 2.810 6 0.198 6 0.037 9 −0.977 2 −0.050 5 −0.003 7 −0.637 0 0.404 7 0.230 9 0.465 9 T＝0.50 s 2.657 4 0.251 2 0.034 4 −0.949 4 −0.059 6 −0.004 0 −0.646 8 0.403 3 0.232 5 0.465 5 T＝0.60 s 2.190 3 0.343 8 0.036 4 −0.735 3 −0.130 6 −0.001 5 −0.654 1 0.402 5 0.244 5 0.470 9 T＝0.70 s 2.026 6 0.344 0 0.035 3 −0.619 7 −0.124 0 −0.003 1 −0.677 0.400 3 0.241 0 0.467 2 T＝0.90 s 1.713 3 0.364 3 0.034 4 −0.429 3 −0.142 6 −0.002 0 −0.721 1 0.403 5 0.246 0 0.472 6 T＝1.00 s 1.416 7 0.430 2 0.028 6 −0.294 7 −0.147 7 −0.003 0 −0.749 4 0.411 8 0.232 4 0.472 9 T＝1.50 s 1.015 6 0.548 5 0.018 0 −0.079 2 −0.155 4 −0.001 9 −0.919 5 0.457 7 0.260 7 0.526 7 
下载: 导出CSV
表 2 震级相关的h值
Table 2 Magnitude-dependent h values
M h M h M h M h 2.8—3.5 1.00 3.6—3.7 1.05 4.0—4.5 1.10 4.6—4.7 1.20 5.0 1.50 5.2 1.50 5.6 2.00 6.4 3.90
下载: 导出CSV
表 3 选用的9次地震的基本信息及记录数
Table 3 Basic information and the number of recordings for the nine events considered
地震编号 发震时刻（北京时间） MS 北纬/° 东经/° 震源深度/km 记录数 年-月-日 时:分:秒 1 2021-05-19 20:05:56 4.6 25.65 99.91 10 50 2 2021-05-21 21:21:25 5.6 25.65 99.92 10 96 3 2021-05-21 21:21:57 4.2（ML） 25.63 99.96 10 50 4 2021-05-21 21:48:34 6.4 25.70 99.88 10 100 5 2021-05-21 21:55:28 5.0 25.67 99.89 9 62 6 2021-05-21 22:31:10 5.2 25.61 99.97 8 97 7 2021-05-21 23:23:34 4.5 25.59 99.98 9 72 8 2021-05-22 02:28:43 4.2（ML） 25.63 99.92 19 37 9 2021-05-22 20:14:36 4.7 25.60 99.92 10 63
下载: 导出CSV
表 4 基于PGA，PGV分别估计的地震破裂方向性参数
Table 4 Rupture directivity parameters estimated based on PGA and PGV，respectively
地震编号 峰值参数 φ/° vr/β k $ {\max C_{\rm{d}}^{0.5} }$ ${\min C_{\rm{d}}^{0.5} } $ ${{\max C_{\rm{d} }^{0.5} }/{\min C_{\rm{d} }^{0.5} } }$ 1 PGA 316.0±11.4 0.93±0.04 0.96±0.03 3.73 0.75 4.95 PGV 309.0±33.6 0.70±0.16 0.86±0.10 1.70 0.80 2.14 2 PGA 296.1±51.4 0.61±0.16 0.68±0.10 1.32 0.83 1.59 PGV 220.6±38.4 0.64±0.10 0.88±0.14 1.56 0.79 1.97 3 PGA 119.0±69.9 0.58±0.17 0.27±0.16 1.31 0.83 1.59 PGV 130.3±104.2 0.64±0.19 0.27±0.15 1.44 0.82 1.75 4 PGA 233.1±58.6 0.63±0.14 0.78±0.14 1.45 0.82 1.78 PGV 167.1±6.8 0.62±0.04 0.82±0.04 1.42 0.81 1.82 5 PGA 247.1±108.1 0.51±0.14 0.71±0.15 1.21 0.82 1.47 PGV 197.4±82.4 0.68±0.18 0.72±0.18 1.50 0.83 1.81 6 PGA 99.4±55.6 0.66±0.17 0.19±0.17 1.55 0.81 1.92 PGV 130.2±64.9 0.70±0.17 0.24±0.17 1.59 0.82 1.94 7 PGA 114.6±41.9 0.63±0.13 0.34±0.14 1.34 0.83 1.61 PGV 97.1±51.4 0.64±0.15 0.22±0.14 1.47 0.82 1.80 8 PGA 149.8±21.9 0.76±0.10 0.87±0.11 1.90 0.80 2.39 PGV 155.6±19.9 0.84±0.12 0.93±0.09 2.43 0.77 3.16 9 PGA 244.4±11.2 0.65±0.08 0.87±0.06 1.57 0.80 1.98 PGV 242.5±18.7 0.70±0.13 0.87±0.08 1.73 0.79 2.18
下载: 导出CSV
-
段梦乔,赵翠萍,周连庆,赵策,左可桢. 2021. 2021年5月21日云南漾濞MS6.4地震序列发震构造[J]. 地球物理学报,64(9):3111–3125. Duan M Q,Zhao C P,Zhou L Q,Zhao C,Zuo K Z. 2021. Seismogenic structure of the 21 May 2021 MS6.4 Yunnan Yangbi earthquake sequence[J]. Chinese Journal of Geophysics,64(9):3111–3125 (in Chinese).
胡进军,谢礼立. 2011. 汶川地震近场加速度基本参数的方向性特征[J]. 地球物理学报,54(10):2581–2589. Hu J J,Xie L L. 2011. Directivity in the basic parameters of the near-field acceleration ground motions during the Wenchuan earthquake[J]. Chinese Journal of Geophysics,54(10):2581–2589 (in Chinese).
雷兴林,王志伟,马胜利,何昌荣. 2021. 关于2021年5月滇西漾濞MS6.4地震序列特征及成因的初步研究[J]. 地震学报,43(3):261–286. Lei X L,Wang Z W,Ma S L,He C R. 2021. A preliminary study on the characteristics and mechanism of the May 2021 MS6.4 Yangbi earthquake sequence,Yunnan,China[J]. Acta Seismologica Sinica,43(3):261–286 (in Chinese).
李大虎,丁志峰,吴萍萍,刘韶,邓菲,张旭,赵航. 2021. 2021年5月21日云南漾濞MS6.4地震震区地壳结构特征与孕震背景[J]. 地球物理学报,64(9):3083–3100. Li D H,Ding Z F,Wu P P,Liu S,Deng F,Zhang X,Zhao H. 2021. The characteristics of crustal structure and seismogenic background of Yangbi MS6.4 earthquake on May 21,2021 in Yunnan Province,China[J]. Chinese Journal of Geophysics,64(9):3083–3100 (in Chinese).
龙锋,祁玉萍,易桂喜,吴微微,王光明,赵小艳,彭关灵. 2021. 2021年5月21日云南漾濞MS6.4地震序列重新定位与发震构造分析[J]. 地球物理学报,64(8):2631–2646. Long F,Qi Y P,Yi G X,Wu W W,Wang G M,Zhao X Y,Peng G L. 2021. Relocation of the MS6.4 Yangbi earthquake sequence on May 21,2021 in Yunnan Province and its seismogenic structure analysis[J]. Chinese Journal of Geophysics,64(8):2631–2646 (in Chinese).
卢永坤,张建国,张方浩,杜浩国,杨黎薇. 2021. 2021年云南漾濞MS6.4地震烈度与震害特征[J]. 地震研究,44(3):429–438. Lu Y K,Zhang J G,Zhang F H,Du H G,Yang L W. 2021. The characteristics of the seismic intensity and damage of the 2021 Yangbi,Yunnan MS6.4 earthquake[J]. Journal of Seismological Research,44(3):429–438 (in Chinese).
苏金波,刘敏,张云鹏,王伟涛,李红谊,杨军,李孝宾,张淼. 2021. 基于深度学习构建2021年5月21日云南漾濞MS6.4地震序列高分辨率地震目录[J]. 地球物理学报,64(8):2647–2656. Su J B,Liu M,Zhang Y P,Wang W T,Li H Y,Yang J,Li X B,Zhang M. 2021. High resolution earthquake catalog building for the 21 May 2021 Yangbi,Yunnan,MS6.4 earthquake sequence using deep-learning phase picker[J]. Chinese Journal of Geophysics,64(8):2647–2656 (in Chinese).
杨九元,温扬茂,许才军. 2021. 2021年5月21日云南漾濞MS6.4地震:一次破裂在隐伏断层上的浅源走滑事件[J]. 地球物理学报,64(9):3101–3110. Yang J Y,Wen Y M,Xu C J. 2021. The 21 May 2021 MS6.4 Yangbi (Yunnan) earthquake:A shallow strike-slip event rupturing in a blind fault[J]. Chinese Journal of Geophysics,64(9):3101–3110 (in Chinese).
岳汉,张勇,盖增喜,王腾,赵里. 2020. 大地震震源破裂模型:从快速响应到联合反演的技术进展及展望[J]. 中国科学:地球科学,50(4):515–537. Yue H,Zhang Y,Ge Z X,Wang T,Zhao L. 2020. Resolving rupture processes of great earthquakes:Reviews and perspective from fast response to joint inversion[J]. Science China Earth Sciences,63(4):492–511. doi: 10.1007/s11430-019-9549-1
Abrahamson N A,Youngs R R. 1992. A stable algorithm for regression analyses using the random effects model[J]. Bull Seismol Soc Am,82(1):505–510. doi: 10.1785/BSSA0820010505
Abrahamson N A,Silva W J. 1997. Empirical response spectral attenuation relations for shallow crustal earthquakes[J]. Seismol Res Lett,68(1):94–127. doi: 10.1785/gssrl.68.1.94
Benioff H. 1955. Mechanism and strain characteristics of the White Wolf fault as indicated by the aftershock sequence[M]//Earthquakes in Kern County, California During 1952. California: State of California Natural Resources, Division of Mines: 199−202.
Ben-Menahem A. 1961. Radiation of seismic surface-waves from finite moving sources[J]. Bull Seismol Soc Am,51(3):401–435. doi: 10.1785/BSSA0510030401
Bernard P,Herrero A,Berge C. 1996. Modeling directivity of heterogeneous earthquake ruptures[J]. Bull Seismol Soc Am,86(4):1149–1160. doi: 10.1785/BSSA0860041149
Boatwright J. 2007. The persistence of directivity in small earthquakes[J]. Bull Seismol Soc Am,97(6):1850–1861. doi: 10.1785/0120050228
Calderoni G,Rovelli A,Ben-Zion Y,Di Giovambattista R. 2015. Along-strike rupture directivity of earthquakes of the 2009 L’Aquila,Central Italy,seismic sequence[J]. Geophys J Int,203(1):399–415. doi: 10.1093/gji/ggv275
Chen J L,Hao J L,Wang Z,Xu T. 2022. The 21 May 2021 MW6.1 Yangbi earthquake:A unilateral rupture event with conjugately distributed aftershocks[J]. Seismol Res Lett,93(3):1382–1399. doi: 10.1785/0220210241
Colavitti L,Lanzano G,Sgobba S,Pacor F,Gallovič F. 2022. Empirical evidence of frequency-dependent directivity effects from small-to-moderate normal fault earthquakes in central Italy[J]. J Geophys Res:Solid Earth,127(6):e2021JB023498.
Convertito V,Caccavale M,De Matteis R,Emolo A,Wald D,Zollo A. 2012. Fault extent estimation for near-real-time ground-shaking map computation purposes[J]. Bull Seismol Soc Am,102(2):661–679. doi: 10.1785/0120100306
Courboulex F,Dujardin A,Vallee M,Delouis B,Sira C,Deschamps A,Honore L,Thouvenot F. 2013. High-frequency directivity effect for an MW4.1 earthquake,widely felt by the population in southeastern France[J]. Bull Seismol Soc Am,103(6):3347–3353. doi: 10.1785/0120130073
Cultrera G,Pacor F,Franceschina G,Emolo A,Cocco M. 2009. Directivity effects for moderate-magnitude earthquakes (MW5.6−6.0) during the 1997 Umbria-Marche sequence,central Italy[J]. Tectonophysics,476(1/2):110–120.
Folesky J,Kummerow J,Shapiro S A,Häring M,Asanuma H. 2016. Rupture directivity of fluid-induced microseismic events:Observations from an enhanced geothermal system[J]. J Geophys Res:Solid Earth,121(11):8034–8047. doi: 10.1002/2016JB013078
Gallovič F. 2016. Modeling velocity recordings of the MW6.0 South Napa,California,earthquake:Unilateral event with weak high-frequency directivity[J]. Seismol Res Lett,87(1):2–14. doi: 10.1785/0220150042
Gong W Z,Ye L L,Qiu Y X,Lay T,Kanamori H. 2022. Rupture directivity of the 2021 MW6.0 Yangbi,Yunnan earthquake[J]. J Geophys Res:Solid Earth,127(9):e2022JB024321. doi: 10.1029/2022JB024321
Joyner W B. 1991. Directivity for nonuniform ruptures[J]. Bull Seismol Soc Am,81(4):1391–1395.
Kane D L,Shearer P M,Goertz-Allmann B P,Vernon F L. 2013. Rupture directivity of small earthquakes at Parkfield[J]. J Geophys Res:Solid Earth,118(1):212–221. doi: 10.1029/2012JB009675
Lengliné O,Got J L. 2011. Rupture directivity of microearthquake sequences near Parkfield,California[J]. Geophys Res Lett,38(8):L08310.
Lin Y Y,Lapusta N. 2018. Microseismicity simulated on asperity–like fault patches:On scaling of seismic moment with duration and seismological estimates of stress drops[J]. Geophys Res Lett,45(16):8145–8155. doi: 10.1029/2018GL078650
McGuire J J. 2004. Estimating finite source properties of small earthquake ruptures[J]. Bull Seismol Soc Am,94(2):377–393. doi: 10.1785/0120030091
McGuire J L,Zhao L,Jordan T H. 2002. Predominance of unilateral rupture for a global catalog of large earthquakes[J]. Bull Seismol Soc Am,92(8):3309–3317. doi: 10.1785/0120010293
Pacor F,Gallovič F,Puglia R,Luzi L,D’Amico M. 2016. Diminishing high-frequency directivity due to a source effect:Empirical evidence from small earthquakes in the Abruzzo region,Italy[J]. Geophys Res Lett,43(10):5000–5008. doi: 10.1002/2016GL068546
Phung V, Atkinson G M, Lau D T. 2004. Characterization of directivity effects observed during 1999 Chi-Chi, Taiwan earthquake[C]//13th World Conference of Earthquake Engineering. Vancouver, BC, Canada: 2740.
Qiang S Y,Wang H W,Wen R Z,Ren Y F,Cui J W. 2023. Characteristics of strong ground motions from four MS≥5.0 earthquakes in the 2021 Yangbi,southwest China,seismic sequence[J]. J Earthq Eng,27(14):3957–3974. doi: 10.1080/13632469.2022.2143941.
Ren Y F,Wang H W,Wen R Z. 2017. Imprint of rupture directivity from ground motions of the 24 August 2016 MW6.2 Central Italy earthquake[J]. Tectonics,36(12):3178–3191. doi: 10.1002/2017TC004673
Ross Z E,Trugman D T,Azizzadenesheli K,Anandkumar A. 2020. Directivity modes of earthquake populations with unsupervised learning[J]. J Geophys Res:Solid Earth,125(2):e2019JB018299. doi: 10.1029/2019JB018299
Ruiz J A,Baumont D,Bernard P,Berge-Thierry C. 2011. Modelling directivity of strong ground motion with a fractal,k−2,kinematic source model[J]. Geophys J Int,186(1):226–244. doi: 10.1111/j.1365-246X.2011.05000.x
Spudich P,Chiou B S J. 2008. Directivity in NGA earthquake ground motions:Analysis using isochrone theory[J]. Earthq Spectra,24(1):279–298. doi: 10.1193/1.2928225
Velasco A A,Ammon C J,Lay T. 1994. Empirical green function deconvolution of broadband surface waves:Rupture directivity of the 1992 Landers,California (MW7.3),earthquake[J]. Bull Seismol Soc Am,84(3):735–750.
Wald D J,Heaton T H,Hudnut K W. 1996. The slip history of the 1994 Northridge,California earthquake determined from strong motion,teleseismic,GPS,and leveling data[J]. Bull Seismol Soc Am,86(1B):S49–S70. doi: 10.1785/BSSA08601B0S49
Wang H W,Ren Y F,Wen R Z,Xu P B. 2019. Breakdown of earthquake self-similar scaling and source rupture directivity in the 2016−2017 central Italy seismic sequence[J]. J Geophys Res:Solid Earth,124(4):3898–3917. doi: 10.1029/2018JB016543
Wang H W,Wen R Z. 2021. Attenuation and basin amplification revealed by the dense ground motions of the 12 July 2020 MS 5.1 Tangshan,China,earthquake[J]. Seismol Res Lett,92(4):2109–2121. doi: 10.1785/0220200400
Wen R Z,Wang H W,Ren Y F. 2015. Rupture directivity from strong-motion recordings of the 2013 Lushan aftershocks[J]. Bull Seismol Soc Am,105(6):3068–3082. doi: 10.1785/0120150100
Yang T,Li B R,Fang L H,Su Y J,Zhong Y S,Yang J Q,Qin M,Xu Y J. 2022. Relocation of the foreshocks and aftershocks of the 2021 MS6.4 Yangbi earthquake sequence,Yunnan,China[J]. J Earth Sci,33(4):892–900. doi: 10.1007/s12583-021-1527-7
Yenier E,Atkinson G M. 2015. Regionally adjustable generic ground-motion prediction equation based on equivalent point-source simulations:Application to central and eastern North America[J]. Bull Seismol Soc Am,105(4):1989–2009. doi: 10.1785/0120140332
Yoshida K. 2019. Prevalence of asymmetrical rupture in small earthquakes and its effect on the estimation of stress drop:A systematic investigation in inland Japan[J]. Geosci Lett,6(1):16–23. doi: 10.1186/s40562-019-0145-z
Yoshida K,Saito T,Emoto K,Urata Y,Sato D. 2019. Rupture directivity,stress drop,and hypocenter migration of small earthquakes in the Yamagata-Fukushima border swarm triggered by upward pore-pressure migration after the 2011 Tohoku-Oki earthquake[J]. Tectonophysics,769(2019):228184.
Zhou Y J,Ren C M,Ghosh A,Meng H R,Fang L H,Yue H,Zhou S Y,Su Y J. 2022. Seismological characterization of the 2021 Yangbi foreshock-mainshock sequence,Yunnan,China:More than a triggered cascade[J]. J Geophys Res:Solid Earth,127(8):e2022JB024534. doi: 10.1029/2022JB024534
-
期刊类型引用(18)
1. 张海洋，毛丰龙，王朝景，刘德强，李博. 2020年唐山古冶5.1级地震前后岩石圈磁场异常变化分析. 地震工程学报. 2024(01): 224-231+240 . 
百度学术
2. 赵慧琴. 南北地震带历史地震与岩石圈磁场空间分布特征分析. 地震地磁观测与研究. 2024(01): 56-62 . 百度学术
3. 张海洋. 2019年丰南4.5级地震前后岩石圈磁场水平矢量动态演化特征. 中国地震. 2024(02): 426-435 . 百度学术
4. 王朝景，李博，苏树朋. 基于多期累积岩石圈磁场变化分析唐山M_S5.1地震震磁异常. 地震研究. 2024(04): 517-527 . 百度学术
5. 梁霄，储飞，徐如刚，孙鸿博，肖伟鹏，王俊. 2014年霍山M_S4.3地震前后重磁场变化特征及机理分析. 地震地质. 2024(05): 1151-1171 . 百度学术
6. 张海洋，李博，刘德强，王朝景，赵慧琴，毛丰龙，苏树朋. 河北及邻区震前岩石圈磁场异常变化分析. 地震工程学报. 2023(01): 191-198 . 百度学术
7. 刘德强. 长清M4.1地震岩石圈磁场异常特征. 大地测量与地球动力学. 2023(04): 425-428 . 百度学术
8. 刘德强，苏树朋，赵慧琴，昝书林. 山西代县M_S3.0地震地磁异常与构造的关系. 华北地震科学. 2023(01): 66-68+74 . 百度学术
9. 刘德强. 华北中部主要活动断裂对岩石圈磁异常分布的影响. 大地测量与地球动力学. 2022(05): 499-504 . 百度学术
10. 刘德强，云萌，王朝景，李博，张海洋，毛丰龙. 有关小震震前数天震中附近岩石圈磁场Z分量异常特征. 华北地震科学. 2022(02): 77-81 . 百度学术
11. 张海洋，苏树朋，赵慧琴. 2022年青海门源6.9级地震前岩石圈磁场异常变化分析. 地震工程学报. 2022(03): 735-743 . 百度学术
12. 刘德强，王朝景，李博，张海洋，毛丰龙. 地磁台站连续观测分钟值连接阈值选择. 地震地磁观测与研究. 2022(02): 63-70 . 百度学术
13. 张海洋，李博，苏树朋，王朝景，刘德强，赵慧琴，毛丰龙. 张家口—渤海地震带岩石圈磁场时空变化特征. 地震. 2022(01): 133-144 . 百度学术
14. 王朝景，李博，刘德强，张海洋，毛丰龙，苏树朋. 河北及周边地区岩石圈磁场分布与震磁特征. 华北地震科学. 2022(04): 83-87+94 . 百度学术
15. 王朝景，刘德强，苏树朋，李博，张海洋，毛丰龙，赵慧琴. 2018年永清M_(W)4.3地震前后岩石圈磁场变化特征. 地震地磁观测与研究. 2021(03): 129-136 . 百度学术
16. 刘德强，王朝景，李博，张海洋，云萌，毛丰龙. 河北唐山M_S5.1地震岩石圈磁场D分量异常分布和应力的关系. 华北地震科学. 2021(04): 75-79 . 百度学术
17. 刘德强，王朝景，苏树朋，李博，赵慧琴，张海洋，毛丰龙. 华北中部岩石圈磁场变化与岩石圈结构间的关系. 地震地磁观测与研究. 2021(05): 97-102 . 百度学术
18. 王朝景，李博，苏树朋，张海洋，刘德强，赵慧琴，毛丰龙. 2018年冀南地区岩石圈磁场变化分析. 高原地震. 2020(03): 9-12+29 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 205
- HTML全文浏览量: 63
- PDF下载量: 61
- 被引次数: 19
