Abstract: The working face is divided into some mining elements along strike, and surface movement and deformation value equals to the superposition of the separate dynamic displacement and deformation values of all these mining elements. A model of predicting dynamic mining subsidence for multi-face mining is made according to the Knothe time function and probability integral method. The cell division method and the conversion between different coordinate systems in the model are studied. An algorithm step of the dynamic prediction model for multi-face mining is designed. The feasibility and practicality of the algorithm are confirmed by comparing of the predicted and measured data and showing the prediction effect of multi-face mining surface subsidence of some mine.