MATLAB-SWPU第七章-地震处理相关知识.ppt

《MATLAB-SWPU第七章-地震处理相关知识.ppt》由会员分享，可在线阅读，更多相关《MATLAB-SWPU第七章-地震处理相关知识.ppt（60页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、7.1 7.1 地震数据格式地震数据格式 常见的地震数据体有常见的地震数据体有SEGSEG格式和格式和NOSEG NOSEG 格式，格式， SEGSEG是（是（ The The Society of Exploration GeophysicistsSociety of Exploration Geophysicists）的缩写，）的缩写， 而而SEGSEG格格式有多种，其中最常用的就是式有多种，其中最常用的就是SEGYSEGY格式地震数据体文件。它有格式地震数据体文件。它有32003200个字节的图头（也叫卷头），个字节的图头（也叫卷头），400400字节的文件头，字节的文件头，240240

2、字节的字节的道头，道数据道头，道数据，.道头，道数据，很多个。图头信道头，道数据，很多个。图头信息我们一般不关心，不太常用。文件头里有些信息可以用到，例息我们一般不关心，不太常用。文件头里有些信息可以用到，例如数据格式了，样品数了，采样间隔了。而道头在实际加载的过如数据格式了，样品数了，采样间隔了。而道头在实际加载的过程中最常用，其中最常用的是程中最常用，其中最常用的是X X、Y Y、道号和线号。但在实际的数、道号和线号。但在实际的数据中它们存放的位置不一定相同，要具体数据体具体分析，可以据中它们存放的位置不一定相同，要具体数据体具体分析，可以用一个用一个SEGYSEGY格式分析器进行转换格式

3、分析器进行转换 7.1 7.1 地震数据格式地震数据格式SEGY数据由文件头和数据体组成数据由文件头和数据体组成文件头总长度为文件头总长度为3600字节，分两部分。字节，分两部分。文件头第一部分文件头第一部分长度：长度： 3200bytes ；组成：；组成： 80bytes*40；特性：；特性：EBCDIC字符集，字符集，参数卡，需要转换为参数卡，需要转换为ASCII码后才能显示。码后才能显示。文件头第二部分文件头第二部分长度：长度：400bytes；数据类型：；数据类型：32位、位、16位的整型；特性：二进制头，位的整型；特性：二进制头，记录数据体信息。记录数据体信息。数据体由多个数据道组成

4、。每道数据分两部分：道头、采样数数据体由多个数据道组成。每道数据分两部分：道头、采样数据，数据类型：据，数据类型：32位的浮点型位的浮点型 ，占，占4个字节。个字节。7.1 7.1 地震数据格式地震数据格式工作站SEGY数据存储格式有两种：IEEE和IBM IEEE和IBM的整型存储与微机格式的存储不同之处，IEEE和IBM的高字节在前、低字节在后，即BigEndian，微机则是低字节在前、高字节在后，即LittleEndian 。文件头格式文件头格式道头格式道头格式例 读取地震数据并绘图nclear allnindata=altreadsegy(mianbomodle.sgy);%读取数据n

5、Tr=length(indata(1,:);%道数nT=length(indata(:,1); %采样点数ns_wplot(indata);%绘剖面例 提取第10道单独显示，并且加入随机噪声nclear allnindata=altreadsegy(mianbomodle.sgy);nTr=length(indata(1,:);nT=length(indata(:,1); ndata=indata(:,10);%提取第10道来显示ntemp=randn(T,1);%生成随机数据ndata=data+temp;ns_wplot(data);n%s_cplot(data);%变面积显示例 在第10

6、道中加入随机噪声nclear allnindata=altreadsegy(mianbomodle.sgy);nTr=length(indata(1,:);nT=length(indata(:,1); ntemp=randn(T,1);%加入随机噪声nindata(:,10)=indata(:,10)+temp;ns_wplot(indata);n%s_cplot(data);%变面积显示7.2 常用处理-滤波nclear allnindata=altreadsegy(mianbomodle.sgy);nTr=length(indata(1,:);nT=length(indata(:,1);

7、ndata=indata(:,10);%提取第10道来显示ntemp=randn(T,1);%生成随机数据ndata=data+temp;nfdata=fft(data);nplot(abs(fdata);例 低通滤波器实现nnfdata=fft(data);nfT=length(fdata);nfor j=100:fT %低通滤波器n fdata(j)=0;nendnplot(abs(fdata);例 低通滤波器实现-反变换nnfdata=fft(data);nfT=length(fdata);nfor j=100:fT %低通滤波器n fdata(j)=0;nendnplot(abs(fd

8、ata);ndata2=ifft(fdata);%逆变换ns_wplot(data);ns_wplot(data2);例 高通滤波器实现nnfdata=fft(data);nfT=length(fdata);nfor j=1:200 %高通滤波器n fdata(j)=0;nendnplot(abs(fdata);例 带通滤波器实现nnfdata=fft(data);nfT=length(fdata);nfor j=1:fT %带通滤波器n if(j300|jltime(i)/tspace-twindow&jltime(i)/tspace-twindow&j= 0 & (nf-fend) = 0

9、 & (nf-fend) = 0 & (nf+fstart-1/tdeta) = 0 & (nf-fend) = 0 & (nf+fstart-1/tdeta) = 0 & (nf-fend) = 0 & (nf+fstart-1/tdeta) = 0 & (nf-fend) = 0 & (nf+fstart-1/tdeta) = 0 & (nf-fend) = 0 & (nf+fstart-1/tdeta) = 0 & (nf-fend) = 0 & (nf+fstart-1/tdeta) = 0 & (nf-fend) = 0 & (nf+fstart-1/tdeta) = 0 & (nf

10、-fend) = 0 & (nf+fstart-1/tdeta) = 0 & (nf-fend) = 0 & (nf+fstart-1/tdeta) = 0);%双门ny2w=xw.*h2w;nsubplot(4,1,4);nplot(linspace(0,N*tdeta,length(h2w),real(ifft(h2w);nstr=sprintf(低截频 %2d Hz, 高截频 %2d Hz 时域：h(t),fstart,fend);ntitle(str);进入夏天，少不了一个热字当头，电扇空调陆续登场，每逢此时，总会进入夏天，少不了一个热字当头，电扇空调陆续登场，每逢此时，总会想起那一把

11、蒲扇。蒲扇，是记忆中的农村，夏季经常用的一件物品。记想起那一把蒲扇。蒲扇，是记忆中的农村，夏季经常用的一件物品。记忆中的故乡，每逢进入夏天，集市上最常见的便是蒲扇、凉席，不论男女老忆中的故乡，每逢进入夏天，集市上最常见的便是蒲扇、凉席，不论男女老少，个个手持一把，忽闪忽闪个不停，嘴里叨叨着少，个个手持一把，忽闪忽闪个不停，嘴里叨叨着“怎么这么热怎么这么热”，于是三，于是三五成群，聚在大树下，或站着，或随即坐在石头上，手持那把扇子，边唠嗑五成群，聚在大树下，或站着，或随即坐在石头上，手持那把扇子，边唠嗑边乘凉。孩子们却在周围跑跑跳跳，热得满头大汗，不时听到边乘凉。孩子们却在周围跑跑跳跳，热得满头

12、大汗，不时听到“强子，别跑强子，别跑了，快来我给你扇扇了，快来我给你扇扇”。孩子们才不听这一套，跑个没完，直到累气喘吁吁，。孩子们才不听这一套，跑个没完，直到累气喘吁吁，这才一跑一踮地围过了，这时母亲总是，好似生气的样子，边扇边训，这才一跑一踮地围过了，这时母亲总是，好似生气的样子，边扇边训，“你你看热的，跑什么？看热的，跑什么？”此时这把蒲扇，是那么凉快，那么的温馨幸福，有母亲此时这把蒲扇，是那么凉快，那么的温馨幸福，有母亲的味道！蒲扇是中国传统工艺品，在我国已有三千年多年的历史。取材的味道！蒲扇是中国传统工艺品，在我国已有三千年多年的历史。取材于棕榈树，制作简单，方便携带，且蒲扇的表面光滑

13、，因而，古人常会在上于棕榈树，制作简单，方便携带，且蒲扇的表面光滑，因而，古人常会在上面作画。古有棕扇、葵扇、蒲扇、蕉扇诸名，实即今日的蒲扇，江浙称之为面作画。古有棕扇、葵扇、蒲扇、蕉扇诸名，实即今日的蒲扇，江浙称之为芭蕉扇。六七十年代，人们最常用的就是这种，似圆非圆，轻巧又便宜的蒲芭蕉扇。六七十年代，人们最常用的就是这种，似圆非圆，轻巧又便宜的蒲扇。蒲扇流传至今，我的记忆中，它跨越了半个世纪，也走过了我们的扇。蒲扇流传至今，我的记忆中，它跨越了半个世纪，也走过了我们的半个人生的轨迹，携带着特有的念想，一年年，一天天，流向长长的时间隧半个人生的轨迹，携带着特有的念想，一年年，一天天，流向长长的时间隧道，袅道，袅