山野莽夫

  • 归档
    • 随笔
    • 建站资源
    • 分享
    • 代码
  • 地球物理学
    • 专业课
    • 概念解释
  • 计算机
  • 互联网
  • 教程
  • 规划
  • 实验室
    • 珍藏的软件
    • 贴吧云签到
    • A1账号自助申请
山野莽夫
小学生的挣扎的点点滴滴
  1. 首页
  2. 地球物理学
  3. 专业课
  4. 地震学
  5. 正文

地震学和地球结构(海洋地壳模型反射波时距曲线)

2016年1月28日 3243点热度 1人点赞 0条评论

 C-1. 编写程序:利用走时和震中距关系(等式 3.3.7, 3.3.8)、及其双曲近似(等式3.3.11 )生成并画出一系列平界面的反射波时距曲线。计算图 3.2-15所示海洋地壳模型的走时,并比较两种方法的计算结果。

  1. 实验原理

    分别运用多个水平层走时曲线精确公式和运用双曲线近似计算海洋地壳模型的每个界面的反射波走时曲线。

    精确公式为

    运用双曲线近似计算公式为


     

二、程序源代码

1.P波

h=[5 2 5 3];i1=zeros(4);

vp=[1.5 5 6.7 8.1];

p=zeros(1,4);

x=zeros(3,41);t=zeros(3,41);

for i=0:40

i1(1)=i;


for j=1:3

p(j)=sind(i1(j))/vp(j);

i1(j+1)=asind(p(j)*vp(j+1));


end


for kk=1:3


for k=1:kk

x(kk,i+1) =2*h(k)*tand(i1(k))+x(kk,i+1);

t(kk,i+1) =2*h(k)/vp(k)/cosd(i1(k))+t(kk,i+1);


end


end

end

x

plot(x(1,:),t(1,:));

hold on

plot(x(2,:),t(2,:));

hold on

plot(x(3,:),t(3,:));

hold on

t1=zeros(1,4);Vn=zeros(1,4);T=zeros(3,41);

for i=1:4

s1=0;s2=0;


for j=1:i

t1(i)=2*h(j)/vp(j)+t1(i);

s1=h(j)*vp(j)+s1;

s2=h(j)/vp(j)+s2;


end

Vn(i)=s1/s2;


end

figure (2);

for k=2:4

for i=0:40

T(k-1,i+1)=sqrt(x(k-1,i+1)*x(k-1,i+1)/Vn(k-1)+t1(k-1)*t1(k-1 ));


end

plot(x(k-1,:),T(k-1,:));

hold on

end

 

2.S波

n=3;

h=[5 2 5 3];i1=zeros(n);

vs=[0 3 4 4.7];

p=zeros(n);

x=zeros(n-1,91);t=zeros(n-1,91);

for i=0:90

i1(2)=i;


for j=2:3

p(j)=sind(i1(j))/vs(j);

i1(j+1)=asind(p(j)*vs(j+1));


end


for kk=2:n


for k=2:kk

x(kk-1,i+1) =2*h(k)*tand(i1(k))+x(kk-1,i+1);

t(kk-1,i+1) =2*h(k)/vs(k)/cosd(i1(k))+t(kk-1,i+1);


end


end


end

plot(x(1,:),t(1,:));

hold on

plot(x(2,:),t(2,:));

hold on

hold on

t1=zeros(1,4);Vn=zeros(1,4);T=zeros(3,91);

for i=2:3

s1=0;s2=0;


for j=2:i

t1(i)=2*h(j)/vs(j)+t1(i);

s1=h(j)*vs(j)+s1;

s2=h(j)/vs(j)+s2;


end

Vn(i)=s1/s2;


end

figure (2)

hold on

for k=3:4


for i=0:90

T(k-2,i+1)=sqrt(x(k-2,i+1)*x(k-2,i+1)/Vn(k-1)+t1(k-1)*t1(k-1));


end

plot(x(k-2,:),T(k-2,:));

end

 

 

三、结果及分析





结果分析:当x较小时,精确公式和双曲线近似的结果是非常近似的。随着x的增大,精确公式和双曲线近似公式计算出来的结果的差距越来越大,可见当x较小时,可以使用双曲线公式近似替代精确公式。随着层数的增加,垂向走时越来越长,当x增大到一定的时候,所有界面的反射波走时曲线会有一个交点,这是因为深度越大波速度越大的原因,
虽然下层反射波传播距离更大,但是时间随着x增大到一定程度会小于上层反射波的走时。

标签: 反射波 地震学 地震学实验 地震学程序 时距曲线 海洋地壳模型
最后更新:2016年1月28日

小菜菜

菜鸟

打赏 点赞
< 上一篇
下一篇 >

文章评论

razz evil exclaim smile redface biggrin eek confused idea lol mad twisted rolleyes wink cool arrow neutral cry mrgreen drooling persevering
取消回复

此站点使用Akismet来减少垃圾评论。了解我们如何处理您的评论数据。

标签聚合
虚拟机 宝塔面板 c语言 模板 地震学程序 ppt wordpress onedrive
最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
曲阜东站--济宁的枢纽高铁站 matlab绘制双对数曲线 最新公式编辑器mathtype分享 阿里云ecs带宽详解 知更鸟begin1.8破解版主题分享 考研英语复试常考话题答题技巧及范文--谈谈你的大学

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2