变网格步长声波方程有限差分数值模拟.pdfVIP

变网格步长声波方程有限差分数值模拟 摘要：有限差分算法是波动方程正演模拟最常用的方法之一。其优点 是包含的地震信息丰富且算法简单容易实现。传统的常网格有限差分 算法在模拟含有低速地表或低／高速夹层时，常常存在不稳定或者频 散现象，为了得到精度较高的模拟结果，就需要减小网格步长，这样 就大大降低了计算效率和浪费了计算机资源。 本文首先介绍了基于声波方程的常规高阶有限差分正演模拟方 法，详细讨论了有限差分正演模拟方法的稳定性和数值频散问题，使 用完全匹配层(PML)边界条件作为有限差分正演模拟的边界条件， 并结合具体的理论模型进行了试算，结果表明完全匹配层吸收边界吸 收效果好，简单容易实现。 为了很好地解决上述问题，本文在前人的研究基础上，发展和改 进了变网格算法，该算法不需要进行波场插值，也不需要用平滑函数 平滑，实现简单，计算效率高且与常网格算法相比能节省大量的计算 机资源。 为了验证该变网格算法的有效性，本文设计了五种不同的地质体 模型，并对这些模型进行了试算，分别从模拟精度、计算效率和内存 需求三个方面与常网格有限差分算法进行了对比分析，结果表明纵向 变网格步长算法明显优于常网格有限差分算法。该算法可以视为对传 统有限差分算法的有效的改进。 图25幅，表6个，参考文献60篇。 关键词：有限差分；数值频散；稳定性条件；变网格；计算时间； PML边界条件 分类号：P631 II Numerical ofAcousticWave with Modeling EquationSpatial Meshes Varying Abstract：Finitedifferenceisoneofthemost methodsforthe popular forward ofwave itsrecord contain modeling equation，sincemay abundantinformationandthe iS toberealized．The algorithmeasy traditionalfinjtedifference hasthe of and algorithmproblemsinstability with surface numerical when withamodel dispersiondealing low-speed mediumor simulationresultswith low／highspeedinterlayer．Toget normal of is method stepsadopted． higherprecision，a decreasing鲥d the time andwastes makes However,it computinglonger computer memory． differencemethodbasedon Theconventi