Research Article  Open Access
Correlation between the Joint Roughness Coefficient and Rock Joint Statistical Parameters at Different Sampling Intervals
Abstract
Joint roughness coefficient (JRC) is a major factor that affects the mechanical properties of rock joints. Statistical methods that are used to calculate the JRC increasingly depend on a sampling interval (Δx). The variation rules of fitting parameters a, b, and b/a at different Δx values were analyzed on the basis of the relationship between the JRC and statistical parameter Z_{2}. The relationship between the fitting parameters a and b was deduced in accordance with the ten standard profiles proposed by Barton. Empirical formulas for the JRC, Z_{2}, and Δx were also established. The estimation accuracy of the JRC was the highest in the analysis of Δx values within 0.1–5.0 mm. JRC tests were conducted through inverse value comparative analysis. Results showed that the outcome calculated using the general formula and the JRC inverse values demonstrate improved agreement and verify the rationality of the general formula. The proposed formula can perform rapid and simple JRC calculation within the Δx range of 0.1–5.0 mm using Z_{2}, thereby indicating favorable application prospects.
1. Introduction
The roughness of a rock joint directly affects the strength, deformation, and seepage characteristics of a rock mass. Barton and Choubey [1] established ten typical profiles through a backcalculation test to quantify the joint roughness coefficient (JRC). The JRC can be visually estimated from the standard profiles that correspond to values within [0, 20] in laboratory tests. Subsequently, many researchers have proposed various improvements to calculate JRC values more precisely than visual assessment [2–5]. To facilitate the application of JRCs in rock engineering practices, a series of new parameters and methods, such as statistical methods [6–10], straight edge method [11, 12], fractal dimension method [13–17], and other quantitative methods, has been introduced to evaluate the JRCs of the ten standard profiles. The statistical parameter method has been extensively used in previous studies, given its convenience. This method mainly includes the first derivative root mean square (Z_{2}), the structure function of the profile (SF), the roughness profile index (R_{P}), the root mean square (RMS) roughness index of the profile, the profile elongation index (δ), and other relevant parameters [18]. Z_{2} is the most commonly used among these statistical parameters.
Tse and Cruden [7] established the regression equations between the JRC and Z_{2} using a 0.5 mm sampling interval (Δx). Yu and Vayssade [8] used three Δx values (i.e., 0.25, 0.5, and 1.0 mm) to describe the effect of the sampling interval on the relationship between the JRC and Z_{2}. These researchers determined that Z_{2} cannot be used without considering the effect of sampling intervals. Since then, the relationship between Z_{2} and Δx has been thoroughly studied. Compared with other statistical parameters, Wu et al. [19] suggested that Z_{2} and SF must be adopted to evaluate the JRC value and propose an empirical relationship between the JRC and Z_{2} at Δx = 0.25 mm, 0.5 mm, and 1 mm. In addition, Zhang et al. [20] proposed an appropriate expression of the Z_{2} and JRC for Δx = 0.5, 1.0, 2.0, and 4.0 mm. Although these tests can be conducted with high accuracy, limitations, such as low accuracy, resolution, and ease of use, still exist.
Recently developed scanning devices can digitize JRCs automatically and more accurately than previous technologies. Tatone and Grasselli [21, 22] developed a roughness evaluation methodology and proposed a new relationship between the JRC and Z_{2} at four Δx values (0.044, 0.25, 0.5, and 1 mm) by analyzing 2D profiles and 3D surface topography using a laser scanner. Furthermore, Song et al. [23] analyzed the variation in Z_{2} using Δx and applied a quadratic polynomial to describe the relationship between Z_{2} and Δx. In summary, strong sensitivity to sampling intervals has been extensively accepted. However, the specific expression of the relationship between Z_{2} and Δx remains unclear.
Some studies have suggested that the relationships between the arbitrary sampling intervals, JRC of natural rock joints, and Z_{2} can be quantitatively characterized on the basis of some statistical relationships, such as power functions [24]. This insight provides a new way of establishing a unified expression. Nevertheless, some drawbacks have been observed in previous research. First, the fitting parameters in the current formula vary for each profile. Second, only the Δx upper bound is defined, and the effect of the Δx lower bound on the joint surface roughness evaluation is ignored.
To overcome the limitations of previous studies, the ten standard profiles presented by Barton are digitally extracted to characterize the relationship between the JRC and Z_{2} at different Δx values. The present study establishes a new formula to estimate JRC values. Moreover, a method for the association model and reliability analysis is proposed. Finally, the results are summarized to derive a conclusion.
2. Association Model Establishment
2.1. Statistical Model of the JRC and Z_{2}
The RMS value of the first deviation formula for profile Z_{2} is expressed as follows:where Δx and L are the sampling interval and horizontal length of the profile (mm), respectively, and M is the total number of sampling intervals. A reconstruction method for the ten standard JRC profiles is introduced through AutoCAD and MATLAB coding. First, the standard profiles are imported into AutoCAD for digitization. The image is saved in JPG format and scaled to realworld dimensions using a 10 cm scale bar. Second, a MATLAB procedure is conducted using equation (1). Finally, the images are inputted to obtain the corresponding Z_{2} values of the ten standard profiles, and the JRC–Z_{2} formula is fitted under different Δx values.
Several studies on the fitting curves of the relationship between the JRC and Z_{2} are summarized in Table 1. The fitting curves for each reference can be expressed aswhere f(Z_{2}) is the function of Z_{2}, and a and b are the fitting parameters. Table 1 summarizes the main findings of the proposed calculation models with the same Δx value. When the sampling interval varies, the fitting equations are also altered. This result is consistent with the findings of Yu and Vayssade [8]. In accordance with these results, the compatibility of six empirical models, namely, , lg(Z_{2}), Z_{2}, tan(Z_{2}), tan^{−1}(Z_{2}), and (Z_{2})^{2} is discussed. To evaluate the fitting effect of the six models, ten sampling intervals (0.1–1 mm) are used to obtain the statistical relationship between the JRC and Z_{2}. The calculation results are listed in Table 2. These results indicate that the correlation coefficients of statistical model fittings 1–4 are 0.98, 0.96, 0.96, and 0.96, respectively. By comparison, the correlation coefficients of statistical models 5 and 6 are less than 0.90, thereby presenting relatively poor fitting results. Therefore, statistical models 1–4 can provide the ideal expressions for the following analysis.


2.2. Formula Derivation
Taking model 1 as an example, the Δx values of 0.1, 0.5, and 1 mm are selected to analyze the effect of sampling intervals on the relationship between the JRC and Z_{2}. A unique bestfit curve is established through each sampling interval (Figure 1). The shapes of the JRC–Z_{2} fitting curves under different Δx values are identical. Moreover, these configurations show that the Δx value mainly reflects the variation in fitting parameters a and b in the JRC–Z_{2} fitting curves.
Section 2.1 suggests that the statistical model between the JRC and Z_{2} and the fitting parameters a and b can be determined, but the variation remains unexplored. A MATLAB calculation result is exported to investigate the effect of Δx variations on a and b in Model 1 (Figure 2). The correlation of the fitting parameters with Δx is represented by the following equation: b = 37.95 + 1.23Δx − 1.23a. The relative coefficient is approximately 0.984, thereby indicating that the different values of a and b in the JRC–Z_{2} statistical model depend on the Δx of the profiles. Therefore, the JRC–Z_{2} fitting curve under different Δx values must be predicted.
The graphical variation in a and b as a function of Δx is illustrated in Figures 3 and 4, respectively. The results for a indicate nonlinear growth and eventual stabilization of the increase in Δx. By contrast, b decreases nonlinearly with the increase in Δx. The resulting equations can be expressed as
Generally, when the joint surface is increasingly smooth, the calculated values for the rock joints gradually decrease. Consequently, the power functions of the current study are formulated because the data points obtained gradually decrease with the increase in Δx.
The correlation coefficient obtained by a and Δx using the power function model is very high (R = 0.987) and can be directly used to describe a. However, the correlation coefficient for b and Δx is relatively poor. Therefore, the variation in the relationship between b/a and Δx is recommended for fitting to establish the relationship between b and Δx (Figure 5). The parameter b/a decreases with the increase in Δx. This outcome is consistent with the power function, and the correlation coefficient is 0.987. Consequently, b/a reflects the variation in the fitting parameters with Δx. Therefore, the empirical formulas in a power law form are expressed as
The resulting equations using Model 1 will be used to estimate the relationship between the JRC and Z_{2} for different sampling intervals and can be expressed as
Furthermore, the statistical model of JRC–Z_{2} can be calculated using Δx.
3. Model Comparison Analysis
A new estimation equation for different sampling intervals is proposed to obtain improved empirical relationships among the hatching JRC, a, b, and Δx. This equation provides a quantitative description and demonstrates favorable application prospects. To verify the equation’s accuracy and rationality, Δx = 0.5 mm is applied to disperse the standard JRC profiles. The Δx value is substituted into the new equation to obtain the statistical model between the JRC and Z_{2}. The fitting curve is calculated and then compared with the previously published values using a sampling interval of 0.5. The comparison results are plotted in Figure 6, in which the fitting curve of the JRC and Z_{2} obtained using the proposed empirical formula agrees well with the findings of the previous studies, thus indicating that the new empirical formula has high applicability.
Following the comparison of prediction results, an error analysis for the JRC values of the ten standard profiles is calculated using the proposed empirical formula and the backcalculation value proposed by Barton and Choubey (Figure 7). The relative errors of both calculations are less than 5%, thereby indicating that the data from the current study are consistent with the previously published values for the 0.5 point spacing. To further verify the applicability of different Δx values, nine typical sampling intervals from 0.05 mm to 10 mm are established to digitize the standard JRC profiles dispersedly, and the corresponding backcalculation values are also listed in Table 3. The effective Δx is defined when Δx agrees with equation (6) in determining the upper and lower bounds of the Δx values.
 
Note: bold values in the table are the JRC values with relative errors greater than 5%. 
The error analysis of the data points and the JRC values with relative errors greater than 5% are listed in Table 3. Within the Δx range of 0.1–5 mm, the calculated JRC values fluctuate in the bounded area and agree with the JRC backcalculation values (Figure 8). The errors exceed 5% when Δx is less than 0.1 mm or greater than 5 mm. This finding is mainly because the size of Δx is directly related to different geometric characteristics. The surface morphology of a joint is mainly divided into three categories, namely, macroscopic geometric contours, surface undulating morphology, and microroughness (Figure 9). When Δx is greater than or equal to Δx_{1}, the curve feature reflects the geometric contour of the rock joint, with the largest level reflecting the geometric and macroscopic shapes of the rock joint. The geometric features of the intermediate fluctuation form are neglected by the Z_{2} parameter. When Δx is less than or equal to Δx_{3}, the curve feature reflects the microscopic rough undulating form of the rock joints, with the minimum level surface roughness reflecting only the minor subsurface geometric characteristics on the valley’s surface. When Δx is between Δx_{1} and Δx_{3}, Z_{2} not only reflects the surface morphology fluctuation degree but also the variation in the geometric contour. The JRC value can be accurately calculated using the Z_{2} formula in the range of Δx_{2} because the macroscopic geometric contours and surface undulations of the rock joint surface morphology are the main factors that affect mechanical properties. Combined with the relative error values displayed in Table 3 and Figure 8, the Δx values within 0.1–5 mm agree with the range of values within Δx_{2}. Therefore, the effective Δx of the new empirical formula is suggested to be within this range. Within this range, the proposed empirical formula can accurately describe the relationship between the JRC and Z_{2} at different Δx values.
Models 2, 3, and 4 are also analyzed, and their empirical formulas are summarized in Table 4. The applicable range of Δx values is from 0.1 mm to 5 mm. Figure 10 illustrates the comparison of the four empirical formulas at different Δx values with high coincidence degree, thus denoting that these empirical formulas can accurately describe the relationship between the JRC and Z_{2} at different Δx values. Moreover, the relationship between the JRC and Z_{2} at different Δx values can be simply and quickly calculated using the proposed formula. The SF, R_{p}, δ, and other statistical parameters can be used to analyze the relationship between the JRC and Δx through the abovementioned method.

(a)
(b)
(c)
(d)
4. Conclusions
In this study, four empirical formulas are proposed to establish the relationship between the JRC and Z_{2} at different Δx values, and the capability is validated using the ten standard profiles. The main conclusions are provided as follows:(1)Ten standard profiles are reconstructed using AutoCAD and MATLAB coding to analyze the relationship between the JRC and Z_{2} at different Δx values. The power function relationship between the fitting parameters and the Δx values shows that the sampling interval has a strong connection to the JRC.(2)Four empirical formulas are proposed to estimate the effect of Δx on JRC estimation, which can bypass the limitation of Δx and quickly obtain a JRC at different Δx values. In comparison with the previous formula, the high consistency of the calculated JRC values of the ten standard profiles and the backcalculation values signifies the high applicability of the new empirical formula.(3)The effective Δx of the new empirical formula is obtained by analyzing the relative errors between the data points and the JRC values and the geometric features of the fluctuation. The results show that the proposed empirical formula can accurately describe the relationship between the JRC and Z_{2} at different Δx values when the Δx values are within 0.1–5 mm.
In summary, this study presents a method that can be adopted to calculate the relationship between 2D JRCs and the statistical parameters at different Δx values conveniently. This method has an important practical value in evaluating the sampling effect on the JRCs of rock joints. However, the proposed empirical formula only applies to 2D profiles. Therefore, a 3D joint morphology must be considered in the future to broaden the applicability of the model.
Data Availability
All the data in figures and tables, which are used to support the findings of this study, are included within the article.
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this article.
Acknowledgments
This study was partially funded by the Zhejiang Provincial Natural Science Foundation of China (Grant no. LQ18D020003), National Natural Science Foundation of China (Grant nos. 41831290, 41472248, and 41572299), and Key Research and Development Projects of Zhejiang Province (Grant no. 2019C03104). This support is gratefully acknowledged.
References
 N. Barton and V. Choubey, “The shear strength of rock joints in theory and practice,” Rock Mechanics, vol. 10, no. 12, pp. 1–54, 1977. View at: Publisher Site  Google Scholar
 A. J. Beer, D. Stead, and J. S. Coggan, “Technical note estimation of the joint roughness coefficient (JRC) by visual comparison,” Rock Mechanics and Rock Engineering, vol. 35, no. 1, pp. 65–74, 2002. View at: Publisher Site  Google Scholar
 G. Grasselli, J. Wirth, and P. Egger, “Quantitative threedimensional description of a rough surface and parameter evolution with shearing,” International Journal of Rock Mechanics and Mining Sciences, vol. 39, no. 6, pp. 789–800, 2002. View at: Publisher Site  Google Scholar
 P. AlamedaHernández, J. JiménezPerálvarez, J. A. Palenzuela et al., “Improvement of the JRC calculation using different parameters obtained through a new survey method applied to rock discontinuities,” Rock Mechanics and Rock Engineering, vol. 47, no. 6, pp. 2047–2060, 2014. View at: Publisher Site  Google Scholar
 C.C. Xia, Z.C. Tang, W.M. Xiao, and Y.L. Song, “New peak shear strength criterion of rock joints based on quantified surface description,” Rock Mechanics and Rock Engineering, vol. 47, no. 2, pp. 387–400, 2014. View at: Publisher Site  Google Scholar
 Y. Gao and L. N. Y. Wong, “A modified correlation between roughness parameter Z_{2} and the JRC,” Rock Mechanics and Rock Engineering, vol. 48, no. 1, pp. 387–396, 2015. View at: Publisher Site  Google Scholar
 R. Tse and D. M. Cruden, “Estimating joint roughness coefficients,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 16, no. 5, pp. 303–307, 1979. View at: Publisher Site  Google Scholar
 X. Yu and B. Vayssade, “Joint profiles and their roughness parameters,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 28, no. 4, pp. 333–336, 1991. View at: Publisher Site  Google Scholar
 Z. Y. Yang, C. C. Di, and K. C. Yen, “The effect of asperity order on the roughness of rock joints,” International Journal of Rock Mechanics and Mining Sciences, vol. 38, no. 5, pp. 745–752, 2001. View at: Publisher Site  Google Scholar
 Z. Y. Yang, S. C. Lo, and C. C. Di, “Reassessing the joint roughness coefficient (JRC) estimation using Z_{2},” Rock Mechanics and Rock Engineering, vol. 34, no. 3, pp. 243–251, 2001. View at: Publisher Site  Google Scholar
 N. Barton and S. Bandis, “Review of predictive capabilities of JRCJCS model in engineering practice,” in Proceedings of the Rock Joints, International Symposium on Rock Joints, pp. 603–610, A A Balkema, Loen, Norway, June 1990. View at: Publisher Site  Google Scholar
 N. Turk, M. J. Gerd, W. R. Dearman, and F. F. Amin, “Characterization of rock joint surfaces by fractal dimension,” in Proceedings of the 28th US Symposium on Rock Mechanics, I. Farmer, Ed., pp. 1223–1236, A A Balkema, Tucson, Rotterdam, June 1987. View at: Google Scholar
 S. Du, Y. Hu, and X. Hu, “Measurement of joint roughness coefficient by using profilograph and roughness ruler,” Journal of Earth Science, vol. 20, no. 5, pp. 890–896, 2009, in Chinese. View at: Publisher Site  Google Scholar
 P. H. S. W. Kulatilake, P. Balasingam, J. Park, and R. Morgan, “Natural rock joint roughness quantification through fractal techniques,” Geotechnical and Geological Engineering, vol. 24, no. 5, pp. 1181–1202, 2006. View at: Publisher Site  Google Scholar
 Y. Li and R. Huang, “Relationship between joint roughness coefficient and fractal dimension of rock fracture surfaces,” International Journal of Rock Mechanics and Mining Sciences, vol. 75, pp. 15–22, 2015. View at: Publisher Site  Google Scholar
 A. Malinverno, “A simple method to estimate the fractal dimension of a selfaffine series,” Geophysical Research Letters, vol. 17, no. 11, pp. 1953–1956, 1990. View at: Publisher Site  Google Scholar
 H. Xie, J.A. Wang, and W.H. Xie, “Fractal effects of surface roughness on the mechanical behavior of rock joints,” Chaos, Solitons & Fractals, vol. 8, no. 2, pp. 221–252, 1997. View at: Publisher Site  Google Scholar
 Y. Li and Y. Zhang, “Quantitative estimation of joint roughness coefficient using statistical parameters,” International Journal of Rock Mechanics and Mining Sciences, vol. 77, pp. 27–35, 2015. View at: Publisher Site  Google Scholar
 Y. X. Wu, Q. S. Liu, and X. Y. Liu, “Study of relationship between joint roughness coefficient and statistical parameters,” Chinese Journal of Rock Mechanics and Engineering, vol. 30, pp. 2593–2598, 2011, in Chinese. View at: Google Scholar
 J. M. Zhang, Z. C. Tang, J. D. Jiang, and T. Zhan, “Research on relationship between statistical parameters and JRC,” Science Technology and Engineering, vol. 15, no. 14, pp. 1–5, 2014, in Chinese. View at: Google Scholar
 B. S. A. Tatone and G. Grasselli, “A new 2D discontinuity roughness parameter and its correlation with JRC,” International Journal of Rock Mechanics and Mining Sciences, vol. 47, no. 8, pp. 1391–1400, 2010. View at: Publisher Site  Google Scholar
 B. S. A. Tatone and G. Grasselli, “An investigation of discontinuity roughness scale dependency using highresolution surface measurements,” Rock Mechanics and Rock Engineering, vol. 46, no. 4, pp. 657–681, 2013. View at: Publisher Site  Google Scholar
 L. B. Song, Q. Jiang, Y. H. Li, C. X. Yang, S. Zhong, and M. Z. Zhang, “Stability of statistical parameters of discontinuities morphology and anisotropy based on different sampling intervals,” Rock and Soil Mechanics, vol. 38, no. 4, pp. 1121–1132, 2017, in Chinese. View at: Google Scholar
 F. T. Sun, C. X. She, and L. T. Wan, “Research on relationship between JRC of Barton’s standard profiles and statistics parameters independent of sampling interval,” Chinese Journal of Rock Mechanics and Engineering, vol. 33, no. 2, pp. 3539–3544, 2014, in Chinese. View at: Google Scholar
Copyright
Copyright © 2019 Man Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.