A Wavelet Comprehensive Threshold Seismic Signal Denoising Method
-
摘要: 小波阈值方法中硬、软阈值方法是地震信号降噪常用方法,但容易造成信号中高频信息丢失导致地震误判和漏判情况发生。小波综合阈值方法继承和发展了硬、软阈值降噪方法的优点,对信号高频部分用硬阈值方法,以提高高频信号能量,对信号低频部分用软阈值方法,提高信号降噪能力的同时保证信号连续性和光滑性。利用噪声信号小波系数小和地震信号小波系数大的特征,进行雷克子波降噪仿真实验和实际地震信号降噪实验。仿真实验表明,小波综合阈值方法降噪后波形MSE值最小,且降噪后与原信号波形最近似,降噪后波形高频部分能量增强且抑制低频部分能量。最后,对实际采集的地震信号进行降噪处理,处理后信号中能量增强被压制,利用处理后的信号可得到地震的初至时间。Abstract: The wavelet threshold method with hard or soft threshold is common for noise reduction of seismic signal, but it is often to cause the loss of high frequency information in the signal to lead in earthquake misjudgment or missing judgment. Taking the advantages of hard threshold and soft threshold denoising method, the wavelet comprehensive threshold improves the high frequency signal energy by the hard threshold treatment, and keeps the signal smoothness and reduce the signal degradation of noise capacity while ensuring signal continuity by lowness of the signal. Based on the small wavelet coefficients of the noise signal and the large wavelet coefficients of the seismic signal, the simulation experiment of the Ricker signal and the noise reduction experiment is carried out. The simulation results show that the MSE of the wavelet comprehensive threshold is the closest to the original signal waveform after noise reduction, and the energy of the high frequency part of the waveform is reduced and the low frequency part is suppressed. Finally, taking the noise reduction of the natural waveform as an example, we found that noise signal energy is suppressed and the first time of the waveform can be obtained.
-
Key words:
- Comprehensive threshold /
- Seismic signal /
- Wavelet method
-
引言
地震活动时起时伏,具有明显的轮回性,但又不是简单的重复,地震发生的不均匀性表现为在一定时间和区域内相对平静或活跃,反映了地震的孕育和发生是区域内应力长期积累、集中、加强并在应力集中区最终释放的结果,因此地震活动会表现出平静—密集的活动特征(王健,2005)。马宗晋等(1982)根据地震活动的增强和平静提出将百年尺度的地震活动韵律称为地震期,在地震期内还有若干个10—20年尺度的地震活跃幕及间隔10、20年的相对平静幕,同时指出华北地区第4地震活跃期始于1815年,目前仍处于第4活跃期。此后,不同学者(汪良谋,1987;陈荣华,1989;朱红彬等,2010)采用不同划分方法对华北地区的地震幕进行了研究,王亚茹等(2014)将华北地区第4活跃期划分为5个地震幕,每个相对平静幕最短16年、最长为22年,目前处于相对平静时段,自1998年张北6.2级地震以来,6级地震平静达到了20年,有随时进入新的地震活跃幕的可能。
首都圈地区是我国的政治、经济、文化中心,在历史上发生多次强烈地震,尤其是1978年唐山大地震和1998年张北地震,都造成了巨大的损失,地震灾害的破坏性和突发性严重威胁着首都圈地区的平稳发展。首都圈地区的地震活动向来受到地震学家的高度重视,很多学者对该区域的地震活动进行了大量的详细分析(刘爱文等,2006;张世民等,2006;胡幸平等,2013)。自2006年文安5.1级地震后,该区域5级地震平静期达到了12年,该区未来发生5级以上地震危险性判断将成为研究的重点。本文对首都圈地区5级以上地震活动进行了研究,分析了历史地震的时空演化特征,探讨了首都圈5级以上地震的发震特点与迁移规律,对首都圈地区的中强地震预测具有一定的意义。
1. 区域背景与资料选取
1.1 区域地质构造背景
首都圈地区(38.5°—41.5°N,114°—120°E)位于华北地块的北部,即张家口-渤海地震带、华北平原地震带和山西地震带的交汇部位,地震活动比较频繁,震中分布范围广,震害严重,据史料记载该区历史上曾发生一系列强震活动,如1679年三河-平谷8级地震以及1976年唐山7.8级地震。该区主要发育一系列NNE—NE向和NW—EW向断裂,这些断裂多集中成带分布,控制了断陷盆地和坳陷的发育,其中NNE—NE向构造带自西向东依次为延庆-怀来构造带、怀柔-北京-涿县构造带、平谷-三河-廊坊构造带以及沧东构造带,NWW向的张家口-渤海构造带从上述4条NNE—NE向构造带的北段斜穿而过,构成了首都圈新生代显著的构造分块特征的基本格局(高文学等,1993)。总体来说,首都圈地区由于东部受到西太平洋板块俯冲带所产生的侧向影响,西部受到青藏高原隆起所引起的深部塑性流动致使浅部物质侧压挤出的影响,形成了一个主压应力为NEE—SWW向的水平应力场(丰成君,2014)。
1.2 资料选取
华北地区第4活跃期始于1815年(马宗晋等,1982),同时首都圈地区进入长达65年的5级地震平静时期,直至1880年发生河北滦县5级地震,因此本文研究时段自1880年开始。本文所用的历史地震目录合并了由Mapsis下载的破坏性地震目录(MS≥4.7)和全国5级以上地震目录,由于历史地震目录根据史料记载综合整理换算而来,其地震三要素与现今地震计得来的地震三要素不可同日而语,关于余震的记载也将不可避免地缺失,为了保持1880年以来地震活动的可对比性与连贯性,下文的统计分析中剔除了1980年前唐山地震的5级以上强余震活动。
2. 首都圈地区5级以上地震时空特征分析
2.1 首都圈5级以上地震活动的时间特征
自1880年至今约130余年,首都圈地区共发生5级以上地震15次(去除唐山地震的余震),其中5.0—5.9级地震10次,6.0—6.9级地震3次,7.0—7.9级地震2次,最大地震为1976年河北唐山7.8级地震。图 1为首都圈地区地震活动时序图,按照地震间隔时间划分为11个(丛)地震,其中1978年唐山地震的多次5级以上强余震与唐山主震归为1组,1991—1999年发生的4次地震由于间隔时间较短划分为同1组(表 1)。统计表明该区地震活动表现出明显的周期性,前10组地震中有8次地震的时间间隔为10±2年,其中最长的为1888年渤海湾7.5级地震与1911年河北蔚县5.9级地震间的间隔,为23年,最短间隔为8年。该区最近一次5级以上地震为2006年文安地震,距今已有12年,这表明该区目前正处于发生5级以上地震的节点上,值得重点关注。
表 1 首都圈地区MS≥5.0地震活动分区表Table 1. Seismicity zones of MS≥5.0 earthquakes in the capital region序号 主震参数 地震间隔/a 分区 时间 纬度/°N 经度/°E 震级/MS 地点 1 1880-09-06 39.70 118.70 5.0 河北滦县 8 东区 2 1888-06-13 38.50 119.00 7.5 渤海湾 23 东区 3 1911-01-25 39.80 114.50 5.9 河北蔚县 12 西区 4 1923-09-14 39.40 115.80 5.5 河北新城 11 中区 5 1934-10-27 39.90 119.20 5.0 河北抚宁 11 东区 6 1945-09-23 39.50 119.00 6.2 河北滦县 12 东区 7 1957-01-01 40.50 115.50 5.0 河北涿鹿 10 西区 8 1967-03-27 38.51 116.50 6.3 河北河间 9 中区 1967-07-28 40.55 115.55 5.4 河北怀来 西区 9 1976-07-28 39.63 118.18 7.8 河北唐山 15 东区 10 1991-05-30 39.50 118.20 5.1 河北唐山 8 东区 1995-10-06 39.80 118.50 5.0 河北唐山 东区 1998-01-10 41.10 114.30 6.2 河北张北 西区 1999-03-11 41.20 114.60 5.6 河北张北 西区 11 2006-07-04 38.90 116.25 5.1 河北文安 中区 2.2 首都圈5级以上地震活动的空间特征
按照断裂走向及地震分布,沿张家口-渤海地震带可把首都圈地区划分为3个分区(图 2),分别称为首都圈西区、中区和东区。其中,东区主要包括唐山老震区及部分渤海地区的平谷-三河-廊坊构造带以及沧东构造带;中区主要包括怀柔-北京-涿县构造带;西区主要包括延庆-怀来构造带。
图 2 首都圈地区分区示意图(断层邓起东等,2003)Figure 2. The zoning map of the capital region表 1分区统计中每组地震中最大地震所处分区表示该组地震主体活动区域,如1967年首都圈中区河间6.3级地震和首都圈西区的怀来5.4级地震,由于中区的地震震级更大,因此该轮地震的主体活动区域为首都圈中区;同样,1991—1999年先后发生4次5级以上地震,鉴于震级最大的张北6.2级地震位于西区,所以此轮地震主体活动区域为首都圈西区。由此,首都圈地区MS≥5.0地震空间分布可以概况为:东区—东区—西区—中区—东区—东区—西区—中区—东区—西区—中区,按照时间可以简单地分为3个轮回:东区—东区—西区—中区—东区—东区、西区—中区—东区、西区—中区,可概括为西区—中区—东区。地震活动表现出明显的自西向东迁移、再跳转到西区开始一个新轮回的特征,具有迁移性、重复性和轮回性。马宗晋等(1992)研究燕山地震带得到其地震活动存在自西向东迁移的定向性,速率为6—12km/a;蒋铭(1985)在研究华北区域第3、第4地震活动期的地震活动特性时发现,华北地区的地震活动出现由西向东整体迁移的特性。
目前首都圈地震活动已经完成了2轮自西向东的轮回,第1个轮回为1911首都圈西区蔚县5.9级地震—1923年中区新城5.5级地震—1934东区抚宁5.0级地震和1945年滦县6.2级地震,历时35年,最大地震为首都圈东区的滦县6.2级地震;第2个轮回为首都圈西区1957年河北涿鹿5.0级地震—1967年中区河间6.3级地震—1976年东区唐山7.8级地震,历时22年,地震震级逐步升高,最大地震为1978年唐山7.8级地震,同样发生在首都圈东区。在首都圈西区发生1998年张北6.2级地震和首都圈中区2006年文安5.1级地震后,第3个轮回已经迁移了半程。目前,距离2006年文安5.1级地震已经过去了12年,根据首都圈MS≥5地震活动保持以10年左右为一周期的特征,首都圈地区正处于5级以上地震发生的时间节点上;根据首都圈地震MS≥5地震自西向东的迁移特征,下个5级以上地震发生在首都圈东区的可能性较大。另外,自1880年以来,首都圈东区发生MS 5.0—5.9地震4次(去除唐山地震的余震)、MS 6.0—6.9地震2次,分别占首都圈5.0—5.9地震的40%以及6.0—6.9级地震的67%,首都圈2次7级以上地震均发生在东区,且前2个轮回的最大地震均发生在首都圈东区,据此推测地震活动最强烈的主体区域集中在东区的可能性较大。
3. 讨论
上述分析可知,自1880年以来首都圈MS≥5.0地震活动表现出一定的规律,一是周期性,地震活动周期表现为约10年左右发生1次(丛)5级以上地震;二是迁移性,地震活动均表现为自西向东定向迁移的特征,这与 马宗晋等(1992)研究燕山地震带得到地震活动自西向东迁移的结果一致;三是轮回性,具体表现为地震活动迁移到首都圈东区后会跳回首都圈西区开始新的轮回,并且这一定向路径是重复的。1880年至今共经历了2个完整的轮回,目前正处于第3个轮回中。在时间上,目前已经到达发生MS≥5.0地震的发震节点;在空间上,地震发生在东区的可能性较大;在震级上,前2个轮回在东区释放的能量最多,第1、2轮回的最大震级分别达到6.2级和7.8级,因此第3轮回中的最大地震可能达到6级左右。然而,这种定性分析的样本量只有2个完整轮回,对其准确程度的判定需要大量样本,但是由于MS≥5.0地震活动的小概率性和长时间性等特性,应该随着地震形势的发展不断总结其规律,并辅助其他的手段方法从规律再推导其机理。
首都圈地震活动表现出的特征可能与其地质构造环境有一定关系。首都圈位于华北地区,受西太平洋自东向西俯冲挤压的作用,同时受到青藏高原北东向的推挤,太平洋板块对大华北地区的推挤作用占主导地位; 印度板块北东向的推挤穿过青藏高原作用于大华北地区的西南角,起辅助作用;二者控制着大华北的新构造运动,控制地震出现成片迁移现象。王绳祖(王绳祖,1993;王绳祖等,1994)根据大陆岩石圈网络状塑性流动与多层构造变形的观点,认为在印度板块的推挤作用下,亚洲中东部大陆存在着不同准周期的塑性流动波,大陆板块驱动力主要通过岩石圈下层(含下地壳和岩石圈地幔)的网络流动和塑性流动波实现其远程传递,并控制板内构造变形和地震活动,造成地震活动沿塑性流动网带的迁移,因此地震活动表现出一定的周期性和迁移性。
4. 结论
根据首都圈地区1880年以来MS≥5.0地震活动规律的研究发现,该区地震活动存在周期性、迁移性和轮回性的发震特性,具体表现为首都圈地区MS≥5.0地震大概间隔10年左右发生1次(丛),从空间上地震有自西向东的迁移特征,且这种定向迁移具有轮回性,首都圈东区地震活动后会跳回首都圈西区开始新的轮回。2006年文安5.1级地震后,首都圈地区5级地震已平静12年,从时间上达到新的发震节点,根据迁移性发生在首都圈东区的可能性较大,其震级可能会达到6级左右。综上所述,首都圈地区地震活动的周期性、迁移性和轮回性可为该区未来中强地震的危险性判断提供一定参考依据。
-
表 1 噪声分类
Table 1. Noise classification
噪声类型 频率范围/Hz 降噪难度 环境噪声 3—20 难 次生噪声 5—30 难 仪器噪声 1—2 易 表 2 仿真降噪后所得结果的SNR和MSE值
Table 2. SNR and MSE from the de-noised signal in simulation
降噪方法 SNR MSE 硬阈值 18.5736 0.2907 软阈值 19.3285 0.2897 小波综合阈值 19.4136 0.2784 -
耿冠世, 俞言祥, 2015.中国西部地区震源破裂尺度与震级的经验关系.震灾防御技术, 10(1):68-76. doi: 10.11899/zzfy20150107 侯跃伟, 赵兵, 田韬, 2015.基于Daubechies小波分析的南京数字化钻孔形变震前变化特征研究.震灾防御技术, 10(2):388-396. doi: 10.11899/zzfy20150219 孔祥茜, 吴继伟, 岳继光, 2005.地震信号小波变换的去噪方法.计算机辅助工程, 14(3):52-56. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jsjfzgc200503012 范涛, 赵兆, 吴海, 2014.矿井瞬变电磁多匝回线电感影响消除及曲线偏移研究.煤炭学报, 39(05):932-940. http://www.oalib.com/paper/4232079 刘霞, 潘洪屏, 高晓春, 2010.基于小波阈值的地震信号去噪处理.科学技术与工程, 10(29):7251-7254. doi: 10.3969/j.issn.1671-1815.2010.29.030 李英, 张淑贞, 许康生, 2006.小波降噪方法在地震信号处理中的应用.西北地震学报, 28(2):159-162. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xbdzxb200602013 曲中党, 吴蔚, 贺日政等, 2015.基于S变换的软阈值滤波在深地震反射数据处理中的应用.地球物理学报, 58(9):3157-3168. doi: 10.6038/cjg20150912 唐守锋, 童敏明, 潘玉祥等, 2011.煤岩破裂微震信号的小波特征能谱系数分析法.仪器仪表学报, 32(7):1521-1527. http://www.cnki.com.cn/Article/CJFDTOTAL-YQXB201107012.htm 魏学强, 袁洪克, 秦晶晶等, 2016.广义S变换地震信号时频分析.震灾防御技术, 11(4):808-813. doi: 10.11899/zzfy20160411 曾宪伟, 赵卫明, 师海阔等, 2010.利用小波包变换对地震信号进行时频分析时小波基函数的选取.地震研究, 33(4):323-328. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzyj201004012 郑作亚, 卢秀山, 李克行, 2007.一类小波基函数的构造及其在测量数据处理中的应用探讨.测绘科学, 32(3):9-11. http://industry.wanfangdata.com.cn/dl/Detail/Periodical?id=Periodical_chkx200703003 Bruni V., Vitulano D., 2006. Wavelet-based signal de-noising via simple singularities approximation. Signal Processing, 86(4):859-876. doi: 10.1016/j.sigpro.2005.06.017 Mousavi S. M., Langston C. A., 2016. Hybrid seismic denoising using higher-order statistics and improved wavelet block thresholding. Bulletin of the Seismological Society of America, 106(4):1380-1393. doi: 10.1785/0120150345 Xia L. M., Zheng W., Liang M. D., 2017. Application of a comprehensive one dimensional wavelet threshold denoising algorithm in blasting signal. In: Balas V. E., Jain L. C., Zhao X. M., Information Technology and Intelligent Transportation Systems. Cham: Springer, 517-527. -