地球物理实验报告Fk 上机实验.docxVIP

Fk 上机实验 ——P马天一、实验目的 应用实验中的程序计算地下层状介质的各种源函数（如爆炸，单力和双力偶）产生的合成地震波形。二、实验内容本次实验使用fk方法正演地震波形，实验中所用的命令如下。 第一步：-Mhk/15/k：表示使用当前目录下的hk文件作为计算格林函数时参考的一维速度模型文件，其中震源深度为15公里，最后的k参数表示hk的文件格式，其格式为： 厚度 横波波速纵横波速比[rho Qs Qp]（密度，S波Q值，P波Q值三选一）-N1024/0.1：表示采样1024个点以0.1s为采样的时间间隔，即最终数据的总共采样时长为102.4s 。其余数字为目标震中距，单位为公里。实验中所用的hk文件内容截图如下第二步： 使用syn程序将格林函数按所需的震源参数整合为地震波形，-M4.5/0/90/0中-M4.5表示地震的震级大小为4.5级，0/90/0为震源参数，即走向0度，倾角90度，滑移角0度的断层（直立走滑断层），-D1表示震源函数持续时间为1s，-A0表示方位角为0度，-OPAS.z中PAS为输出的文件名，-G后面加的是格林函数第一分量的名字。第三步：上诉命令是将两个不同层状地层模型分别得到的.z文件用pssac绘制到同一张图上 hk文件： hk2文件： 从上面两个地层模型可以看出来，第一个层状地层模型地震波的速度随着深度的增加而增大，而第二个层状地层模型含有一个低速层。最后可以得到下面的震中距-时间图像：图一、两个不同层状地层模型的波形图分析：其中黑色线表示第一个模型的数据，红色线表示第二个模型的数据，按发震时刻对齐。由图可见在震中距较小的位置两个模型的到时相差不大，几乎是一致。随着震中距的增加，含低速层模型的到时相对于第一个模型要长一些，而且两者之间的到时差越来越大。第四步：将在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 在原底图上绘制第二个底图，-X2.1i表示第二张相对于第一张右移错开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 在原底图上绘制第三个底图，-X2.1i表示第三张相对于第二张右移错开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 将提取的不同方位角的波形数据绘制到底图上最后得到的图像为：图三、不同方位角的波形图（同模型同震中距）分析：从图像上可以看出波形存在明显的四象限分布特征。特别的在90度，180度，270度可以发现明显的不同之处，可以作为进一步的观察思量。