《反演与层析成像》实验 - 图文

本科生实验报告

实验课程 地球物理反演与层析成像 学院名称 地球物理学院

专业名称 地球物理学 学生姓名 学生学号 指导教师

实验地点 5417 实验成绩

二〇-四年 十二 月 二〇一五年一月

学生（签名）： 年 月 日 学生实验 心得 指导 教师 评语 成绩评定： 指导教师（签名）： 年 月 日

实验一 射线追踪-建立观测系统

一、 实验原理

1.1 初至拾取

初至拾取也就是对初至波的到达时间的记录，这是进行层析反演计算的基础数据。数据的好坏直接关系到层析成像的效果好坏。如果是进行实际资料计算时，初至拾取则为野外采集所得数据；若为理论模型的计算，则是在理论模型上进行正演计算所得的初至波旅行时。 1.2 建立初始模型

初始模型的建立，就是按照实际模型的地表起伏给定一假设的速度场。为方便起见，这个速度场可以是从地表开始以某一初始速度以相同的梯度往下递增。这一初始速度可以同先验知识给出，又或者从初至拾取中获得。炮－道比较接近时，可以把地震波看作直线传播，已知炮检距及初至时间就可以算出地表速度。初始模型的深度一定要够大，能够让射线自由地传播，不会出现遇到模型底部强迫反射的现象。 1.2 射线追踪

射线追踪的准确与否是影响层析成像的关键。鉴于日本科学家Aszkawar提出的LTI算法的高效率和高精度，本文采用该算法进行旅行时计算及射线追踪。LTI算法基于Fermat原理,即地震波沿着一条传播时间最短的路径进行传播。 该方法把模型离散成均匀的正方形单元，旅行时和射线路径的确定只与单元边界上的点有关。假设单元边界上任一点的旅行时可由该边界上相邻两个离散点的旅行时线性插值得到。如图1所示为一匀速正方形单元，单元边界平行于坐标轴，A,B是二个旅行时己知的点，要求射线穿过A,B边界到达D点的最小旅行时及射线路径。设射线从C点通过，C点旅行时可用线性内差公式由A,B二点的值表示。D点的旅行时为C点旅行时与波在C, D间直线传播时间之和。然后根据Fermat原理，就可求出C点的位置（即射线路径）和D点的最小旅行时。该方法也分成向前和向后处理，向前处理只计算各单元边界上的节点的旅行时；向后处理根据Fermat原理追踪射线路径．确定的射线路径不是单元边界上离散节点的连线，而是穿过单元边界上正好满足最小旅行时条件的那一点的线。这样，在一个单元内射线路径为直线，同时边界上的折射角随入射角连续变化。所以，

该方法具有较高的精度和计算速度。

图1匀速正方形单元旅行时线性插值几何关系 A、B 为已知点，C是内插点，D是待求点

二、 程序源代码 #include\

#include\struct data { double begin; double end; double slope; double length;

double transform[12][9]; double transform2[108]; int n;

double time; struct point {

double x; double y; }point[100]; }ray[144];

struct ww {

double x; double y;

}d[100];