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

《成都理工大学层析成像实验报告.docx》由会员分享，可在线阅读，更多相关《成都理工大学层析成像实验报告.docx（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、地球物理层析成像实验报告姓 名：XXXXXX、学 号： 20100XXXXXX指导老师：XXXXX电 话：XXXXXXXXXX实习一:Matlab初识与概述部分及编程1 .实验要求：计算下述函数在自变量采样点集上的值: -x2-4x-3)/2,-3x-1y = -x2+1,-1 x 1(f+4x_3)/2,lx=-3 & x(i) 1 & x(i) t t=v(l,i+l);imax=i+l;%else% imax 1=0;endend imax end实习二、射线追踪1 .实验内容：模型为一三层模型，范围lOOmXIOOm。该模型含有三个水平层，上下两个 层速度为3000m/s,中间为一高速

2、层，速度为4500m/s,高速层厚度范围是 5060m。激发点有19个,坐标依次为是El(0、E2(0叫、E3(0,15),，E10(0,95)接收点有 21 个，坐标分别为 Rl(100, 0)、R2 (100,5)3(100,10), R21(100,100)。将此模型网格化为10X10个单元格.根据LTI法追踪出19X21 条射线的路径和旅行时。2 .实验源程序：#include #include double conv(double Wa,double Re口,double Con,int n_f,int m_f); #defineA5定义道数# define B 900定义角度个数#

3、 define C 61定义子波长度# define D 70定义采样点个数#define pi 3.14159265定义 pi 常量 void mainQ(double h=80.0/v=2000.0/dx=15.0,dt=0.002, f=30.0;定义第一层厚度, 第一层速度，检波点间隔，采样间隔,子波频率double 2=0.05*仇/180.0;定义初始角度double lsjt,w;各种数组的定义int reA; 定义反讲系数序列double WaC+DlLReC+DlonC+Dl;定义子波数组,反射系数数 组,褶积结果数组double angleBjeceiveBLtA;定义角度

4、数组,接收点数组，褶积数组 double dataAD,SlLRA;定义地震记录数组，炮点数组，检波点数组 FILE *fp,*fpL*fp2,*fp3,*fp4;fp = fopen(子波.txtJw+”)；fpl二fopen(反射系数,segy” Jwb+”)；fp2=fopen(接收点.txt，”w+”)；fp3二fopen( 地震记录.segy,rJwb+”)；fp4=fopen(”接受时间.txt，”w+”)；forint i=0;iC;i+) w=pi*pi*f*f*i-30)*(i-30)*dt*dt;Wai=l-2.0*w)*exp(-1.0*w; fprintf(fpj% 叭

5、n，Wai); fclose(fp);S0=0.0;h=80.0;for(i=0;iA;i+)Ri=95.0+i*dx;开始射线追踪for(i=0;iB;i+) anglei=a*(i+l);receivei=2,0*h*tan(anglei);printf(H%lfn,l,receivei);for(int j=O;jA;j+)if(fabs(receivei-Rj0,1fprintf(fp2receive i);tj=receive i/sin (angle i)/v;rej=(int)tj/dt);printf(%d %dn,J+l,rej);fprintf(fp4/%lf %dn，tj

6、Lrej);fclose(fp2);fclose(fp4);for(i=0;iA;i+)for(int j=O;j D;j+ifQ=rei)Rej=0.2;else Rej=O.O;lt=float(Rej);fwrite(</sizeof(float)ll/fpl)；褶积函数的调用conv(Wa,ReonCD);for(j=0;jD;j+) datai j=Conj+30;ls=float(datai jJ;fwrite(&lslsizeof(float),l/fp3);fclose(fpl);fclose(fp3);定义褶积函数double conv(double Wa,double Re,double Conjnt n_fjnt m_f) int ij;int N=n_f+m_f-l;double sum;for(i=n_f;iN;i+)Rei=0.0;for(i=m_f;iN;i+)Wai=0.0;for(i=0;iN;i+)sum=0.0;forQ=0;j=i;j+) sum+=Rej*Wai-j; Coni=sum;)return Coni; )3 .运行结果(1) 地震记录:反射系数:(2)1 0.001.00o. 100. 10-