A High-Efficiency Spectral Element Method Based on CFS-PML for GPR Numerical Simulation and Reverse Time Migration

被引:15
作者
Wang, Xun [1 ,2 ,3 ]
Yu, Tianxiao [1 ,2 ,3 ]
Feng, Deshan [1 ,2 ,3 ]
Ding, Siyuan [1 ,2 ,3 ]
Li, Bingchao [1 ,2 ,3 ]
Liu, Yuxin [1 ,2 ,3 ]
Feng, Zheng [1 ,2 ,3 ]
机构
[1] Cent South Univ, Key Lab Metallogen Predict Nonferrous Met & Geol E, Minist Educ, Changsha 410083, Peoples R China
[2] Key Lab Nonferrous Resources & Geol Hazard Detect, Changsha 410083, Peoples R China
[3] Cent South Univ, Sch Geosci & InfoPhys, Changsha 410083, Peoples R China
基金
中国国家自然科学基金;
关键词
Time-domain analysis; Numerical models; Finite element analysis; Numerical analysis; Computational modeling; Mathematical models; Numerical simulation; Complex frequency shifted perfectly matched layer (CFS-PML); ground penetrating radar (GPR); per-element GPU parallel framework; spectral element method (SEM); GROUND-PENETRATING RADAR; PERFECTLY MATCHED LAYER; SEISMIC-WAVE PROPAGATION; DOMAIN METHOD; EVANESCENT WAVES; SPLIT-FIELD; FDTD; BOUNDARY; ABSORPTION; PERFORMANCE;
D O I
10.1109/JSTARS.2023.3234199
中图分类号
TM [电工技术]; TN [电子技术、通信技术];
学科分类号
0808 ; 0809 ;
摘要
Improving the accuracy and efficiency of the numerical simulation of ground penetrating radar (GPR) becomes a pressing need with the rapidly increased amount of inversion data and the growing demand for migration imaging quality. In this article, we present a numerical spectral element time-domain (SETD) simulation procedure for GPR forward modeling and further apply it to the reverse time migration (RTM) with complex geoelectric models. This approach takes into account the flexibility of the finite element methods and the high precision of the spectral methods. Meanwhile, in this procedure, the complex frequency shifted perfectly matched layer (CFS-PML) is loaded to effectively suppress the echo at the truncated boundary, and the per-element GPU parallel framework used can achieve up to 5.7788 times the efficiency compared with the CPU calculation. The experiments on SETD spatial convergence and CFS-PML optimal parameter selection showed that, under the same degree of freedom, the SETD offered substantially better accuracy compared with the traditional FETD. The experiments on RTM of different profiles with different orders of SETD via a complex geoelectric model verify the universality of the algorithm. The results indicate that the RTM imaging effect has been significantly improved with the increase of SETD order. It fully proves the great potential of efficient and high-precision SETD simulation algorithm in the RTM imaging direction and shows certain guiding significance for underground target structure exploration.
引用
收藏
页码:1232 / 1243
页数:12
相关论文
共 55 条
[1]   Evanescent waves in PML's:: Origin of the numerical reflection in wave-structure interaction problems [J].
Bérenger, JP .
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, 1999, 47 (10) :1497-1503
[2]   Improved PML for the FDTD solution of wave-structure interaction problems [J].
Berenger, JP .
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, 1997, 45 (03) :466-473
[3]   Application of the CFSPML to the absorption of evanescent waves in waveguides [J].
Bérenger, JP .
IEEE MICROWAVE AND WIRELESS COMPONENTS LETTERS, 2002, 12 (06) :218-220
[4]   A PERFECTLY MATCHED LAYER FOR THE ABSORPTION OF ELECTROMAGNETIC-WAVES [J].
BERENGER, JP .
JOURNAL OF COMPUTATIONAL PHYSICS, 1994, 114 (02) :185-200
[5]   Reverse-Time Migration from Rugged Topography to Image Ground-Penetrating Radar Data in Complex Environments [J].
Bradford, John H. ;
Privette, Janna ;
Wilkins, David ;
Ford, Richard .
ENGINEERING, 2018, 4 (05) :661-666
[6]   Dedalus: A flexible framework for numerical simulations with spectral methods [J].
Burns, Keaton J. ;
Vasil, Geoffrey M. ;
Oishi, Jeffrey S. ;
Lecoanet, Daniel ;
Brown, Benjamin P. .
PHYSICAL REVIEW RESEARCH, 2020, 2 (02)
[7]   The application of finite-difference time-domain modelling for the assessment of GPR in magnetically lossy materials [J].
Cassidy, Nigel J. ;
Millington, Tim M. .
JOURNAL OF APPLIED GEOPHYSICS, 2009, 67 (04) :296-308
[8]  
CLAYTON R, 1977, B SEISMOL SOC AM, V67, P1529
[9]  
Cohen G., 2002, SCIENTIF COMPUT
[10]   Performance of regular PML, CFS-PML, and second-order PML for waveguide problems [J].
Correia, Davi ;
Jin, Jian-Ming .
MICROWAVE AND OPTICAL TECHNOLOGY LETTERS, 2006, 48 (10) :2121-2126