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

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