基于数字化观测资料的强地震短期预测新方法的研究.pdfVIP

基于数字化观测资料的强地震短期预测新方法研究 刘希强I沈萍2李红1山长仑1季爱东1张萍3蔡明军4 1山东省地震局济南250011：2中国地震局地球物理研究所北京100081；3辽宁省地震局辽宁 10031；4云南省地震局云南650041) 埘地震波记录中是否含有区域孕震信息．含有多少异常信息的识别程度取决于两个方面的研 究进展：其一是地震观测分辨率的不断提高；其二是现代信号处理技术和方法的深入开发和应用。 随着数字化地震观测技术的提高和台网密度的增加、地震波理论的发展、现代信号处理方法和技 术的提高．依靠直接来自震源和介质的数字化地震波信息，基于地震波参数进行地震预测的问题 又有了新的契机。地震波信号的时频分布既与震源的辐射图像有关，又依赖于传播介质的吸收性 质。地震在孕育过程中．震源区介质性质发生改变的同时，具有”滤波器”作用的介质对地震波的 吸收能力也会发生变化。本文提f{{了定量刻画非稳态地震波信号变化特征的指标和方法。因而， 对相同射线路径或近似射线路径j：的地震波记录开展时频特征研究，可深入地研究其物理内涵， 探索出具有物理基础的强震预测方法。 1．时频能量衰减因子及其偏离度方法 小波变换适用于分析固定比例带宽的非平稳信号，而线调频小波变换可以分析具有固定不变 带宽和非固定比例带宽的非平稳信号(Mann。et a1．，1995)，是小波变换的一种推广，具有潜在的应 用前景。本文基于高斯线调频小波变换(刘希强，2002)提出进行强地震短期预测的时频能量衰减 因子及其偏离度方法。主要思路如下：(1)分别选取P波、s波和包括P波，s波及面波在内的整个 波列，然后分别对数字化地震波数据进行高斯线调频连续小波变换；(2)构造高斯线调频小波变换 CWT(a， 的时间一尺度图序列。尺度图定义为小波变换((w丁(a-b))模值的平方，即：E(a，方)一I 6)I 2。(3)求不同尺度下尺度图对应信号的时频能量衰减因子分布特征。选取一定尺度范围和相 同)由线性拟合方法得到．称为高斯线调频小波变换时频能量衰减因子，其拟合误差称为小波变换 时频能量衰减因子线性偏离度且(下文同)。 2．归一化小波包分频能量比方法 (11对地震波信号进行n层小波包分解，分别提取第n层从低频段到高频段2n个频段成分 Snl④… 的小波变换系数；(2)对小波包分解系数进行重构，提取各频带范围的信号，即s—SnOo r 2dt= 1／2”]，[1／2”．1／2”1]，[(2～1)／2”，1]。(3)求各频带信号的总能量，即：E0=『Isnl(￡)I ：I。小I 雨 ／三r—一 E0『2。(5)建立有、 构造归一化特征向量一=EE,,0／E，E。l／E，…，EⅢ”-f1)／E]，其中Ei／∑J 无中强震孕育条件下的异常特征向量场和正常特征向量场，研究归一化分频能量比(下文同)的 时间和空间变化特征，寻求分频能量的异常指标，用于地震预测。 3．1999年11月29日辽宁岫岩5．4级地震震例研究 日5．4级地震震中为中心，选取营口台、岫岩台和丹东台记录的主震周围200公里范围3级以上 地震i31次；2．5—2．9级地震160次；震前1．0—2．4级地震169次地震波形记录作为研究对象。 从研究地震事件中再选取信噪比较高的记录资料。最后确定能够用于分析的营171台地震事件数是 ·260· 341次、岫岩台335次、丹东台183次。距离主震及其余震分布区由近及远的这三个台站依次是营 口台、岫岩台和丹东台。从三个台站震前地震波参数异常明显程度来看，营口台最明显，岫岩和丹 东依次减弱；震前异常的阶段异常过程表现出相似