matlabswpu 地震处理相关知识.pptxVIP

第七章 地震处理相关知识;7.1 地震数据格式;7.1 地震数据格式;7.1 地震数据格式;;;;;;;例 读取地震数据并绘图;;例 提取第10道单独显示，并且加入随机噪声;;例 在第10道中加入随机噪声;7.2 常用处理-滤波;;例 低通滤波器实现;;例 低通滤波器实现-反变换;;例 高通滤波器实现;;;例 带通滤波器实现;;;常用处理-子波提取;;常用处理-子波提取;常用处理-希尔伯特变换;;常用处理-能量谱;;7.3 地震属性提取-能量;;属性提取-FFT分析;属性提取-FFT分析;;7.4 生成理论地震模型;nRandRange=0.2; ntrand=-nRandRange+rand(1,length(st))*nRandRange*2;%随机数组 str=sprintf('随机信号 n(t) 范围 %.1g到%.1g',-nRandRange,nRandRange); %xt=st+ntrand; xt=st; nfft=2^ceil(log(length(xt))/log(2));%产生大于xt的最接近2的n次方的整数 xw=fft(xt,nfft); plot(nt,bt); title('地震子波 b(t)'); axis([0 0.12 -1 1]);; subplot(2,2,1); plot(ntrand); title(str); subplot(2,2,2); stem(kesai); title('反射系数序列 \xi(t)'); subplot(2,2,3); plot(linspace(0,N*tdeta,length(xt)),xt); title('x(t)'); axis([0 0.12 -1 1]); subplot(2,2,4); plot(abs(xw)); title('|X(W)|');;7.5 相关计算;subplot(2,1,1); plot(ntX,xt); title('x(n)'); subplot(2,1,2); plot(ntY,yt); title('y(n)'); figure(201); subplot(3,1,1); plot(linspace(0,MX*tdetaX,length(rxy)),rxy); title('\bf\fontsize{12}{\gamma}_{xy}(\tau)'); subplot(3,1,2); plot(linspace(0,MX*tdetaX,length(ryx)),ryx); title('\bf\fontsize{12}{\gamma}_{yx}(\tau)'); subplot(3,1,3); plot(linspace(0,MX*tdetaX,length(ryx)),r_xy); title('\bf\fontsize{12}{\gamma}_{xy}(-\tau)');;7.6 复杂的滤波程序;计算时间域的信号;计算时间域的信号;频率域滤波后结果;频率域滤波后结果;设计的滤波器频域和时间域结果;设计的滤波器频域和时间域结果;时域褶积滤波与原始信号在频域进行比较;不同截频滤波效果比较;不同截频滤波效果比较;不同截频滤波效果比较;对不同频率的H(W)求傅里叶反变换;对不同频率的H(W)求傅里叶反变换;对不同频率的H(W)求傅里叶反变换