长白山火山地震频谱特征研究.pdfVIP

科技论坛 ·65 · 长 白山火山地震频谱特征研究 周文斌 (吉林省地震局敦化地震台，吉林 敦化 133700) 摘 要 ：讨论 了应用快速傅丽叶变换将时间域的地震波信号转换成频率域频谱信亏的基本原理，对火山地震 的类型进行 了分类阐 述。通过对长白山天池火山地震监测数据进行频谱分析 ，确定了长白山近年来发生的火山地震均属于构造型火山地震 ，认为 目前长白山 天池火山正处于火山活动 的早期，尚不具备喷发的危险。 关键词 ：火山；长白山；频谱 ；构造地震 1方法原理 快速傅立叶变换 (FFT)的作用是要是用于将一个有限数字信号 从时间域转换到频率域。什么是时间域呢?简单的可以这样理解 ，就 是一个信号是与时间相关的，用平面坐标来说就是横坐标是以时间 t，纵坐标是时间t上的一个函数值 f(x)，我们叫这个鬲数 f(x)为时间 域。当一个信号经过FFT变换之后 ，将变成频率域。什么是频率域 呢?同样可以这样理解，它其实就是一个用频率与相位来标识的一 个函数。像F(x)=Asin(03x+ )，在数学中A表示幅值，而 03表示频 率 ， 表示相位。我们假设 Ffx)一个频率域函数 ，那么其频率谱就是 JF(x)l= R【(x)+I(x)]“Ix)，其中R与 I表示取实部与虚部函数。需要说明 的一点 ，那就是F丌’变换是一个复数变换 ，所 以F(x)是一个复数函 数。 现代火山地震监测一般采用宽频带数字地震仪，所获取的地震 其 中high—frequency earthquake为构造型火 山地震 LP earthquakes为长周期火 山地震 波形数据实质上就是离散的时间域信号。通过快速傅立叶变换 Harmonictremor是谐 波型火山颤动 (FFT)就可以将这些时间域的地震信号转换成频谱。 图 11985年 9月 9日哥仑 比亚记录到的一张地震 图 2火山地震类型分类 根据地震波形进行火山地震分类的方法主要就是频谱分析。频 谱分析方法主要反映整个观测波形的频谱特性。其缺点是，对浅源 地震 ，随着观测距离的增加 ，地震面波的发育可能对整个频谱产生 明显的影响，致使频谱明 向低频方向偏移。为了防止因为这种偏 移而导致错误的分析结果，一般采用多台频谱综合进行分析。 按照火 山地震的成因或频谱特征，火 山地震可以分为 四大类 ， 即：构造型火 山地震 (Tectonicevent)、混合型火 山地震 (Hybrid event)、长周期火山地震(Longperiodevent)、火山颤动(Tremor)(刘 国明等 ，2006)。 构造型火山地震的发震机制类似于普通的构造地震，是发生在 火山区的由构造应力场或岩浆活动引起岩石破裂产生的地震，也可 因围岩的剪切破裂或者断层运动而产生。其P、S波清晰，波形以短 图2 5个ML3．0级以上火山地震的谱图像 周期成分为主。其优势频率通过在 5-15Hz之间。与典型构造地震 不同的是，构造型火山地震通常以震群的形式发生。由于岩浆可以 在很大的范围内从很深的位置上升 ，因此有些构造型火 山地震 的震 群往往可以持续很长时间。构造型火山地震的急剧增加通常是火山 活动的早期信号。这种类型火山地震也被称为A类火山地震。 长周期火山地震(LP)是在地下应力或者温度梯度作用下，岩浆 活动或者其它流体活动引发的一种地震波形 ，其特点是频率低，P 波明显，很少有 s波 ，而且其波形周期明显比构造型火山地震长，主 频一般小于 2Hz，5Hz以上的频率成人较少。这种火山地震也称为B 类火 山地震 。 混合型火山地震属于二次触发的地震，第一次构造地震产生的 地震波先到达地面被接收到(短周期成分 )，向下传播的地震波又触 发 了下面的岩浆 ，引发 了热液活动 ，产生 的地震