# Scale Dependence of Waviness and Unevenness of Natural Rock Joints through Fractal Analysis

Authors: Yingchun Li; Shengyue Sun; Hongwei Yang
Journal: Geofluids
 (2020)

---

## Abstract

The scale dependence of surface roughness is critical in characterising the hydromechanical properties of field-scale rock joints but is still not well understood, particularly when different orders of roughness are considered. We experimentally reveal the scale dependence of two-order roughness, i.e., waviness and unevenness through fractal parameters using the triangular prism surface area method (TPM). The surfaces of three natural joints of granite with the same dimension of1000mm×1000mmare digitised using a 3D laser scanner at three different measurement resolutions. Waviness and unevenness are quantitatively separated by considering the area variation of joint surface as grid size changes. The corresponding fractal dimensions of waviness and unevenness in sampling window sizes ranging from 100mm×100mm to 1000mm×1000mm at an interval of 100mm×100mm are determined. We find that both the fractal dimensions of waviness and unevenness vary as the window size increases. No obvious stationarity threshold has been found for the three rock joint samples, indicating the surface roughness of natural rock joints should be quantified at the scale of the rock mass in the field.


---

## Body

## 1. Introduction

Rock masses contain a large body of joints. The mechanical and hydraulic behaviours of rock joints highly affect the hydromechanical properties of rock masses. Due to geological processes, rock joints with natural surfaces occur over a broad scale from millimeters to kilometers. Accurate description of joint roughness at the relevant scale is crucial for predicting the hydromechanical coupling of the rock mass.The surfaces of natural rock joints exhibit varying degrees of roughness. Roughness refers to the inherent unevenness and waviness of a rock joint surface relative to its mean plane [1]. Initially, waviness and unevenness represent large-scale undulations observed in the field and small-scale roughness sampled in the laboratory, respectively [2, 3]. Large-scale undulations dominate joint dilatancy since they are too large to be sheared off, and small-scale roughness affects joint shear strength as it is usually damaged under shear. The surface of a laboratory-sized rock joint also exhibits two-order asperities, i.e., first-order waviness and second-order unevenness [4–8]. Waviness with comparatively larger wavelength and amplitude primarily contributes to dilation, whereas unevenness of a smaller asperity size is sheared and damaged, providing shear resistance to the shear movement. That is to say, a rock joint surface is characterised by two-order roughness at various scales [9]. Although numerous empirical and statistical approaches have been proposed to quantify the roughness of rock joints [10–15], they have rarely taken into account the two-order roughness of a joint surface that plays distinct roles in the mechanical and hydraulic behaviours of rock joints [7, 16, 17].The roughness of a natural rock joint surface depends on the scale of examination, which is referred to as scale effect. Bandis et al. [18] reported that the value of JRC (joint roughness coefficient) decreased as the rock joint size increased, i.e., negative scale effect. On the other hand, conflicting results including positive and no scale effects have been observed [19–22]. By examining the morphological characteristics of the surface roughness of a large-scale rock joint replica (1000mm×1000mm), Fardin et al. [23] also stated that there was a stationarity threshold beyond which the scale dependency of surface roughness vanished, i.e., the roughness remained unvaried once the scale exceeded the size of the stationarity threshold. Due to these controversial findings, the nature of how scale affects the surface roughness remains enigmatic.Fractal theory [24, 25] has been successfully applied to characterise the roughness of rock joints at varying scales. Many approaches to estimate the fractal dimension of a rock joint profile have been proposed, including ruler length [26], box counting [27], variogram [28, 29], spectral [30, 31], roughness length [9, 23, 32, 33], and line scaling [34]. The triangular prism surface area method (TPM) [35–37], project covering method (PCM) [38–40], and cubic covering method (CCM) [41] are shown to be applicable for determining the fractal geometry of a three-dimensional joint surface. However, few of them have considered the individual fractal dimensions of waviness and unevenness since a universal single value was commonly assumed.In this paper, we examine the fractal characteristics of waviness and unevenness of three natural granite rock joints dimensioned up to1000mm×1000mm. We find that each-order roughness possesses individual fractal dimension at varying sizes from 100mm×100mm to 1000mm×1000mm. The waviness and unevenness of a rock joint surface are separated by considering the surface area variation as grid size changes. The fractal dimension of each-order roughness is calculated using TPM (triangular prism surface area method). Evident scale dependency of fractal dimension of each-order roughness has been observed. However, the stationarity threshold of the joint surface roughness is possibly absent.


## 2. Data Acquisition

We scanned and reconstructed the three-dimensional surfaces of three natural granite joints 70 (labeledS1, S2, and S3, respectively) sourced from a quarry in Fujian Province, China (Figure 1).Figure 1
Digitisation of a1000mm×1000mm rock joint surface.The dimension of each rock joint surface was1000mm×1000mm. The rock joint surfaces were initially covered by a very thin layer of dust. We carefully cleaned the dust with wipes to avoid damaging the surface roughness. Visual observation suggested that the granite joints are light gray and unweathered.When the rock joint surfaces were naturally dried in the laboratory, a Creaform MetraSCAN 3D 750 system was employed to digitise the joint surface at three measurement resolutions with point spacings being 0.5 mm, 1.0 mm, and 2.0 mm, respectively. The optical scanning system consists of a HandyPROBE for scanning, a C-Track sensor to locate the position of the HandyPROBE, a C-Track controller for data acquisition, and a laptop for image processing and display (Figure1). It ideally can acquire the topographic information of an object up to several meters at a minimum point spacing of 0.05 mm. During data acquisition, we scanned the rock joint surface region by region very slowly and carefully to fully capture the morphological properties of the rock joint surfaces. The digitised surface of the rock joint was visualised simultaneously over scanning through the laptop monitor, which ensured that each surface was thoroughly reconstructed without small empty areas remained. We attempted to obtain a more detailed joint surface at the point spacing of 0.2 mm but failed due to memory limitation of the laptop. The graphic processing software, Geomagic Studio, was employed to coordinate the data acquired through the scanner. The graphic processing software, PolyWorks, converted the format of the data imported from Geomagic Studio to the format that is readable by the data processing software, MATLAB.


## 3. Fractal Dimensions of Waviness and Unevenness

We used the well-established triangular prism surface area method (TPM) [35, 37] to estimate the fractal dimensions of waviness and unevenness. The principle of TPM is that the true surface area of a joint surface is measurable once the heights of all points on the joint surface above a base reference plane are established. For a square grid with a side length of δ (Figure 2), the elevation at the center of the grid cell (h0) is determined by the elevations of its four points:
(1)h0=14hi,j+hi,j+1+hi+1,j+hi+1,j+1,where hi,j, hi,j+1, hi+1,j, and hi+1,j+1 are the elevations of the four points, respectively (Figure 2).Figure 2
Illustration of the triangular prism surface area method (after Clarke [35]). δ is the grid size, S1, S2, S3, and S4 represent the areas of four triangles; i,j, i,j+1, and (i+1, j+1) are, respectively, the coordinates of the four points of the triangular prism, h0 denotes the elevation at the center of the grid cell, and hi,j, hi+1,j, hi,j+1, and hi+1,j+1 are the elevations of the four points, respectively.The area of one of the triangles,S1, is as follows:
(2)S1=l1l1−a1l1−b1l1−c1,where l1=1/2a1+b1+c1,(3)a1=hi,j−hi,j+12+δ2,b1=hi,j−h02+12δ2,c1=hi,j+1−h02+12δ2.Similarly, the areas of the other three triangles, i.e.,S2, S3, and S4, are calculated, respectively. The true area of a joint surface in a given grid cell sized of δ×δ is as follows:
(4)Si,j=S1+S2+S3+S4.The joint surface area is as follows:(5)Sδ=∑i,j=1NδSi,j,where Nδ denotes the number of total grid cells. The joint surface area is a function of grid size (δ) by [37]:
(6)Sδ=Aδ2−D,where D is the fractal dimension of a joint surface, and A is a coefficient. Note that the original approach of Clarke [35] estimated the fractal dimension (D) through the relationship between joint surface area Sδ and grid size square (δ2), i.e., Sδ=Aδ22−D. However, Equation (6) using grid size is proven to be mathematically correct and experimentally reliable [36, 37, 42] as the use of the grid size square underestimated fractal dimension [43].To double-logarithmise Equation (6), we have:
(7)lnSδ=lnA+2−Dlnδ,where D and A are estimated from the slope and intercept of the lnSδ−lnδ plot, respectively.To investigate the scale dependency of surface roughness, the fractal dimensions of waviness and unevenness of the three rock joint samples at varying window sizes from100mm×100mm to 1000mm×1000mm are estimated (Figure 3). The square window of different sizes is selected from the central part of a rock joint surface. We first plot the relationship between lnSδ and lnδ based on Equations (1), (2), (4), (5), (6), and (7). The joint surface area is computed at varying grid sizes through Equation (5). Table 1 shows the grid size used in calculating the joint surface area at varying window sizes of different measurement resolutions from 100mm×100mm to 1000mm×1000mm at the interval of 100mm×100mm. As illustrated in Figure 4, the principle of grid size determination is to ensure that the side length of the sampling window is divisible by the grid size (δ) that is a multiple of the point spacing. Figure 5 demonstrates the double-logarithmic relationship between surface area lnSδ and grid size lnδ of the three rock joint samples at the dimension of 800mm×800mm under the resolution of point spacing at 1.0 mm. The surface areas of the three joint samples are calculated through TPM with each grid size ranging from 1 mm to 400 mm, i.e., δ=400,200,160,100,80,50,40,32,25,20,16,10,8,5,4,2,1×1mm. Waviness and unevenness are separated by considering the area variation of a joint surface at varying grid sizes. Specifically, as the grid size decreases, the joint surface area increases to approximately the real surface area. When the grid size exceeds 30 mm, the slope of the lnSδ−lnδ plot decreases remarkably. Under this circumstance, the joint surface area is primarily contributed by waviness, whereas the surface area of unevenness is excluded. For all the three rock joint samples, the slopes of the lnSδ−lnδ curves vary noticeably at the grid size of 30 mm at which waviness and unevenness are separated. Figure 6 illustrates the decomposition of a rock joint surface into profiles of waviness and unevenness.Figure 3
Illustration of the square windows of different sizes froma=100mm×100mm to j=1000mm×1000mm chosen from the central part of a rock joint surface.Table 1
Grid sizes used in TMP calculation for estimating the fractal dimensions of the rock joint samples of varying window sizes.
Window size (mm×mm)Point spacing (ω) (mm)Grid size (δ) (mm)100×1000.5100,50,40,25,20,10,8,5,4,2,1×0.51.050,25,20,10,5,4,2,1×1.02.050,25,20,10,5,4,2,1×2.0200×2000.5200,100,80,50,40,25,20,16,10,8,5,4,2,1×0.51.0100,50,40,25,20,10,8,5,4,2,1×1.02.050,25,20,10,5,4,2,1×2.0300×3000.5300,200,150,120,100,75,60,50,40,30,25,24,20,15,12,10,8,6,5,4,3,2,1×0.51.0150,100,75,60,50,30,25,20,15,12,10,6,5,4,3,2,1×1.02.075,50,30,25,15,10,6,5,3,2,1×2.0400×4000.5400,200,160,100,80,50,40,32,25,20,16,10,8,5,4,2,1×0.51.0200,100,80,50,40,25,20,16,10,8,5,4,2,1×1.02.0100,50,40,25,20,10,8,5,4,2,1×2.0500×5000.5500,250,200,125,100,50,40,25,20,10,8,5,4,2,1×0.51.0250,125,100,50,25,20,10,5,4,2,1×1.02.0125,50,25,10,5,2,1×2.0600×6000.5600,400,300,240,200,150,120,100,80,75,60,50,48,40,30,25,24,20,16,15,12,10,8,6,5,4,3,2,1×0.51.0300,200,150,120,100,75,60,50,40,30,25,24,20,15,12,10,8,6,5,4,3,2,1×1.02.0150,100,75,60,50,30,25,20,15,12,10,6,5,4,3,2,1×2.0700×7000.5700,350,280,200,175,140,100,70,56,50,40,35,28,25,20,14,10,8,7,5,4,2,1×0.51.0350,175,140,100,70,50,35,28,25,20,14,10,7,5,4,2,1×1.02.0175,70,50,35,25,14,10,7,5,2,1×2.0800×8000.5800,400,320,200,160,100,80,64,50,32,40,25,20,16,10,8,5,4,2,1×0.51.0400,200,160,100,80,50,40,32,25,20,16,10,8,5,4,2,1×1.02.0200,100,80,50,40,25,20,16,10,8,5,4,2,1×2.0900×9000.5900,600,450,360,300,225,200,180,150,120,100,90,75,72,60,50,45,40,36,30,25,24,20,18,15,12,10,9,8,6,5,4,3,2,1×0.51.0450,300,225,180,150,100,90,75,60,50,45,36,30,25,20,18,15,12,10,9,6,5,4,3,2,1×1.02.0225,150,90,75,50,45,30,25,18,15,10,9,6,5,3,2,1×2.01000×10000.51000,500,400,250,200,125,100,80,50,40,25,20,16,10,8,5,4,2,1×0.51.0500,250,200,125,100,50,40,25,20,10,8,5,4,2,1×1.02.0250,125,100,50,25,20,10,5,4,2,1×2.0Figure 4
Illustration of grid size (δ) determination in TPM calculation. The rationale is that the side length of sampling window is divisible by the grid size (δ) that is a multiple of the point spacing.Figure 5
Double-logarithmic relationship between joint surface areaSδ and grid size δ. Dw and Du denote the fractal dimensions of waviness and unevenness, respectively.Figure 6
Decomposition of a rock joint surface into waviness and unevenness.The unevenness is acquired by subtracting the waviness from the whole joint surface. The fractal dimensions of waviness (Dw) and unevenness (Du) of a rock joint surface are determined from the two slopes of each bilinear curve, respectively (Figure 4). Actually, similar bilinearity of the lnSδ−lnδ plot of tension-induced rock joint surfaces has been reported by several researchers [40, 41]. They found that the rock joint surface has nonuniversal fractal dimensions, depending on the measurement scale. However, the nature of the two-order fractal dimensions that are explained above was not unveiled by the authors.


## 4. Results

### 4.1. Scale Effect

Figure7 shows that the fractal dimensions of two-order roughness are scale-dependent for the three rock joint samples digitised at three measurement resolutions. Tables 2–4 list the fractal dimensions of waviness and unevenness of rock joint samples S1 to S3 at varying sizes, respectively. For joint sample S1 at a fixed point spacing, the fractal dimension of waviness is the highest at the sampling window of 100mm×100mm, followed by a decrease once the sampling window grows to 200mm×200mm. As the side length of the sampling window increases to 400 mm, the fractal dimension of waviness peaks with a value smaller than that at the sampling window of 100mm×100mm. When the side length of the sampling window increases from 400 mm to 1000 mm, the fractal dimension of waviness generally decreases with slight fluctuations at the side lengths of 700 mm and 900 mm. The fractal dimension of unevenness of joint sample S1 at a certain point spacing, however, is the smallest at the window size of 100mm×100mm. The fractal dimension of unevenness rises continuously to a peak value as the sampling window size increases to 400mm×400mm, followed by an overall decrease as the sampling window size is increased to the maximum value of 1000mm×1000mm.Figure 7
Fractal dimensions of waviness and unevenness of three rock joints of three resolutions at varying window sizes.Dw and Du are fractal dimensions of waviness and unevenness, respectively. ω denotes the point spacing.
(a)
Two-order fractal dimensions of surface roughness of joint sampleS1(b)
Two-order fractal dimensions of surface roughness of joint sampleS2(c)
Two-order fractal dimensions of surface roughness of joint sampleS3Table 2
Fractal dimensions of waviness and unevenness of rock joint sampleS1 at varying window sizes.
Point spacing (mm)Fractal estimationWindow size (mm×mm)100×100200×200300×300400×400500×500600×600700×700800×800900×9001000×10000.5Dw2.0086772.0040982.0065422.0066242.0058972.0045572.0048992.003482.0038382.0034Aw103604092093180166200258300369600503400653400827600964300R2w0.9990.91160.96210.96710.95610.95990.96940.91790.95430.9056Du2.0084632.0091242.010932.012042.011442.01182.011112.010922.010922.01064Au103504161094500169400263500378800513900669900846800988400R2u0.99060.98730.98930.99070.99040.99090.99080.99060.97830.98861.0Dw2.010182.0048892.005172.0071592.0060782.0045092.0048892.0035462.0038092.003378Aw104204106092560166600258600369500503400653600827500964500R2w0.9990.93650.92640.98180.96860.96410.9820.93670.95890.9172Du2.0078412.0080212.0091092.010762.010162.010412.0097712.0096092.0095382.009236Au10330415009371016880026250037730051190667400843900984600R2w0.99440.99430.99160.99310.99310.9930.99410.9940.99340.99212.0Dw2.0075282.0045912.0063112.0067962.006362.0045212.0051162.0036272.0041352.003681Aw103204100093040166300259000369600503800653800828500966100R2w0.9990.9450.94120.97660.95180.95460.96880.91640.96360.9025Du2.0063112.0064092.0073472.0093432.0085712.0085852.0078492.0078212.0075672.00734Au102804128093470168000261000375200508800663700838800978600R2u0.98790.98940.98910.98430.9930.98570.99480.98780.99030.9905Note:Dw and Du are fractal dimensions of waviness and unevenness, respectively. Aw and Au are the coefficients during linear correlations for estimating Dw and Du, respectively. R2w and R2u represent the coefficients of determination during linear correlations for estimating Dw and Du, respectively.Table 3
Fractal dimensions of waviness and unevenness of rock joint sampleS2 at varying window sizes.
Point spacing (mm)Fractal estimationWindow size (mm×mm)100×100200×200300×300400×400500×500600×600700×700800×800900×90001000×10000.5Dw2.0020652.0046582.0040252.0042222.0044052.0032382.0029182.0030862.0034752.003312Aw101404092091850163500255903669004987065180082650963700R2w0.9990.9710.93430.97280.97560.9810.94420.94490.96640.9589Du2.010162.0098582.010472.010132.0099152.010422.010152.010752.010752.01063Au104104162093790166600260703756002.0101566720084620987600R2u0.98760.98930.98930.98860.98580.989510900.98930.97190.98821.0Dw2.0026962.0043372.0046872.0041362.0039122.0033242.0029572.0030512.0032872.003409Aw101504086092090163400255503671004989065170082590964300R2w0.9990.96520.96540.9540.96230.97770.93990.94650.96040.9578Du2.0080482.0084272.0089472.0088062.0086112.0090512.0087962.0088592.0093292.009198Au103504143093410166000259803742005090066470084320983800R2w0.98490.99190.99020.99240.99070.99280.99340.99320.99210.99082.0Dw2.0039452.0052022.0032652.0047572.0035712.0033082.0028962.0032252.003472.003585Aw101904100091580163900255203670004987065230082640956300R2w0.9990.94990.93620.94180.93840.97710.95470.93870.95380.9417Du2.007172.0072152.0073442.0070522.0068962.007162.0068282.0069522.007332.007264Au103104129092930165100258103721005059066080083810977700R2u0.97990.98950.98850.99190.99270.98860.99270.99090.98940.9899Note:Dw and Du are fractal dimensions of waviness and unevenness, respectively. Aw and Au are the coefficients during linear correlation for estimating Dw and Du, respectively. R2w and R2u represent the coefficients of determination during linear correlations for estimating Dw and Du, respectively.Table 4
Fractal dimensions of waviness and unevenness of rock joint sampleS3 at varying window sizes.
Point spacing (mm)Fractal estimationWindow size (mm×mm)100×100200×200300×300400×400500×500600×600700×700800×800900×9001000×10000.5Dw2.0051992.002592.0043662.004452.0036022.0024772.0025432.002572.0021842.002377Aw101204095091870163500254600364700496700649000820100957800R2w0.9990.96850.92460.96340.95820.91010.92520.90420.91580.9196Du2.010582.010562.011372.010852.010562.011032.010692.010682.010882.01061Au103904158093880166900260600375100510300666400843300984100R2u0.97390.98590.98830.98640.98460.98830.98970.98940.97040.98581.0Dw2.0042822.0032352.0032862.0031042.0036232.0020342.0025572.0028962.0019132.002271Aw10180406009143016260025460364100496800647700819100955500R2w0.9990.90920.91330.90280.93210.90130.91240.9290.91020.9049Du2.008662.0091182.0098792.0092582.0089362.0093832.0090352.0091482.0093122.0091Au103404141093450166100259300373400507900663500839900979800R2w0.98310.98870.98990.98930.990.99270.99290.99270.99150.99012.0Dw2.005642.0038152.0040732.003722.003212.0021922.0029692.0033132.0023312.003075Aw102304069091750163100254300364400497600648500820600963300R2w0.9990.90130.91090.90690.90650.90910.9290.91840.93470.9201Du2.0070112.0068352.0073412.0069062.0069792.0072312.0069032.0071052.0071912.007291Au102804112092820164900257700371000504600659400834700977300R2u0.990.9890.99180.99220.99220.9910.99060.99240.99070.9921Note:Dw and Du are fractal dimensions of waviness and unevenness, respectively. Aw and Au are the coefficients during linear correlation for estimating Dw and Du, respectively. R2w and R2u represent the coefficients of determination during linear correlations for estimating Dw and Du, respectively.For rock joint sampleS2 under a certain measurement resolution, the fractal dimension of waviness fluctuates slightly as the window side length increases from 100 mm to 400 mm, preceding a gradual decrease when the window size grows to 700mm×700mm. As the sampling window size increases from 700mm×700mm to 100mm×1000mm, the fractal dimension of waviness increases marginally. The fractal dimension of unevenness seemingly exhibits no general tendency. The magnitude of the fractal dimension of unevenness roughly levels off with several unremarkable fluctuations at different sampling window sizes.For rock joint sampleS3 at a fixed point spacing, the fractal dimension of waviness is maximum at the window size of 100mm×100mm and then predominantly declines as the window side length increases to 1000 mm with slight fluctuations. The variation of fractal dimension of unevenness resembles that of rock joint sample S2 without noticeable trend.To quantify the variation of fractal dimensions of waviness and unevenness as the window size is enlarged, the percent error relative to the value of window size of100mm×100mm is calculated as follows:
(8)δi=Di−D100D100×100%,where δi and Di represent the percent error and fractal dimension of waviness or unevenness at a window size between 200mm×200mm to 1000mm×1000mm, respectively. D100 is the fractal dimension of waviness or unevenness at the window size of 100mm×100mm.Figure8 presents the percent errors of waviness and unevenness of the three rock joint samples at window sizes from 200mm×200mm to 1000mm×1000mm at three measurement resolutions. Generally, the effect of window size on the fractal dimension of waviness is more pronounced than that of unevenness. Particularly for rock joint samples S2 and S3, the percent errors of fractal dimension of unevenness are lower than 0.1%. The variations of fractal dimensions of both waviness and unevenness of rock joint sample S1 are overall larger than those of rock joint samples S2 and S3 at the same window size and resolution. For rock joint sample S1, the percent errors of the fractal dimension of waviness at varying window sizes under the resolution of 1.0 mm are unanimously the highest, and 2.0 mm the lowest. The percent error of the fractal dimension of unevenness under the resolution of 0.5 mm is the highest, and 2.0 mm the lowest except at the window size of 400mm×400mm where the percent error of the fractal dimension of unevenness under the point spacing of 1.0 mm is marginally smaller than that under 2.0 mm. For rock joint sample S2, the percent errors of the fractal dimension of waviness at varying window sizes under the resolution of 0.5 mm generally are the highest, and 2.0 mm the lowest except at the window side lengths of 300 mm and 800 mm. The percent errors of the fractal dimension of unevenness of rock joint samples S2 and S3 at the resolution of 1.0 mm are the highest.Figure 8
Effect of window size on the fractal dimensions of waviness and unevenness of three rock joints of three resolutions. Percent errors are relative to the values of window size of100mm×100mm. Dw and Du are fractal dimensions of waviness and unevenness. ω denotes the point spacing.
(a)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS1(b)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS2(c)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS3Figures7 and 8 show that the fractal dimension of each-order roughness varies from 2.001 to 2.014, and the percent errors of fractal dimensions of waviness and unevenness of the three rock joint samples are numerically small, which are less than 1%. One may draw the conclusion that the scale effect of the fractal dimensions of waviness and unevenness could be neglected. Actually, the low values of percent errors result from the low values of the fractal dimensions (Figure 7) which are common for the rough surfaces of naturally formed rock joints [39, 41, 44]. Many studies reported that the fractal geometry of the surface roughness of three-dimensional rock joints is slightly larger than 2.0 [38–41]. Zhou and Xie [41] showed that the fractal dimensions of the surface roughness of tension-induced rock joints of varying degrees of roughness are all smaller than 2.07. Similarly, several rock joints collected from in situ also exhibited surface roughness of fractal dimensions around 2.05, whereas the JRC (joint roughness coefficient) values [10] of these joint surfaces reached as high as 14.0 [38]. That is to say, a naturally surfaced rock joint possesses a fractal dimension varying in a very narrow band. A small variation in fractal dimension likely leads to noticeable change of surface roughness [26, 45].The low values of the percent errors of fractal dimensions do not necessarily mean that the variation of surface roughness is negligible as window size change, because the widely used indicators of surface roughness such asJRC and asperity slope can be mathematically related with fractal dimension through certain relationships [45, 46]. These relationships commonly involve scaling coefficients of quite high values [45].To illustrate the effect of fractal dimension variation on the joint roughness change, the well-established relationship linking fractal dimension (D) and JRC is employed [26]:
(9)JRC=−0.87804+37.7844D−10.015−16.9304D−10.0152.Equation (9) was originally proposed to estimate JRC through the fractal dimension (D) of a two-dimensional joint profile. Since the relationship between JRC and the fractal dimension (D) of a three-dimensional joint surface is unavailable, the above formulation is directly adopted by extending the two-dimension to three-dimension by replacing (D−1) with (D−2). Additionally, our purpose is not to quantify JRC through fractal dimension (D), but to demonstrate the significant joint roughness variation due to the small change of fractal dimension. Considering the two-order roughness separation, JRC values of waviness and unevenness (JRCw and JRCu) are, respectively,
(10)JRCw=−0.87804+37.7844Dw−20.015−16.9304Dw−20.0152,JRCu=−0.87804+37.7844Du−20.015−16.9304Du−20.0152.Figure9 shows the effect of window size on the percent errors of JRC values of waviness and unevenness (JRCw and JRCu) of the three rock joint samples of different measurement resolutions. For joint samples S1 to S3, JRC values of waviness (JRCw) exhibit strong scale dependency without general trend. The maximum percent errors of JRCw of joint samples S1, S2, and S3 occur on window sizes of 1000mm×1000mm, 200mm×200mm, and 600mm×600mm, respectively, suggesting the randomness of the scale effect. Additionally, similar to the fractal dimension, the effect of window size on JRCw is much more pronounced than on JRCu. For joint sample S1, the percent errors of JRCw are generally appreciably larger than those of JRCu with a maximum value smaller than 40%. For joint samples S2 and S3, the percent errors of JRCw can be as high as 130%, whereas the percent errors of JRCu are all lower than 15% with a majority less than 10%, which indicates that the scale effect of JRCu may be insignificant.Figure 9
Effect of window size on theJRC values of waviness and unevenness of three rock joints of three resolutions. Percent errors are relative to the values of window size of 100mm×100mm. JRCw and JRCu are the JRC values of waviness and unevenness, respectively. ω denotes the point spacing.
(a)
Percent errors of theJRC values of waviness and unevenness of joint sample S1(b)
Percent errors of theJRC values of waviness and unevenness of joint sample S2(c)
Percent errors of theJRC values of waviness and unevenness of joint sample S3


### 4.2. Effect of Measurement Resolution

We digitised the surface of each rock joint sample using three different resolutions. Figure10 shows the percent errors of the fractal dimensions of waviness and unevenness of the three rock joint samples at resolutions of 1.0 mm and 2.0 mm, respectively, relative to the point spacing of 0.5 mm. Both the fractal dimensions of waviness and unevenness are dependent of measurement resolution. The fractal dimension of unevenness is much more sensitive to the measurement resolution compared with that of waviness. For all the three samples, the fractal dimension of unevenness is the largest at the highest resolution and vice versa. The percent error of fractal dimension of unevenness under the point spacing of 2.0 mm is roughly two times that under point spacing of 0.5 mm. Nevertheless, the fractal dimension of waviness of all the three samples is seemingly unaffected by the resolution. In this study, waviness are asperities in wavelength longer than 30 mm that is substantially greater than the prescribed point spacings. The fractal dimension of waviness is theoretically independent on the resolution of point spacing less than 30 mm. In many cases, the difference of fractal dimension of waviness of the three samples is negligible (Figures 8 and 11). The minor discrepancies mainly originated from the fact that it is impossible to acquire exactly the same coordinates of a certain point due to the error caused by the manual operation of the HandyPROBE. For the example, in Figure 12, the coordinates (400 mm, 800 mm) of a point of a profile are perhaps digitised as (402 mm, 802 mm), (399 mm, 801 mm), and (402 mm, 798 mm) during the scanning with point spacings of 0.5 mm, 1.0 mm, and 2.0 mm, respectively.Figure 10
Effect of resolution on the fractal dimensions of waviness and unevenness of three rock joints of varying sampling window sizes. Percent errors are calculated relative to the values of resolution of 0.5 mm.Dw and Du are fractal dimensions of waviness and unevenness, respectively. ω denotes the point spacing.
(a)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS1(b)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS2(c)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS3Figure 11
Effect of resolution on theJRC values of waviness and unevenness of three rock joint samples of varying sampling window sizes. Percent errors are calculated relative to the values of resolution of 0.5 mm. JRCw and JRCu are the JRC values of waviness and unevenness, respectively. ω denotes the point spacing.
(a)
Percent errors of theJRC values of waviness and unevenness of joint sample S1(b)
Percent errors of theJRC values of waviness and unevenness of joint sample S2(c)
Percent errors of theJRC values of waviness and unevenness of joint sample S3Figure 12
Illustration of the coordinate deviation due to manual operation.Figures8 and 10 show that both window size and measurement resolution affect the fractal dimensions of waviness and unevenness. The measurement resolution has no less influence than window size on the variation of fractal dimension. Previous studies [18–23] have reported positive, negative, and no scale effects when large-scale rock joint surfaces are measured. The controversial results may be stemmed from the combinative effects of sample size and measurement resolution. Most of the studies on large-scale rock joints in the literature [21, 23, 47] used inconsistent measurement resolutions to digitise rock joint surfaces of varying scales. In these studies, the point spacing was increased as larger areas (or lengths) of the rock joint surface were examined. Therefore, the change in roughness with increased rock joint size reported in this work was also possibly a result of the varying measurement resolution besides the variation in roughness [22].


## 4.1. Scale Effect

Figure7 shows that the fractal dimensions of two-order roughness are scale-dependent for the three rock joint samples digitised at three measurement resolutions. Tables 2–4 list the fractal dimensions of waviness and unevenness of rock joint samples S1 to S3 at varying sizes, respectively. For joint sample S1 at a fixed point spacing, the fractal dimension of waviness is the highest at the sampling window of 100mm×100mm, followed by a decrease once the sampling window grows to 200mm×200mm. As the side length of the sampling window increases to 400 mm, the fractal dimension of waviness peaks with a value smaller than that at the sampling window of 100mm×100mm. When the side length of the sampling window increases from 400 mm to 1000 mm, the fractal dimension of waviness generally decreases with slight fluctuations at the side lengths of 700 mm and 900 mm. The fractal dimension of unevenness of joint sample S1 at a certain point spacing, however, is the smallest at the window size of 100mm×100mm. The fractal dimension of unevenness rises continuously to a peak value as the sampling window size increases to 400mm×400mm, followed by an overall decrease as the sampling window size is increased to the maximum value of 1000mm×1000mm.Figure 7
Fractal dimensions of waviness and unevenness of three rock joints of three resolutions at varying window sizes.Dw and Du are fractal dimensions of waviness and unevenness, respectively. ω denotes the point spacing.
(a)
Two-order fractal dimensions of surface roughness of joint sampleS1(b)
Two-order fractal dimensions of surface roughness of joint sampleS2(c)
Two-order fractal dimensions of surface roughness of joint sampleS3Table 2
Fractal dimensions of waviness and unevenness of rock joint sampleS1 at varying window sizes.
Point spacing (mm)Fractal estimationWindow size (mm×mm)100×100200×200300×300400×400500×500600×600700×700800×800900×9001000×10000.5Dw2.0086772.0040982.0065422.0066242.0058972.0045572.0048992.003482.0038382.0034Aw103604092093180166200258300369600503400653400827600964300R2w0.9990.91160.96210.96710.95610.95990.96940.91790.95430.9056Du2.0084632.0091242.010932.012042.011442.01182.011112.010922.010922.01064Au103504161094500169400263500378800513900669900846800988400R2u0.99060.98730.98930.99070.99040.99090.99080.99060.97830.98861.0Dw2.010182.0048892.005172.0071592.0060782.0045092.0048892.0035462.0038092.003378Aw104204106092560166600258600369500503400653600827500964500R2w0.9990.93650.92640.98180.96860.96410.9820.93670.95890.9172Du2.0078412.0080212.0091092.010762.010162.010412.0097712.0096092.0095382.009236Au10330415009371016880026250037730051190667400843900984600R2w0.99440.99430.99160.99310.99310.9930.99410.9940.99340.99212.0Dw2.0075282.0045912.0063112.0067962.006362.0045212.0051162.0036272.0041352.003681Aw103204100093040166300259000369600503800653800828500966100R2w0.9990.9450.94120.97660.95180.95460.96880.91640.96360.9025Du2.0063112.0064092.0073472.0093432.0085712.0085852.0078492.0078212.0075672.00734Au102804128093470168000261000375200508800663700838800978600R2u0.98790.98940.98910.98430.9930.98570.99480.98780.99030.9905Note:Dw and Du are fractal dimensions of waviness and unevenness, respectively. Aw and Au are the coefficients during linear correlations for estimating Dw and Du, respectively. R2w and R2u represent the coefficients of determination during linear correlations for estimating Dw and Du, respectively.Table 3
Fractal dimensions of waviness and unevenness of rock joint sampleS2 at varying window sizes.
Point spacing (mm)Fractal estimationWindow size (mm×mm)100×100200×200300×300400×400500×500600×600700×700800×800900×90001000×10000.5Dw2.0020652.0046582.0040252.0042222.0044052.0032382.0029182.0030862.0034752.003312Aw101404092091850163500255903669004987065180082650963700R2w0.9990.9710.93430.97280.97560.9810.94420.94490.96640.9589Du2.010162.0098582.010472.010132.0099152.010422.010152.010752.010752.01063Au104104162093790166600260703756002.0101566720084620987600R2u0.98760.98930.98930.98860.98580.989510900.98930.97190.98821.0Dw2.0026962.0043372.0046872.0041362.0039122.0033242.0029572.0030512.0032872.003409Aw101504086092090163400255503671004989065170082590964300R2w0.9990.96520.96540.9540.96230.97770.93990.94650.96040.9578Du2.0080482.0084272.0089472.0088062.0086112.0090512.0087962.0088592.0093292.009198Au103504143093410166000259803742005090066470084320983800R2w0.98490.99190.99020.99240.99070.99280.99340.99320.99210.99082.0Dw2.0039452.0052022.0032652.0047572.0035712.0033082.0028962.0032252.003472.003585Aw101904100091580163900255203670004987065230082640956300R2w0.9990.94990.93620.94180.93840.97710.95470.93870.95380.9417Du2.007172.0072152.0073442.0070522.0068962.007162.0068282.0069522.007332.007264Au103104129092930165100258103721005059066080083810977700R2u0.97990.98950.98850.99190.99270.98860.99270.99090.98940.9899Note:Dw and Du are fractal dimensions of waviness and unevenness, respectively. Aw and Au are the coefficients during linear correlation for estimating Dw and Du, respectively. R2w and R2u represent the coefficients of determination during linear correlations for estimating Dw and Du, respectively.Table 4
Fractal dimensions of waviness and unevenness of rock joint sampleS3 at varying window sizes.
Point spacing (mm)Fractal estimationWindow size (mm×mm)100×100200×200300×300400×400500×500600×600700×700800×800900×9001000×10000.5Dw2.0051992.002592.0043662.004452.0036022.0024772.0025432.002572.0021842.002377Aw101204095091870163500254600364700496700649000820100957800R2w0.9990.96850.92460.96340.95820.91010.92520.90420.91580.9196Du2.010582.010562.011372.010852.010562.011032.010692.010682.010882.01061Au103904158093880166900260600375100510300666400843300984100R2u0.97390.98590.98830.98640.98460.98830.98970.98940.97040.98581.0Dw2.0042822.0032352.0032862.0031042.0036232.0020342.0025572.0028962.0019132.002271Aw10180406009143016260025460364100496800647700819100955500R2w0.9990.90920.91330.90280.93210.90130.91240.9290.91020.9049Du2.008662.0091182.0098792.0092582.0089362.0093832.0090352.0091482.0093122.0091Au103404141093450166100259300373400507900663500839900979800R2w0.98310.98870.98990.98930.990.99270.99290.99270.99150.99012.0Dw2.005642.0038152.0040732.003722.003212.0021922.0029692.0033132.0023312.003075Aw102304069091750163100254300364400497600648500820600963300R2w0.9990.90130.91090.90690.90650.90910.9290.91840.93470.9201Du2.0070112.0068352.0073412.0069062.0069792.0072312.0069032.0071052.0071912.007291Au102804112092820164900257700371000504600659400834700977300R2u0.990.9890.99180.99220.99220.9910.99060.99240.99070.9921Note:Dw and Du are fractal dimensions of waviness and unevenness, respectively. Aw and Au are the coefficients during linear correlation for estimating Dw and Du, respectively. R2w and R2u represent the coefficients of determination during linear correlations for estimating Dw and Du, respectively.For rock joint sampleS2 under a certain measurement resolution, the fractal dimension of waviness fluctuates slightly as the window side length increases from 100 mm to 400 mm, preceding a gradual decrease when the window size grows to 700mm×700mm. As the sampling window size increases from 700mm×700mm to 100mm×1000mm, the fractal dimension of waviness increases marginally. The fractal dimension of unevenness seemingly exhibits no general tendency. The magnitude of the fractal dimension of unevenness roughly levels off with several unremarkable fluctuations at different sampling window sizes.For rock joint sampleS3 at a fixed point spacing, the fractal dimension of waviness is maximum at the window size of 100mm×100mm and then predominantly declines as the window side length increases to 1000 mm with slight fluctuations. The variation of fractal dimension of unevenness resembles that of rock joint sample S2 without noticeable trend.To quantify the variation of fractal dimensions of waviness and unevenness as the window size is enlarged, the percent error relative to the value of window size of100mm×100mm is calculated as follows:
(8)δi=Di−D100D100×100%,where δi and Di represent the percent error and fractal dimension of waviness or unevenness at a window size between 200mm×200mm to 1000mm×1000mm, respectively. D100 is the fractal dimension of waviness or unevenness at the window size of 100mm×100mm.Figure8 presents the percent errors of waviness and unevenness of the three rock joint samples at window sizes from 200mm×200mm to 1000mm×1000mm at three measurement resolutions. Generally, the effect of window size on the fractal dimension of waviness is more pronounced than that of unevenness. Particularly for rock joint samples S2 and S3, the percent errors of fractal dimension of unevenness are lower than 0.1%. The variations of fractal dimensions of both waviness and unevenness of rock joint sample S1 are overall larger than those of rock joint samples S2 and S3 at the same window size and resolution. For rock joint sample S1, the percent errors of the fractal dimension of waviness at varying window sizes under the resolution of 1.0 mm are unanimously the highest, and 2.0 mm the lowest. The percent error of the fractal dimension of unevenness under the resolution of 0.5 mm is the highest, and 2.0 mm the lowest except at the window size of 400mm×400mm where the percent error of the fractal dimension of unevenness under the point spacing of 1.0 mm is marginally smaller than that under 2.0 mm. For rock joint sample S2, the percent errors of the fractal dimension of waviness at varying window sizes under the resolution of 0.5 mm generally are the highest, and 2.0 mm the lowest except at the window side lengths of 300 mm and 800 mm. The percent errors of the fractal dimension of unevenness of rock joint samples S2 and S3 at the resolution of 1.0 mm are the highest.Figure 8
Effect of window size on the fractal dimensions of waviness and unevenness of three rock joints of three resolutions. Percent errors are relative to the values of window size of100mm×100mm. Dw and Du are fractal dimensions of waviness and unevenness. ω denotes the point spacing.
(a)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS1(b)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS2(c)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS3Figures7 and 8 show that the fractal dimension of each-order roughness varies from 2.001 to 2.014, and the percent errors of fractal dimensions of waviness and unevenness of the three rock joint samples are numerically small, which are less than 1%. One may draw the conclusion that the scale effect of the fractal dimensions of waviness and unevenness could be neglected. Actually, the low values of percent errors result from the low values of the fractal dimensions (Figure 7) which are common for the rough surfaces of naturally formed rock joints [39, 41, 44]. Many studies reported that the fractal geometry of the surface roughness of three-dimensional rock joints is slightly larger than 2.0 [38–41]. Zhou and Xie [41] showed that the fractal dimensions of the surface roughness of tension-induced rock joints of varying degrees of roughness are all smaller than 2.07. Similarly, several rock joints collected from in situ also exhibited surface roughness of fractal dimensions around 2.05, whereas the JRC (joint roughness coefficient) values [10] of these joint surfaces reached as high as 14.0 [38]. That is to say, a naturally surfaced rock joint possesses a fractal dimension varying in a very narrow band. A small variation in fractal dimension likely leads to noticeable change of surface roughness [26, 45].The low values of the percent errors of fractal dimensions do not necessarily mean that the variation of surface roughness is negligible as window size change, because the widely used indicators of surface roughness such asJRC and asperity slope can be mathematically related with fractal dimension through certain relationships [45, 46]. These relationships commonly involve scaling coefficients of quite high values [45].To illustrate the effect of fractal dimension variation on the joint roughness change, the well-established relationship linking fractal dimension (D) and JRC is employed [26]:
(9)JRC=−0.87804+37.7844D−10.015−16.9304D−10.0152.Equation (9) was originally proposed to estimate JRC through the fractal dimension (D) of a two-dimensional joint profile. Since the relationship between JRC and the fractal dimension (D) of a three-dimensional joint surface is unavailable, the above formulation is directly adopted by extending the two-dimension to three-dimension by replacing (D−1) with (D−2). Additionally, our purpose is not to quantify JRC through fractal dimension (D), but to demonstrate the significant joint roughness variation due to the small change of fractal dimension. Considering the two-order roughness separation, JRC values of waviness and unevenness (JRCw and JRCu) are, respectively,
(10)JRCw=−0.87804+37.7844Dw−20.015−16.9304Dw−20.0152,JRCu=−0.87804+37.7844Du−20.015−16.9304Du−20.0152.Figure9 shows the effect of window size on the percent errors of JRC values of waviness and unevenness (JRCw and JRCu) of the three rock joint samples of different measurement resolutions. For joint samples S1 to S3, JRC values of waviness (JRCw) exhibit strong scale dependency without general trend. The maximum percent errors of JRCw of joint samples S1, S2, and S3 occur on window sizes of 1000mm×1000mm, 200mm×200mm, and 600mm×600mm, respectively, suggesting the randomness of the scale effect. Additionally, similar to the fractal dimension, the effect of window size on JRCw is much more pronounced than on JRCu. For joint sample S1, the percent errors of JRCw are generally appreciably larger than those of JRCu with a maximum value smaller than 40%. For joint samples S2 and S3, the percent errors of JRCw can be as high as 130%, whereas the percent errors of JRCu are all lower than 15% with a majority less than 10%, which indicates that the scale effect of JRCu may be insignificant.Figure 9
Effect of window size on theJRC values of waviness and unevenness of three rock joints of three resolutions. Percent errors are relative to the values of window size of 100mm×100mm. JRCw and JRCu are the JRC values of waviness and unevenness, respectively. ω denotes the point spacing.
(a)
Percent errors of theJRC values of waviness and unevenness of joint sample S1(b)
Percent errors of theJRC values of waviness and unevenness of joint sample S2(c)
Percent errors of theJRC values of waviness and unevenness of joint sample S3


## 4.2. Effect of Measurement Resolution

We digitised the surface of each rock joint sample using three different resolutions. Figure10 shows the percent errors of the fractal dimensions of waviness and unevenness of the three rock joint samples at resolutions of 1.0 mm and 2.0 mm, respectively, relative to the point spacing of 0.5 mm. Both the fractal dimensions of waviness and unevenness are dependent of measurement resolution. The fractal dimension of unevenness is much more sensitive to the measurement resolution compared with that of waviness. For all the three samples, the fractal dimension of unevenness is the largest at the highest resolution and vice versa. The percent error of fractal dimension of unevenness under the point spacing of 2.0 mm is roughly two times that under point spacing of 0.5 mm. Nevertheless, the fractal dimension of waviness of all the three samples is seemingly unaffected by the resolution. In this study, waviness are asperities in wavelength longer than 30 mm that is substantially greater than the prescribed point spacings. The fractal dimension of waviness is theoretically independent on the resolution of point spacing less than 30 mm. In many cases, the difference of fractal dimension of waviness of the three samples is negligible (Figures 8 and 11). The minor discrepancies mainly originated from the fact that it is impossible to acquire exactly the same coordinates of a certain point due to the error caused by the manual operation of the HandyPROBE. For the example, in Figure 12, the coordinates (400 mm, 800 mm) of a point of a profile are perhaps digitised as (402 mm, 802 mm), (399 mm, 801 mm), and (402 mm, 798 mm) during the scanning with point spacings of 0.5 mm, 1.0 mm, and 2.0 mm, respectively.Figure 10
Effect of resolution on the fractal dimensions of waviness and unevenness of three rock joints of varying sampling window sizes. Percent errors are calculated relative to the values of resolution of 0.5 mm.Dw and Du are fractal dimensions of waviness and unevenness, respectively. ω denotes the point spacing.
(a)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS1(b)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS2(c)
Percent errors of fractal dimensions of waviness and unevenness of joint sampleS3Figure 11
Effect of resolution on theJRC values of waviness and unevenness of three rock joint samples of varying sampling window sizes. Percent errors are calculated relative to the values of resolution of 0.5 mm. JRCw and JRCu are the JRC values of waviness and unevenness, respectively. ω denotes the point spacing.
(a)
Percent errors of theJRC values of waviness and unevenness of joint sample S1(b)
Percent errors of theJRC values of waviness and unevenness of joint sample S2(c)
Percent errors of theJRC values of waviness and unevenness of joint sample S3Figure 12
Illustration of the coordinate deviation due to manual operation.Figures8 and 10 show that both window size and measurement resolution affect the fractal dimensions of waviness and unevenness. The measurement resolution has no less influence than window size on the variation of fractal dimension. Previous studies [18–23] have reported positive, negative, and no scale effects when large-scale rock joint surfaces are measured. The controversial results may be stemmed from the combinative effects of sample size and measurement resolution. Most of the studies on large-scale rock joints in the literature [21, 23, 47] used inconsistent measurement resolutions to digitise rock joint surfaces of varying scales. In these studies, the point spacing was increased as larger areas (or lengths) of the rock joint surface were examined. Therefore, the change in roughness with increased rock joint size reported in this work was also possibly a result of the varying measurement resolution besides the variation in roughness [22].


## 5. Discussion and Implications

We explored the fractal characteristics of two-order asperities of natural rock joints in window sizes from100mm×100mm to 1000mm×1000mm with three different resolutions.Both the fractal dimensions of waviness and unevenness vary as the window size increases.Nevertheless, no apparent stationarity threshold has been found for the three natural granite joints. The findings contrast with the results reported in Fardin et al. [23]. Using the roughness length method [33], Fardin et al. [23] studied the fractal dimension of the surface roughness of a rock joint up to the window size of 1000 mm ×1000 mm.They claimed that for sampling windows having a size larger than a stationarity threshold of roughly500mm×500mm, the fractal dimension remained almost unvaried and can be considered as reliable estimation. The discrepancy between Fardin et al. [23] and our study possibly resulted from the following reasons. First, a single fractal dimension, i.e., the fractal dimension of unevenness, was estimated by Fardin et al. [23]. In the calculation of the roughness length method, the local trend, the asperities with long wavelength (which is termed waviness in this study) are excluded to avoid overestimating roughness in small windows [32, 33]. Second, the fractal dimension of unevenness with relatively small and similar wavelength and amplitude probably keeps almost constant as the sampling data reaches a certain large volume since fractal dimension is essentially determined based on statistical consideration. For example, the fractal dimensions of unevenness of rock joint samples S2 and S3 in this study (Figures 7(b) and 7(c)) appear unchanged when the window size exceeds 400mm×400mm. Moreover, the surface of merely one rock joint was studied by Fardin et al. [23], based on which affirmative conclusions cannot be drawn without examining the fractal features of several more large-scale rock joints with varying surface roughness.Quantification of the surface roughness of a natural rock joint is critical for predicting the mechanical and hydraulic properties of a rock mass. Particularly, the shear behaviour of rock joints is strongly affected by surface roughness. Under a low normal stress (relative to rock strength), a rock joint fails in shear due to asperity dilation. In this case, both waviness and unevenness contribute to the degree of dilatancy. When the normal stress is high, joint shear behaviour is mainly controlled by the damage or degradation of waviness since unevenness is much more easily sheared off. The permeability of rock joints under shearing is closely associated with the dilation and degradation of asperities [48–53] . That is to say, waviness and unevenness of the surface roughness of a natural rock joint should be separately considered for accurately predicting the hydromechanical behaviour of a rock mass. Currently, surface roughness of rock joints at field scales is commonly evaluated from laboratory experimentation on small joint samples through scaling laws [18, 54]. These laws are unable to, respectively, take into account the variations of waviness and unevenness at different scales.Additionally, our findings based on three natural rock joints sized1000mm×1000mm show that the fractal dimensions of waviness and unevenness seemingly vary without a universal trend as the joint sample size changes. In other words, waviness and unevenness of rock joints should be quantified at the scales of rock masses in the field mainly due to the random nature of asperity distribution along a joint surface [55].


## 6. Conclusions

We investigated the scale effect of surface roughness of natural rock joints using fractal approach. Three natural granite joints dimensioned of1000×1000mm were digitised and reconstructed at three different measurement resolutions. The fractal characteristics of two-order roughness, i.e., waviness and unevenness, were separately quantified through the classic triangular prism surface area method (TPM). We found that each-order roughness of a natural rock joint in window sizes varying from 100mm×100mm to 1000mm×1000mm owns individual fractal dimension. Although both the fractal dimensions of waviness and unevenness are scale-dependent, no noticeable stationarity threshold of scale effect has been found primarily due to the randomness of roughness distribution. Additionally, the measurement resolution has remarkable influence on the fractal dimension of unevenness, whereas its effect on that of waviness is negligible. Surface roughness quantification plays a key role in predicting the hydromechanical behaviour of rock masses. Our findings suggest that waviness and unevenness should be separately characterised at the field scale of the rock mass with an appropriate consistent measurement resolution. The conclusions are drawn by examining three natural rock joints with the dimension up to 1000mm×1000mm. The existence of a stationarity threshold of a larger value remains questionable due to the absence of experimental data. One may argue that the surface roughness of a natural rock joint likely exhibits more than two-order roughness.However, roughness characterisation with two-order roughness is sufficient for the purpose of accurately estimating the mechanical and hydraulic properties of rock joints. Further studies are to investigate the anisotropy of the fractal features of surface roughness since the hydromechanical behaviours of rock joints are strongly direction-dependent.



---

*Source: 8818815-2020-08-25.xml*