地震层析成像概论讲稿.docxVIP

《地震层析成像概论》大作业张义蜜，2012260301272016-01-041简述用于地震走时成像方法中的射线追踪算法及原理。1.1打靶法地震射线的打靶法原理简单，它是由射线的运动学方程(1式)定义的初值问题，可用来进行完全射线的追踪计算(存在速度界面时采用斯奈尔定律)。在寻找炮点和检波器点间的地震射线的两点问题是一个反演问题，其中未知参量是射线的初始方向，函数取最小值的原则是使射线终点到检波器间的距离为最小。这类算法的主要挑战是反演问题的非线性行为，且这种非线性行为随介质的复杂程度急剧增加(见图1所示)。（1）图2是打靶法原理图1.1.1近(旁)轴射线追踪射线理论的一个还未提及的重要领域是近轴射线近似，它被广泛地用来进行数据预测。1.1.2完全非线性打靶算法在弱非均匀介质中，非线性迭代打靶算法解决边值问题是十分有效的，但当介质复杂程度增加时则变得不稳定。图3说明了这种情况，展示出具有不同复杂程度的两个速度模型中射线倾角与射线终点到检波器间距离的关系。给定非线性迭代打靶算法易出的错误，完全非线性算法至少是值得研究的。然而，相关文献较少，或许是近来太多的基于网格节点和波前算法出现就是用来克服上述局限。1.2弯曲(调整)法射线追踪调整算法的原理就是不断循环调整连接炮点与检波器间任意初始路径直至成为真实射线路径为止(即：满足费马原理所说的稳定时间路径，见图4所示)。通常的做法是推导出可以循环求解的射线运动学方程的边值形式。1.2.1伪弯曲法伪弯曲法原理上与上述弯曲算法相似，但不对射线方程直接求解。该算法较早提出之一是Um和Thurber(1987)，主要基于射线可由一组线性插值点表示。给定某些任意路径，其目的是不断调整每个节点的位置，从而满足射线方程。利用确定射线路径的法方向，然后直接运用费马原理而有效地实现上述算法。??????Zhao等(1992)修改了Um和Thurber(1987)的三点扰动算法，且容许速度界面的存在。这种情况下，定义射线的点序列包括射线穿过每个界面上的点。当更新至界面上这些点时，将沿界面扰动(此时，连接界面两侧的点不动)直到满足斯奈尔定律为止。Koketsu和Sekine(1998)在球坐标系下提出了相似的算法。1.2.2其它弯曲算法除了弯曲和伪弯曲算法外，还有其他类似的弯曲算法也是值得一提的。Prothero等(1988)提出了基于函数极小值的三维弯曲算法。利用大量的搜寻找到炮点与检波器点间最小时间圆弧作为初始射线路径。利用三角谐函数的求和对射线进行扰动，去寻找产生走时最小的振幅系数。Debski和Ando(2004)开发了一种所谓“光谱射线追踪器”，其与Prothero等(1988)算法相近，不同的是射线由Chebyshev多项式定义。弯曲问题可形成函数求极小值问题，其中Chebyshev多项式分项系数变成可调整的变量直至两点间走时为最小。其中不是利用线性化方法去寻求最小值，而是利用遗传基因算法去产生和选择多项式分项系数。射线追踪中大多数求解边值问题可归类为打靶法和弯曲法。然而，其它的算法也是存在的，较著名的是基于结构扰动理论(Cerveny,2001)。在这类算法中，假设在参考介质中已知两点间的射线路径，其目的是要通过对参考介质的扰动确定真实介质中相应两点的射线路径。1.3基于网格(节点)波前扩展的算法炮点与检波器间的射线追踪的另一类选择算法是计算由网格节点定义介质中所有节点上的走时和波前模拟。完整的走时场隐含了波前面是时间[即：T(x)的等值线]和所有可能路径(由?T定义)的函数1.3.1快速行进法(FastMarchingMethod)当波前自身相交(也就是说多值路径存在)时，初至波前自身包含绞缠或不连续，通常与后续波前相交。由于通常只计算初至波走时，这些信息被丢弃。为了保证计算的稳定性，计算新走时时，一般不用不连续介质两面不同的差分节点。克服这一问题方法之一是：求解粘性程函方程，其中对不连续进行了平滑。这种平滑解的局限是仅限于初至波波前的计算。业已证明，上述粘性局限解同样可利用熵满足算子求解方程(18)得到，其中估计?T时考虑了方向性1.3.2最短路径算法最短路径算法是在网格化节点速度场中计算所有节点走时的另一种流行方法不需求解差分方程，而是由网格节点的连线作为具有走时的射线段。采用类似Dijkstra算法来求取给定点到网格上所有节点的最短路径。根据费马原理两点间最短路径(最小走时路径)相应于实际射线路径。最短路径网格通常是由单元或中心节点来参数化。Nakanishi和Yamaguchi(1986)把速度场参数化为由常速度单元组成，而节点定义在单元边界上(参见图26a)。该算法的优势在于每对节点间的走时可以很容易地计算(t?d?s)，这里d是两节点间的距离，s是含节点单元的波慢度。计算精度可通