GeoInvexam勘查技术与工程培训讲义.ppt

Prof. Dr. Junxing Cao Prof. Dr. Junxing Cao : Geophysical Inversion * 实验：线性方程组解算 SIRT，CG/LSQR，SVD * 地球物理反演概论课程结业考察题 1. 地震走时层析成像实验（勘查技术与工程专业） 单元数：9 ? 12 单元边长：3.0 ? 5.0m 参数：速度（km/s) 实验一：单边激发单边接收 (左发右收) 数据：12 ? 12 实验二：四边激发，四边接收 (每边设置4个激发点，在其它三边所有接收点接收；上下两边的激发点位置参见紫色圈叉图标所示；直立两边的激发点位置参见蓝色星状图标所示) 数据: 4?(4?12)+2?(4 ? 9）+ 2 ?(4 ?12)+ 4?(4 ? 9) 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 5.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 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 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 2.0 2.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 3.0 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? * 内容与要求 一、实验内容 1. 直射线正演： 使用直射线追踪方法计算走时的正演；分块均匀模型 2. 反演（矩阵方程求解）： 分别使用SIRT 和CG/LSQR两种方法 3. 编程语言： Matlab or Fortran 二、 报告内容： 格式： 1. 报告题目：地震走时层析成像实验 ——《地球物理反演概论》课程结业报告（副标题） 2. 学号_______ 姓名 _______ 专业 ________手机_______ _ 内容： （1） 程序功能说明；（ 2 ) 可编译程序源代码及附带数据； （3）实验结果，图及文字说明（至少包含五幅图） 三、提交 时间与方式： 1. 提交时间：4月30日前 2. 提交 方式：电子邮件（caojx@cdut.edu.cn） 3.打包文件标题《地球物理反演概论》-勘查技术与工程-学号-姓名 地球物理反演概论课程结业考察题 2. 密度异常体反演实验（勘查技术与工程专业） R=10m △?=0.25g/cm3 H=50m 地面 R=10m △?=0.20g/cm3 H=50m Model 1 H=70m H=90m R=10m △?=0.25g/cm3 R=10m △?=0.15g/cm3 Model 2 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?30m? ?20m? * 内容与要求 一、实验内容 1. 单异常体二维正、反演； 2.多异常体二维正、反演 3. 编程语言： Matlab or Fortran 说明: (1) 图示为异常体的 几何和物性参数及球心埋深；（2）数据观测点置于地面，间距2.5m, 共80观测点（数据），以（中间）异常体地面投影为中心对称配置；（3）离散单元尺度5?5m , 离散单元个数 40(水平)?30（垂直） ，注意成像范围与扩边。(4) 使用分块均匀模型; (5)分别使用SIRT 和CG/LSQR两种方法解算矩阵方程。 二、 报告内容： 格式： 1. 标题：1. 报告题目：密度异常反演实验 ——《地球物