山野莽夫

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

地震学和地球结构(震源追踪)

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

C-2. (a) 编写程序:震源在地表,选择震源一定范围的入射角,追踪穿过一个地球模型的射线。图形的输出显示震源、射线路径、地球表面、核幔边界以及内核-外核边界。

(b) 用PREM模型追踪震源地表的地震射线,表示出上地幔间断面和地核对射线的影响。

(c) 该程序生成一个走时曲线图,能根据这个图给出上地幔间断面吗?

一、实验原理

运用地球的射线参数公式求取射线路径的最低点。

运用圆的几何角度关系,和角度公式计算每层的入射点

走时计算公式

读入prem模型的速度和深度数据,计算出一组数据,运用matlab绘图。

二、程序源代码

射线追踪程序

s=xlsread('perm.xlsx');

r=s(:,1);

v=s(:,2);

for i=1:26

zhong=r(i,1);r(i,1)=r(54-i,1);r(54-i,1)=zhong;

zhong=v(i,1);v(i,1)=v(54-i,1);v(54-i,1)=zhong;

end

for jiao=0:0.01:0.1

% jiao=input('请输入入射角');%入射角

p=r(1)*sin(jiao)./v(1);

N=size(r);

for i=2:N

k=v(i)*p./r(i);

if k>=1;

if k==1

M=i

Q=0;break;

else if r(i)==r(i-1)

M=i-1

R=r(M,1);

V=v(M,1);

y=R/V;

Q=p/(R*(y^2-p^2))

break;

else Xl=(v(i)-v(i-1))./(r(i)-r(i-1));

r(i)=(Xl*r(i-1)-v(i-1))./(Xl-1/p);

v(i)=r(i)/p;

M=i

Q=0

break;

end

end

else continue;

end

end

r(M,1)=int64(r(M,1));

t=0;

T=0;

zb1=zeros(1,6371);

zb2=zeros(1,6371);

for i=M:-1:2

K=r(i-1)-r(i);

for j=0:K

R=r(i)+j;

V=v(i)+j*(v(i-1)-v(i))./(r(i-1)-r(i));

y=R/V;

deta=p/(R*(y^2-p^2)^0.5);

T=deta+T;

zb1(j+t+1)=R*sin(T);

zb2(j+t+1)=R*cos(T);

end

t=K+1+t;

end

for i=1:6370

if zb1(i)==0 break

end

end

zb=zeros(2,i-1);

phi=i-1;

for j=1:i-1

zb(1,j)=zb1(j);

zb(2,j)=zb2(j);

end

zb3=zeros(2,2*phi);

for j=1:phi

zb3(1,phi+j)=zb(1,j);

zb3(1,phi+1-j)=-zb(1,j);

zb3(2,phi+j)=zb(2,j);

zb3(2,phi+1-j)=zb(2,j);

end

plot(zb3(1,:),zb3(2,:));hold on;

len=[6371 3630 1221.5];%边界

x=-6371:6371;

for i=1:3

y=sqrt(len(i).^2-x.^2);

plot(x,y);hold on;

end

end

走时曲线程序

s=xlsread('perm.xlsx');

r=s(:,1);

v=s(:,2);

Sj=zeros(2,50);

for i=1:26

zhong=r(i,1);r(i,1)=r(54-i,1);r(54-i,1)=zhong;

zhong=v(i,1);v(i,1)=v(54-i,1);v(54-i,1)=zhong;

end

xishur=r;

xishuv=v;

for g=1:50

r=xishur;

v=xishuv;

jiao=(51-g)*0.02;

p=6371*sin(jiao)/1450;

N=size(r);

for i=2:N

k=v(i)*p./r(i);

if k>=1;

if k==1

M=i;

Q=0;break;

else if r(i)==r(i-1)

M=i-1;

R=r(M,1);

V=v(M,1);

y=R/V;

Q=p/(R*(y^2-p^2));

break;

else Xl=(v(i)-v(i-1))./(r(i)-r(i-1));

r(i)=(Xl*r(i-1)-v(i-1))./(Xl-1/p);

v(i)=r(i)/p;

M=i;

Q=0;

break;

end

end

else continue;

end

end

r(M,1)=int64(r(M,1));

T=0;

D=0;

for i=M:-1:2

K=r(i-1)-r(i)

for j=0:K

R=r(i)+j;

V=v(i)+j*(v(i-1)-v(i))./(r(i-1)-r(i));

y=R/V;

deta=p/(R*(y^2-p^2)^0.5);

time=y^2/(R*(y^2-p^2)^0.5);

D=deta+D;

T=time+T;

end

end

Sj(1,g)=D/pi*180;

Sj(2,g)=T;

end

plot(Sj(1,:),Sj(2,:));

  1. 结果及分析

结果分析,此图在入射角为0.04的情况下得出。因为速度是连续变化的,所以绘制出图线为曲线。因为上地幔间断面为速度间断面,所以使射线弯折。入射角越小,地震波穿越地球越深。当入射角为0度的时候,地震波射线直接穿过地核,到达对称位置,入射角越大,地震波穿越深度越浅,深度很小时就到达最低点后返回。上地幔间断面和核幔边界因为速度突然减小,所以图线会下凹。核幔边界则会偏向地核。在内外核边界速度增大,本身图线应该下凹,但是下凹程度减小,近似为直线。走时曲线也为曲线,在弯折处为速度突变面。所以第一个弯折出为上地幔间断面。

标签: PERM模型 地震学 地震学实验 地震学程序 震源追踪
最后更新: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来减少垃圾评论。了解我们如何处理您的评论数据。

联系方式

QQ群 | TG群 | 邮箱

最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
探讨一下谷歌云的路线问题 三度体球体磁场正演 Deepin(linux)操作系统如何使用windows远程桌面 QQ截图配置方法 地震学和地球结构(震源追踪) excel中逆序一列的数据
标签聚合
宝塔面板 onedrive c语言 地震学程序 虚拟机 ppt wordpress 模板
最近评论
小菜菜 发布于 7 个月前(11月24日) 这玩意已经废了,成收割工具了,不能再用了。
eamon 发布于 7 个月前(11月07日) 我一年不用了才发现这个休眠管理费每月15,一共扣了我135元,然后我消费还消费不了,我宁愿消费掉也不...
magic 发布于 12 个月前(07月03日) 请问账号不注销会有什么影响吗?
magic 发布于 12 个月前(07月01日) 我想问一下 如果不注销账号就留着会怎么样
qwp6601 发布于 1 年前(06月04日) 有没有方法改为bing

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2