The Xishan coal mine area in Beijing, China has a long history of mining. Many landslide hazards, in addition to collapses and ground fractures, have occurred in this area. This study used multi-temporal satellite images to extract this region’s deformation information, identify landslides and analyze the deformation evolution process of these landslides. Taking the Anzigou ditch as an example, we investigate the “Quarry–Landslide–Mudslide” disaster chain model. We found that the landslide evolution process is closely related to the geological conditions, and usually goes through four stages: initial deformation, slope front swelling and collapsing, rear part connecting and rupturing, and landslide creeping. The surface deformation can be identified and tracked by high-resolution optical images and InSAR monitoring. Under the combined effects of rainfall and topographic conditions, medium and large landslides may occur and trigger a “Quarry–Landslide–Mudflow” disaster chain. The identification and analysis of these landslide hazards and the disaster chain help with geological disaster prevention, and provide reference for early identification and research of similar disasters.

Ground-based synthetic aperture radar interferometry (GB-InSAR) has the characteristics of high precision, high temporal resolution, and high spatial resolution, and is widely used in highwall deformation monitoring. The traditional GB-InSAR real-time processing method is to process the whole data set or group in time sequence. This type of method takes up a lot of computer memory, has low efficiency, cannot meet the timeliness of slope monitoring, and cannot perform deformation prediction and disaster warning forecasting. In response to this problem, this paper proposes a GB-InSAR time series processing method based on the LSTM (long short-term memory) model. First, according to the early monitoring data of GBSAR equipment, the time series InSAR method (PS-InSAR, SBAS, etc.) is used to obtain the initial deformation information. According to the deformation calculated in the previous stage and the atmospheric environmental parameters monitored, the LSTM model is used to predict the deformation and atmospheric delay at the next time. The phase is removed from the interference phase, and finally the residual phase is unwrapped using the spatial domain unwrapping algorithm to solve the residual deformation. The predicted deformation and the residual deformation are added to obtain the deformation amount at the current moment. This method only needs to process the difference map at the current moment, which greatly saves time series processing time and can realize the prediction of deformation variables. The reliability of the proposed method is verified by ground-based SAR monitoring data of the Guangyuan landslide in Sichuan Province.