A parallelization of the wavenumber integration acoustic modelling package OASES

Abstract

The wavenumber integration seismo-elastic model OASES can simulate the wave propagation in layered media, consisting of rough interfaces and elastic and porous layers. Range dependence is achieved by coupling vertical sections of layers, or cuts, together using the spectral super-element method. Only the specific frequencies, receiver depths, and offsets of interest need to be calculated when using the wavenumber integration technique. However, as pulse propagation at higher frequencies is simulated, denser sampling of frequencies must be used. For complex media with many layers, and many vertical sections, the computation time quickly escalates. By exploiting that each frequency response can be calculated independently, a parallelization of the OASES package has been implemented and is presented here. This makes otherwise computationally unfeasible or unpractical simulations feasible. In this work, the parallel OASES package is applied to, and benchmarked on, several acoustic and seismic problems. The increased computation capacity is used to simulate and image the full wave field of several cases, reducing the computation time in one of the cases from 1.5 years to 5 h.