fk上机实验报告作业4

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

1、Fk上机实验报告廖一帆PB007155一.方法简介该程序用于计算层状介质下的各种源函数(如爆炸，单力和双力偶)产牛的合成地震波形。二.实验内容本次实验使用fk方法正演地震波形，实验中所用的命令如下：./fk.pl-MhV15/k-N1024/0.105101520253035404550556065707580第一个Mhk/15/k,表示以hk为模型并且文件名是hk,深度为15公里，文件是很快类型。-N1024表示数据数目有1024个。实验屮所用的hk文件内容截图如下:5.53.181.73060010・53.6

2、41・73160016・03.871・73160090.04.501.733900这些数据分别表示厚度横波波速纵横波速比WhoQsQp](密度，S波Q值，P波Q值三选一)-N1024/0表示以0.1s为采样间隔，共采样1024个点，即最终数据的总时长为102.4s其余数字为目标震中距，单位为公里。下面开始将hk文件变位地震波形图syn-M4.5/0/90/0-DI-AO-OPAS.z-Ghk_15/3O.grn.Osyn程序将格林函数按所需的震源参数整合为地震波形，・M4.5/0/90/0中•M4.5表示地震的震

3、级大小为4.5级，0/90/0为震源参数，即走向0度，倾角90度，滑移0度的断层，・D表示震源函数持续吋间，・A0表示方位角为0度，-OPAS.z中PAS为输出的文件名最终的输出文件类型是.乙，・G后面加的是格林函数第一分量的名字。将两个不同层状地层模型分别得到的n文件用pssac绘制到同一张图上实验中所用的命令为：psbasemap・R0/50/0/100-B10:ntime(s),,:/10:,,distance(km)n:WSne-JX6i・K-P-Y2i>sac.pspssachk_15/*.z-R0/5

4、0/0/100-Ekt-3-JX6i・K-O・Ml»sac.pspssachk2_15/*.z-R0/50/0/100-Ekt-3-JX6i-K-O-Ml-Wl/255/0/0»sac.pshk2文件是将hk文件中的第二行中的10.5公里地层分为一个和原来一样的5.5公里地层，以及一个5公里厚，横波波速为3.4的地层。即第二个层状地层模型含有一个低速层。第一句先画了底图，第二句语句用黑色线条描绘第一个模型的数据，第二个语句用红色线条描绘第二个模型。最后可以得到下面的图像：(EM)①0U5S一p1000102030

5、4050time(s)图一、两个不同层状地层模型的波形图图像按照发震时刻对齐，在震中距较小的时刻，到时儿乎完全一样没有区别,但是当震中距逐渐增大，含低速层即模型二的到时逐渐增大，模型一与模型二的到时z差越来越大。将同一个模型（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

6、-JX2i/6ihk_15/.z»dist.ps将z分量数据绘制到底图上psbasemap-R0/50/0/100-B100/10-JX2i⑹-K-O・X2.1i»dist.pspssac-R0/50/0/100-Ekt-3・M.5・K-O-JX2i/6ihk_15/.t»dist.ps将t分量数据加到底图上psbasemap-R0/50/0/100-B100/10-JX2i/6i-K-O・X2.1i»dist.pspssac-R0/50/0/100-Ekt-3-M.5-K-O-JX2i/6ihk_15/.r»d

7、ist.ps将「分量数据加到底图上最后得到的图像如下：9080706050403020100图二、同一hk模型下得到的三分量地震数据波形图图中左边第一张图为z分量数据，第二张图为t分量数据，第三张图为I•分量数据。将同一模型（hk模型）、同一震中距（30km）、不同方位角的波形图按发震时刻对齐，绘制波形图实验中所用的命令如下：gawk'BEGIN{for（i=0;i<360;i+=5）printHshaz.sh“,i}‘

8、sh从方位角为0度开始到360度结束，每5度提取波形数据psbasemap-R0/20/0/

9、360-B100/10-JX6i-K-P-Y2i>sacl.ps绘制底图pssac-R0/20/0/360-Eat-3-M.1-K-O-JX6ia*.z»sacl.ps将提取的不同方位角的波形数据绘制到底图上，其中・E选项后面的第一个参数决定sac文件的Y轴排列方式，此处为a,表示按方位角排列。最后得到的图像为：翱芻需翱据駅芻需翱黑翱芻話翱据翱黑需O33^33332222