反射地震记录道的形成.doc

《工程与环境物探》课程实验报告 一个反射地震记录道的形成 姓 名： xxx 学 号：xxxx 解决的问题： 假设在地面以下半无限空间内有层弹性介质，则共有个波阻抗界面，第i层的p波速度为v[i],密度为p[i],第i层反射界面上的反射系数用R[i]表示 已知模型如图所示。 地面 p1 v1 h1 R1 p2 v2 h2 R2 p3 v3 h3 R3 p4 v4 h4无限延伸 图 水平层状介质模型示意图 2．R[i]=(z[i+1]-z[i])/( z[i+1]+z[i] 2.选取适当的子波，本论文选取Ricker子波 雷克子波表达式： 3.解决思路，方法： 主导思想：振幅反射系数序列与子波的褶积构成一个地震记录道 如果地下实际介质存在N个反射界面的话，地面可以接收到每一个界面上的反射波，于是一个实际地震道上记录了N个反射波。根据实际钻井和测井资料可知，地下介质近似为层状介质，并存在大量的薄层。这表明每个层界面都可能产生反射，成为一反射面。由此可见，一个实际地震记录道就是由这些无数多个反射子波（地震子波）组成的复合振动。在地面上接收到N个反射波总和的一个地震道记录X(t)可表示为一个反射记录道是地层反射系数序列和地震子波的褶积（卷积）结果。于是反射道可表示为： 4.程序代码（模块）及说明 #include stdio.h #include math.h #define pai 3.1415926 #define n 51 //n表示子波长度// #define m 128 //m表示反射序列长度// void main() { double v[]={2000.0,2500.0,3000.0,3500.0}; //速度参数// double p[]={2.2,2.4,2.6,2.8}; //密度参数// double h[]={100.,125.,150.}; //层厚度参数// double wv[n]={0.0}; //子波数据赋初值// double Rf[m]={0.0}; //反射系数// double r[3]={0.}; //波阻抗// double con[n+m-1]={0.}; //合成记录// double z[4]={0.}; //波阻抗// double f=40.0; //频率// double dt=0.002; //采样间隔// double t; int i,j; FILE *fp; fp=fopen(反射道.txt,w); //==========子波求取==========// for(i=0;in/2;i++) { t=i*dt; wv[25+i]=(1-2*(pai*f*t)*(pai*f*t))*exp(-(pai*f*t)*(pai*f*t)); wv[25-i]=wv[25+i]; } for(i=0;in;i++) { //printf(%9.8f\n,wv[i]); } //=========求反射系数=======// for(i=0;i3;i++) { z[i]=p[i]*v[i]; z[i+1]=p[i+1]*v[i+1]; r[i]=(z[i+1]-z[i])/(z[i+1]+z[i]); //printf(%f\t,r[i]); } for(i=0;im;i++) { if(fabs(i*dt-h[0]/v[0])1e-4) {Rf[i]=r[0];} else if(fabs(i*dt-(h[0]/v[0]+h[1]/v[1]))1e-4) {Rf[i]=r[1];} else if(fabs(i*dt-(h[0]/v[0]+h[1]/v[1]+h[2]/v[2]))1e-4) {Rf[i]=r[2];} else Rf[i]=0.0; } for(i=0;im;i++) { printf(%f\t,Rf[i]); } //======