利用重叠加窗方式恢复地面运动.pdfVIP

地 震 地 磁 观 测 与 研 究 第36 卷 第2 期 SEISMOLOGICAL AND GEOMAGNETIC Vol .36 N o . 2 20 15 年 4 月 OBSERVATION AND RESEARCH Apr. 20 15 doi: 10. 3969/j. issn. 1003-3246. 2015. 02. 007 利用重叠加窗方式恢复地面运动 高 亮 杨维辉 孙 炜 邵 帅 （中国黑龙江 164 155 德都地震台） 摘要 在频率域，利用地震记录反卷积恢复地面运动，由于FFT 运算导致的谱泄漏与系 统传递函数对高频与低频的放大效果会导致地面运动恢复失败，本研究采用重叠加窗方 式，可避免单一加窗的时域失真，并减小频域谱泄漏。结果表明，恢复地面运动取得了 较好效果。 关键词 谱泄漏；地面运动恢复；系统传递函数；重叠加窗；时域 0 引言 对地震观测的基本要求是，准确地记录与研究地震发生时由地震波传递引起的地面 质点运动。从模拟窄带地震计到现今宽频带、大动态的数字地震台网，观测仪器不能直 接无畸变地记录整个地面运动过程。早在20 世纪初地震仪器理论形成时，科学家已经 意识到此问题。受当时客观条件所限，该问题研究进展较为缓慢。到了现代，伴随着地 震学的快速发展，开展了震源机制、地震波传递动力学特征与运动学特征及提取前兆信 息等研究工作，如果不考虑去除仪器自身特性对地面运动的畸变，得出的最大振幅、拐 角频率、纵横波振幅比等参数可能就是错误的。因此，如何通过地震记录图去实现真实 地面质点运动的复原很有必要。地震学家Aki 及Richard 归纳4 种去除仪器响应的方法， 其中两种方法用于时间域，另外两种需要变换到频率域，然后反变换复原到时间域。时 间域的处理方法中，金星等（2004 ）对强震动记录数据作了较好地处理，并提出递推滤波 公式，可以实时仿真地震速度，并用于地震记录。在频率域方法方面，频率域处理效率高， 缺点是无法实时处理，李鸿吉（1992 ）利用“引用周期”的方法，避免反卷积过程产生的高 频与低频失真，并利用现代控制理论中的解状态方程，验证反卷积结果，取得较好的效果； 刘瑞丰等（1997 ）及何少林等（1997 ）分别沿用傅里叶变换与构造记录系统效果等价稳定滤 波器的办法，对长周期与中长周期的地震记录进行仿真，也取得成功。 1 频率域中恢复记录方法及缺陷分析 在频率域中，对地震记录进行地面运动恢复，基本原理是，基于线性时不变系统的 时域特性，输出信号等于系统仪器响应与输入信号的卷积y (n) h(n) x (n) ，对其进行 FFT 运算，转换到频率域，其频率域表达式为Y(w) H (w) X (w) ，则X (w) Y(w) / H (w) ，然 后对X (w) 进行IFFT ，转换回时域，得到x (n) IFFT (Y (w) / H (w)) 。此为频率域中恢复记 作者简介：高亮（1982—），男，湖北监利人，研究生，工程师，2009 年毕业于中国地质大学研究生院，现从事 地震监测与分析工作 本文收到日期：20 14-12-26 第2 期 高 亮等：利用重叠加窗方式恢复地面运动 37 录的一般方法，本质上是一个反卷积过程。在实际地震记录数据转换处理中，直接利用 该方法处理数据，往往会出现一些问题。这是因为，目前数字地震观测系统频率域响应 H (w) ，其振幅谱在高频段与低频段均趋近于零，则1/ H (w)趋近于无穷大。而对于需要进 行地面运动复原处理的地震记录来说，数据是有限长的，实际上等同于对无限长数据截断， 天然加上一个矩形窗，在对地震记录进行FFT 运算时，在高低频段，特别是低频端会 产生谱泄露，频谱分析中，矩形窗的旁瓣峰值幅度是 13 dB ，与主瓣相比，最大处仅