SAC上机实验报告作业范文7.docVIP

SAC实验报告 P昊天 一.数据描述 本次上机实验选取的地震数据为亚洲部分台站所记录的甘肃定西地震的数据，一共有来自19个不同台站的36道数据，台站分布如图1所示。震中位于东经104.26°，北纬34.51°，震源深度8.0km，发震时刻为格林尼治时间2013年7月21日23时45分。 图1：包含震源机制解的台站分布图。其中黑色实心三角表示台站的位置，震源球表示震中位置及震源机制解。 二.实验内容 1.绘图前的预处理： 使用TauP软件将prem模型下的P波理论到时标注的SAC文件中，写入头文件t0中； 使用SAC命令ppk，人工标注P波到时，写入头文件t1； 使用SAC标注P波初动方向，命令为P，写入头文件t2，“+1”表示初动向上，“-1”表示初动方向向下”表示初动方向无法判断。 2.绘图 震中距—走时曲线图（图2） 所用脚本： ls *BHZ* |gawk {print r,$1;print lh gcarc t1;} |sac2000 |gawk NF==3{if($1==gcarc)printf(%f\t,$3);if($1==t1)printf(%f\n,$3)} Ptxt psbasemap -R0/4000/0/600 -B10:distance(km):/200:time(s):WSne -JX6i -K -P -Y2i sac.ps pssac *BHZ* -R0/4000/0/600 -Edt-3 -JX6i -K -O -M.1 -Vsac.ps psxy PB11007.txt -R0/4000/0/600 -JX6i -K -O -Sx.5 sac.ps cat PB11007.txt |sort -n -k1 |psxy -R -JX -K -O sac.ps 图2：震中距---走时曲线图。其中黑色线为手工标注的P波到时曲线，红色线为prem模型下P波的理论到时曲线，绿色线为ak135模型下的P波理论到时曲线。 包含P波初动的震源机制解对应的震源球（图3） 图3：包含P波初动的震源机制解对应的震源球。其中“X”表示初动方向向下，“O”表示初动方向向上。对于不确定的点，可能和三.结果讨论 1.从走时曲线可以看出，在2. 震源球的效果主要由手动标记P波初动方向的准确性决定。从整体来看吻合较好，个别导致四.对课程的建议 对于参数建议