Modelling ground-penetrating radar wave propagation using graphics processor unit parallel implementation of the symplectic Euler method
H. Fang, J. Lei, J. Zhang, J. An and F. Wang
Journal name: Near Surface Geophysics
Issue: Vol 17, No 4, August 2019 pp. 417 - 425
Info: Article, PDF ( 1.17Mb )
Price: € 30
Inversion of ground-penetrating radar signals requires accurate and efficient for-ward modelling. The symplectic Euler method promises good results when simulating ground-penetrating radar wave propagation in substructures, but its computational efficiency is limited by the same Courant–Friedrichs–Lewy stability condition as the finite-difference time-domain method. A two-dimensional graphics processor unit–accelerated parallel symplectic Euler algorithm is used to simulate ground-penetrating radar wave propagation. We compared the reflection waveforms as well as the sim-ulation time of the complex underground structure models simulated by the parallel symplectic Euler method with traditional finite-difference time-domain method. Re-sults show that the parallel symplectic Euler algorithm achieves the same level of accuracy as the standard finite-difference time-domain method. Moreover, it signifi-cantly improves the computational efficiency, as the calculation speed is improved by more than 21 times. We verify the performance of the proposed algorithm through a map of the single-track radar data for a three-layered pavement model and a simula-tion wiggle map for a structural damage pavement model. This provides a theoretical basis for accurately interpreting ground-penetrating radar detection data and efficient forward modelling for the next step of inversion imaging.