地震子波波形显示及一维地震合成记录精要.docxVIP

地震勘探原理实验一 地震子波波形显示及一维地震合成记录 姓名： 学号：专业：地球物理勘察技术与工程 级 一、实验目的 认识子波，对子波的波形有直观的认识。（名词：零相位子波，混合相位子波，最小相位子波；了解子波的分辨率与频宽的关系；） 利用褶积公式合成一维地震记录。 二、实验步骤 雷克子波 零相位子波 （最小相位子波） n= m1/m2为最大波峰m1和最大波谷m2之比 钟型子波 xw为初相m为时间域主波峰与次波峰之比 w(t)=exp(-2*Fm^2*t^2*ln(n))*sin(T-2*pi*Fm*t) n=m1/m2 最大相位子波 （最大相位子波请同学们自己查找相关文献完成，非必须完成） 其中代表子波的中心频率， t =i*dt，dt为时间采样间隔，i为时间离散点序号; 这里可以为 = 10，25，40，100 Hz等，采样间隔dt=0.002秒，i为0～256; 根据公式编程实现不同频率的零相位子波的波形显示； 不同中心频率的零相位子波图 = 25: = 100: 其地质模型为： 震源 检波器 第一层介质r(100)=1.0 第二层介质r(200)=-0.7 第三层介质r(300)=0.5 第四层介质r(400)=0.4 第五层介质r(450)=0.6 设计反射系数(n=512)，n为地层深度，其中，为第一层介质深度；，为第二层介质深度；，为第三层介质深度；，为第四层介质深度；，为第五层介质深度；其它为0。 地震波在介质中传播，当到达介质分界面时，发生反射和透射，反射波被检波器接受，生成地震记录。反射系数表示地震波在两层介质分界面的能量重新分配，如r(100)=1.0，表示地震波入射到分界面时，只有一种波，反射纵波（或反射横波）。反射系数不为1.0时，表示当地震波入射到分界面时，产生两种反射波。反射系数为正，表示反射波相位与入射波相位相差2π；反射系数为负，表示反射波相位与入射波相位相差π，存在半波损失。 应用褶积公式合成一维地震记录，并图形显示； 应用褶积公式求f（n）的程序为： #includestdio.h #includemath.h #define PI 3.1415926 #define FM 100 void main() { double fac(double x[],double y[],double z[],int m,int n); FILE *fp; int i,j,x; double W,dt=0.002,t,a[256]; double b[512]={0}; double r[512]={0}; r[100]=1.0; r[200]=-0.7; r[300]=0.5; r[400]=0.4; r[450]=0.6; fp=fopen(Date.txt,w+); printf(please input x:\n); scanf(%d,x); for(i=0;i256;i++) { t=i*dt; if(x==1) W=exp(-2*FM*FM*t*t*log(2))*sin(2*PI*FM*t); else if (x==2) W=(1-2*pow(PI*FM*t,2))*exp(-pow(PI*FM*t,2)); else if (x==3) W=exp(-FM*FM*t*t*log(2))*cos(2*PI*FM*t+PI/4); a[i]=W; } fac(r,a,b,512,256); for(j=0;j512;k++) { fprintf(fp,%f\n,b[j]); } } double fac(double x[],double y[],double z[],int m,int n) { int i,j; for(i=0;i=m+n-1;i++) { double sum=0.0; for(j=0;j=m;j++) { if(i-j0i-j=256) sum+=x[j]*y[i-j]; } z[i]=sum; } } 三、实验结果 一维反射系数序列图形显示为： 零相位子波与反射系数褶积后的地震记录图形显示： = 25: = 100: 最小相位子波与反射系数褶积后的地震记录图形显示： =