地球物理上机实验报告范文4.docVIP

对地震数据的处理( SAC、GMT使用） P 张盛毅 实验要求 1、绘制2、绘制3、绘制实验数据说明 此次研究的地震是发生于北京时间2013年07月22日07时45分（）在甘肃省定西市岷县、漳县交界（北纬34.51度，东经104.26度）的6.6级地震，其震源深度8km。 （网上查到的却是20km？） 数据来源于由记录得到的共，基本分布在地震事件的四周。实验所用数据为BHZ分量数据。 度9.43 34.54 IU MAJO 138.20 36.55 00,10 IU MAKZ 81.98 46.81 00,10 IU TATO 121.50 24.97 00,10 实验内容 计算到时，并绘制理论到时曲线 使用taup_setsac用于将震相到时信息写入SAC头端t0中，其命令格式如下： taup_setsac -mod prem -ph P-0,S-1,PcP-8,ScS-9 -evdpkm *.U 命令说明 -mod?prem?：指定地球模型为PREM模型，也可以指定其他模型，如ak135； -ph?P-0,S-1,PcP-8,ScS-9?：将P波到时标记到T0中，S到时标记到T1中，PcP标记到T8，ScS标记到T9。与此同时，该命令还会修改相应的KTn，可根据需要选取不同参数。 -evdpkm?：该选项设置深度的单位为km（由于历史的原因，该命令默认SAC文件中的震源深度单位为米，而现在的SAC文件中的震源深度单位为千米。） *.U?:?需要标记理论到时的SAC文件，此次用的为*BHZ*文件。 上机实验：taup_setsac -mod prem -ph P -evdpkm *BHZ* 在sac下，使用ppk命令手动标注p波到时t1。 所用的sh文件中的命令如下 ①psbasemap -R0/30/0/600 -B10:distance(deg):/200:time(s):WSne -JX6i -K -P -Y2i sac.ps 用于绘制底图，gmt命令。 ②taup_curve -mod prem -h 8 -ph P 利用走时计算软件TauP计算prem模型下的p波理论到时，并保存到tauo_curve.gmt文件中，第一列为震中距（角度值），第二列为理论到时（秒） ③cat taup_curve.gmt |gawk NR1{print $1,$2}|sort -g |psxy -R0/30/0/600 -JX6i -K -O -W2p/255/0/0 sac.ps 将理论到时曲线画在底图上 ④pssac *BHZ* -R -Edt-3 -JX -K -O -M.1 -V sac.ps 将地震记录绘制在底图上，通过-Edt -3进行到时校正，得到相对于发震时刻的波形图。-M.1表示振幅的0.1倍。 ⑤ls *BHZ* |gawk {print r $1;print lh gcarc t1} gcarc_t1 读取震中距及手动标注P波初动时刻。 ⑥echo q gcarc_t1 在gcarc_t1末尾写入q，表示退出。 ⑦cat gcarc_t1 |sac |gawk NF==3{if($1==gcarc)printf(%f\t,$3);if($1==t1)printf(%f\n,$3)}|sort -g sac.txt 实际到时写入gcarc_1文件中。 ⑧psxy sac.txt -R -JX -K -O -W3 sac.ps 绘制实际到时曲线 理论到时曲线与实际到时曲线在distance为10°到25°之间符合的比较好，但其它位置拟合的效果不太好。红线为prim模型。 绘制台站位置及震源球 gmt pscoast -R60/140/0/65 -JM6i -K -P -B10 -Y2i -Dh -A200 loc.ps 绘制台站以及地震所在位置的底图 ls *BHZ* |gawk {print r $1;print lh stlo stla} stlo_stla.txt 提取台站位置信息 echo q stlo_stla.txt cat stlo_stla.txt |sac |gawk NF==3{if($1==stlo)printf(%f\t,$3);if($1==stla)printf(%f\n,$3)} |gmt psxy -R60/140/0/65 -JM6i -K -O -St.5 -G0/0/0 loc.ps 绘制台站，用黑三角表示 gmt psmeca cmt_1 -R60/140/0/65 -JM6i -K -O –Sc1 -Gred loc.ps 在地震位置绘制震源球，cmt_1数