地震波分析软件SAC2000.pptxVIP

地震波分析软件SAC2000主讲：张朝玉 2007年3月6日内容目录 SAC软件简介； SAC基本操作命令； 实验任务； 附录；SAC软件简介 SAC(Seismic Analysis Code)是美国加州大学Lawrence Livermore国家实验室(LLNL)开发研制的，用于研究连续信号，特别是时间序列数据的通用程序。 SAC是基于UNIX平台下的波形分析软件，主要针对SAC (一种IRIS和全球数字地震台网通用的二进制压缩记录格式)记录格式波形的读入和分析。IRIS：Incorporated Research Institutions of Seismology Sac数据文件格式 数据文件由两部分组成：一个固定长度的 头段和位于其后的一个或两个数据段；头段：对字母数字码的sac文件，头段储存在前30行，后面紧跟着数据段。头段变量介绍头段变量介绍头段变量介绍头段变量的数据类型 SAC的基本功能； 分析功能：算术运算，傅氏变换，频谱估计技术，无限脉冲响应(IIR)和有限脉冲响应(FIR)滤波，信号叠加，数据抽取，数值内插，相关分析，地震震相读取等功能； 图形功能：简述图形设备和每个图形命令。说明 SAC安装在SUN服务器上，IP为：5，机器名为：SUN2500。 为大家开设了学习帐号，user: stu11401，password: stu2005。 大家可以用telnet登陆，但是不支持图形功能，建议使用如下软件登陆：x-win32，EXCEED， F-Secure SSH登陆，方便绘图。SAC基本操作命令1. SAC的 进入：sac 推出：quit2. 产生sac数据文件； 使用命令：FUNCGEN 语法：FG {type},{DELTA v},{NPTS n},{BEGIN v} Type: IMPULSE, STEP, BOXCAR, TRIANGLE, SINE {v1 v2}, LINE {v1 v2}, QUADRATIC {v1 v2 v3}, CUBIC {v1 v2 v3 v4}, RANDOM {v1 v2}, SEISMOGRAM; user: fg sine 2 delta 0.1 npts 2003. 文件读入命令：READ，READALPHA，READSDD(CSS)，READHDR，READSP； user: r p0400.p0600.sac注释：命令大小写都可以，桔色字母表示这个命令的 简写。内存内最多能读入100 个文件， 选用more 参数（ 如r more file11,file22）可以在保持内存文件下，读入新的文件。但若无more 参量，第二次读入的文件将刷新内存中读入的数据。4. 文件写出命令： 输出sac 或ASC波形文件 WRITE {SAC | ALPHA} {OVER | filelist} WRITEHDR，WRITESDD ， WRITESP {RLIM | AMPH | RL | IM | AM | PH} { OVER | filelist} user: w mine.sac5. 查看头文件：LISTHDR； 改变头文件内容：CHNHDR user: lh p0400.p0600.sac user: ch LEVEN TRUE6. 文件打开和关闭命令； OAPF: 打开字母数字型的震相读取文件 CAPF: 关闭字母数字型的震相读取文件 OHPF: 打开HYPO格式的震相读取文件 CHPF: 关闭HYPO格式的震相读取文件 WHPF: 写内容到HYPO格式的震相读取文件中 user: oapf tiger关于字母数字震相文件的格式7.图形控制命令 开始绘图到一个或多个图形设备上：BEGIN DEVICES devicesdevices是下面之中的一个或多个：TERMINAL, SGF, XWINDOWS, SUNWINDOWSENDDEVICES devicesERASE 注：目前工作站仅支持xwindows图形窗口8. 绘图命令； PLOT：在一个窗口绘制一条曲线； PLOT1：在一个窗口绘制多条曲线，共用x轴，而每条曲线有自己的y轴； PLOT2：在一个窗口绘制多条曲线（重叠）； PLOTSP：给频谱文件绘图； user: read p1800.p2000.sac p1900.p2100.sac p2000.p2200.sac user: p {p1 p2} user: fft user: psp9. 滤波命令；FFT：对时间序列数据作傅氏变换IFFT：对时间序列数据作傅氏反变换BANDPASS {BUTTER|BESSEL|C1|C2}, {CORNERS v1 v2}, {NPOLE