地震 matlab,MATLAB在地震信号处理中的应用实例.docx

MATLAB在地震信号处理中的应用实例

谭雨文I刘国明2

(I.丰满地震台，吉林 丰满132108 : 2长白山天池火山监测站，吉林 安图133613)

摘要 实现数抿的可视化是地震监测和硏究工作中最重要的问题之一。本文介绍了 MATLAB在信号处理 工作中的主要特点，结合实例以代码的形式介绍了几个台站工作中最常遇到的数据可视化编程方法，通过 分析指出MATLAB实现地卷数据可视化和其它实现数据可视化方法相比具有方便灵活，图形质虽高的特 点,

关键词：MATLAB:可视化；数据处理；地丧

中图分类号P315. 6文献标志码：A 文章编号:1674-8565 (2011) 03-0061-06

0引言

U前我国已经基本完成了由模拟地震信号 监测到数字化地震观测的转变。数字化地震监 测的大力开展给地震研究人员提供了更广阔的 研究平台，也加快了地震研究工作向更深层次 的发展，随着计算机技术的普及和数字信号处 理技术的发展，对地震台站产出的地震和前兆 数据进行有效的加工处理已成为可能"

MATLAB是20世纪70年代发展起来的一 门新的计算机高级语言。MATLAB语言更直观' 更简洁的代码代替了传统的C和FORTRAN等 其它语言的冗长代码，给用户提供了最简洁最 直观的开发环境：该语言主要有以下特点：d 语言简洁紧凑，使用方便灵活，有相当丰富有 效的库函数；②既具有结构化的控制语句，也 可以实现面向对象编程；③具有强大的图形处 理功能；④有功能强大的工具箱，工具箱包含 了最新的数据处理方法，可以方便快捷地实现 数据处理的*的；⑤应用MATLAB绘出的图形 具有可编辑性且图形质量高，能满足出版要求：

基于以上特点.MATLAB正在越来越广泛 地应用于地震观测研究实践中。

1 MATLAB用于地震频谱分析

I. 1地震频谱分析方法简介

所谓地震频谱分析是通过快速傅丽叶变换 (FFT)获得一个地震事件的优势频率.从而实 现信号从时间域到频率域的转化。频谱分析的 方法可以获得地震事件的优势频率，从而判断 地震的类型，在火山地震类型识别和爆破分析 中有非常广泛的应用，

MATLAB信号处理工具箱函数FFT用于快 速傅丽叶变换，可以快速实现地震频谱分析, 具体实现的时候需要将地震记录原始数据通过 EDAS-IAS软件或者NEWSEIS等地震分析软何 将原始地震数据转换成ASCII码，并将ASCII 码文件中的数据头文件去掉，使数据文件变成

收稿日期：2011-06-10 修订日期：2011-07-11

作者简介：谭雨文(1966-).男.汉族吉林省舒兰市人.1990年毕业于长春地质学院地球物理系.T-程师.现主要从事地後分 析与预报研究工作，

可被MATLAB识别的形式。图1是一个 MATLAB实现地震频谱分析的实例，所用的数 据是长口山天池火山监测站记录到的一个火山 地震.通过FFT变换得到了该火山地震的优势 频率为4H乙

1. 2 MATLAB实现地震频谱分析实例代码 地震频谱分析MATLAB代码：

fid=fopen (*200412172303. dat) ; %读取三 分向地震波形数据

a = fscanf

(fid ‘％门;

n = length

(a) 13;

for i = 1 ： n

i1= (i-

11*3+1： i2 = il + 1; i3打2+1;

uz (i)

=a (il) ;(?小

ue (i) =a (i2) ; un

(i3) =a (i3) ; %提取垂直向波形数据uz end

subplot (2, 1. 1):

plot (uz) ; %绘垂直向波形图 xlabel ('Ils1);

ylabel (* Magnitude I count*);

subplot (2, 1. 2):

dt=O. 02; %设置采样间隔为0.Q2

x=0： dt： (n-1) * dt；

n2 = n!2;

s = fft (uz) ; %实现快速傅氏变换

df = l. 01 (n * dt) ； %归一化处理 f

=0： df： (n2-1) *df；

plot (f, abs (s (1: n2) ) I max (abs (s

(1: n2))));绘归一化振幅谱

xlabel CflHz*);

ylabel ('Spectra');

Vs

0.8

0.6

O.4

0.2

O5

1O15

f7Hz

图1地震信号的原始波形和频率谱

20

25

Fig. 1 Original signal and its spectra

2地震月频次图实例

2. I地震月频次图简介

地震月频次图是指以月份为横坐标，以每月 发生的地震数量为纵坐标，画出的统计图形。地 震)］频次图简明的反应地震数量在时间尺度上 的变化，在地震分析