地球自由振荡的保结构模拟

大地震等诸多激励均能激发全球自由振荡现象,通常表现为驻波形式的全球整体振荡.现有的地震波数值模拟方法多为非保结构方法,无法压制长时程计算中的积累误差.本文采用优化的三阶辛格式谱元法,对地球自由振荡及全球尺度的地震波传播进行了长时程模拟.通过与传统的基于Newmark算法的谱元方法结果对比分析,明确验证了本文所得优化的三阶辛格式谱元法在模拟地球自由振荡等大规模长时程问题上的优越性和准确性.上述进展在方法论层面为今后探测、刻画全球尺度地球非均匀结构的驻波数值方法奠定了部分基础,并为相关研究领域提供了新的选择.     

求解声波方程的辛RKN格式

将声波方程变换至Hamiltion体系,构造了适用于高效声波模拟的二阶显式辛Runge-Kutta-Nyström（RKN）格式,运用根数理论得到此格式的阶条件方程组. 针对两个自由度的辛条件方程组,根据三次项截断误差最小原理得到一种误差最小辛格式;通过分析声波的时间演进方程的稳定性,选择不同的辛系数使演进方程更稳定,并得到了另一种更为稳定辛格式;在频散关系分析中,选择使数值频散最小的辛系数,得到第三种最小频散辛格式. 在理论分析中,这组辛RKN格式相比常见格式在精度控制、数值频散压制以及稳定性提升等方面均具有明显优势;在数值实验中,通过具体算例验证了理论分析的正确性.     

正则化预处理迭代算法在频率域声波模拟中的应用

高精度及高效频率域声波数值模拟的关键在于高效求解声波方程经离散化后得到的大型稀疏线性方程组.该方程组系数矩阵具有很强的稀疏性,非对称性和非正定性等特征,常用的迭代算法难以准确、高效地求解.为了改善数值模拟迭代算法的收敛性与稳定性,在算法基础上添加预条件算子是求解该类方程的常用方案.本文基于以上思路,引入正则化技术来构造合适的预条件算子,提出正则化预条件迭代算法,以加速求解方程组.通过包含有均匀介质和高非均匀度介质(Marmousi)模型的数值模拟实验结果表明:与单独使用迭代算法相比,本文提出的正则化预条件迭代算法在计算量方面仅多了一次矩阵-矢量相乘,内存消耗未增加;同时,基于该算法的数值模拟结果能够满足精度要求,较单独使用迭代法能够有效改善收敛性质,加快收敛速度;而且,在二维模型算例下,与LU分解算法相比,基于该算法的内存消耗大幅下降.