High-efficiency Numerical Simulation of Electromagnetic Response of Complex Geoelectric Structures Based on GPU Parallel Symplectic Algorithm

For forward modeling of large size and fine structures, the numerical accuracy and computational efficiency are not high due to the stability conditions and dense grid number. In this paper, the symplectic Euler algorithm, surface conformal technique and GPU acceleration technique were combined to establish a precise and high-efficiency numerical model of electromagnetic wave propagation in complex geoelectric structures. And we realized the refined and efficient calculation of the electromagnetic response of arbitrary shape underground target. The results show that the accuracy and efficiency of ground penetrating radar (GPR) forward modeling are greatly improved, which provides a theoretical basis for accurately interpreting GPR detection data, and provides an accurate and efficient forward modeling for the inversion imaging.