When limit-state functions are highly nonlinear, traditional reliability methods, such as the first-order and second-order reliability methods, are not accurate. Monte Carlo simulation (MCS), on the other hand, is accurate if a sufficient sample size is used but is computationally intensive. This research proposes a new system reliability method that combines MCS and the Kriging method with improved accuracy and efficiency. Accurate surrogate models are created for limit-state functions with minimal variance in the estimate of the system reliability, thereby producing high accuracy for the system reliability prediction. Instead of employing global optimization, this method uses MCS samples from which training points for the surrogate models are selected. By considering the autocorrelation of a surrogate model, this method captures the more accurate contribution of each MCS sample to the uncertainty in the estimate of the serial system reliability and therefore chooses training points efficiently. Good accuracy and efficiency are demonstrated by four examples.