地震反射道

工程物探VB课程报告学年学期：2016~2017学年第1学期院别：土木工程学院专业：勘查技术与工程专业方向：岩土方向班级：勘查1401学生：任霄瑶学号：3140604120指导教师：江凡一、课设目的认识学习地震反射记录道二、课设内容1.解决思路：当地下介质近似为层状介质时，并存在大量的薄层，这表明每个层界面都可以产生反射并成为一个反射面。由于在层界面的两侧速度和密度存在一定的差异，地震波到达层界面之后会出现反射、折射等现象，因此可以计算得到不为零的反射系数，由此可见，一个实际地震记录道就是由这些无数多个子波组成的复合震动。显然，振动的幅值大小与界面的反射系数成正比。利用地震子波与反射系数做卷积，就得到了单道地震记录道。在地面上接收到N个反射波总和的一个地震道记录X(t)可表示为一个反射记录道是地层反射系数序列tR和地震子波tb的褶积（卷积）结果。于是反射道可表示为：2.求每个反射界面的反射系数，即求得反射系数序列波阻抗z[i]=p[i]*v[i]反射系数R[i]=(z[i+1]-z[i])/(z[i+1]+z[i]3.选取适当的子波.4.程序说明：Dimlayer(100)AsDoubleDimntAsIntegerDimv1AsDoubleDimv2AsDoubleDimrou1AsDoubleDimrou2AsDoubleDimflagAsIntegerDimtralen(2000)AsDoubleDimw(2000)AsDoubleDimfAsDoublettntnNnbRbRt1(X）DimdtAsDoubleDimltwAsIntegerDimrefAsDoubleConstPI=3.1415926535PrivateSubForm_Load()flag=0Picture1.AutoRedraw=TruePicture1.Scale(0,0)-(99,99)DimoridepAsIntegeroridep=50Fori=0To99layer(i)=oridepNextiFori=0To99Picture1.Line(i,layer(i))-(i+1,layer(i+1))Nextiv1=1000#v2=2000#rou1=1.1rou2=1.2ref=(v2*rou2-v1*rou1)/(v2*rou2+v1*rou1)f=100dt=0.001ltw=Int((4#/f/dt))+1Callricker(w(),f,dt,ltw)rickerPicture3.AutoRedraw=TruePicture3.Scale(-20,0)-(20,40)Fori=0Toltw-1Picture3.Line(w(i),i)-(w(i+1),i+1),vbRedNextitralen(Int(oridep/v1*1000)+1)=refnt=120CallConvolution(w,ltw,tralen,nt)Picture2.AutoRedraw=TruePicture2.Scale(-15,0)-(1000,120)Fori=0To99Forj=0Tont-1Picture2.Line(tralen(j)+i*10,j)-(tralen(j+1)+i*10,j+1),vbBlueNextjNextiEndSubPrivateSubPicture1_MouseDown(ButtonAsInteger,ShiftAsInteger,xAsSingle,yAsSingle)flag=1DimpointXAsIntegerDimpointYAsIntegerpointX=Int(x)+1pointY=Int(y)layer(pointX)=pointYPicture1.ClsFori=0To98Picture1.Line(i,layer(i))-(i+1,layer(i+1))NextiEndSubPrivateSubPicture1_MouseMove(ButtonAsInteger,ShiftAsInteger,xAsSingle,yAsSingle)If(flag=1)ThenDimpointXAsIntegerDimpointYAsIntegerpointX=Int(x)+1pointY=Int(y)layer(pointX)=pointYPicture1.ClsFori=0To98Picture1.Line(i,layer(i))-(i+1,layer(i+1))NextiEndIfEndSubPrivateSubPicture1_MouseUp(ButtonAsInteger,ShiftAsInteger,xAsSingle,yAsSingle)flag=0Picture2.ClsFori=0To99Forj=0Tonttralen(j)=0Nextjtralen(Int(layer(i)/v1*1000)+1)=refCallConvolution(w,ltw,tralen,nt)Forj=0Tont-1Picture2.Line(tralen(j)+i*10,j)-(tralen(j+1)+i*10,j+1),vbBlueNextjNextiEndSubFunctionConvolution(x()AsDouble,lenXAsInteger,y()AsDouble,lenYAsInteger)Dimx2(5000)AsDouble'实际用到长度lenN=lenX+2*lenYDimy2(5000)AsDouble'实际用到长度lenYDimz(5000)AsDouble'实际用到长度lenX+lenY-1Fori=0TolenY-1y2(i)=y(lenY-i-1)NextiFori=0TolenX-1x2(i+lenY)=x(i)NextiFori=0TolenX+lenY-2Forj=0TolenY-1z(i)=z(i)+x2(i+j+1)*y2(j)NextjNextiDimkAsIntegerk=Int(lenX/2)Fori=0TolenY-1y(i)=z(i+k)NextiEndFunctionFunctionricker(w()AsDouble,fAsDouble,dtAsDouble,ltwAsInteger)Fori=0Toltw-1tmp=(i-ltw/2)*dttmp=-PI*PI*f*f*tmp*tmpw(i)=20#*Exp(tmp)*(1+2*tmp)NextiEndFunction3.图形解释：①地层深度Z=50，V1=x=1000，ρ1=a=1.1，V2=y=2000，ρ2=b=1.2②地层深度Z=50，V1=x=1000，ρ1=a=1.1，V2=y=2000，ρ2=b=1.2③地层深度Z=50，V1=x=800，ρ1=a=1.1，V2=y=1500，ρ2=b=1.2④地层深度Z=59，V1=x=1000，ρ1=a=1.2，V2=y=2000，ρ2=b=1.4三、课设总结地震子波由震源激发，在地层中传播，因为在沉积地层中，每层介质的物理性质不相同，从而使得地震波的传播速度也不相同。当地震波传播到两层介质的分界面时，会发生反射，由于每层介质的反射系数不同，所以反射波的能量也不相同，检波器接收到不同时刻的、不同能量的反射波，形成一个地震记录。