Theoretical Study on Radial Distribution Laws of Rock Mass Damage Factors under Decoupled Charge Blasting

In this paper, the radial distribution laws of damage factors under decoupled charge blasting are studied for the optimization design of blasting parameters. -rough defining the critical radial decoupling coefficient, the damage zone around the borehole is partitioned and the characteristics are described. Based on the damage factor defined by Taylor’s effective elastic modulus, the formulas of the radial distribution laws of damage factors are derived by the attenuation law of stress wave and the theory of thickwalled cylinder, respectively, which are then superposed to obtain the formula under the combined action of explosion stress wave and quasistatic gas. Experimental verification indicates that the theoretical values, which have a good correlation with the test data and are of high accuracy, can characterize the radial distribution laws of damage factors and estimate the damage range. When a radial decoupling coefficient is less than the critical value, the attenuation rate of damage factors firstly increases and then decreases with the increase of distance, and a serious damage zone is caused. Conversely, it decreases gradually, and the serious damage zone is not caused.-erefore, on the premise of stable detonation, it is necessary to apply an appropriate radial decoupling coefficient which is larger than the critical value to smooth or presplit blasting.


Introduction
During the process of formation and evolution of rock mass, there are a lot of initial damages such as microcracks and microcavities in the interior [1,2]. Microcracks are activated, grown, and coalesced gradually to form macrocracks under the load, which finally causes the physical and mechanical properties of rock mass to be lowered and even be failed [3]. Hence, rock mass failure is a process of the accumulation and evolution of damage [4]. Blasting action can inevitably cause different degrees and ranges of damage to the reserved rock mass. Due to the irreversible process of damage, the gradual accumulation and evolution cause rock mass to be fractured, which endangers the integrity and stability of rock mass. erefore, the study on blasting damage is always one of the primary research topics in the field of rock engineering [5], especially damage control in the bedrock blasting excavation of an underground radioactive waste repository and a nuclear power plant [6,7].
At present, studies on blasting damage in rock mass have mainly focused on constitutive models [4,[8][9][10], evolution mechanisms and characteristics [11][12][13][14], tests and safety assessments [15][16][17][18][19], control technologies [20][21][22][23], and so on. Among them, the damage range and degree are two extremely important contents of damage characteristics. However, up to now, relatively few studies have been done on the radial distribution laws of the damage factors within the blasting damage range. Carrying out the study can not only identify the radial distribution laws of damage factors and the interrelation between damage factors and blasting parameters but also estimate the damage range and further control damage.
Some scholars have studied the radial distribution laws of the blasting damage factors. Dai et al. [24] found that the damage factors around a single hole decreased in an approximate power function law with the increase of distance and fitted a relation between damage factors and scaled distance through the experimental data. Suo [25], through ultradynamic strain tests of coal samples, put forward that the volumetric strain accorded with the attenuation law of a power function and applied Taylor's and Grady's damage models to derive the distribution function of the macroscopic damage degree generated by preweakening the blasting of coal seams. Tang et al. [26] derived the distribution laws of damage factors based on the relation that elastic modulus was proportional to transmission frequency of stress wave as well as the relation of blasting dominant frequency to charge mass and to the distance from the explosion center. Meanwhile, five damage zones were partitioned. Pan et al. [27] found that damage factors presented an attenuation law of exponential function with distance and fitted out related parameters through numerical simulation. Although the abovementioned research achievements enrich the theory of blasting damage to some extent, there are the following shortages to be further improved. Empirical formulas fitted out by model experiments and numerical experiments still need theoretical supports. e existing theoretical analysis only considers the contributions of explosion stress wave to damage evolution, without considering the contributions of the quasistatic gas. When the strain analysis is applied, only the volumetric strain is considered, and the equivalent strain is ignored. erefore, it is necessary to conduct a theoretical analysis on the radial distribution laws of damage factors under the combined action of explosion stress wave and quasistatic gas, which can provide a theoretical guidance for optimizing blasting parameters and carrying out model experiments and numerical experiments.
Based on the effective elastic modulus of mesodamage mechanics, the radial distribution laws of damage factors under the combined action of explosion stress wave and quasistatic gas are derived by the attenuation law of explosion stress wave and the change law of elastic stress field in the thick-walled cylinder. e correctness and reliability is validated by experimental data so as to lay a theoretical foundation for the optimization design of blasting parameters.

Definition of Damage Factor.
At present, the definition of damage factor has not had a unified form. Applying macroscopic material parameters to characterize damage degree has been the common method. According to macroscopic material parameters, the damage factor of rock mass, D, can be expressed as where E is the elasticity modulus before blasting, c p is the wave velocity before blasting, E is the equivalent elasticity modulus after blasting, c p is the wave velocity after blasting, and η is the reduction rate of the wave velocity.
In the strict sense, D > 0 represents the fact that rock mass is damaged. However, considering the influence of blasting damage on the properties of rock mass, the code SL47-94 stipulates that rock mass can be affected by blasting damage when η > 10% and the corresponding damage threshold is 0.19.

Partition of Damage
Zone. Generally, due to strong impact, rock mass blasting with the coupled charge or smaller decoupling coefficient charge forms the crushed zone, the fracture zone, and the vibration zone with the borehole as the center. In engineering, the decoupled charge is mostly used to weaken shock effect and eliminate the crushed zone so as to reduce the blasting effect on the reserved rock mass. at is to say, the purpose is to avoid producing the serious damage zone of "cracking without breaking" from the perspective of damage. e decoupled charge structure usually contains the radial and axial decoupled charge structure, so applying the volumetric decoupling coefficient can describe the radial and axial decoupling coefficients. e smooth blasting technology is used in the perimeter controlled blasting. e radial and axial decoupled charge structure is usually adopted in the borehole. Particularly, it is more likely to form local fracture and damage around the charge segment. It is obvious that the influence of radial decoupling coefficients around the charge segment on the radial distribution laws of damage factors should be studied in advance or emphatically under a certain volumetric decoupling coefficient. For this purpose, the maximum radial decoupling coefficient of "cracking without breaking" caused by rock mass blasting, K 0 , is defined as the critical radial decoupling coefficient whose value is related with properties of rock mass and explosive. When the radial decoupling coefficient is greater than K 0 , namely, K ≥ K 0 , rock mass around the borehole can be partitioned into three zones: the general damage zone, the slight damage zone, and the undamaged zone in turn, as shown in Figure 1(a). When K < K 0 , rock mass around the borehole can be partitioned into four zones: the serious damage zone, the general damage zone, the slight damage zone, and the undamaged zone in turn, as shown in Figure 1 e abovementioned rock mass damage in each zone, which represents growth or coalescence of microcracks, results from the combined action of explosion stress wave and quasistatic gas. e characteristics of rock mass damage in each zone can be defined. In the serious damage zone, rock mass discontinuities and rock blocks are all fractured, but do not fully form intersecting blocks; that is, blasting cracks in rock mass have not been fully coalesced. In the general damage zone, rock mass discontinuities are fractured while rock blocks are not; that is, microcracks in rock blocks have not been grown to form blasting cracks. In the slight damage zone, rock mass discontinuities are locally fractured, but growth of fractured cracks is extremely limited and microcracks in rock blocks are not grown. In the undamaged zone, pressure decreases with the increase of distance from the borehole and is not enough to make microcracks grow, so rock mass is not damaged.

Analysis and Assumptions.
In general, there exist two different actions under the blasting load, that is, the action of explosion stress wave and the action of quasistatic gas, and both cause different degrees of damage to rock mass. Although the radial and axial decoupled charge structure is used in smooth or presplit blasting, which can weaken the explosion impact greatly, damage caused by the action of explosion stress wave can provide favorable conditions for the action of quasistatic gas so that damage caused by the action of quasistatic gas cannot be ignored. erefore, it is necessary to study the radial distribution laws of damage factors under the combined action of the two. e abovementioned action processes are sequential but complex, so it is difficult to separate them completely. To simplify the analysis, when deriving the analytical solution of damage factors under the combined action of the two, the action of explosion stress wave and quasistatic gas should be considered separately in advance and then be superposed together. Besides, the blasting damage is a function of time and space and it is complex to consider the continuous time effect when studying the radial distribution laws of damage factors. erefore, the distribution laws of the maximum damage factor in two discrete time points selected from two periods of explosion stress wave and quasistatic gas are only focused on. Pressure of time points acted by explosion stress wave takes the stress peak value and that acted by quasistatic gas takes the quasistatic pressure value.
In smooth or presplit blasting, the cracks grow directionally along the connecting line of the borehole centers and the growth of cracks in surrounding rock is restrained due to the stress concentration effect between the adjacent holes. e stress concentration zone along the connecting line of the borehole centers is relatively smaller than the damage zone of surrounding rock, so the effect on damage of surrounding rock can be ignored. In addition, the damage characteristics in surrounding rock are generally focused on in smooth blasting. e stress concentration effect under the influence of adjacent holes is mainly on the connecting line of the borehole centers, while the influence of the stress concentration effect on surrounding rock is less, especially on the right back of boreholes. erefore, smooth or presplit blasting can be simplified as single-hole blasting for analytical analysis.
It is easy to obtain initial parameters of homogeneous material before blasting, while the change laws of material parameters after blasting is unclear. And for a mass point, the wave velocity changes with time after blasting. Because the propagation velocity of stress wave is faster than the growth velocity of microcracks and wave velocity in the equation σ � ρ r cυ is still the original wave velocity before blasting, c p , it is difficult to derive the relation between wave velocity and stress after blasting by using this equation. e essence of material damage is nucleation, growth, and connection of microcracks, while the mechanical description of microcracks is exactly made by mesodamage mechanics to reveal the essence of material damage. On the one hand, mesodamage mechanics ignores the microscopic physical process in which damage is too complex and avoids the heavy calculation of statistical mechanics. On the other hand, it contains geometrical and physical characteristics of mesodamage of different materials and provides a clear physical background for the damage variables and the evolution process of damage. Besides, the average idea is adopted in mesodamage mechanics and the results of study on damage mechanism of mesoscopic structure are reflected on the behavior description of macromechanics [28,29]. Shock and Vibration erefore, the interrelation between microcracks and the radial distribution laws of damage factors can be communicated from the perspective of effective elastic modulus of mesodamage mechanics. e explosive is usually in the center of a borehole. Detonation gas expands around the explosive with the same high speed and acts on the borehole wall. In addition, rock mass is generally in the state of initial ground stress which is very small relative to the blasting load. Meanwhile, rock mass under the blasting load shows a quasibrittle characteristic. For smooth blasting holes, the amplitude value of explosive stress around the charge segment is the largest and reduces gradually in the process of the axial transmission to the air segment. When the borehole is filled with detonation gas, quasistatic pressure is finally formed so as to make damage around the charge segment greater than that around the air interval segment [24]. erefore, from the most unfavorable perspective, the radial distribution laws of damage factors within the range of the charge segment under decoupled charge blasting are only considered when the linear charge density is fixed. To highlight the focus and simplify the analysis, the following assumptions are made: (1) e influence of initial ground stress is not taken into account. (2) Distribution of borehole pressure is uniform.
(3) Rock mass is elastic material containing microcracks with uniformly random distribution. (4) e action of explosion stress wave and quasistatic gas are considered separately.

Derivation of eoretical Formulas.
According to the theory of mesodamage mechanics, the calculation methods of the effective elastic modulus mainly include Taylor's method without considering the interaction among microcracks, the methods considering the weak interaction among microcracks, such as the self-consistent method, the generalized self-consistent method, Mori-Tanaka's method, and the differential method, and the statistical mesomechanics method considering the strong interaction among microcracks. For most quasibrittle materials, the influence of the interaction among microcracks on the constitutive relation is not obvious, and meanwhile, there is no simple and effective method to calculate the interaction effect among microcracks [28,29]. erefore, this paper uses the relatively simple Taylor method to analyze rock mass with microcracks having the same size and uniformly random distribution of orientation and location. e effective elastic modulus of rock mass with microcracks is [8] where E is the effective elastic modulus, E is the initial elastic modulus, N is the number of microcracks in per unit volume, and a is the average radius of activated microcracks. For the three-dimensional problem, A � 16(1 − μ 2 ) (10 − 3μ)/(45(2 − μ)), where μ is the Poisson ratio and m � 3. For the two-dimensional problem, A � π and m � 2. e average radius of microcracks in rock mass can be calculated by Grady's expression [30]: where K IC is the fracture toughness, ρ r is the rock mass density, c p is the P-wave velocity, and _ ε max is the maximum of the volumetric strain rate. e abovementioned effective elastic modulus is the homogenization modulus value based on the statistical law. It can be approximately applied to the change laws of the average modulus in an annulus formed by the microelements with the distance of r from the borehole center; namely, E(r) � f (N(r)).
e number of microcracks caused by blasting is directly proportional to the peak value of explosion stress wave [31]; namely, where k d is the proportional coefficient under the action of explosion stress wave, p 0d is the peak value of pressure on the borehole wall under the action of explosion stress wave, r b is the borehole radius, r is the scaled distance, r is the distance from the borehole center, and α is the attenuation index. For the intact rock mass, α where μ d is the dynamic Poisson's ratio, while for the fractured rock mass, the value may be greater than 2.0. us, the damage factor under the action of explosion stress wave can be obtained from where for the three-dimensional problem, A d � 16(1 − μ 2 d ) (10 − 3μ d )/(45(2 − μ d )) and for the two-dimensional problem, A d � π.
Based on the thick-walled cylinder theory, the stress change law under the action of quasistatic gas is where p 0j is the pressure on the borehole wall under the action of quasistatic gas. It is assumed that the average radius of microcracks is invariable. In addition, Poisson's ratio and the number of activated microcracks under explosion stress wave are different from those under quasistatic gas, so A and k also change. us, the damage factor under the action of quasistatic gas can be obtained from where k j is the proportional coefficient under the action of quasistatic gas. For the three-dimensional problem, A j � 16(1 − μ 2 j )(10 − 3μ j )/(45(2 − μ j )), where μ j is the static Poisson's ratio. For the two-dimensional problem, Blasting damage is the result of the combined action of explosion stress wave and quasistatic gas, but their contributions to damage are different at different stages. In fact, the two actions intertwine and couple with each other. erefore, the damage factor under the combined action of the two can be expressed by [27] In conclusion, with the consideration of the combined action of explosion stress wave and quasistatic gas, the damage factor around a borehole is where the proportional coefficients, k d and k j , can be determined by experiments or analogy methods and the pressure on the borehole wall under decoupled charge, p 0d and p 0j , can be determined according to the expansion laws of detonation gas and detonation pressure. In particular, for the radial decoupled and axial continuous charge structure, the detonation pressure takes an average; namely, p m � 1/8ρ e D 2 m , where ρ e is the explosive density and D m is the detonation velocity. Rock mass can be dealt with as a plane strain problem. For the radial decoupled and axial interval charge structure, the detonation pressure is p m � 1/4ρ e D 2 m . Rock mass near the charge segment is constrained sideward by rock mass near the air segment and when simplified as a plane stress problem, there are some discrepancies with the reality, so it is more practical to consider it as a three-dimensional problem. Pressure on the borehole wall under explosion stress wave usually adopts the quasistatic pressure to multiply by the pressure enhancement coefficient which is determined by the radial decoupling coefficient values. When the radial decoupling coefficient is less than the critical value, the blasting effect belongs to the strong impact and the bigger pressure enhancement coefficient is chosen. Conversely, it belongs to the weak impact and the smaller pressure enhancement coefficient is chosen. In general, the pressure enhancement coefficients can be taken from 2 to 11. rough analysis, it is found that the broad trend of distribution curves under the radial decoupled charge is basically close, but the change rate of distribution curves and the maximum of damage factors are different so as to cause different damage ranges. A comparison of the distribution laws of damage factors under three decoupling coefficients is made. e material parameters are given in the second row of Table 1, and blasting parameters are given in the second row of Table 2. e curves of radial distribution laws are shown in Figure 2.
Damage factors under decoupled charge blasting decrease in a power function and finally tend to zero with the increase of distance from the center of the borehole, as shown in Figure 2.
is is mainly because the action of explosion stress wave is prominent at the initial stage and the contribution to the blasting damage is greater. With the decrease of explosion stress wave, the action of quasistatic gas is increasingly prominent. e theoretical formula gives full consideration to the combined action of explosion stress wave and quasistatic gas, which can characterize the radial distribution laws of damage factors and reflect the action process of rock mass blasting. e attenuation rate of damage factors decreases gradually when K ≥ K 0 and firstly increases and then decreases when K < K 0 , as shown in Figure 2. e main reason is that damage caused by explosion stress wave is smaller and the conditions offered for the action of quasistatic gas are limited when K ≥ K 0 , while damage caused by explosion stress wave is larger and the conditions offered for the action of quasistatic gas are superior when K < K 0 .
In Figure 2, when K ≥ K 0 , the explosion impact is greatly weakened and the serious damage zone is not generally produced. A serious damage zone whose largest boundary is located in the first turning point of the distribution curve is generated when K < K 0 . With other things being equal, the damage range decreases gradually with the increase of radial decoupling coefficients. When 0.19 is taken as a damage threshold value, the damage range of the decoupling coefficient of 1.31 is about 3.57 times that of the decoupling coefficient of 1.68. e difference of the damage range is related with explosive properties, charge structures and rock mass properties, and so on.
From the perspective of the size of the damage range and whether the serious damage zone is produced or not, the appropriate radial decoupling coefficient which is greater than the critical value is chosen in the field charge structure by explosive parameters and rock mass properties to protect the integrity and stability of surrounding rock.

Model Experiment.
Two cylinder models numbered A and B, respectively, are poured by using fine stone concrete. e diameter is 100 cm and the height is 50 cm. e curing period under the natural condition is more than 28 days. Meanwhile, the cubic standard samples of 15 cm are made. After curing, the physical and mechanical parameters are tested, which are listed in the third row of Table 1. It is observed that this concrete is equivalent to rock mass with low strength. A borehole in the center of the cylinder is drilled, as shown in Figure 3. e blasting parameters are listed in the third row of Table 2.
In the process of the experiments, the models are erected to test P-wave velocity of concrete by applying the opposite test method so as to avoid the second disturbance and the influence of water injection during the process of the coring test. Besides, damage tested by the opposite test method contains the comprehensive damage around the air and charge segment, which can test the maximum degree and range of damage. It is shown from the distribution characteristics of pressure in the borehole that damage around the charge segment is more serious than that around the air segment. erefore, damage around the charge segment is often a key concern and the opposite test results can reflect damage around the charge segment, which shows that the Shock and Vibration 5 test method is reasonable and feasible. Four mutually vertical P-wave test lines, named successively from I to IV in a clockwise direction, are laid out on the surface of the models. ere are eight test points in each line to verify each other, and then the average of P-wave velocity is taken to reduce the error caused by randomness of concrete properties. e area close to the borehole is densely arranged to reflect the trend of damage change obviously, while the area far away from the borehole is sparsely arranged. erefore, considering the radius of the concrete models and the radius of the test probes, the distance of the test points from the center of the borehole are 4 cm, 6 cm, 8 cm, 12 cm, 17 cm, 23 cm, 31 cm, and 40 cm in turn, as shown in Figure 3.
In the models numbered A and B, the radial decoupling coefficients of 1.25 and 2.0 are, respectively, used to carry out the blasting experiments. e explosive is No. 2 rock emulsion explosive with the charge amount of 15 g, which are detonated by an instantaneous electric detonator, as shown in Figures 4(a) and 4(b). e sandy clay is used to stem the borehole, as shown in Figure 4(c). ere are some axial air intervals near the charge segment. e condition after blasting is shown in Figure 4(d). And then the P-wave velocity of the models is tested by using the opposite test method and the RSM-SY5(T) sonic detector, as demonstrated in Figure 4(e). Before testing the P-wave velocity of the models, firstly, the area of the test points is polished to smoothness. Secondly, auxiliary circles with the same diameter as the probe are drawn with the test points as the center to ensure that the transmitting probe and the receiving probe are in alignment. irdly, the coupling agent is used to make the probe fit the concrete surface closely. e tested curves of P-wave are displayed in Figure 5.
e damage factors are calculated by bringing the tested results of P-wave velocity into (1). In addition, the coefficient A, the average radius of cracks, and the borehole pressure are calculated with the parameters in Tables 1 and 2. e proportional coefficient, k, is fitted by the blasting experiments. e theoretical curves are calculated by bringing the  [24] 2000 --3580 20.0 --  abovementioned parameters into (9). e tested values and the calculated curves of damage factors under the two radial decoupling coefficients are shown in Figure 6. In Figure 6, the test values of damage factors under the two radial decoupling coefficients decrease nonlinearly with the increase of the distance from the center of the borehole. When the radial decoupling coefficient is larger, the attenuation rate is faster and the maximum damage factor and the damage range around the borehole are smaller. Conversely, due to the appearance of the serious damage zone under the smaller radial decoupling coefficients, the attenuation rate is slower and the maximum damage factor and the damage range around the borehole are larger. e radial decoupling coefficient has a great influence on the damage range. When 0.19 is taken as a damage threshold value, the damage range of the decoupling coefficient of 1.25 is about 5 times that of the decoupling coefficient of 2.0, as shown in Figure 6. us, it is necessary to choose an appropriate radial decoupling coefficient for controlling the damage range from the perspective of the charge structure. e theoretical formula reflecting the combined action of explosion stress wave and quasistatic gas is adopted to calculate damage factors based on the model parameters. It is found that the correlation is high and the theoretical formula can describe the radial distribution laws of damage factors around the borehole so as to estimate the damage range, as displayed in Figure 6. Besides, when the radial decoupling coefficients change, it can also reflect the influence of the load characteristics in the borehole on damage; namely, when the radial decoupling coefficient is larger than the critical value, the impact effect of the blasting load on the borehole wall is greatly weakened.

Cited Model Experiment.
In the experiment [24], cement mortar models, whose cement-sand ratio was designed as 1 : 2.5 and water-cement ratio was 0.4, were made. e physical and mechanical parameters are presented in the fourth row of Shock and Vibration Table 1. e single-hole blasting is carried out in the models, and the blasting parameters adopted are listed in the fourth row of Table 2. e fitting relation between damage factors and the scaled distance is shown by the dashed line in Figure 7. A curve is obtained from the theoretical formula derived in this paper, as shown by the solid line in Figure 7.
On the basis that stress and strain around the borehole decrease in a power function with the increase of distance from the blasting source, Dai et al. [24] used the power function, D(r) � A(r/r b ) − B , where A and B reflect coefficients, such as rock mass properties, charge structure, detonation conditions, and so on, to fit out a quantitative relation which has a good consistency with the change laws of stress and strain in (4).
Obtained from the derived theoretical formula in this paper which can reflect the combined action of explosion stress wave and quasistatic gas, the radial distribution laws of damage factors have a good correlation with the tested values, as shown in Figure 7. Meanwhile, the correlation coefficient is higher than that fitted by a power function because the attenuation index under the action of explosion stress wave is different from that under the action of quasistatic gas. e same attenuation index was used in [24], while the two different attenuation indexes are considered in the theoretical formula derived in this paper.
It is thought that the attenuation laws of damage factors have a good consistency with the change laws of stress and strain [24]. Damage of rock mass is the result of the actions of explosion stress and strain from the aspect of mechanical mechanism. Based on the effective elastic modulus and the change law of stress, this paper derives a theoretical formula under the combined action of explosion stress wave and quasistatic gas. e formula can describe the distribution laws of damage factors and estimate the damage range, which can not only provide an important reference for studying the damage characteristics, but also reveal the blasting damage mechanism. erefore, the formula has a great theoretical significance.

Conclusion
e formula of the radial distribution laws of damage factors caused by decoupled charge blasting under the combined action of explosion stress wave and quasistatic gas is derived and verified. e following conclusions are drawn: (1) e derived formula can describe the radial distribution laws of damage factors, embody the influence of the radial decoupling coefficients on damage factors, and communicate the interrelation between damage factors and blasting parameters. Meanwhile, the formula is of high accuracy and wide practicability. It can be conveniently applied to engineering to guide the optimization design of blasting parameters (2) e attenuation rate of damage factors decreases gradually when K ≥ K 0 and first increases and then decreases when K < K 0 . In addition, the attenuation rate is in positive correlation with the radial decoupling coefficients at the initial stage. But with other things being equal, the damage range is in negative correlation with the radial decoupling coefficients  [24], and empirical formula curve [24] of damage factors under decoupled charge blasting.

Shock and Vibration
(3) e formula of the radial distribution laws of damage factors is derived from the perspective of mechanics. Parameters in the formula have definite physical meanings. And the formula has a relatively simple form and a strict theoretical basis. Besides, the mechanical mechanism of blasting damage can be revealed. erefore, it can enrich the damage theory of rock mass blasting (4) In engineering, an appropriate radial decoupling coefficient which is larger than the critical value should be adopted to the charge structure of smooth or presplit blasting on the premise of stable detonation so as to reduce the disturbance of the blasting action to surrounding rock and protect the integrity and stability of rock mass.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.