山野莽夫

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

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

2016年1月28日 2896点热度 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来减少垃圾评论。了解我们如何处理您的评论数据。

标签聚合
宝塔面板 ppt 虚拟机 模板 onedrive 地震学程序 c语言 wordpress
最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
关于全球付globalcash万事达虚拟信用卡的安全和使用问题的提示和下车方法 Deepin+宝塔面版配置Hexo静态博客 地震学和地球结构(海洋地壳模型反射波时距曲线) 阿里云国际站跃上云端新春采购季活动新老用户都可领代金券(文末有福利)! 关于 MP3音频文件格式详细解析

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2