Industry-scale hydrogen is mainly produced by steam methane reforming (SMR), which uses natural gas as the feedstock and fuel and co-produces CO₂. This study aims to numerically evaluate hydrogen production by SMR under various reacting conditions. Unlike the previous studies with limited scenarios, the performance of SMR is continuously evaluated in a high-dimensional input-parameter space. The SMR plant including a combustor, a reformer, and a water-gas shifter is modeled in Aspen HYSYS software. The four key parameters, including methane fraction of the feedstock, reformer pressure and temperature, and shifter temperature, are treated uncertain and 50 samples are drawn from a four-dimensional parameter space defined by their ranges. Each sample is input to HYSYS model and mass ratio of each component in product streams is obtained as the output variables. Based on the 50 pairs of input-output data, response surfaces of the outputs are developed to surrogate HYSYS models. The fast response surface models are then used to calculate global sensitivity indices and evaluate SMR processes. Results show the reformer performance is controlled by temperature rather than pressure, and a temperature higher than 900 ${}^{\circ }$C can maximize the reaction rate. The water-gas shifting reaction is inhibited in the reformer but significantly enhanced in the shifter. Hydrogen is mainly produced in the reformer while the major function of the shifter is to convert CO to nontoxic CO₂.