The purpose of this study is to evaluate the rock mass stability around the tunnels by three-dimensional (3-D) numerical modeling. Complex geological and tunnel features have been included in the model using the software package 3DECTM. Based on the available information on stratigraphy, geological structures, and mechanical properties of rock masses and discontinuities from the underground mine, 3-D stress analyses were performed on various cases with different K0 (lateral stress ratio) values, constitutive models, as well as support systems. Relations between these factors and the distributions of stress, displacement, and failure zones around the tunnels were obtained and discussed. Finally, comparisons between the field deformation measurements and numerical simulations are used as validations for the simulations, helping to determine the appropriate K0 value and material constitutive models for the underground mine.