基于PML吸收边界和辛算法的地震波场模拟

摘要

本文在辛算法研究的基础上,结合褶积微分算子及通量校正传输(FCT)等数值模拟技术,采用完全匹配层(PML)吸收边界对来二维介质中的地震波场进行了数值模拟.地震波场的数值模拟结果表明,辛算法不仅具有保持体系结构的特性,并且具有长时间跟踪能力,虽然在计算时增加了运算量,但具有很强的数值模拟稳定性;采用了褶积微分算子及FCT 数值模拟技术和PML 吸收边界,突出了空间微分的局部属性,压制网格频散,可以对地震波场特征进行更精确的数值模拟;利用辛算法对地震波场进行模拟,不仅可以计算各种反射波,而且可以对复杂构造所产生的绕射波波场特征进行数值模拟.