Dynamics model of the rocket and its debris were established, and the quaternion algorithm was used to compute the carrier attitude. A prediction method of safety control zone based on infill-sampling Kriging model was proposed, which combined the characteristics of Monte Carlo and Kriging surrogate models. Then a flow chart of prediction was provided furthermore, and one booster rocket′s debris safety control zone was simulated. Simulation results show that, compared with the Monte Carlo method without losing computational accuracy, the proposed method has equivalent computational accuracy and higher computational efficiency:it can meet the rapid iteration requirement for engineering design. While compared with the conventional extreme deviation overlying method, it can reduce the area of safety control zone dramatically and has relatively good engineering application value.