禹州郭村煤矿采空区稳定性问题,是南水北调中线一期工程重大工程技术问题。基于概率积分法沉陷理论,运用采空区在自重下残余变形极限值计算公式,以郭村煤矿采空区剖面4为基础,应用FLAC3 D软件对采空区无注浆处理时渠道开挖前和施工运行中两种工况建立数值模型,研究其变形规律。计算出郭村煤矿采空区竖直残余变形极限值为331 mm ,水平残余变形极限值为82.8 mm。研究表明:采空区上方两端地表处的水平位移值最大,其值受渠道开挖影响很小;地表竖直位移总体呈随采空区埋深增大而减小,随离渠道水平距离增大而减小的趋势;地表变形受渠道的影响相对较大,这将为采空区地面开挖等提供理论参考。