工程与环境物探实验报告.docVIP

工程与环境物探教程 实验报告 一个反射地震纪录道的形成 专 业：勘查技术与工程 学 号：XXX 姓 名： XXXXXXX 解决的问题 假设在地面以下半空间内有n+1层弹性介质，则有n个波阻抗界面，在第I 层的p波速度为v[i]，密度为p[i]，第i层反射界面上反射系数用R[i]表示，进行实验一个反射地震纪录道的形成。 已知模型如图所示， p1=2000kg/m3 v1=800m/s p2=1600kg/m3 v1=1000m/s p3=2600kg/m3 v1=2000m/s p4=1800kg/m3 v1=900m/s p5=3200kg/m3 v1=2300m/s h1=h2=h3=h3=50m 2．计算要求 1.求每个反射界面的反射系数 波阻抗： Z[i]=D[i]*v[i] — D[i]为第i层密度，v[i]为第i层速度 反射系数： R[i]=(z[i+1]-z[i]/(z[i+1]+z[i] 整理后反射系数表达式如下： R[i]=(V[i+1]*D[i+1]-V[i]*D[i])/(V[i+1]*D[i+1]+V[i]*D[i]) 实验选取选取Ricker（雷克）子波： z[i]=(1-π*f*i*π*f*i)*exp(-(π*f*i)*(π*f*i）); 2.思路方法 主导思想：地震记录道由振幅反射系数序列与子波的褶积构成 多层反射界面，在地面接收的反射波总和的一个地震到记录X（t）可表示为地层反射系数序列R[i]和地震子波z[i]的褶积（卷积）的结果，公式表示如下： X(t)=n=1 4.程序代码 # includestdio.h # includemath.h # includestring.h # define pi 3.1416 int Wave(float,float ); //定义子波 int Reflect(float,float); //定义反射系数 int Convolution(float,float,float ); //定义合成记录 int main() { float d=0.004,f=30.,h=50; / /设定采样间隔 频率 厚度 if(Wave(f,d)!=1) printf(error); if(Reflect(h,d)!=1) printf(error); if(Convolution(f,d,h)!=1) printf(error); return 1; } /***********子波的求取***********/ int Wave(float f,float d) { FILE *fpw; if((fpw=fopen(Wave.text,w))==NULL)printf(error); int i; float Wa[100],t; for(i=0;i50;i++) { t=(50-i)*d; Wa[i]=(float)(1-(pi*pi*f*f*t*t))*(float)exp(-2*pi*pi*f*f*t*t); } for(i=50;i100;i++) { t=(i-50)*d; Wa[i]=(float)(1-(pi*pi*f*f*t*t))*(float)exp(-2*pi*pi*f*f*t*t); } for(i=0;i100;i++) { fprintf(fpw,%f\n,Wa[i]); } fclose(fpw); return 1; } /*******反射序列的求取*******/ int Reflect(float h,float d) { FILE *fpr; fpr=fopen(Rflct.text,w); int i,c=0; float v[5]={800,1000,2000,900,2300}; float p[5]={2000,1600,2600,1800,3200}; float Re[50]; for(i=0;i50;i++) Re[i]=0; for(i=0;i4;i++) { c+=(int)(h/(v[i]*d)); Re[c]=(float)(v[i+1]*p[i+1]-v[i]*p[i])/(v[i+1]*p[i+1]+v[i]*p[i]); } for(i=0;i50;i++) { fprintf(fpr,%f\n,