地球物理上机实验报告范文3.docxVIP

地球物理上机实验报告一.数据描述 本次上机实验选取的地震数据为亚洲部分台站所记录的甘肃定西地震的数据，一共有来自19个不同台站的36道数据，台站分布如图1所示。震中位于东经104.26°，北纬34.51°，震源深度8.0km，发震时刻为格林尼治时间2013年7月21日23时45分。图1：包含震源机制解的台站分布图。其中黑色实心三角表示台站的位置，震源球表示震中位置及震源机制解。二.实验内容1.绘图前对数据的预处理： （1）使用TauP软件将prem模型下的P波理论到时标注的SAC文件中，写入头文件t0中； （2）使用SAC命令ppk，人工标注P波到时，写入头文件t1； （3）使用SAC标注P波初动方向，命令为P，写入头文件t2，“+1”表示初动向上，“-1”表示初动方向向下。2.绘图 （1）震中距—走时曲线图（图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)} PB11007.txtpsbasemap -R0/4000/0/600 -B10:distance(km):/200:time(s):WSne -JX6i -K -P -Y2i sac.pspssac *BHZ* -R0/4000/0/600 -Edt-3 -JX6i -K -O -M.1 -Vsac.pspsxy PB11007.txt -R0/4000/0/600 -JX6i -K -O -Sx.5 sac.pscat PB11007.txt |sort -n -k1 |psxy -R -JX -K -O sac.ps图2：震中距—走时曲线图。其中黑色线为手工标注的P波到时曲线，红色线为prem模型下P波的理论到时曲线，绿色线为ak135模型下的P波理论到时曲线。（2）包含P波初动的震源机制解对应的震源球（图3）图3：包含P波初动的震源机制解对应的震源球。其中“X”表示初动方向向下，“O”表示初动方向向上。三.结果讨论1.从走时曲线（图2）可以看出，理论到时曲线和手动标注到时曲线大致相符，但是在震中距较小和较大的地区，理论到时比手动标注的到时偏大，可能是模型中的速度设定偏低，或者是手动标注存在的误差的原因。2. 包含P波初动的震源机制解对应的震源球（图3）大致结果还行，少数的几个点不很理想，可能是由于手动标注P波初动存在的误差，或者在出射角估计以及投影半径计算的误差而导致的。