Abstract

Rock quality designation (RQD) is widely adopted as a fundamental tool in characterizing rock masses since it was devised in 1964. Since the conventional RQD calculation is limited due to its dependence on the selected threshold, previous research introduced generalized RQD to adequately reflect the anisotropy and scale effect of RQD. However, the influence of the joint development inside rock mass on generalized RQD remains unclear. The objective of this work is to investigate characteristics of the generalized RQD in view of different development degrees of discontinuities in rock mass, including spacing (density) and trace length. Three-dimensional fracture network modelling is employed to simulate the actual rock mass of open-pit iron mine in China. Virtual scanlines are set to obtain RQD values in different directions. The results primarily show that the generalized RQD should be introduced to calculate the RQD with different thresholds to fully reflect the anisotropy of rock mass. The optimal threshold can be obtained based on an anisotropic coefficient, which is defined by (RQDmax-RQDmin). It is also indicated that the fracture spacing has a great influence on both the anisotropy of RQD and the selection of the optimal threshold. The optimal threshold of the generalized RQD increases with the increase in the fracture spacing. In addition, the scale effect of RQD in different models is discussed by changing the length of the scanlines. The longer the scanlines we set, the more stable RQD value can be obtained in the model. It is recommended to fit much longer scanline to get realistic RQD in heavily fractured rock mass.

1. Introduction

The presence of discontinuities plays an essential role in the strength and stability of rock masses, which controls the damage characteristics of rock mass [14]. Rock quality designation (RQD) is an important parameter in evaluating the quality of rock masses using drill core cataloging. This parameter was initially proposed by [5]. It provides a quite simple description of rock mass quality qualitatively [6]. At present, it has been widely employed as the foundational parameter in the rock mass classification system such as Rock Mass Rating (RMR) and the Rock Quality Tunneling Index (Q-system) [711]. RQD is internationally adopted in areas of mines, rock mechanics, and engineering geology [12, 13]. The value of RQD is expressed as follows [14], where is all pieces of core with length greater than 100 mm and X is the total length of the drilling borehole.

However, three limitations in the aforementioned calculation principles are observed when obtaining RQD. First, the “threshold value of 10 cm” is used as the only criterion for judging rock quality. Whether this threshold can be used as the most reasonable threshold value for distinguishing rock mass has been questioned by many scholars [1519]. Second, due to the geological survey level and engineering cost, the value of RQD can be different for a given location when measured from cores with different drilling orientations [20, 21]. The conventional RQD cannot fully reflect the heterogeneity and anisotropy of the rock mass [22, 23]. Third, traditional methods often directly calculate the size of RQD based on the length of each section of the core in the borehole, ignoring the scale effect of RQD [2426].

To overcome the aforementioned limitations, many domestic and foreign scholars have conducted relevant research. The authors in [27] found that RQD changes with the size of the scanline; that is, RQD has scale effect. Thus, the authors proposed corresponding empirical formulas as follows:where λ is the density of structural surface (average number of structural surfaces per meter of scanline) and t is the selected threshold for calculation. Harrison believed that a reasonable selection of thresholds can expand the range of RQD change and proposed a corresponding equation [16]. Based on the relationship between threshold and RQD change, the generalized RQD was suggested in light of rock mass anisotropy [28]. The authors in [29] investigated the generalized RQD through in situ cataloging of discontinuities and 2D network simulation technology. The results pointed out the limitations in using 0.1 m as the only threshold. The authors in [30] conducted the three-dimensional fracturing network method to simulate rock slope of the Baihetan hydropower station. By arranging numerous scanlines, it is indicated that the generalized RQD has obvious spatial and scale effects.

The RQD index is influenced by different parameters which mainly include joint orientation, joint spacing, fracture frequency, and joint roughness [31]. Although predecessors have conducted extensive research on the limitations of the traditional method to calculate the generalized RQD, how the development of the discontinuities in rock mass will affect the characteristics of the generalized RQD remains unclear. Therefore, this study adopts 3D fracture network simulation [32, 33] to generate numerous models of an actual rock mass by considering different jointing degrees. A series of virtual scanlines were arranged in different directions in each model. The influence of joint development on the anisotropy of the generalized RQD is analyzed. The scale effect of generalized RQD is also studied by changing the length of the scanline, which provides a certain theoretical basis for selecting a reasonable drilling depth of rock projects.

2. Model Establishment of an Actual Rock Mass

2.1. Study Area and Geological Setting

Practical investigations on RQD measurements were carried out at Anqian open-pit iron mine, which is located in southwest of Anshan, in Liaoning Province, China. The strata of this area were precipitated in medium and high metamorphic series with NW-SE trend. The major rock types at the open-pit slope are ferruginous quartzite, schist, and phyllite of the Anshan Group, which belongs to Precambrian.

We adopted the noncontact sampling method to conduct geological logging at Xudonggou stope using the ShapeMetriX3D system [34]. Four key discontinuity parameters of 258 discontinuities, i.e., orientation, number of sets, set spacing, and mean trace length, were collected and analyzed. Three sets of discontinuities intersect the rock mass of the stope: a bedding and two major joints, according to the variation range of dip angle (see Figure 1).

2.2. Model Parameters

Gu [35] divided the joints of rock mass into five levels according to the scale of structures: I, II, III, IV, and V. Among them, the IV-grade structural planes extend from tens of centimeters to 20–30 m, with a large number distributed randomly, thereby affecting the integrity and mechanical properties of the rock mass. This issue is the subject of statistical analysis and simulation of structural planes. Studies have found that when the number of structural planes is constant, the development degree of discontinuities is mainly related to their spacing and ductility, and the occurrence and distribution type of structural planes have a minimal, even negligible, effect on it [36].

Regarding the quantitative description of structural plane spacing and ductility, the International Society for Rock Mechanics [37] proposed the ductility of structural planes into five levels (<1, 1–3, 3–10, 10–20, and >20; m) and the spacing into 7 levels (<20, 20–60, 60–200, 200–600, 600–2000, 2000–6000, and >6000; mm). Based on the preceding classification and the scale of the case study, this work only selects seven representative structural plane spacing values and five representative structural ductility ranges, which are tabulated in Table 1.

The discontinuities inside the rock mass are mainly divided into occurrence, ductility, and density (spacing). The accuracy of the above elements greatly determines the accuracy of the 3D joint network model [38]. Three sets of discontinuities at Xudonggou stope are developed in the rock mass. The distribution of occurrence follows the Fisher distribution [39], which contains only one parameter, the Fisher constant K. The spacing and trace length distribution obey the negative exponential distribution [40]. The specific parameters of 3D fracture network simulation are listed in Table 2. According to the result of [41], the value of CP3 is 2.0 for a uniform distribution of fracture orientation. For most fracture geometries, CP3 will vary between 1.0 and 3.0. Therefore, the value of CP3 is set to 2.0 here. Meanwhile, the density of the discontinuities can be calculated from the spacing presented in Table 1. The solution method is based on the work of [42]; and the calculation process is not mentioned in this study.

2.3. Model Generation and Scanline Layout

The existence of discontinuities causes the rock mass to have an obvious scale effect. The parameters of rock mass change with the sampling scale [43]. With the increase in the sample to a certain critical value, the parameters of rock mass tend to be stable. The sampling size corresponding to the critical value is the representative elementary volume (REV). The authors in [44] introduced the index of the blockiness degree of fractured rock masses. The result indicated that the REV is 4–8 times of the spacing between discontinuities and no more than 8 times of the spacing. According to the selected maximum spacing of 6 m, the scale of the study area is set 50 m × 50 m × 50 m.

A total of 35 rock mass models with different structural plane development levels are generated using 3-Dimensional Distinct Element Codes (3DEC, Itasca Consulting [45]), as illustrated in Figure 2. To fully reflect the development of the discontinuities inside the rock mass, 5 cut planes are selected in each model. 18 scanlines are arranged on each cut plane at intervals of 10°. The directions of scanlines cover 0°–360°. Due to the symmetry of polar coordinates, we only calculate the scanlines in the range of 0°–180°. The selection of rock mass section and the layout of scanlines are displayed in Figure 3. According to the intersection of the internal structural plane of the model and the survey line, the program can automatically calculate the line density in the survey line direction. Combining the Priest–Hudson formula (equation (2), we can obtain the RQD value of the model in a certain direction.

3. Anisotropy of the Generalized RQD

Because of the discontinuous nature of rock masses, they are anisotropic, inhomogeneous [46, 47]. The values of RQD in different locations of the same rock mass should be inconsistent. That is, the distribution of RQD also has anisotropy. This section introduces the concept of the generalized RQD to analyze the anisotropy of RQD based on the degree of joints development.

3.1. Generalized RQD

Selecting a reasonable threshold can expand the variation range of RQD, which plays an important role in studying the heterogeneity and anisotropy of rock mass. Therefore, the concept of generalized RQD was introduced to maximize the variation range of the conventional RQD. For any threshold t, the percentage of the ratio of the sum of the length of the core length not less than t along a certain scanline (borehole) to the length of the scanline (borehole) is defined as the generalized RQD [48].

According to the arrangement in Figure 3, the scanlines of the selected section are set in each model with different thresholds. The average value can be calculated based on the measurement results from those scanlines with the same direction of each section. It is regarded as the RQD inside the rock mass in a certain direction. All calculation results are drawn into a 2D polar diagram. As shown in Figure 4, the angle in the figure indicates the direction of the scanline and the polar radius indicates the magnitude of RQD. Due to constraints on length, Figure 4 only lists the RQD polar coordinate plane diagrams of nine models composed of extremely dense, medium dense, and extremely wide spacing with extremely low, medium, and high ductility, respectively.

3.2. Determination of Optimal Threshold

Figure 4 describes the RQD anisotropy under different structural plane development levels and different threshold conditions. As the distance remains the same, the polar coordinate plot changes within a narrow range with the ductility changes from 1–3 m to 20–30 m. As the internal ductility remains the same, however, the polar coordinate plot changes significantly with the spacing changes from 0.2 m to 6 m. Hence, it is indicated that the RQD anisotropy is more sensitive to the spacing of the discontinuities. In addition, the polar coordinate plot in different models changes significantly with different thresholds, which shows that the threshold also has a greater impact on the RQD anisotropy. Therefore, the distance between the discontinuities within the model and the selected threshold are the main factors that affect the calculation results of RQD.

Since the anisotropy of the rock mass leads to the anisotropy of RQD, RQD anisotropy can reflect the anisotropy of rock mass. A reasonable threshold should be selected to maximize the variation range of RQD, which can reflect the anisotropy characteristics of rock mass. In this study, we define the optimal threshold of RQD as it can make the degree of RQD anisotropy reach the maximum in the model. Taking the model with a spacing of 1 m and a ductility of 8–14 m as an example, we set the thresholds from 0.1 m to 1.0 m. The degree of the RQD anisotropy can be characterized as the range of RQD (RQDmax-RQDmin). If the range of RQD reaches the maximum, the corresponding threshold can be considered as the optimal threshold.

It is illustrated in Figure 5 that the RQD range increases first with the increase in the threshold and then decreases when reaches a certain critical value. In this model, the RQD range reaches the maximum with the threshold of 0.7 m, which means the optimal threshold is 0.7 m. Therefore, the conventional RQD owes obvious limitation because the calculation with the only threshold of 0.1 m cannot adequately characterize the anisotropy of rock mass.

Furthermore, the optimal thresholds are not the same in different models. The degree of the discontinuity’s development will affect the selection of the optimal threshold. Figure 6 presents the relationship between the optimal threshold and the development of the discontinuities (including spacing and ductility) among 35 different models. It shows that the optimal threshold of RQD increases with the increase in the spacing between the discontinuities, but the ductility has no significant effect on it. With the distance of 0.2–0.5 m, the threshold does not change with the ductility of the discontinuity. As the distance increases to 0.8–1.5 m, the optimal threshold changes slightly with the ductility. When the distance reaches 3–6 m, the optimal threshold starts to fluctuate with the change of the ductility.

Therefore, the influence of the ductility of the discontinuities on the optimal threshold of RQD is mainly related to the spacing of the discontinuities, which becomes remarkable with the increase in the spacing. Consequently, the distance between the discontinuities is the major factor for both the selection of the optimal threshold and the anisotropy of RQD.

4. Scale Effect of the Generalized RQD

Engineers usually calculate the RQD according to the length of the core in the borehole, ignoring the impact of the drilling depth on the RQD value. In fact, the true RQD of rock mass at any drilling depth (the length of the scanline) cannot be obtained [49, 50]. Thus, the issue is discussed using 3D fracture network simulation. The scale effect of the generalized RQD can be analyzed by changing the length of the scanlines.

Generally, engineering geological drillings are vertical because of the geological environment and drilling technology. Therefore, the scanlines are arranged as boreholes in the generated rock mass model, which are perpendicular to the xoy plane. A total of 121 scanlines are shown in Figure 7. The scanline is arranged every 5 m to get enough data.

Priest and Hudson [27] suggested a theoretical equation related to the calculation of the generalized RQD (equation (3)), with respect to the hypothesis of negative exponential distribution of joint spacing. Equation (2) was proposed with regard to the effect of the scanline length. Another possible more realistic expression for the generalized RQD was also suggested by [51]. Compared with equation (3), the calculation of equation (4) is concerned more about the influence of short scanline on the RQD result. In these equations, λ is the density of the discontinuities, t is the calculation threshold, and L is the length of the scanline.

Five typical models with intervals of 0.2, 0.5, 0.8, and 1 m were selected as examples. The values of the generalized RQD with different lengths of the scanline are calculated according to equation (4) and the optimal thresholds of each model. The variation curves of the generalized RQD are plotted in Figure 8. It is assumed that the true RQD can be obtained under the full length of scanline, i.e., 50 m in this study. Then, the rate of RQD change is brought forward to compare the generalized RQD under each scanline length in different models with that corresponding to the scanline length of 50 m. The relationship between the rate of RQD change and the length of the scanline is shown in Figure 9.

The various models indicate that the generalized RQD converges with the increase in the scanline length, which has an obvious scale effect. The RQD result can represent the rock mass quality more accurately under the longer scanline. However, the convergence of the generalized RQD in different models is not consistent. When the spacing is 0.2 m, the RQD changes from large to small and then to converge as the length of the scanline increases. As the distance between the discontinuities increases, the change of RQD increases first and then decreases to convergence with the length of the scanline. The ductility of the discontinuities also affects the convergence of the RQD. The rate of the RQD change in each model increases as the ductility increases. It is indicated that the convergence length of scanline is positively associated with the ductility of the discontinuities.

For those rock masses without a recognizable structural pattern, it is difficult to determine the real RQD. When the number of scanlines is large enough, the average value of RQD can be close to the real RQD of rock mass infinitely. Therefore, the average RQD obtained from 121 scanlines with full length (50 m in this study) is assumed to be the real RQD. Subsequently, the average value obtained from different lengths of scanlines, which is recorded as RQD#, is regarded as the undetermined RQD of the rock mass. This study defines the variation rate of RQD as the difference between RQD# and RQD. With respect to the engineering requirement, it is considered to meet the research requirements when the variation rate of RQD is not greater than 2%, which is described in equation (5). We defined the RQD# in a certain scanline satisfied the inequality as the convergence RQD value, and the corresponding scanline length is the convergence length of the scanline.

Table 3 tabulates the shortest convergence length of the scanline and the convergence RQD value that meet equation (5) requirements of each model. The degree of joint development distinctly affects the convergence length of the scanline. The convergence length generally exhibits a slight increase with the increase in the spacing. However, the increase is not uniform in each model. Compared with the spacing, the ductility greatly influences the convergence length of the scanline. With the same spacing, the convergence length in different models increases with the increase in ductility.

According to the combined results of Figures 7 and 8 and Table 3, the generalized RQD has an obvious scale effect. The RQD result is more stable with respect to the longer measuring line, and it can represent the quality of the entire rock mass more accurately. In addition, a longer scanline is needed to obtain a reliable RQD value for a heavily fractured rock mass.

5. Conclusions

This study used 3D structural plane network simulation technology to generate 35 rock mass models with different structural plane spacing and ductility, arranged survey lines inside the model instead of drilling holes, obtained enough values of the generalized RQD, and studied the relationship between the generalized RQD characteristics and the development of discontinuities. The specific conclusions are as follows:(1)The traditional method using 0.1 m as the only threshold for RQD has significant limitations. It cannot fully reflect the anisotropy and scale effect of the rock mass. The generalized RQD concept is introduced, and the optimal threshold is obtained based on the degree of RQD anisotropy. The results show that the degree of joint development affects the selection of the optimal threshold. The optimal threshold increases with the increase in the spacing, while of the discontinuities have a minimal influence on the selection of the optimal threshold. However, the influence degree of the ductility shows an increasing trend with the increase in the spacing.(2)The degree of joint development has a significant effect on RQD anisotropy [52]. When the spacing is larger, the rock mass has higher anisotropy, while the ductility has a minimal effect on RQD anisotropy, which can even be ignored. Therefore, if the discontinuities with a larger spacing are encountered in the practical engineering, the number of drillings should be increased so that the RQD calculation can fully reflect the anisotropy of the rock mass.(3)The generalized RQD has an obvious scale effect. The longer the scanline, the more stable the RQD obtained, and the more accurately it can represent the quality of the rock mass. In addition, when the discontinuities are highly developed in the rock mass, the scale effect of RQD is more obvious, and a longer scanline is needed to obtain a real RQD. Therefore, the sufficient length of the scanline is necessary to obtain stable and reliable RQD. Particularly, in heavily fractured rock masses, the length of the scanline should be appropriately increased to ensure the accuracy of RQD [53].

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 in the publication of this article.

Acknowledgments

This study was partially funded by the Zhejiang Provincial Natural Science Foundation of China (Grant no. LQ18D020001) and National Natural Science Foundation of China (Grant no. 42007254), and this support is gratefully acknowledged.