Abstract: Due to the characteristics of propofol of high time-varying, and complex compartment model, the traditional method of nonlinear mixed effects modeling (NONMEM) has miscellaneous of variables and plenty of artificial factors in the estimation of propofol. This study was aimed to build a propofol prediction model based on the differential evolution (DE) algorithm and grey model. DE was used to optimize the pa-rameter of multi-variable grey model (MGM) and to build a model of prediction of the plasma concentration of propofol based on the grey model. It was compared with the results of NONMEM algorithm. In conclusion, the median performance error (MDPE) of DE-MGM was -4.6%, while the result of NONMEM is -12.13%. The median absolute performance error (MDAPE) of GA-BP neural network is 13.19%, while that of NONMEM is 23.12%. The experimental results suggest that the new method is suitable to determine the short half-life of anesthesia drug propofol with higher accuracy.