论文部分内容阅读
数值方法求解波动方程一直都是地球物理学界和计算数学界重要的研究课题。数值方法模拟地震波传播现象,可以为天然地震预测和地震波勘探提供理论依据。现有的数值方法使用非常细化的网格,可以得到很好的地震波传播模拟效果,但是采用非常细化的网格计算时,必然会增加计算时间和计算存储量。研究高效的数值方法模拟地震波传播现象,成为求解地震波动方程的关键问题。由于地下介质分布比较复杂,介质之间存在着波速变化,给地震波传播模拟造成了很大困难。高分辨率有限体积法可以得到波在不连续介质中传播的高分辨率的数值解,本文使用高分辨率有限体积法求解声波方程和弹性波方程,构造了波传播的算子,并对波传播算子进行修正限制;数值算例表明:该方法不但可以得到高分辨率的波场快照,而且有效地控制了数值振荡。通过数值解和精确解比较、频散分析,证明了高分辨率有限体积法模拟震波传播的有效性。为了提高地震波传播模拟的计算效率,本文提出了使用自适应网格细化算法结合高分辨率有限体积法求解地震声波和弹性波方程。首先,使用高分辨率有限体积法得到基网格上的数值解,接着,对数值解进行截断误差估计,根据预先设定的阈值,确定需要进行网格细化的区域,将得到的细化网格嵌套在基网格上。自适应网格细化算法自动生成层层嵌套的细化网格,时间步长与空间步长同时细化,采用树状数据结构存储网格点。本文分别使用自适应网格和非自适应网格求解地震声波和弹性波方程,并对计算时间和计算存储量进行比较,结果表明自适应网格细化算法的高效性。提升算法得到的第二代小波变换相对于小波变换运算速度更快,并且保持了多尺度分解的特性。本文利用第二代小波的快速变换和多尺度分解性质,构造了多尺度层层嵌套的自适应网格求解二维声波方程,网格点的空间导数使用有限差分法得到,并给出了声波数值传播的波场快照、计算精度的分析和计算效率的比较。数值实验结果,说明了基于第二代小波的自适应网格求解声波方程的可行性。