《地球物理层析成像实验报告》由会员分享，可在线阅读，更多相关《地球物理层析成像实验报告（24页珍藏版）》请在金锄头文库上搜索。
1、本 科 生 实 验 报 告实 验 课 程 地 球 物 理 层 析 成 像 学 院 名 称 地 球 物 理 学 院 专 业 名 称 勘 查 技 术 与 工 程 学 生 姓 名 学 生 学 号 指 导 教 师 曹 俊 兴 实 验 地 点 5417 实 验 成 绩 二一五 年 三 月 二一五 年 四 月学生实验 心得在学习了地球物理层析成像之后，收获了很多专业知识，比如学会了利用层析成像的手段反演出地下地质体的异常，同时也学会了利用我们的专业知识解决不同的地质问题。程序语言作为一种工具一方面起到了辅助作用，另一方面我们也学会了一种思维方式，如何设计程序，如何用程序解决我们的复杂问题。在今后的学习工作
2、当中，进一步拓宽思路，勇于创新，能够获得更多的知识。学生（签名）： 2015 年 4 月 28 日指导教师评语成绩评定：指导教师（签名）：年 月 日地震走时层析成像实验地球物理正反演概论课程结业报告学号：201205060423姓名：马力衡专业：勘查技术与工程手机：18011445091摘要运用 C 语言程序，正演得到地震走时和射线在传播过程中经过离散化处理单元格内的距离。通过反演程序反演出地下异常速度值，将反演所得速度值成图与原始速度成图进行比较，得出结论。离散化处理模型建立，单边激发,四边激发直射线正演，单边激发，四边激发反演异常值，用代数重建算法迭代慢度矩阵对单边，四边激发进行迭代。关键
3、词：离散化 反演 迭代 第 1 章 地震走时层析成像实验1.1 实验内容1.1.1 直射线正演：使用直射线追踪方法计算走时的正演；分块均匀模型1.1.1.1 单边激发正演程序：#include #include void main()int v129;int m,n，i，j;FILE *fp0;fp0=fopen(速度.txt,r);for(i=0;i=y0)X0=x0;Y0=y_0;y_1=xlij*x1+bi;if(y_1=y0)Y1=y_1;X1=x1;elsex_1=(y0-bi)/xlij;if(x_1=x0)Y1=y0;X1=x_1;elsex_1=(y1-bi)/xlij;X1=
4、x_1;Y1=y1;flijnm=sqrt(X1-X0)*(X1-X0)+(Y1-Y0)*(Y1-Y0);elsey_1=xlij*x1+bi;if(y_1=y0)X1=x1;Y1=y_1;x_0=(y0-bi)/xlij;if(x_0=x0)X0=x_0;Y0=y0;elsex_0=(y1-bi)/xlij;X0=x_0;Y0=y1;flijnm=sqrt(X1-X0)*(X1-X0)+(Y1-Y0)*(Y1-Y0);elsex_0=(y0-bi)/xlij;if(x_0=x0)X0=x_0;Y0=y0;x_1=(y1-bi)/xlij;X1=x_1;Y1=y1;flijnm=sqrt(X1
5、-X0)*(X1-X0)+(Y1-Y0)*(Y1-Y0);else flijnm=0.0;fprintf(fp_ds,%f ,flijnm);Timeij+=flijnm/vnm;/每个单元格的走时；printf(%d ,vnm);/射线传播的总时间 fprintf(fp_ds,n);/以上判定射线存在的单元格情况/FILE *fp1;fp1=fopen(走时.txt,w);for(i=0;i#include#includedouble funfpo1(double xl,double b,FILE *fp,FILE *fp1)int i,j,m,n; double fl129;double
6、X0,Y0; double X1,Y1; double x0,x1,y0,y1; double x_0,x_1,y_0,y_1;double v129,Time=0.0;for(i=0;i=y0)X0=x0; Y0=y_0;y_1=xl*x1+b; if(y_1=y0)Y1=y_1;X1=x1;elsex_1=(y0-b)/xl;if(x_1=x0)Y1=y0; X1=x_1;elsex_1=(y1-b)/xl;X1=x_1; Y1=y1;flnm=sqrt(X1-X0)*(X1-X0)+(Y1-Y0)*(Y1-Y0);elsey_1=xl*x1+b;if(y_1=y0)X1=x1; Y1=y
7、_1;x_0=(y0-b)/xl;if(x_0=x0)X0=x_0; Y0=y0; elsex_0=(y1-b)/xl;X0=x_0; Y0=y1;flnm=sqrt(X1-X0)*(X1-X0)+(Y1-Y0)*(Y1-Y0);elsex_0=(y0-b)/xl;if(x_0=x0)X0=x_0; Y0=y0;x_1=(y1-b)/xl;X1=x_1; Y1=y1;flnm=sqrt(X1-X0)*(X1-X0)+(Y1-Y0)*(Y1-Y0);fprintf(fp,%f ,flnm);Time+=flnm/vnm;fprintf(fp,n);fprintf(fp1,%fn ,Time);r
8、eturn 0;void main() int i,j;/以下是读取原始速度值/FILE *fp=fopen(每个单元格的距离.txt,w),*fp1=fopen(时间.txt,w);double x_jf4=0.0,y_jf4;/左侧激发点坐标double x_js4=0.0,y_js4;/右侧激发点的坐标double X_jf4,Y_begin4=0.0;/上下激发点坐标double X_end4,Y_js4=0.0;/下顶下激发点坐标for(i=0;i#includevoid main()double r0, d01212129, a1,a2,Time1212，md129; int N，
9、i,j,n,m;double M129;double wucha1212=0.0;/反演误差FILE *fp7,*fp8,*fp9,*fp10; fp7=fopen(每一小格的距离.dat,r); fp8=fopen(走时.txt,r); fp10=fopen(不为零的个数.txt,w);for(i=0;i12;i+) for(j=0;j12;j+) for(n=0;n12;n+)for(m=0;m9;m+)fscanf(fp7,%lf , /*for(i=0;i12;i+) for(j=0;j12;j+) for(n=0;n12;n+)for(m=0;m9;m+)printf(%5.2lfn
10、,d0ijnm); */for(i=0;i12;i+) for(j=0;j12;j+) fscanf(fp8,%lf,/*for(i=0;i12;i+) for(j=0;j12;j+) printf(%5.6lfn,Timeij);*/for(n=0;n12;n+)for(m=0;m9;m+) Mnm=0.0; mdnm=0.0;for(n=0;n12;n+)for(m=0;m9;m+) for(i=0;i12;i+) for(j=0;j12;j+) if(d0ijnm!=0)Mnm+=1; /算出系数矩阵中每一列不为零值的数的个数fprintf(fp10,%fn,Mnm); N=0; /迭代
11、次数 10000 次 /计算误差/while(N10000) for(i=0;i12;i+)for(j=0;j12;j+) r0=0.0; for(n=0;n12;n+)for(m=0;m9;m+)r0+=d0ijnm*mdnm;wuchaij=(Timeij-r0); /printf(%f,wuchaij); /重建算法迭代慢度矩阵 /for(m=0;m12;m+)for(n=0;n9;n+) a1=0.0; a2=0.0;for(i=0;i12;i+) for(j=0;j12;j+) a1+=d0ijmn*d0ijmn; a2+=(wuchaij*d0ijmn); mdmn+=a2/(a1
12、*Mmn); N+; /以 txt 格式输出最后数据/ fp9=fopen(最终结果.txt,w); for(i=0;i12;i+)for(j=0;j9;j+) fprintf(fp9,%lf ,1.0/mdij); fprintf(fp9,n); 数据；10 次2.784751 3.280478 3.261611 2.454797 3.711088 2.493625 3.307422 3.441376 2.771530 3.145613 3.042454 2.956347 2.564969 3.632866 2.761524 3.229051 3.030151 3.015298 3.0241
13、39 3.211870 3.507230 3.014290 3.717386 2.895093 3.202817 2.976226 2.986879 2.889704 3.105180 3.377056 3.284153 3.382822 2.999874 3.057579 3.020697 2.938980 2.886266 2.913581 2.910847 3.054787 3.280013 3.066888 3.096661 3.088917 2.967574 2.805710 2.954641 2.924597 2.918879 3.156665 3.081835 3.116152 3.125369 2.938478 2.843217 2.995171 2.971224 2.908500 3.053036 3.1266
