Abstract: Accurate prediction and controlling of the tunnel deformation is the key point in ensuring the tunnel con? struction safety. Aiming at the insufficiency of current time series prediction method of tunnel surrounding rock deformation, a time series prediction method of the tunnel deformation based on multivariable Gauss Process (GP) -Differential Evolution Algorithm (DE) is proposed. According to the results of tunnel automation monitoring, the multivariable phase space is reconstructed, and the input dimensions are reduced by principal component analysis. On this basis, the GP-DE model is used to predict the tunnel deformation. Taking the Gaoligou tunnel in Jilin province as an example, the surrounding rock displacement on the vault crown is predicted, and the prediction results are compared with that of BP neural network and SVM model. The results show that the GP-DE model of multi-variable time series has higher prediction accuracy, and the predicted value is in good agreement with the measured value,proving that it is an effective method for tunnel displacement prediction.