Fk上机实验报告作业4.docxVIP

Fk上机实验报告PB007155廖一帆一．方法简介该程序用于计算层状介质下的各种源函数（如爆炸，单力和双力偶）产生的合成地震波形。二．实验内容本次实验使用fk方法正演地震波形，实验中所用的命令如下：./fk.pl -Mhk/15/k -N1024/0.1 05 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80第一个Mhk/15/k，表示以hk为模型并且文件名是hk，深度为15公里，文件是很快类型。-N1024表示数据数目有1024个。实验中所用的hk文件内容截图如下：这些数据分别表示厚度 横波波速 纵横波速比 [rho Qs Qp]（密度，S波Q值，P波Q值三选一）-N1024/0.1：表示以0.1s为采样间隔，共采样1024个点，即最终数据的总时长为102.4s其余数字为目标震中距，单位为公里。下面开始将hk文件变位地震波形图syn -M4.5/0/90/0 -D1 -A0 -OPAS.z -Ghk_15/30.grn.0syn程序将格林函数按所需的震源参数整合为地震波形，-M4.5/0/90/0中-M4.5表示地震的震级大小为4.5级，0/90/0为震源参数，即走向0度，倾角90度，滑移0度的断层，-D表示震源函数持续时间，-A0表示方位角为0度，-OPAS.z中PAS为输出的文件名最终的输出文件类型是.z，-G后面加的是格林函数第一分量的名字。将两个不同层状地层模型分别得到的.z文件用pssac绘制到同一张图上实验中所用的命令为：psbasemap -R0/50/0/100 -B10:time(s):/10:distance(km):WSne -JX6i -K -P -Y2i sac.pspssac hk_15/*.z -R0/50/0/100 -Ekt-3 -JX6i -K -O -M1sac.pspssac hk2_15/*.z -R0/50/0/100 -Ekt-3 -JX6i -K -O -M1 -W1/255/0/0sac.pshk2文件是将hk文件中的第二行中的10.5公里地层分为一个和原来一样的5.5公里地层，以及一个5公里厚，横波波速为3.4的地层。即第二个层状地层模型含有一个低速层。第一句先画了底图，第二句语句用黑色线条描绘第一个模型的数据，第二个语句用红色线条描绘第二个模型。最后可以得到下面的图像： 图一、两个不同层状地层模型的波形图图像按照发震时刻对齐，在震中距较小的时刻，到时几乎完全一样没有区别，但是当震中距逐渐增大，含低速层即模型二的到时逐渐增大，模型一与模型二的到时之差越来越大。将同一个模型（hk模型）下得到的三分量地震数据绘制到同一个图中实验过程中的命令如下：psbasemap -R0/50/0/100 -B100/10 -JX2i/6i -K -P -Y2i dist.ps:这个语句是绘制底图绘制底图 下面三个语句是绘制三分量图pssac -R0/50/0/100 -Ekt-3 -M.5 -K -O -JX2i/6i hk_15/.z dist.ps将z分量数据绘制到底图上psbasemap -R0/50/0/100 -B100/10 -JX2i/6i -K -O -X2.1i dist.ps pssac -R0/50/0/100 -Ekt-3 -M.5 -K -O -JX2i/6i hk_15/.t dist.ps将t分量数据加到底图上psbasemap -R0/50/0/100 -B100/10 -JX2i/6i -K -O -X2.1i dist.ps pssac -R0/50/0/100 -Ekt-3 -M.5 -K -O -JX2i/6i hk_15/.r dist.ps将r分量数据加到底图上最后得到的图像如下：图二、同一hk模型下得到的三分量地震数据波形图图中左边第一张图为z分量数据，第二张图为t分量数据，第三张图为r分量数据。将同一模型（hk模型）、同一震中距（30km）、不同方位角的波形图按发震时刻对齐，绘制波形图实验中所用的命令如下：gawk BEGIN{for(i=0;i360;i+=5) print sh az.sh,i} |sh从方位角为0度开始到360度结束，每5度提取波形数据psbasemap -R0/20/0/360 -B100/10 -JX6i -K -P -Y2i sac1.ps绘制底图pssac -R0/20/0/360 -Eat-3 -M.1 -K -O -JX6i a*.zsac1.ps将提取的不同方位角的波形数据绘制到底图上，其中-E选项后面的第一个参数决定sac文件的Y轴排列方式，此处为a，表示按方位角排列。最后得到的图像为：图三、同一hk模型、同一震中距、不同方位角的波形图从图像上可以看出波形存在明