Abstract

The mechanical properties and permeability evolution of sand-infilled rock joints during the shear process is an important issue in rock engineering, such as it pertains to hydraulic fractures filled with proppant. Shear can disrupt the preexisting hydraulic and mechanical equilibrium conditions, thus affecting fluid flow. In this study, we simulate the shear behavior of rock joints with variable roughness and sand infilling thickness using the discrete element code PFC2D. Rock joint roughness is evaluated by the joint roughness coefficient (JRC), and sand infilling thickness is evaluated by a thickness ratio (i.e., ratio of infill thickness to rock height) ranging from 0.02 to 0.20. The results show that peak shear strength decreases with the thickness ratio in a relation that can be expressed by a hyperbolic function. We also measure the permeability evolution during shearing and find that the permeability of infilled rock joints increases with both the thickness ratio and JRC.

1. Introduction

The mechanical and hydraulic properties of rock joints are of great significance in rock engineering, including near-field modeling of activated faults [13], petroleum, shale, and geothermal reservoirs [49]. Shear stimulation of existing fracture networks is considered a method for permeability enhancement in geothermal systems [1012]. Similar to fault zones (Figures 1(a) and 1(b)) with continuous infilled gouge layers [13, 14], weathered joints (Figure 1(c)) are generally filled with undeformed minerals, which can significantly influence the mechanical properties and permeability of rock masses [15]. The shear behavior and permeability of filled discontinuities in rock are determined by several parameters, including the infill thickness and joint surface roughness [1618]. Although the influence of joint roughness and infill thickness on rock mechanics has been studied previously, the coupled effect of these parameters on the mechanical behavior and permeability evolution of filled discontinuities is rarely reported.

Shear failure occurs when a frictional threshold value is reached. A sand-infilled rock joint can dramatically reduce the shear failure strength of a rock mass, and the infilled materials can dominate the shear behavior of intact rock due to its relatively low frictional properties [1922]. The infill thickness plays a vital role in estimating the mechanical behavior of a rock mass. Specifically, a thicker infill implies a smaller peak shear strength [23, 24]. Several models have been proposed to predict the peak shear strength of infilled joints by considering the ratio of infill thickness (t) to the height of an idealized (i.e., regular or planar) joint wall asperity (a) [2529]. In these models, shear strength decreases with t prior to reaching the critical value t/a. However, natural rock joints commonly have undulating surfaces, and the resistance of the roughness may pose a significant difference to that of a hypothetical smooth surface. Hence, it is difficult to estimate the shear strength of natural joints because a reasonable t/a ratio is not easily determined.

In order to characterize natural joint roughness, Barton [30] introduced the joint roughness coefficient (JRC), which varies from 0 to 20, corresponding to the smoothest to the roughest joint surfaces. Based on the JRC, Barton and Choubey [31] established an empirical law of rock joint friction to estimate the shear strength of rock joints. Numerous subsequent studies have been conducted to investigate shear strength criteria [3235], mechanical behavior [3639], and porosity and/or permeability evolution of infilled joints [4043]. Coupled shear-flow experiments of a single rock joint have also been performed to study the effect of joint shear deformation on the permeability evolution of rock joints [4452]. Specifically, Ye et al. conducted a series of novel shear stimulation experiments to probe the fundamental principle of permeability enhancement in rock joints. The experimental results demonstrate that dilatant shear deformation and cracks propagation are two crucial and integral mechanisms for permeability creation in rock joints.

Fluid flow through a single rough fracture varies with shear displacement. More specifically, the permeability of a rock mass decreases with loading stress and increases with shear-induced dilation. However, research on the coupled effects of JRCs and infill thickness on permeability evolution is underreported in the literature. In this study, we numerically simulate the shearing process of infilled rock joints using particle flow code (PFC) software. The paper is organized as follows: the distinct element numerical method and material parameters are introduced in Section 2, the simulation procedure is illustrated in Section 3, and results and discussion are presented in Section 4. Finally, several important conclusions are drawn in Section 5.

2. Numerical Method and Material Parameters

2.1. Description of DEM and PFC2D

The distinct element method (DEM) is based on solving Newton’s second law of motion for an assembly of particles with a predefined constitutive model applied to the grain-grain contacts [53]. Particle flow code 2D (PFC2D) is a two dimensional program based on DEM theory [54, 55]. In PFC2D, particles are idealized as rigid particles that interact with one another at their contacts.

The calculation cycle in PFC2D is a time-stepping algorithm, and the contact force is updated by the force-displacement law. When particles come into contact, the contact forces are calculated as a function of relative displacements and specified stiffness. The linear relationships can be expressed as follows: where represents the total normal force, is the normal stiffness, is the total normal displacement, is the unit normal vector to the contact plane, is the relative shear force, and and are the shear stiffness and relative shear displacement, respectively.

For sand-infilled discontinuities, we only consider the frictional parts of particles. Thus, in the direct shear numerical simulation, shear slip occurs when the iteratively checked contact force exceeds the maximum allowable shear contact force. where represents the friction coefficient at the contact.

2.2. Contact Bond Model

The numerically simulated specimens consist of two groups, intact rock and infill, represented as yellow and red particles, respectively, in Figures 2(a) and 2(b). An intact rock with joints is represented as a heterogeneous material of cemented grains. A real cement is present in sedimentary rock, while granular interlock in crystalline rock (e.g., granite) can be considered the notional cement. A parallel bond provides the mechanical behavior of a cement-like material substance between two contacting particles [56]. The parallel bond component acts in parallel with the linear component and establishes an elastic interaction between the pieces (Figure 2(c)). Parallel bonds, which can transmit both force and moment between pieces, have been successfully employed in estimating the mechanical behavior of rocks [5759].

2.3. Permeability Enhancement Model

A sigmoidal logarithmic function is employed to illustrate the permeability versus shear displacement curve [60]. where is the maximum increment in permeability, is the permeability change, is the shear displacement, is a constant related to the intrinsic property of the materials, and and are shear displacements at which 5% and 95% of the total permeability enhancement has occurred, respectively.

According to equation (3), the relative difference in permeability is

Combining equations (3) and (4), we obtain

Equation (4) can be rearranged as where is the initial permeability, is the shear strain, and are shear strains at which 5% and 95% of the total permeability has occurred, respectively, and equals 11.0.

2.4. Material Parameters

As shown in Figure 2(a), the numerically simulated granite specimen consists of upper and the lower halves and the infilled material (assumed as nonbonded to which a zero bond strength is assigned). The elastic behavior of the granular assembly is represented by assigning an effective modulus and normal-to-shear stiffness ratio to the grain-to-grain contacts. The average, maximum, and minimum ball radius were , , and , respectively. Each ball had at least three contacts with other neighboring balls, and balls without contacts were algorithmically identified as floaters. Different material parameters, including the normal-to-shear stiffness ratio, calibrated Young’s modulus, and uniaxial compressive strength (UCS) are listed in Table 1.

3. Numerical Simulation Design

The direct shear configuration consists of two sets of densely packed particles of nonuniform size and a gouge sandwiched between them (Figure 3(a)). The two sets represent the upper and lower rocks, respectively. In the DEM modelling, 4811 particles were packed into a rectangular box with a length of 60 mm and height of 40 mm. The profiles of the infilled assembly were attributed by fitting five typical joint roughness curves [31]. In the numerical simulation of the shear tests, five groups of specimens (i.e., T1 to T5) were compacted, and five typical JRCs (i.e., 0–2, 4–6, 8–10, 14–16, and 18–20) and five infill thickness ratios ranging from 0.02 to 0.2 were considered (Figure 3(b)). The simulation parameters are listed in Table 2. Figure 4 shows several assembled specimens with variable roughness and infill thickness.

The simulated specimens were compacted and tested as follows: (1) generation of rock specimens, (2) identification of the roughed infill with a certain thickness and JRC, (3) removal of floating particles, (4) application of a linear parallel contact bond to the intact rock and fill, (5) application of a constant normal stress of 20 MPa on the wall using the servo control algorithm in all numerical runs, (6) shearing of the upper half using a shear rate of 0.1 mm/s, and (7) termination of the shearing process when a shear strain of 0.10 is reached. During the shearing process, mechanical and physical parameters (e.g., shear stress and porosity) were recorded in real time.

4. Results and Discussion

Our numerical simulation results indicate that evolution of the shear strength and permeability of rock joints is strongly related to the infill thickness ratio (t/h) and JRC. Section 4.1 describes the effect of the thickness ratio on the shear behavior of infilled joints, and Section 4.2 demonstrates the coupled effects of the thickness ratio and JRC on the permeability evolution.

4.1. Effect of t/h on the Shear Behavior of Infilled Joints
4.1.1. Peak and Residual Shear Strength

Figure 5 illustrates shear stress-shear strain curves of infilled specimens with various thickness ratios and JRCs. Shear behavior near the peak shear strength is shown in the magnified plots, i.e., Figures 5(a-2)–5(e-2). Prior to the peak shear strength, elastic and yield stages increase quickly, and the shear stress-shear strain curves undergo a gradual decline after reaching the peak strength.

As can been seen in Figures 5(a-2)–5(e-2), a greater infill thickness ratio implies a higher peak shear strain. Similarly, the residual shear strength increases with the infill thickness ratio. However, the residual shear strengths are similar after a relative high thickness ratio is exceeded (e.g., ). The shear strength is therefore predominantly governed by the infill and trends toward a stable value [26].

Peak and residual shear strengths are plotted in Figures 6(a) and 6(b), respectively, and the simulated data are listed in Table 3. A higher thickness ratio clearly implies a smaller peak shear strength, which is consistent with previous studies [22, 23, 61]. A critical thickness of 0.05 divides the peak shear strength versus thickness ratio curves into two stages. Prior to the critical value, the peak shear strength drops sharply with increasing thickness ratio, while a slight decrease can be identified afterwards. Interestingly, a higher residual shear strength, corresponding to a greater thickness ratio, is observed before reaching a critical thickness ratio of 0.10.

4.1.2. Peak Shear Strength Model

Infilled material in joints can reduce the peak shear strength of rock masses. More specifically, a thicker infill implies a smaller shear strength [23, 24]. The relationship between shear strength and thickness ratio can be expressed by the following function [26]: where is the normalized shear strength of the infilled or unfilled rock joints, is the variation of the stress (MPa), is the normal load (MPa), is the thickness of the infill (m), and refers to the height of the rock (m) and is a constant. and are constants that depend on the normal load and surface roughness, and their magnitudes can be obtained by a regression of the simulated data.

According to the hyperbolic function, Indraratna et al. [24] developed a conceptual normalized shear strength model for infilled joints that can be expressed as where is the friction angle of intact rock, is the ratio of infill thickness to the critical thickness (), is the friction angle of the infilling material, and , , and are constants. Curves fitted according to equation (8) and key regressive parameters are presented in Figure 6(a). The mechanical parameters are listed in Table 3.

4.2. Coupled Effects of t/h and JRC on the Permeability Evolution of Infilled Joints
4.2.1. Permeability Evolution

Dilation or compaction usually occur during the direct shear test [31, 62] and can have a substantial effect on the permeability evolution of an infilled rock mass, which can be directly measured in physical experiments by monitoring change in specimen layer thickness. Although permeability and porosity are closely related, the evolution of porosity is difficult to observe during real experiments.

Assuming that the rock matrix is assumed impermeable and fluid only flows through the porous infills within the void spaces of the fracture. Permeability is normally enhanced during shear slip and dilation. To illustrate the relationship between permeability and shear strain, we employed an existing model to predict the permeability of infilled joints by iterating the derived porosity. The model proposed by Segall and Rice [63] using the test data from Marone et al. [46] indicates that permeability evolution can be estimated by either porosity or specimens thickness [64, 65]. An estimation of the permeability can be expressed as follows: where is the relative difference in permeability, is the relative difference in specimen thickness, and is the relative difference in porosity.

According to equations (9) and (10), the relative difference in permeability (for convenience, termed as permeability) can be obtained. The evolution of permeability is illustrated in Figure 7. It is clear that permeability increases approximately linearly with increasing shear strain before the yield stage is reached. However, permeability tends to maintain a stable value after the critical shear strain is attained. The thickness ratio has a key effect on the permeability of infilled rock joints. More specifically, permeability increases with the thickness ratio. The simulated data are included in Table 4.

As shown in Figure 7, permeability increases rapidly with shear strain before the yield stage is reached, during which the permeability trends to a constant value. The relation between the relative difference in permeability and shear displacement can be represented by a sigmoidal logarithmic function as shown in equation (6). The calibrated curves are plotted in Figure 7. The regressive curves show good agreement with the simulated results. In addition, the increasing stage, yield stage, and stable stage are also well illustrated.

4.2.2. Effects of t/h and JRC on Permeability Evolution

Figure 8 and Table 4 show the effects of the JRC and thickness ratio on the relative difference in permeability. It is clear that permeability increases with both t/h and JRC. In the permeability versus thickness ratio curves, permeability increases rapidly at first; however, the growth ratio then decreases gradually after the peak permeability is reached. This result is consistent with the findings by Marone et al. [46].

Permeability increases with both the thickness ratio and JRC, which may be explained by the following two reasons. Firstly, loose media in joints can be assumed as tectonically crushed rock materials or the product of the decomposition or joint weathering, and the bond strength between joint surfaces is zero [23, 24]. The particles of infilled materials are unbonded but can still resist one another through surface friction, which means that these particles would be rearranged throughout the entire test. This rearrangement is affected by the applied normal stress and joint roughness. The rearrangement activity would be more intense with a thicker infill, which may result in a larger porosity variation as particles mechanically degrade one another. Secondly, as can be seen from Figure 8(b), a higher JRC corresponding to greater permeability due to greater joint dilation can be expected with a relatively higher JRC.

5. Conclusions

In this paper, the coupled effects of JRC and sand-infilled thickness on the shear behavior and permeability evolution were studied by conducting numerical simulations using PFC2D. Five groups of JRCs were adopted to represent joints with different roughness, and five groups of infill thickness ratios were used to study the impact of infill thickness on the shear behavior and permeability evolution. Several important conclusions can be drawn. (1)Peak and residual shear strength were determined in the simulated direct shear tests, and the effects of the infill thickness are significant. The peak shear strength decreases with the thickness ratio in the form of a hyperbolic function. A critical thickness ratio of 0.05 divides the thickness ratio-peak shear strength curves into two stages, i.e., sharp drop stage and slight decrease stage. However, the residual shear strength increases with the thickness ratio before a critical thickness ratio is reached, after which point the residual shear strength trends toward a constant value(2)The permeability evolution during a direct shear test can be estimated using the relative porosity difference. Permeability is enhanced by increasing shear strain, and a sigmoidal logarithmic function can be employed to illustrate their relationship(3)The thickness ratio and JRC have coupled effects on the permeability evolution of infilled rock joints. More specifically, permeability increases with both the thickness ratio (t/h) and JRC. The permeability versus thickness ratio curves show that permeability increases significantly at first; however, the growth ratio decreases gradually when the thickness ratio reaches a critical value

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 is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

Hongwei Zhang contributed equally to this work.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (grant numbers 51604267, 51674242), the State Key Laboratory of Coal Resources and Safe Mining, CUMT (grant number SKLCRSM15X05), and the Fundamental Research Funds for the Central Universities (grant number 2017XKZD06).