论文部分内容阅读
地下介质中波的传播情况的研究,对于地质勘探具有十分重要的应用价值,为了更好的模拟地下波场,需要不断地发展和完善解决问题的方法.而波动方程正演作为正问题研究的重要分支,在众多领域有着重要的应用.例如,我们可以用弹性动力学中的各种波动方程来描述地震波在地下介质的传播,并且通过波动方程的正演模拟来获取地震资料.对于波动方程正问题而言,它面临着构造的显格式对步长有一定的限制、数值解不适定等困难,并且实际应用中会有相当大的计算量和存储量.因此,波动方程正演的研究具有重要的理论意义及广泛的应用价值.对于二维波动方程的研究为更高维度波动方程的研究奠定了理论基础. 首先,本文介绍正演模拟的主要方法,包括有限差分法、有限元法以及伪谱法.其中,对波动方程的初边值问题以及有限差分法三种常用差分格式进行简单介绍,包括古典显格式、古典隐格式以及VonNeumann格式;对于有限元法,介绍其找近似解的主要六步过程;对于伪谱法,介绍中心差分和时间导数,以及傅里叶变换和空间导数. 其次,以二维波动方程作为数学模型.然后应用显式格式与隐式格式的有限差分法对模型进行离散,构造十一点隐式格式和十五点隐式格式,并且采用不同的方式来处理边界条件,然后验证算法的收敛性.接下来,使用MATLAB对构造的显式差分计算格式进行数值模拟,分析模拟效果,说明第二种和第三种边界条件处理要优于一般的边界条件处理方法. 最后,在二维波动方程的基础上,研究三维波动方程的数学模型.其中对三维波动方程进行离散,构造十五点隐式格式和二十一点隐式格式,并采用不同的方式来处理边界条件,然后验证算法的收敛性.接下来,使用MATLAB对构造的显式差分计算格式进行数值模拟,分析模拟效果,说明在三种模型下,虚拟辅助边界层法处理要优于其它两种方法.