A novel multigrid based preconditioner for heterogeneous Helmholtz problems