背景噪声中地震信号的可识别判据及量化验算
2023-11-01
来源:乌哈旅游
第34卷 第s/6期 2013年 12月 地震地磁观测与研究 SEISM0L0GICAL AND GE0MAGNETIC OBSERVAT10N AND RESEARCH Vo1.34 No.5/6 Dec. 2013 doi:10.3969/j.issn.1003—3246.2013.os/o6.024 背景噪声中地震信号的可识别 判据及量化验算 崔庆谷 樊丽云 李万金 (中国昆明650224云南省地震局) 摘要 背景噪声与地震信号的强度或周期差异决定地震震相是否可从视觉上有效识别。为寻找 平稳背景噪声中地震信号的视觉可识别判据,研究能反映两者频率差异及强度差异的参量用于上 述两种信号对比,得到背景噪声中地震信号视觉可识别条件的量化判据,并用实际观测资料及仿 真数据进行验证。在此基础上,针对功率谱计算过程中的校正环节进行讨论,利用Albuquerque实 验室给出的测试实例阐明USGS标准中功率谱计算的规范步骤及细节。 关键词 背景噪声;地震信号;视觉可识别判据;振幅谱密度;功率谱密度 0 概述 背景噪声、地震信号、正弦信号是地震观测中经常面对且必须就某一特征进行比较的信 号。在地震定位和大震速报中,从背景噪声中识别地震信号并标注震相是分析人员的日常工 作,只有当地震波强度大于背景噪声或者地震波的周期与背景噪声相比存在较大差异时,分析 人员才能从背景噪声中识别震相并进行标注。因此,震相识别的过程就是对地震波和背景噪 声的强度(或周期)进行对比并做出判断的过程。长期以来,这个过程一直依靠人的视觉完成, 过程中隐含对地震信号视觉可识别条件的模糊判据。随着数字化技术的推广应用,各种数据 处理方法广泛用于震相的自动识别和定位(沈萍,2002;李山友,2006),视觉可识别已经不再是 地震震相识别的必要条件。但原始的、未经处理的地震记录数据依然是地震工作者在第一时 间直接面对和处理的对象,有必要对地震信号的视觉可识别条件做进一步研究。 本文在理论分析的基础上得到功率谱密度与振幅谱密度之间的转换关系,并结合USGS 测试标准(Charles R Hutt,1999),用实例展示功率谱密度与振幅谱密度量化计算中的细节问 题。在此基础上,研究背景噪声中地震信号的可识别判据,并用真实背景噪声及地震信号为例 加以证实。 量化对比关系的建立 能量信号和功率信号的强度对比不能直接根据振幅谱密度、功率谱密度的数值大小得 结论,而是必须建立两者之间的量化关系。由于地震信号的强度指标是振幅谱密度(安艺敬 作者简介:崔庆谷(1967一),男,高级工程师,2003年中国地震局地球物理研究所固体地球物理专业毕业,获理学博 士学位,曾从事地震计设计及地震计性能研究工作,现从事地震监测工作 基金项目:云南省自然科学基金面上项目(2008CD220)资助 本文收到日期:2013 04 1 7 第5/6期 崔庆谷等:背景噪声中地震信号的可识别判据及量化验算 145 一,1987;彼得・鲍曼,2006)。如果将功率谱密度转换为对应的振幅谱密度就可以将背景噪声 持续时间为T的离散地震信号X(kAt)可以用傅里叶级数表示(Jens Havskov,2007) 与地震信号进行对比。 x(kAt)一∑F ×ejX2 ̄rXnXk/N (1) 其中,N为持续时间T内的采样点数(对于分段DFT,分段重叠部分必须重复计人N),k一0, 1,2,3,…,T/At,N—T/At,At为采样时间间隔,F 为傅里叶谱系数。如果只考虑正频率, A 一F 称为地震信号的振幅谱 A 兰 ∑ (尼△£)8 。积 为消除时间窗长度T对谱的影响,引入振幅谱密度A: 一(2) T×A 一At∑x(kAt)e (3) 令FDF 一∑x(kAt)e 川 为z(kAt)的离散傅里叶变换,则 A 一At×F FT (4) 对于持续时间无限长的背景噪声,在实际计算中必须首先进行截断,假设截断长度等于地震信 号的持续时间T,采样间隔△ 与地震信号相同。则背景噪声信号z ( )截断后的表达式为 xT(t)一 ㈤ 其功率谱密度为(Jens Havskov,2007) P 一 ×I F .r I z (6) 考虑到样本系列的个数N—T×F ,其中F 为采样率,则式(6)变为 P 一—N 一× F× I FDFT” I (7) 对比式(4)、式(7),可以得到功率谱密度与振幅谱密度的转换关系为 A 一√ 截断长度的差异而改变。 (8) 在没有近场强干扰的情况下,背景噪声是平稳随机信号,可任意截取一段有限长度的背景噪声 进行计算。由于式(7)计算过程中已经扣除截断长度的影响,得到的功率谱密度数值不会因为 2转换关系的检验 为了检验功率谱密度与振幅谱密度的转换关系式(8),分别用式(4)、式(8)计算一段人工 合成信号sin(27rX 8Ox )+2sin(2zX140x )+randn(size(t))的振幅谱密度,其中t为时间, randn(size(t))39持续时间为 、方差为1.0的伪随机白噪声,假设该信号的单位为ITI,计算结 果见图1。由图可见,由功率谱密度转换得到的振幅谱密度与按定义计算的结果一致,证明上 述转换方法是正确的。 地震地磁观测与研究 3 USGS标准中功率谱密度计算校正环节 根据功率谱密度的定义式(7),背景噪声信号的离散傅里叶变换F FT必须经过采样率F 校正、窗长度(样本数N)校正才能得到功率谱密度。在实际计算过程中,功率谱估计的方法 很多,不同资料所采用的校正环节不尽相同、得到的计算结果也存在差异,由此造成数据交换 和成果共享的困难,有必要用通行的计算标准进行统一。 此结果证明,Peterson算法中已经做了窗长度及采样率校正,与式(7)的定义相同,是标准 Welch周期图方法,可以用Matlab平台中的Welch函数实现。 -r.N l10 10‘l0 10。” 1O \量 船10n 按 0 5 i0 15 2O 25 30 35 40 45 50 f/H 工 10 I 5 图1 功率谱密度和振幅谱密度转换关系检验 Fig.1 Diagram for comparing the results of am plitude spectral density from two arithme一 10 tic form as 餐 / /… / ; . . __ 、 。 4地震信号可识别条件模型检验 4.1微小地方震信号识别 以云南省易门地震台2009年5月16日18点记录的真实地噪声和小地震事件为研究对 象,分别计算地震事件的振幅谱密度A 及背景噪声振幅谱密度A z,结果见图3(a);两种信号 的振幅谱密度比为R—A /A ,见图3(b);实际记录数据见图3(c)。 从图3(a)中可以看出,该地震事件的振幅谱密度A 在O一2.5 Hz、4.8 9.6 Hz、10.5 13.6 Hz期间的值明显大于地背景噪声振幅谱密度A 。其中,4.8—9.6 Hz之间的S波强度 高于背景噪声强度,容易从背景噪声中识别。而1O.5 13.6 Hz频段的P波及O 2.5 Hz频 段的尾波强度稍弱,但依然可以从背景噪声中识别出来[图3(c)]。进一步计算可见,在R>2 时,地震波能够从背景噪声中识别出来,见图3(b)。 第5/6期 崔庆谷等:背景噪声中地震信号的可识别判据及量化验算 147 3 姜 墨2 × 嚣 娄。 蝼 _们6 考3 。 一3 1171 1173 1175 t/s ll77 1179 1181 图3 易门台记录地震波及地噪声的振幅谱密度(a)、振幅谱密度比R(b)及真实记录时程(c) Fig.3 Diagrams of Amplitude spectral density of the background noise and the seismic event get from Yimen station,Yunnan province(a),ratio of the tWO Amplitude spectral density R (b)。and the record data(c) 为进一步证实振幅谱密度比作为信号可识别指标的有效性,在地震波强度保持不变时,将 图3中地背景噪声强度提高到原噪声的2.0和4.0倍,再次计算两种信号的振幅谱密度比,结 果见图4(a),实际记录数据见图4中(b)、(c)、(d)图。 从图4(a)中可以看出,随着背景噪声强度增大,地震事件与背景噪声的振幅谱密度比R 逐渐减小。对比图4(a)、(c),当背景噪声强度增加2.0倍时,P波已经无法识别,此时P波对 应频段(10.5—13.6 Hz)的谱密度比R≤2.0,但S波频段(4.8—9.6 Hz)对应的谱密度比R 的峰值依然大于2.0,对应记录中的S波依然可以识别。 对比图4(a)、图4(d),当背景噪声强度增加为4.0倍时,P波、S波均无法识别,此时谱密 度比R值小于2.0。 6 4 们l 2 0 砉5 0 5 O ,肫 10 15 酱一s 专5 0 簧一s 图4地方小震与背景噪声振幅谱密度比R随背景噪声强度倍数(1.0、2.0、4.0倍)的变化曲线(a)、 1.0倍背景噪声时的记录时程(b)、2.0倍背景噪声时的记录时程(c)、4.0倍背景噪声时的记录时 程(d) Fig.4 Amplitude spectral density ratios of local evenbto background noise varying with strength(10. .2.0,4.0 multiple)of the background noise 4.2近震信号识别 以云南省勐腊地震台2010年12月5日21时记录的背景噪声和近震事件(S—P- ̄300 s) 为研究对象,分别计算地震事件的振幅谱密度与1.0倍、4.0倍背景噪声振幅谱密度比,得到2 148 地震地磁观测与研究 条振幅谱密度比R随频率变化曲线[图5(a)]、实际记录数据见图5(b)、图5(c)。 1 8 Q=1 O 2 10 10 2 砉1 一 交。 1 1080 11O0 l120 f/Hz 量 l140 t/s l160 l180 1200 0 一 厂 小 Z 倒 80 I100 112O 1140 t/s 1160 l18O 1200 O 图5近震与背景噪声振幅谱密度比R随背景噪声强度倍数(1.0、4.0倍)的变化曲线(a)、1.0倍背景 噪声时的记录数据(b)、4.0倍背景噪声时的记录数据(c) Fig.5 Amplitude spectral density ratios of near—field event—to—background noise varying with strength (1.0.2.0.4.0 multiple)of the background noise 从图5(a)可见,振幅谱密度比R存在2个峰值区,分别代表S波和P波所在频段,随着背景 噪声强度增大,事件波形与背景噪声的振幅谱密度比R逐渐减小。对比图5(a)、图5(c),当背景 噪声强度增加为4.0倍时,P波无法识别。P波对应频段(2.0 5.0 Hz)的谱密度比R≤2.0,但S 波频段(1.o 2.0 Hz)对应的谱密度比R依然大于2.0,对应记录中的s波依然可以识别。 4.3远震信号识别 以云南省勐腊地震台2010年12月1日05时记录的真实地噪声和远震事件为研究对象, 分别计算地震事件的振幅谱密度与1.0倍、4.0倍、8.5倍背景噪声的振幅谱密度比,得到3条 振幅谱密度比R随频率变化曲线[图6(a)],实际记录数据见图6(b)、图6(c)、图6(d)。 :4 昌 0 世 一4 一 10 8 曼 2 墨4 0 .遗 蚓0 d d 230O 2 7u0 t/s 3100 3500 2300 2700 Us 3100 350。 图6远震波与背景噪声振幅谱密度比R随背景噪声强度倍数(1.0、4.0倍)的变化曲线(a)、3.0倍 背景噪声时的记录数据(b)、8.5倍背景噪声时的记录数据(c) Fig.6 Amplitude spectral density ratios of tele—seismic event—to—background noise varying with strength(1.0,2.0,4.0 multiple)of the background noise 从图6(a)可见,该远震对应的振幅谱密度比R存在2个峰值区,分别代表体波(1.0— 2.0 Hz)和Rlayleigh面波(5.0 90 s)所在频段,随着背景噪声强度增大,振幅谱密度比R逐 渐减小。对比图6(a)、图6(d),当背景噪声强度增加为8.5倍时,体波已经无法识别,体波对 应频段(1.O 2.0 Hz)的谱密度比R≤2.0,但面波频段的谱密度比R依然大于2.0,对应记录 中的面波依然可以有效识别。 综合以上分析结果,选取更多震例进行计算,结果表明:无论对地方震、近震或远震,地震 第5/6期 崔庆谷等:背景噪声中地震信号的可识别判据及量化验算 149 波与背景噪声振幅谱密度比R作为地震信号的视觉可识别判据是有效的,当R≥2.0时,对应 频段信粤可以从背景噪声中有效识别;当R<2.0时,对应频段信号难以肉眼识别。 5 结论 对于被噪声所包围的特定地震波图像,虽然能够依靠视觉对地震波是否可识别做出判断, 但难以对具体的判断条件进行精确描述。虽然频域RMS一直作为地震波可识别判据用于地 震台站和台网的监测能力计算,但由于RMS值来源于平稳随机噪声的功率谱,代表某一频带 宽度(可变)内平稳随机信号的平均强度,因此将RMS值作为地震信号的可识别条件与部分 通行标准(安艺敬一,1987;彼得・鲍曼,2006)不一致。 振幅谱密度能够反映两种信号的强度差异,也能够以强度随频率变化的方式反映两者的 周期差异,该参数~直作为地震波强度的参数出现在教材中(安艺敬一,1987;彼得・鲍曼, 2006)。利用功率谱密度与振幅谱密度的转换关系,可以实现地震波与背景噪声的量化对比。 实际和模拟的记录数据证明,振幅谱密度比是背景噪声中地震信号可识别的一种可靠判据。 参考文献 李山有,朱海燕,武东坡,等.基于振动和瞬时频率的震相自动识别方法[J].世界地震工程,2006,22(4):1—4. 沈萍,郑治真,刘希强,等.小震综合识别研究EJ].地震学报,2002,24(2):169—175. 安艺敬一,P G理查兹.定量地震学——理论和方法・第一卷[M].北京:地震出版社,1987:552—554. 彼得・鲍曼.新地震观测实践手册第一卷[M].北京:地震出版社,2006:94一100. Charles R Hutt.Standards for Seismometer Testing(A Progress Report)l D 1.U S Geological Survey,Albuquerque Seismo— logical Laboratory,Kirtland AFB-East,Albuquerque,New Mexico,USA,1990:8—10. Jens Havskov,Gerardo Alguacil(著),赵仲和赵建和(译)、地震观测技术与仪器[M].北京:地震出版社,2007:118—123. Peter W Rodgers.Frequency Limits For Seismometers As Determined From Signal—To—Noise Ratios.Part 1.The Electro— magnetic Seismometer J I.BSSA,1992,82(2):1 071—1 098. Peter W Rodgers.Frequency Limits For Seismometers As Determined Ftom Signal—To—Noise Ratios.Part 2.The Feedback Seismometer l J 1.BSSA,1992,82(2):1 099—1 123. Peterson J.Observations and modeling of Seismic background noise[M].U S Geol Survey Open-File Report,1993:93—320. Criterion for distinguishing seismic signal from steady background noise visually:its quantificatiOn and verification Cui Qinggu,Fan Liyun and Li Wanjin (Earthquake Administration of Yunnan Province,Kunming 650224,China) Abstract Small seismic event waves and weak teleseismic signals are always buried in background noise and cannot be identified visually;the difference both in strength and frequency between background noise and seismic signal decides whether the seismic signal can be distinguished from steady background noise.For searching a parameter which can indicate the difference both in strength and frequency between background noise and seismic signa1.a criterion for distinguishing seismic signa1 from steady background noise visually are worked out and then demonstrated by practica1 data and its simulation.Furthermore,which kind of amending steps should be used in calculation of power spectral density are argued and then the normal way of USGS standard for power spectral density calculation is detailed and demonstrated bv an example of Albuquerque seismological laboratory. Key words:background noise,seismic signal,visual distinguishable criterion,amplitude spectral density,power spectral density