山野莽夫

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

水平多层介质反射波时距曲线地震检波器组合方向频率特性曲线

2016年2月2日 3964点热度 0人点赞 0条评论

 一、实验目的  
         通过编制C语言程序,了解水平层状介质情况下反射波时距曲线、以及均方根速度和平均速度近似的特点和差异。理解地震检波器组合的方向特性曲线、频率特性曲线。  
         二、主要内容  
         1、二层参数:V2 = 1500m/s,h2 = 700m。用精确方程计算反射波时距曲线,用均方根速度近似的双曲线方程计算反射波时距曲线,用平均速度近似的双曲线方程计算反射波时距曲线,将三种计算结果绘制在同一张图上,并对计算结果进行分析  
         2、编制地震检波器简单线性组合的方向特性曲线C或matlab计算程序,将检波器组合个数分别为2,3,4的方向特性曲线绘制在同一张图上(参考图见附图1),并对计算结果进行分析。  
         3、编制地震检波器简单线性组合的频率特性曲线C或matlab计算程序,将组内相邻检波器时差分别为0.01、0.002、0.005的频率特性曲线绘制在同一张图上(其中检波器组合个数为7,频率范围为0~200Hz,参考图见附图2),并对计算结果进行分析。  
         三、实验报告内容  
         本实验结束后要求写出实验报告,主要包括以下内容:  
         1、  实验目的及实验内容  
         2、  基本原理阐述  
         3、  实验结果分析  
         4、  程序及其附图  
         5、  实验体会  

一、试验目的    

了解水平层状介质情况下反射波时距曲线、以及均方根速度和平均速度近似的特点和差异。理解地震检波器组合的方向特性曲线、频率特性曲线。

二、实验原理    

1.水平多层反射波时距曲线

.对于一个n层水平水平状模型,一射线倾斜入射,在第n层A点反射后到达观测点S的时间为

由斯奈尔定理可以确定每层的入射角度和透射角度和反射角度。


炮检距可以表示为


使用这两个公式分别计算,绘出时距曲线。

当x较小的时候,可以用平均速度和均方根速度来等效走时曲线,




2.组合方向特性曲线    

当波的主频固定的时候,组合后的波的振幅至于波的射线方向有关,我们称与射线方向有关的组合振幅特性为组合方向特性。归一化的的组合振幅特性表达式为

3.组合频率特性曲线    

将波传播方向改为时差,即可将

三、实验内容    

1、二层参数:V2 = 1500m/s,h2 = 700m。用精确方程计算反射波时距曲线,用均方根速度近似的双曲线方程计算反射波时距曲线,用平均速度近似的双曲线方程计算反射波时距曲线,将三种计算结果绘制在同一张图上,并对计算结果进行分析

2、编制地震检波器简单线性组合的方向特性曲线C或matlab计算程序,将检波器组合个数分别为2,3,4的方向特性曲线绘制在同一张图上(参考图见附图1),并对计算结果进行分析。

3、编制地震检波器简单线性组合的频率特性曲线C或matlab计算程序,将组内相邻检波器时差分别为0.01、0.002、0.005的频率特性曲线绘制在同一张图上(其中检波器组合个数为7,频率范围为0~200Hz,参考图见附图2),并对计算结果进行分析。

四、实验结果    

1.水平多层时距曲线

2.组合检波方向特性曲线

  1. 组合检波频率特性曲线    

五、结果分析

1.水平多层反射波时距曲线

水平多层反射波使用三个公式计算时距曲线,用均方根速度和平均速度还是和精确公式还是有细微的误差,不过误差并不大,当炮检距比较小的时候,误差非常小,近似为双曲线,可以用这两个速度近似代替。平均速度计算的最大,平方根速度计算的次之,精确公式的最小。

2.组合方向特性曲线

组合方向特性曲线有明显规律性,对于高视速度的规则波,近乎垂直到达个检波点,相邻检波器之间的时差趋于零,达到最大值。组合数目多少对特性曲线也有影响,组合数目增加,通放带边界点左移,压制带内的极值降低。高视速度的波可以落入通放带,,组合输出波的振幅得到加强,低视速度的波组合后受到压制。

3. 频率特性曲线

从频率特性曲线上看,波的视速度趋于无穷时,组合后对所有频率成分都没有滤波作用,随着△t的增大,频率特性曲线通放带和压制带越来越明显,并且通放带变窄,表明组合有频率滤波作用,对于高频成分有压制作用。

六、源代码    

1.水平多层时距曲线(C语言)

#include<stdio.h>

#include<math.h>

#include<stdlib.h>

#define
dzkt5
"file.dat"

 

int main()

{

    float x[41], t[41],t1[41],t2[41], p, h[100], v[100],v1,v2,s1,t0,s2,s3,s4,ni,ii[100],d1=0.0,d2=0.0;int i, j, m;

    FILE *fp;

    printf("请输入层数\n");

    scanf("%d", &m);

    //printf("请依次输入每层的层厚和波速,用回车或空格隔开");

    for (i = 0;i<m;i++)

    {

        printf("请输入第%d层层厚和波速,空格或回车隔开\n",i+1);

        scanf("%f%f", &h[i], &v[i]);

    }

    fp=fopen (dzkt5, "a+");//文件打开

    for (i = 0;i<=40;i++)

    {

        t[i] = 0.0;

        x[i] = 0.0;

        s1 = 0.0;//初值

        s2 = 0.0;

        s3 = 0.0;//初值

        s4 = 0.0;//初值

        t0 = 0.0;

                for (j = 0;j<m;j++)

        {

            p = sin(i / 180.0 * 3.141592653) / v[j];

            

                ii[j] = asin(v[j]/v[0]*sin(i/180.0*3.141592653));

                                            t[i] = 2 *h[j] /v[j]/ cos(ii[j])+t[i];

                x[i] = 2 * h[j] * tan(ii[j])+x[i];

        
 

    
 

            s1 = h[j] / v[j] * v[j] * v[j] + s1;//均方根速度

            s2 = h[j] / v[j] + s2;

            t0 = 2 * h[j] / v[j] + t0;

            s3 = h[j] + s3;

            s4 = h[j] / v[j] /+ s4;

            

        }//累加计算

        

        v1 = pow((s1 / s2),0.5);//均方根速度

        v2 = s3 / s2;//平均速度

        

        t1[i] = sqrt(t0*t0 + x[i] * x[i] / v1 / v1);//均方根速度计算时间

        t2[i] = sqrt(t0*t0 + x[i] * x[i] / v2 / v2);//平均速度计算时间

        

        printf("%.6f\t%.6f\t%.6f\t%.6f\n", x[i], t[i],t1[i],t2[i]);

        fprintf(fp, "%.6f\t%.6f\t%.6f\t%.6f\n", x[i], t[i],t1[i],t2[i]);//存储到文件中

    }

    
 

        fclose(fp);//文件关闭

    return 0;

 

}

2.组合检波方向特性曲线

T=0:0.001:1;

n=2;

fai=(sin(n*3.14*T))./(n*sin(3.14*T));

plot(T,abs(fai),'r');

hold on;

n=3;

fai=(sin(n*3.14*T))./(n*sin(3.14*T));

plot(T,abs(fai),'--');

hold on;

n=4;

fai=(sin(n*3.14*T))./(n*sin(3.14*T));

 

plot(T,abs(fai),':k');

hold off;

 

3.组合检波频率特性曲线

n=7;

f=0:1:200;

t=0.01;

fai=(sin(n*3.14*f*t))./(n*sin(3.14*f*t));

plot(f,abs(fai),'r');

hold on;

t=0.002;

fai=(sin(n*3.14*f*t))./(n*sin(3.14*f*t));

plot(f,abs(fai),'--');

hold on;

t=0.005;

fai=(sin(n*3.14*f*t))./(n*sin(3.14*f*t));

plot(f,abs(fai),'k:');

hold off

 

 

 

这个实验首先做了水平多层反射波时距曲线的绘制,水平多层反射波水平单界面反射波有着相似之处,我在这个任务中编制了了一个通用型程序,这个程序适合更多层底层模型反射波时距曲线的分析。另外本次试验做了组合检波的方向特性曲线和频率特性曲线,这两种曲线分析了组合检波是一个滤波器,有滤波作用,可以压制噪声,增强有用信号,提高信噪比。这个小程序用matlab编译,代码比较简单。

 

 

 

学生(签名):

2015 年 月 日

 

 

 

 

 

 

 

成绩评定:

指导教师(签名):

年 月 日


 

标签: c语言 C语言小程序 地震勘探 方向特性曲线 时距曲线 检波器 检波器组合 水平多层 频率特性曲线
最后更新:2016年2月2日

小菜菜

菜鸟

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

文章评论

  • 勒煜鑫

    水平层折射波时距方程C语言编程可以分享一下吗

    2022年6月8日
    回复
  • 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工具分享
    考研英语复试常考话题答题技巧及范文--谈谈你的大学 ubuntu16.04开启root账户登录 简单网站源码(应付建站公选课作业) 莱卡lycamobile购买、激活、充值、关闭语音信箱、使用、复活等相关事宜 明朝皇帝 盘点一下云服务大厂CDN和对象存储的免费额度
    标签聚合
    虚拟机 onedrive 地震学程序 c语言 ppt wordpress 模板 宝塔面板
    最近评论
    小菜菜 发布于 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