Input parameters, such as rock mass strength parameters and deformation modulus, considered in the design of underground openings involve some uncertainty. The current uncertainty in these parameters is due to the inherent variability of these parameters. To quantify these parameters and design underground openings, the statistical methods must be utilized. In this research, a statistical method was used to define the GSI of rock mass (Geological Strength Index), block volume (Vb), and joint conditions (Jc). Using the GSI distribution function obtained from field data and intact rock strength characteristics, the statistical distribution functions of rock mass parameters were defined using the Monte Carlo method. The statistical analysis of the stability in Azad-pumped storage powerhouse cavern was carried out through the point estimate method. The appropriate support system was suggested according to … Journal Papers