地震波波动方程数值模拟方法