MATLAB-SWPUMATLAB-SWPU第七章第七章- -地震地震处理相关知识处理相关知识7.1 7.1 地震数据格式地震数据格式常见的地震数据体有常见的地震数据体有SEGSEG格式和格式和NOSEG NOSEG 格式，格式， SEG SEG是（是（ The The Society of Exploration GeophysicistsSociety of Exploration Geophysicists）的缩写，）的缩写， 而而SEGSEG格式格式有多种，其中最常用的就是有多种，其中最常用的就是SEGYSEGY格式地震数据体文件。它有格式地震数据体文件。它有32003200个字节的图头（也叫卷头），个字节的图头（也叫卷头），400400字节的文件头，字节的文件头，240240字节的道头，字节的道头，道数据道数据，.道头，道数据，很多个。图头信息我们道头，道数据，很多个。图头信息我们一般不关心，不太常用。文件头里有些信息可以用到，例如数据一般不关心，不太常用。文件头里有些信息可以用到，例如数据格式了，样品数了，采样间隔了。而道头在实际加载的过程中最格式了，样品数了，采样间隔了。而道头在实际加载的过程中最常用，其中最常用的是常用，其中最常用的是X X、Y Y、道号和线号。但在实际的数据中它、道号和线号。但在实际的数据中它们存放的位置不一定相同，要具体数据体具体分析，可以用一个们存放的位置不一定相同，要具体数据体具体分析，可以用一个SEGYSEGY格式分析器进行转换格式分析器进行转换 7.1 7.1 地震数据格式地震数据格式SEGY数据由文件头和数据体组成数据由文件头和数据体组成文件头总长度为文件头总长度为3600字节，分两部分。字节，分两部分。文件头第一部分文件头第一部分长度：长度：3200bytes；组成：；组成：80bytes*40；特性：；特性：EBCDIC字符集，字符集，参数卡，需要转换为参数卡，需要转换为ASCII码后才能显示。码后才能显示。文件头第二部分文件头第二部分长度：长度：400bytes；数据类型：；数据类型：32位、位、16位的整型；特性：二进制头，位的整型；特性：二进制头，记录数据体信息。记录数据体信息。数据体由多个数据道组成。每道数据分两部分：道头、采样数数据体由多个数据道组成。每道数据分两部分：道头、采样数据，数据类型：据，数据类型：32位的浮点型位的浮点型，占，占4个字节。个字节。7.1 7.1 地震数据格式地震数据格式工作站SEGY数据存储格式有两种：IEEE和IBMIEEE和IBM的整型存储与微机格式的存储不同之处，IEEE和IBM的高字节在前、低字节在后，即BigEndian，微机则是低字节在前、高字节在后，即LittleEndian。文件头格式文件头格式道头格式道头格式例读取地震数据并绘图nclearallnindata=altreadsegy(mianbomodle.sgy);%读取数据nTr=length(indata(1,:);%道数nT=length(indata(:,1);%采样点数ns_wplot(indata);%绘剖面例提取第10道单独显示，并且加入随机噪声nclearallnindata=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道中加入随机噪声nclearallnindata=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常用处理-滤波nclearallnindata=altreadsegy(mianbomodle.sgy);nTr=length(indata(1,:);nT=length(indata(:,1);ndata=indata(:,10);%提取第10道来显示ntemp=randn(T,1);%生成随机数据ndata=data+temp;nfdata=fft(data);nplot(abs(fdata);例低通滤波器实现nnfdata=fft(data);nfT=length(fdata);nforj=100:fT%低通滤波器nfdata(j)=0;nendnplot(abs(fdata);例低通滤波器实现-反变换nnfdata=fft(data);nfT=length(fdata);nforj=100:fT%低通滤波器nfdata(j)=0;nendnplot(abs(fdata);ndata2=ifft(fdata);%逆变换ns_wplot(data);ns_wplot(data2);例高通滤波器实现nnfdata=fft(data);nfT=length(fdata);nforj=1:200%高通滤波器nfdata(j)=0;nendnplot(abs(fdata);例带通滤波器实现nnfdata=fft(data);nfT=length(fdata);nforj=1:fT%带通滤波器nif(j300|jltime(i)/tspace-twindow&jltime(i)/tspace-twindow&j=0&(nf-fend)=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-fend)=0&(nf+fstart-1/tdeta)=0&(nf-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(低截频%2dHz,高截频%2dHz时域：h(t),fstart,fend);ntitle(str);进入夏天，少不了一个热字当头，电扇空调陆续登场，每逢此时，总会想起进入夏天，少不了一个热字当头，电扇空调陆续登场，每逢此时，总会想起那一把蒲扇。蒲扇，是记忆中的农村，夏季经常用的一件物品。记忆中的故那一把蒲扇。蒲扇，是记忆中的农村，夏季经常用的一件物品。记忆中的故乡，每逢进入夏天，集市上最常见的便是蒲扇、凉席，不论男女老少，个个手持乡，每逢进入夏天，集市上最常见的便是蒲扇、凉席，不论男女老少，个个手持一把，忽闪忽闪个不停，嘴里叨叨着一把，忽闪忽闪个不停，嘴里叨叨着“怎么这么热怎么这么热”，于是三五成群，聚在大树，于是三五成群，聚在大树下，或站着，或随即坐在石头上，手持那把扇子，边唠嗑边乘凉。孩子们却在周下，或站着，或随即坐在石头上，手持那把扇子，边唠嗑边乘凉。孩子们却在周围跑跑跳跳，热得满头大汗，不时听到围跑跑跳跳，热得满头大汗，不时听到“强子，别跑了，快来我给你扇扇强子，别跑了，快来我给你扇扇”。孩。孩子们才不听这一套，跑个没完，直到累气喘吁吁，这才一跑一踮地围过了，这时子们才不听这一套，跑个没完，直到累气喘吁吁，这才一跑一踮地围过了，这时母亲总是，好似生气的样子，边扇边训，母亲总是，好似生气的样子，边扇边训，“你看热的，跑什么？你看热的，跑什么？”此时这把蒲扇，此时这把蒲扇，是那么凉快，那么的温馨幸福，有母亲的味道！蒲扇是中国传统工艺品，在是那么凉快，那么的温馨幸福，有母亲的味道！蒲扇是中国传统工艺品，在我国已有三千年多年的历史。取材于棕榈树，制作简单，方便携带，且蒲扇的表我国已有三千年多年的历史。取材于棕榈树，制作简单，方便携带，且蒲扇的表面光滑，因而，古人常会在上面作画。古有棕扇、葵扇、蒲扇、蕉扇诸名，实即面光滑，因而，古人常会在上面作画。古有棕扇、葵扇、蒲扇、蕉扇诸名，实即今日的蒲扇，江浙称之为芭蕉扇。六七十年代，人们最常用的就是这种，似圆非今日的蒲扇，江浙称之为芭蕉扇。六七十年代，人们最常用的就是这种，似圆非圆，轻巧又便宜的蒲扇。蒲扇流传至今，我的记忆中，它跨越了半个世纪，圆，轻巧又便宜的蒲扇。蒲扇流传至今，我的记忆中，它跨越了半个世纪，也走过了我们的半个人生的轨迹，携带着特有的念想，一年年，一天天，流向长也走过了我们的半个人生的轨迹，携带着特有的念想，一年年，一天天，流向长长的时间隧道，袅长的时间隧道，袅结束