GPU-Accelerated Transient Electro-Magnetic Modelling
Y. Liu, M. Hagdorn and A. Ziolkowski
Event name: 81st EAGE Conference and Exhibition 2019
Session: Poster: Electromagnetic Methods A
Publication date: 03 June 2019
Info: Extended abstract, PDF ( 486.59Kb )
Price: € 20
We present a GPU-accelerated transient electromagnetic modelling method implemented in CUDA C where the numerical efficiency can be improved by a factor of over 100 compared with the serial C code. The method uses the Fourier pseudo-spectral method to solve the spatial derivatives, and uses the rapid expansion method (REM) to solve the temporal evolution of the field by a summation of Chebyshev polynomials. The results are free of numerical dispersion and accurate to the Nyquist in both space and time. We show that if the model consists only of 1-D or 2-D structures, the computational effort to solve a 3-D field can be reduced by up to two orders of magnitude without loss of accuracy. We demonstrate the accuracy of the code by a 2.5-D modelling example in a VTI anisotropic half space. The significant improvement in numerical efficiency generalizes the use of REM in large-scale, time-domain EM modelling or inversion problems.