《地球物理层析成像实验报告》.doc

本科生实验报告 实验课程地球物理层析成像学院名称地球物理学院专业名称勘查技术与工程学生姓名学生学号指导教师曹俊兴实验地点5417实验成绩二〇一五 年 三 月二〇一五 年 四 月 学生实验 心得 在学习了地球物理层析成像之后，收获了很多专业知识，比如学会了利用层析成像的手段反演出地下地质体的异常，同时也学会了利用我们的专业知识解决不同的地质问题。程序语言作为一种工具一方面起到了辅助作用，另一方面我们也学会了一种思维方式，如何设计程序，如何用程序解决我们的复杂问题。在今后的学习工作当中，进一步拓宽思路，勇于创新，能够获得更多的知识。 学生（签名）：2015年 4 月 28 日 指导 教师 评语成绩评定： 指导教师（签名）：年 月 日地震走时层析成像实验 ——《地球物理正反演概论》课程结业报告 学号：0423 姓名：马力衡 专业：勘查技术与工程 手机： 摘要 运用C语言程序，正演得到地震走时和射线在传播过程中经过离散化处理单元格内的距离。通过反演程序反演出地下异常速度值将反演所得速度值成图与原始速度成图进行比较得出结论#include stdio.h #include math.h void main() { int v[12][9]; int m,n，i，j; FILE *fp0; fp0=fopen(速度.txt,r); for(i=0;i12;i++) { for(j=0;j9;j++)fscanf(fp0,%d,v[i][j]); } double b[12];//截距; double xl[12][12]; //斜率; double y_jf[12],y_js[12];//激发点与接收点的纵坐标 for(i=0;i12;i++) { y_jf[i]=1.5+3.0*i;//激发点点坐标的方程 for(j=0;j12;j++) {y_js[j]=1.5+3.0*j;//接收点坐标的方程xl[i][j]=(y_js[j]-y_jf[i])/(45-0);//斜率printf(%f\n,xl[i][j]); } } for(i=0;i12;i++) { b[i]=1.5+i*3;//每条射线截距 } //以上在求射线的斜率和射线在纵轴上的截距 // double ft_t=0.0;//每一格的时间; double fl[12][12][12][9]; //每一格射线的长度; double Time[12][12];//每条射线的时间; double X0,Y0;//第一个点坐标; double X1,Y1;//第二个点坐标; double x_0,x_1,y_0,y_1; //判定的 x，y; double x0,x1,y0,y1;//小格的边界; FILE *fp_ds; fp_ds=fopen(每一小格的距离.dat,w); for(i=0;i12;i++) { for(j=0;j12;j++) {Time[i][j]=0.0;for(n=0;n12;n++){for(m=0;m9;m++){fl[i][j][n][m]=0;x0=5*m;x1=x0+5;y0=3*n;y1=y0+3;y_0=xl[i][j]*x0+b[i];if(y_0=y1y_0=y0){X0=x0;Y0=y_0;y_1=xl[i][j]*x1+b[i];if((y_1=y1)(y_1=y0)){Y1=y_1;X1=x1;}else{x_1=(y0-b[i])/xl[i][j];if(x_1=x1x_1=x0){Y1=y0;X1=x_1;}else{x_1=(y1-b[i])/xl[i][j];X1=x_1;Y1=y1;}}fl[i][j][n][m]=sqrt((X1-X0)*(X1-X0)+(Y1-Y0)*(Y1-Y0));}else{y_1=xl[i][j]*x1+b[i];if(y_1=y1y_1=y0){X1=x1;Y1=y_1;x_0=(y0-b[i])/xl[i][j];if(x_0=x1x_0=x0){X0=x_0;Y0=y0;}else{x_0=(y1-b[i])/xl[i][j];X0=x_0;Y0=y1;}fl[i][j][n][m]=sqrt((X1-X0)*(X1-X0)+(Y1-Y0)*(Y1-Y0));}else{x_0=(y0-b[i])/xl[i][j];if(x_0=x1x_0=x0){X0=x_0;Y0=y0;x_1=(y1-b[i])/xl[i][j];X1=x_1;Y1=y1;fl[i][j][n][m]=sqrt((X1-X0)*(X1-X0)+(Y1-Y0)*(Y1-Y0));}else fl[