Abstract

Hydromechanical coupling in rock masses is an important issue for many rock mechanics and hydrogeology applications. The change of a water-bearing state will induce the fracture of the intact rocks and further accelerate the shear slip instability of the sheared surface. To investigate the weakening effect of water content on the mechanical properties of a rock mass, laboratory direct shear tests combined with three-dimensional analysis of sheared surfaces were carried out on sandstone samples with different water contents. The variogram parameters, sill and range, were applied to quantify the morphology of shear fracture surfaces, to reflect the shear failure process of the intact rock, and to provide a basis for resliding instability of jointed rock. It was determined that the sill represents the height of the fluctuation body in the fracture surface and the range represents the single fluctuation body and may reflect the frequency of fluctuations. The test results revealed that the increase in water content had a clear weakening effect on the shear strength and deformation behavior of rock, especially under saturated conditions. Moreover, the distribution of water in the samples directly affected the crack initiation and propagation and characteristics of the fracture morphology.

1. Introduction

Hydromechanical coupling in rock masses is an important issue for many rock mechanics and hydrogeology applications [17]. In mining and civil engineering, the redistribution of the stress field during the excavation of tunnels and underground caverns leads to the initiation and propagation of cracks, which can cause dramatic changes in the permeability of a rock mass [8, 9]. Hence, it is very important in engineering practice to study the failure and instability mechanisms of rock masses under the hydromechanical coupling action of groundwater [10]. Xiao et al. evaluated the equivalent porous media (EPM) and the relationship between force and deformation, and the flow rate of fluid through a fracture is especially critical [11]. Heiland and Raab reported that the permeability evolution of low-porosity sandstone is similar to that of crystalline rock under triaxial compression [12]. In the prefracture instability region, permeability decreases with compression stress and increases after expansion. When the volume of the specimen returns to the initial state, the permeability does not recover. Wang et al. reported that the pore pressure gradient field has an important influence on the propagation direction of the main tensile fracture [13]. Geertsma found that an increase in local pore water pressure can promote crack growth, while the overall increase in pore water pressure may inhibit the crack growth, especially for tensile cracks [14]. In the field of fluid injection research, it is considered that pore water pressure may induce fault slip and seismic events [1517]. However, the hydromechanical coupling mechanisms still are poorly understood [18, 19].

It is worth noting that the effect of water content, which is the intermediate link of the rock mass changing from a dry state to a pore water pressure state, often has been neglected. It is well known that an increase in water content will reduce the strength of a rock mass and will shorten the time before instability [20]. Meanwhile, the failure of rock is quite commonly induced by compressive and shear coupling loading [2126]. Xu et al. reported that the tensile and shear damage of the weakest parts of a rock slope are the main factors that trigger failure of the rock slope [27]. Liu et al. developed a new mesoscopic shear apparatus and studied the crack propagation path and damage mechanisms under compression-shear loading [28, 29]. Thus, it is very meaningful to study the shear failure process of rocks with different water contents.

The aim of the present work was to study the shear instability process and secondary slip risk after shear failure, based on the evolution of water infiltration into the rock. Compression-shear tests of sandstone with different water contents were carried out. The effect of water content on the mechanical behavior and on the morphological characteristics of the sheared fracture surface was analyzed. To further assess the risk of continued slip instability after intact rock sheared, the variogram function was applied to quantify the anisotropy of postshearing 3D geometry of the fracture surface.

2. Experimental Procedure

2.1. Sample Preparation

The sandstone samples were obtained from the Three Gorges Reservoir, Chongqing, China. The samples were comprised primarily of quartz, feldspar, chert, and muscovite with a grain size distribution of 0.1 to 0.5 mm. Drilling cores without obvious fractures were selected and were cut into cubes with dimensions of approximately . Young’s modulus of the samples was 6.79 GPa, the Poisson ratio was 0.26, the uniaxial compressive strength was 81.04 MPa, and the density was 2.32 g/cm3. A picture of some of the sandstone samples is presented in Figure 1.

2.2. Test Design

The tests were divided into three groups, which were taken from the same intact rock block with the same normal stress and identical stress paths. Three water contents (0%, 50%, and 100%) were selected to be tested, and the detailed water content levels were determined in the following manner [30]: (1)Firstly, a sample was dried in a drying oven at 105°C for 48 hours. It then was cooled to room temperature in a dry vessel. The water content () of the sample then was considered to be 0%(2)In the saturation process, a dry sample was immersed in a tank filled with purified water for 1 hour in a vacuum and then was weighed. The operation was repeated every hour until the weight of the sample remained unchanged. The sample was considered to be saturated () in this study(3)A saturated sample was placed in a laboratory with excellent ventilation. The sample was weighed every hour. The mass change values were used to calculate the water content according to where (%) is the water content of the specimen and (g) and (g) are the wet and dry masses of the sample, respectively.

Because of inaccuracy in these measurements, a water content value between 49% and 51% was assumed to be 50%. Figure 2 shows the loading status of the specimens. In the figure, is a load to produce normal stress and is a load to produce shear stress. The normal stress () that was applied to samples was considered 2 MPa. To digitize the morphology of the surface, the system projects a series of structured white light fringe patterns onto the surface of the object (see Figure 3). Two CCD cameras were employed to capture images of these patterns automatically. From these image pairs, high-density 3D point cloud data from the surface were calculated accurately [31]. Finally, digitized surface data were reconstructed at an interval of 0.2 mm.

3. Results

3.1. Shear Mechanical Behavior

In view of the consistent trend of the shear stress-shear strain curves under different conditions [30, 33], curves under saturated condition () (Figure 4(a)) were analyzed. In the present study, a statistical comparison and analysis of the mechanical parameters were made at the peak shear stress under different conditions (Figures 4(b)4(d)). It can be observed that as the shear strain increased during the initial stage, the shear stress and the normal strain increased in concavity form. This was due to the fact that the normal strain during the initial stage of loading was mainly associated with closure of the microcracks, rather than to dilation resulting from microcracking and microfracture closure [34]. Then, the shear stress continued to rise after a period of fluctuation until shear failure occurred. With increasing water content, the mechanical parameters at the peak shear stress showed a decreasing trend (Figures 4(b)4(d)). It should be concluded that the weakening effect of water on the mechanical properties of rock includes physical chemistry phenomena such as exfoliation, hydration, solution, oxidation, and abrasion [35] and mechanical action such as pore pressure [13]. When the water content reached 50%, the shear stress decreased by about 6%. However, when the water content reached 100%, the peak shear stress decreased by 20%. This indicated that the effect of mechanical action on shear strength was greater and the action of water physical chemistry just lowers the energy required for microcracking [36], whereas pore pressure reduces the effective stress. Under the condition where the sample was saturated with water, this resulted locally in the formation of pore water pressure in the sample and promoted the further initiation and expansion of microcracks, thus greatly reducing the peak shear stress because the seepage velocity of the water was low.

Meanwhile, peak shear strain increased when the water content reached 50% and decreased when the water content reached 100%. Peak normal strain decreased as the water content was increased. This further indicated that water softened the rock and pore pressure reduced the fracture toughness of tension cracks [13]. When the water content reached 100%, the tension cracks propagated further under the combined action of the local pore water pressure and weakening of the fracture toughness of the specimens without further expansion deformation which resulted in a significant decrease in the peak shear strain and peak normal strain.

This section may be divided into subheadings. It should provide a concise and precise description of the experimental results and their interpretation as well as the experimental conclusions that can be drawn.

3.2. Anisotropic Characteristics of the Shear Fracture Surface

The morphology of the fracture surface not only reflects the failure process [37] but also plays a dominant role in the mechanical and hydraulic behavior of discontinuous rock masses. The mechanical responses of a rock joint, such as its peak shear strength, are all dependent on the shear direction [38, 39]. Geostatistical methods are used to describe the anisotropy of joint surfaces because they are designed to characterize spatial variations [40]. In the present study, the variogram was employed to quantify the anisotropic characteristics of the shear fracture and was defined as where is the number of pairs of data whose lag is , is the height at point , and is the height at a radial distance in a direction from .

Here, a spherical model was used to fit the variogram function (Figure 5) and was defined to be where is the nugget effect, is the sill, is the range, and is the lag distance [41].

For the denser data, the nugget effect was zero [42]. The parameters of the variogram used in the present paper include the range and the sill.

The parameter distribution curve for the variogram at a water content of 100% is presented in Figure 6. It can be observed that the fracture surface is undulating along the shear direction, but the fracture face is relatively flat along the vertical direction of the shear (Figures 6(a) and 6(b)). Correspondingly, the maximum sill is along the shear direction and the maximum of range is perpendicular to the shear direction (Figures 6(c) and 6(d)). Therefore, it follows that the sill represents the height of the fluctuation body in a fracture surface and the range represents the single fluctuation body and can reflect the frequency of fluctuation.

4. Discussion

The morphology of the fracture surface not only reflects the failure process [37] but also plays a dominant role in the mechanical and hydraulic behavior of discontinuous rock masses. In the literature, the morphology of the fracture surface was measured based mainly on 2D and 3D parameters. The 2D parameters are mainly concentrated in the JRC proposed by Barton and Choubey [43]. The 3D parameters are mainly concentrated in the fractal dimension [4446]. However, the 2D parameters are limited by excessive data [47], while the 3D parameters are limited by poor directionality [48]. Wang et al. [39] proposed that based on the morphology of sheared fracture surfaces, the process of crack initiation and propagation during shear failures can be obtained and pointed that two parameters ( and ) in the variogram were effectively to quantify the anisotropy of fracture surfaces: the sill () can reflect the height of the fluctuation body in the fracture surface and the range () can reflect the single fluctuation body and its fluctuation frequency.

To quantify the anisotropy of the shear fracture surface, the anisotropic coefficient () was adopted [49], which is defined as the ratio of physical and mechanical parameters in the weakest direction to those in the strongest direction. Curves of the parameters of the variogram function of the shear fracture surface under different water contents are shown in Figure 7. It can be seen that both the maximum values of the sill and range parameters change with the same trend, decreasing when the water content was 50% and increasing when the water was 100% (see Figures 7(a) and 7(b)). This was due to the fact that the sample had a uniform internal stress distribution under dry or saturated conditions and the crack propagation path was consistent with the stress direction. However, when the water content was 50%, under compression-shear stress, the water in the sample penetrated from the interior into the surrounding area. Due to restriction of the seepage channels, the water content is too high, or a higher pore water pressure is generated locally, which results in an uneven distribution of the stress field inside the sample. In consequence, crack propagation at the center of the sample is more complicated than that in the surrounding area and leads to a decrease in the maximum values of sill and range. Similarly, when the water content was 0% and 100% (see Figures 7(c) and 7(d)), the directionality of the shear fracture surface parameters was more evident, resulting in a reduction in the anisotropy coefficient ( and ) and more obvious anisotropy characteristics.

5. Conclusions

The morphology of a fracture surface not only reflects the failure process but also plays a dominant role in the mechanical and hydraulic behavior of discontinuous rock masses. In the present investigation, direct shear tests were carried out on sandstone to investigate the effect of water content on the shear failure mechanism of the rock and the anisotropy characteristics of the shear fracture surface. The morphology of the surface was digitized using a 3D scanner, and the anisotropy of the shear fracture surface was calculated by employing the variogram. The following conclusions were obtained: (1)According to the results of the present study, the increase in water content has a significant weakening effect on the shear strength and deformation of rock. When the water content is increased to 50%, the shear stress was decreased by about 6%. However, when the water content reaches 100%, the peak shear stress decreased by 20%(2)The weakening effect of water on the mechanical properties of rock included physical chemistry changes (exfoliation, hydration, solution, oxidation, and abrasion) and mechanical action (due to pore pressure). The effect of mechanical action on shear strength was greater(3)It is reasonable that the variogram can be applied to describe the anisotropic characteristics of the shear fracture surface. When the sample was in an unsaturated state, the water content will redistribute under external load, which directly affects the crack propagation path and hence the morphology of the fracture surface(4)The relationship between the weakening effect of the water content on the mechanical properties of the rock and its fracture surface morphology was inconsistent. However, both the mechanical properties and fracture surface morphology are related to the microcracking mechanism

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by the National Natural Science Foundation of China (grant number 51704187); Shandong Province Key Research and Development Plan, China (grant number 2019GSF109068); Shandong Natural Science Foundation, China (grant numbers ZR2017BEE054 and ZR2019QEE029); Postdoctoral Innovation Project of Shandong Province, China (grant number 201703024); and Postdoctoral Application Research Project of Qingdao, China (grant number 2017-267). The authors are thankful for all of the support for this basic research.