Snapshot wavefield decomposition for heterogeneous velocity media