Geomechanical characterization of a heterogenous rock mass using geological and laboratory test results: a case study of the Niobec Mine, Quebec (Canada)

Abstract To conduct a successful geomechanical characterization of rock masses, an appropriate interpretation of lithological heterogeneity should be attained by considering both the geological and geomechanical data. In order to clarify the reliability and applicability of geological surveys for rock mechanics purposes, a geomechanical characterization study is conducted on the heterogeneous rock mass of Niobec Mine (Quebec, Canada), by considering the characteristics of its various identified lithological units. The results of previous field and laboratory test campaigns were used to quantify the variability associated to intact rock geomechanical parameters for the different present lithological units. The interpretation of geomechanical similarities between the lithological units resulted in determination of three main rock units (carbonatite, syenite, and carbonatite-syenite units). Geomechanical parameters of these rock units and their associated variabilities are utilized for stochastic estimation of geomechanical parameters of the heterogeneous rock mass using the Monte Carlo Simulation method. A comparison is also made between the results of probabilistic and deterministic analyses to highlight the presence of intrinsic variability associated with the heterogeneous rock mass properties. The results indicated that, for the case of Niobec Mine, the carbonatite-syenite rock unit could be considered as a valid representative of the entire rock mass geology since it offers an appropriate geomechanical approximation of all the present lithological units at the mine site, in terms of both the magnitude and dispersion of the strength and deformability parameters. Article Highlights Evaluating the reliability and applicability of geological survey outcomes for rock mechanics purposes. A geomechanical characterization study is conducted on the heterogeneous rock mass by considering the various identified rock lithotypes. The geomechanical parameters of intact units and their associated variabilities are used to stochastically estimate the geomechanical parameters of the heterogeneous rock mass by employing the Monte Carlo Simulation. A comparison is also made between the results of probabilistic and deterministic geomechanical analyses. The results indicate that, in the case of Niobec Mine, the combined syenite-carbonatite rock unit could be considered as a valid representative of the entire rock mass.


Introduction
Site characterization is an essential preliminary phase for implementing a successful rock mechanics program in any underground mining activity. As part of an underground mining plan, site characterization facilitates the subsequent geomechanical classifications by determining the geological settings and lithological characteristics of the area in which the mining activity is taking place. Aside from determining the geological settings and/or hydrogeological characteristics of the mining environment, site characterization contributes in estimation of the strength and deformability parameters of the numerous lithological units identified within the rock mass at the mine site [1,2]. In fact, prior to any geomechanical investigation on rock mass, it is essential to characterize the intact rock in terms of lithological variability, by considering the mechanical properties of the key constituent lithological units [3][4][5].
From a geological perspective, various intact lithological units can be distinguished based on their lithological and petrographic differences. In some cases, numerous units could be identified within the rock having slight differences in mineral assemblage or alteration intensity [5,6]. This approach of dealing with numerous identified lithological units may not be desirable in a geomechanical perspective since it could impose unnecessary complexities in rock mass classification/characterizations or numerical stability analysis. On the other hand, assuming only strong or only weak intact rock quality would either be too conservative or result in overestimation of rock mass strength [5,7]. Therefore, a typical solution to deal with lithological heterogeneity of rock in a geomechanical perspective, could be through defining a periodic presence of strong and weak rock units with their corresponding geomechanical properties.
Lithological heterogeneity should be considered when estimating the geomechanical parameters of rock masses. Variable lithological compositions along with the presence of geological structures would result in anisotropy and heterogeneity in rock mass properties. Heterogeneity and anisotropy in rock mass geomechanical characteristics, originate from the presence of multiple rock formations having horizons with different alternations [2,5]. The existing uncertainty in the strength and deformability parameters of a rock mass stems from the inherent variability associated with the rock formation process as well as a general lack of knowledge throughout the characterization process. The inevitable presence of such uncertainties in intact rock properties, complicates the process of rock mass characterization for a reliable estimate of rock mass properties [6,8].
Different methods for estimating the strength and deformability parameters have been used in geomechanical characterization of rock masses. However, employment of the empirical Hoek-Brown failure criterion in conjunction with the Geological Strength Index (GSI) classification system has been reported to be the most common method for estimation of rock mass properties [9][10][11]. Conventional deterministic application of the above-mentioned approach, is not capable of addressing the intrinsic variable nature of rock masses since the inherent variability of the parameters are not taken into account. Probabilistic estimation of rock mass properties on the other hand, would be capable of incorporating the inherent variability of geomechanical parameters and depict a more realistic picture of the rock mass geomechanical behavior [6,[11][12][13][14][15].
The application of probabilistic methods in rock mass geomechanical characterization have been studied in many researches at different domains of geotechnical engineering such as rock slope stability e.g. [16][17][18][19] and underground rock mechanics such as tunnel stability [20][21][22][23], stope stability [24][25][26] and pillar stability [27][28][29]. However, few studies used probabilistic approaches to quantify the impact of lithological heterogeneity identified by geological surveys and laboratory tests on the uncertainties of the rock mass geomechanical characterization. As one of the pioneering studies, Kim and Gao [30] utilized the Monte Carlo Simulation (MCS) method to estimate the geomechanical parameters of a heterogeneous rock mass and demonstrated that the best-fitted statistical distribution for the mechanical properties of a basalt would be the third type asymptotic distribution of the smallest values. Hoek et al. [31] focused on characterization and determination of engineering properties of tectonically undisturbed sedimentary rock deposits (such as molassic formations) and recommended the use of empirical charts to estimate the GSI and m i values of heterogeneous rock masses. Sari et al. [6], employed the MCS method, to incorporate the uncertainties associated with the intact rock strength and discontinuity parameters into a spread sheet model using statistical distributions, to quantify the variability of rock mass properties for Ankara andesite (Turkey). Tziallas et al. [4] determined the mechanical properties of heterogeneous rocks in the laboratory and proposed a methodology for predicting the rock mass strength of flysch formations consisting of siltstone-sandstone alternations in different proportions. Pepe et al. [5], used GSI classification system in conjunction with the Hoek-Brown failure criterion for geomechanical characterization of the highly heterogeneous Sanremo flysch formation located in Western Italy. Grenon et al. [32] evaluated two successive field and laboratory campaigns to investigate intact rock properties of an underground mining project in the Canadian arctic. The obtained results were compared to the pre-defined target levels of confidence associated with different stages of a mining project. More recently, Contreras et al. [33] applied the Bayesian methodology to estimate the intact rock strength parameters of the Hoek-Brown failure criterion, through the analysis of data from compression and tension tests; they explained the essential differences between frequentist and Bayesian statistics in quantifying the inevitable uncertainty in experimentally determined rock mechanics parameters.
Beside the accomplished researches, further studies are yet required to clarify how lithological and petrographic heterogeneities of a rock mass should be interpreted to provide a useful input for the geomechanical characterization. Identification of numerous lithological units within a rock mass (as the output of geological site characterizations) will not necessarily provide reliable data for rock mechanics experts since the geomechanical parameters of the identified lithological units may not show considerable differences. In fact, the best practice to deal with numerous identified lithological units within rocks should be developed and characterized for rock mechanics experts, not only to increase the efficiency but also to decrease the confusion caused by the presence of dozens of different lithological units, in geomechanical characterization process of heterogeneous rock masses.
To this end, this study uses the results of geological surveys to conduct a geomechanical characterization of rock, in order to evaluate the reliability and applicability of geological survey results for rock mechanics and rock engineering purposes. Geomechanical characterization is conducted on the intact rock specimens collected from the Niobec Mine (Quebec, Canada). Lithological and geomechanical properties of each identified lithological units at the mine, would be brought together from the results of previously conducted field and laboratory test campaigns. lithological similarities between the various constituent lithological units along with the expert judgment would be considered to identify main representative rock units at the Niobec Mine. Subsequently, the geomechanical parameters of each rock unit would be determined and their associated variabilities would be quantified and used to estimate the geomechanical parameters of the entire heterogeneous rock mass. For this purpose, MCS method will be employed to provide a stochastic estimation of strength and deformability parameters of rock mass using Hoek-Brown failure criterion in conjunction with GSI classification method. A comparison will also be made between the results of probabilistic and deterministic analyses to highlight the presence of intrinsic variabilities associated to the heterogeneous rock mass properties. Finally, a valid representative rock unit will be chosen over the others providing the most reasonable and reliable approximation of the geomechanical characteristics for the entire rock mass at the Niobec Mine.

Niobec Mine
The Niobec underground mine as part of the Saint-Honoré alkaline complex is located 13 km northwest of the city of Saguenay (Chicoutimi) within the limits of the municipality of Saint-Honoré, Quebec, Canada (Fig. 1).
The local geology of the region consists of an elliptical-shaped late Precambrian/early Cambrian intrusion core of alkaline silicates (syenite) and dolomitic/calcitic carbonatites [34]. The intrusion is overlain by an average of 75 m of Paleozoic limestone and dolomite of Trenton group [35]. Carbonatite in the intrusion core may have a fine to coarse grain dolomitic or calcitic composition. Carbonatite is brecciated or foliated in texture and contains variable amounts of accessory minerals (such as biotite, magnetite, apatite, pyrite, chlorite and pyrochlore). The foliation is often defined by the alignment of accessory minerals which are constituted from 1 to 10% of the total mineral assemblage [34]. The brecciated carbonatite contains 15% to 90% of syenite fragments from a few centimeters to approximately 10 m in size, showing total or partial alteration to chlorite. Furthermore, the hematitic alteration of carbonatite is also observed with variable intensities. The intrusion core is surrounded by massive or brecciated syenite having different degrees of alteration to chlorite [35]. The zones of economic mineralization in the deposit correspond to irregular shaped sub-vertical lenses oriented towards east-west direction with variable thickness of 73-171 m [34].

Geological description of the rock mass
As mentioned above, the core of the intrusion at the Saint-Honoré alkaline complex is formed by syenite and carbonatite and the whole deposit is surrounded by syenite [34]. By means of the thin section analysis and macroscopic images, the identified lithological units at the Niobec Mine were grouped into three main rock units based on their mineral compositions and content of syenite fragments namely carbonatite, syenite and brecciated carbonatitesyenite. Carbonatite which is competent in term of quality, can contain various accessory minerals such as magnetite, biotite, apatite, chlorite and pyrite whose concentrations vary from 1 to 10%. These minerals are often present as foliation in carbonatite without affecting the resistance of the unit. Carbonatite can be massive to brecciated and can be altered to hematite in different degrees. Carbonatite is brecciated when it contains from 15 to 90% of syenite fragments. Partial or total alteration to chlorite can occur, resulting in resistance reduction in the brecciated carbonatite. The resistance of syenite varies depending on its degree of alteration [34,35].

Research Article
SN Applied Sciences (2021) 3:640 | https://doi.org/10.1007/s42452-021-04617-1 A lithological classification system has been developed based on the rock's composition, grain size and level of accessory minerals to standardize the information regarding the geological description of drill core samples [37]. The various rocks forming the alkaline complex are codified according to their main mineral compositions and categorized into different lithological units. Accordingly, the first category was distinguished by dominance of carbonatite in its composition, while the second category was distinguished by the dominance of syenite. Furthermore, the third category included lithologies associated with particular cases or any other lithological units excluded from the first two categories. The codes which were used to designate the different lithological units, start with the letter "S" if the rock contains between 50 and 95% of syenite fragments and with the letter "C" if contains less than 5% of syenite fragments. Moreover, if the rock contains between 5 and 49% of syenite in its composition, the suffix "S" was added to the type of carbonatite. The rock was designated by the code "Si" or "Sa" if it contains more than 95% of syenite and respectively when not altered and altered in chlorite. Table 1 presents the geological codes established and used by the Niobec Mine authorities to describe the different lithological units identified within the rock mass. A detailed description of the codes attributed to each lithological unit utilized in Table 1, is provided by Lajoie [35].

The geomechanical parameters of the intact rock
Over the past years, various mapping and drilling campaigns have conducted laboratory tests on core samples extracted from different locations within the rock mass at the Niobec Mine [35,[38][39][40][41][42][43][44] to determine the physical and mechanical properties of various existing intact rock units. Table 2 presents the type and number of laboratory tests conducted by different campaigns at the Niobec Mine. Detailed information on the methodologies and specifications of each measuring campaign is provided by Lajoie [35]. The laboratory test results for all the campaigns were compiled and gathered into a database to facilitate the analysis of data. The generated dataset was also revised by Itasca [45] as part of a feasibility study conducted on the mine operation to assess the possibility of applying block caving method. The results were firstly grouped based on similarities of the identified lithological units and secondly based on the depth at which the samples were taken. The analysis of the results indicated that the variation of geomechanical properties, is mainly caused by the lithological differences rather than the depth. In fact, review of the laboratory test results demonstrated that the intact geomechanical properties show large variations in the average values calculated for different  (Tables 3  and 4). Considering the intrinsic nature of the rock materials, there is an uncertainty associated with each of the intact rock parameters, especially when compared to manmade materials. This uncertainty is generated through the process of formation and continuous modification of rock over its geologic history and results in both micro and macro scale property variations from one spatial location to another [8]. The variability associated with the intact  Inherent variability in uniaxial compressive strength (σ c ) and tensile strength (σ t ) stems from the variability of index properties and petrographic characteristics of rock. According to Langford and Diederichs [8], index properties (e.g., the total porosity, density, water content and durability) should be consistent for the rock specimens taken from the same formation (lithological context). While, variations in petrographic characteristics develop different patterns of micro cracks within each rock sample, resulting in variability of uniaxial compressive and tensile strength values. Variability in petrographic characteristics is defined by the variation in grain size and shape, texture (degree of interlocking), micro fractures (grain boundaries, mineral cleavages), nature of cement, degree of chemical alteration and anisotropy (orientation of the cracks and the condition of stress distribution). Similarly, variability in the deformability parameters such as Young's Modulus and the Poisson's ratio is generated due to the variation in parameters such as porosity and degree of jointing, water content and the vibration effects from blasting [8].
In order to quantify the existing variability associated to the intact strength and deformability parameters at the Niobec Mine, lithological units were examined in different categories. Firstly, the strength and deformability parameters of the carbonatite units containing few accessory minerals (C5, C3N and CCA) were compared with those of the foliated units and containing many accessory minerals (C3C, C3A, C3NA, C3B and C3D). Secondly, the strength and deformability parameters of calcific units (C3C, C3D and CCA), were compared to those of dolomitic (C3A, C3B, C5, C3N and C3NA) and syenitic units. Finally, third comparison was made between the strengths and deformability parameters of northern units (C3N, C3NA and C3NT) and all other units [37]. The average values of Young's modulus and Poisson's ratio as well as the uniaxial compressive and tensile strength for each lithological unit obtained from the previous test campaigns are presented respectively in Tables 3 and 4.
Studying the obtained results specified that the largest variations belong to the carbonatitic lithological units. These variations in strength and deformability parameters showed a significant dependence on the degree of alteration.
Box-plot diagrams (Figs. 4 and 5) illustrate the dispersion of strength and deformability test results. Box-plot diagrams are able to compare the series of obtained results for the main lithological units and identify the outliers to be excluded for further calculations. It was indicated that the syenite has lower quality than carbonatite in terms of geomechanical properties. The average values of strength and deformability parameters of carbonatite and   The parameters of Mohr-Coulomb and Hoek-Brown failure criteria for the intact rock were also determined from the results of the uniaxial, and triaxial compression tests carried out by the previous campaigns. The experimental results were first grouped together in order to determine the parameters such as c, φ, m i and σ ci for each of the present lithological units (Table 7). Accordingly, linear regression lines were fitted to the graphs generated from the results of the laboratory tests. The failure envelope curves were subsequently plotted for each of the identified lithological units.
The intact rock failure parameters such as the material constant m i are the parameters mostly defined by the rock type. Variability of such parameters is normally caused by variability in mineral content (lithological heterogeneity), the degree of "particle interlocking", foliation and grain size (texture) [8,46]. According to Table 6, the intact failure parameter values of different carbonatite units, are relatively constant except for the units either with high alteration degrees or with high presence of accessory minerals. Moreover, based on the results of triaxial tests, intact specimens with a more uniform lithological context and  less concentration of micro cracks, showed higher strength than the specimens with high concentration of accessory minerals, or with a large amount of syenite fragments altered in chlorite. The test results are then combined in order to compare the failure parameters of the carbonatite and syenite units. Figure 6 presents the obtained Mohr-Coulomb and Hoek-Brown failure envelopes for the groups of carbonatites and syenite rock units while Fig. 7 illustrates the obtained Mohr-Coulomb and Hoek-Brown failure envelopes for the carbonatite-syenite rock unit. The average values of the failure parameters for the three mentioned rock units are reported in Table 8.
According to the presented results in Table 7, the estimated values of failure parameters for carbonatite is slightly higher than syenite, but fairly close to the values obtained for the carbonatite-syenite rock unit.
Even though a slight difference was observed between the calculated strength and deformability parameters of carbonatite and syenite intact units, the carbonatitesyenite provided a quit fair approximation of the two units. Therefore, assuming a homogeneous intact rock containing a random combination of syenite and carbonatite units could provide an agreeable geomechanical description of intact rock for the case of Niobec Mine. Table 9 summarizes the average values to be considered in characterizing the properties of the intact rock at the Niobec mine.

Geomechanical characterization of the rock mass
It was shown in the previous section that in the scale of intact rock, a homogenous carbonatite-syenite rock rather than considering only a carbonatitic or syenitic intact rock provides a more realistic approximation for the existing rock at the Niobec Mine. This assumption is aligned with the findings of local site characterization programs conducted during the deposit exploitation which emphasized that the nature of rock mass at the Niobec Mine is too complex to be discretized geomechanically hence the observed irregularities in the quality of intact rock units should be considered to occur in local scales [34,35]. In this part, the same hypothesis is examined in the rock mass scale to see whether the results would support the consideration of a unified carbonatite-syenite rock mass unit for the Niobec Mine or not. Geomechanical characterization of rock mass, can be accomplished through using the obtained results

Probabilistic estimation of rock mass geomechanical parameters
Hoek-Brown failure criterion is used in conjunction with the Geological Strength Index (GSI) to estimate the geomechanical parameters of the syenitic, carbonatitic  Table 10). The mean GSI value and the corresponded standard deviation of the rock mass were adopted from the study conducted by Lavoie [37]. Based on the concentration of discontinuities and the number of joint sets identified within the carbonatite-syenite unit, the rock mass was recognized to have a blocky structure with rough and slightly altered joint surface condition. Therefore, to assign an appropriate value of GSI to the rock mass at the Niobec Mine, the calculated value of RMR 89 as equal to 71 was used (as reported by Golder [34] and Lajoie [35]) and the GSI was calculated as 74 based on the conversion formula provided by Hoek and Brown [3]. It should be noted that estimating the GSI value from RMR 89 might be unreliable, especially for poor quality rock masses and for rocks with lithological peculiarities that cannot be easily incorporated in the RMR classification [9]. However, in this case, due to the high quality of rock mass and lack of GSI data resulted from direct GSI estimations, the above-mentioned conversion formula is used.
Since the reported RMR 89 value was estimated for the entire rock mass and no GSI value was available separately for different units within the rock mass, the same GSI value of 74 was considered for each syenite and carbonatite units for further geomechanical calculations (Table 10). These values are considered as random inputs when calculating the geomechanical parameters.
The Kolmogorov -Smirnov statistical test was applied on the above-mentioned parameters to assign the best-fitted probability distributions function (PDF) to each parameter. According to the results of Akaike Information Criterion (AIC) goodness of fit test, Weibull and Normal distributions were determined to be the most appropriate statistical distributions for the parameters. However, for the sake of simplicity, normal PDFs were chosen to be assigned to the parameters. Assuming normal distribution function to describe the random characteristics of the geomechanical parameters is consistent with the suggestions made by many studies e.g. [6, 20, 22, 26, and 47].
To eliminate generation of negative and/or false values, the normal PDFs were truncated by using the actual reported minimum and maximum values. Figures 8 a-c illustrate examples of the truncated PDF plots for intact rock material properties of the carbonatite-syenite rock mass.
Subsequently, the MCS method was applied to calculate the Hoek-Brown strength and deformability parameters of each rock mass type through Eqs. 1-6 [9].
where D represents the degree of disturbance of the rock mass (ranging from 0 to 1 for undisturbed in-situ rock masses to highly disturbed rock masses) [9]. Since the rock 3 )  mass is assumed to be undisturbed by controlled blasting, the parameter D is considered equal to 0 for this study. Each output parameter was generated from 10,000 iterations of randomly selected combinations of input parameters in accordance with their assigned PDF using Latin hypercube sampling algorithm (LHS). The advantage of LHS is that it provides smoother resulting PDFs with fewer iterations by using stratified sampling models. For the sake of simplicity, all the input parameters were assumed as independent variables. Although the dependence of the Hoek-Brown parameters could jeopardize the probabilistic estimation of the output parameters, Sari et al., [48] stated that the dependence between the Hoek-Brown parameters does not have a significant effect on the estimation of strength and deformability parameters using probabilistic simulations.
By the aim of MCS, the variability associated to the output geomechanical parameters were quantified and the Kolmogorov-Smirnov test determined the best-fitted distribution function for each parameter. The mean, standard deviation and best-fitted PDF of strength and deformability parameters for each rock mass type are reported in Table 11. Figures 9, 10, 11a-f respectively illustrate the obtained truncated PDF plots of output geomechanical parameters for syenite, carbonatite and carbonatite-syenite rock mas ses.
The results of rock mass geomechanical characterization indicated that, similar to the intact properties, the obtained values for syenite and carbonatite rock masses show slight differences and the carbonatite rock mass quality was estimated to be higher than syenite. However, estimated geomechanical parameters for the carbonatitesyenite rock mass were fairly close to the values calculated only for the carbonatite.

Deterministic estimation of rock mass geomechanical parameters
The RocData 5.0 software [49], was used to provide a deterministic estimation of the carbonatite-syenite rock mass parameters through the calculated mean values of the input parameters (Table 10). Similarly, Eqs. 1-6 were used by the software and the factor D was considered equal to 0 [50]. The failure envelopes of the Hoek-Brown criterion, traced for the calculated parameters, are presented in Fig. 12. Table 12 presented the average estimated values of strength and deformability parameters for the carbonatite-syenite rock mass.

Discussion on the geomechanical characterization results
Studying the results of geological and geomechanical characterizations, emphasizes the necessity of adopting an accurate interpretation of data considering both lithological and geomechanical descriptions of a rock mass. In fact, lithological heterogeneity shouldn't provide a misleading indication on the presence of geomechanical variety throughout that rock mass. In this regard, studying the field and laboratory tests results of different test campaigns at the Niobec Mine, revealed that lithological heterogeneity rather than depth was responsible for the variation in geomechanical properties of the three main identified rock units. Studying the results proved that the largest variations in geomechanical parameters belong to the carbonatitic lithological units either with maximum separation from the syenitic lithological units or with intense alteration to chlorite. In fact, the resulted variability in strength and deformability parameters of intact samples was determined to be significantly defined by the degree of carbonatite alteration. Furthermore, geomechanical characterization of intact samples identified the carbonatitic units with a better quality than syenitic units. Determination of the Mohr-Coulomb and Hoek-Brown failure parameters for each constituent lithological unit and subsequently each rock unit (Tables 7 and 8), revealed that the strength of different carbonatitic lithological units was relatively constant except for the units either with intense alteration or with high presence of accessory minerals. In fact, according to the results of previously conducted triaxial tests, it was concluded that high concentrations of accessory minerals within the intact samples, presence of a contact between two different lithological units and presence of large quantities of altered syenite fragments to chlorite could result in reduction of rock strength compared to equivalent samples with more uniform lithological context and less concentration of micro cracks.
However, in geomechanical perspective, the observed difference between the calculated strength and deformability parameters of carbonatite and syenite units, wouldn't be sufficient enough to justify the assumption of considering them as distinguished units within the intact rock; however, the carbonatite-syenite rock unit provided a quite reasonable approximation of geomechanical properties of the two aforementioned units. Besides, findings of local site characterization programs during the deposit exploitation phases, demonstrated that the nature of rock mass at the Niobec Mine is too complex to be discretized geomechanically and the captured irregularities in the quality of intact rock units are considered to occur in local scales [34,35]. Therefore, assuming a homogeneous intact rock containing a random combination of syenite and carbonatite rock units could provide a valid geomechanical description of intact rock for the case of Niobec Mine.
In the rock mass scale, estimation of geomechanical parameters separately for syenite, carbonatite and carbonatite-syenite rock units also determined the quality of carbonatite rock mass to be slightly higher than syenite but fairly close to the quality of the carbonatite-syenite rock mass. Even the dispersions of the estimated strength and deformability parameters, were obtained to be very similar between carbonatite and carbonatite-syenite rock masses (Table 11). The obtained results of rock mass  parameters for the particular case of carbonatitesyenite rock mass, proved the significant presence of variability associated to each parameter and inability of conventional deterministic approaches to address them entirely (Table 12). Conventional deterministic approaches cannot capture the variability of rock mass properties since they assign a single mean value to the parameters instead of defining them through probability distribution functions (which consider the mean and standard deviation values). This simplification not only ignores the inherent variable nature of rock mass, but also by assuming a better rock mass quality in terms of strength and stiffness produces unreliable and unclear results that could be misleading in subsequent rock mechanics analyses. It is hence recommended to incorporate uncertainties associated with the variability of rock mass parameters into consideration to ensure that the actual behavior of rock mass is properly reflected.

Conclusion
The results of previously conducted geomechanical and geological field and laboratory tests at the Niobec Mine (Quebec, Canada) were combined to characterize the heterogeneous rock mass by considering the lithological and mechanical properties of identified lithological units' constituents. The aim of this study was to find a reasonable agreement between the geomechanical parameters in relation to the extensive lithological variability for describing the rock mass properties. The results of intact rock characterization indicated that the carbonatite-syenite rock unit could be considered as an appropriate representative lithology to define the rock mass geomechanical properties. Furthermore, estimated rock mass geomechanical parameters for syenitic, carbonatitic and carbonatitic-syenitic units also indicated that considering the carbonatite-syenite rock mass instead of trying to distinguish the syenite and the carbonatite as separate units provide a reasonable and  reliable approximation of the geomechanical quality of rock mass at the Niobec Mine. Moreover, consideration of the carbonatite-syenite rock unit to represent the entire rock mass provides a good agreement between both the geomechanical and geological perspectives. Finally, the use probabilistic approaches instead of conventional deterministic methods in rock mass geomechanical characterization programs are highly recommended since a more realistic portray of the intrinsic nature of rock materials is depicted by considering the inherent variability associated to the geomechanical parameters.
Acknowledgements The authors would like to acknowledge the funding received by a grant from Natural Sciences and Engineering Research Council of Canada (NSERC) (CRDPJ 469057-14) for this study. We also thank the personnel at the Niobec Mine for their assistance.

Fundings
The authors have no relevant financial or non-financial interests to disclose This study was funded by a grant from Natural Sciences and Engineering Research Council of Canada (NSERC) (CRDPJ 469,057-14).

Conflict of interest
The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.