针对复杂采空区上建设用地问题,将概率积分法和FLAC3D数值模拟相结合,以概率积分法计算所得地表残余沉降量为已知结果,综合岩石试验结果和工程类比经验,建立多煤层复杂采空区的地质模型,进行反复试算、调整,得出研究区的地表移动变形和垂向应力。同时提出了一种基于FLAC3D的采空区上建筑物稳定性计算分析的方法和程序,认为三维数值模拟能较好地解决采空区上建筑物的稳定性问题,弥补了概率积分法不能反映地质体应力的缺点。