小波分析在地震趋势预测中的应用.docVIP

小波分析在地震趋势预测中的应用 邵辉成1) 杜兴信1) 金学申2) 杜长娥1)1) 中国西安 710068 陕西省地震局2）中国石家庄 050021 河北省地震局 　 摘要 本文应用Morlet小波对近年来华北地区几个引人注目的地震能量序列进行了动态周期分析，发现地震活动既有长期较稳定的活动周期，也存在有一定时限的短期活动周期。并用动态优势周期对未来地震活动趋势进行了估计，结果表明，地震都发生在“特征值”高值附近。 关键词：小波分析 预测 周期 特征值 引言 地震活动具有一定的周期或准周期性，这已为许多事实所证实，一般可用Fourier变换、极大熵谱分析等得到地震活动周期，并根据优势周期，考虑相位，进行地震预测(彭美煊，1996；李学新等，1993)。 但地震的孕育过程是一非线性过程，地震时间序列也是非平稳序列，它在不同时刻携带着不同的地震孕育信息。Fourier变换、最大熵谱分析等方法都是建立在稳态信号基础上的频谱分析方法，仅能给出信号总体所包含的各种频率成分。窗口Fourier 变换(Gabor变换)虽具有一定时频域特征，但它的时—频窗口大小固定不变，也只适合分析所有特征尺度大致相同的各种过程。因此，这些方法不适合用来分析象地震序列那样包含多尺度的信号过程(郑治真等，1996；赵松年等，1997)。 小波分析是80年代发展起来的可同时进行时间域和频率域分析的方法，它能揭示信号不同时刻的频率特征，在信号分析、语音合成、图象识别、地震勘探、震相识别以及气候分析等方面取得了具有科学意义和应用价值的重要成果(刘太中等，1995；刘希强等，1998)，本文尝试了用小波变换来研究地震时间序列中包含的多层次时间周期结构，并讨论其在地震活动预测中的应用，为地震预测提供一种新的思路。 1小波分析及小波函数的选取 具有有限能量的函数f(x)的小波变换以函数簇??()为积分核的积分变换是(即对函数进行伸缩和平移所得的一组函数簇)(赵松年等，1997) ?f)(a,b)=Wf(a,b) ?????????????????? = f(t) ?a,b(t)dt ?????????????????? =f(t)?()dt???? a0，f?L(R) ?????（1） 这样就可把一维函数f(x)变换成a和b的二维函数，其中a是尺度参数，b是定位参数，函数?称为小波。随尺度参数a的变化，函数f(x)的变换方式可看成是通过一个尺度固定的虑波过程，这样以来，尺度因子就可以赋于映射标度一定的物理意义，可以把尺度参数与周期联系起来(周晓兰等，1995；Laura et al.，1995)。通过变换尺度参数a和平移参数b可得到不同时刻的时频域结构特征(杜兴信，1997)。 要从信号中提取有用的信息，应恰当地选择或构造合适的小波函数(高静怀等，1996)，Morlet小波在时频域中有最好的分辩率，是地球物理过程和湍流分析研究中常用的小波。为此，本文采用它对地震能量时间序列进行小波变换。Morlet小波可近似表示为 ( ?(t)= ? ee?????????? ??5??????????????? (2) 是振荡角频率，其实质是一个受高斯包络调制的谐波。(这里 2 小波分析在地震预测中的应用 地震活动周期是孕震信息的一种表现形式，我们可用小波变换来分析地震序列，得到不同时刻地震序列所包含的周期成分和优势周期，并以此来对地震趋势进行预测。但地震活动优势周期有一定的时限性，而地震活动优势周期在什么时间发生突变从地震序列本身目前无法预知，因此我们假定：当前地震活动格局还将维持一段时间，即目前的优势周期在下次地震之前不会发生突变。 地震活动的强弱包括地震所释放的能量，同时包括地震发生的频次，但考虑到资料的均一性，取某一地区一定时间段内地震序列能量平方根，组成时间序列f(t) ?????????????????????????????????????? (3) 其中E是震级为Ms所对应的地震能量(log(E)=1.5Ms+11.8)。 利用小波分析，就可得到不同时刻地震能量序列的优势周期。可以看出，这里主要分析的是地震能量释放的周期。根据周期的意义，某一地震后，如存在一稳定的优势周期，作为预测,那该地震后，周期整数倍所对应的时刻，将是下次地震发生的预测时间，也就是对地震后这些时刻进行了一次预测，但由于地震序列的离散性和随机性，把落在单位时间内预报的次数称为特征值(郑治真等，1993)。下面将分析特征值与地震发生时间之间的关系。 依据以上思路，我们对华北地震区近年来几个显著地震事件(表1)分别进行讨论。 表1、资料的选取 地区 资料时间段(年) 震级下限(Ms) 预测震级(Ms) 采样间隔 陕西地区 1900～1993 3.5 4.5 1年 大同、张北地区