Simulation method for determining the optimal backfill height in multi-level goaf areas