Full text loading...
-
GPU-Accelerated Transient Electro-Magnetic Modelling
- Publisher: European Association of Geoscientists & Engineers
- Source: Conference Proceedings, 81st EAGE Conference and Exhibition 2019, Jun 2019, Volume 2019, p.1 - 5
Abstract
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.