SAC上机实验报告作业范文2.docxVIP

SAC上机实验报告P乔垒数据描述这次实验是学习使用sac和gmt等软件，并且学习脚本的简单编写与运行。选取的地震事件是甘肃定西的地震。一共选取了亚洲19台站的36道数据。台站基本分布在地震事件的四周。实验所用数据居委宽频带的Z分量数据。震中位置为（104.26E，34.51N），震源深度为8km。发震时刻为格林尼治时间2013年7月21日23时45分55秒999毫秒。下表为台站信息：台网台站经度/度纬度/度仪器编号HKHKPS114.1422.28ICBJT116.1740.0200,10ICENH109.4930.2800,10ICHIA119.7449.1700,10ICKMI102,7425.1200,10ICLSA91.1329.7000,10ICMDJ129.5944.6200,10ICQIZ109.8419.0300,10ICSSE121.1931.0900,10ICWMQ87.7043.8100,10ICXAN108.9234.0300,10IIAAK74.4942.6400,10IIKURK78.6250.7200,10IINIL73.2733.6500,10IUCHTO98.9418.8100,10IUKBL69.4334.54IUMAJO138.2036.5500,10IUMAKZ81.9846.8100,10IUTATO121.5024.9700,10实验中所用台站的具体位置信息内容的图一中。实验内容第一部分：计算不同台站的理论到时，绘制理论到时与时间到时曲线Step1：使用Taup命令中的taup_setsac，计算prem模型中的p波理论到时并写入头文件t0中Step2：使用ppk命令手动标注p波到时t1脚本为:psbasemap -R0/4000/0/600 –B500:distance(km):/200:time(s):WSne -JX6i/9i -K -P -Y2i deg_time.ps#绘制地图taup_curve -mod prem -h 8 -ph P#计算prem模型下的p波理论到时，并保存到tauo_curve.gmt文件中，第一列为震中距（角度值），第二列为理论到时（秒）cat taup_curve.gmt|gawk NR1{print $1,$2}|sort -g |psxy -R -JX -P -K -O -W3/255/0/0 deg_time.ps#将理论到时曲线画在底图上，其中要用sort命令将数据按震中距从大到小排序，使所画曲线有意义pssac *BHZ* -R -Edt-3 -JX -K -O -M.1 -V deg_time.ps#将地震记录绘制在底图上，-Edt -3表示按sac头文件中o变量的值校正波形的时间坐标，时间的发生时刻相对于参考时刻的秒数。参考时刻设为发震时刻，因此校正后即表示相对于地震时刻的波形图。即发震时刻为0.M表征振幅。ls *BHZ* |gawk {print r $1;print lh gcarc t1} gcarc_t1#读取震中距及手动标注P波初动时刻echo q 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 deg_time.txt#将实际到时写入gcarc-t1文件psxy deg_time.txt -R -JX -K -O -W3 deg_time.ps#绘制实际到时曲线成图为：图一 理论与实际到时曲线图由图可知，理论到时曲线与实际到时曲线在中间的位置符合的比较好，但是在近震或者远震位置拟合的效果不好。第二部分：将台站位置以及地震位置和震源机制解绘制在底图上直接绘制将台站以及震源球绘制在地图上脚本为：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.txtcat 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 -