Abstract: Heavy rainfall is an important weather that causes flood disasters, and it is also one of the most important natural disasters in our country. With development of the high-resolution numerical models, numerical weather prediction has been the main method for heavy rainfall forecasting. However, the accuracy of the numerical prediction depends on the completeness of the atmospheric motion equations, the accuracy of the initial state, the reasonability of the physical process, and the robustness of the calculation method. As the atmosphere is a nonlinear chaotic system, the small errors in these aspects will cause significant uncertainty in the forecast results. Therefore, the improvement of the rainstorm numerical prediction is closely related to the development of data assimilation, parameterization of physical processes, and ensemble prediction, especially for the role of parameterization schemes of cloud microphysical processes that produce precipitation in numerical models. In addition, in order to improve and perfect the numerical model, the investigation of the evaluation method of the forecast results can also not be ignored and is a crucial part of the numerical prediction. This review describes the development of several key numerical weather prediction techniques. The four-dimensional ensemble variational assimilation method, the microphysics parameterization scheme, and the stochastic kinetic energy backscatter method of ensemble prediction model perturbation are highlighted. An evaluation method of model results based on kinetic energy spectrum analysis is also proposed. Finally, the future research directions in these aspects are summarized.