基于卷积滤波的谱元法在长时程波场模拟中的应用
Application of spectral element method based on convolution filtering in long-term wavefield modeling作者机构：中国科学院大学计算地球动力学重点实验室北京 100049
出 版 物：《地球物理学报》 (Chinese Journal of Geophysics)
年 卷 期：2024年第67卷第5期
页 面：1832-1838页
核心收录：
学科分类：070801[理学-固体地球物理学] 07[理学] 0708[理学-地球物理学]
基 金：41725017:国家自然科学基金 2020YFA0713400:国家重点研发计划
主 题：长时程波场模拟 卷积滤波 谱元法 CFL条件 数值稳定性
摘 要：地震波波场数值模拟是深入认识复杂地下介质的重要手段.谱元法兼顾了有限元方法的传统优势,并提高了数值精度和稳定性,已成为目前地震波波场数值模拟中最常用的方法之一.随着超大规模并行计算的普及和发展,高分辨模拟所需的网格尺寸越来越小,网格规模越来越大,导致长时程模拟的大步长需求越来越高,成为倍受关注的研究热点.区别于传统的保结构算法,本文采用时间频散变换对时间离散引起的误差进行补偿,并引入空间滤波用于打破CFL(Courant-Friedrichs-Lewy)条件对时间采样间隔的限制.为了克服谱元法中因为非均匀空间采样点引起的滤波频谱计算难题,本文采用空间卷积滤波替换传统频率域滤波,有效的降低了算法的计算需求.该滤波方法也适用于复杂起伏模型的波场模拟.将本文方法用于保结构算法,并针对四阶辛Nyström方法进行了测试,显著增大了四阶辛Nyström方法的时间采样步长.新方法可以有效的用于地球自由振荡的模拟以及其他需要长时程数值模拟的工作中.