An emerging kind of prediction model for sequential data with multiple time series is proposed. Because D-vine copula provides more flexibility in dependence modeling, accounting for conditional dependence, asymmetries, and tail dependence, it is employed to describe sequential dependence between variables in the sample data. A D-vine model with the form of a time window is created to fit the correlation of variables well. To describe the randomness dynamically, Bayesian theory is also applied. As an application, a detailed modeling of prediction of abnormal events in a chemical process is given. Statistics (e.g., mean, variance, skewness, kurtosis, confidence interval, etc.) of the posterior predictive distribution are obtained by Markov chain Monte Carlo simulation. It is shown that the model created in this paper achieves a prediction performance better than that of some other system identification methods, e.g, autoregressive moving average model and back propagation neural network.