扩展矩形盒波前法与扩展波前法在初至波计算中的对比研究

基于Ｖｉｄａｌｅ“扩展矩形盒波前法”的思想，通过有限差分求解程函方程，利用ＭＡＴＬＡＢ编程计算了二维速度模型的初至时间。理论研究和模拟结果表明：扩展矩形盒波前法计算初至时间不会像射线追踪法一样出现盲区，具有较高的计算精度，可大大降低计算量；但在处理复杂结构模型和强速度界面时会出现不稳定现象，导致所计算的走时并非真正的初至时间。为了进一步提高计算精度、算法的稳定性和解决地震波传播的因果性问题，给出了结合局部算法的扩展波前法结果。由于扩展波前法结合局部算法遵循了因果性，适应复杂结构模型和强速度界面，但计算的大量时间要用于寻找波前上的最小初至时间，效率较低。