Abstract

A new investigation method is proposed for recording large-sized joint profiles and making statistical analyses of the joint roughness coefficient (JRC) values of the 10–300 cm sized profiles. The mechanical hand profilograph is used for joint roughness measurement due to its advantage of easy operation and high accuracy in recording joint traces. Based on the proposed method, it provides sufficient samples from various positions on the large joint profile, which allows the statistical evaluation of JRC values. A neutrosophic number (NN) is employed for revealing determinate and/or indeterminate information as it consists of determinate and indeterminate parts. Due to the uncertainty of JRC in the real world, NN is chosen to represent the JRC value, which is not only random but also a fuzzy indefinite parameter. The neutrosophic function is used to analyze and express the scale effect of joint surface roughness, and its derivative is used to describe the changing trend of the scale effect. The results show that the JRC value of the joint profile is related to the scale and has a negative effect on the surface roughness of the rock joint. The indeterminate information about the scale effect on joint roughness is described by the neutrosophic functions, and the derivative indicated that the JRC values of small samples are more sensitive than those of large-sized examples. When the length of the sample exceeds the stationarity limit of 80 cm, the roughness appears to be almost scale independent.

1. Introduction

Rock joint roughness plays an important role in hydraulic and mechanical characteristics of discontinuous rock masses [14]. In the past few decades, considerable efforts have been devoted to the investigation of the basic properties of joint roughness. In practice, it is found that roughness is usually scale dependent, and the scale effect on roughness is very important for engineering design [5]. This delivers an urgent need to more accurately determine the relationship between roughness and sample size.

To study the scale effect, the roughness of samples of different sizes should be measured from the same joint surface. Therefore, the method used to measure the roughness should provide a quick and effective measurement of the joint profile in length ranging from a few centimeters to meters. Generally, the joint roughness measurement methods can be classified into “contact methods” and “noncontact methods” [4]. Due to their advantages of easy operation and low cost, the contact method is usually used to measure the topography of rock discontinuities. Weissbach [6] noted that the most reliable method of characterizing a surface in three dimensions remains a combination of a set of linear profile measurements using a stylus instrument. Mogilevskaya [7] measured the macroscopic relief of a joint surface with a VNIIG needle-tip meter. Micallef and Williams [8] used an H-shaped rock profiler to measure locally common smooth and irregular natural rock surfaces. Alameda-Hernández et al. [9] applied a needle profilometer and image-processing method for fast and economical digitization of rock surface profiles. Although these methods can obtain high-resolution and high-precision joint profiles, the allowable size of the measured joint is usually very limited (≤1 m), because the digitization of a large joint is very time-consuming and inconvenient to operate in situ. Recently, the noncontacted instruments have emerged as an alternative for quickly and accurately measuring joint surfaces in both laboratory and situ environments [2]. To investigate the scale-dependency of surface roughness, Fardin et al. [10] digitized a natural rock joint surface replica in the laboratory using a laser scanner with high accuracy and resolution. Feng et al. [11] employed a total station for making high in situ noncontact measurements of fracture roughness. This method was proved to have a similar accuracy level as that of laser scanners. Haneberg [12] introduced a method to obtain outcrop-scale directional roughness profiles based on three-dimensional photogrammetric or laser scanner point clouds. Tatone and Grasselli [2] used a laser scanning system to digitize the surface roughness of large-scale rock fractures in situ. The noncontact methods greatly improved the speed and accuracy of roughness measurements. However, the accuracy and precision of measurement are affected by operation conditions such as surface covering, light requirement, and color distortions caused by rock components [13]. In addition, the measurement noise involved in the data is another fundamental limitation [14].

The surface roughness of rock joints varies with the sample scale and location. The joint surface roughness varies significantly with position, even when measured on the same joint and along the same direction [15]. Ye et al. [16] and Liu et al. [17] studied the roughness of natural joint profiles, and the results showed that the surface roughness differs with positions of the large joint. Roughness heterogeneity affects the results of scale-dependent investigations, especially when joint samples are insufficient. Barton [18] concluded that sampling bias and sampling disturbance influenced the reliability of the previous findings of scale effect on joint roughness. Some researchers have found that statistical analysis of roughness helps reveal the scale effect more accurately. Du et al. [19] analyzed the relationship between the asperity amplitudes of joint profiles and sample size, and a negative scale effect was observed. However, Hsiung et al. [20] stated that the average value of JRC statistical data is difficult to represent the global roughness of a given joint due to its variation. Chen et al. [21] and Wang et al. [22] applied the geological statistics method to investigate the anisotropy and the scale effect of joint roughness. Li et al. [23] studied the fractal dimensions of waviness and unevenness in sampling window sizes ranging from 100 mm × 100 mm to 1000 mm × 1000 mm and found that there is no obvious stationarity threshold for rock joint samples. Yan et al. [24] found that the joint roughness scale effect depends on shear directions, and it tends to be stable when the scale of the study increases to a certain value. Due to the irregular and inconsistent properties of the joint surface, it is still difficult to establish a function for quantitatively expressing the statistical results of roughness to accommodate different sample sizes.

To address these issues, this study represents a new investigation method to analyze the scale dependency of rock joint surface roughness. The rest of the paper is organized as follows. Section 2 contains the literature review on some limitations in previous studies on scale effect. Section 3 contains the measurement method and description of the scale effect of joint surface by neutrosophic functions. Section 4 contains the application of the proposed method in a case study. Section 5 contains the conclusions.

2. Limitations in Previous Studies on Scale Effect

First, the lack of appropriate devices hinders accurate and quick measurement of field-scale joint surfaces. It is widely recognized that joint roughness consists of a large-scale waviness component and a small-scale unevenness component [25]. Investigations of rock joint surface roughness were carried out based on small samples in previous studies. Recently, the optical instrument has been proved to be a good option for measuring large joint surfaces in situ [13, 26, 27]. Yet there are disadvantages, such as the cost of the equipment, relevant software, and costs of training. In rock engineering practices, conventional contact profiling methods are still used in joint roughness measurements. These inexpensive methods can provide repeatable measures of joint traces that are ready to be digitized and analyzed. The simple profile gauge and 10 cm long steel stylus comb are habitually used in the field for estimating the small-scale joint roughness [28]. But these devices/apparatus (e.g., profile comb, electrical profilograph, and needle profilometer) cannot measure the joint profiles with several meters or even tens of meters, in length. As a result, large-scale roughness components are rarely measured and ignore the mechanical impact of field-scale discontinuities [10].

Second, the sampling bias and disturbance affect the reliability of the previous findings concerning some of the apparent scale effects [18]. Numerous potential scale effects are evident in rock mechanics. Unfortunately, conflicting results about scale effects on rock joint surface roughness were produced in previous studies. Indeed, a lot of confusion was caused by the difficulties in obtaining representative samples. Joint surface roughness varies not only from joint to joint but also with scale and the measurement locations [29].

The following example was used to show the JRC values over a 100 cm long specimen. Here, the JRC values of different sized profiles were calculated based on the method proposed by Zhang et al. [30]. The detailed roughness determination process is introduced as follows.

First, assume that a horizontal reference line crosses the selected profile. The average vertical distance between this horizontal reference line and the points on the digitized profile is obtained bywhere L is the digitized profile length and x and y are the x- and y-coordinates of the points on the profile.

Then, the modified root mean square is determined as follows:

Based on the logistic correlation, the JRC value can be determined aswithwhere λ is the roughness index.

This procedure was repeated to determine the JRC evaluation for all profiles.

In Figure 1, a digitized profile was separated using different sampling lengths (10 cm, 20 cm, and 50 cm). Their JRC values are plotted in Figure 2. The JRC values changed with sample length and showed a negative scale effect. In addition, the JRC values change because of the positions of the joint samples. For example, the JRC values of the profiles in length of 10 cm ranged from 7.80 to 20 and from 5.26 to 12.90 for the 20 cm long joint samples. The JRC values of smaller samples showed a higher disparity/variation. Although the JRC mean values of small-sized samples were higher than the larger ones, some individual small samples showed lower roughness than the larger samples. Furthermore, as the surface roughness of joint samples between the adjacent sampling lengths was neglected, the characteristics of joint roughness cannot be considered comprehensively. For example, there were two 50 cm long samples in the range of 0–50 cm and 50–100 cm in Figure 1, but the roughness of the joint profile was in the range of 25–75 cm. Based on these two samples, it was impossible to obtain a JRC value that could truly reflect the integrated characteristic of the 50 cm long joints. Also, only one joint sample could be acquired when the sample size was larger than 50 cm. Thus, researchers may arbitrarily select the sampling regions from a large joint surface. Moreover, inadequate small-sized samples cannot comprehensively represent the real characteristic of joint roughness. It is therefore of urgent need to have generally adopted methods for obtaining representative samples.

3. A New Survey Method

3.1. Profile Recording Equipment

The profiles of rock joints were measured using mechanical hand profilograph (Figure 3). The profilograph is composed of several parts, including feeler, drawing pen, balance block, fixed board, bubble levels, and drawing paper. The top of the feeler, which is 80 mm long, is linked to the end of the balance block via a spring. The bottom consists of a steel needle. A small steel ball is mounted on the feeler point that keeps in contact with the rock joint surface during the measuring process. The drawing pen is fixed on the balance block and synchronized with it, moving up and down depending on the geometry of the joint surface. The fixed board is used to adjust and underlay the drawing paper. The board and paper must be kept perpendicular to the joint surface during measurement to avoid the “sudden jump” plotting problem when imitating the rock joint profile. Meanwhile, the balance block-feeler-drawing pen should also be perpendicular to the fixed board (l1l2, l2l3, l1l4, l1l3, and l2//l4 in Figure 3).

To measure the profiles in lengths of 200 cm, 300 cm, 400 cm, 500 cm, 600 cm, 700 cm, 800 cm, 900 cm, and 1000 cm, stainless steel boards made of aluminum alloy were attached. This board, which was 2 mm thick, 100 cm, or 50 cm in length and 20 cm in width (Figure 4(a) and Figure 4(b)), held the underlying drawing paper in place, providing stiffness and lightweight. Two lines of mounting holes, 5 mm in diameter, were set along the length of the board at an interval of 10 cm, beginning 5 cm from the left end. The distance from the center of the mounting hole to the near margin was 2 cm. The front and back layers are mutually staggered and connected to each other with bolts crossing the mounting holes. For example, a 200 cm long fixed board was assembled as shown in Figures 4(c) and 4(d). The fixed board was then carefully laid on the unfolded drawing paper (Figure 4(e)). The drawing paper was tightly rolled around the fixed board and held place with clamps. This is vital as the paper for recording joint profiles should be completely flat without any wrinkles. By performing this process, the profiles of large-scale rock joints can be copied accurately on the drawing paper. These profiles were then digitized using the grayscale image processing technique. More details of the profilograph can be found in [3]. Finally, the joint roughness was determined according to the digitized profile data. The sampling resolution has a significant impact on the joint roughness measurement [4]. In this study, the joint profiles were digitized with a sampling interval (SI) of 0.5 mm, which is a very common sampling resolution used in previous studies on joint roughness determination.

3.2. Profile Sampling Process

Figures 2 and 3 show that the joint roughness of the profiles differs depending on the sampling locations. These samples, obtained by equally dividing the original joint profile, were inadequate yet. Consequently, they cannot accurately represent the real roughness characteristics. The sampling procedure is shown in Figure 5. First, it is assumed that the number of joint profiles in length of l is n. The overlap length is and the remaining length is :

A large joint profile (L = 800 cm) was considered to show how to extract the small-sized samples (l = 200 cm). Seven samples (n = 7) are required to be obtained from the original joint profile. The overlap of two adjacent samples and the length of the remaining part were calculated using (5) and (6), where the overlapped length  = 100 cm and the length of the remaining part k = 100 cm. Although it is unlike the independent samples obtained by equally dividing the original joint profile, it exposited more refined and integrated roughness characteristics of the rock joint. Samples from different positions of the rock joint were obtained as far as possible. Each sampling occurrence is independent, and thus, the statistical regularities of different sized samples could be analyzed.

3.3. Neutrosophic Function with NNs for Investigating the Scale Dependency of Joint Roughness

Smarandache [3235] proposed NN under the indeterminate environments. It consisted of a certain part and an uncertain part. The former is given as α, and the latter is expressed as βI. Then, NN can be expressed as μ = α + βI for α, β ∈ R. Here, I denotes indeterminacy, and R represents the real numbers. NN indicates a changeable interval number corresponding to different indeterminate ranges of I ∈ [I, I+]. Since NN is flexible and convenient in expressing indeterminate information in indeterminate situations, it has been widely used to represent determinate and/or indeterminate information in many areas.

For any two NNs μ1 = α1 + β1I and μ2 = α2 + β2I for μ1, μ2 ∈ μ, and I ∈ [I, I+], we developed some basic operational laws as follows:

According to the roughness data, it is preferred to use the mean value of all same sized joint samples to represent the joint roughness under this scale. The standard deviations have not been taken into consideration yet, which can quantify the amount of variation or dispersion of a set of joint roughness values. In Figure 1, the JRC values of the joint samples in a length of 10 cm are more discrete than the 20 cm samples. Thus, the standard value of 10 cm long samples should be larger than those of 20 cm in length. The JRC values ranging from JRCmean − σ (lower limit) to JRCmean + σ (upper limit) were chosen to describe the joint roughness properties. The range of JRC values reflected the majority JRC values for the samples of the same size. Generally, the upper and lower limits decreased with the length of rock joint samples increasing. Their relationships were expressed by the following logarithmic functions:where JRCdown and JRCup are the predicted values of both lower and upper limits, and m1, m2, m3, and m4 are the fitting coefficients.

Thus, the neutrosophic function of JRC for I ∈ [0, 0.5] is expressed aswhere m1 and m2 are the fitting coefficients of the lower bound values and nm1 and nm2 are twice the difference of the fitting coefficients between the lower and upper bound values.

According to (9), the variation trend of scale effect can be expressed by the derivative of JRC(l). The expression of the derivative DJRC is

4. Application

In this study, the natural rock joint (Figure 6) was selected from a bedding surface located in Changshan County, Zhejiang Province, China. This slate rock was slightly weathered and had a relatively smooth, planar surface. An 800 cm long joint profile was measured using the mechanical hand profilograph (Figure 7(a)). The obtained joint profile was digitized using the grayscale image processing technique and shown in Figure 7(b).

In order to determine a reasonable sample number, the roughness of joint profiles 10 cm, 50 cm, and 100 cm long was calculated based on the relationship between the surface amplitude asperities and JRC values [35]. The approximate value of field-scale JRC is approximately determined from the measured values of amplitudes and lengths of joint profiles. The number of samples with the same size changed from 10 to 200 at an interval of 10. The mean value and standard deviation of the samples in length of 10 cm, 50 cm, and 100 cm are given in Figure 8, which shows that the statistical results (i.e., mean value and standard deviation) tend to be consistent until the sample number reaches 100. Therefore, a sample number of 100 was taken for the joint profiles of each size in this study.

The statistical results of the JRC values of the joint profiles with the same length were separately analyzed. The frequency distributions and cumulative probabilities of the 10 cm, 20 cm, 50 cm, and 100 cm long joint samples are shown in Figure 9. They are represented graphically by characteristic bell-shaped curves. It was found that when the sample length increased, the JRC mean value and standard deviation decreased. The mean value JRCmean and standard deviation σ of JRC values of the 10 cm long samples were 10.56 and 5.34. Values ranging from JRCmean − σ (lower limit) to JRCmean + σ (upper limit) were chosen to describe the roughness characteristics. The lower and upper limits of 10 cm long samples are 5.22 and 15.91. As shown in Figure 9, 75%, 63%, 67%, and 67% of the data of the samples in different lengths distributed within the interval [JRCmean − σ, JRCmean + σ]. The lower limits, upper limits, and JRC mean values of the joint sample in length from 10 cm to 300 cm were plotted in Figure 10. A decreasing trend was found in both the upper and lower limits as sample size progressively increased, which fit the exponential functions. Thus, m1 and m2 were 19.59 and −2.26 and m3 and m4 were 6.61 and −0.43 (Figure 10). Then,

Thus, the NN function JRC (l) with I ∈ [0, 0.5] for describing the JRC values of different sized joint samples was

Based on (12), the derivative of JRC (L) for indicating the changing rates in different sizes was expressed as

The derivative DJRC shows the rate of change of the JRC values relative to the change of joint profile length. DJRC = 0 means that JRC is almost constant around this length. As shown in Figure 11, the derivative values increased as the sample length increased. DJRC of joint profile in small length changed rapidly, but it tended to be almost constant when the length of the joint profile exceeds a certain value, called the stationarity limit. It indicated that JRC values of small samples were more sensitive than those of larger ones. For samples smaller than this, the joint roughness was scale dependent. Beyond this limit, however, the roughness behavior was scale independent. The sample length was assumed to reach the stationarity limit once DJRC equal to 0.05. Thus, here, the 80 cm sample length represented the stationarity limit. When the sample sizes exceeded 80 cm, and the roughness behavior was practically scale independent.

5. Conclusions and Future Work

Improving the characterization of scale effect on rock joint roughness is an important step for predicting the mechanical behavior of the rock fractures. So far, people have paid considerable efforts to explore the joint roughness scale effect. Several limitations exist in previous studies: (1) the lack of appropriate devices for accurately and quickly measuring large joint surfaces in situ, (2) the representative of test sample was not fully considered, and (3) the statistical data of joint roughness fails to be comprehensively reflected. To address these problems, this study suggested a new method as follows.

A simple mechanical hand profilograph was adopted in this study for the measurement of large rock joint profiles. For measuring the joint profiles several meters long, a fixed board was employed, consisting of two fixed-sized, stainless steel boards that were conveniently assembled.

The roughness along different sections of a joint surface varied significantly. Consequently, a new sampling process was proposed for obtaining sufficient rock joint samples. The relationship between the number of samples and the statistical mean value and standard deviation were considered for determining a reasonable number of joint samples for statistical purposes.

The JRC statistical analysis was used to study the scale dependency of joint roughness. The neutrosophic functions JRC(l) were proposed for expressing the indeterminate information about the scale effect on joint roughness. The upper and lower bounds of JRC values reduced with the sample size, showing a negative scale effect on joint roughness. Based on the neutrosophic functions JRC(l), the variation trend of scale effect was indicated by the derivative of JRC(l).

The scale effect of rock joints was quantitatively expressed by neutrosophic functions, and the derivative indicated that the JRC values of small samples are more sensitive than those of large-sized examples. The stationary limit of joint roughness is 80 cm, when the length of the samples exceeded it, and the roughness behavior was almost scale independent.

In our future work, the direct shear test will be performed for the verification of JRC estimation results of different sized joint specimens. In addition, the applications of the neutrosophic functions in other domains, such as intelligent manufacturing, machine learning, and data mining, will be conducted.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Rui Yong proposed the main idea of this paper. Shaonan Tan and Jun Ye collected the natural rock joint samples and measured different sized joint surfaces. Shigui Du worked on joint roughness assessment and statistical analysis on the JRC values.

Acknowledgments

The study was funded by the National Natural Science Foundation of China (nos. 41502300, 41427802, and 41572299) and Zhejiang Provincial Natural Science Foundation (no. LQ16D020001).