山野莽夫

  • 归档
    • 随笔
    • 建站资源
    • 分享
    • 代码
  • 地球物理学
    • 专业课
    • 概念解释
  • 计算机
  • 互联网
  • 教程
  • 规划
  • 实验室
    • 珍藏的软件
    • 贴吧云签到
    • A1账号自助申请
山野莽夫
小学生的挣扎的点点滴滴
  1. 首页
  2. 地球物理学
  3. 专业课
  4. 地震学
  5. 正文

地震分析软件SAC

2016年5月1日 4462点热度 1人点赞 0条评论

C-1. 针对IzmitMw7.4和 DuzceMw7.1地震,写出完成以下功能的命令操作

# lunch sac                                                                     

# read the vertical-component data recorded by station BV

# plot the data                                                                         

# read all three component data recorded at station BV                                  

# plot multiple data in one page                                                        

# zoom in at 0 - 7 s relative to the waveform reference time (typically the origin time)

# you have to plot it again to reflect the change                                       

# zoom in at 2 s before P, and 4 s after S arrivals                                     

# enforce that y axis has the same scale                                                

# you have to plot it again to reflect the change                                       

# quit sac

C-2. 把SAC数据写回磁盘,写出完成以下功能的命令操作

cp DATA/*.BV.*z* CYC.BP # make a backup of the vertical data                                          

sac # lunch sac                                                                     

read DATA/*.YJ.BV.EHZ.SAC# read the data file YJ.BV.EHZ.SAC                                             

plot # plot it out                                                                 

remean # remove the mean                                                             

w DATA/*BV*Z*.SAC # overwrite the data back to the disk (using the original name YJ.BV.EHZ.SAC) 

w DATA/temp.SAC # Write the data back to the disk using a new name temp.sac                   

quit #quit sac                                                                     

rm DATA/temp.SAC # clean up what just created.                                                 

cd DATA     ls # restore the original data

C-3.列出SAC头文件,写出完成以下功能的命令操作

sac # launch sac again                           

read DATA/*.YJ.BV.EHZ.SAC # read the data YJ.BV.EHZ.SAC                 

LISTHDR # list the SAC header   

LISTHDR evla evlo stla stlo # list the SAC header evla evlo stla stlo

C-4.对地震记录进行旋转,写出完成以下功能的命令操作

DATA/YJ.BV.EHN.SAC # read the north component seismogram                                                  

CHNHDR EVLA 0 EVLO 90# change the incident and azimuth of the instrument to be north and horizontal         

w YJ.BV.EHN.SAC # save the header (the above three lines are not needed if the headers are already in) 

r DATA/YJ.BV.EHE.SAC # read the east component seismogram                                                   

CHNHDR EVLA 90 EVLO 90# change the incident and azimuth of the instrument to be east and horizontal          

w YJ.BV.EHE.SAC # save the header (the above three lines are not needed if the headers are already in) 

w YJ.BV.*Z.SAC # read two horizontal component seismograms

rotate to gcp  # rotate to the "great circle path"

w *.*.BV.*            # save the rotated data (first radial, second transverse)

r *.*.BV.*        # read three-component seismograms (vertical, radial and transverse)

xlim pdw              # zoom in

p1ot                      # plot them out, try to identify Love and Rayleigh waves

C-5.设计不同的滤波器计算FFT,写出完成以下功能的命令操作

r *.PKD.*      # read the transverse-component data recorded at station PKD

 rmean                 # remove the mean

 bp  2:8  # apply two-way 2-8 Hz band-pass-filter to the data

 w  BK.PKD.HHT.SAC.bp   # save the band-pass-filtered data as BK.PKD.HHT.SAC.bp

 r     BK.PKD.HHT.SAC.bp           # read the original data again

 p1ot                           # plot the band-pass-filtered data

 cut on         # cut the data

 r    BK.PKD.HHT.SAC.bp    # read the data again

 rmean                 # remove the mean

 taper type hanning width 0.05       # apply a hanning window to the 0.05 width of the data

 fft wmean amph               # compute FFT

 keepam               # keep only the amplitude information (and throw the phase information)

 p1ot                    # plot the amplitude spectrum

 loglog                 # use the log-log plot

 xlim on        # set the x-axis plot range (Frequency)

 p1ot                           # plot again in log-log scale

C-6.震相拾取,写出完成以下功能的命令操作

cut on# cut the data

r *.BV.*.SAC# read three component data recorded at station BV

p1ot # plot the data

ch a undef      # remove the original P and S arrival pick

apk VALIDATION ON     # use the auto picker

p1ot # see the auto picker result

ppk  om         # use the manual picker

ppk v3           # pick all 3 component together with the same time, turn the annoying bell off

wh 4 #to save the picked arrival, you can use the command "wh". But this command could   mess up the data if you have used the cut command. So be careful  

C-7.编写一个SAC宏文件,完成以下功能

sac # launch sac

read *.BV.*.SAC # read all three component data recorded at station BV                    

mul  4 5 6# multiply data 4,5 ,6 respectivly

fft wmean amph # Fourier transform

quit # quit sac

 

标签: SAC 地震分析软件
最后更新:2016年5月1日

小菜菜

菜鸟

打赏 点赞
< 上一篇
下一篇 >

文章评论

razz evil exclaim smile redface biggrin eek confused idea lol mad twisted rolleyes wink cool arrow neutral cry mrgreen drooling persevering
取消回复

此站点使用Akismet来减少垃圾评论。了解我们如何处理您的评论数据。

联系方式

QQ群 | TG群 | 邮箱

最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
matlab绘制双对数曲线 视速度的解释 Visual Studio项目上传到码云 C语言简单文件读取代码(使用fscanf) 致青春毕业季幻灯片模板的补充 地震子波显示及合成地震记录
标签聚合
ppt 地震学程序 宝塔面板 wordpress c语言 虚拟机 onedrive 模板
最近评论
小菜菜 发布于 7 个月前(11月24日) 这玩意已经废了,成收割工具了,不能再用了。
eamon 发布于 8 个月前(11月07日) 我一年不用了才发现这个休眠管理费每月15,一共扣了我135元,然后我消费还消费不了,我宁愿消费掉也不...
magic 发布于 12 个月前(07月03日) 请问账号不注销会有什么影响吗?
magic 发布于 12 个月前(07月01日) 我想问一下 如果不注销账号就留着会怎么样
qwp6601 发布于 1 年前(06月04日) 有没有方法改为bing

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2