基于Wigner-Ville分布的地磁前兆异常提取方法及其在川滇地区的应用

地震 ›› 2014, Vol. 34 ›› Issue (4): 88-99.

基于Wigner-Ville分布的地磁前兆异常提取方法及其在川滇地区的应用

何康¹, 晏锐², 王春³, 郑海刚¹, 李军辉¹, 方震¹

1.安徽省地震预报研究中心,安徽合肥230031;
2.中国地震台网中心,北京100045;
3.黑龙江省地震局,黑龙江哈尔滨150090

收稿日期:2014-03-11 修回日期:2014-07-23 发布日期:2020-09-11
作者简介:何康(1982 ),男,湖北大悟人,硕士研究生毕业,主要从事地震电磁学和地震预报等研究。
基金资助:
“十二五”国家科技支撑项目( 2012BAK19B02- 03)和中国地震局2013年度震情跟踪青年课题(20 13020403)资助

Extracting Geomagnetic Precursory Anomaly based on Wigner-Ville Distribution and Application in the Sichuan- Yunnan Region

HE Kang¹, YAN Rui², WANG Chun³, ZHENG Hai-gang¹, LI Jun-hui¹, FANG Zhen¹

1. Earthquake Administration of Anhui province,Hefei 230031,China;
2. China Earthquake Networks Center, Bejjing 100045,China;
3. Earthquake Administration of Heilongjiang province, Harbin 150090,China

Received:2014-03-11 Revised:2014-07-23 Published:2020-09-11

摘要/Abstract

摘要： 为了能够从强背景场中提取较弱的震磁信息,采用平滑伪Wigner-Ville 分布方法对川滇地区的地磁场Z分量数据开展时频分析。通过研究Z分量滤波后主要周期为2~6 h部分中的不同周期分量变化,发现其中2.4h的周期分量受到外源场影响较大,而4.4h的周期分量在多次显著地震前均出现明显的高值异常。

关键词: 地磁Z分量, Wigner-Ville分布, 前兆异常, 川滇地区, 芦山地震

Abstract: To extract weak seismo-magnetic information from strong background field, we carry out time- frequency analysis of geomagnetic Z-component in Sichuan—Yunnan region using Smoothed and Pseudo Wigner- Ville Distribution( SPWVD). By studying different periodic components of the filtered Z-component , which main periods are 2~6 h,we find that the 2.4 h periodic component is easy to be affected by external field, and the 4.4 h periodic component shows obvious high value of anomalies before many significant earth- quakes in the region.

Key words: Geomagnetic Z-component, Wigner-Ville distribution, Precursory anomaly, Sichuan-Yunnan region, Lushan earthquake

中图分类号:

P315.7

何康, 晏锐, 王春, 郑海刚, 李军辉, 方震. 基于Wigner-Ville分布的地磁前兆异常提取方法及其在川滇地区的应用[J]. 地震, 2014, 34(4): 88-99.

HE Kang, YAN Rui, WANG Chun, ZHENG Hai-gang, LI Jun-hui, FANG Zhen. Extracting Geomagnetic Precursory Anomaly based on Wigner-Ville Distribution and Application in the Sichuan- Yunnan Region[J]. EARTHQUAKE, 2014, 34(4): 88-99.

参考文献

[1] 徐文耀. 地磁学[M]. 北京:地震出版社,2003.
[2] 陈斌、自然正交分量方法在地震地磁监测中的应用[J].地震研究,2011,34(4):466-469.
[3] 田文通,辛长江,薛锋,等.甘肃不同台站地磁垂直分量变化特征分析[J]. 防灾科技学院学报,2011,13(4): 35-40.
[4] 何宇飞,杨冬梅,朱荣,等.四分位法在地磁日变化观测数据分析中的应用[J].地震地磁观测与研究,2011, 32(1): 65-70.
[5] DunlopIS, Menk F W, Hansen H J, et al. A muhistations study of long period geomagnetic pulsations at cusp and boundary layers atitudes[J]. J Atmos Terr Phys, 1994,56: 667-679.
[6] Kivelson M G,Russell C T.Introduction to Space Physics[ M]. London: Cambridge University Press. 1965.
[7] 胡昌华,张军波,夏军,等.基于Matlab的系统分析与设计一 -小波分析[M].西安:西安电子科技大学出版社,1999.
[8] 刘本永. 非平稳信号分析导论[M].北京:国防工业出版社,2006.
[9] Cohen L.Time frequency analysis[J]. IEEE Signal Processing, 1996,16(1): 22-28.
[10] 张贤达,保铮.非平稳信号分析与处理[M].北京:国防工业出版社,1998.
[11] Wigner E.On the quantum correction for thermodynamic equilibrium[J]. Phys Rev, 1932, 40:749-759.
[12] Ville J.Theorie et application de la notion de signal analytique[J] . Cables et T ransmission, 1948,2A: 61-74.
[13] T ACM Claasen, W FG Mecklenbrauker.The Wigner distribution-a tool for time frequency signal analysis, Part IT: discrete time signals[J]. Phillips J Res, 1980, 35: 276-300.
[14] Matz G,Hlawatsch. W igner distributions( nearly) everywhere: time- frequency analysis of signals,systems, random processes,signal spaces and frames[J]. Signal Processing, 2003, 83(7): 1355-1 378.
[15] J B Wu,J Chen, P Zhong.Time- frequency- -based source separation technique for elimination of cross- terms in Wigner distribution[J]. IEE Electronics, 2003, 39(5): 475-477.
[16] 张梅军,李曙光,孙启亮.四种典型时频分析方法研究与仿真[].仪器仪表与分析监测,2010, 3:25-27.
[17] 潘建华. SPWD在脉内混合信号时频分析中的应用[J].舰船电子对抗[J]. 2005, 28(5): 41-44.
[18] 吕久明,罗景青. PWVD在调频信号时频分析中的应用[J].弹道学报[J]. 2007. 19(2): 93-96.
[19] 张宁,林春生,庞学亮,等.复杂海况下舰船微弱信号的时频辨识研究[J]. 兵工学报,2009,30(6): 834-838.
[20] 夏洪瑞, 张怿平,张乃勋. FIR与PWVD串联滤波消除随机噪声[J].石油地球物理勘探,2009, 44(1): 33-37.
[21] 候树杰、时频分析方法对比及S变换在地震数据处理中的应用[J]. 油气地质与采收率,2011, 18(3): 44-46.
[22] 何康,晏锐,郑海刚,等,伪魏格纳一维勒分布在地磁时频分析中的应用[J]. 中国地震,2013, 29(1): 157-167.
[23] 杜斌,安明智,林洋,等.磁通门磁力仪在野外安装的实施方案[J]. 地震地磁观测与研究,2008,29(6): 154-160.
[24] 王晓美,滕云田,王喜珍,等. GM4型磁通门磁力仪性能检验方法[J].地震地磁观测与研究,2008,29(6): 88-94.
[25] 吴琼,王喜珍,杨冬梅,等. FHDZ- M15地磁组合观测系统OVERHAUSER探头对FGE探头的影响初探[J].地震地磁观测与研究,2008, 29(6): 161-169.
[26] 葛哲学, 陈仲生. MATLAB时频分析技术及其应用[M].北京:人民邮电出版社,2006.
[27] WuX Y, Liu T Y. Time frequency analysis on wigner- vlle distribution of seismie signal based on time-frequency rearrangement[J]. Oil Geophysical Prospecting, 2009, 44: 86-91.
[28] 康国发. 地磁场H和Z分量月均值的频谱分析[J].云南大学学报(自然科学版),1990, 12(1):88-92.
[29] 蒋骏,李胜乐,张雁滨,等.地震前兆信息处理与软件系统[M]. 北京:地震出版社,2000.
[30] 詹志佳, 高金田.频谱分析在地磁与震磁研究中的应用[J].地震地磁观测与研究,1992,13(4):30-37.
[31] 徐文耀. 地球电磁现象物理学[M].合肥:中国科学技术大学出版社,2009.
[32] 林云芳. 我国震磁关系研究的新进展[].地震,1990,6: 26-34.
[33] 徐文耀. 中国地磁学研究的进展[J].地球物理学报,1994, 37(增刊): 284-294.
[34] Huang Q H.One possible generation mec hanism of co- seismic electric signals[J]. roceedings of the Japan Academy, 2002, 78(7B): 173-178.
[35] 马占虎,杜学彬,谭大诚.地震电性变化及其物理机理初步讨论[J]. 西北地震学报,2004, 26(3):234-239.
[36] 陆阳泉,钱家栋,刘建毅.大型花岗岩标本缓慢膨胀破裂过程中电阻率和声发射前兆特征的实验研究[J].西北地震学报,1990,12(2): 35-41.
[37] 刘煜洲,刘因,王寅生,等.岩石破裂时电磁辐射的影响因素和机理[J].地震学报,1997, 19(4):418-425.
[38] Ren H X, Chen X F, Huang Q H.Numerical simulation of coseismic electromagnetic fields associated with seismic waves due to finite faulting in porous media[J]. GeophysJ Int, 2012,188: 925-944.
[39] 赵超,申俊峰,张自力,等,含磁铁矿岩石热电实验研究及其地震地质意义[J]. 矿物岩石,2012,32(1): 16-20.
[40] 钱家栋,陈有发,等.地震地电学译文集[M].北京:地震出版社,1989.
[41] 国家地震局预测预防司.电磁学分析预报方法[M].北京:地震出版社,1998.
[42] 杨涛,刘庆生,付媛媛,等.震磁效应研究及进展[J]. 地震地磁观测与研究,2004, 25(6): 63-71.
[43] 李延峰,朱胜伟,杨广华、格尔木地磁Z分量静日变化特征与前兆异常分析[J].高原地震,2013, 25(1): 35-41.
[44] 曾小苹,郑吉盎,王墨談,等.震前特大地磁异常及其短临预警意义[J].中国工程科学,2011, 4: 48-53.
[45] 胡久常,刘伟,郭敏瑞,等.汶川8.0级特大地震前的地磁垂直分量日变“双低点”异常[J]. 地震学报,2009, 31(5): 589-593.
[46] 王武星,丁鉴海,余素荣,等.汶川Ms8.0地震前地磁短临异常与强震预测探索[J].地震学报,2009,31(2): 172-179.