Based on the catalog data of earthquakes with Ms ≥ 2.5 in the Longmenshan fault zone from January 2012 to September 2021, we establish an earthquake time interval series grouped by earthquake magnitude and then use the SARIMA model to predict the series in different periods. By analyzing the fitting effect of the models, the optimal model parameters of different magnitude sequences and the corresponding period values are obtained. Among them, the adjusted R2 values of each model with Ms ≥ 2.5 and Ms ≥ 3.0 sequences are more than 0.86, up to 0.911; the short-time prediction effects are good, and the values of predicted RMSE are 10.686 and 8.800. The prediction results of the models show that the overall trend of the subsequent earthquake time interval in the Longmenshan fault zone is stable, and the prediction results of the Ms ≥ 3.0 sequence have a weak fluctuating growth trend; that is, the number of earthquakes with the Ms ≥ 3.0 in this area will decrease slightly, and the seismicity will decrease in a period of time. The analysis results and method can provide a scientific basis for earthquake risk management and a feasible way to predict earthquake occurrence times.