华北地震活动性论文 小波变换论文

华北地震活动性 论文 政研论文下载论文大学下载论文大学下载关于长拳的论文浙大论文封面下载 小波变换论文 华北地震活动性论文 小波变换论文 运用小波变换方法研究华北地区的地震活动性 摘要:运用小波变换方法，以不同的小波时间尺度 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 了华北地区1480,2009年M?5.0和1970,2009年M?4.0地震的贝尼奥夫应变数据，得到了不同时间尺度的地震活跃期和平静期。分析表明:目前华北地区从1480年以来的M?5.0地震活动刚过活跃期，处于应力积累阶段，而1970年以来M?4.0地震存在7年左右的活动周期，目前正在进入4级以上地震的活跃期。 关键词:贝尼奥夫应变;小波变换;地震活跃期;地震平静期;华北地区 0 引言 表征地震活动性的各种参数从不同角度对地震活动的特征进行了完整、系统的描述，成为地震预测预报的重要参考依据，但是地震活动在时间轴上的分布是否有规律可循，多年来一直是个有争议的问题。马宗晋等(1982)、张国民等(1987，2000)、张培震等(2003)对我国的地震活动作了分析，发现在不同的时期，地震活动显示了平静和活跃变化;同一地震区内的地震活动具有某种程度的自相似性，在时间分布上具有某种“似周期性”分布特征;在不同时期地震发生的频度和强度有较明显的不同(时振梁等，1997)。Benioff(1951)曾把1900,1950年全球的地震活动分为5个时段，指出地震活动存在地震活跃期和相对平静期;马宗晋等(1982)指出地震活动强弱随时间分段现象在我国大陆地区是相当普遍的，而且强弱变化的时间 范围长短不一;张国民和张培震(2000)、傅征祥等(2000)讨论了华北地区M?7.0强震的时间分布，发现在每一个地震活跃期内大地震的积累频度曲线随时间的变化大致呈指数函数形式;万永革等(2003)也分析了中国大陆地区M?5.0地震活动的周期性，得到了不同时间尺度的地震平静期和地震活跃期;刘虹和苏有锦(2007)研究了1970年以来云南地区MS?7.0地震的各种活动参数，发现震源区地震活动水平在大地震前后有很大的差异;苏有锦等(2001)研究了20世纪以来云南地区MS?5.0地震的活跃期与平静期的交替过程。总之，活跃期和平静期的划分对地震预测预报，特别是强震的中短期分析预报甚至短临预报工作具有十分重要的意义，是我国地震预测研究中一个不可缺少的环节(马宏生等，2002，2005)。 华北地区作为中国内陆的地震多发区，一直备受地震学家的关注。如何从不同尺度上划分活跃期和平静期，并在大尺度活跃期(或平静期)中提取小尺度的活跃期和平静期，进一步为地震预测预报提供依据,在传统的Mt图上难以找到地震活动演化过程突变点。小波分析方法具有能够根据分析对象不同而自动调整有关参数的“自适应性”和能够根据观测对象不同而自动“调焦”的特性，因而被广泛应用于地震研究的各个领域(刘希强等，1998;吴忠良，1999;邵辉成等，2000;徐义贤，王家映，2000)。笔者采用小波变换方法，从不同时间尺度和震级尺度研究了华北地区的地震在时间上的分布规律，以期为地震预测预报提供参考。 1 小波变换的原理 一个地震时间序列f(t)是由不同强度的地震构成的，在t时刻地震的强度可以表示为 f(t)~ta.(1) 其中a为奇异性指数或标度指数。奇异性是指t?0时a的数值(特别是a<0)决定了f(t)?0的速度。标度性是指若变量t变成λt，则f(t)和f(λt)是自相似的，即[KH*1] f(λt)=(λt)a=λata=λaf(t).(2) 式(2)常称为标度律。自相似的标度律是多种尺度(即无特征尺度)现象的特征。地震资料的时间序列含有多种不同尺度。直观地看f(t)是杂乱无章的，但小波变换好比显微镜，调节其中的放大倍数就可以清楚地看出f(t)在各个层次上的变化趋势。小波变换Tg是将f(t)分解成具有局部特性的小波g(t): Tg(a,t0)=1a??[KG-1*4/9]?-?f(t)g(t-t0a)dt.(3) 其中，1ag(t-t0a)小波是将具有局部特性的小波函数g(t)通过平移t0和放大(放大倍数为1,a)而构成的。参数a具有时间的量纲，也称为小波尺度。 设θ(t)为具有低通性质的平滑函数，以它的一阶、二阶导数作为小波对f(t)作小波变换，可以证明(Logan，1977) 式中，*表示卷积。也就是说，用φ1(t)、φ2(t)对f(t)作小波变换， 分别相当于f(t)被θ(t)平滑后对t求一阶或二阶导数。因此，对某一固 定a值，[KG*4]f(t)*θa(t)的拐点既是Wφ[HT4.5]1[HT5]f(a,τ)的极值点， 又是Wφ[HT4.5]2[HT5]f(a,τ)的过零点， 由此可检测地震活动信号的急剧变化之处(Mallat，1991;Mallat，Hwang，1992)。 本研究中采用具有低通滤波性质的平滑高斯函数，根据它的一阶导数、二阶导数作为小波基函数进行突变点分析。小波函数的表达式为 选择该函数是因为它对称、可微、可积，时域、频域都是高斯型且呈平方型指数衰减，在时、频两域均具有很好的局域性。运用小波变换可以准确地识别出突变点，运用高斯函数作为小波基的小波变换相当于对原来的信号做平滑处理，而运用高斯函数的一阶导数作为小波基的小波变换可精确地找到信号所对应的突变点:若脉冲向上则表示信号增大，反之表示信号减弱(万永革，2007)。如果使用傅立叶变换，就只能从频率轴来分析该信号是由哪些频率的波组成的，它只是将这些信息铺开到整个频率轴上。但是地震资料f(t)组成的序列并不是平稳的，而是随时间不断变化，傅立叶变换并不能提取f(t)中的奇异性和突变点的信息。 2 地震资料选取 应变与构造应力和地震活动是联系紧密的物理量，在地震分期研究中较为常用(夏浩明，1987)。为了表征某个地区在不同时期的地震活动特征，根据该地区的应变积累和释放变化，可将强震活动划分为若干周期，每一周期包括一个应变能释放高潮和低潮。一个地区 的地震能量释放是局部状态的表征，由于地震辐射能量的平方根与地震的应变释放成正比，在研究中常用应变释放来表征地震能量释放。根据历史地震目录可以推算出地震应变释放的大小。 本文采用贝尼奥夫曲线(Benioff，1951)来表征华北地区地震应变能的释放。计算贝尼奥夫应变时采用时振梁等(1997)的公式: εb(t)=[KG-*3]?10(4.8+1.5M)/2.[JX-*2] (7) 式中，εb(t)为以t时刻为中心的0.1 a之内的贝尼奥夫应变积累，即某一地区0.1 a内所有大于M级地震对应变释放的贡献;M为震级。地震目录资料的选取对研究结果有很大的影响。华北地区的地震资料在1480年以后是比较完整的(黄玮琼等，1994)，因此笔者结合《中国历史强震目录(公元前23世纪—公元1911年)》(中国地震局震害防御司，1995)和中国地震地震台网中心地震目录选取了1480~2009年华北地区(110?,120?E，34?,42?N)1 449个M?4.0的地震进行研究(图1)。 3 地震资料的小波分析 贝尼奥夫应变曲线表述了一个地区地震能量积累的特性，通过研究贝尼奥夫应变曲线可以反映该区内地震能量的释放变化过程，但是用贝尼奥夫应变不能明显地划分出不同时间尺度和层次上地震活动的平静期和活跃期。笔者分别对华北地区1480,2009年M?5.0的地震和1970,2009年M?4.0的地震进行计算，得到了两组贝尼奥 夫应变数据，然后分别以高斯函数及其一阶导数为小波基进行小波变换，分析不同尺度的平静期和活跃期。 图2为研究区内1480,2009年M?5.0地震的贝尼奥夫应变曲线及其相应的小波变换。图2a最上端一条是贝尼奥夫应变去均值的原始资料曲线,下边为不同时间尺度以高斯函数为小波基的小波变换曲线;图2b是不同时间尺度以高斯函数一阶导数为小波基的小波变换,图中数字是对应的小波时间尺度。从小波变换曲线能够看出，总体上华北地区自1480年以来M?5.0地震大约呈300 a的周期性变化，无论采用哪种小波尺度，这种周期性都较为明显。在过去的530多年时间内华北地区大致经历了两个活动期，前一活跃期高峰在1668年郯城8.5级地震前后，后一活跃期高峰在1976年唐山7.8级地震前后。从小波变换的曲线形态来看，目前华北地区处于应变积累阶段，且处于应变较低状态。由于历史地震资料主要都是来源于《中国历史地震目录(公元前23世纪—公元1911年)》，所以笔者对这些地震资料的分析结果跟张国民和张培震(2000)、马宏生等(2005)结果一致，说明小波变换的分析结果是可信的。华北地区刚刚经历强震活跃期，不会有巨大地震或群集性中强地震发生。从以高斯一阶导数为小波基的分析还能看出，自1800年以后，当小波时间尺度小于5 a时，有4个明显的突变点对应4个活跃期，第一个活跃期对应1830年6月12日河北磁县7.5级地震，第二个对应1888年6月13日渤海湾7.5级地震，第三个对应1937年8月1日山东菏泽7.0级地震， 第四个对应1966年3月8日河北7.2级地震和1978年7月28日唐山7.8级地震，当小波时间尺度大于5 a时，这两个地震对应的曲线突变点融合为一个突变点。 从小波变换曲线还可以看出:? 随着时间尺度的增加，极值点逐渐变少，且相邻的极值点逐渐融合成较大的极值点;? 各极值点的幅值反映了以同尺度高斯函数为小波基变化对应点局部梯度的大小，极值点的正负反映了滤波后贝尼奥夫曲线的边沿方向特性——正极大值点对应上升沿，负极大值点对应下降沿;? 极值点随尺度变化的幅值大小反映了曲线边沿变化特性——幅值较大反映滤波后贝尼奥夫曲线较剧烈的变化趋势(较陡的边沿)，较小的幅值对应滤波后贝尼奥夫曲线较平缓的变化趋势(较平的边沿)，由此可得到地震活跃期的开始年代(地震活跃突变点);? 当尺度增大时极值点消失或成对结合成新的极值点，这反映了信号中隐含的某种“嵌套结构”。 运用小波变换能够较好地分离和刻画不同时间尺度的地震活动性，探测序列突变。地震活动性的强弱，是相对于不同时间尺度而言，因此，地震活动的活跃和平静分期的问题只有建立在具体的时间尺度上才有意义。不同时间尺度上的地震活动，可能与地壳块体的整体应力应变、不均匀性，以及与邻区相互作用有密切的关系。笔者选用华北地区1970,2009年的M?4.0地震，是因为这个时间段内M?4.0的地震已被非常准确、系统、完整地记录下来，且M?4.0地震对中 强震的发生有很大的影响，选取这个震级更符合实际的地震能量的释放演化过程。 从图3可以看出，自1970年以来，在0.5,5 a的时间尺度上华北地区M?4.0地震一共有5个活跃期，大致呈现出7 a左右的周期，而且其强度随时间在变化，随着小波尺度的增加，地震活跃期突跳更加平滑。第一个活跃期对应的是1976年河北唐山7.8级地震前后，第二个对应1983年山东菏泽5.9级地震前后，第三个对应1989和1991年山西省大同6.1、5.8级地震前后，第四个对应1998年河北张北6.2级地震前后，第五个对应2004年河北邢台5.0级和2006年河北文安5.1级地震前后，但是这个活跃期比较弱，当小波时间尺度增加到4 a时已经很不明显。根据曲线分析，1998,2009年研究区内缺少4级以上地震，在5 a以下的时间尺度上，研究区内目前应该正在进入4级以上地震的活跃期。 从图3还可以看出，每个活跃期的出现总是具有双峰特点的，当小波尺度小于3 a时，这个特点非常明显，1970年以来一直很有规律性，且以高斯一阶导数函数为小波基时比以高斯函数为小波基时更加明显，这也说明在分析突变时，小波基的选择对结果有一定的影响，选择合适的小波基函数可以更好地突出地震活动的规律。 4 结论与讨论 (1)小波分析方法能得出不同时间尺度的地震活动特征，反映地震演化过程。运用该方法不仅可以将隐含在序列中的各种随时间变 化的周期显现出来，还能确定地震活动突变点的位置，该位置对应于地震活跃期和平静期的转折处。 (2)华北地震活动在不同的时间尺度上表现出不同的特征，华北地区M?5.0地震活动出现两次高峰期，周期约为300 a，且刚刚经历了1个活跃期，处于应力积累的平静阶段，该地区目前不会有巨大地震或群集性中强地震发生。1970年以来，研究区M?4.0地震活动的周期大约为7 a左右，而且其强度随时间在变化，目前正在进入4级以上地震的活跃期。因此，对于地震活跃或平静期的讨论，只有建立在具体的时间尺度上才有意义。 对中国地震局地球物理研究所万永革研究员和陈化然研究员给予的指导表示衷心的感谢～ 参考文献: 傅征祥，刘杰，刘桂萍，等.2000.中国大陆成组强震活动水平的三分法及其应用,J,.地震，20(增刊):34-38. 黄玮琼，李文香，曹学锋.1994.中围地震资料完整性研究之二,J,.地震学报，16(4):423-431. 刘虹，苏有锦.2007.云南地区7级大震震源区地震活动特征研究,J,.地震研究，30(1):18-23. 刘希强，周蕙兰，郑治真，等.1998.基于小波包变换的弱震相识别方法,J,.地震学报，20(4):373-380. 马宗晋，傅征祥，张郢珍，等.1982.1966,1976年中国九大地震,M,.北京:地震出版社. 马宏生，刘杰，张国民，等.2002.用应变积累释放模型研究中国大陆地块分区地震活动,J,.地震学报，24(6):569-578. 马宏生，刘杰，张晓东，等.2005.基于G—R关系的应变积累释放模型研究中国大陆强震的分区活动,J,.地震学报 27(4):355-366. 苏有锦，李忠华，刘祖萌，等.2001.20世纪云南地区MS?5.0级地震活动的基本特征,J,.地震研究，24(1):1-9. 邵辉成，杜兴信，金学申，等.2000.小波分析在地震趋势预测中的应用,J,.中国地震，16(1):48-52. 时振梁，汪良谋，傅征祥，等.1997.中国大陆中长期强震危险性预测方法研究,M,.北京:海洋出版社. 吴忠良.1999.鉴别地震活动周期性和随机性的广义Rydelek Sacks检验,J,.地震学报，21(5):487-494. 万永革，齐福荣，孟晓春，等.2003.中国大陆及华北地区地震资料的小波分析,J,.大地测量与地球动力学，23(4):28-33. 万永革.2007.数字信号处理的matlab实现,M,.北京:科学出版社. 夏浩明.1987.大华北M?4.0地震阶级性分区特征及其意义,J,.地球物理学报，30(3):281-291. 徐义贤，王家映.2000.基于连续小波变换的大地电磁信号谱估计 方法,J,.地球物理学报，43(5):677,683. 张国民.1987.我国大陆强震活动的韵律特征,J,.地震地质，9 (2):27-37. 张国民，张培震.2000.大陆强震机理与预测中期学术进展,J,. 中国基础科学，(10):4-10. 张培震，邓起东，张国民，等.2003.中国大陆的强震活动与活动 地块,J,.中国科学(D辑)，33(增刊):12-20. 中国地震局震害防御司.1995.中国历史强震目录(公元前23世 纪—公元1911年),M,.北京:地震出版社. Benioff H.1951.Crustal strain characteristics derived from earthquake sequences,J,.Trans Amer Geophy Union，32(4): 203-208. Logan B.1977.Information in the zero-crossirig of bandpass signals,J,.Bell Syst Tech J，56，510-532. Mallat S.1991.Zero-crossing of a wavelet transform,J,.IEEE Trans，IT-37(1):1019-1033. Mallat S，Hwang W L .1992.Singularity detection and processing with wavelets,J,.IEEE Trans AMS，IT-38:617-643. Research on the Characteristics of Seismic Activity in North China Using Wavelet Transform ZHANG Xiaotao1，2，FENG Xiangdong3，ZHANG Yuegang3， ZHANG Shuangfeng1，YANG Yaqiong3，HAN Liping1 (1.Central Seismic station of Handan，Earthquake Administration of Hebei Province，Handan 056001，Hebei，China) (2.Jilin University，Changchun 130026，Jilin，China) (3.Earthquake Administration of Hebei Province，Shijiazhuang 050021，Hebei，China) Abstract:Referring to the M?5.0 earthquakes from 1500 to 2000 and the M?4.0 earthquakes from 1970 to 2009 in North China，we analyze the Benioff strain using wavelet transform in various scales.Analysis shows that on a centuryscale the M?5.0 earthquakeactivity period has just finished and now the stress is accumulating.And since 1970 there has existed activity cycle of about 7 years for the M?4.0 earthquakes;now the M?4.0 significant earthquakeactivity period is beginning in North China. Key words:Benioff strain;wavelet transform;seismic active period;seismic quiet period;North China