Abstract

The shear strength parameters (i.e., the internal friction coefficient and cohesion ) are very important in rock engineering, especially for the stability analysis and reinforcement design of slopes and underground caverns. In this paper, a probabilistic method, Copula-based method, is proposed for estimating the shear strength parameters of rock mass. The optimal Copula functions between rock mass quality and , and for the marbles are established based on the correlation analyses of the results of 12 sets of in situ tests in the exploration adits of Jinping I-Stage Hydropower Station. Although the Copula functions are derived from the in situ tests for the marbles, they can be extended to be applied to other types of rock mass with similar geological and mechanical properties. For another 9 sets of in situ tests as an extensional application, by comparison with the results from Hoek-Brown criterion, the estimated values of and from the Copula-based method achieve better accuracy. Therefore, the proposed Copula-based method is an effective tool in estimating rock strength parameters.

1. Introduction

The shear strength parameters including friction coefficient (the tangent of friction angle ) and cohesion are the most fundamental mechanical parameters for rock mass. They are of very high importance for stability evaluation and support design of rock mass in rock engineering, especially for landslide risk mitigation. Different approaches including laboratory tests, in situ tests, back analysis, and empirical equations can be used to estimate these strength parameters. In the laboratory tests, scale effects are involved since rock specimens are of small sizes [1, 2]. The back analysis method is also recommended by [3, 4] the fact that it might be more reliable than laboratory or in situ test method. Some empirical equations between shear strength parameters and geological indices are widely used due to simplicity and practicality [513]. The empirical equations based on famous Hoek-Brown criterion are [13] where is a reduction of the rock material constant and obtained by where is the Geological Strength Index and is a disturbance factor, depending on the degree of disturbance by which the rock mass has been subjected to blast damage and stress relaxation. The varies from zero for undisturbed in situ rock mass to 1 for very disturbed rock mass. and in (1) can be calculated by

In (1), is uniaxial compressive strength, , and is the upper-limit value of confining stress discussed by Hoek et al. [13].

It is important to take full advantage of in situ test data in order to establish the correlativity between geological indices and mechanical parameters of rock mass. With regard to the correlativity of the above variables, an effective way is to establish their joint distribution function. The common statistical methods for obtaining joint distribution function of variables are generally based on two assumptions; that is, the variables have the same marginal distributions and follow the normal distribution [14]. However, rock mass parameters do not always follow the normal distribution especially as the number of samples is small [15]. With respect to small samples, normal information spread estimation method (NISEM) [15] is suitable. Meanwhile, Copula theory [1621] is a fairly effective method to establish the joint distribution function of nonnormal distribution variables. The advantages of Copula theory are that the marginal distributions and dependence structure of variables can be studied separately and the variables can follow any distributions. The Copula theory has been widely used in many fields such as insurance [16], finance [17], hydrology [1820], and geotechnical/geological engineering [21].

In this paper, the correlations between the rock mass quality indices (i.e., rock quality designation (RQD), rock mass quality (), and Rock Block Index (RBI)) and the shear strength parameters ( and ) are investigated based on the results of in situ direct shear tests at Jinping I-Stage Hydropower Station in China. It shows that presents the strongest correlation with the shear strength parameters. Then, the joint distribution functions of and , and are proposed based on Copula theory. By means of the functions and the given guarantee rates, the shear strength parameters can be estimated. Compared with those estimated by Hoek-Brown criterion, the estimated values using the proposed Copula-based method are closer to the in situ tests ones.

2. Correlations between Rock Mass Quality Indices and Shear Strength

2.1. In Situ Test Data

The Jinping Hydropower Station is located on the middle reach of the Yalong River, in Sichuan, China. It has a double-curvature arch dam with a height of 305 m and a total installed capacity of 3300 MW [22, 23]. The Hydropower Station is located in the east edge of the Tibet Plateau tectonic belt. Due to the rapid crust uplift during the Quaternary, the steep slopes exist on both sides of the valley deeply cut by Yalong River. The stratum in the dam zone consists of a series of metamorphic rocks of middle of Triassic period, of which the great majority are marbles with grey or white and fine-grained texture. The main mineral compositions of the marble are calcite and dolomite. The uniaxial compressive strength of the fresh marble is about 60–75 MPa.

In order to investigate the shear strength parameters of the representative marbles, 12 sets of in situ direct shear tests with five similar specimens in each set were performed in the exploration adits of Jinping I-Stage Hydropower Station. Figure 1 shows the sketch of in situ direct shear test which was performed using flat stacking method [24]. The side length of each cube specimen is 1 m which contains in general two sets of joints (bedding joints and cross joints; see Figure 1). RQD and RBI [25, 26] of each specimen can be obtained using scanline field mapping method. The scanline used for measuring RQD was parallel to shear direction. The shear direction was parallel to the bedding joint along which the shear deformation happened. The five similar specimens in a set were tested at different normal stress from 2 to 10 MPa. The shear strength parameters ( and ) of each set are obtained by linear regression according to Mohr-Coulomb criterion. A good fitting straight line between the peak shear stress and normal stress can be derived from each set test data for which the of fitting degree is greater than 0.9 (see Figure 2). The test results are shown in Table 1, in which the rock mass quality indices are the average of each set samples, respectively.

2.2. Correlation Measure

Pearson linear correlation coefficient [27], Spearman’s rank correlation coefficient [28], and Kendall’s rank correlation coefficient [29] are widely used to measure how relevant the variables are. They can be calculated by where is the number of observations, and are the observed values of variables ( and ), and are sample mean, and are sample variance, and are the rank of and , if and , and if and ; otherwise, , , .

Table 2 shows the results of the , , and of test data. There are clearly positive correlations (the same variation trend) between three rock mass quality indices (RQD, RBI, and ) and . The correlation coefficients between and are much greater than those between RQD and , RBI and . There is clearly negative correlation (changes in the opposite trend) between the three indices and . The correlation between RBI and is slightly better than and , while RQD and have the worst correlation. In summary, has the best correlation with shear strength parameters ( and ). In fact, RQD and RBI evaluate the rock mass quality only by the integrality of rock core. While many other factors, such as joint surface roughness and degree of weathering for dominantly oriented joint sets, joint water reduction, stress reduction factor, and volumetric joint count, are fully involved in -system, the -system, therefore, provides a comprehensive description on rock mass quality [30]. The shear strength parameters of jointed rock mass are generally controlled by complex geological conditions and mechanical behaviors of rock material and so on. The best correlation between and the shear strength parameters in this statistical study also indicates that -system associated with multiple factors of influence on and of jointed rocks more precisely depicts the quality of rock mass.

The correlation measure indicates that, for the in situ tests presented in this paper, the marbles with higher rock mass quality indices show higher friction angle and lower cohesion. Many researchers had also found that, for rocks and soils, there was a negative correlation between and [31, 32]. If this negative correlation is true between and , it is easily understood that a variable has a positive correlation with one of shear strength parameters while a negative correlation must appear with the other based on the pure mathematical logic. As shown in Figure 1, there are some intact rock segments on the sheared bedding joint in the case of relatively low value due to cross joint, which might explain why higher has higher and lower . The shear failure occurs along the weak bedding joint with a relatively low value on condition that the tested specimen has a high value (good rock mass quality). When a low value appears, these intact rock block segments on the shear surface must be thus cut through before reaching peak strength, which can increase the cohesion of the failure surface.

3. Copula-Based Method

3.1. Copula Theory

Considering multivariate random variables with marginal distribution functions as , where is the number of random variables and is the value of random variable (), is the marginal distribution function and is the probability when . The joint distribution function of the multivariate random variables can been expressed as [33]

The Copula function is to connect multivariate probability distributions to their marginal probability distributions. Thus, the multivariate joint distribution is expressed in terms of its marginal distribution function and the associated dependence structure as where , called Copula function, is a uniquely determined mapping whenever are continuous and captures the essential features of the dependence between random variables. is essentially the joint distribution of multivariate random variables of through their marginal distributions [33].

The joint distribution function established using Copula-based method has two major advantages: the marginal distribution and dependence structure can be studied separately, and the variables can follow any distributions. Therefore, the problem of determining reduces to determining marginal distribution of variables and Copula function .

In order to clearly demonstrate the calculation procedure using the present Copula-based method, the calculation flow chart is shown in Figure 3. The detailed procedures are as follows.

Step 1. The correlation between target variables and relevant variables (, RQD, and RBI) is analyzed according to the correlation coefficients obtained by (4) or else. The best relevant variable (i.e., ) can be selected, as shown in Section 2.2.

Step 2. Utilize suitable functional forms and methods (such as normal information spread estimation method (NISEM)) to establish marginal distribution function of variables which could be accepted according to corresponding test (e.g., Kolmogorov-Smirnov (K-S) test). It is addressed in Section 3.2.

Step 3. Introduce suitable multivariable joint distribution functional forms (i.e., Copula functions) and confirm the optimal Copula function indicated by a minimum AIC value. Such is stated in Section 3.3.

Step 4. Estimate the values of target variables using the Copula method on condition that a guarantee rate and the best relevant variable have been given, which is presented in Section 4.

3.2. Verification of Marginal Distribution Functions

For small samples, the normal information spread estimation method (NISEM) [15] is an effective way to generate probability density functions. The main difference between NISEM and the available classical methods [15, 34] is that the variables are not assumed to follow a distribution (such as normal, log normal, and exponential distribution) in advance. It has been proved that the NISEM used for generating the probability density functions of geotechnical parameters with small samples is more accurate and effective than the other classical methods [15]. The probability density function of random variables by using NISEM is expressed as where is the window width of standard normal information spread function and is the observed value of variable . By the approaching principle of normal information diffusion, can be calculated by where and are the maximal and minimum observed values of variable . The value of is related to , as shown in Table 3 [35].

For the variables in Table 1, while . Therefore, the estimation of probability density functions can be obtained from (7). The values of and (parameters of ) are shown in Table 4.

By means of the probability density functions (see (7)), the marginal distribution functions can be established as

Kolmogorov-Smirnov (K-S) test commonly used to check whether a set of data follows a certain distribution [36] is adopted to test the established marginal distribution functions. The procedure of K-S test involves forming the cumulative frequency distribution and computing the test statistic : The cumulative frequency distribution is defined as where is the indicator function, equal to 1 if and equal to 0 otherwise.

The established marginal distribution functions will be accepted if , otherwise they will be denied. Where is the critical value of a given significance level (usually [37]), the can be obtained using the table established by [38] when and have been determined. Here, the significance level is 5% and is 12, so based on the relevant table in [38]. The K-S test results are shown in Table 5. Obviously, all the established marginal distribution functions can be accepted.

3.3. Determination of the Optimal Copula Functions

In order to obtain the optimal Copula functions, the common method of the goodness-of-fit is adopted. Therefore, two steps are taken to determine the optimal Copula function: (1) selecting the suitably qualified Copula functions according to the characteristics of marginal distribution functions and dependent structure of variables and (2) determining the parameters of Copula functions selected.

With respect to two-dimensional Copula functions, the parameter of Copula functions can be obtained by means of Kendall’s rank correlation coefficient , which has been calculated in (4). The relation between and is [39] where and are the marginal distribution functional values on condition that the values of two variables are equal to and (i.e., and ), while is the introduced Copula functional value with the parameter .

Akaike information criteria (AIC) [40] are commonly used for goodness-of-fit test. AIC are based on entropy principle and are a simple and effective method for goodness-of-fit test [41]. And the optimal Copula is the one which has the minimum AIC value; the value of AIC can be calculated as where is the number of parameters of Copula function ( for two-dimensional Copula function) and is the Copula functional value. is the empirical probability functional value; it can be obtained as where is the number of joint observations under the condition of as well as .

In order to visualize the dependence structure of the measured data () and to obtain an insight on a suitable Copula, the measured data in original space should be transformed into the standard uniform random vector . Tang et al. [21] defined the random vector by adopting the empirical distributions of measured data. and are expressed as where (or , ) denotes the rank of (or , ) among (or , ) in an ascending order. are the values of variables () or () corresponding to () or ().

The scatter plots of empirical cumulative distribution of and are shown in Figure 4. and and are substantially symmetric about the diagonal lines. Although the data points ( and ) in right Figure 4(b) seem to be not symmetric as well as those in left Figure 4(a), most of them are located around diagonal line (i.e., symmetrical line) in Figure 4(a). Archimedean Copulas have a wide range of application [1921, 42] due to several competitive advantages: (1) the ease with which they can be constructed; (2) the great variety of families of copulas which belong to this class; and (3) the many nice properties possessed by members of this class [39]. As results, the symmetrical Archimedean Copulas, Nelsen No. 1, Nelsen No. 2, Nelsen No. 4, Nelsen No. 5, Nelsen No. 12, Nelsen No. 14, Nelsen No. 15, and Nelsen No. 18 [39] are selected to model the positive correlation between and , while Nelsen No. 1, Nelsen No. 2, Nelsen No. 5, Nelsen No. 7, Nelsen No. 8, and Nelsen No. 15 [39] are selected to model the negative correlation between and . By means of AIC calculated by (13), the goodness-of-fit test is used to determine the optimal Copula functions. The chosen symmetrical Archimedean Copula functions and the values of AIC are shown in Table 6. It shows that Nelsen No. 1 has the minimum AIC value for establishing the joint distribution function of and , while Nelsen No. 2 shows the minimum for that of and . Therefore, Nelsen No. 1 and Nelsen No. 2 can be used as the optimal Copula functions to construct the joint distribution functions of and , respectively.

4. Application to Estimating Shear Strength Parameters

4.1. Guarantee Rate

Instability would happen if the strength parameters of rock mass are overvalued, while economic waste could be caused if they are undervalued. Assuming that the estimated values of shear strength parameters are and , the actual values are and ; the guarantee rates and are then defined as the probability of and , respectively, which can be expressed as

The estimated value is overvalued if guarantee rates are low, while being undervalued if high. In consequence, the estimated values can be evaluated by means of the guarantee rates . On the other hand, the estimated and with a certain guarantee rates can be calculated by solving the conditional probability of the optimal Copula functions. Based on the optimal Copula functions of and established in Section 3.3, and can be calculated aswhere is the measured value of .

For direct shear test, in general, the predominant slope method and the least square method [43] are used to statistically process the test data. Using the least square method, the best linear fitting relations of normal stress and shear stress for each class of rock mass quality are obtained, and the slope of fitting straight line is , and the intercept at axial is based on Mohr-Coulomb criterion. The predominant slope method is obviously different from the least square method although they are all based on Mohr-Coulomb criterion and linear fitting method. For the predominant slope method, the upper- and lower-limit values of cohesion should be conducted (Figure 5); the lower one is often suggested. Considering the obvious differences of values in Table 1, 1–6 sets are grouped as Class II while 7–12 sets being grouped as Class III in rock mass quality. Their scatter plots of normal stress versus shear stress are shown in Figure 5. The estimated results using the two methods are shown in Table 7.

As the prescribed standard [44] for the hydroelectric engineering in China, the suggested design parameters must have a guarantee rate of more than 0.8. Therefore, it could be suggested that the of more than 0.95 are too high, while those of less than 0.20 are too low. Relative to the 12 sets test results in Table 1, the guarantee rates and (by (17a) and (17b)) associated with the estimated and (in Table 7) are shown in Table 8. Most of them are not acceptable: all are too low using least square method, which means that the estimated are overvalued; the are larger than 0.95 for 7–12 sets using predominant slope method, which means that the estimated are too undervalued. The data points must be located around the fitted straight line as far as possible on account that the least square method focuses on the best linear fitting of normal stress and shear stress (Figure 5), which can cause the relatively low guarantee rates defined as (17a) and (17b), while the data pairs of and are above the lower-limit fitting line (Figure 5) when the lower-limit values are adopted using the predominant slope method, which can lead to the large guarantee rates . Therefore, the guarantee rate of the value is low by the least square method but high by the predominant slope method.

4.2. Estimated Shear Strength Parameters Using Copula-Based Method

The Copula-based fitting formulas established above are based on the test data of the 12 sets before; therefore, these tested shear strength parameters cannot be used to examine the extensional suitability of the Copula function obtained. The additional in situ direct shear tests are thus performed for the other 9 sets of the marbles with similar properties to check the application of Copula-based method to estimating shear strength parameters of rock mass. The and of the tested samples are also investigated at the same time. The test results are shown in Table 9.

The estimated and using Copula-based method have the guarantee rate of 0.8. The Hoek-Brown criterion (see (1) to (3), the disturbance factor ), as a comparative method, is also adopted to calculate shear strength parameters (in Figure 6). The indicator to assess the predicted suitability of Copula-based method and Hoek-Brown criterion is the root mean square errors (RMSE) obtained by where and are obtained from the in situ shear test and and are derived from Copula-based method or Hoek-Brown criterion. Clearly, the smaller the RMSE is, the more suitable the estimated result is.

Figure 6 indicates that the estimated shear strength parameters by Copula-based method are closer to in situ test values than those by Hoek-Brown criterion for the rock mass specimens in Table 9, which means that the proposed Copula-based method would be a more effective tool in estimating shear strength parameters of the marbles involved in this paper. The Hoek-Brown criterion with a wide range of applications is empirically established on the basis of the test results of a large number of different types of rock mass; the calculated accuracies of and are however seriously affected by the uncertain test conditions of (see (1)) and undetermined disturbance factor and so forth (see (1)–(3)) [45, 46]. However, the optimal Copula-based fitting functions established in the paper are on the basis of probabilistic analysis on the in situ direct shear test results of the 12 sets of similar marbles under the same test conditions, which should show better estimations of shear strength parameters of the rock mass with similar geological and mechanical properties.

5. Conclusion

Using the proposed Copula-based method, the shear strength parameters with a given can be estimated provided that has been obtained. The calculated guarantee rate can be also used to analyze the suitability of the estimated values from the statistical methods for direct shear test data (such as predominant slope method and the least square method). For the marbles tested in this paper, compared with Hoek-Brown criterion, the Copula-based method is a more effective tool in estimating shear strength parameters.

The paper only considers the relationships between two of the parameters. Actually, the Copula functions could be of multidimensions and can be used to evaluate the parameters related with multiple variables. Furthermore, the Copula functions can be used for reliability analysis of estimated rock mass parameters.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This study is supported by the National Natural Science Foundation of China (41172243 and 41472245) and the Fundamental Research Funds for the Central Universities (CDJZR12205501 and 106112014CDJZR200009).