matlab segy文件,利用matlab实现segy格式数据的读写研究和分析.ppt

利用matlab实现segy格式数据的读写研究和分析

2.SEGY数据格式 地震数据一般以地震道为单位进行组织,采用SEG-Y文件格式存储。SEG-Y格式是由SEG (Society of Exploration Geophysicists)提出的标准磁带数据格式之一,它是石油勘探行业地震数据的最为普遍的格式之一。 标准SEG-Y文件一般包括两部分, 1、文件头总长度为3600字节,分两部分。 (1)文件头第一部分 长度: 3200bytes ;组成: 80bytes*40;特性 :EBCDIC字符集,参数卡,需要转换为ASCII码后才能显示。 (2)文件头第二部分 长度:400bytes;数据类型:32位、16位的整型;特性:二进制头,记录数据体信息。 2、数据体由多个数据道组成。每道数据分两部分:道头、采样数据。 (1)数据道道头 长度:240 bytes;数据类型:32位、16位的整型;特性:记录采样点数、采样间隔、CDP号、XLine号、Line号以及坐标信息等。 (2)采样数据 长度:采样数*采样点字节数; 5.数据介绍 文件11102.sgy是一个浅层的单边放炮,一边接收,SEGY格式的地震数据。总共有24道,采样数是4096,采样间隔是100us,采样数据类型是4个字节IEEE Float(32 bit)。 读取成图 读取彩图 原始数据图 原始数据图 7.结论认识 通过以上比较,程序读取成图能够真实的反应原始数据图,证实了该方法可有效读取原始数据。 地震勘探采集的地震资料经处理后多以SEGY格式记录。因而用来进行解释的地震数据也基本来源于SEGY格式文件。要了解探区的地质构造,准确估算探区石油储量!,确定井位,都需要对大量的地震数据进行精确的解释, 这样,准确读取SEGY格式地震文件,便显得至关重要。 利用matlab实现SEGY格式数据的读写只是用来理论研究,不能用于实际生产,要想在实际生产中运用,还必须调试成C语言程序或其他程序。 六.程序存在的问题 ⑴ 此程序只对于标准的SEGY格式数据能够进行读写; ⑵ 程序没有解决地震数据二维与三维的区别; ⑶ 没有形成一个比较智能系统化的函数。 课题名称: 利用matlab实现SEGY格式数据的读写研究与分析 专业班级: 物探11102 学生姓名: 马文博 指导老师:毛宁波教授 报告 提纲 一、课题来源 二、研究的目的和意义 三、主要研究内容 四、国内外现状和发展趋势 五、技术思路及成果图 六、程序存在的问题 一、课题来源 利用matlab实现SEGY格式数据的读写研究与分析 ——源自导师的研究项目 二、研究的目的和意义 1、学会使用matlab软件。 2、深入了解segy格式数据的具体内容 目的 意义 1、理论与实践相结合。 2、通过编程,提高动手能力 3、用matlab编辑segy格式数据读写程序 三、主要研究内容 (1)检索有关的文献,收集相关资料,完成开题报告; (2)学会使用matlab软件; (3)深入了解segy格式数据的具体内容; (4)用matlab编辑segy格式数据读写程序; (5)归纳总结并撰写毕业论文 四、国内外现状和发展趋势 SEG Y数据交换格式SEG Y rev 0自1975年推出以来在地球物理界得到了广泛的应用,并引发了一系列的革新。随着三维数据采集技术以及高速度、大容量记录媒体的应用,SEG Y rev 0已经不能满足现阶段数据采集、处理及存储的需求,因此,SEG技术标准委员会(the SEG Technical Standards Committee)推出了新的格式标准——-sEG Y rev 1.0。 1930 年 SEG (勘探地球物理家学会) 1975 年 SEG Y数据交换格式SEG Y rev 0推出 之 后 SEG Y rev 1.0推出 五,技术路线及成果图 matlab软件 安装和使用软件 系统的了解matlab语言 简单的编程 segy数据格式 掌握segy数据的具体内容 利用matlab编程实现segy数据的读写 利用读取的数据进行简单的处理 成果图 1,技术思路 注:此图是我做C语言课程设计时指导老师陈义群老师所给. 400字节二进制文件头 字节 描述 3201-3204 作业标识号 3205-3208 测线号,对3-D叠后数据而言,它将典型地包含纵向测线(In-line)号 3209-3212 卷号 3213-3214 每个道集的数据道数,叠前数据强制要求 3215-3216

  • 0
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是一个简单的使用卷积神经网络恢复segy格式地震数据Matlab代码: ``` %% 加载数据 data = readSegy('data.segy'); [ntr,ns] = size(data); %% 数据预处理 data = (data-mean(data(:))) ./ std(data(:)); data = reshape(data,[ntr,1,ns]); %% 定义卷积神经网络模型 inputSize = [ntr 1 ns]; numFilters = 32; layers = [imageInputLayer(inputSize) convolution2dLayer([1,5],numFilters) reluLayer maxPooling2dLayer([1,2],'Stride',[1,2]) convolution2dLayer([1,5],numFilters*2) reluLayer maxPooling2dLayer([1,2],'Stride',[1,2]) convolution2dLayer([1,5],numFilters*4) reluLayer transposedConv2dLayer([1,4],numFilters*2,'Stride',[1,2],'Cropping',[0,1,0,1]) convolution2dLayer([1,5],numFilters*2) reluLayer transposedConv2dLayer([1,4],numFilters,'Stride',[1,2],'Cropping',[0,1,0,1]) convolution2dLayer([1,5],1) regressionLayer]; %% 训练模型 options = trainingOptions('sgdm','MaxEpochs',10,'InitialLearnRate',0.01); net = trainNetwork(data,data,layers,options); %% 使用模型恢复数据 recon = predict(net,data); recon = reshape(recon,[ntr,ns]); recon = (recon.*std(data(:))) + mean(data(:)); %% 显示恢复的数据 figure;imagesc(data(:,:,1));colormap(gray);title('原始数据'); figure;imagesc(recon(:,:,1));colormap(gray);title('恢复数据'); ``` 其中 `readSegy` 函数可以从segy格式文件取地震数据,需要自己实现或者从网络上下载。这个代码使用一个简单的卷积神经网络模型来恢复数据,可以根据具体情况进行调整和改进。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值