摘要
程函方程和输运方程被广泛应用于地震波的模拟和反演问题,但在实际应用中由于点震源附近的程函方程和输运方程存在奇点,在变速介质中无法直接求解.本文采用参数化的方法将走时和振幅的奇异部分分离出来,剩下的平滑部分满足一个新的程函方程和输运方程.求解这个改造后的程函和输运方程,采用Lax-Friedrichs哈密顿函数构造迭代系统,同时利用WENO高阶差分格式离散化计算,可以提高走时和振幅计算精度.此方法可采用雅可比迭代方法实现GPU加速计算.最后,用实例证明本文方法是旅行时和振幅的高精度近似解,在计算效率方面达到50倍的加速比.