快速高精度的频率空间域声波数值模拟方法研究

频率空间域地震波正演优势明显:无累计误差;多炮反演时,只需几个频率就能得到比拟时间域全带宽反演的结果;容易实现时间域难以精确实现的吸收衰减数值模拟.频率域波动方程被离散化而形成稀疏、非对称、不定、复系数、大规模的系数矩阵,需要求解大规模的线性方程组,最大计算瓶颈是需要海量的计算机内存、导致计算量庞大.采用嵌套剖分法存储稀疏矩阵,极大地减少内存需求和计算量.提出了频率空间域声波多炮模拟快速计算流程,提高了模拟的计算速度.数值试验证明嵌套剖分法具有节省存储量、计算速度快、计算效率高等优势,为后续的频率空间域声波全波形反演奠定很好的基础.

快速精确的二维频率空间域弹性波数值模拟

二维频率空间域的数值模拟方法具有以下的优势:多炮模拟时,计算成本比时间域方法低;无累计误差;在地震反演中处理多震源模拟时,只需要有限的几个频率就可以得到好的反演结果.差分离散化形成的稀疏系数矩阵,需要求解一个巨大规模的线性方程组,最大瓶颈是需要海量的计算机内存,导致计算量庞大.本文在前人研究的基础上,采用嵌套剖分网格排序法,极大限度减少对计算机内存的需求,从而减少了计算量.针对弹性波数值模拟的特征,提出二维频率空间域弹性波多炮模拟的快速计算流程.数值模拟试验证明使用嵌套剖分排序法的弹性波多炮数值模拟比压缩存储法具有节省存储量、计算效率高等优势,为后续的二维频率空间域弹性波全波形反演奠定了很好的基础.