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 条
[21]  
Komatitsch D, 1999, INT J NUMER METH ENG, V45, P1139, DOI 10.1002/(SICI)1097-0207(19990730)45:9<1139::AID-NME617>3.0.CO
[22]  
2-T
[23]   Simulations of ground motion in the Los Angeles basin based upon the spectral-element method [J].
Komatitsch, D ;
Liu, QY ;
Tromp, J ;
Suss, P ;
Stidham, C ;
Shaw, JH .
BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA, 2004, 94 (01) :187-206
[24]   Spectral-element simulations of global seismic wave propagation - I. Validation [J].
Komatitsch, D ;
Tromp, J .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2002, 149 (02) :390-412
[25]   Spectral-element simulations of global seismic wave propagation - II. Three-dimensional models, oceans, rotation and self-gravitation [J].
Komatitsch, D ;
Tromp, J .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2002, 150 (01) :303-318
[26]  
Komatitsch D, 2005, GEOPHYS MONOGR SER, V157, P205, DOI 10.1029/156GM13
[27]   Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers [J].
Kuzuoglu, M ;
Mittra, R .
IEEE MICROWAVE AND GUIDED WAVE LETTERS, 1996, 6 (12) :447-449
[28]   A 3-d spectral-element method using mixed-order curl conforming vector basis functions for electromagnetic fields [J].
Lee, JH ;
Xiao, T ;
Liu, QH .
IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, 2006, 54 (01) :437-444
[29]   A 3-D Discontinuous Spectral Element Time-Domain Method for Maxwell's Equations [J].
Lee, Joon-Ho ;
Chen, Jiefu ;
Liu, Qing Huo .
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, 2009, 57 (09) :2666-2674
[30]   GPR simulation based on complex frequency shifted recursive integration PML boundary of 3D high order FDTD [J].
Li, Jing ;
Zeng, Zhaofa ;
Huang, Ling ;
Liu, Fengshan .
COMPUTERS & GEOSCIENCES, 2012, 49 :121-130