In this paper, we make long-term predictions based on numbers of current confirmed cases, accumulative dead cases of COVID-19 in different regions in China by modeling approach. Firstly, we use the SIRD epidemic model (S-Susceptible, I-Infected, R-Recovered, D-Dead) which is a non-autonomous dynamic system with incubation time delay to study the evolution of the COVID-19 in Wuhan City, Hubei Province and China Mainland. According to the data in the early stage issued by the National Health Commission of China, we can accurately estimate the parameters of the model, and then accurately predict the evolution of the COVID-19 there. From the analysis of the issued data, we find that the cure rates in Wuhan City, Hubei Province and China Mainland are the approximately linear increasing functions of time t and their death rates are the piecewisely decreasing functions. These can be estimated by finite difference method. Secondly, we use the delayed SIRD epidemic model to study the evolution of the COVID-19 in the Hubei Province outside Wuhan City. We find that its cure rate is an approximately linear increasing function and its death rate is nearly a constant. Thirdly, we use the delayed SIR epidemic model (S-Susceptible, I-Infected, R-Removed) to predict those of Beijing, Shanghai, Zhejiang and Anhui Provinces. We find that their cure rates are the approximately linear increasing functions and their death rates are the small constants. The results indicate that it is possible to make accurate long-term predictions for numbers of current confirmed, accumulative dead cases of COVID-19 by modeling. In this paper the results indicate we can accurately obtain and predict the turning points, the end time and the maximum numbers of the current infected and dead cases of the COVID-19 in China. In spite of our simple method and small data, it is rather effective in the long-term prediction of the COVID-19.