sac上机实验报告作业范文7

《sac上机实验报告作业范文7》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、6(r70°80°90°100°110°120°130°140°60650°40°30°2(r10°0°DSAC实验报告PB11007115许昊天一.数据描述本次上机实验选取的地震数据为亚洲部分台站所记录的甘肃定西地震的数据，一共有来自19个不同台站的36道数据，台站分布如图1所示。震中位于东经104.26°,北纬34.51°,震源深度8.0km,发震时刻为格林尼治时间2013年7月21日23时45分。60°50°40°30°20°10°0°60°70°80°90°100°110°120°130°140°图1：包含震源机制

2、解的台站分布图。其中黑色实心三角表示台站的位置，震源球表示震中位置及震源机制解。二.实验内容1•绘图前的预处理:1）使用TauP软件将prem模型下的P波理论到时标注的SAC文件中，写入头文件tO中；2）使用SAC命令ppk,人工标注P波到时，写入头文件变量tl；3）使用SAC标注P波初动方向，命令为P,写入头文件变量t2,以“+1”表示初动向上，以表示初动方向向下，“0”表示初动方向无法判断。2•绘图1）震中距一走时曲线图（图2）所用脚本：Is*BHZ*

3、gawk'{print$l;print〃lhgcarcsac200

4、0

5、gawk'NF==3{if（$l==z/gcarcz/）printf$3）;if（$1二二〃tl〃）printf（〃％f〃，$3）}'>PB11007115.txtpsbasemap-R0/4000/0/600-BIO:"distance（km）〃：/200:"time（s）〃：WSne-JX6i-K-P一Y2i>sac・pspssac*BIIZ*-R0/4000/0/600-Edt-3-JX6i-K-0-M.1-V>>sac.pspsxyPB11007.txt-R0/4000/0/600-JX6i-K-0-Sx.

6、5»sac.ps2005001000150020002500300035004000distance（km）图2：震中距一-走时曲线图。其中黑色线为手工标注的P波到时曲线，红色线为“cm模型下P波的理论到时曲线，绿色线为akl35模型下的P波理论到时曲线。2）包含P波初动的震源机制解对应的震源球（图3）□图3：包含P波初动的震源机制解对应的震源球。其中“X”表示初动方向向下，“0”表示初动方向向上。对于不确定的点，可能同时出现“0”和X。三•结果讨论1•从走时曲线可以看出，理论曲线和手动标注曲线整体吻合较好，在震屮距较大的

7、区域存在一定误差。这应是由手动标记的误差和理论模型的误差两方面原因共同导致的。2.震源球的效果主要由手动标记P波初动方向的准确性决定。从整体來看吻合较好，个别波形的比较难以分辨前期背景噪声和P波初动的分界，导致判断不准。四•对课程的建议个人感觉对于GMT以及Linux下常用操作还不是很熟悉，上课的时候往往就只是跟着老师敲代码，对于程序具体细节、参数的理解不够。建议老师可以增加布置课后编程作业，让大家的编程基本功更扎实。