MATLABSWPU第七章地震处理相关知识.pptVIP

第七章 地震处理相关知识 7.1 地震数据格式 7.2 常用处理 7.3 地震属性提取 7.1 地震数据格式 7.1 地震数据格式 7.1 地震数据格式 文件头格式 道头格式 例 读取地震数据并绘图 clear all indata=altreadsegy(‘mianbomodle.sgy’);%读取数据 Tr=length(indata(1,:));%道数 T=length(indata(:,1)); %采样点数 s_wplot(indata);%绘剖面 例 提取第10道单独显示，并且加入随机噪声 clear all indata=altreadsegy(mianbomodle.sgy); Tr=length(indata(1,:)); T=length(indata(:,1)); data=indata(:,10);%提取第10道来显示 temp=randn(T,1);%生成随机数据 data=data+temp; s_wplot(data); %s_cplot(data);%变面积显示 例 在第10道中加入随机噪声 clear all indata=altreadsegy(mianbomodle.sgy); Tr=length(indata(1,:)); T=length(indata(:,1)); temp=randn(T,1);%加入随机噪声 indata(:,10)=indata(:,10)+temp; s_wplot(indata); %s_cplot(data);%变面积显示 7.2 常用处理-滤波 clear all indata=altreadsegy(mianbomodle.sgy); Tr=length(indata(1,:)); T=length(indata(:,1)); data=indata(:,10);%提取第10道来显示 temp=randn(T,1);%生成随机数据 data=data+temp; fdata=fft(data); plot(abs(fdata)); 例 低通滤波器实现 ……… fdata=fft(data); fT=length(fdata); for j=100:fT %低通滤波器fdata(j)=0; end plot(abs(fdata)); 例 低通滤波器实现-反变换 ……… fdata=fft(data); fT=length(fdata); for j=100:fT %低通滤波器fdata(j)=0; end plot(abs(fdata)); data2=ifft(fdata);%逆变换 s_wplot(data); s_wplot(data2); 例 高通滤波器实现 ……… fdata=fft(data); fT=length(fdata); for j=1:200 %高通滤波器fdata(j)=0; end plot(abs(fdata)); 例 带通滤波器实现 ……… fdata=fft(data); fT=length(fdata); for j=1:fT %带通滤波器if(j300||j40)fdata(j)=0;end end plot(abs(fdata)); 常用处理-子波提取 clear all read_segy_file(‘mig.sgy’)； seismic=s_data; sdata=s_seismic2wavelet(seismic,{color,blue},{type,zero-phase},{wlength,80}); s_wplot(sdata); 常用处理-子波提取 clear all read_segy_file(mig.sgy); seismic=s_data; bt=s_seismic2wavelet(seismic,{color,blue},{type,zero-phase},{wlength,80}); kesai=zeros(1,200); kesai(12)=-0.2;kesai(33)=0.4; kesai(37)=-0.15; kesai(52)=-0.5; kesai(72)=0.35;kesai(87)=-0.1; plot(bt.traces); st=conv(bt.traces,kesai); s_wplot(st); 常用处理-希尔伯特变换 clear all seismic=read_segy_file(mig.sgy); sdata=s_hilbert(seismic); s_cplot(sdata); s_cplot(seismic); 常用处理-能量谱 clear all seismic=read_segy_file