必威体育精装版《地震学》实验报告.doc

本科生实验报告 实验课程 地震学 学院名称 地球物理学院 专业名称 地球物理 学生姓名 学生学号 指导教师 程先琼 实验地点 5417 实验成绩 二〇一五年 五 月 二〇一五年六月 填写说明 适用于本科生所有的实验报告（印制实验报告册除外）； 专业填写为专业全称，有专业方向的用小括号标明； 格式要求： 用A4纸双面打印（封面双面打印）或在A4大小纸上用蓝黑色水笔书写。 打印排版：正文用宋体小四号，1.5倍行距，页边距采取默认形式（上下2.54cm，左右2.54cm，页眉1.5cm，页脚1.75cm）。字符间距为默认值（缩放100%，间距：标准）；页码用小五号字底端居中。 具体要求： 题目（二号黑体居中）； 摘要（“摘要”二字用小二号黑体居中，隔行书写摘要的文字部分，小4号宋体）； 关键词（隔行顶格书写“关键词”三字，提炼3-5个关键词，用分号隔开，小4号黑体)； 正文部分采用三级标题； 第1章 ××(小二号黑体居中，段前0.5行) 1.1 ×××××小三号黑体×××××（段前、段后0.5行） 1.1.1小四号黑体（段前、段后0.5行） 参考文献（黑体小二号居中，段前0.5行），参考文献用五号宋体，参照《参考文献著录规则（GB/T 7714－2005）》。 学生实验 心得 在第一个实验中我们了解这个网站而该网站是一个获取全球数字地震台网地震波记录及相关波形记录处理软件的重要网站。在第二次实验中我们其实为一次正演模型，即模拟P波入射到核幔边界后的位移反、折射系数与入射角的关系。通过第三次实验我们计算了理论的时距曲线，在近震范围内，可以忽略地球曲率。由于地球介质物性在水平方向的变化远小于其垂直方向变化，因此许多研究中，为简单起见，将地球介质简化成横向均匀各向同性的弹性介质。最后一次实验为计算一维P波的速度，这次实验与上两次实验不同的是，本次实验是反演计算而上两次的实验内容为正演。而这次用到的方法为Herglotz－Wiechert反演方法。 学生（签名）： 年 月 日 指导 教师 评语 成绩评定： 指导教师（签名）： 年 月 日 实验二 计算核幔边界P波反射、 透射系数 如果在核—幔界面（近似当成平面）以上，下地幔的介质密度及P、S波的速度分别为：5.5×103kg/m3,13.7km/s,7.2km/s;外地核的介质密度及P、S波速度分别为： 9.9×103kg/m3,8.0km/s,0.0km/s。用你熟悉的语言，设计一程序计算Ｐ波以任意角度入射时，界面的地震波垂直向位移的反、透射系数，并绘图表示反、透射系数与入射角的关系。 实验原理 程序源代码 clear all; clc; %??D?2?êy x=1:1:90; vp1=13.7; vs1=7.2; vp2=8.0; vs2=0; Density1=5500; Density2=9900; Mu1=vs1^2*Density1; Namd1=vp1^2*Density1-2*Mu1; Namd2=vp2^2*Density2; a1=1:1:90; p=zeros(90); a2=zeros(90); b1=zeros(90); n1=zeros(90); n2=zeros(90); n=zeros(90); E=zeros(90); F=zeros(90); G=zeros(90); H=zeros(90); I=zeros(90); Rp=zeros(90); Rs=zeros(90); Tp=zeros(90); rp=zeros(90); rs=zeros(90); tp=zeros(90); for i=1:1:90; p(i)=sin(a1(i)*pi/180)/vp1; n1(i)=sqrt(1/(vp1)^2-(p(i))