Characterization of rock joint surface anisotropy considering the contribution ratios of undulations in different directions

Anisotropy in rock joint is strongly dependent on undulating surface morphology. Recent research of the morphology showed the parameter can express the different types of anisotropic characteristics of the joint surface separately. This report aims to analyze the common characteristic of the anisotropic distribution and exhibit the anisotropic variation trend. The joint morphology function consists of two morphology functions of regular plane in orthogonal directions, and the anisotropic variation determined by the contribution ratios of the two morphology. The roughness weight ratio in orthogonal direction of joint surface is used as an index to describe the anisotropic variation behavior, which proposes the anisotropic variation coefficient (AVC). On this basis, it is divided into 5 levels from strong anisotropic to isotropic. According to the assumption of anisotropic arc distribution, the anisotropic analytic function is derived and the agreement between the deduced curves and measured data therefore suggests the possibility of defining the morphology anisotropy through the index AVC. Finally, we verify the characteristic of three natural rock joints, and prove the proposed function can reflect the anisotropic distribution trend. The new index can be used to describe the anisotropic variation behaviour of rock joint surfaces.


Scientific Reports
| (2020) 10:17117 | https://doi.org/10.1038/s41598-020-74229-z www.nature.com/scientificreports/ estimated. Therefore, a thorough investigation in the commonality of anisotropic characteristics between different joint surfaces needs to be carried out. Quantitative characterization of anisotropy is part of safety assessments for rock engineering projects. Based on the study by Belem et al. [21][22][23] , the degree of anisotropy can be categorised by the ratio of the small half-axis and the large half-axis of the apparent anisotropy ellipse into five types: (i) surface morphology is anisotropic; (ii) surface is more anisotropic than isotropic; (iii) surface is more isotropic than anisotropic; (iv) surface morphology is considered as isotropic; (v) surface morphology is isotropic. The anisotropy is simplified to the ratio of the maximum directional roughness to the minimum directional roughness (2D or 3D) in any direction 13 . Unfortunately, this ratio has ignored the influence of surface morphology in other directions on the anisotropy. To clearly show the degree of anisotropy, the factor of the composition of joint surfaces should also be considered.
In this study, we propose a new method for constructing joint surfaces with contribution ratio based on the investigation of the anisotropic distribution in existing researches. Furthermore, a new index anisotropic variation coefficient (AVC) is defined for quantitatively describing the degree of variation in anisotropy. Moreover, according to the behaviour of anisotropy comprehensively described by the proposed method, an anisotropic analytic function is derived for analysing the variation of anisotropy. Finally, the proposed AVC and anisotropic analytic function are applied in natural rock joints, and the accuracy is proved. In doing so, the anisotropic variation behaviour of rock joint surfaces might be effectively assessed.

Methodology
Investigation of the anisotropy of rock joint surfaces. The anisotropy of a joint surface is illustrated by plotting the roughness parameters obtained in different directions on polar coordinates [24][25][26][27] . Theoretically, the shape of the anisotropy of the joint surface distributed on the polar diagram is a circle when the roughness parameters in all directions are the same, and the shape becomes an irregular arc when the roughness parameter changes to one or some directions. As shown in Supplementary Table S2 online, the distribution variation, which is controlled by the values of the roughness parameters and its direction, is clearly reflected in the anisotropic distribution maps. Although the selected roughness parameters and their dimension descriptions are different, there still is a specific variation trend, that all distribution features of anisotropy are basically orthogonal. In order to analyse this trend, we normalize the anisotropic distribution maps in Supplementary Table S2 online, and plot them in a polar diagram which have uniform orthogonal direction (see Supplementary Fig. S1 online). It is shown the distribution shape present "∞" when the roughness in 90 • −270 • direction is zero. It means that the joint surface undulates only towards the 0 • −180 • direction. To better distinguish these two orthogonal directions, the 0 • −180 • direction with greater roughness is defined as the dominant direction, and the corresponding 90 • −270 • direction is the disadvantaged direction. when the roughness in disadvantaged direction gradually increases, the distribution shape like a "peanut shell" gradually expanded. Finally, the distribution changes to an "O" shape when the roughness is similar in both directions, and the joint surface is almost isotropic.
This phenomenon indicates that the complex joint surface morphology affected by geological movements in different directions may be explained by the superposition of two undulating surfaces in orthogonal directions, and the anisotropic characteristics can be approximatively determined by the morphological changes of the joint surface in these two directions. Therefore, we begin to study the variation characteristics of anisotropy by constructing the joint surface in orthogonal directions, which are helpful for determining the common of the anisotropic distribution in rock joints.

Construction of anisotropic joint surface.
To determine the common of the anisotropic distribution, we need to construct the rock joints with different undulating surfaces. Due to the complication of actual joint surfaces, it is necessary to consider the regular plane that undulating in one direction. As a special case in the study of joint surface, its morphology can be considered to be extended by the same profiles along the undulation direction (Fig. 1a). Thus, we first need to obtain the morphological function of a random profile. Yang et al. proposed the morphology of a joint profile can be approximated satisfactorily by using summation of several www.nature.com/scientificreports/ simple sine (or cosine) waves (Fig. 1b), and its morphological function z(x) varies periodically in x is defined with the magnitude a, the frequency b and the phase c as follows 28 : Based on the above discussion, the three-dimensional description of joint surface is carried out by extending the two-dimensional morphological function of the x-axis direction in the coordinate system by n units along the y-axis, and its three-dimensional morphological function Z(x) is given by Eq. (2). Similar result of Z(y) is also obtained and given by Eq. (3). In Fig. 2, both images are plotted as a function of x and y.
The actual rock joint is not regular as the above surface morphology, and thus we cannot directly establish the function of spatial position. However, the commonality between these regular or irregular planes is their anisotropy. According to the analysis of the anisotropic distribution in "Methodology" section, the actual joint surface can be considered as formed by the superposition of the same regular planes in orthogonal directions, and its morphology depends on the contributions of these two directions to the undulation. Consequently, an irregular plane can be mathematically described by the two deterministic components of regular plane in orthogonal directions, and the morphological function is given by: Here Z x and Z y represent the morphological functions of the regular plane in the x and y directions, respectively; k 1 and k 2 represent the contribution ratios of the x and y directions to the surface undulation, respectively. When k 1 = 1 or k 1 = 0 , the joint surface exhibits a regular undulating morphology; when 0 < k 1 < 1, the joint surface exhibits an irregular undulating morphology (Fig. 3). Therefore, the joint surfaces with different undulating morphology can be characterized by the contribution ratios.

Description of anisotropic variation behavior.
To identify the anisotropic variation behaviour, the contribution ratio is first selected to characterize the joint surface with different undulating morphology. Then, by means of the roughness parameters, we describe the morphological characteristics in different orientations. Finally, the anisotropic distributions are drawn on polar coordinates. Considering that the k 1 has the same contribution relationship in the ranges of [0, 0.5] and [0.5, 1] , except the undulating directions. Moreover, the x and y directions defined in this study are relative because the undulations of joint surfaces in these two directions are the same. Therefore, when k 1 ∈ [0, 0.5] , it can fully characterize the degree of undulating contributions in different directions. In the process of selecting the roughness parameters, we determine to use the directional parameters to better reflect the directionality of the morphology. In addition, considering the influence of the 2D and 3D parameters on anisotropy, we choose θ * max (C + 1) 2D and θ * max (C + 1) to measure morphological characteristics for more comprehensively analysing the anisotropic variation behaviour of joint surfaces 29,30 . Their calculation are as follows: www.nature.com/scientificreports/ where L 0 is the maximum potential contact length ratio, A 0 is the maximum potential contact area ratio, θ * max is the maximum apparent dip angle in the shear direction, and C is the roughness fitting coefficient.
According to the proposed construction method, MATLAB is used to construct six different joint surfaces whose k 1 are 0, 0.1, 0.2, 0.3, 0.4 and 0.5 respectively, with each size being 100 mm × 100 mm and height variation being less than 10 mm. Then, by cropping the square sampling windows into circles, eighteen 2D profiles of the same nominal length are extracted at angular increments of 10° (see Supplementary Fig. S2 online), and the 2D parameter θ * max (C + 1) 2D and 3D parameter θ * max (C + 1) are used to measure the morphological characteristics. The anisotropic distribution from Table 1 shows that when the k 1 ∈ [0,0.5], the behaviours of anisotropy reflected by different directional parameters agree with the approximate results, in which the variations are all from anisotropic to isotropic. Among them, the anisotropic distribution reflected by 2D parameter has obvious fluctuation compared to 3D parameter. This is because the local geometry of the joint surface described by 2D parameter ignores the continuity of the morphology between different directions, which indicates that the anisotropy described by 3D parameter is more accurate.

Analysis
Index AVC for describing the variation of anisotropy. It is verified from the results in "Methodology" section that the anisotropic variation is determined by the joint morphology in orthogonal directions. Based on this, we define an index anisotropic variation coefficient (AVC), which quantifies the degree of anisotropic variation by calculating the weight ratio of roughness in the orthogonal direction. For one of the orthogonal directions, however, the joint surface has different roughness in its forward and reverse direction. Take the direction of 0 • −180 • as an example, the roughness parameter in the forward direction ( 0 • → 180 • ) is not same as the reverse direction ( 180 • → 0 • ). So, the sum of those two roughness parameter values is taken as the weight of roughness to anisotropy in this direction. Then the AVC is determined as follows: where the P 90 and P 270 represent the forward and reverse roughness parameter in the direction of 90 • −270 • , P 0 and P 180 represent the forward and reverse roughness parameter in the direction of 0 • −180 • .
According to the definition of AVC, we calculate the degree of anisotropic variation at different contribution ratios in Table 1. The results show that the smaller the AVC is, the more prominent the anisotropic variation behaviour is. To better quantify the degree of anisotropic variation, we refer to the anisotropic classification proposed by Belem et al. 21 and provide a new anisotropic classification based on AVC which are shown in Table 2.
Anisotropic analytic function. To explain the anisotropic variation, we assume that the values of roughness parameters in different directions fall on an arc. Therefore, the anisotropic distribution can be interpreted by the AVC as follows: When the AVC is close to 0, its distribution shape is approximately tangent to two circles, indicating that the joint surface undulates in only one direction (Fig. 4a); when the AVC in the range of (0,1), its distribution shape approximately changes to two intersecting circles. And the undulation of the vertical joint surface gradually increases as the AVC increases (Fig. 4b,c); when the AVC is close to 1, its distribution shape changes to a circular, indicating that the roughness parameters in each direction are similar (Fig. 4d). This phenomenon is closely related to the roughness parameters in the orthogonal directions. Therefore, we divide the anisotropic distribution into four parts and select arc segments with the range of (0, π/2) for anisotropic variation analysis (the red parts in Fig. 4). www.nature.com/scientificreports/ According to the basic equation of the circular curve, we set ∆x as the distance from the centre of the circle to the origin of the coordinates, the value of P 0 and P 90 are the polar radiuses which correspond to the polar angle θ are 0° and 90°, respectively (see Supplementary Fig. S3 online). Then, the roughness parameter P (ρ, θ) for any direction in the range of (0°, 90°) can be calculated by the following relationship:  www.nature.com/scientificreports/ let K 1 = P 90 P 0 , then, ∆x can be expressed as: 2 and substitute into Eq. (6), we obtain the following relationship: The above function explains the anisotropic variation of one arc segment in the range of 0, π 2 . Hence, the anisotropic distribution function, which can be used to explain the entire anisotropic distribution, is obtained by calculating each arc segment and given by: where P is the value of the roughness parameter in the direction of 0 • −180 • . When θ ∈ (0, π/2) ∪ 3π 2, 2π , the value of P is P 0 . When θ ∈ π 2, 3π 2 , the value of P is P 180 ; K ρ represents the calculated value related to the K 1 in the corresponding angle range.
Using the above function, the anisotropic reanalysis of the above six sets of joint surfaces with different contribution ratios is performed, and the anisotropic distributions are listed in Table 1 (the red curves). As a comparison of prediction effects, Fig. 5 gives an error analysis of the deduced curves and measured data, and the average relative error − δ is calculated as follows: where M represents the number of roughness parameters, S is the measured roughness parameters, and L is the values of deduced curves. The results show that the deduced curves are almost identical to the measured data and the average relative errors between the two types of values range from a minimum of 3% to a maximum of 8%.

Application and discussion
Actual joint surface measurements. To verify the anisotropic variation behaviour of actual joint surfaces, we first need to collect natural joints. After field investigations, there are three open-pit mines with different lithology (sandstone, tuff, and limestone) in Shaoxing City, Zhejiang Province, China, and the exposed joints Table 2. Classification of anisotropic variation levels. www.nature.com/scientificreports/ are well preserved (see Supplementary Fig. S4 online). Then we select a set of 1000 mm × 1000 mm joints from each mine as the research object and use a portable laser scanner (MetraSCAN 3D, Creaform, Canada) with the scanning accuracy of 0.5 mm to obtain the 3D morphological information of its surface. Finally, three digitized surfaces labeled as S1, T1 and L1 are obtained, which are shown in Supplementary Fig. S5 online.

Range of the variation in AVC
Description and analysis of anisotropy. The original anisotropic behaviours are obtained by using the 3D parameter θ * max (C + 1) . However, due to the complicated distribution in actual joint surfaces, it is difficult to distinguish the locations of the orthogonal directions. To this end, we have developed an adaptive algorithm that can automatically match the roughness parameters of the orthogonal direction for anisotropic analysis. It mainly implements by using the anisotropic analytic function to obtain the deduced curves in each vertical direction on the joint surface, and then compares with the measured data. Among them, the vertical direction with minimum relative error is the location of orthogonal directions. This algorithm is applied to the above three joints, and the calculation results listed in Table 3 indicates that there are obvious orthogonal distribution characteristics in the S1, T1 and L1 joints. In addition, the relative errors of actual anisotropic distribution and analytical anisotropic distribution stable within 10%. Therefore, the anisotropic analytic function can effectively reflect the anisotropic variation behaviour of rock joints.
Extraction and statistics of index AVC. From the anisotropy distribution of actual joints, the orthogonal direction is not specified in the directions of 0 • −180 • and 90 • −270 • . To better distinguish the orthogonal direction, the direction with smaller weight of roughness is determined as the inferior direction, and the forward and reverse roughness parameters in this direction are defined as P F Inf and P R Inf , respectively. On the contrary, the direction with greater weight of roughness is determined as the advantageous direction, and the forward and reverse roughness parameters in this direction are defined as P F Adv and P R Adv , respectively. Finally, the AVC can be further expressed as follows: Therefore, we calculate the AVC of each joint according to the above formula. The results show that the AVC values of S1, T1 and L1 joints are 0.77, 0.63 and 0.98, respectively. Compared with the classification of anisotropic variation levels (Table 2), the S1 joint is presents weak anisotropic variation, the T1 joint presents medium anisotropic variation, and the L1 joint is close to isotropic distribution.

Conclusions
The method for superimposing regular surfaces in orthogonal directions is used to construct the different undulating morphology of the joint surfaces and to characterize the anisotropic behaviour, with index AVC being proposed to quantify the degree of anisotropic variation. Based on the assumption of anisotropic arc distribution, the anisotropic distribution function is derived and the capability of the anisotropic variation behaviour are also validated. Some concluding remarks are presented below: A specific distribution law is investigated in the anisotropy of the joint surface. Such a law shows that a rock joint surface can be simplified to the superposition of regular surfaces in orthogonal directions, and the anisotropy is described by the contribution ratios in these two directions. Six joint surfaces with different contribution ratios are constructed base on the proposed construction method of joint surfaces. The behaviours of anisotropy expressed by 2D and 3D parameters indicate that the contribution ratio can effectively reflect the process of joint surfaces from anisotropy to isotropy. Figure 5. Error analysis of deduced curves and measured data.
Scientific Reports | (2020) 10:17117 | https://doi.org/10.1038/s41598-020-74229-z www.nature.com/scientificreports/ The AVC is proposed by defining the roughness weight ratio in the orthogonal direction of joint surfaces. According to the size of AVC, we divide the level of anisotropic variation into five different intervals to describe the degree of anisotropy. Additionally, we derive an anisotropic analytic function based on the assumption of anisotropic arc distribution, and verify it can be used as an effective method for analyzing anisotropic characteristics.
By measuring the anisotropy of three natural rock joints, it is found that there are obvious orthogonal distribution characteristics. The anisotropic analytic function is used to determine the locations of orthogonal directions. The results show that the inferior directions of S1, T1 and L1 joint are 100 • −280 • , 140 • −320 • and 170 • −350 • , and their advantageous directions are 10 • −190 • , 50 • −230 • and 80 • −260 • , respectively. According to the extracted AVC value of each joint, the S1 joint is presents weak anisotropic variation, the T1 joint presents medium anisotropic variation, and the L1 joint is close to isotropic distribution.
In summary, this study provides a new method for accurately judging the strength and weakness of joint mechanical properties and a new index to quantify the degree of anisotropy, which can be used to evaluate the stability of engineering rock mass. However, it is worth noting that AVC is proposed based on the anisotropic orthogonal distribution reflected by mathematical statistics method. Whether the other two methods (empirical value method and fractal dimension method) have the same distribution law and then apply AVC to determine the anisotropic variation behavior needs to be further verified. Additionally, the anisotropic distribution function is derived based on a circular curve equation and thus it should be further studied with an elliptical curve. Furthermore, the scale effect of AVC will be studied in the future. License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.