山野莽夫

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

地震勘探时距曲线绘制

2016年5月23日 5912点热度 2人点赞 0条评论


试验要求:

一、实验目的

通过编制和运行C或Matlab程序,进一步理解地震反射波、折射波和直达波以及绕射波时距曲线方程及其特点。

二、主要内容

1、假设地下有一个水平地质界面,埋藏深度100m,地层是均匀的各向同性介质,地震波纵波传播速度为1500m/s,请编制C或Matlab程序,绘制出中点放炮共炮点的反射波、直达波和折射波时距曲线(参考附图1)。Matlab程序编制中,炮检距:-500m~500m,间距10m;临界角取35°。

2、假设上述水平界面模型中界面埋藏深度变为1000m,上述时距曲线有何变化,绘出相应图形,并从理论上加以分析;

3、假设地下有一个倾斜界面,地层倾角为20°,地层是均匀的各向同性介质,地震波纵波传播速度为1500m/s,请编制C或Matlab程序,绘制出端点放炮共炮点的反射波、直达波和折射波时距曲线(参考附图2)。C或Matlab程序编制中,炮点处界面的法向深度为100m;炮检距:-500m~500m,间距10m;临界角取35°。

4、假设地下有一个绕射点,埋藏深度为1000m,绕射点以上地震波传播速度为1500m/s,编制C或Matlab程序计算并绘制出绕射波时距曲线图。C或Matlab程序编制中,炮检距:-1000m~2000m,间距20m;绕射点在地面的投影点距离炮点200m(参考附图3)

三、实验报告内容

本实验结束后要求写出实验报告,主要包括以下内容:

  1. 实验目的及内容
  2. 基本原理阐述
  3. 实验结果分析
  4. 程序及其附图
  5. 实验体会


实验报告:

一、实验目的

通过编制和运行Matlab程序,并绘制出时距方程曲线。进一步理解地震反射波、折射波和直达波以及绕射波时距曲线方程及其特点。

二、实验原理

水平界面的反射波时距方程为

可见这是一个原点为对称的标准双曲线。

直达波的时距方程为

折射波的时距方程为 故折射波的时距曲线为一条直线,斜率为下层速度的倒数。

倾斜界面的反射波时距方程为 故也为标准双曲线。

直达波方程为

上倾折射波方程为

绕射波时距曲线方程为

故绕射波时距曲线也双曲线。

三、试验内容

1、水平界面时距曲线(埋深100米)

假设地下有一个水平地质界面,埋藏深度100m,地层是均匀的各向同性介质,地震波纵波传播速度为1500m/s,请编制C或Matlab程序,绘制出中点放炮共炮点的反射波、直达波和折射波时距曲线(参考附图1)。Matlab程序编制中,炮检距:-500m~500m,间距10m;临界角取35°。

2、水平界面时距曲线(埋深1000米)

假设上述水平界面模型中界面埋藏深度变为1000m,上述时距曲线有何变化,绘出相应图形,并从理论上加以分析;

3、倾斜界面时距曲线

假设地下有一个倾斜界面,地层倾角为20°,地层是均匀的各向同性介质,地震波纵波传播速度为1500m/s,请编制C或Matlab程序,绘制出端点放炮共炮点的反射波、直达波和折射波时距曲线(参考附图2)。C或Matlab程序编制中,炮点处界面的法向深度为100m;炮检距:-500m~500m,间距10m;临界角取35°。

4、绕射波时距曲线

假设地下有一个绕射点,埋藏深度为1000m,绕射点以上地震波传播速度为1500m/s,编制C或Matlab程序计算并绘制出绕射波时距曲线图。C或Matlab程序编制中,炮检距:-1000m~2000m,间距20m;绕射点在地面的投影点距离炮点200m(参考附图3)

四、试验结果

1.水平界面反射波直达波折射波时距曲线(埋深100米)

2.水平界面反射波直达波折射波时距曲线(埋深1000米)


3.倾斜层反射波直达波折射波时距曲线

4.绕射波时距曲线


五、结果分析

1.水平界面的时距曲线为以炮点对称的标准双曲线,双曲线的极小点为法相反射时间,直达波时距曲线是反射波时距曲线的渐近线。水平界面折射波的时距曲线是一条直线,直线的斜率是下层速度的倒数,下层速度越高,折射波斜率越小,折射波时距曲线越平缓,反之越陡。盲区内不存在折射波,所以折射波的时距曲线可以说是认为的把曲线延长到坐标轴。

2.把埋深为100米和埋深1000米的水平层时距曲线作对比,可以发现埋深越大,反射波直达波折射波时距曲线越平缓。

3.单倾斜层反射波时距曲线是一个标准的双曲线,双曲线的极值点的值是最小时间,反射波时距曲线极值电位于激发点上倾一侧,双曲线渐近线不再是直达波时距曲线,斜率仍然是下层速度的倒数。折射波仍然是一条直线。

4.绕射波时距曲线仍为双曲线,移动激发点绕射波的时距曲线形状不变,随着激发点原理坐标原点,绕射波极小点时间增大。

六、源代码(matlab)

1.水平界面反射波直达波折射波时距曲线(100米和1000米埋深)

h=input('input h=');%埋深

V1=input('input V1=');%波速

i=input('input i='); %入射角度

x1=input('input x1=');%坐标左断点

x2=input('input x2='); %坐标右端点

x=x1:x2;

t=1/V1*sqrt(power(x,2)+4*h*h);%反射波

plot(x,t,'r');

hold on;

t1=abs(x/V1);%直达波

plot(x,t1,'k--');

hold on

V2=V1/sind(i);

t2=(abs(x)/V2+2*h*cosd(i)/V1);%折射波

plot(x,t2,'b:');

hold on

2.倾斜层反射波直达波折射波时距曲线

fai=20;%倾角

h=100;%埋深

v1=1500;%波速

xm=-2*h*sind(fai);

x=-500:500;

t=1/v1*sqrt(power((x+500),2)+4*h*h-2*(x+500)*xm);%反射波

plot(x,t,'r');

hold on

t1=abs(x+500)/v;%直达波

plot(x,t1,'k--');

hold on

t2=(x+500)*cosd(fai)*sind(i)/v1+(2*h+(x+500)*sind(fai))/v1*cosd(i);折射波

hold on

plot(x,t2,'b:');

3.绕射波时距曲线

h=input('input h='); 埋深

v=input('input v='); 波速

d=input('input d='); 距离

x1=input('input x1=');左端点

x2=input('input x2='); 右端点

x=x1:20:x2;

% x=-1000:20:2000;

% v=1500;d=200;h=1000;d=200;

t=sqrt(d*d+h*h)/v+sqrt(power((x-d),2)+h*h)/v;

plot(x,t);%走时

本次试验主要做的是水平界面、倾斜界面反射波直达波和折射波的走时曲线,还做了绕射波的反射波时距曲线,绘制出了图线,并且做了分析,

认识了时距曲线的特征。本次试验使用matlab编程,matlab在编程绘图方和C语言相比有着自己的有优点,而且程序比简单,适合小数据量的分析。

学生(签名):

2015 年 月 日

成绩评定:

指导教师(签名):

年 月 日

标签: 反射波 地震勘探 折射波 时距曲线 直达波 绕射波
最后更新:2016年5月23日

小菜菜

菜鸟

打赏 点赞
< 上一篇

文章评论

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工具分享
四川电信校园宽带拨号问题 贪吃蛇C语言代码 数据结构(c语言版)所有源代码 浦科特固态硬盘m7vc 宝塔docker管理器的使用(以docker搭建临时邮件forsaken-mail为例) 地震学和地球结构(震源追踪)
标签聚合
宝塔面板 onedrive 模板 地震学程序 虚拟机 c语言 ppt wordpress
最近评论
小菜菜 发布于 6 个月前(11月24日) 这玩意已经废了,成收割工具了,不能再用了。
eamon 发布于 7 个月前(11月07日) 我一年不用了才发现这个休眠管理费每月15,一共扣了我135元,然后我消费还消费不了,我宁愿消费掉也不...
magic 发布于 11 个月前(07月03日) 请问账号不注销会有什么影响吗?
magic 发布于 11 个月前(07月01日) 我想问一下 如果不注销账号就留着会怎么样
qwp6601 发布于 12 个月前(06月04日) 有没有方法改为bing

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2