Abstract: In order to reliably predict the solar wind in the near earth environment, the deep learning multi-step prediction method is adopted to predict the solar wind speed at the Lagrange point 1 (L1) between the sun and the Earth in the next 24, 48, 72, and 96 hours. The SDO image data are used to extract the area of the coronal hole and other input features from the NASA OMNIWEB data set to form the multivariable time series as the input information for the solar wind speed prediction. In the multi-step prediction method, the output is no longer the solar wind speed at a certain moment in the future, but the solar wind speed for several hours around that moment at the same time. Three deep learning prediction models are used respectively: the single LSTM model, the encoder-decoder LSTM model and the CNN-LSTM encoder-decoder model as the base model, and the multi-step output data are included for the comparative analysis. Experiments show that, compared with the single-step output prediction, the multi-step prediction improves the correlation coefficient and the forecast accuracy of the solar wind speed forecasts in the next 24, 48, 72, and 96 hours. For the 24-hour prediction, the best experimental results are achieved: the correlation coefficient with the observed data is 0.789, the RMSE is 62.469 km/s, and the MAE is 46.476 km/s.