数字化台网的近震震相自动识别.PDFVIP

维普资讯 第 24卷 第 1期 西 北 地 震 学 报 Vo1．24 N仉 1 2002年 3月 N0RTHW ESrERN SEISM ()LOG】CALJOURNAL March，ZOOZ 数字化 台网的近震震相 自动识别 朱元清 ，佟玉霞，于海英 ，宋俊高 (上海 市地 震局 ，上海 200062) 摘要：首次提 出丁数字化 台网对 于近震震相 自动 识 别 的 “波形变化值 增 长”算法和 “无 后续震相 ”判捂．前者根据波形 变干匕值来描述近震波形的变化 ，用谊值 的增长捕捉地 震波振幅或频率的突变；后者在识别到初动后考察后续波形的变化情况，如果没有其 它震相 ，初动识别无效．谊判据使算法在初动震相振 幅过小识别不清时，不会误将后 续振幅较大 的震相识别为初动震相．该方法经实例检验表 明精度符舍地震速报 的要 求 ，并鳊制 了相应 的程序． 关t词 ：曩相 自动识暑I；波形变化懂增长；地曩波形 ；近震 中田分类号 ：P315．63 文献标识码：A 文章编号 ：1ooo-0844(2oo2)oi-0008-08 0 引言 在从模拟地震台网向数字化地震台网、地震 台阵发展的过程 中，地震震相的 自动识 0成为 一 项重要课题．所谓震相 自动识别 ，就是当地震波到达记 录系统时，实时地震数据处理系统能 够将各个台站地震波形的各个震相的到时 自动识别 出来．其准确程度直接关系到数字化地震 台网、台阵的智能化 ，它的成熟和完善必将对地震应急和快速响应起到重要作用． 对这一问题 国内外已有不少地震工作者进行了多方面的研究，也取得了一定进展 ，但仍然 是世界地震观测领域中尚未根本解决 的难题．本文从地震波形 的基本特征 出发，对近震 P渡 震相和 s渡震相的 自动识别作进一步研究． 1 震相 自动识别的理论基础 地震发生后 ，各 台站记录中地震波形的初动振幅和频率总是明显地高于脉动，后续震相的 振幅或频率 比前一个震相同样有大幅度 的变化．这就是震相最普遍也是最显著的特征 ]． 一 般定义震 中距小于 1000km 的地震为近震 ，其波形具有 以下特点 ： (1)震 中距小于 100km 的地方震 ，F渡 的周期约为 0．05～O．2S，S渡 的周期约为 0．1～ 0．5s，S渡 的振幅 明显地大于 波 的振幅． (2)震 中距大于 IOOkm 小于 160km 时．初动震相仍然是 p渡 ，它 的周期约为 0．1～1s， S波 的周期约为 0．5～2s，S波 的振幅 明显地大于 F渡 的振幅． (3)震中距大于 16okm时 ，P 、S超前于 p、S．P波成为初动震相 ，周期为 0．5～3s；S 渡周期为 2～5s．震相序列为 P、p…S S，后续震相在振幅上明显大于 P震相． 收稿 日期 ：2001—03-13 作者简舟 ·朱元清 (1953--)，男(汉族)，研究员，从事固体地球物理学和地震观涮技术工作 维普资讯 6 西 北 地 震 学 报 第 24卷 1．1 算法基础 根据地震波的性质和数字化地震记录 ，我们设计如下 的算法 ： (1)在给定采样率后，由数字地震计拾取的地震连续信号经过间隔采样形成数字地震波 形数据 ，它相 当于一种低通滤波． 设如图 1是一段数字化地震波形 ，n、b^ d f是相邻 的 6个采样 点， 、b、C、d。、 、f。为各采样 点在纵轴上 的对应 ，1 dl 点， 、bz、C、dz、 、f2是各点在横轴上 的对应 点．若采样率给 I 定 ，相邻两点间的时间差必然是一个确定值 ，则 (。一b)、(6_一 C)、( —d。)、(d．一 )和 ( 一f-)