本申请公开一种适用于煤矿采空区覆岩及地表变形的数值模拟方法及装置，其中，方法包括：基于预设模型试验，拟合不同荷载条件下垮落带岩体蠕变参数随荷载的变化方程；构建采空区围岩、覆岩以及地表的三维计算模型，从而模拟工作面的开挖过程，得到工作面开挖后研究区域的重分布应力场；进而读取垮落带范围内每一个计算单元的荷载，并代入垮落带岩体蠕变参数随荷载的变化方程，计算垮落带范围内每一个计算单元的蠕变参数；基于蠕变参数，模拟垮落带岩体在荷载作用下的蠕变过程，得到煤矿采空区覆岩及地表长期变形模拟结果。由此，解决了相关技术中未能考虑垮落带不同区域的岩体蠕变特性的差异，降低煤矿采空区长期变形数值模拟结果的准确性等问题。