Analysis of GPR Wave Propagation Using CUDA-Implemented Conformal Symplectic Partitioned Runge-Kutta Method

Accurate forward modeling is of great significance for improving the accuracy and speed of inversion. For forward modeling of large sizes and fine structures, numerical accuracy and computational efficiency are not high, due to the stability conditions and the dense grid number. In this paper, the symplectic partitioned Runge-Kutta (SPRK) method, surface conformal technique, and graphics processor unit (GPU) acceleration technique are combined to establish a precise and efficient numerical model of electromagnetic wave propagation in complex geoelectric structures, with the goal of realizing a refined and efficient calculation of the electromagnetic response of an arbitrarily shaped underground target. The results show that the accuracy and efficiency of ground-penetrating radar (GPR) forward modeling are greatly improved when using our algorithm. This provides a theoretical basis for accurately interpreting GPR detection data and accurate and efficient forward modeling for the next step of inversion imaging.