The spatial distribution of design parameters greatly affects tunnel behavior during and after construction, as well as in the long-term temporal responses. However, the tunnel design parameters commonly used in numerical modeling tend to be representative or average values of global-scale properties. Furthermore, the uncertainty and spatial variation of the design parameters increase as the tunnel scale increases. Consequently, the probability of failure also increases. In order to achieve structural stability in large-section tunnels, the design framework must take into consideration the quantitative effect of design parameter variations on tunnel behavior. Therefore, this paper suggests a statistical approach to numerical modeling to explore the effect of spatially distributed design parameters in a circular tunnel. Also, the effect of spatial variation in the lining strength is studied in this paper. The numerical results suggest that the deformation around the tunnel increases with an increase in the variation of the design parameters.