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, elastic and porous layers. The wavenumber integration technique requires only the interesting frequencies, receiver depths, and offsets, to be calculated. As we seek to model pulse propagation at higher frequencies, denser sampling of frequencies must be applied. 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 achieved and is presented here. This makes otherwise unfeasible or unpractical simulations computationally achievable. 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 case from 1.5 years to 5 hours.The wavenumber integration seismo-elastic model OASES can simulate the wave propagation in layered media, consisting of rough interfaces, elastic and porous layers. The wavenumber integration technique requires only the interesting frequencies, receiver depths, and offsets, to be calculated. As we seek to model pulse propagation at higher frequencies, denser sampling of frequencies must be applied. 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 achieved and is presented here. This makes otherwise unfeasible or unpractical simulations computationally achievable. 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 case from 1.5 years to 5 hours.