三维弹性波方程的改进近似解析离散化方法及波场模拟

【摘要】： NAD类算法是一类新的求解波动方程的有效算法,该类算法在时间层上做截断的Taylor展开,利用偏微分方程将关于时间的高阶偏导数转化为关于空间的高阶偏导数,然后使用网格点处的函数值及其梯度共同逼近空间高阶偏导数。理论和数值实验都表明这类算法可以更有效地保留波场信息,压制数值频散,提高计算的精度和效率。 本文分别在理论分析和数值计算两个方面系统地研究了NAD类算法,主要成果有以下几个方面:改进了三维ONADM方法,得到更好的计算效果,并给出它的稳定性条件;将二维的RNADM推广到三维,完成该算法在三维情形下的精度分析,给出了三维RNADM方法的稳定性条件;发展了SNADM方法,给出了三维SNADM方法的计算精度并对SNADM方法进行了频散分析。 和经典的4阶LWC格式相比:NAD类算法在粗网格下可以有效的压制数值频散,得到清晰的波场快照;在达到消除数值频散要求时,4阶LWC格式的计算时间是ONADM计算时间的559倍,是RNADM方法计算时间的1222倍,是SNADM方法计算时间的1955倍;而此时ONADM方法存储量只有4阶LWC存储量的3.3%,RNADM方法和SNADM方法的存储量只有4阶LWC格式存储量的5.8‰。