山野莽夫

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

斯奈尔定律和Zoeppritz方程

2016年1月30日 7442点热度 2人点赞 2条评论

 一、实验目的

1.利用密度、上下界面的纵横波速度通过斯奈尔定律求出该界面的0°~90°入射角的反射角度和透射角度。

2.利用Zoeppritz方程绘制反射和透射能量系数曲线。

二、实验原理

1. 斯奈尔定律反映了在弹性分界面上入射波和反射波和透射波射线之间的角度关系。


当入射波是p波和sv波的时候还会波形转换,分别产生反射透射sv波和反射透射p波,当SH波入射时,无波形转换,只有反射SH波和透射SH波。所有波都会满足斯奈尔定律。通过斯奈尔定律,在我们知道相邻两层波形参数和入射角度之后我们就可以求得入射波的反射p波sv波和透射p波sv波的角度。

2. 入射波遇到弹性分界面,波就会产生反射波和透射波,当然这些波的角度都符合斯奈尔定理。这样的弹性分界面将会由一个波变成多个波,波的能量也会发生变化,zoeppritz方程就可以确定波的转换和能量分配关系。本实验报告主要讨论p波入射。

P波入射的zoeppritz方程如下


这样就可以确定入射波和反射波透射波及转换波的振幅比,

为反射P波P11的反射系数;为反射S波P1S1的反射系数;

为透射P波P12的透射系数;为透射S波P1S2的透射系数;

得到这些数据以后直接用matlab绘图。

三、实验内容和步骤

模型如下

 

用斯奈尔定律解出0~90度对应入射角的反射p波和反射s波的角度,和透射p波及透射s波的角度,分别代入zoeppritz方程计算出反射透射系数并汇出能量系数曲线。

四、实验结果


 

五、实验分析:

随着透射角度不断增大,反射p波的反射系数不断减小,透射p波的透射系数不断增大,转换波反射sv波的反射系数变化也不大在0度和90度的时候都为0。转换波透射sv波的的透射系数变化也不大,在0度和90度的时候都为0;透射系数总大于0,所以透射总是存在。

六、源代码(本次实验使用matlab)

 

vp1=3300;vs1=1585;p1=2.4;vp2=3100;vs2=1989;p2=2.24;

s=zeros(4,90);

d=zeros(4,91);


for i=0:90

a1=sind(i);

a2=vs1/vp1*sind(i);

b1=vp2/vp1*sind(i);

b2=vs2/vp1*sind(i);%斯奈尔定律求反射透射角度

b=[-a1;sqrt(1-a1*a1);sind(2*i);2*a2*a2-1];%列向量B

A=zeros(4);

A(1,1)=a1;

A(1,2)=sqrt(1-a2*a2);

A(1,3)=-b1;

A(1,4)=-sqrt(1-b2*b2);

A(2,1)=sqrt(1-a1*a1);

A(2,2)=-a2;

A(2,3)=sqrt(1-b1*b1);

A(2,4)=-b2;

A(3,1)=sind(2*i);

A(3,2)=vp1/vs1*(1-2*a2*a2);

A(3,3)=p2/p1*vs2*vs2*vp1/vs1/vs1/vp2*2*b1*sqrt(1-b1*b1);

A(3,4)=p2/p1*vp1*vs2/vs1/vs1*(1-2*b2*b2);

A(4,1)=1-2*a2*a2;

A(4,2)=-vs1/vp1*sind(2*i);

A(4,3)=-p2/p1*vp2/vp1*(1-2*b2*b2);

A(4,4)=p2/p1*vs2/vp1*2*b2*sqrt(1-b2*b2);%矩阵A赋值

x=linsolve(A,b);%求解线性方程组

d(1:4,(i+1))=x;


end

g=d';

plot(0:90,g);%绘图

 

 

 

学生实验 心得

通过做这个实验,认识到了斯奈尔定律用来计算一定角度入射波的反射P波和反射SV波的反射角度,和透射P波和透射SV波的投射角度。并用Zoeppritz方程编写了程序来计算不同角度入射波反射能量分配和透射能量分配关系。并绘制了图线,清楚的认识到了P波入射,在界面上反射和透射关系。

 

 

 

学生(签名):

2015 年 月 日

指导

教师

评语

 

 

 

 

 

 

 

成绩评定:

指导教师(签名):

年 月 日


 

标签: zoeppritz方程 反射系数 地震勘探 振幅比 斯奈尔定律 能量分配 转换波 透射系数
最后更新:2016年1月30日

小菜菜

菜鸟

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

文章评论

  • 猫野

    您好,请问一下,如果是SV波入射时,又该如何求解左布里兹方程呢

    2021年5月4日
    回复
    • 小菜菜

      @猫野 抱歉,多年不学了,忘干净了

      2021年5月10日
      回复
  • 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工具分享
    认识测井软件--CifLog 地震断层P轴和T轴计算 关于 星愿浏览器Twinkstar便携版制作方法以及便携版成品下载地址 cloudflare简单使用测评以及如何自定义节点(例如使用香港节点) 瑞士免费不限速网盘Degoo注册就送105G免费空间(可做任务无限扩容)
    标签聚合
    虚拟机 地震学程序 wordpress c语言 宝塔面板 onedrive ppt 模板
    最近评论
    小菜菜 发布于 7 个月前(11月24日) 这玩意已经废了,成收割工具了,不能再用了。
    eamon 发布于 8 个月前(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