地震子波波形显示.pdfVIP

地震子波波形显示 及一维地震记录合成 一、实验目的 1、 认识地震子波（以雷克子波为例） ，对子波有直观的认识。 2 、 利用线性褶积公式合成一维地震记录。 二、实验内容 1、雷克子波： 2 2 2 f / t 零相位子波源： w t e m cos2 f m t 程序： fm=30; r=3; t=0.002; for n=1:200 w(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*cos(2*pi*fm*t*n); end figure(1),plot(w); 图像（ 1）： 图（1） 2 2 2 f / t 最小相位子波： w t e m sin 2 fm t 程序： fm=10; r=3; t=0.002; for n=1:200 w(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*sin(2*pi*fm*t*n); end figure(1),plot(w); 图像： 图（2 ） 其中 f m 代表子波的中心频率， 代表子波宽度 ,随着 的增大，子波能量后 移，当 =7 时，最小相位子波可视为混合相位子波。 因而若将最小相位子波中 =7 ，则可以间接地得到混合相位子波的图形为： 图（3 ） 从图（ 3 ）中可以看出，经过改变 的值后，整个雷克子波的能量发生了明 显的后移； 从一定程度上可以反映最小能量子波、 混合能量子波、 最大能量子波 的图像上的区别是能量集中区域上的不同。 2 、利用线性褶积公式合成一维地震记录 （1）利用零相位雷克子波作为震源合成一维地震记录： 程序： fm=30; r=3; t=0.002; for n=1:200 w(n)=exp(-(2*pi*fm/r)^2*(t*n)^2)*cos(2*pi*fm*t*n); end for j=1:500 r(j)=0; end r(100)=1.0; r(200)=-0.7; r(300)=0.5; r(400)=0.4; r(500)=0.6; 对地层界面的反射系数进行设定 for j=1:699 f(j)=0; for i=1:500 if(j-i0j-i200)f(j)=f(j)+r(i)*w(j-i); 褶积的主要过程 end end end figure(2),plot(f); 图像： 图（4 ） 在这个一维地震记录图像上可以看出，除了 100 、200 、300 、400 、500 这些存在反射系数的点上发生了变化， 在这些点的周围也发生了些