中国地震台网观测数据统计_梁姗姗
区域地震台网历史震相数据整合与共享服务
SEISMOLOGICAL AND GEOMAGNETICOBSERV ATION AND RESEARCH第41卷 第5期2020年 10月Vol.41 No. 5Oct. 2020地震地磁观测与研究doi: 10. 3969/j. issn. 1003-3246. 2020. 05. 0250 引言地震台网震相报告是地震台网观测的重要产出成果，也是地震学研究领域的重要资料之一（代光辉等，2019），内容包括震相到时及地震发震时刻、震中位置、震源深度和震级等一系列地震事件相关参数，可以为地壳速度结构、震源参数等地球科学研究工作提供必不可少的数据支持（赵荣国，1993；梁珊珊，2015）。
我国是多地震灾害国家，也是世界上开展地震观测较早、地震研究普及且深入的国家之一。
新中国成立后，特别是1966年邢台地震后，党和国家高度关注地震监测预报工作，区域地震台网迅速发展，到1976年底，中国大陆地区29个省、自治区、直辖市均建成区域地震台网。
经过多年的不懈努力，我国区域地震台网积累了大量观测数据，由于时间跨度较长，观测资料的存储方式差别较大。
早期，震相数据填写在纸介质的“地震卡片”上；1975年以后，震相数据被编辑为“地震观测报告”；2000—2007年，进入模拟记录和数字记录并行阶段，纸介质的地震观测报告和电子版数据存放在各省地震局。
然而，由于纸介质震相数据自然老化日趋严重，保管技术不够完善，部分纸介质出现以下问题：被水浸泡后无法分离；受潮严重，字迹模糊；纸介质发黄变脆，无法重复查阅等。
因此，部分珍贵资料面临无法系统和深层挖掘利用的危险（许建生等，2008）。
随着现代信息技术的发展与普及，纸介质的历史资料已经无法满足当前计算机处理的需求，急需进行数字化转换。
为了方便珍贵历史资料的使用，中国地震局地球物理研究所2014年起组织各省地区域地震台网历史震相数据整合与共享服务刘 伟1） 王庆良1） 王丽艳2） 柴旭超1） 王文青1） 朱飞鸿1）1）中国西安710054中国地震局第二监测中心2）中国贵阳550001贵州省地震局摘要 区域地震台网震相数据是区域地震台网产出的重要成果，是开展地球科学研究的重要资料。
２０２０年新疆伽师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 LV 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 gw＠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 ys 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 ys 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 qu 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 qu 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 qu 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 qu 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 qu 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 mm 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 gu 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 gu 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 yz 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 pl 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 de 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 cf 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 gu 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 pa 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 gi o n ．K e yw 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 ps ://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 pl 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 pl 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 eT e c h n i qu 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 pl 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 pl 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 qu 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 ys 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 ea 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 gu 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 ys Ｇ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 po 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 qu 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 ya 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 qu a k e E n g i n e e r i n g Jo 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 ji 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 es 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 rc 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 nY 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 nd s o u r ce 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 gM 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 og i c a l S o c i e t y o fA 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 fbＧ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 ci 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 wm 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 nd a t aＧd r i ve n(T b D D):ac a s es t u d y of t h e２０２１Y a ng 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,２０２１,６４(９):３１２６Ｇ３１３４．(本文编辑:贾源源)２９４１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年。
雷波断裂带全新世活动证据
第４６卷 第１期２０２４年２月地 震 地 质ＳＥＩＳＭＯＬＯＧＹＡＮＤＧＥＯＬＯＧＹＶｏｌ．４６，Ｎｏ．１Ｆｅｂ．，２０２４ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２５３－４９６７．２０２４．０１．００９张国霞，孙浩越，李伟，等．２０２４．雷波断裂带全新世活动证据［Ｊ］．地震地质，４６（１）：１４１—１６１．ＺＨＡＮＧＧｕｏ ｘｉａ，ＳＵＮＨａｏ ｙｕｅ，ＬＩＷｅｉ，ｅｔａｌ．２０２４．ＥｖｉｄｅｎｃｅｆｏｒｔｈｅｈｏｌｏｃｅｎｅａｃｔｉｖｉｔｙｏｆｔｈｅＬｅｉｂｏｆａｕｌｔｚｏｎｅ［Ｊ］．ＳｅｉｓｍｏｌｏｇｙａｎｄＧｅｏｌｏｇｙ，４６（１）：１４１—１６１．雷波断裂带全新世活动证据张国霞１） 孙浩越１） 李 伟１，２） 孙 稳１，３）１）中国地震局地质研究所，地震动力学国家重点实验室，北京 １０００２９２）中国地震局第二监测中心，西安 ７１００５４３）中铁工程设计咨询集团有限公司，北京 １０００５５摘 要 荥经－马边－盐津构造带是青藏高原东南缘与四川盆地之间的重要边界构造带，虽然晚第四纪以来地震活动频繁，但构造带内各断裂的活动性尚不明确，对认识和评价区域现今地壳变形模式和地震危险性带来了很大的不确定性。
特别是构造带南段的雷波断裂带附近曾发生过多次６级以上强震，但该断裂带的研究程度很低，其最新活动时代的厘定还缺乏可靠依据。
为解决这一问题，文中基于高分辨率遥感影像解译和野外地质地貌调查，在雷波断裂带的北支和南支断裂上分别开展了古地震探槽研究工作。
基于古地震事件识别标志，在２条分支断裂上分别揭示了３次和７次古地震事件，相关放射性碳样品的测年结果限定了北支断裂３次古地震事件的发震时间分别为２１１９０—２０５９０ＢＣ（ＥＰ１）、２０５５０—１２１２０ＢＣ（ＥＰ２）和１２０９０ＢＣ之后（ＥＰ３），而南支断裂最新的２次古地震事件的发震时间为９２７０—５０４０ＢＣ（ＥＬ６）和５０００ＢＣ之后（ＥＬ７）。
古地震探槽研究结果表明，雷波断裂带在全新世以来曾发生过多次地表破裂型强震事件，为全新世活动断裂；其北支和南支断裂揭露的事件互不相同，指示雷波断裂带的分支断裂为独立的发震构造。
辽宁地震前兆台网观测数据跟踪分析探讨
辽宁地震前兆台网观测数据跟踪分析探讨杨振鹏;孙素梅;赵倩;康晓珊;张琪【摘要】地震前兆台网跟踪分析工作的展开,推进前兆台网从观测向观测、应用并重的转变,强化观测数据质量产出,发挥其观测数据在震情监视、地震会商中的作用,使监测与预报更好的结合,有利于服务地震监测和科学研究.统计了2014年以来的前兆跟踪分析不明原因数据资料和预报研究中心地震会商数据资料,通过对比研究,试图探索跟踪分析新的数据记录方式.【期刊名称】《防灾减灾学报》【年(卷),期】2018(034)002【总页数】5页(P46-50)【关键词】地震前兆台网;数据跟踪分析;不明原因异常【作者】杨振鹏;孙素梅;赵倩;康晓珊;张琪【作者单位】辽宁省地震局,辽宁沈阳 110034;辽宁省地震局,辽宁沈阳 110034;辽宁省地震局,辽宁沈阳 110034;辽宁省地震局,辽宁沈阳 110034;辽宁省地震局,辽宁沈阳 110034【正文语种】中文【中图分类】P315.60 引言近年来中国地震局领导高度重视震情工作，指出震情工作是地震工作的核心，一定要推进震情工作的研究。
在监测预报司的领导下，推进地震前兆台网日常工作重心从观测向观测、应用并重转变[1]。
2014年8月开始正式开展前兆台网数据的日常跟踪分析业务，前兆跟踪分析工作的展开，强化观测数据质量产出，加强数据处理分析工作，发挥其按照观测数据在震情监视、会商中的作用，使监测与预报更好的结合，有利于服务地震监测和科学研究。
跟踪分析工作纳入国家台网年度评比，辽宁省前兆台网近年取得较好的名次。
自跟踪分析工作开展以来，已经更新到二代版本，数据异常跟踪分析工作有条不紊的进行，目前辽宁区域前兆数据跟踪分析工作已经参与到预报中心会商中来，借着预报分析改革的东风，成立前兆数据组，并在会商中参与讨论。
但因跟踪分析和地震会商侧重点不同，在会商的作用有限，所以试图探索和讨论跟踪分析新的数据记录方式。
1 辽宁前兆数据跟踪分析不明原因异常2014年国家前兆台网推出“前兆数据跟踪分析软件”在各省区域中心展开，以其作为基本平台，地方前兆台站作为基本操作单元，分析记录形变、电磁、地下流体观测数据曲线变化，形成分析产品，录入前兆数据库。
基于GIS的贵港市雷暴日时空变化特征
基于GIS的贵港市雷暴日时空变化特征杨宇川;梁珊珊;梁文希【摘要】利用贵港市及周边邻近的气象站近30a(1981-2010年)的观测资料,采取统计学方法及GIS技术,对贵港市雷暴日的时空分布特征进行分析研究.结果表明:贵港市近30a年平均雷暴日数为73.4d,最多的年雷暴日数是100d(1983年),最少的年雷暴日数是57d(2009年),年平均雷暴日数呈明显减少的趋势;全年各月均有雷暴发生,各月雷暴日数差异大,呈单峰型分布,峰值出现在7月份;雷暴日大多集中在汛期(4-9月),雷暴出现频率达89.5％;其中6-8月为全年雷暴高峰月.季节上,雷暴主要出现在春夏两季.在空间上,贵港市年雷暴日数自北向南逐渐增多;北部山区的年雷暴日数多在50-70d,中部多在70-80d,南部为80-90d.该研究可为贵港市雷暴的监测预警、评估以及开展防雷检测防灾工作提供了科学的参考依据.【期刊名称】《气象研究与应用》【年(卷),期】2017(038)003【总页数】3页(P122-124)【关键词】贵港;GIS;雷暴;时空【作者】杨宇川;梁珊珊;梁文希【作者单位】贵港市气象局,广西贵港537100;贵港市气象局,广西贵港537100;梧州市气象局,广西梧州543000【正文语种】中文【中图分类】P49雷电是发生在积雨云中的强烈的发电现象，雷电的发生会伴随着瞬时高电压、大电流、强电磁辐射等有巨大破坏作用的物理效应，随着现代社会经济和科技的迅速发展，雷电灾害所造成的人员伤亡和经济财产损失越来越严重［1－2］。
因此研究雷电活动的气候分布，在国民济和国防建设中涉及防雷减灾问题上有重要的应价值，张敏锋［3］等利用经验正交函数（EOF）和主值数分析对我国年平均雷暴日的气候特征及异常变进行了研究，指出近30a来我国大部分地区（除东北地区外）雷暴频数在波动中减少。
马明等［4］对中国雷暴活动特征的时空分布分析表明，中国的雷暴活动表现为在夏季有明显的北进，而冬季则完全南退的季节变化特点。
中国地震台网观测数据统计_梁姗姗
1985 and 2011 for events shallower than 50 km and JB travel-time curves for surface focus
44
地震地磁观测与研究
36 卷
3 震级测定类型及统计
中国地震台网观测报告给出记录事件的振幅和周期值，及相应的震级值，以便读者
到 2011 年底，中国地震台网积累大量震相数据，日常资料分析包括以下震相： ① 地 方 或 区 域 震 相：P、S、Pg、Pn、Sg、Sn； ② 远 震 或 极 远 震 震 相：P、PP、PKP、 PKP2、pP、pPKP、PKS、S、SS、sP、sS、PcP、PcS、ScP、ScS、sPKP、SKS、 SKKS、Pdif。记录分析的各种震相达 5 557 699 条（表 2），大部分为清晰的直达波、核 幔反射波、地表反射波震相和少量地核穿透波震相。中国地震台网记录的各种震相数据 为分析、解释震相的成因和物理意义，为利用各种震相来测定震源基本参数、震源的力 学性质及探讨地球内部构造等研究工作，提供了坚实的数据基础。
10
8
2004-12-26T00:58:51.1 3.15°N 95.79°E
40
8.9
2005-03-28T16:09:34.7 2.03°N 97.05°E
34
8.6
2006-04-20T23:24:59.3 60.97°N 167.48°E
32
8.3
2006-11-15T11:14:11.9 46.69°N 153.27°E
目前，中国地震台网中心编辑出版的《中国数字地震台网观测报告》（月刊）汇编 107 个台站数据。该报告主要提供给中国地震台站、各区域地震台网和国内的科研单位、大 学及研究人员使用。中国地震台网观测报告是人机交互处理，结合地震震相记录特征， 对每一个震相都进行反复调整核对，并对查出问题进行分析、纠正，将每日发生的地震 进行精细分析，处理精度较高。因此，本文应用的统计数据资料完整性、可靠性较高。
玉树MS7.1地震地表破裂带及其同震位移分布
野外地 质地 貌 研 究 表 明甘 孜一 玉 树 断 裂 带 近 ５万 年 来 的 平 均 左 旋 滑 动 速 率 约 为 １ ２
ｍｍ ／ ［ ａ ，晚第 四纪 以来 玉树两侧 （ 柯 段 和 当江 段） 平滑 动 速率 约 为 ７ 邓 水 ．３ｍｍ ／ ｌ 。 ａ３ ＿ ］ １３ ７ ８年 在玉树 县西北 侧隆 宝镇发 生 不小 于 ７级 强震［ ，１ ９ ６ ８ ６年在 玉 树 县南 侧 洛 须至 巴塘 ］
第 ３ Ｏ卷
第 ３期
地
震
Ｖ ０．３ １ ０，Ｎ ｏ ．３
２１ ０ ０年 ７月
ＥＡＲＴＨＱＵＡＫＥ
Ｊ １ ，２ ｌ ｕ． ０Ｏ
玉树 ７ １地 震 地 表 破 裂 ． 带 及 其 同震 位 移 分 布
张军 龙 陈长 云 胡朝 忠 杨 攀 新 ， ， ，
熊仁伟 李 智敏 任 金 卫 ， ，
１ 甘 孜一 玉树 断 裂 带 地 质概 况
印度 板块 以约 ４ ０ｍｍ／ ａ的速率 与欧 亚 板块 的北 向汇 聚造 就 了青藏 高原 广 泛 隆升 和众
＊ 收 稿 日期 ： ０ ０ ０ — ５ ２ １ — ５０ ；修 改 回 日期 ：２ １ — ５ １ ０ ００—５
基 金 项 目 ：国家 重 点 基 础研 究 发 展 计 划 项 目（０ ８ Ｂ ２ ７ ） 助 ２０ Ｃ ４ ５ ０ 资 作 者 简 介 ： 军 龙 （ ９ ３） 张 １ ７ 一 ，男 ，山西 阳城 人 ， 研 究 员 ， ０ ５年 获 硕士 学 位 ，主要 从 事 地 震 地 质 等 研究 。 副 ２０
据 中国地震 台 网测 定 ，２ １ ０ ０年 ４ １ 月 ４日７时 ４ ９分 ３ 秒 （ ７ 北京 时 间） 青海 省玉树 县 附 ，
广西测震台网台站地方性震级偏差分析
广西测震台网台站地方性震级偏差分析
牟剑英;向巍;李攀攀;徐宁
【期刊名称】《地震地磁观测与研究》
【年(卷),期】2022(43)2
【摘要】选用广西数字化测震台网2008—2020年记录的广西及邻区1251个M_(L)≥2.0地震事件震相资料,采用震级残差统计方法,获取广西51个地震台单台震级与台网平均震级的平均偏差和标准差,分析观测仪器类型、台基场地响应和震中距等因素对震级偏差的影响,进而提出台基校正值,修正区域量规函数。
结果发现,南宁、天等、北海等15个地震台站的震级偏大,桂林、贺州、天峨等16个地震台站的震级偏小,崇左、灵山、平果等20个地震台站震级适中;在单台震级共14784个样本中,震级偏差集中在-0.4—0.4,呈正态分布;震级偏差受台基场地响应和震中距因素影响明显,当场地响应值大于1.0时震级偏大,小于1.0时震级偏小,等于1.0时震级较为适合,当震中距在70 km范围内时单台震级偏小,当震中距在70—300 km时单台震级较为适合,当震中距在300—550 km时,单台震级稍有偏大。
修正后的区域量规函数和台基校正值,可有效提高广西测震台网测定震级的准确性。
【总页数】9页(P18-26)
【作者】牟剑英;向巍;李攀攀;徐宁
【作者单位】广西壮族自治区地震局
【正文语种】中文
【中图分类】P31
【相关文献】
1.福建测震台网单台震级偏差与台站场地响应关系的初步分析
2.内蒙古数字测震台网单台近震震级偏差研究
3.云南测震台网台站场地响应与ML震级测定的讨论
4.江西测震台网新、旧震级国家标准测定地方性震级的对比分析
5.河北测震台网单台震级偏差分析
因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
0 0 20 40 60 80 100 120 140 160 180
ᡰሏ/Ń
0 0 20 40 60 80 100 120 140 160 180 ᡰሏ/Ń
图 2 中国地震台网观测报告主要震相走时曲线及表面源 JB 模型走时曲线
Fig.2 Travel-time picks collected by the report of National Seismological Network of China between
ᡌ/s ᡌ/s
3600
3600
3200
3200
2800
2800
2400
2000
1600
1200
800
SCS SCP
400 PCP
SS SKKS
S
SKS diff
2400
2000
1600
1200
800
SCS SCP
400 PCP
SS
S PP
P
SKKS
SKS SKP PKP
Pdiff
1985 and 2011 for events shallower than 50 km and JB travel-time curves for surface focus
44
地震地磁观测与研究
36 卷
3 震级测定类型及统计
中国地震台网观测报告给出记录事件的振幅和周期值，及相应的震级值，以便读者
9
8
2007-01-13T04:23:22.7 46.29°N 153.80°E
10
8.1
2007-08-15T23:40:57.8 13.40°S 76.60°W
39
8.2
2007-09-12T11:10:19.3 4.90°S 101.41°E
25
8.6
2007-09-12T23:48:56.9 2.85°S 100.67°E
0 引言
地震台网观测报告是地震学研究领域的重要资料之一，其中包含各种震相到时并提 供一定精度的地震事件相关参数（如发震时刻、震中位置、震源深度和震级等）并且可以 为速度结构，震源物理等相关研究提供基础性数据资料。故地震观测报告是进行地球科 学研究工作必不可少的观测资料。
中国最早的地震观测报告见于 20 世纪初。19 世纪末，法国天主教会在上海徐家汇 设立观象台；1904 年增设地震观测项目，同年开始汇编地震报告（国家地震局科技监测 司，2001）。随后几十年，中国开始自行设立地震台，老一辈科学家在极其艰苦条件下， 悉心收集、整理地震资料。随着时代变迁，台址和名称有所变动，但整理汇编的地震报 告基本是连续的。全国范围的地震台网较正规的地震报告，始于 1954 年开始（赵荣国， 1988）。之后，地震观测技术迅速发展，国际地震机构定期出版地震报告、地震目录及相 关资料，按照国际惯例，部分台站地震资料进行国际交换，促进国际交流与合作，为国 际科技合作服务。
目前，中国地震台网中心编辑出版的《中国数字地震台网观测报告》（月刊）汇编 107 个台站数据。该报告主要提供给中国地震台站、各区域地震台网和国内的科研单位、大 学及研究人员使用。中国地震台网观测报告是人机交互处理，结合地震震相记录特征， 对每一个震相都进行反复调整核对，并对查出问题进行分析、纠正，将每日发生的地震 进行精细分析，处理精度较高。因此，本文应用的统计数据资料完整性、可靠性较高。
2001-11-14T09:26:09.8 35.92°N 90.53°E
11
8.2
2002-11-03T22:12:41.0 63.50°N 147.40°W
4
8.1
2003-09-25T19:50:05.0 41.79°N 143.71°E
16
8.2
2004-12-23T14:59:04.4 49.30°S 161.30°E
20
8.7
参考地点
中国青海省（昆仑山口） 美国阿拉斯加州中部 日本北海道 麦夸里岛以北 印尼苏门答腊岛北部西岸远海 印尼苏门答腊岛北部 俄罗斯西伯利亚东部 千岛群岛 千岛群岛 秘鲁沿岸近海 印尼苏门答腊岛南部
印尼苏门答腊岛南部
四川汶川县 萨摩亚群岛 智利中部沿岸近海 日本本州东岸近海
ᄬ
10000
9000 8000 7000 6000 5000
作者简介：梁姗姗（1989—），女，助理工程师，主要从事数字地震资料处理与解释工作 基金项目：留学回国人员科技项目择优资助 本文收到日期：2015-02-14
42
地震地磁观测与研究
36 卷
供了宝贵的研究资料。由地震频度分布图（图 1）可知，2001 年后，中国数字地震台网记 录地震事件数目大幅提高，记录 7 万多次地震事件。自 2001 年以来，全球地震活动水平 明显偏高，尤其是 8 级以上强震多发，2001—2011 年共发生 16 次 8 级以上地震（表 1）， 尤其是 2004 年 12 月 26 日苏门答腊 8.9 级地震和 2008 年 5 月 12 日汶川 8.0 级地震，较 强余震频发。
给出震级的离散程度。
地震台数字化改造之前，直接应用模拟位移记录测定各种震级。随着我国数字地震
观测技术的发展，“十五”期间，中国地震局完成“中国数字地震观测网络”项目建设，将
模拟地震观测系统进行数字化改造，并新建一批数字化地震台站（刘瑞丰等，2008）。由 于各数字地震台站观测仪器的类型和频带不同，观测地动速度，而非地动位移，所以，
2 震相统计
震相分析是地震学的基础工作之一，主要任务是从地震图上读取、测定与地震参数 相关的各种震相读数。地震台网提供的震相读数，作为原始数据，在地学研究领域具有 广泛的学术意义和使用价值（赵荣国，1988）。
第4期
梁姗姗等：中国地震台网观测数据统计
43
1985 年—2011 年，中国地震台网观测报告提供各类地震波震相，包括地方震、区域震、 远震以及深震，地震报告不包括非天然地震。2000 年以前，我国大部分台站为模拟地震 记录，记录仪器主要为短周期地震仪和中长周期地震仪，通过分析模拟记录来识别震相 并进行震级计算。自 2001 年起，国家数字地震台网投入使用，为了保证观测数据的连续 性与完整性，分析人员根据短周期和中长周期地震仪特性，将宽频带速度型数字记录仿 真成标准的模拟记录，如：长周期（763）、中长周期（SK）和短周期（DD-1）模拟地震记录， 分析人员在仿真后的波形记录图上识别震相和计算震级等。
中国地震台网观测数据统计
梁姗姗 黄志斌
（中国北京 100045 中国地震台网中心）
摘要 收集了由中国地震台网中心编辑出版的中国地震台网观测报告，对报告中提供的 震相和震级数据进行统计，并详细介绍中国地震台网测定的震相类型、震相数量，以及 观测报告中测定的各种震级标度和震级计算公式。本文目的在于迎合地震监测技术，为 地震研究和预报工作提供高质量、高水平的地震观测报告。 关键词 中国地震台网；震相；震级
第 36 卷 第 4 期 2015 年 8 月
地震地磁观测与研究 SEISMOLOGICAL AND GEOMAGNETIC
OBSERVATION AND RESEARCH
doi: 10. 3969/j. issn. 1003-3246. 2015. 04. 008
Vol.36 No. 4 Aug. 2015
到 2011 年底，中国地震台网积累大量震相数据，日常资料分析包括以下震相： ① 地 方 或 区 域 震 相：P、S、Pg、Pn、Sg、Sn； ② 远 震 或 极 远 震 震 相：P、PP、PKP、 PKP2、pP、pPKP、PKS、S、SS、sP、sS、PcP、PcS、ScP、ScS、sPKP、SKS、 SKKS、Pdif。记录分析的各种震相达 5 557 699 条（表 2），大部分为清晰的直达波、核 幔反射波、地表反射波震相和少量地核穿透波震相。中国地震台网记录的各种震相数据 为分析、解释震相的成因和物理意义，为利用各种震相来测定震源基本参数、震源的力 学性质及探讨地球内部构造等研究工作，提供了坚实的数据基础。
震相 数目
震相 数目
P 1764992
PKS 18899
Pg 133186
S 758627
Pn 144006
Sg 136532
Pdiff 4833
Sn 82379
pP 496685
sS 218551
sP 407462
ScS 135259
PcP 163064
ScP 83576
PcS 84007
SS 233841
目 前， 中 国 地 震 台 网 中 心（CENC）承 担 中 国 地 震 台 网 观 测 报 告 的 处 理 与 编 辑 工 作， 利用有人值守单台记录资料，编辑和出版了《中国数字地震台网观测报告》和《中国地震台 站观测报告》，并提供在线和离线服务。
1 资料来源
研 究 使 用 的 观 测 资 料 主 要 来 自《中 国 地 震 年 报 》和《中 国 数 字 地 震 台 网 观 测 报 告 》 （1985 年 1 月—2011 年 12 月）。中国地震台网观测数据记录了全球范围内 86 944 次地震， 时间跨度较大、空间分布广。数字化技术采用以来，随着大量数字化台站建设及模拟台 站改造，台站观测水平和分析能力逐步提高，地震观测数据大量增加，为地震学研究提
15
8.2
2008-05-12T06:27:59.5 31.01°N 103.42°E
14
8.2
2009-09-29T17:48:10.7 14.80°S 171.57°W
16
8.1
2010-02-27T06:34:14.0 35.80°S 72.80°W
35
8.8
2011-03-11T05:46:19.0 38.10°N 142.50°E