时频分析在地震信号初至估计中的应用研究.pdfVIP

时频分析在地震信号初至估计中的应用 许朝阳 范万春 季国华 刘俊民 (西北核技术研究所 西安710024) 摘要时频分析旨在构造一种时间和频率的联合密度函数，继而揭示信号所含频率分量及其演 化特性，是传统Fourier分析的新进展。本文通过对地震信号噪声特性的分析，研究了基于联合时 频谱图三重相关的地震信号初至的估计方法，并给出了实验结果和结论。 关键词 时频分析 三重相关 时频谱图 初至估计 1 引言 一般地，地震学把地震信号的P波到时作为地震信号的初至。地震信号的初至估计作为 地震信号检测的重要组成部分，对于地震事件的定位意义重大。目前，地震信号的检测包括 人工检测和利用计算机自动检测。人工检测就是地震分析人员根据所掌握的对地震波形了解 的先验知识，从直观的地震图中挑选出可供分析的事件震相段并加以分析。由于地震信号和 台站本底噪声构成的时间序列属于时变非平稳信号，一方面，震中距较远，震级较小的地震 在台站记录的信号信噪比可能很低，另一方面，台站附近以及记录系统本身也可能出现远远 高于地震噪声本底的干扰。基于这些因素，人工检测往往无法准确地进行地震信号的初至估 计，而且有时无法判断地震信号的初至。 近年来兴起的联合时频分析技术，是一种刻画非平稳信号的有效手段。时频分析技术对 信号同时在时域和频域进行分析，不仅能够揭示信号中包含的频率分量，而且能够表征每一 分量随时间的演化过程。本文通过对地震信号的噪声特性的分析，研究了基于联合时频谱图 三重相关的地震信号初至的估计方法，并给出了实验结果和结论。 2地震信号及其噪声的特性分析 地震信号是指地震波的振动信号，噪声指环境噪声和信道噪声等。研究地震信号的检测 技术，首先需要了解地震信号记录中的噪声情况[1]。 2．1 环境噪声 环境噪声从大的方面可以分为短周期的人为干扰和海洋、气象因素引起的2秒以上的较 长周期噪声两种。人为干扰主要是交通和生产活动所致；气象、海洋、河流、湖泊等因素引 起的周期在2～8秒的振动称为脉动，其中最为主要的是由于海浪等面波从海岸传至内陆造成 的。 一37— 2．2信道噪声 除环境噪声外，另一类噪声主要是在信号传输过程中出现的。这类噪声是一种随机过程， 噪声与传输的信号之间没有相关性，大致可以分为两类：一类是各态历经的，噪声的瞬时值 是随机的，但是有确定的平均特性，这类噪声是由大量的短促脉冲的叠加形成的，这些脉冲 的产生是随机的，但其统计特性具有一定的规律性，通常称为起伏噪声。另一类噪声是非各 态历经的，可分为规则的与无规则的两种，前者具有特征性的波形，如电源噪声、内燃机点 火系统的感应等；后者是指个别的脉冲干扰，如某些天电干扰、电器干扰等。 2．3设备噪声 由于地震记录系统一般包括机械装置、机电转换、电信号放大、模数转换等多个组成部 分，每个部分都多多少少对系统本身引入自身的噪声，并且这些部分又是前后衔接的，各部 分噪声有一个累加甚至是放大的过程，最后总和的噪声对地震信号的影响往往不可忽视。 可见，地震信号的背景噪声是随机的，它在地震记录上表现为杂乱无章的波形，分布于 整个记录中，频带很宽，传播速度不确定，且无一定的传播方向；但遵循某种统计规律，在 地震勘探中存在的随机干扰噪声是具有各态历经性质的平稳随机过程[z]。 由于地层的吸收和衰减作用，地震信号的高频能量随着传播距离的增加迅速减弱，而随 机噪声在整个时间域中的分布几乎是平稳的，其特性描述如下： (1)随机噪声的谱不是白化的； (2)随机噪声在二维频谱中是分布于全局的； (3)有色随机噪声的道与道之间的互相关不等于零，只要它的平均振幅大于信号，就会 严重地影响到剩余静校正及速度分析的精度，这时噪声和信号很难分离； (4)随机噪声在时间域中的分布几乎是平稳的。 图l地震信号噪声特性 图2地震信号有用分量的特性 图1和图2分别显示了采样率为40Hz的天然地震信号垂向记录的噪声和有用分量的幅 频响应曲线。 通过对图1和图2的分析，可以得出如下结论：地震信号的噪声分量基本符合随机噪声 的特性，其幅频响应均匀地分布在o～14Hz之间；地震信号为非对称概率分布，地震