地震层析成像的正演实验精要.docVIP

1-1-1正演的速度模型图 图1-1-2分块均匀的模型  1-2正演的后的走时图  图1-3 反演前后速度对比图  图1-5-a第0炮，第5接收点的数据  图1-5-b-1正演第1炮，第8接收点（0为开始的激发点，0开始的接收点） 图1-5-b-2 与1-5-b-2对应的验证图形（附注：由于本人u盘被病毒入侵，导致本人做得CAD图丢失，此图引用廖松杰同学的CAD图像，但是1-5-b-2由本人程序自己得出，特此说明。） 图 1-5-c四边放炮，四边接收左方第2激发， 图1-5-b单边接收第0炮， 第25接收的r图像 第5接收点的r数据图  正反演的程序 单边放炮单边接收： #includestdio.h #includestdlib.h #includemath.h void fun1(int n,double R[144][108],double t[12][12]); void fun2(double k,double o,double t[12][12],double R[144][108],int m,int n);//k为斜率，o为炮点坐标，相当于截距；j； double fun(double x1,double y1,double x2,double y2); void main() { FILE *fp; int i,j,m,n,l,f; double c[12],d[12],K[12][12],r[12][9]={0},v[12][9]={0.0},t[12][12]={0.0},u,o,w,R[144][108]={0.0},k,v2[144]={0}; float v1; //*******************************************************************************************************// for(i=0;i12;i++)//第i行 { for(j=0;j9;j++)//第j列 { v[i][j]=3000; } } v[2][2]=5000.0; v[3][2]=5000.0; v[8][5]=2000.0; v[8][6]=2000.0; for(i=0;i12;i++)//第i行 { for(j=0;j9;j++)//第j列 { v2[j*12+i]=v[i][j]; } } fp=fopen(速度,wb); for(i=0;i12;i++) { for(j=0;j9;j++) { v1=v[i][j]; fwrite(v1,sizeof(float),1,fp); } }fclose(fp); //*******************************************************************************************************// //************************************* 计算各点的斜率***************************************************// for(i=0;i12;i++) { //printf(第%d炮的斜率\n,i); c[i]=(i+0.5)*5.0;//***************************************************************激发点********// for(j=0;j12;j++) { d[j]=(j+0.5)*5;//***********************************************************接收点********// K[i][j]=(d[j]-c[i])/(9.0*3);//**********************************************K斜率*********// //printf(K[%d][%d]=%f\n,i,j,K[i][j]); }//printf(\n); } //**********************************************************************************************************//