In this study, the probability integral method, Synthetic Aperture Radar Interferometry (InSAR), and the Okada dislocation model were collaboratively used to analyze deformation in the Alhada Pb-Zn mine. The predicted deformation values of the subsidence centers in three subsidence areas were 107 mm, 120 mm, and 83 mm, respectively, as predicted using the probability integral method. The coherent scatterer InSAR technique was used to analyze the time-series deformation of the mining area, and the same subsidence center locations and similar deformation values were observed. The Okada dislocation model was used to invert the optimal parameters of the underground-mining ore body causing the surface subsidence, indicating that the surface subsidence is mainly caused by the mining of ore bodies in the 888 and 848 middle sections. We further simulated ground deformation using the multi-source Okada model. The results showed that the predicted and modeled deformations are highly correlated with the observed deformation. Through the analysis and comparison of the InSAR results, it was concluded that the three subsidence areas do not threaten the stability of the main buildings in the mining area. Using theoretical, observational, and modeling methods, the development and evolution of the subsidence area in mines can be established, which could provide basic data for subsidence control work and guarantee mine production safety.