Experimental and Numerical Studies on the Mechanical Properties of Columnar Jointed Basalt

Abstract

Columnar jointed basalts (CJB) are featured by the joint network of discontinuities, which is hard for geotechnical design and construction. Therefore, investigating the mechanical behaviors of CJB is significant for the long-term use of the engineering structures. In this paper, field research studies, such as the rigid bearing plate tests and acoustic tests, were employed, to obtain the deformation modulus and the acoustic velocity of columnar jointed basalts, and the formula expressed by deformation modulus and acoustic velocity was then established, which would be a useful guidance for the site operation. Based on the monitoring and testing data of CJB, several numerical simulation models with different joint angles and weak thickness were built to further discuss the mechanical behaviors of CJB. The numerical simulation results show that the joint angle of 30° ∼ 60° presents weak antipressure abilities, for the shearing slipping force and interformational sliding are remarkable in them. Besides, the small material property difference is helpful for the compressive strength of rock masses. After adopting the most unfavorable joint angle and width to model the tunneling process in CJB, it suggests that it is easily damaged along the directions of rock joints in CJB when constructing in tunnels.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

The columnar jointed basalts (CJBs) are particular geological structures formed during the cooling of basalt, which usually cut basalt into regular or irregular prisms. Due to the internal implicit joints and unloading cracks in CJB, the mechanical properties of rock masses are complicated. It was discovered that during the excavation of CJB, the special failure patterns and mechanical responses normally occurred [1–4]. In China, the columnar jointed basalts are widely distributed at several regions, which create great challenges for the design and construction of engineering projects, such as the tunnels, mines, and dams [5], for their special jointing patterns [6].

Therefore, investigating the mechanical behaviors of CJB has become urgent matters in recent years due to their particular damage modes and mechanical behaviors during the CJB tunneling and excavation, which poses adverse effects on project stability. For example, in tunnel projects, during the excavation process of underground tunnels in CJB, collapses may easily occur due to their unfavorable mechanical properties [7–10]; in dam projects, the dam foundation bears substantial hydraulic pressure during the operation period, but rock mass materials of the dam foundation are normally prone to failure due to excavation and mechanical loading, especially in CJB geological conditions [11, 12]. Hence, evaluating the mechanical properties and construction stability of CJB is necessary and essential [13]. After a review of related studies, it shows that there are several approaches to implement related work, such as field surveys [14], experimental tests, and numerical simulations [15–17]. Most of the scholars combined both methods to carry out the scientific research studies.

Several field investigations were employed in the areas where columnar jointed basalts occur [18]. The affecting factors of mechanical properties of columnar jointed basalts primarily contain the distribution of internal joints, geological environment [19–22], pore water pressure [23], and climatic conditions. A uniform distribution survey line method was put forward to reveal the structure, deformation, and failure characteristics from formation to their exposed state of CJB [24].

Experimental tests contain two types which are in situ tests and laboratory tests [25]. In situ tests are the most normal and conventional way in determining the parameters of rock masses. In situ testing methods, such as field triaxial tests, plate bearing tests, microseismic monitoring technologies, and acoustic tests [26], have been widely employed, and a number of research achievements are reported. Carrying out field tests shows that rock mass time-dependent relaxation, cracking, and collapse are dramatic features in CJB [27]. Laboratory tests on rock mass and physical modeling tests have been conducted as well [28–30], to evaluate mechanical behavior and failure mechanisms [31, 32]. Limited by the scale of laboratory tests, field tests become an effective way to confirm mechanical properties of rock masses. Observation window methods, 3D scanning methods, and P-wave tests were employed to study the features of internal implicit joints of CJB at drainage tunnels, and it suggested that the CJB is heterogeneous, showing the characteristics of mosaic structure and mainly quadrilateral in shape, and the joints are distributed with distinct tendency, which could influence the P-wave anisotropy of CJB. Besides, the artificial CJB specimens were made to conduct the uniaxial compression tests, and the anisotropic deformation strength was discussed. Scanning electron microscope (SEM) tests were also adopted to analyze the mechanical behaviors and microscopic failure of rock masses, which indicated that the confining pressure is essential for the mechanical behaviors of the rock masses [33].

Numerical simulation methods are quick and effective approaches to investigate the mechanical behaviors of CJB. A novel geometric modeling method based on the topological properties of field columnar basalt was employed, which provided a method to estimate rock mass modulus using numerical simulation. 3DEC was used to simulate the uniaxial compression tests, which indicated that the elastic-plastic and elastic brittle behaviors are different with the presence of confining pressure [15, 34, 35]. Chen et al. [36] adopted PFC^3D to investigate the effect of the joint inclination on the mechanical behavior and associated cracking process and displacement. A modified constrained centroid Voronoi tessellation (CCVT) algorithm was put forward to analyze the influence of columnar jointed structure on the mechanical properties of CJB. Zhou et al. [37] addressed a statistical damage model based on the Weibull distribution considering the joint orientation associated with Jaeger’s and modified Hoek–Brown failure criteria for CJB, and the discrete element method of PFC was employed to verify the model proposed. A computational interaction algorithm with high numerical accuracy and efficiency under the loading conditions of plasticity, damage, and coupled elastoplastic damage was developed.

As a consequence, the objective of this paper is to investigate the mechanical behaviors of CJB. Firstly, field research studies were conducted in the engineering project, where the CJB was widely distributed. The in situ research studies mainly contained the rigid bearing plate tests and acoustic tests on columnar jointed basalts, to study the particular mechanical properties of CJB. Secondly, based on the field testing data, the formula between deformation modulus and acoustic velocity was established. For the further discussion of the mechanical properties of CJB, several numerical simulation models with different joint angles and widths were built. After confirming the most unfavorable joint angle and width of CJB, an engineering example about tunneling in CJB with the most adverse geological conditions was modeled.

2. Engineering Background

The Baihetan hydropower station contains a super-high arch dam, whose height is 289 m, and it is built in the southwest of China. Furthermore, it is also the second largest hydropower station around the world, whose installed capacity is 16000 MW. Figure 1(a) depicts the topographical range of Baihetan arch dam, and it can be seen that Baihetan arch dam is constructed in an asymmetrical V-shaped valley, whose left slope is gentle and the right slope is steep. The width of the river valley is from 449 m to 534 m at the normal pool level. Figure 1(b) illustrates the geological condition along the dam axis, which shows that the rock masses at the dam site mainly contain Permian Emeishan basalt and breccia lava with 11 layers of basalt, i.e., P₂β₁-P₂β₁₁. Besides, the faults and shear belts are in wide distribution, F₁₇ is the largest fault region in the left slope, the F₁₄ and F₁₆ extend across the riverbed close to the right slope. The interlayer shear belts are widely developed on the top of each basalt layer, for example, C₃ and C_3-1 are located at the left slope from the elevation of 720 m to 760 m, C₃, C_3-1, C₄, and C₅ in the right side of the dam are distributed from the elevation of 640m to 780m. Some intrastratal shear belts, such as Ls₃₃₁, Ls₃₃₁₈, Ls₃₃₁₉, Rs₃₃₁, and Rs₃₃₆, mainly reveal in the third basalt stratum, which are ${\text{P}}_2{\beta }_3^{2-2}$ and ${\text{P}}_2{\beta }_3^3$. The range of dam foundation on the left bank of the Baihetan hydropower station refers to the elevation of the left bank is 600∼834 m, the range of dam foundation on the right bank is 600∼834 m, the range of dam foundation on the riverbed refers to the elevation of both river banks is below 600, and the lowest elevation of the dam foundation is 545 m.

[figures omitted; refer to PDF]

Figure 2 shows the field images of different types of rock masses in Baihetan arch dam, there are three kinds of rock masses, i.e., columnar jointed basalt (CJB), breccia lava, and cryptic basalt. The columnar jointed basalt (CJB) is looser than other two rock masses; thus, the anisotropy is more significant. However, the columnar jointed basalt (CJB) is a common primary structure in basalt, and it extensively spreads in the Baihetan hydropower station. The CJB can be divided into three types, which mainly conclude Category-I CJB, Category-II CJB, and Category-III CJB, and they are depicted in Figure 1(b) [14]. The columnar jointed basalts are primarily located in ${\text{P}}_2{\beta }_2^2$, ${\text{P}}_2{\beta }_2^3$, ${\text{P}}_2{\beta }_3^2$, ${\text{P}}_2{\beta }_3^3$, ${\text{P}}_2{\beta }_4^1$, ${\text{P}}_2{\beta }_6^1$, ${\text{P}}_2{\beta }_7^1$, and ${\text{P}}_2{\beta }_8^2$, and they are nonuniform with different columnar size and length. Meanwhile, in the riverbed and both slopes, the CJB is widely formed in the rock strata ${\text{P}}_2{\beta }_3^2$ and ${\text{P}}_2{\beta }_3^3$, and it is vital to investigate the rock masses at the foundation surface. In view of the adverse geological conditions in the Baihetan dam foundation, during the construction, the excavation stability must be taken into consideration. Figure 3 shows the layout of the field testing positions at the dam foundation, it illustrates that cryptic basalt, breccia lava, and CJB are distributed in different elevations, and several testing points are arranged at suitable positions. In practical engineering, at least one group and three points should be arranged in each group for different lithology, different strata, and different parts. Generally speaking, no less than five groups (a total of 15 points) should be appropriate, to fully reflect the situation of different lithologies with different integrity degrees. In this paper, the test of 30 group CJB rigid pressure plates on the left bank foundation of the Baihetan hydropower station is taken as an example.

[figures omitted; refer to PDF]

[figure omitted; refer to PDF]

As for the excavation process, the excavation of the slope at the abutment of the left bank dam foundation began in September 2013 at the elevation of 834.0 m, and the anchor cable construction of the foundation surface above the elevation of the dam foundation was completed in November 2016. The excavation of the slope at the abutment area of the right bank foundation began in October 2014 at the elevation of 834.0 m, and the excavation of the right bank foundation reached the elevation of 600m in September 2016. Figure 4 shows the sketch of dam foundation excavation, and the dynamic loading is applied upon the protecting zone. According to the existing research achievements, the thickness of the protecting zone is confirmed associated with the geometry of blasting-induced damage zone, and the damage zone is defined as the rock region where has been influenced by the blast loading excavation [38]. The thickness of dam foundation protecting zone of some constructed large-scale projects is listed in Figure 4 [39–42]. Therefore, based on the comprehensive comparison, the dam foundation of the riverbed began to be excavated by reserving a 5 m protective layer in October 2016. The protective layer of the riverbed was fully excavated in January 2017, and the excavation was completed by March 10. In the excavation process, the noncolumnar jointed rock mass is excavated according to a ladder section per 10 m. In the excavation of dam foundation, the phenomenon of unloading relaxation and fracture of column-jointed basalt was relatively obvious, so the excavation of column-jointed basalt rock mass was carried out according to a ladder section every 5 m.

[figure omitted; refer to PDF]

When acoustic waves from different measuring points were monitored from the original acoustic record and travel to the two receiving transducers, the distance between the two receiving transducers was divided by the time difference, and then the acoustic velocity of the pore wall rock mass in this section can be obtained. The formulas are as follows:$$\begin{array}{}\text{(2)}& \Delta t=t_1-t_2,\text{(3)}& v_p=\frac{L}{\Delta t}\times {10}^6,\end{array}$$where t₁ and t₂ are the arriving time to two receiving transducers of the same excitation signal, respectively; Δt is the time difference; L is the distance between two receiving transducers, and L = 0.2 m, which is a fixed value; $v_p$ is the speed of the acoustic wave.

3.3. Experimental Results of Rigid Bearing Plate and Acoustic Tests

According to the results of rigid bearing plate and acoustic tests, the deformation modulus, elastic modulus, and the acoustic velocity were monitored, as listed in Table 1. Besides, the acoustic waves propagate in different types of rock masses, the propagation velocity could be reduced to a vary degree, thus defining the attenuation rate η, to estimate the quality of rock masses, and it could be expressed by the following formula:$$\begin{array}{}\text{(4)}& \eta =1-\frac{v_w}{v_0},\end{array}$$where $v_0$ is the incident velocity of acoustic wave and $v_w$ is the transit velocity of acoustic wave in the rock masses. Normally, when η is less than 10%, it is considered that the excavation has little or no influence on the rock masses; when the value of η is between 10% and 15%, the rock masses are considered slightly affected by the excavation; when η is larger than 15%, it is believed that the excavation has influence on the rock mass quality.

Table 1

Results of field tests at the left bank of the Baihetan hydropower station [43].

Monitoring position	Elevation (m)	Elastic modulus Ee (GPa)	Deformation modulus E0 (GPa)	$v_0$ (m/s)	$v_0$ (m/s)	η (%)
EA2-1	653.0	15.2	9.8	4900	4350	11.2
EA2-2	651.5	15.1	9.1	4900	4541	7.3
EA2-3	653.0	7.1	4.3	4350	3685	15.3
EC1-1	651.2	5.6	3.5	4900	4147	15.4
EC1-2	651.4	5.0	2.9	4350	3742	14.0
EC2-1	651.3	6.2	3.5	4900	3781	22.8
EB1-1	656.2	2.4	1.4	4350	3492	19.7
EB1-2	656.2	1.8	1.0	4350	3278	24.6
EB2-1	656.5	6.9	3.2	4900	3620	26.1
EL.632-1	632.0	17.7	13.2	4900	4739	3.3
EL.632-2	632.0	13.2	8.7	4900	4659	4.9
EL.632-3	632.0	19.3	12.2	4900	4577	6.6
EL.637-1	637.0	6.0	4.4	4900	3580	26.9
EL.637-2	638.0	6.8	3.4	4900	4413	9.9
EL.637-3	636.0	13.1	8.2	4900	4332	11.6
EL.647-1	647.0	12.3	7.4	4900	4757	2.9
EL.647-2	646.0	17.0	9.8	4900	4440	9.4
EL.647-3	648.0	9.5	5.4	4900	3937	19.7
E11-1	592.0	11.4	8.5	4900	4833	1.4
E11-2	592.0	6.4	4.4	4900	4464	8.9
E11-3	592.0	7.1	4.5	4900	4375	10.7
E11-4	587.0	29.7	23.7	4900	5125	0
E11-5	587.0	12.6	9.8	4900	4802	2.0
E11-6	587.0	6.8	5.1	4900	4569	6.8

In Table 1, there is a great dispersion of the data obtained from field experiments, which indicates that the mechanical behaviors of CJB are significantly different. According to the data in Table 1, the diagram of elastic modulus and deformation modulus varying with elevation is depicted and shown in Figure 7. Figure 7(a) shows the data points’ distribution of elastic modulus and deformation modulus in different elevation. It illustrates that the mechanical behaviors of CJB are different even in the same elevation. Figure 7(b) shows the bar graph, which represents the mean value and data standard deviation of elastic modulus and deformation modulus of CJB in the same elevation. From the figure, it indicates that the mechanical property differences in low elevation are larger than those in high elevation.

[figures omitted; refer to PDF]

The elevation from 650 m to 660 m is the region of columnar jointed basalts, thus conducting further analysis in this area. From Figure 3, it can be seen that there are nine testing points in the area of A2, B1, B2, C1, and C2: three points are in A2 area, two points are in area B1, one point is in area B2, two points are in area C1, and one point is in area C2. In the upstream of the test area, the intraformational faulted belt LS₃₃₁₉ was observed in the area of A1 and A2, the occurrence is N10°E, SE∠25°, and the bandwidth is 5∼15 cm, as shown in Figure 1. In addition to EA2-1 and EA2-2 points located in the footwall of LS₃₃₁₉, the other seven points were located in the upper wall of LS₃₃₁₉. Under its influence, the deformation modulus values of EA2-1 and EA2-2 points are 9.8 GPa and 9.1 GPa, respectively, with an average of 9.5 GPa. The deformation modulus values of the remaining seven points ranged from 1.0 GPa to 4.3 GPa, with an average of 2.8 GPa. During the process of dam foundation excavation, if dislocation interfaces are distributed in the protecting zone, the rock masses could be easily affected. The rock mass relaxation phenomenon in the upper wall of the dislocation interfaces is generally more serious than those in the footwall. Therefore, in this area, the intraformational faulted belt LS₃₃₁₉ indeed has great effects on the rock mass deformation modulus.

The deformation modulus values of the three points in the area A are the highest, which are significantly different from each other due to the influence of the rock structure. The deformation modulus values in B area are minimum, whose values are from 1.7 GPa to 3.17 GPa. The deformation modulus value in C area is from 2.9 GPa to 4.0 GPa, which corresponds to the real situation. Meanwhile, combined with the data in Table 1, it can be seen that the acoustic velocity in the area of CJB ranges from 3278 m/s to 5125 m/s, with an average of 4207.2 m/s. The acoustic velocity attenuation rate of η in CJB is within 0 and 26.9%, and the proportion of η less than 10% is 43.3%; between 10% and 15% is 13.3%; and larger than 15% is 43.3%. Therefore, the CJB is sensitive to the elevation changing.

Figure 8(a) shows the results of in situ tests of CJB at different stages. In 2015, nine groups of tests were conducted on CJB at the elevation of 650∼660 m. With the downward excavation of the dam foundation, nine groups of tests were added to CJB at the elevation of 630∼650 m in 2016. Finally, in 2017, six groups of tests were added to CJB at the elevation of 585∼595 m and 563∼570 m, respectively, so that the number of left bank foundation testing points was increased to thirty. Compared with the fitting results using all testing data, the fitting curve in 2015 (N = 9) is located above it, and under the same acoustic velocity, the deformation modulus is larger. The curve in 2016 takes the acoustic velocity of 4.5 km/s as the cutoff point, and when $v_p$ is less than 4.5 km/s, the deformation modulus is larger, while it is smaller when $v_p$ is greater than 4.5 km/s. In 2017, the curve takes the P-wave velocity of 4.8 km/s as the cutoff point, and when $v_p$ is less than 4.8 km/s, the deformation modulus is smaller, while it is larger when it is larger than 4.8 km/s. Considering that a small number of testing points will produce great uncertainty to the fitting result, the parameter selection of the foundation surface in the later stage should be determined according to the comparison relation with more number of testing points.

[figures omitted; refer to PDF]

There are generally two forms of E₀-$v_p$ function, i.e., power function and exponential function, as shown in Figure 8(b). The power function can be expressed by $E_0=0.01398e^{1.40806v_p}$, and the exponential function can be fitted by $E_0={{\displaystyle 0.000674v_p}}^{6.244}$. When $v_p$ is lower than 4.0 km/s or higher than 4.75 km/s, the exponential function is higher than that of the power function. When $v_p$ is between 4.0 km/s and 4.75 km/s, the power function is slightly higher than that of the exponential function, and there is almost no difference between them. In addition, when the same data were used to establish the power function and exponential function, the difference was within 3%.

The different mechanical properties at different elevation of CJB have been discussed above, and the fitting functions of deformation modulus and acoustic velocity are proposed, which are valuable for the determining parameters of CJB by rigid bearing plate tests or acoustic tests. Therefore, in the next section, different types of CJB with different joint angles and thickness of weak intercalations are further discussed by numerical simulation methods.

4. Numerical Simulation

4.1. Model Establishing and Property Assignment

In this section, the numerical simulation software is employed to deal with the problems of columnar jointed basalt, and several numerical models were established with different angles and widths of joints in rock mass, to analyze the strength and failure modes of irregular CJB. As Figure 9 shows, the cubical CJB samples with 50 mm width and 50 mm height were modeled, and there are seven angles of the joint, i.e., β = 0°, 15°, 30°, 45°, 60°, 75°, and 90°. Besides, different thickness of weak plane in CJB, such as 1.25 mm, 2.5 mm, 3.75 mm, 5 mm, 6.25 mm, and 10 mm, was also built under the most unfavorable angle of joints in CJB, as shown in Figure 9.

[figure omitted; refer to PDF]

The constitutive model named ubiquitous-joint model can be adopted to simulate the presence of an orientation of the weak plane in Mohr–Coulomb material where shear failure can be simulated. The compressive strength of the rock mass is a function of the material and joint properties, formed by the direction of the compressive stress and its projection onto the plane of weakness. Therefore, the rock mass of CJB is assigned by the Mohr–Coulomb, and the joint in CJB is assigned ubiquitous-joint properties. The properties adopted in this paper for rock mass and joints are listed in Table 2.

Table 2

Mohr–Coulomb and ubiquitous-joint properties used in numerical simulation.

Material	Bulk modulus (GPa)	Cohesion (kPa)	Dilation (°)	Friction angle (°)	Tension limit (kPa)
CJB	10	7000	0	40	1000
Joint	—	300	0	29	50

Based on the analysis results of numerical simulation of uniaxial compression, the adverse joint angle and thickness are confirmed; thus, the excavation stability of the drainage tunnel can be further discussed. The materials properties are assigned as the parameters listed in Table 2. As shown in Figure 10, the dimension of the model is 50mm in length (L), 50 in width (W), and 25 m in depth (h). According to the principle of tunnel symmetry, the tunnel is set as a quarter section drainage tunnel, and the diameter of the tunnel is 4 m (D). The width of joints in CJB is 10 cm, and the excavation length of each step is 5 m.

[figure omitted; refer to PDF]

Figure 12 shows the cutting plane contour of the plastically deforming area in CJB with different joint angles calculated by FLAC^3D, it indicates that the joint angle from 15° to 75° shows a complicated stress state, which tends to be easily damaged. When the joint angle is 0° or 90°, the core of rock sample is mostly bearing compression and lack of shear slipping, and they are bearing stress states in the compression process. However, with the increase of joint angle, the element in CJB bears different types of stress, especially in the joint angle of 45° to 75°.

[figure omitted; refer to PDF]

Figure 13 shows the stress-strain curves of CJB with different joint thickness, and the values of joint thickness are 0 mm, 1.25 mm, 2.5 mm, 3.75 mm, 5 mm, 6.25 mm, and 10 mm, respectively. With the increase of thickness of weak intercalation, the stress strength gradually decreases and then there was a tiny increase, which indicates that the small material property difference is helpful for the compressive strength of rock masses. When the thickness of weak intercalation increases from 0 mm, the stress strength decreases sharply, because the mechanical properties of weak intercalation are in poor quality, and it has a great influence on the overall strength of rock mass. When the thickness of weak intercalation increases to a certain value (3.75 mm), the stress strength of CJB appears a tiny increase, because the rock specimen gradually becomes intact, and it is advantageous for the stress strength of CJB.

[figure omitted; refer to PDF]

Figure 14 shows the cutting plane contour of the plastically deforming area in CJB with different joint thickness calculated by FLAC^3D, and it suggests that the stress state in each thickness of joint in CJB is variety during the compression process. By comparison, the thickness of 6.25 mm and 10 mm is in few stress state, which indicates that the stress strength in these cases is higher than others, and this result is corresponding to the stress-strain curve depicted before.

[figure omitted; refer to PDF]

Figure 15 shows the curves of maximum displacement in each excavation step in CJB and intact rock masses. It illustrates that the CJB is adverse geological conditions for tunneling, as it is easily damaged along the directions of rock joints in them. In Figure 15(a), to compare the different mechanical properties of CJB, tunneling in intact rock masses is modeled as well, and the maximum displacement of at the top of tunnel is monitored. The figure suggests that the maximum displacement in CJB is larger than in intact rock masses. Following the excavation step, the vault crown settlement gradually increases and then tends to be stable. Within the excavation length of 20 m, the displacement variation value is in high level. The reasons for this are mainly the initial moment of tunneling is unstable, ground disturbance is remarkable, besides, the number of joints in CJB is larger than the later excavation zone, as shown in Figure 15(b), which easily causes rock slipping and instable failure.

[figures omitted; refer to PDF]

5. Conclusion

Detailed in situ tests of columnar jointed basalts (CJBs) in the Baihetan hydropower station dam foundation have been conducted, and the field tests mainly include rigid bearing plate tests and acoustic tests. Based on the monitoring data from the tests, the relationship between deformation modulus and acoustic velocity was expressed by function expressions. To further investigate the mechanical behaviors of CJB with different joint angles and widths, the numerical simulation models were then established, and the material parameters of rock masses were assigned related to the data from field tests. According to the results of numerical simulation, the most adverse joint angle and width were proposed. Finally, adopting the most unfavorable joint angle and width to model the tunneling process in CJB, it suggests that it is easily damaged along the directions of rock joints in CJB when constructing in tunnels. The conclusions are as follows:

(1) The discreteness of data monitored from the field tests is large, and the mechanical behaviors of CJB are obviously different from those in the same elevation, the mechanical property differences in low elevation are larger than those in high elevation of Baihetan hydropower station, which means that the joints state and distribution patterns in CJB have great influence on the mechanical behaviors of rock masses.

(2) The power function $E_0=0.01398e^{1.40806v_p}$ and exponential function $E_0={{\displaystyle 0.000674v_p}}^{6.244}$ can be well employed to express the function relationship between deformation modulus and acoustic velocity, which are valuable for the determining parameters of CJB by rigid bearing plate tests or acoustic tests.

(3) The CJB with different joint angles performs distinct mechanical properties, and the angles of 30°∼60° present weak antipressure abilities, because the shear slipping force in CJB with 0° joint is the smallest, while is more remarkable in CJB with the angle of 30°∼60°, which may occur instable failure.

(4) When the thickness of weak intercalation increases from 0 mm, the stress strength decreases sharply, because the mechanical properties of weak intercalation are in poor quality, and it has a great influence on the overall strength of rock masses. When the thickness of weak intercalation increases to a certain value, the stress strength of CJB appears a tiny increase, because the rock specimen gradually becomes intact, and it is advantageous for the stress strength of CJB.

(5) When tunneling in columnar jointed basalts, it is easily damaged along the directions of rock joints; thus, in the related constructions, it is significantly important to take some measures to prevent the damage of the joints.

Acknowledgments

The authors are grateful for the financial supports from the Department of Communications of Zhejiang Province Science Foundation of China (grant no. 2020046) and the State Key Program of the National Natural Science of China (grant no. 52038008).

References

[1] Y. Xia, C. Zhang, H. Zhou, "Structural characteristics of columnar jointed basalt in drainage tunnel of Baihetan hydropower station and its influence on the behavior of P-wave anisotropy," Engineering Geology, vol. 264,DOI: 10.1016/j.enggeo.2019.105304, 2020.

[2] P. Lin, J. Shi, P. Wei, Q. Fan, Z. Wang, "Shallow unloading deformation analysis on Baihetan super-high arch dam foundation," Bulletin of Engineering Geology and the Environment, vol. 78 no. 8, pp. 5551-5568, DOI: 10.1007/s10064-019-01484-4, 2019.

[3] Q.-X. Meng, H.-L. Wang, W.-Y. Xu, Y.-L. Chen, "Numerical homogenization study on the effects of columnar jointed structure on the mechanical properties of rock mass," International Journal of Rock Mechanics and Mining Sciences, vol. 124,DOI: 10.1016/j.ijrmms.2019.104127, 2019.

[4] G. Feng, X. C. Wang, M. Wang, Y. Kang, "Experimental investigation of thermal cycling effect on fracture characteristics of granite in a geothermal-energy reservoir," Engineering Fracture Mechanics, vol. 235,DOI: 10.1016/j.engfracmech.2020.107180, 2020.

[5] X. Wang, C. Liu, S. Chen, L. Chen, K. Li, N. Liu, "Impact of coal sector’s de-capacity policy on coal price," Applied Energy, vol. 265,DOI: 10.1016/j.apenergy.2020.114802, 2020.

[6] Y. H. Hatzor, X.-T. Feng, S. Li, G. Yagoda-Biran, Q. Jiang, L. Hu, "Tunnel reinforcement in columnar jointed basalts: the role of rock mass anisotropy," Tunnelling and Underground Space Technology, vol. 46,DOI: 10.1016/j.tust.2014.10.008, 2015.

[7] H. Bao, K. Zhang, C. Yan, H. Lan, F. Wu, H. Zheng, "Excavation damaged zone division and time-dependency deformation prediction: a case study of excavated rock mass at Xiaowan hydropower station," Engineering Geology, vol. 272,DOI: 10.1016/j.enggeo.2020.105668, 2020.

[8] H. Ji, J. C. Zhang, W. Y. Xu, "Experimental investigation of the anisotropic mechanical properties of a columnar jointed rock mass: observations from laboratory-based physical modelling," Rock Mechanics and Rock Engineering, vol. 50 no. 7, pp. 1919-1931, DOI: 10.1007/s00603-017-1192-4, 2017.

[9] N. Xu, J. Wu, F. Dai, "Comprehensive evaluation of the stability of the leftbank slope at the Baihetan hydropower station in southwest China," Bulletin of Engineering Geology and the Environment, vol. no. 4, 2017.

[10] Y.-X. Xiao, X.-T. Feng, B.-R. Chen, G.-L. Feng, Z.-B. Yao, L.-X. Hu, "Excavation-induced microseismicity in the columnar jointed basalt of an underground hydropower station," International Journal of Rock Mechanics and Mining Sciences, vol. 97, pp. 99-109, DOI: 10.1016/j.ijrmms.2017.04.012, 2017.

[11] F. Wu, J. Liu, T. Liu, H. Zhuang, C. Yan, "A method for assessment of excavation damaged zone (EDZ) of a rock mass and its application to a dam foundation case," Engineering Geology, vol. 104 no. 3-4, pp. 254-262, DOI: 10.1016/j.enggeo.2008.11.005, 2009.

[12] Y. L. Bao, J. L. Wan, "The innovation on foundation excavation for Laxiwa double arch dam," Journal of Northwest Hydroelectric Power, vol. 22 no. 5, pp. 79-81, 2006. in Chinese

[13] X. Shen, X. Niu, W. Lu, "Rock mass utilization for the foundation surfaces of high arch dams in medium or high geo-stress regions: a review," Bulletin of Engineering Geology and the Environment, vol. 76 no. 2, pp. 795-813, DOI: 10.1007/s10064-016-0892-4, 2016.

[14] J. Cui, Q. Jiang, X. Feng, "Insights into statistical structural characteristics and deformation properties of columnar jointed basalts: field investigation in the Baihetan dam base, China," Bulletin of Engineering Geology and the Environment, vol. 77 no. 2, pp. 775-790, DOI: 10.1007/s10064-017-1029-0, 2017.

[15] M. Espada, J. Muralha, J. V. Lemos, "Safety analysis of the left bank excavation slopes of baihetan arch dam foundation using a discrete element model," Rock Mechanics and Rock Engineering, vol. 51 no. 8, pp. 2597-2615, DOI: 10.1007/s00603-018-1416-2, 2018.

[16] J. Wang, Y. Zhang, Z. Qin, S. Song, P. Lin, "Analysis method of water inrush for tunnels with damaged water-resisting rock mass based on finite element method-smooth particle hydrodynamics coupling," Computers and Geotechnics, vol. 126,DOI: 10.1016/j.compgeo.2020.103725, 2020.

[17] W. Han, G. Li, Z. Sun, H. Luan, C. Liu, X. Wu, "Numerical investigation of a foundation pit supported by a composite soil nailing structure," Symmetry, vol. 12 no. 2,DOI: 10.3390/sym12020252, 2020.

[18] Q. Fan, Z. Wang, J. Xu, M. Zhou, Q. Jiang, G. Li, "Study on deformation and control measures of columnar jointed basalt for baihetan super-high arch dam foundation," Rock Mechanics and Rock Engineering, vol. 51 no. 8, pp. 2569-2595, DOI: 10.1007/s00603-017-1378-9, 2017.

[19] B. Chen, S. Zhang, Y. Li, Z. Li, H. Zhou, "Physical simulation study of crack propagation and instability information discrimination of rock-like materials with faults," Arabian Journal of Geosciences, vol. 13 no. 18,DOI: 10.1007/s12517-020-05966-8, 2020.

[20] J. Chen, J. Zhao, S. Zhang, Y. Zhang, F. Yang, M. Li, "An experimental and analytical research on the evolution of mining cracks in deep floor rock mass," Pure and Applied Geophysics, vol. 177 no. 11,DOI: 10.1007/s00024-020-02550-9, 2020.

[21] J. Xu, A. Haque, W. Gong, "Experimental study on the bearing mechanisms of rock-socketed piles in soft rock based on micro X-ray CT analysis," Rock Mechanics and Rock Engineering, vol. 53 no. 8, pp. 3395-3416, DOI: 10.1007/s00603-020-02121-3, 2020.

[22] H. Y. Pan, D. W. Yin, N. Jiang, Z. G. Xia, "Crack initiation behaviors of granite specimens containing crossing-double-flaws with different lengths under uniaxial loading," Advances in Civil Engineering, vol. 2020,DOI: 10.1155/2020/8871335, 2020.

[23] D. Liu, Z. Gu, R. Liang, "Impacts of pore-throat system on fractal characterization of tight sandstones," Geofluids, vol. 2020,DOI: 10.1155/2020/4941501, 2020.

[24] Q. Jiang, X.-T. Feng, Y. H. Hatzor, X.-J. Hao, S.-J. Li, "Mechanical anisotropy of columnar jointed basalts: an example from the Baihetan hydropower station, China," Engineering Geology, vol. 175, pp. 35-45, DOI: 10.1016/j.enggeo.2014.03.019, 2014.

[25] C. Zhu, M. He, M. Karakus, X. Cui, Z. Tao, "Investigating toppling failure mechanism of anti-dip layered slope due to excavation by physical modelling," Rock Mechanics and Rock Engineering, vol. 53 no. 11,DOI: 10.1007/s00603-020-02207-y, 2020.

[26] J. Yang, J. Dai, C. Yao, S. Jiang, C. Zhou, Q. Jiang, "Estimation of rock mass properties in excavation damage zones of rock slopes based on the Hoek-Brown criterion and acoustic testing," International Journal of Rock Mechanics and Mining Sciences, vol. 126,DOI: 10.1016/j.ijrmms.2019.104192, 2020.

[27] Q. Fan, X. Feng, W. Weng, Y. Fan, Q. Jiang, "Unloading performances and stabilizing practices for columnar jointed basalt: a case study of Baihetan hydropower station," Journal of Rock Mechanics and Geotechnical Engineering, vol. 9 no. 6, pp. 1041-1053, DOI: 10.1016/j.jrmge.2017.07.003, 2017.

[28] D. Zhong, Y. Liu, L. Cheng, Q. Yang, Y. Chen, "Study of unloading relaxation for excavation based on unbalanced force and its application in baihetan arch dam," Rock Mechanics and Rock Engineering, vol. 52 no. 6,DOI: 10.1007/s00603-018-1653-4, 2018.

[29] D. Sow, C. Carvajal, P. Breul, "Modeling the spatial variability of the shear strength of discontinuities of rock masses: application to a dam rock mass," Engineering Geology, vol. 220, pp. 133-143, DOI: 10.1016/j.enggeo.2017.01.023, 2017.

[30] L. Yan, W. Xu, R. Wang, Q. Meng, "Numerical simulation of the anisotropic properties of a columnar jointed rock mass under triaxial compression," Engineering Computations, vol. 35 no. 4, pp. 1788-1804, DOI: 10.1108/EC-07-2017-0240, 2018.

[31] X. Wang, N. Wu, L. Hong, Y. Yan, "Influence of joint angle on the instability failure characteristics and AE evolution law of underground caverns," European Journal of Environmental and Civil Engineering,DOI: 10.1080/19648189.20201837252, 2020.

[32] C. Zhu, Z. G. Tao, S. Yang, S. Zhao, "V shaped gully method for controlling rockfall of high-steep slope in China," Bulletin of Engineering Geology and the Environment, vol. 78 no. 4, pp. 2731-2747, 2019.

[33] C. Jia, W. Xu, S. Wang, R. Wang, J. Yu, "Experimental analysis and modeling of the mechanical behavior of breccia lava in the dam foundation of the Baihetan hydropower project," Bulletin of Engineering Geology and the Environment, vol. 78 no. 4,DOI: 10.1007/s10064-018-1228-3, 2018.

[34] A. Shi, Y. Wei, J. Wu, D. Ren, M. Tang, "Study on the shear deformation of intralayer shear bands at the Baihetan hydropower station dam foundation," Bulletin of Engineering Geology and the Environment, vol. 79 no. 7,DOI: 10.1007/s10064-020-01799-7, 2020.

[35] C. Jin, S. Li, J. Liu, "Anisotropic mechanical behaviors of columnar jointed basalt under compression," Bulletin of Engineering Geology and the Environment, vol. 77 no. 1, pp. 317-330, DOI: 10.1007/s10064-016-0942-y, 2016.

[36] M. Chen, S.-Q. Yang, P. G. Ranjith, Y.-C. Zhang, "Cracking behavior of rock containing non-persistent joints with various joints inclinations," Theoretical and Applied Fracture Mechanics, vol. 109,DOI: 10.1016/j.tafmec.2020.102701, 2020.

[37] C. Zhou, M. Karakus, C. Xu, J. Shen, "A new damage model accounting the effect of joint orientation for the jointed rock mass," Arabian Journal of Geosciences, vol. 13 no. 7,DOI: 10.1007/s12517-020-5274-3, 2020.

[38] Y. Zhang, W. Lu, M. Chen, P. Yan, Y. Hu, "Dam foundation excavation techniques in China: a review," Journal of Rock Mechanics and Geotechnical Engineering, vol. 5 no. 6, pp. 460-467, DOI: 10.1016/j.jrmge.2013.08.002, 2013.

[39] X. Wang, Y. Liu, Z. Tao, W. Wang, Q. Yang, "Study on the failure process and nonlinear safety of high arch dam and foundation based on geomechanical model test," Engineering Failure Analysis, vol. 116,DOI: 10.1016/j.engfailanal.2020.104704, 2020.

[40] Y. Y. Luo, S. Z. Shen, "Application of controlled blasting technology in dam foundation excavation of Longtan hydropower station," Hongshui River, vol. 22 no. 3, pp. 43-45, 2003. in Chinese

[41] L. Zhao, M. Tao, G. Wu, "The excavation technique at the trenches of the dam abutment on left bank of the Xiluodu hydropower station on the Jinshajiang river," pp. 203-210, . in Chinese

[42] S. J. Guan, X. C. Yuan, "Excavation construction and supervision for dam foundation of Jinping dam," Yangtze River, vol. 40 no. 18, pp. 61-63, 2009. in Chinese

[43] S. Zhu, "Study on relaxation property of rock mass excavation on left bank dam foundation of Baihetan hydropower project," , 2019. Master’s thesis