3D adaptive internal multiples modeling based on globally optimised Fourier finite-difference method. (27th July 2021)