叶 佩，李庆春

1.中国国土资源航空物探遥感中心，北京 100083

2.长安大学地质工程与测绘学院，西安 710054

0 引言

射线理论及射线方法是研究地震波传播理论的重要途径之一，射线理论也经常被用于研究不均匀介质和复杂构造中的地震波传播问题。利用旅行时信息进行层析成像最常用的正演方法就是射线追踪，该方法还被广泛用于偏移、地震波反演及其他正演模拟中，计算地震波的旅行时和波前［1－11］在实际应用中得到了不断的发展与完善。旅行时线性插值（linear travel－time interpolation，LTI）射线追踪法是Asawaka［4］提出来的一种基于线性假设的网格单元扩展方法，该方法已广泛应用于层析成像，具有计算速度快、精度高、易于实现等优点。近年来，国内很多学者对该方法做了大量的研究。如：张建中等［6］推导出了旅行时非线性插值及变速介质中旅行时插值计算公式；为了改善该方法在近场旅行时的精度，降低累计误差，聂建新等［12］提出了旅行时二次／线性联合插值（简称 QLTI）；彭直兴等［13］在震源的近场采用非线性插值计算旅行时，在远场仍用线性插值的计算方法，改善了传统的LTI法因线性插值引起的误差；为了改善旅行时的计算精度和计算效率，孙章庆等［14－15］利用窄带技术来模拟波前，并在局部采用线性插值计算走时，随后，又将该方法扩展到复杂地表走时计算。由于传统的LTI法在进行旅行时计算时射线方向考虑不全，计算的节点旅行时不一定最小，导致追踪的射线路径不一定满足最小旅行时。……