摘要
利用希尔伯特-黄变换方法对2013年芦山MS7.0地震前甘肃省平凉台、四川省成都台和云南省元谋台地电场观测资料进了分析研究。为增强对研究结果的说服力,还对平凉、成都和元谋地电场观测台站在正常时段变化情况,与地震发生前夕的处理结果进行了对比,对成都地电场观测台站在芦山MS7.0和汶川MS8.0地震前的变化情况也进行了对比分析,并得出:1希尔伯特-黄变换方法能够对地电场观测资料进行有效处理,处理结果简单直观;2在芦山地震前1~2个月地电场时频图出现能量增强现象,出现变化的周期范围约为90~1400h,能量最小变化幅度有3个数量级,最大变化幅度达到6个数量级;3通过上述三个地电场观测结果的对比,以及成都地电场观测结果的类比方式,大大提高了分析研究结果的可信度。在讨论部分,对地电场观测中能够观测到震前异常变化的可能性,从实验和数值模拟两个方面进行了相关论述。
In this paper,Hilbert-Huang Transformation method is introduced to deal with geo-electric field data observed by the stations around the epicenter of 2013 Lushan MS7.0earthquake.These stations include the stations of Gansu Pingliang,Sichuan Chengdu and Yunnan Yuanmou.We compare the results of the three geo-electric field stations between the normal time with the eve of the Lushan MS7.0earthquake,the results of Chengdu geoelectric field station about the Lushan Ms7.0and Wenchuan Ms8.0earthquakes are also analyzed.We conclude that:(1)HHT method is effective for the data analysis observed by the geo-electric field station,the time-frequency spectrum may be simple and intuitive;(2)About 1~2 months before the Lushan MS7.0earthquake,the energy of geo-electric field enhanced,with period range of approximately 90~1400hours,the minimum amplitude of energy variation is about 3orders of magnitude,the largest ranges up to about 6orders;(3)The credibility of results obtained are greatly enhanced through comparison obtained from the three geo-electric field stations,and analogy obtained from the Chengdu geo-electric field station.In the discussion section of results,we discussed the?situation of seismic electromagnetic signals based on the experimental and numerical simulation method in the domestic and abroad in recent years,these results obtained by them can provide theoretical support to our results in this paper.
引文
[1]Varotsos P,Alexopoulous K,Nomicos K.Seismic electric currents[J].Prakt Akad Athenon,1981,56:277-286.
[2]Varotsos P,Alexopoulous K.Physical properties of the variations of the electric field of the earth preceding earthquakes,I[J].Tectophysics,1984a,110:73-98.
[3]黄清华.地震电磁观测研究简述[J].国际地震动态,2005,323(11):2-5.
[4]汤吉,詹艳,王立凤,等.汶川地震强余震的电磁同震效应[J].地球物理学报,2010,53(3):526-634.
[5]傅承义,陈运泰,祁贵仲.地球物理学基础[M].北京:科学出版社,1985,447.
[6]徐文耀.Sq发电机电流的逐日变化和Sq指数[J].地球物理学报,1992,35(6):676-683.
[7]黄清华,刘涛.新岛台地电场的潮汐响应与地震[J].地球物理学报,2006,49(6):1 745-1 754.
[8]谭大诚,赵家骝,席继楼,等.潮汐地电场特征及机理研究[J].地球物理学报,2010,53(3):544-555.
[9]范莹莹,杜学彬,Jacques Zlotnicki,等.汶川MS8.0大震前的电磁现象[J].地球物理学报,2010,53(12):2 887-2 898.
[10]谭大诚,赵家骝,席继楼,等.青藏高原中强地震前的地电场变异及构成解析[J].地球物理学报,2012,55(3):875-885.
[11]田山,张磊,王建国,等.汶川、玉树大地震前的地电场异常[J].地球物理学进展,2012,27(3):0 878-0 887.
[12]马钦忠,冯志生,宋治平,等.崇明与南京台震前地电场变化异常分析[J].地震学报,2004,26(3):304-312.
[13]田山,王建国,徐学恭,等.大地电场观测地震前兆异常提取技术研究[J].地震学报,2009,31(4):424-431.
[14]钱家栋,顾左文,赵家骝,等.《地震台站观测环境技术要求第2部分:电磁观测》(GB/T 1953.2-2004).北京:地震出版社,2004.
[15]杜学彬,席继楼,谭大诚,等.《地震台站建设规范地电台站第2部分:地电场台站》(DB/T 18.2–2006).北京:地震出版社,2006.
[16]石春香,罗奇峰.时程信号的Hilbert-Huang变换与小波分析[J].地震学报,2003,25(4):398-405.
[17]杨培杰,印兴耀,张广智.希尔伯特-黄变换地震信号时频分析与属性提取[J].地球物理学进展,2007,22(5):1 585-1 590.
[18]张立,傅容珊,周挚,等.基于HHT提取重力固体潮的地震前兆信息[J].地震学报,2007,29(3):222-226.
[19]Huang N E,Zheng Shen,Long S R,et al.The empirical mode decomposition and the Hilbert spectrum for non-liner and non-stationary time series analyses[J].Proc R Sco Lond,1998,545:903-995.
[20]Huang N E,Zheng Shen,Steven R L.A new view of nolinear water waves:the Hilbert spectrum[J].Annu Rev Fluid Mech,1999,31:417-457.
[21]Huang N E,Wu M C,Long S R,et al.A confidence limit for the empirical mode decomposition and Hilbert spectral analyses[J].Proceeding of the Royal Society of London,2003,459:2 317-2 345.
[22]Rilling G,Flandrin P,Gonalves P,et al.“Bivariate Empirical Mode Decomposition”,Signal Processing Letters[J].IEEE,2007,14(12):936-939.
[23]安张辉,元丽华,李宁,等.HHT方法在地电场数据处理中的应用[J].地球物理学进展,2010,25(2):525-532.
[24]安张辉,杜学彬,元丽华,等.HHT方法在受城市轨道交通干扰地电场观测数据中的应用[J].地震学报,2011,33(2):243-251.
[25]高永新,胡恒山.水平分层孔隙介质中点源激发的震电波场的数值模拟及分析[J].地球物理学报,2009,52(8):2 093-2 104.
[26]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].Geophys J Int,2012,188(3):925-944.
[27]Ikeya M,Takaki S,Matsumoto H,et al.Pulsed charge model of fault behavior producing seismic electric signals(SES)[J].J Circuits Syst Comput,1997,7:153-164.
[28]Huang Q H.Rethinking earthquake-related DC-ULF electromagnetic phenomena:towards a physics-based approach[J].Nat Hazards Earth Syst Sci,2011,11:2 941-2 949.
[29]Varotsos P.The physics of seismic electric signals[M].Terra Pub,Tokyo,2005.
[30]黄清华,林玉峰.地震电信号选择性数值模拟及可能影响因素[J].地球物理学报,2010,53(3):535-543.
[31]Huang Q H,Lin Y F.Selectivity of seismic electric signal(SES)of the 2000Izu earthquake swarm:A 3-D FEM numerical simulation model[J].Proc Japan Acad,2010a,86:257-264.