三峡重点监视区构造地震与矿震时频谱特征分析.pdfVIP

第 31卷 第4期 地 震 地 质 Vo1．31．No．4 2009年 12月 SEISMOLOGY AND GEOLOGY 1)ec．．2009 三峡重点监视区构造地震 与矿震时频谱特征分析 张丽芬 ’’、 廖武林 曾夏生。 钟羽云 ’ 1)中国地震局地震研究所，武汉 430071 2)中国地震局地球物理研究所，北 京 100081 3)湖北中南勘察基础工程有 限公司，武汉 430081 4)浙江省地震 局，杭州 31O013 摘 要 基于三峡数字遥测地震台网记录的数字化波形资料 ，利用谱图时频分析方法 ，研 究了三峡 地 区构造地震和矿震波谱 的时频特征，发现二者的时频谱特征存在着明显 的差异 ：1)构造地震波主 频高于矿震 ；2)构造地震波的频率成分比矿震丰富，能量密度频谱随频率轴展布较宽，前者带宽足后 者的4倍 ；3)矿震能量强度 的峰值在整个波列中出现的时问较构造地震早 ，且能量衰减较快 。 关键词 构造地震 矿震 时频分析 三峡 中图分类号 ：P315．3 1 文献标识码 ：A 文章编号 ：0253—4967(2009)04—0699—08 O 引言 长江三峡工程库区的地震类型较多，大致可 以分为构造地震 、水库诱发地震 、矿震、浅表卸 荷型地震等。然而该地区的地震多属于小震、微震，地震发生后很难快速判识地震类型 ，又由于 不同类型地震序列的特征和震级上限不同，这就给地震趋势预测带来了一定的困扰 。随着三峡 水库后期蓄水的提前 ，今后2～5a内，甚至更长时期 内三峡水库因蓄水而引起的库区地震灾害 的危险性还难以预料 ，所以库区地震类型的识别对于地震趋势的预测研究具有十分重要的科学 意义和现实意义。 许多研究者借助于地震信号这一有利工具对地震类型识别进行过不同程度的研究，曹肃朝 (1993)和林怀存 (1990)分别利用河北北部和鲁中南地震台网的资料，对构造地震和塌陷地震 进行了分析，发现二者的震相特征存在着明显的不同。张萍等 (2005)根据辽宁数字地震台网 记录的爆破 、矿震与地震的数字资料 ，采用波谱分析的方法，发现爆破 、矿震与地震的纵横波拐 角频率及最大谱值等存在差异。和雪松等 (2006)利用小波包识别地震和矿震 ，赵永等 (1995) 从地震记录图的角度对北京地 地震 、爆破和矿震的特征进行波形识别，刘希强等 (2005)选取 Morlet小波作为基函数 ，利用能量衰减因子 (。)的方法研究辽宁省抚顺和北镇台记录的矿震和 地震的a值变化特征及其随震级变化的特点 。 [收稿 日期] 2009—01—07收稿 ，2009—11—09改回。 [基金项 目] “十一五”国家科 技 支撑 计划 项 目(2006BAC01B040202)、国家 自然科 学 基金 重点 项 日 中国地震局地震研究所所长基金 (IS200756052)共同资助 。 700 地 震 地 质 但鉴于波谱特征与震源和接收台站之问的射线路径介质的非均匀性的相关性 ，所得出的判 别指标往往具有一定的区域局限性。所以本文想利用三峡遥测地震 台网提供的数字地震资料 ， 通过能在时问、频率、能量域联合反映信号特征的非平稳时间序列的时频分析方法 ，对三峡地 区 构造地震和矿震的识别做一些初步研究。 1 时频分析方法 时频分析方法 (胡 昌华 ，2002)是指用时间和频率的联合函数来表示非平稳信号 ，并对其进 行分析和处理的一种方法。常用的线性时频表示有短时傅立叶变换 (STFT)、小波变换 (wT)和 戈勃 (Gabor)展开等。而非线性时频表示则是二次型的能量分布表示 ，常用的非线性时频表示 有 wV(Wigner—Ville)时频分布、马根诺一希尔 (Margenau—Hil1)时频分布 、乔伊一威廉斯 (Choi． Willimas)时频分布等。