石油物探2025，Vol.64Issue(4)：639-652,14.DOI:10.12431/issn.1000-1441.2024.0040
复杂地质条件的间断有限元地震波数值模拟及GPU加速
A discontinuous Galerkin finite-element method of seismic modeling in complex media using GPUs
摘要
Abstract
The discontinuous Galerkin finite-element method(DGFEM),which is a high-order finite element method adapting to complex surface conditions,has attracted extensive attention.Based on triangular unstructured meshes and local Lax-Friedrichs flux,the matrix forms of DGFEM calculation using elastic,viscoelastic,and poroelastic wave equations are established,and the general calculation format for single wave field components is developed,which improves the scalability of DGFEM programming.Based on this format,the procedure to construct a universal CUDA kernel is developed,which can be easily extended for more complex media and 3D cases,and the CPU+GPU parallel computing framework of 2D DGFEM is established.The results of a theoretical model and a complex mountain model reveal that the general calculation format and CUDA kernels constructed in this paper can accurately simulate P-waves,S-waves,and slow P-waves described by using acoustic,elastic,viscoelastic,and poroelastic wave equations.Compared to single-core CPU simulation,the speedup ratio of 2D DGFEM elastic-wave GPU calculation is about 100 on the average.Meanwhile,the simulation time for elastic,viscoelastic,and poroelastic waves is approximately 1.7,2.3,and 3.0 times that of acoustic wave simulation,respectively.This result can be used to guide multi-process load balancing in the simulation of complex coupled media.关键词
间断Galerkin有限元方法/弹性波/粘弹性波/孔隙弹性波/数值模拟/GPU并行计算Key words
discontinuous Galerkin finite-element method/elastic wave/viscoelastic wave/poroelastic wave/numerical simulation/GPU parallel computing分类
能源科技引用本文复制引用
韩德超,刘卫华,张春丽,袁媛,白鹏..复杂地质条件的间断有限元地震波数值模拟及GPU加速[J].石油物探,2025,64(4):639-652,14.基金项目
国家自然科学基金企业创新发展联合基金(U19B6003-04)资助.This research is financially supported by the NSFC Enterprise Innovation and Development Joint Funds(Grant No.U19B6003-04). （U19B6003-04）
