一种快速求解程函方程的间断伽辽金方法  

A discontinuous Galerkin method to solve eikonal equation fast

作　　者：王静男 孙章庆[1] 孙建国[1] 崔栋 韩复兴[1] 高正辉[1] WANG Jingnan;SUN Zhangqing;SUN Jianguo;CUI Dong;HAN Fuxing;GAO Zhenghui(College of Geo-exploration Science and Technology,Jilin University,Changchun 130026,China;Research Institute of Petroleum Exploration&Development,Department of Geophysics,PetroChina,Beijing 100083,China)

机构地区：[1]吉林大学地球探测科学与技术学院,长春130026 [2]中国石油勘探开发研究院油气地球物理研究所,北京100083

出　　处：《世界地质》2021年第4期935-944,共10页World Geology

基　　金：国家重点研发计划课题(2019YFC0312004);国家自然科学基金项目(42004127、42074150)资助。

摘　　要：作为一种在勘探地震学里近期发展起来的有限元类数值求解程函方程的算法,间断伽辽金方法通过拓展基函数空间提高数值精度,但存在计算效率低的问题。为了解决该问题,提出一种并行快速扫描间断伽辽金方法。算法在快速扫描过程中对节点的计算顺序进行Cuthill-McKee变换,将能够并行计算的节点分为不同任务集合,进而实现不同任务集的并行计算。数值结果的分析表明:并行快速扫描间断伽辽金方法在保证计算精度的前提下,能稳定地提升计算效率,并且还能适应各种复杂速度模型。As a recently developed finite element numerical algorithm for solving eikonal equations in exploration seismology,discontinuous Galerkin method can be used to improve numerical accuracy by extending basis function space,but there exists a problem of low computational efficiency.To solve this problem,for discontinuous Galerkin method obtained by combining fast sweeping method,Cuthill-McKee transformation of calculating order is performed in fast sweeping process and nodes which can be calculated at the same time are divided into different task sets,so as to realize the parallel calculation of task sets.The numerical analysis and examples show that the algorithm improved calculating efficiency while taking account of accuracy,and adapted in conditions with complex velocity models.

关 键 词：程函方程 间断伽辽金 快速扫描法 Cuthill-McKee排序 并行算法 

分 类 号：P631.4[天文地球—地质矿产勘探]