成都理工大学层析成像实验报告

地球物理层析成像实验报告姓名、XXXXXX学号20100XXXXXX指导老师XXXXX电话XXXXXXXXXX实习一:初识与概述部分及编程实验要求

1.Matlab计算下述函数在自变量采样点集上的值:\-X2-4X-3/2,-3X-1y=2-x+1,-1x1编写脚本程序，分别创建采样向量，函数计算相应的函数值向量，并利用函数绘图―f+4x_3/2,lx3绘制曲线源程序f2plo.tMx,yatlab x-yAyi=-xi2+1;elseclear,clc x=-3:

0.01:3;A⑴八yi=-xi2+4*xi-3/2;N=lengthx;y=zerosl,N;fori=1:N ifendxi=-3xi-1yi=-x2-4*xi

3.程序运行结果end plotx,y-3/2;elseif xi—1xi1PJ-igure1III国If S3IFile EditView InsertTools DesktopWindow Help.按要求完成下述练习4编写函数返回给定向量中最大元素的下标值换言之，返回的满足不小于任何⑴maxi,imax=maxiv将你的学号作为测试数据每一位对应一个向量元素进行测试imax vimaxvfunctionimax=funvimax=l;fori=1:119ifvl i+lt t=vl,i+l;imax=i+l;%else%imax1=0;endend imaxend实习

二、射线追踪实验内容

1.模型为一三层模型，范围该模型含有三个水平层，上下两个层速度为中间为一高速层，速度为高速层厚度范围是激发点〜lOOmXIOOm有个,坐标依次为是⑸、叫、，接收点有个，坐3000m/s,4500m/s,5060m标分别为、将此模型网格化为19El0E20E30,15,…E100,9521个单元格.根据法追踪出条射线的路径和旅行时Rl100,0R2100,5^3100,10,…,R21100,10010X10实验源程序

2.LTI19X21口,〃定义道数#includestdio.h#includemath.h double convdouble Wa[],double Re〃定义角度个数double Con[],int n_f,int m_f;#defineA5〃定义子波长度#define B900〃定义采样点个数#define C61〃定义常量#define D70#define pi

3.14159265pi voidmainQ〃定义第一层厚度,第一层速度，检波点间隔，采样间隔,子波频率〃//double h=

80.0v=

2000.0dx=

15.0,dt=

0.002,f=

30.0;仇/〃定义初始角度double2=

0.05*

180.0;〃各种数组的定义〃double lsjt,w;〃定义反讲系数序列・・・〃定义子波数组,反射系数数组,褶int re[A];积结果数组double Wa[C+D lLRe[C+D l]£on[C+D l];〃定义角度数组,接收点数组，褶积数组〃定义地震记录数组，炮点数组，检波点数组double angle[B]jeceive[BLt[A];doubledata[A][D],S[lLR[A];FILE子波；*fp,*fpL*fp2,*fp3,*fp4;二反射系数；fp=fopen.txtJw+”接收点，；fpl fopen,segy”Jwb+”二地震记录；fp2=fopen.txt”w+””接受时间，,r；fp3fopen.segy Jwb+”fp4=fopen.txt”w+”for[int i=0;iC;i++叭，{w=pi*pi*f*f*[i-30*i-30*dt*dt;Wa[i]=[l-

2.0*w*exp[-

1.0]*w];fprintffpj%n Wa[i];}fclosefp;S

[0]=

0.0;h=

80.0;fori=0;iA;i++〃开始射线追踪{R[i]=

95.0+i*dx;}fori=0;iB;i++{angle[i]=a*i+l;receive[i]=2,0*h*tanangle[i];H,l,printf%lf\n receive[i];forint j=O;jA;j++{iffabsreceive[i]-R[j]]0,1]{fprintffp2receive[i];t[j]=receive[i]/sin angle[i]/v;re[j]=int[t[j]/dt;,,，,printf%d%d\n J+l re[j];fprintffp4/%lf%d\n t[jLre[j];}}}fclosefp2;fclosefp4;fori=0;iA;i++{forint j=O;j D;j++]{ifQ==re[i]Re[j]=

0.2;else Re[j]=O.O;；lt=floatRe[j];〃褶积函数的调用〃/lfwritelt sizeoffloatl/fpl}convWa,Re£onCD;forj=0;jD;j++{data[i][j]=Con[j+30];ls=floatdata[i][j]J;l,/fwritels sizeoffloatlfp3;}}fclosefpl;〃定义褶积函数fclosefp3;},doubleconvdoubleWa[]double Re[],double Con[]jnt n_fjnt m_f{int ij;int N=n_f+m_f-l;double sum;fori=n_f;iN;i++Re[i]=

0.0;fori=m_f;iN;i++Wa[i]=

0.0;fori=0;iN;i++sum=

0.0;forQ=0;j=i;j++sum+=Re[j]*Wa[i-j];Con[i]=sum;return Con[i];运行结果

3.地震记录:

1120.

00.

000.10-■o.101反射系数:。