地震matlab频域分析,MATLAB在地震数据分析中的应用

MATLAB在地震数据分析中的应用 中国矿业大学信息与电气工程学院 薛慧 顾玲玲中国矿业大学计算机学院 谢晓宇 摘要：奉文主要介绍运用M^TLAB，将宽带滤波嚣的精八仿真到窄带滤波器上，得到和窄带滤波嚣一样的波形，较为精确地恢复了地面运动。 关键词：MATLAB仿真位移泣形恢复 地震数据分析是地震预报的基础环节，可以提高地震分析的准确度，在地震观测过程中起到了关键的作用。MATLAB数字滤波器为地震分析提供了优秀的工具．对于去除地震观测波形中的各种干扰和地面位移波形的恢复表现得特别优秀。 一、基本原理 在地震记录处理中位移波形的恢复一直作为一个基本问题加以研究：地震仪的传递函数有多重零点，加上地震仪的非线性仪器噪声地震仪和地震记录系统有限的分辨力，真实地面位移只能被 近似地恢复。 在线性系统中．一个系统可以用它的系统函数或脉冲响应来表示‰圳@x(1J 式中．对于地震观测系统来讲．x(t)为输入信号即地面运动．Y (t)为系统的输出，即地震记录数据．h(t)为系统的脉冲响应。在频率域内，该式可以表示为 式中，H(u)、Y(∞)、x(u)分别为h(t)、Y(t)、x(0的傅里叶变换： 设一个频带很宽的线性系统．系统甬数为H(u)；另一个频带很窄的系统．系统甬数为H’fu】，频谱为Y’(∞)。对于同样的输入Xful有 y(ca)=H(t^，)x(∞)．Y’(∞)=H’(uJx(u)得到．HI掣 Y’(u)=H(u) 变换到时域得到频带较窄系统的输出y’(t)。由此可得，若知道了宽频带和窄频带的传递函数，理论上可以从宽频带系统的输出推测出频带窄系统的输出。这样就可以在记录某种信号时采用宽频带记录，然后仿真到其他各种窄频带的记录仪器上对信号进行分析。 二、仿真处理 1仿真要求 在做仿真处理时，要求原始记录系统频带宽，具有稳定明确的传递函数，被仿真的地震系统传递函数也应明确稳定。 2传递雨数的变换 对于具体地震记录系统，要实现计算机仿真首先要将模拟域的传递函数变换为离散域中的传递函数。在z域中地震记录系统传递函数的幅频响应首先计算从连续域到离散域的变换．然后计算在离散域中的幅频响应。从连续域到离散域可采用双线性变换法。条件是离散域频率响应和模拟域频率响应一致。其公式为： 7’ 卜z—l I+z—t 双线性变换引入了一个非线性频率畸变，即双线性变换把在无限区间内的模拟滤波器的频率响应映射到数字滤波器中，是从直流到fo的有限段，且fO=fs／2(fs为采样频率)：实际常用的几种地震仪是频率范围有限的滤渡器．IjCDSN的采样率是仪器最大频率的两倍：因此．在记录仪的主频段内可得到最佳的双线性变换效果．不 同的只是在高频部分产生频率畸变； 3 MATLAB仿真结果分析 下面举例给出将宽频带系统的输出仿真到窄频带输出及地面 运动恢复的方法。 设计一个巴特沃斯模拟宽频带带通滤波器．设计指标为：通带频率：1-40Hz．低端阻带边界0．2Hz，高端阻带边界42 5Hz．通带最大衰减ldB，阻带最小衰减20dB。再设计一个窄频带带通滤波器．设计指标为通带频率：5-7Hz，低频段过渡带宽1 5Hz．高频段过渡带宽0 5Hz．通带最大衰减ldB．阻带最小衰减20dB。假设一个 信号。x(O=sinI二Ⅳ·0．5cos20n／,-0 5sin38nt fs=50Hz。 仿真结果 运用宽带仪器的输出、宽带和窄带仪器的传递甬数得到的窄 带仪器的输出(仿真)．并与原窄带仪器的输出进行比较 在上面的仿真过程中．设计的频带较宽．几乎包含了所有输入信号的频率。输入信号经过宽带仪器输出后波形基本