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

/才声Z乂庶 本 科生实验报告 实验课程 地球物理层析成像 学院名称地球物理学院 学院名称 地球物理学院 专业名称 勘查技术与工程 学生姓名 学生学号 指导教师 曹俊兴 实验地点 5417 实验成绩 二。一五年三月二。一五年四月 二。一五年三月 二。一五年四月 在学习了地球物理层析成像之后，收获了很多专业知识，比如 学 生 实 验 心学会了利用层析成像的手段反演出地下地质体的异常，同时也学会 了利用我们的专业知识解决不同的地质问题。程序语言作为一种工 具一方面起到了辅助作用，另一方面我们也学会了一种思维方式， 如何设计程序，如何用程序解决我们的复杂问题。在今后的学习工 作当中，进一步拓宽思路，勇于创新，能够获得更多的知识。 学 生 实 验 心 得 学生（签名）： 指 导 教 师 评 语成绩评定：指导教师（签名）：年 月 日2015年4月 28日 指 导 教 师 评 语 成绩评定： 指导教师（签名）： 年 月 日 地震走时层析成像实验 《地球物理正反演概论》课程结业报告 学号：201205060423 姓名：马力衡 专业：勘查技术与工程 手机摘要 运用c语言程序，正演得到地震走时和射线在传播过程中经过离散化处理单元格 内的距离。通过反演程序反演出地下异常速度值，将反演所得速度值成图与原始 速度成图进行比较，得出结论。离散化处理模型建立，单边激发，四边激发直射 线正演，单边激发，四边激发反演异常值，用代数重建算法迭代慢度矩阵对单边, 四边激发进行迭代。 关键词：离散化反演迭代 第1章地震走时层析成像实验 *3.0 * 3.0 3.0 夜) 3.0 -% 3.0 *3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 ； *3.0 3.0 5.0 3.0 3.0 3.0 3.0 3.0 3-0 i *30 3.0 5.0 3.0 3.0 3.0 3.0 3.0 3.0 *30 30 3.0 3.0 3.0 3.0 3.0 3.0 3.0, 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3Q *3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 *3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 1 3.00 *3.0 3.0 3.0 3.0 3.0 2.0 2.0 3.0 3-01 ?3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 *3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 ? 3.0 3.0 ? 3.0 L? 3.0 3.0 ? 3.0 ? 3.0 ? ? 实验内容 单元数：9x12 单元边长：3.0x5. 0m 参数:速度＜ km/s) 实验一：单边激发单边接收 (左发右收)数据：12 x 12 实验二：四边激发,四边接收 (每边设置4个激发点.在其它 三边所有接收点接收：上下两边 的激发点位置参见紫色匿叉图标 所示；直立两边的激发点位置参 见蓝色星状图标所示) 数据：4x(4x12)+2x(4 x 9 )+ 2 x(4 x12)+ 4x(4 x 9) 直射线正演： 使用直射线追踪方法计算走时的正演；分块均匀模型 单边激发正演 程序： Sinclude stdio. h Sinclude math. h void main() ( int v[12] [9]; int m, n, i, j; FILE *fpO; fpO=fopen(“速度.txt, r); for(i=0;i12;i++) ( for(j=0;j9;j++) fscanf (fpO, %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]=L 5+3. o*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);〃斜率 printfC%f\n, xl[i][j]); } for(i=0;i12;i++) ( b[i]=L5+i*3;〃每条射线截距 ) 〃以上在求射线的斜率和射线在纵轴上的截距〃 double ft_t=O. 0; 〃每一格的时间； double fl[12] [12] [12] [9]; 〃每一格射线的长度; double Time[12] [12]; 〃每条射线的时间； double X0,Y0; 〃第一个点坐标； double XI,