２０２０年新疆伽师MS６．４地震震源区微震检测与目录完备性分析

第４５卷㊀第６期２０２３年１１月
地㊀震㊀工㊀程㊀学㊀报
C H I N A E A R T H Q U A K EE N G I N E E R I N GJ O U R N A L
V o l ．４５㊀N o ．６
N o v e m b e r ,２０２３
㊀㊀收稿日期:２０２１Ｇ１１Ｇ２２
㊀㊀基金项目:地震科技星火计划项目(X H ２１３７０３Y )
㊀㊀第一作者简介:梁姗姗,高级工程师,主要从事地震定位和震源机制反演等研究工作.E Ｇm a i l :l i a n g s h a n s h a n ＠s e i s ．a c ．c n .㊀㊀通信作者:张广伟,副研究员,主要从事地震定位和震源机制反演研究.E Ｇm a i l :z h a n g g
w＠m a i l ．u s t c ．e d c ．c n .梁姗姗,张广伟,邹立晔,等．２０２０年新疆伽师M S ６．４地震震源区微震检测与目录完备性分析[J ]．地震工程学报,２０２３,４５(６)
:１４８５Ｇ１４９２．D O I :１０．２００００/j
．１０００Ｇ０８４４．２０２１１１２２００１L I A N GS h a n s h a n ,Z HA N GG u a n g w e i ,Z O UL i y e ,e t a l ．M i c r o s e i s m i c d e t e c t i o n a n d c a t a l o g u e c o m p l e t e n e s s a n a l y
s i s o f t h e f o c a l a r e a o f t h e ２０２０J i a s h i M S ６．４e a r t h q u a k e [J ]．C h i n aE a r t h q u a k eE n g i n e e r i n g J o u r n a l ,２０２３,４５(６):１４８５Ｇ１４９２．D O I :１０．２００００/j ．１０００Ｇ０８４４．２０２１１１２２００１
２０２０年新疆伽师M S ６．４地震震源区
微震检测与目录完备性分析
梁姗姗１,张广伟２,邹立晔１,刘艳琼１,姬运达１
(１．中国地震台网中心,北京１０００４５;
２．应急管理部国家自然灾害防治研究院,北京１０００８５
)摘要:中强地震发生后,
地震检测因受到尾波的干扰可能会遗漏部分微震事件,影响地震目录的完备性.文章利用波形模板匹配方法对２０２０年新疆伽师M S ６．４地震序列开展微震检测,相比原始的中国地震台网中心统一地震目录,新检测出１７５６个微震事件,地震数量增加了１．３倍.基于检测后的余震目录计算最小完备震级为M L １．２,地震活动性b 值为０．７６,较原始目录的M L １．６和０．７７均有所降低.通过伽师震源区地震序列活动特征分析,结果表明前震序列在主震前短时间内(前３６小时)
出现地震活动的密集增强,相应的b 值显示为低值;主震发生后地震序列完备震级较高,随着时间的推移,完备震级缓慢降低并趋于稳定,并且呈周期性的波动.本研究提高了伽师震源区地震目录的完备性,为精细化描述该地区地震序列时空演化特征提供了关键数据基础.
关键词:伽师M S ６．４震源区;微震检测;最小完备震级;地震活动性中图分类号:P ３１５．２㊀㊀㊀㊀㊀㊀文献标志码:A㊀㊀㊀文章编号:１０００Ｇ０８４４(２０２３)０６－１４８５－０８
D O I :１０．２００００/j
．１０００Ｇ０８４４．２０２１１１２２００１M i c r o s e i s m i c d e t e c t i o na n d c a t a l o g u e c o m p l e t e n e s s a n a l y
s i s o f t h e f o c a l a r e a o f t h e ２０２０J i a s h i M S ６．４e a r t h q
u a k e L I A N GS h a n s h a n １,Z H A N G G u a n g w e i ２,Z O U L i y e １,L I U Y a n q i o n g １,
J IY u n d a １
(１．C h i n aE a r t h q u a k eN e t w o r k sC e n t e r ,B e i j i n g １
０００４５,C h i n a ;２．N a t i o n a l I n s t i t u t e o f N a t u r a lH a z a r d s ,M i n i s t r y o f E m e r g e n c y M a n a g e m e n t ,B e i j i n g １
０００８５,C h i n a )A b s t r a c t :M a n y m i c r o Ｇe v e n t sm a y b em i s s e d i n e a r t h q
u a k e d e t e c t i o n d u e t o t h e i n t e r f e r e n c e o f c o Ｇd aw a v e s a f t e rm o d e r a t e Ｇs t r o n g e a r t h q u a k e s ,t h u s a f f e c t i n g t h e c o m p l e t e n e s s o f e a r t h q
u a k e c a t a Ｇl o g u e ．I n t h i s p a p e r ,m i c r o s e i s m i c d e t e c t i o no f t h e J i a s h i M S ６．４e a r t h q u a k e s e q
u e n c ew a s c a r r i e d o u t u s i n g t h e m a t c h e df i l t e r m e t h o d ．１７５６m i c r o s e i s m i ce v e n t sw e r ed e t e c t e d ,w h i c h w a s１．３t i m e sm o r e t h a n t h e c a t a l o g u e f r o m C h i n aE a r t h q
u a k eN e t w o r k sC e n t e r (C E N C )．T h em i n i m u m
m a g n i t u d e o f c o m p l e t e n e s s o f t h e d e t e c t e da f t e r s h o c kc a t a l o g
u e i s M L １．２a n d t h e b v a l u e o f s e i s Ｇm i c i t y i s ０．７６,w h i c ha r e l o w e r t h a nt h o s eo f t h eo r i g i n a l c a t a l o g
u e (M L １．６a n d b ＝０．７７)．T h e c h a r a c t e r i s t i c s o f s e i s m i cs e q u e n c ea c t i v i t i e s i nJ i a s h i f o c a l a r e a w e r ea n a l y
z e d ,a n dt h er e s u l t s s h o wt h a t t h e f o r e s h o c k s e q u e n c e i n t e n s i v e l y s t r e n g t h e n sw i t h i n a s h o r t t e r m (t h e f i r s t ３６h o u r s )b e f o r e t h em a i n s h o c k a n d t h e c o r r e s p o n d i n g b v a l u e i s l o w．T h em a g n i t u d e o f c o m p
l e t e n e s s o f t h e s e i s m i c s e q u e n c e i s h i g h a f t e r t h em a i n s h o c k ,t h e n i t s l o w l y d
e c r e a s e s a n d t e n d s t o b e s t a b l ew i t h p e r i o d i c
f l u c t u a t i o n s ．T h e s t u d y i m p r o v e s t h e c o m p l e t e n e s so f t h e e a r t h q u a k e c a t a l o g
u e i nJ i a s h i f o c a l a r e a ,t h u s p r o v i d i n g k e y d a t a f o r a d e t a i l e d d e s c r i p t i o n o f t h e s p
a c e Ｇt i m e e v o l u t i o n c h a r a c t e r Ｇi s t i c s o f t h e e a r t h q u a k e s e q u e n c e i n t h e r e g
i o n ．K e y
w o r d s :f o c a l a r e a o f J i a s h i M S ６．４e a r t h q u a k e ;m i c r o s e i s m i c d e t e c t i o n ;m i n i m u m m a g n i t u d e o f c o m p l e t e n e s s ;s e i s m i c i t y
０㊀引言
地震目录是研究地震活动性㊁地震发震过程㊁地震预测以及地震危险性分析的重要基础资料.地震目录的完备性直接关系到研究结果的科学性和可靠性.已有研究表明,强震发生后,由于地震波形受主震和强余震尾波的干扰,使波形互相叠加,导致很难清晰识别叠加地震事件的震相,因此主震后短时间会遗漏相当比例的余震.近年来,随着模板匹配技术的发展和应用,这一问题得到了很好地解决[
１Ｇ６
].２０２０年１月１９日新疆喀什地区伽师县发生
M S ６．４地震(３９．８４ʎN ,７７．２１ʎE ).M S ６．４地震发生后,震中附近又相继发生１次M S ５．２余震和多次
M S ４．０以上余震.此前一天,即１月１８日０时０５
分,距离M S ６．４地震震中不远处曾发生１次M S ５．４前震.此次伽师M S ６．４地震序列致使１人死亡㊁２人轻伤,４０００余间房屋受损,部分水库㊁道路㊁桥梁等基础设施遭到破坏,直接经济损失达１６．２亿元
(h t t p
s ://b a i k e ．b a i d u ．c o m /).伽师M S ６．４地震位于天山西南麓柯坪塔格推覆构造最前缘的柯坪断裂附近,余震序列位于接近平行的E W 向柯坪断裂和奥兹格尔他乌断裂之间.由于受印度板块向欧亚板块俯冲碰撞作用的影响,伽师震源区附近历史地震活动频繁,１９９６年以来曾发生多次M S ６．０以上地震(图１),其中最大的１次为１９９６年３月１９日的M S ６．９地震
.
五角星表示M S ６．０以上地震(１９７０年１月１日至今),圆圈表示模板事件(２０２０年１月１７日 ２０２１年１月３１日)
图１㊀研究区域内M S ȡ６．０历史地震和模板事件分布
F i g ．１㊀D i s t r i b u t i o no f t h eh i s t o r i c a l a n d t e m p
l a t e e v e n t s o f M S ȡ６．０i n t h e r e s e a r c ha r e a ６８４１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年
㊀㊀据中国地震台网中心统一地震编目显示(２０２０年１月１７日０时至１月３１日２４时,北京时),此次地震序列共记录到地震１４７７次,含１７１次单台记录结果.由于伽师震源区附近地震台站较为稀疏,特别是震中北部和南部地区台站空白,造成台站最大孔隙角较大(图１),因此很难应用传统的地震定位方法来确定更多小地震的目录信息.因此,为获得伽师震源区更为完整的地震事件资料,本文选取震中１２０k m范围内４个台站,利用２０２０年１月１７日至１月３１日(北京时)中国地震台网记录的伽师震源区波形资料,以初始地震定位目录中信噪比较高的７３３个地震作为模板事件,采用匹配定位方法,检测微小震事件[７],并补充到初始地震目录中,提高地震序列的完备性,进而利用M C值和b值,对此次地震序列的时间演化特征加以分析研究.
１㊀微地震检测
本研究使用２０２０年１月１７日至２０２０年１月３１日(北京时)新疆区域台网的波形资料,选择震中距在１２０k m以内的共４个台站的波形(八盘水磨B P M,西克尔X K R,喀什K S H,阿图什A T S),部分台站分布如图１.选用全国统一地震编目的目录中相同时间段内M Lȡ１．０㊁信噪比较高的共７３３个地震作为模板事件(图１).
本研究采用匹配定位方法(M a t c ha n dL o c a t e,简称M&L)[８]对常规产出的地震目录中所遗漏的地震事件进行微地震检测.这一方法与波形模板匹配方法相类似[９],二者均利用模板事件波形与可能的微震信号做互相关叠加,从而检测出微震事件.因为在叠加之前需要对模板事件周围三维空间进行扫描搜索,来计算同一台站上可能存在的事件位置与参考位置之间的相对走时差,并据此对互相关波形进行矫正叠加,所以相较于当今广泛流行的微震检测方法,M&L法具有以下优点:(１)可检测到事件的震级下限更小;(２)可以检测到距离模板事件较远的微事件;(３)对速度模型的依赖性较小;(４)可以给出高精度微事件位置信息.在采用匹配定位方法进行遗漏地震事件检测过程中,根据叠加后相关波形的平均相关系数(C C)㊁信噪比(S N R)两部分来确定检测所需的标准阈值;当检测到微震事件后,则根据台站记录中参考震相的振幅来确定震级.
微震检测之前,我们需要对连续波形进行预处理,包括去均值㊁去线性趋势等.微震检测时,地震波的滤波频段设为２~８H z,将长周期噪声干扰去除.在三维网格空间内,以参考模板事件的震源位置作为搜索中心,对微震的位置进行搜索,其经度方向和纬度方向的搜索范围为０．０５ʎ,搜索间隔为０．０１ʎ;深度方向的搜索范围为２k m,搜索间距为０．１k m.计算模板事件与可能的微震事件之间的走时差,然后利用该走时差对互相关波形进行走时校正并叠加.模板事件波形与连续波形之间的滑动互相关以S g震相作为参考震相.地震模板匹配中,计算S g的理论到时和慢度参数所用一维区域速度模型,与崔仁胜等[１０]对２０２０年伽师地震序列进行地震定位和震源机制求解的相同;模板波形的互相关窗口为S g波理论到时的前１s和后３s.本研究在微震检测过程中,将互相关系数C C１㊁C C２和波形信噪比经验性地设置为０．３５㊁０．３和１０倍背景相关系数[１１].
通过模板检测计算和人工复核,伽师震源区２０２０年１月１７日 ２０２０年１月３１日时间段内,共检测到３６２２个地震事件(图２).其中包括７３３个模板事件的自检测结果,以及与中国地震台网中心统一编目已定位目录中相同的地震事件１１３３个和新检测地震事件１７５６个,所检测出来的地震事件个数为传统方法已定位地震个数的１．３倍,新检测出来的地震事件个数大于已定位地震事件的个数.且新检测地震目录比原始地震目录在数量上有较大的提升[图２(a)],原始地震目录的震级范围在M L１．０~２．０之间,而新地震目录震级范围在M L０．５~１．５之间,补充了更多小震级事件.原始地震目录通常存在小震级地震事件缺失,其原因在于背景噪声对地震波形的干扰,以及大地震后能量快速释放而引起的后续地震事件波形的互相叠加.从时间分布上可以看出[图２(b)],主震后１天内遗漏事件较多,遗漏频次随着时间推移,逐渐趋于稳定.需要说明的是,尽管有原始目录中１５８个地震事件没有检测出来,这些事件主要集中在主震后２天内,且震级较小,主要分布在M L０．５~１．５之间,产生的原因可能是:原始目录中部分事件采用２个或者３个台站的数据进行计算,而本研究为了获取准确度和精度更高的震源参数位置,微震检测采用４个台站的数据进行定位,因此少于４个台站记录的微震事件会被忽略.
７８４１
第４５卷第６期㊀㊀㊀㊀㊀梁姗姗,等:２０２０年新疆伽师M S６．４地震震源区微震检测与目录完备性分析㊀㊀㊀㊀㊀
图２㊀不同目录震级Ｇ频次和时间Ｇ频次分布
F i g．２㊀D i a g r a mo fm a g n i t u d eＧf r e q u e n c y a n d t i m eＧf r e q u e n c y o f d i f f e r e n t c a t a l o g u e s
㊀㊀图３展示了微震事件检测实例.模板为２０２０年１月１８日０时１７分的M L２．０地震事件,而检测到的事件对应为２０２０年１月１８日０时５２分的M L１．３事件.地震模板事件波形与M&L所检测到的微震信号的叠加平均相关系数为０．５７３１;
由于归一化互相关确定的互相关值,当两个地震事件的震源时间函数和辐射类型一致时,其互相关值为１.但是,对于实际的地震波形,
即使上述两
红色波形为模板地震事件波形,地震时间为２０２０Ｇ０１Ｇ１８,００:１７:１９;灰色波形为用于识别微震事件的连续
波形,新的微震事件识别结果,其时间为２０２０Ｇ０１Ｇ１８,００:５２:４３,震级为M L１．３;左侧大写字母表示台站代
码和分量;右侧数字表示互相关系数
图３㊀微震识别模板匹配扫描结果实例
F i g．３㊀E x a m p l e o fm i c r o s e i s m i c r e c o g n i t i o n r e s u l tw i t h t h em a t c h e d f i l t e rm e t h o d
者一致,噪声的出现也会使得互相关值小于１.模
板事件波形之间的平均背景互相关值通常在０．０３
左右,当采用９倍的平均背景互相关值作为地震检
测的阈值,其C C值在０．２左右,图３展示的本研究
检测结果C C值为０．５７３１,表明其相关性较高,成
功检测到１次微震事件.
８８４１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年
我们将本研究检测到的遗漏地震,以及自检测地震事件和原始目录中没有检测出来的地震事件组合在一起,得到了更加完整㊁丰富的伽师地震序列目录(图４).从检测事件的空间分布上可以看出,原始目录在平面上分布较分散,很难明显分辨地震的展布方向,而本研究所得目录的地震分布更加集中,展布形态更加清晰.由图４可知,２０２０年伽师地震序列主要呈两个优势方向展布,分别为NWW向和S S W向,其中大多数余震分布在NWW向,与S S W 向余震呈不对称的共轭特征分布.此外,主要余
震
图４㊀原始目录和本研究所得目录空间分布
F i g．４㊀S p a t i a l d i s t r i b u t i o no f t h e e v e n t s f r o mt h e o r i g i n a l c a t a l o g u e a n d t h e o n e o b t a i n e d i n t h i s s t u d y
序列以S向零散分布一些余震.
图５给出了本研究所得目录的地震序列的震
级Ｇ时间和震源深度分布.伽师地震序列MＧt图显
示[图５(a)],前震序列和余震序列均呈现衰减趋
势,无论是地震的震级还是发生频次都在相应地减
小,但前震序列活动频度远小于余震序列活动频度
.
图５㊀地震序列震级及深度演化图
F i g．５㊀M a g n i t u d e a n dd e p t he v o l u t i o nm a p o f e a r t h q u a k e s e q u e n c e
９８４１
第４５卷第６期㊀㊀㊀㊀㊀梁姗姗,等:２０２０年新疆伽师M S６．４地震震源区微震检测与目录完备性分析㊀㊀㊀㊀㊀
１月１８日随着０时５分M L ５．６地震的发生,
震源区地震活动明显增强,在１９日下午６时,主震发生前３小时出现短暂的地震活动平静,２０２１年云南漾濞
M S ６．４地震的前震序列也出现类似的特征[１２]
,该特征为前震序列的典型特征,可能是深部大断裂发生破裂之前快速蠕动的临震表现.从震源深度分布可以看出[图５(b )
],前震深度和余震深度均在约５~１０k m 和１５~２０k m 处存在明显的优势发震层.２㊀地震目录完备性和序列活动特征分析
完备震级M C (M a g n i t u d eo fC o m p
l e t e n e s s )是评估地震台网监测能力的一种方法,并且通过确定完备震级,使用该完备震级以上的地震目录进行地震序列活动性特征讨论时,可以确保研究结果的合
理性和可靠性[
１３Ｇ１４
].地震频次随震级变化而变化,并呈曲线分布,其频次峰值所对应的震级即为M C
值,代表了目录的完备性.若M C 值估算得太高,则有可能导致未能充分使用地震数据,而未能获得完整的地震信息;但若M C 估算值比真实值低,则又有可能使得后续分析所用数据并不完整,导致模型参
数值计算的错误结果,进而得出错误结论[
１５
].本文使用最大曲率法MA X C (M a x i m u m C u r v a t u r e
T e c h n i q
u e )对比分析了模板匹配前后地震目录的最小完整震级M C ;作为一种常用的快速㊁直接㊁稳健
的方法,该方法不需要进行任何参数拟合.同时,该方法也能够在地震目录数据相关性较小时,计算出合理的M C 值.以往研究表明,b 值反映地震活动和区域应力变化,b 值的升降变化不仅直接反映了区域内不同震级地震的比例关系,而且体现了区域内地下介质的应力状态变化.一般而言,应力水平升高,则b 值降低;b 值升高,
则应力水平降低.本研究采用Z M A P 程序包[１
６]
中的最大似然法[１７
]
计算b 值.为了对匹配定位检测所得地震目录的质量加以
评估,我们分别分析了新旧目录各自的完备性.检测后地震目录的完备震级由台网中心统一地震目录的M L １．６,降至检测目录的M L １．２(图６).这再次表明:匹配定位方法可有效提高微小地震的检测能
力,并降低完备震级M C .继而以各自的最小完整
震级为起始震级,分别计算其b 值,从而获得检测前后的伽师震源区的b 值.检测前的b 值为０．７７,误差为ʃ０．０４;而检测后的b 值则为０．７６,
误差为ʃ０．０２.比较遗漏事件地震补充前后b 值结果,b 值数值与误差均降低,这是由于检测后b 值的计算地震完备度较高
.
图６㊀匹配定位检测前后最小完备震级以及b 值对比
F i g ．６㊀C o m p a r i s o nb e t w e e n t h em i n i m u m m a g n i t u d e o f c o m p
l e t e n e s s a n d t h e b v a l u e sb e f o r e a n da f t e r d e t e c t i o n ㊀㊀基于本研究所得最新目录,
我们以天为单位,对２０２０年１月１７日 １月３１日研究时间段的完备震级和b 值进行统计,以分析伽师震源区的完备震级
M C 和b 值的时变演化特征(
图７),其中每个时窗内地震事件的最小样本数为５０个.如图７(a )所示,主震后的短时间内,完备震级较高,最高值达M L ２．２;
这可能是因为主震尾波干扰对后续波形造成了严重影响.随着时间的推移,完备震级慢慢变小,并逐渐趋于稳定,并且呈周期性的波动;产生这种现象的原因可能是主震波形能量逐渐衰减,后续波形受主震影响较小,M C 周期性波动与白天㊁夜晚噪声水平不同有关.由于大部分地震都发生在２０k m 或更浅,
在分析b 值随时间变化演化特征时,为了分清b 值的变化来源,我们参照刘雁冰等[１８]
分
析汶川地震b 值随时间变化的方法来分析演化特征.因此,计算b 值随时间变化只考虑２０k m 或更浅的区域,之后同样以天为单位,分析余震序列b 值
随时间演化的特征[图７(b )].在时间段内地震序０９４１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年
列发生过程中,b 值总体低于１．２,变化范围在０．６~
１．２之间,M S ６．４主震前b 值约为０．６左右,
存在低值的现象,可能显示应力不断积累的过程,符合大地震前的变化规律;主震发生后,b 值起伏明显,可能与震后余震区断层上较为剧烈的应力调整有关[
１９]
.当b 值降低时可能反映了局部的差应力增大,因为很多余震都在震后时间段内发生,而b 值的缓慢恢复,则反映了差应力逐渐下降,并低于震前水平.这种b 值随时间的变化趋势可能反映了伽师震源区应力积累和释放的演化过程
.
图７㊀地震序列完备震级和b 值演化图
F i g ．７㊀M a g n i t u d e o f c o m p
l e t e n e s s a n d b Ｇv a l u e e v o l u t i o n m a p o f s e i s m i c s e q
u e n c e ３㊀结论
本文通过模板匹配方法,对２０２０年１月１９日
新疆伽师M S ６．４地震后１５天的连续波形数据,开展遗漏地震检测,共检测出１７５６个新地震事件,约为中国地震台网中心目录给出的已定位地震数量１３０６的１．３倍.并将检测出来的地震目录完备震级
进行对比,使完备震级M C 由M L １．６减小到M L １．２,地震活动性b 值由检测前的０．７７ʃ０．０４减小到
０．７６ʃ０．０２.
利用本研究所得到地震目录计算得到的伽师震源区的完备震级和b 值随时间变化显示,震源区的完备震级在主震后迅速上升,最高值达M L ２．２,然后降低并逐渐稳定在M L １．０左右;主震发生前地震序列的b 值较低,表明此时震源区应力状态处于较高水平;而主震发生后的b 值随时间的推移而迅速
增大.
本研究所获得的地震目录为２０２０年伽师地震
序列增加了更多的样本,为精细化描述该地区地震序列时空演化特征提供了关键数据基础,为该区域未来地震活动性判断提供了科学的参考资料.另外,在更加完整的目录前提下进一步计算出来的b 值可以为后续分析伽师地区地震活动性以及强震危险性都具有现实的意义.本研究表明,模板匹配方法可有效提高遗漏地震的检测率及地震目录的完备性;同时,模板检测自动化的计算模式节省了大量的人为工作量,在中强地震发生后能够高效而快速地识别出微震事件,为震后应急决策救援提供关键的地震序列目录数据.参考文献(R e f e r e n c e s
)[１]㊀侯金欣,王宝善．２０１４年鲁甸M S ６．５地震前后地震活动性[J ]．
地球物理学报,２０１７,６０(４):１４４６Ｇ１４５６．
H O UJ i n x i n ,WA N GB a o s h a n ．T e m p o r a l e v o l u t i o no f s e i s m i c i Ｇt y b e f o r e a n d a f t e r t h e ２０１４L u d i a n M S ６．５e a r t h q u a k e [J ]．C h i Ｇn e s e J o u r n a l o fG e o p h y
s i c s ,２０１７,６０(４):１４４６Ｇ１４５６．[２]㊀谭毅培,邓莉,曹井泉,等．２０１５年河北滦县震群发震机理分析
[J ]．地球物理学报,２０１６,５９(１１):４１１３Ｇ４１２５．T A N Y i p e i ,D E N G L i ,C A O J i n g q u a n ,e t a l ．S e i s m o l o g i c a l m e c h a n i s ma n a l y s i so f ２０１５L u a n x i a ns w a r m ,H e b e iP r o v i n c e [J ]．C h i n e s e J o u r n a l o fG e o p h y s i c s ,２０１６,５９(１１):４１１３Ｇ４１２５．[３]㊀尹欣欣,杨立明,赵林林,等．九寨沟M ７．０级地震余震目录完
备性研究[J ]．地球物理学进展,２０２０,３５(２):４７５Ｇ４７９．Y I NX i n x i n ,Y A N GL i m i n g ,Z HA OL i n l i n ,e t a l ．M i s s i n g e
a r t h Ｇq u a k e s d e t e c t i o na n d c o m p l e t e n e s s o f e a r t h q u a k e c a t a l o g
u e s o f t h e ２０１７J i u z h a i g o u M ７．０e a r t h q u a k e [J ]．P r o g r e s s i nG e o p h y
s Ｇi c s ,２０２０,３５(２):４７５Ｇ４７９．
[４]㊀王鹏,侯金欣,吴朋．２０１７年九寨沟M S ７．０地震序列活动特征
分析[J ]．中国地震,２０１７,３３(４):４５３Ｇ４６２．
WA N GP e n g ,H O U J i n x i n ,WU P e n g ．T e m p
o r a le v o l u t i o no f t h e s e i s m i c i t y o ft h e２０１７J i u z h a i g o u M S ７．０e a r t h q u a k es e Ｇq u e n c e [J ]．E a r t h q
u a k eR e s e a r c hi n C h i n a ,２０１７,３３(４):４５３Ｇ４６２．
[５]㊀胡涛,陈继锋,杨立明,等．青海茫崖５．０级地震震源区微震检
测与序列活动特征分析[J ]．地震工程学报,２０２０,４２(４):９４１Ｇ
９４７．
HU T a o ,C H E N J i f e n g ,Y A N G L i m i n g
,e ta l ．D e t e c t i o n o f m i s s i n g e a r t h q u a k e s a n d a c t i v i t y c h a r a c t e r i s t i c s o f t h eM a n g y
a M ５．０e a r t h q u a k es e q u e n c ei n Q i n g h a i [J ]．C h i n a E a r t h q
u a k e E n g i n e e r i n g J
o u r n a l ,２０２０,４２(４):９４１Ｇ９４７．[６]㊀徐志国,梁姗姗,刘敬光,等．２０２０年云南巧家M S ５．０地震微震
检测及发震构造初步探讨[J ]．中国地震,２０２１,３７(２):３９０Ｇ３９９．
X UZ h i g u o ,L I A N GS h a n s h a n ,L I UJ i n g g u a n g
,e t a l ．M i c r o s e i s Ｇm i c d e t e c t i o na n d s e i s m o g e n i c s t r u c t u r e a n a l y s i s o f t h eQ i a o j
i a １９４１第４５卷第６期㊀㊀㊀㊀㊀梁姗姗,等:２０２０年新疆伽师M S ６．４地震震源区微震检测与目录完备性分析㊀㊀㊀㊀㊀
M S５．０e a r t h q u a k e i nY u n n a nP r o v i n c e i n２０２０[J]．E a r t h q u a k e R e s e a r c h i nC h i n a,２０２１,３７(２):３９０Ｇ３９９．
[７]㊀Z H A N G M,W E NLX．A n e f f e c t i v em e t h o d f o r s m a l l e v e n t d eＧt e c t i o n:m a t c ha n d l o c a t e(M&L)[J]．G e o p h y s i c a l J o u r n a l I nＧt e r n a t i o n a l,２０１５,２００(３):１５２３Ｇ１５３７．
[８]㊀张淼．地震定位和检测[D]．合肥:中国科学技术大学,２０１５．Z H A N G M i a o．E a r t h q u a k el o c a t i o na n dd e t e c t i o n[D]．H e f e i: U n i v e r s i t y o f S c i e n c e a n dT e c h n o l o g y o fC h i n a,２０１５．[９]㊀P E N GZG,Z H A OP．M i g r a t i o no f e a r l y a f t e r s h o c k s f o l l o w i n g t h e２００４P a r k f i e l de a r t h q u a k e[J]．N a t u r eG e o s c i e n c e,２００９,２(１２):８７７Ｇ８８１．
[１０]㊀崔仁胜,赵翠萍,周连庆,等．２０２０年１月１９日新疆伽师M S６．４地震序列的活动特征和发震构造[J]．地震地质,２０２１,
４３(２):３２９Ｇ３４４．
C U IR e n s h e n g,Z H A O C u i p i n g,Z HO U L i a n q i n g,e ta l．S e i sＧ
m i c i t y f e a t u r e a n d s e i s m o g e n i c f a u l t o f t h e M S６．４e a r t h q u a k e
s e q u e n c e o n J a n u a r y１９,２０２０i n J i a s h i,X i n j i a n g[J]．S e i s m o l oＧ
g y a n dG e o l o g y,２０２１,４３(２):３２９Ｇ３４４．
[１１]㊀郭铁龙,高原．剪切波分裂揭示的青藏高原上地壳地震各向异性基本特征[J]．地球物理学报,２０２０,６３(３):１０８５Ｇ１１０３．
G U O T i e l o n g,G A O Y u a n．S e i s m i ca n i s o t r o p y i nt h eu p p e r
c r u s tw i t h i nT i b e t a nP l a t e a u r e v e a l e db y s h e a rＧw a v e s p l i t t i n g
[J]．C h i n e s e J o u r n a l o fG e o p h y s i c s,２０２０,６３(３):１０８５Ｇ１１０３．[１２]㊀段梦乔,赵翠萍,周连庆,等．２０２１年５月２１日云南漾濞M S６．４地震序列发震构造[J]．地球物理学报,２０２１,６４(９):
３１１１Ｇ３１２５．
D U A N M e n g q i a o,Z H A O C u i p i n g,Z HO U L i a n q i n g,e ta l．
S e i s m o g e n i cs t r u c t u r eo ft h e２１M a y２０２１M S６．４Y u n n a n
Y a n g b i e a r t h q u a k e s e q u e n c e[J]．C h i n e s e J o u r n a l o fG e o p h y sＧ
i c s,２０２１,６４(９):３１１１Ｇ３１２５．
[１３]㊀张广伟．云南地区地震的重新定位及b值研究[J]．中国地震,２０１６,３２(１):５４Ｇ６２．
Z HA N G G u a n g w e i．R e s e a r c ho ne a r t h q u a k e r e l o c a t i o na n d bＧ
v a l u e i n t h eY u n n a na r e a[J]．E a r t h q u a k eR e s e a r c h i nC h i n a,
２０１６,３２(１):５４Ｇ６２．
[１４]㊀徐志国,梁姗姗,盛书中,等．２０１９年四川长宁M S６．０地震序列重定位和震源特征分析[J]．地震学报,２０２０,４２(４):３７７Ｇ
３９１,５０９．
X UZ h i g u o,L I A N GS h a n s h a n,S H E N GS h u z h o n g,e t a l．R e l oＧ
c a t i o n a n
d s o u r c
e c h a r a c t e r i s t i c s o ft h e２０１９C h a n g n i n g
M S６．０e a r t h q u a k es e q u e n c e[J]．A c t a S e i s m o l o g i c a S i n i c a,２０２０,４２(４):３７７Ｇ３９１,５０９．
[１５]㊀解孟雨,孟令媛．川东南长宁地区地震活动及b值演化特征[J]．中国地震,２０２１,３７(２):４９４Ｇ５０７．
X I E M e n g y u,M E N GL i n g y u a n．S e i s m i c i t y a n d e v o l u t i o n c h a rＧ
a c t e r i s t i c s o f bＧv a l u e so fC h a n g n i n g a r e a i ns o u t h e a s t e r nr eＧ
g i o no fS i c h u a n B a s i n[J]．E a r t h q u a k e R e s e a r c hi n C h i n a,
２０２１,３７(２):４９４Ｇ５０７．
[１６]㊀W I E M E RS．As o f t w a r e p a c k a g e t o a n a l y z e s e i s m i c i t y:Z MA P [J]．S e i s m o l o g i c a lR e s e a r c hL e t t e r s,２００１,７２(３):３７３Ｇ３８２．[１７]㊀S H IY L,B O L T B A．T h es t a n d a r de r r o ro f t h em a g n i t u d eＧ
f r e q u e n c y b v a l u e[J]．B u l l e t i no f t h eS e i s m o l o
g i c a l S o c i e t y o f
A m e r i c a,１９８２,７２(５):１６７７Ｇ１６８７．
[１８]㊀刘雁冰,裴顺平．汶川地震前后b值的时空变化及构造意义[J]．地球物理学报,２０１７,６０(６):２１０４Ｇ２１１２．
L I U Y a n b i n g,P E I S h u n p i n g．T e m p o r a l a n d s p a t i a l v a r i a t i o n o f
bＧv a l u e b e f o r e a n d a f t e rW e n c h u a n e a r t h q u a k e a n d i t s t e c t o n i c
i m p l i c a t i o n[J]．C h i n e s eJ o u r n a lo fG e o p h y s i c s,２０１７,６０(６):
２１０４Ｇ２１１２．
[１９]㊀姜丛,蒋长胜,尹凤玲,等．基于数据驱动的时间序列b值计算新方法(T b D D):以２０２１年云南漾濞M S６．４地震序列为例
[J]．地球物理学报,２０２１,６４(９):３１２６Ｇ３１３４．
J I A N GC o n g,J I A N GC h a n g s h e n g,Y I NF e n g l i n g,e t a l．An e w
m e t h o df o rc a l c u l a t i n g bＧv a l u e o ft i m es e q u e n c e b a s e d o n
d a t aＧd r i v
e n(T b D D):ac a s es t u d y o
f t h e２０２１Y a n
g b i M S６．４
e a r t h q u a k e s e q u e n c e i n Y u n n a n[J]．C h i n e s eJ o u r n a lo fG e oＧ
p h y s i c s,２０２１,６４(９):３１２６Ｇ３１３４．
(本文编辑:贾源源)
２９４１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年。