# Influence of Parallel Joint Spacing and Rock Size on Rock Bulk Modulus

Authors: Tao Zhao; Gaojian Hu; Tao Wang; Huan Zhang
Journal: Advances in Civil Engineering
 (2022)

---

## Abstract

The size effect on the bulk modulus of rocks has been reported by previous studies. The accuracy in selecting the parameter for the rock mechanics analysis determines further accuracy of the calculation results. Moreover, given the influence of the size effect and joint spacing, the rock bulk modulus often changes. Thus, it is essential to examine the size effect on the bulk modulus. This study elucidated the influence of rock size and parallel joint spacing on the bulk modulus using the regression analysis and 12 sets of numerical plans. The results demonstrated that the bulk modulus decreased with an increase in rock size, and the curve represents an exponential function. The bulk modulus linearly increased as the parallel joint spacing increased. Furthermore, the characteristic size of the bulk modulus linearly decreased as the parallel joint spacing increased. In contrast, the characteristic bulk modulus linearly increased as the parallel joint spacing increased. The specific forms of these relationships were also elucidated in this study.


---

## Body

## 1. Introduction

The rock bulk modulus reflects the ability of a rock to resist deformation under external loads, using which one can characterize the mechanical characteristics of rocks. The bulk modulus is a comprehensive response of the rock mineral composition, fluid, pore, and structure under the action of the original internal environment, such as the formation pressure and temperature [1]. These unique characteristics make the bulk modulus an important parameter to study.The joints in a rock fundamentally affect its bulk modulus. For instance, parallel joints can exist in sedimentary rocks, thereby significantly affecting the bulk modulus [2]. At the same time, the number of parallel joint spacings (PJS) in the rock causes a difference in the bulk modulus. For instance, Zhao et al. [3] have used the finite element method to study the bulk modulus of rock. Moreover, some scholars have studied the effects of different factors on bulk modulus. For instance, Blake and Faulkner [4] have conducted seepage tests on granite to examine the effect of bulk modulus on fracturing permeability. Furthermore, Davarpanah et al. [5] studied the relationship between the modulus ratio and bulk modulus. Liu and Zhang [6] conducted an indoor triaxial creep test and found that the absolute value of the bulk modulus decreases with an increase in creep stress. Scholars have also addressed the calculation methods for bulk modulus. For instance, Shen et al. [7] studied the relationship between the tensile failure strain of rock and the initial fracture density spacing. They ultimately established a calculation method for the rock bulk modulus under tensile conditions. Zhao et al. [8] studied the failure and cracking characteristics of a broken rock mass cut using joints. Li et al. [9] obtained the volumetric fracture and crack propagation law of rocks through a dynamic load test. Overall, research literature extensively examined bulk modulus, but only a few studies have considered the size and PJS, while the relationship between bulk modulus and PJS is poorly understood.The rock has a side effect, and the change in rock size also affects the rock bulk modulus. For instance, the change in the size of some sedimentary rocks rich in joints fundamentally affects the change in the rock bulk modulus [10]. In this regard, Pang et al. [11] demonstrated that the bulk modulus of coal increases with a decrease in particle size. Jiang [12] reported that the size of the rock particles affects the rock pore structure, thereby triggering changes in the rock bulk modulus. Zhang and Yang [13] quantified the bulk modulus of montmorillonite by measuring changes in the length and density of the sample. Some other scholars have studied the effect of irregular particles on the bulk modulus. For instance, Kerimov et al. [14] have elucidated the effects of irregularly shaped particles and particle size on the porosity, permeability, and bulk modulus of granular porous media. Researchers also examined the effect of size on bulk modulus, but the elements of PJS, as well as the relationship between bulk modulus and size, are both poorly understood.The mechanical parameters of rocks fundamentally vary with the size of the jointed rock and eventually tend toward a stable value, which is defined as the representative essential volume (REV). Some researchers have previously proposed various methods to evaluate the REV, whereas Hu and Ma [15] used a realistic failure process analysis (RFPA) to study the characteristic size. Ying et al. [16] established a method based on volume rupture strength (P32) and a statistical test method to estimate rock mass REV. Wu et al. [17] studied the effect of size on bulk modulus, while taking into account the effect of the model location, and reported that the REV size was 18 m. Liu et al. [18] elucidated the size effect of the defective rock mass strength through uniaxial and confining pressure tests and obtained a REV size of 5 m × 10 m. Hu et al. [19] obtained the relationship between the characteristic size of the rock elastic modulus and the PJS. Overall, although the bulk modulus has been extensively studied, only a few researchers investigated the relationship between the characteristic size of the bulk modulus (CSBM) and PJS and established a model of the CSBM and PJS.This study established 12 numerical models to investigate the influence of the PJS and rock size on rockK. Within this research aim, (1) the corresponding stress-strain curves were analyzed, (2) the relationship between K and PJS, as well as the relationship between K and size were both established. Finally, (3) a model of the CSBM and PJS and a model of the rock characteristic bulk modulus (CBM) and PJS were established.


## 2. Numerical Simulation Plans

This study primarily focuses on two aspects: (1) elucidating the influence of PJS on rock K, with PJS of 10, 20, 30, 40, and 50 mm and (2) elucidating the effect of rock size with PJS on rock K, with the rock sizes of 100, 200, 400, 600, 800, 1,000, and 1,200 mm. Table 1 summarizes the research plans of this study based on [15].Table 1
Research plans [15].
Numerical simulationPJS (mm)Research plansPlan 1Plan 2Plan 3Plan 4Plan 5Plan 6Plan 7l = 100 mml = 200 mml = 400 mml = 600 mml = 800 mml = 1000 mml = 1200 mmPlan 8s = 1010 × 10010 × 20010 × 40010 × 60010 × 80010 × 100010 × 1200Plan 9s = 2020 × 10020 × 20020 × 40020 × 60020 × 80020 × 100020 × 1200Plan 10s = 3030 × 10030 × 20030 × 40030 × 60030 × 80030 × 100030 × 1200Plan 11s = 4040 × 10040 × 20040 × 40040 × 60040 × 80040 × 100040 × 1200Plan 12s = 5050 × 10050 × 20050 × 40050 × 60050 × 80050 × 100050 × 1200s is the parallel joint spacing, and l is the rock size.This study uses RFPA as the simulation software. The boundary conditions and rock mechanical parameters used in the numerical simulation are summarized in referred to [15].The roughness coefficient of the joint applied for the numerical simulations was 3, its elastic modulus was 1.5 MPa, its compressive strength was 2 MPa, its Poisson’s ratio was 0.3, and its internal friction angle was 30°.


## 3. Numerical Results and Analysis

### 3.1. Stress-Strain Curve Analysis

The research analysis included plotting of the stress-strain curves in the plans 1 to 7 (Figure1). Furthermore, the stress-strain curves in plans 8 to 12 were plotted as well (Figure 2).Figure 1
Stress-strain curves of rocks of different sizes. (a)l = 100 mm, (b) l = 200 mm, (c) l = 400 mm, (d) l = 600 mm, (e) l = 800 mm, (f) l = 1000 mm, and (g) l = 1200 mm.
(a)(b)(c)(d)(e)(f)(g)Figure 2
Stress-strain curves of rocks with different PJS. (a)s = 10 mm, (b) s = 20 mm, (c) s = 30 mm, (d) s = 40 mm, and (e) s = 50 mm.
(a)(b)(c)(d)(e)Figure1 illustrates the effect of the PJS on the compressive strength of rocks with different rock sizes. As seen in Figures 1(a)–1(g), the laws of the stress-strain curves were similar as the rock size increased from 100 to 1,200 mm. The strain of the rock gradually increased with the increase of the stress under the pressure effect. Moreover, it exhibited a linear elastic failure, which was further exacerbated by a plastic failure. Figure 1(a) is shown as an example, where the rock size is 100 mm. As seen, the compressive strength of the rock gradually increased, as the PJS increased from 10 to 50 mm.Figure2 displays the effect of the rock size on the compressive strength of the rock at different PJS. As seen from Figures 2(a)–2(e), when the PJS increased from 10 to 50 mm, the laws of rock stress-strain curves were also similar. Like in Figure 2, the rock was also destroyed by the elastic deformation first, which was further exacerbated by plastic deformation. Figure 2(e) is shown as an example with a PJS of 50 mm. As seen, the compressive strength of the rock gradually decreased with the increase of the rock size from 100 mm to 1,200 mm, thereby somewhat manifesting the size effect. Notably, the law of the curve is consistent with the law from [15], thereby indicating that the size effect of the rock was the same when the parameters of the joints changed.The elastic modulus fundamentally reflects the proportional relationship between the stress and strain of a material during the elastic deformation stage. Figures1 and 2 shows the obtained values of the elastic modulus and Poisson’s ratio for each working condition, as summarized in Tables 2 and 3.Table 2
Values of elastic modulus.
Numerical plansRock size (mm)Elastic modulus (GPa)Plan 8Plan 9Plan 10Plan 11Plan 12s = 10 mms = 20 mms = 30 mms = 40 mms = 50 mmPlan 1l = 1001.7502.9573.2363.5193.338Plan 2l = 2000.8641.5421.6811.7691.607Plan 3l = 4000.4380.7340.8590.9110.834Plan 4l = 6000.2900.4880.5630.5700.557Plan 5l = 8000.2180.3800.4060.4390.418Plan 6l = 10000.1720.2980.3270.3330.344Plan 7l = 12000.1450.2430.2720.2970.273Table 3
Values of Poisson’s ratio.
Numerical plansRock size (mm)Poisson’s ratioPlan 8Plan 9Plan 10Plan 11Plan 12s = 10s = 20s = 30s = 40s = 50Plan 11000.30000.20490.22000.23690.2930Plan 22000.23260.16670.14290.16540.2396Plan 34000.15230.20000.21110.27540.2861Plan 46000.17540.18750.18670.21180.2486Plan 58000.18750.19640.20970.21770.2638Plan 610000.10530.17500.23400.26890.2673Plan 712000.15860.16670.16070.19440.3013Bulk modulus is a relatively stable material constant. Fundamentally, there is a relationship between bulk modulus K, elastic modulusE, and Poisson’s ratio ν: K=E/3×1−2ν. Furthermore, the bulk modulus K of each working condition was solved according to this equation and the values from Tables 2 and 3, as shown in Table 4.Table 4
Values of bulk modulus.
Numerical plansRock size (mm)Bulk modulus (GPa)Plan 8Plan 9Plan 10Plan 11Plan 12s = 10s = 20s = 30s = 40s = 50Plan 11001.4581.6701.9262.2292.687Plan 22000.5380.7710.7840.8811.028Plan 34000.2100.4080.4950.6760.650Plan 46000.1490.2600.2990.3300.369Plan 58000.1170.2090.2330.2590.295Plan 610000.0730.1530.2050.2400.246Plan 712000.0710.1220.1340.1620.229


### 3.2. Influence of PJS onK

Table3 summarizes the statistical data, which revealed no clear relationship between Poisson’s ratio, rock size, and PJS, and the data obtained are relatively discrete.The statistical data, shown in Table4, demonstrate that K gradually increased as the PJS increased. A data point plot of K and PJS for each size was drawn, and the corresponding curve was fitted, as shown in Figure 3.Figure 3
Fitting curves of bulk modulus and PJS.As seen in Figure3, when the rock size was fixed and unchanged, K was affected by the PJS and increased with the increasing PJS. This pattern of variation remained the same even when the rock size varied. When the PJS was fixed and unchanged, K was affected by the rock size and decreased with an increase in the rock size. This finding indicates that K of the rock with parallel joints exhibited a positive correlation with the PJS and a negative correlation with the rock size. To illustrate this relationship in detail, see the formulae for the regression curves in Table 5.Table 5
Fitting relationships betweenK and PJS.
Rock size (mm)Fitting formulaFitting coefficient (R2)100Ks=0.027s+1.0890.976200Ks=0.013s+0.4730.927400Ks=0.012s+0.1430.907600Ks=0.005s+0.1280.917800Ks=0.004s+0.1000.9161000Ks=0.004s+0.0530.9071200Ks=0.004s+0.0360.942The fitting formula from Table5 was further used. The relationship indicated that K and PJS exhibited a linear relationship, and the mathematical model for K and PJS was formalized as(1)Ks=as+b,where K(s) is K when PJS is s (GPa) and s is PJS (mm); a and b are the parameters.The values of the parametersa and b from Table 5 are listed in Table 6 according to (1). The values from Table 6 were used to draw the fitting curves of the rock size and the parameters a and b (see Figure 4).Table 6
Values ofa and b under different rock sizes.
ParameterValuesl = 100 mml = 200 mml = 400 mml = 600 mml = 800 mml = 1000 mml = 1200 mmA0.0270.0130.0120.0050.0040.0040.004B1.0890.4730.1430.1280.1000.0530.036Figure 4
Fitting curve diagrams. (a) Parametera; (b) parameter b.
(a)(b)Figure4 shows that the parameters a and b exhibited a power function relationship to s. Thus, we established the following formulae:(2)a=1.241l−0.833,(3)b=381.501l−1.272.From equations (1)–(3), we obtained a special relational formula for K and PJS:(4)Ks=1.241sl−0.833+381.501l−1.272.(3) was used to quantify K, thereby providing a special relationship between K and PJS. Overall, it is applicable for solving K on a two-dimensional plane. In particular, for a known rock size, K can be obtained when the PJS was determined.


### 3.3. Influence of the Size Effect of K

A data point plot ofK and the rock size under each PJS was drawn according to the values from Table 4. The corresponding curve was fitted, as shown in Figure 5.Figure 5
Fitting curves of bulk modulus and rock size.Figure5 shows that when the PJS was fixed and unchanged, K was affected by the rock size and decreased with the increase in size. This pattern of variation remained the same even if PJS differed. Furthermore, when the rock size was fixed and unchanged, K was affected by rock size and increased accordingly. This phenomenon suggests that, for a rock with parallel joints, its K would have a positive correlation with the PJS, but a negative correlation with the rock size. For details of this relationship, see the formula of the regression curves summarized in Table 7.Table 7
Fitting relationships betweenK and rock size.
PJS (mm)Fitting formulaFitting coefficient (R2)10Kl=0.110+3.254e−0.0093l0.99520Kl=0.191+3.321e−0.0094l0.98830Kl=0.243+4.517e−0.0111l0.97940Kl=0.293+5.132e−0.0126l0.96350Kl=0.326+6.507e−0.0129l0.983The fitting formula in Table7 was used to infer an exponential relationship between K and the rock size. The mathematical model for K and rock size was proposed:(5)Kl=d+fe−gl,where K(l) is K at the size l and the unit is GPa; d, f, and g are parameters.The values of the parametersd, f, and g in Table 7 are listed in Table 8 according to (5). Furthermore, the values from Table 8 were used to draw the fitting curves of the rock size (d, f, and g are drawn), as shown in Figure 6.Table 8
Values ofd, f, and g.
ParametersValuess = 10 mms = 20 mms = 30 mms = 40 mms = 50 mmd0.10960.19050.24310.29270.3255f3.25413.32064.51745.13176.5067g0.00930.00940.01110.01260.0129Figure 6
Fitting curves of parameters and PJS. (a) Parameterd; (b) parameter f; (c) parameter g.
(a)(b)(c)Figure6 shows that the parameters d, f, and g all were linearly associated to s, thereby laying the foundation for the following formulae:(6)d=0.005s+0.072,(7)f=0.083s+2.051,(8)g=1.048×10−4s+0.008.Furthermore, equations (5)–(8) were used to obtain a special relational formula for K and the rock size:(9)Kl=0.0083s+2.051e−1.048×10−4s+0.008l+0.005s+0.072.(9) quantifies K, while also providing a special relationship between K and the rock size. Moreover, it can be used for solving K on a two-dimensional plane. Thus, K can be obtained when PJS is known and when the rock size is determined.


### 3.4. Relationships of CSBM, CBM, and PJS

#### 3.4.1. Derived Formula of CSBM

The size effect of K can be characterized by the characteristic size of the bulk modulus (CSBM). In particular, Liang et al. [20] previously provided the method for quantifying the characteristic size; one can solve the CSBM by(10)k=fge−gl,k≤γ,l≥lngf−lnγg,where γ is the acceptable absolute value of the inclination.


#### 3.4.2. Relationship of CSBM and PJS

The CSBM was solved, as summarized in Table9, when the PJS was 10, 20, 30, 40, and 50 mm. The regression curves for the CSBM and PJS are shown in Figure 7.Table 9
Relationship between CSBM and PJS.
PJS (mm)1020304050Characteristic size of bulk modulus (mm)690.66683.52622.58570.2574.1Figure 7
Fitting curve of CSBM and PJS.Figure7 shows that, as the PJS increased from 10 to 50 mm, the CSBM decreased from 690.66 mm to 574.1 mm. Moreover, a linear relationship between CSBM and PJS was discerned, where the slope of the curve was found to be negative. Thus, according to the fitting curve, the following specific relationship was formalized:(11)Bs=−3.464s+732.144,where B(s) is the characteristic size of bulk modulus (unit: mm).(11) quantifies the CSBM, thereby providing a special relationship between the CSBM and PJS, which can be used for solving the CSBM in a two-dimensional plane. Note that, in field applications, the CSBM can be generally obtained when the PJS is measured.


#### 3.4.3. Relationship of CBM and PJS

The value of CSBM was substituted into (5), and the characteristic bulk modulus (CBM) of rocks with different PJS values is described in Table 10.Table 10
Relationship between CBM and PJS.
PJS (mm)1020304050Characteristic bulk modulus (GPa)0.1150.1960.2480.2970.329Figure8 shows that the CBM increases gradually with the increase of the PJS, thereby exhibiting a linear relationship. On this basis, the special relation was obtained:(12)Kws=0.005s+0.078,where Kws is the CBM (unit: GPa).Figure 8
Fitting curve of CBM and PJS.(12) quantifies the CBM and provides a special relationship between the CBM and PJS. It can be used for solving the CBM on a two-dimensional plane. In field applications, the CBM can be obtained when the PJS is measured.


### 3.5. Verification Analysis

To evaluate the accuracy of our results, we analyzed the applicability of (5). The verification was conducted using the data of elastic moduli for different rock sizes from Figure 3.11 (page 44) in reference [21], as shown in Table 11. From these elastic modulus, the corresponding bulk modulus were calculated and summarized in Table 11. According to Table 11, we obtained a scatter plot and relationship curve of the bulk modulus and rock size (see Figure 9).Table 11
Elastic modulus and bulk modulus with different rock sizes.
Rock size2 m4 m8 m12 m16 mElastic modulus (GPa)11.388.067.667.397.29Bulk modulus (GPa)2.04841.45081.37881.33021.3122Figure 9
Fitting curves of bulk modulus with different sizes.The relationship between the bulk modulus and different sizes is shown in Figure9 and was formalized as follows:(13)Kl=1.34+4.37e−0.91l,where K(l) (GPa) is K of the rock when the rock size is l and l (m) is the rock size.The function type of (13) generally conformed to the mathematical model, introduced in (5). It can be therefore concluded that the numerical simulation and experimental conclusions were consistent. Moreover, the evaluation analysis indicates that the mathematical model, proposed in (5), is applicable for the solution of the bulk modulus with respect to the size.


## 3.1. Stress-Strain Curve Analysis

The research analysis included plotting of the stress-strain curves in the plans 1 to 7 (Figure1). Furthermore, the stress-strain curves in plans 8 to 12 were plotted as well (Figure 2).Figure 1
Stress-strain curves of rocks of different sizes. (a)l = 100 mm, (b) l = 200 mm, (c) l = 400 mm, (d) l = 600 mm, (e) l = 800 mm, (f) l = 1000 mm, and (g) l = 1200 mm.
(a)(b)(c)(d)(e)(f)(g)Figure 2
Stress-strain curves of rocks with different PJS. (a)s = 10 mm, (b) s = 20 mm, (c) s = 30 mm, (d) s = 40 mm, and (e) s = 50 mm.
(a)(b)(c)(d)(e)Figure1 illustrates the effect of the PJS on the compressive strength of rocks with different rock sizes. As seen in Figures 1(a)–1(g), the laws of the stress-strain curves were similar as the rock size increased from 100 to 1,200 mm. The strain of the rock gradually increased with the increase of the stress under the pressure effect. Moreover, it exhibited a linear elastic failure, which was further exacerbated by a plastic failure. Figure 1(a) is shown as an example, where the rock size is 100 mm. As seen, the compressive strength of the rock gradually increased, as the PJS increased from 10 to 50 mm.Figure2 displays the effect of the rock size on the compressive strength of the rock at different PJS. As seen from Figures 2(a)–2(e), when the PJS increased from 10 to 50 mm, the laws of rock stress-strain curves were also similar. Like in Figure 2, the rock was also destroyed by the elastic deformation first, which was further exacerbated by plastic deformation. Figure 2(e) is shown as an example with a PJS of 50 mm. As seen, the compressive strength of the rock gradually decreased with the increase of the rock size from 100 mm to 1,200 mm, thereby somewhat manifesting the size effect. Notably, the law of the curve is consistent with the law from [15], thereby indicating that the size effect of the rock was the same when the parameters of the joints changed.The elastic modulus fundamentally reflects the proportional relationship between the stress and strain of a material during the elastic deformation stage. Figures1 and 2 shows the obtained values of the elastic modulus and Poisson’s ratio for each working condition, as summarized in Tables 2 and 3.Table 2
Values of elastic modulus.
Numerical plansRock size (mm)Elastic modulus (GPa)Plan 8Plan 9Plan 10Plan 11Plan 12s = 10 mms = 20 mms = 30 mms = 40 mms = 50 mmPlan 1l = 1001.7502.9573.2363.5193.338Plan 2l = 2000.8641.5421.6811.7691.607Plan 3l = 4000.4380.7340.8590.9110.834Plan 4l = 6000.2900.4880.5630.5700.557Plan 5l = 8000.2180.3800.4060.4390.418Plan 6l = 10000.1720.2980.3270.3330.344Plan 7l = 12000.1450.2430.2720.2970.273Table 3
Values of Poisson’s ratio.
Numerical plansRock size (mm)Poisson’s ratioPlan 8Plan 9Plan 10Plan 11Plan 12s = 10s = 20s = 30s = 40s = 50Plan 11000.30000.20490.22000.23690.2930Plan 22000.23260.16670.14290.16540.2396Plan 34000.15230.20000.21110.27540.2861Plan 46000.17540.18750.18670.21180.2486Plan 58000.18750.19640.20970.21770.2638Plan 610000.10530.17500.23400.26890.2673Plan 712000.15860.16670.16070.19440.3013Bulk modulus is a relatively stable material constant. Fundamentally, there is a relationship between bulk modulus K, elastic modulusE, and Poisson’s ratio ν: K=E/3×1−2ν. Furthermore, the bulk modulus K of each working condition was solved according to this equation and the values from Tables 2 and 3, as shown in Table 4.Table 4
Values of bulk modulus.
Numerical plansRock size (mm)Bulk modulus (GPa)Plan 8Plan 9Plan 10Plan 11Plan 12s = 10s = 20s = 30s = 40s = 50Plan 11001.4581.6701.9262.2292.687Plan 22000.5380.7710.7840.8811.028Plan 34000.2100.4080.4950.6760.650Plan 46000.1490.2600.2990.3300.369Plan 58000.1170.2090.2330.2590.295Plan 610000.0730.1530.2050.2400.246Plan 712000.0710.1220.1340.1620.229


## 3.2. Influence of PJS onK

Table3 summarizes the statistical data, which revealed no clear relationship between Poisson’s ratio, rock size, and PJS, and the data obtained are relatively discrete.The statistical data, shown in Table4, demonstrate that K gradually increased as the PJS increased. A data point plot of K and PJS for each size was drawn, and the corresponding curve was fitted, as shown in Figure 3.Figure 3
Fitting curves of bulk modulus and PJS.As seen in Figure3, when the rock size was fixed and unchanged, K was affected by the PJS and increased with the increasing PJS. This pattern of variation remained the same even when the rock size varied. When the PJS was fixed and unchanged, K was affected by the rock size and decreased with an increase in the rock size. This finding indicates that K of the rock with parallel joints exhibited a positive correlation with the PJS and a negative correlation with the rock size. To illustrate this relationship in detail, see the formulae for the regression curves in Table 5.Table 5
Fitting relationships betweenK and PJS.
Rock size (mm)Fitting formulaFitting coefficient (R2)100Ks=0.027s+1.0890.976200Ks=0.013s+0.4730.927400Ks=0.012s+0.1430.907600Ks=0.005s+0.1280.917800Ks=0.004s+0.1000.9161000Ks=0.004s+0.0530.9071200Ks=0.004s+0.0360.942The fitting formula from Table5 was further used. The relationship indicated that K and PJS exhibited a linear relationship, and the mathematical model for K and PJS was formalized as(1)Ks=as+b,where K(s) is K when PJS is s (GPa) and s is PJS (mm); a and b are the parameters.The values of the parametersa and b from Table 5 are listed in Table 6 according to (1). The values from Table 6 were used to draw the fitting curves of the rock size and the parameters a and b (see Figure 4).Table 6
Values ofa and b under different rock sizes.
ParameterValuesl = 100 mml = 200 mml = 400 mml = 600 mml = 800 mml = 1000 mml = 1200 mmA0.0270.0130.0120.0050.0040.0040.004B1.0890.4730.1430.1280.1000.0530.036Figure 4
Fitting curve diagrams. (a) Parametera; (b) parameter b.
(a)(b)Figure4 shows that the parameters a and b exhibited a power function relationship to s. Thus, we established the following formulae:(2)a=1.241l−0.833,(3)b=381.501l−1.272.From equations (1)–(3), we obtained a special relational formula for K and PJS:(4)Ks=1.241sl−0.833+381.501l−1.272.(3) was used to quantify K, thereby providing a special relationship between K and PJS. Overall, it is applicable for solving K on a two-dimensional plane. In particular, for a known rock size, K can be obtained when the PJS was determined.


## 3.3. Influence of the Size Effect of K

A data point plot ofK and the rock size under each PJS was drawn according to the values from Table 4. The corresponding curve was fitted, as shown in Figure 5.Figure 5
Fitting curves of bulk modulus and rock size.Figure5 shows that when the PJS was fixed and unchanged, K was affected by the rock size and decreased with the increase in size. This pattern of variation remained the same even if PJS differed. Furthermore, when the rock size was fixed and unchanged, K was affected by rock size and increased accordingly. This phenomenon suggests that, for a rock with parallel joints, its K would have a positive correlation with the PJS, but a negative correlation with the rock size. For details of this relationship, see the formula of the regression curves summarized in Table 7.Table 7
Fitting relationships betweenK and rock size.
PJS (mm)Fitting formulaFitting coefficient (R2)10Kl=0.110+3.254e−0.0093l0.99520Kl=0.191+3.321e−0.0094l0.98830Kl=0.243+4.517e−0.0111l0.97940Kl=0.293+5.132e−0.0126l0.96350Kl=0.326+6.507e−0.0129l0.983The fitting formula in Table7 was used to infer an exponential relationship between K and the rock size. The mathematical model for K and rock size was proposed:(5)Kl=d+fe−gl,where K(l) is K at the size l and the unit is GPa; d, f, and g are parameters.The values of the parametersd, f, and g in Table 7 are listed in Table 8 according to (5). Furthermore, the values from Table 8 were used to draw the fitting curves of the rock size (d, f, and g are drawn), as shown in Figure 6.Table 8
Values ofd, f, and g.
ParametersValuess = 10 mms = 20 mms = 30 mms = 40 mms = 50 mmd0.10960.19050.24310.29270.3255f3.25413.32064.51745.13176.5067g0.00930.00940.01110.01260.0129Figure 6
Fitting curves of parameters and PJS. (a) Parameterd; (b) parameter f; (c) parameter g.
(a)(b)(c)Figure6 shows that the parameters d, f, and g all were linearly associated to s, thereby laying the foundation for the following formulae:(6)d=0.005s+0.072,(7)f=0.083s+2.051,(8)g=1.048×10−4s+0.008.Furthermore, equations (5)–(8) were used to obtain a special relational formula for K and the rock size:(9)Kl=0.0083s+2.051e−1.048×10−4s+0.008l+0.005s+0.072.(9) quantifies K, while also providing a special relationship between K and the rock size. Moreover, it can be used for solving K on a two-dimensional plane. Thus, K can be obtained when PJS is known and when the rock size is determined.


## 3.4. Relationships of CSBM, CBM, and PJS

### 3.4.1. Derived Formula of CSBM

The size effect of K can be characterized by the characteristic size of the bulk modulus (CSBM). In particular, Liang et al. [20] previously provided the method for quantifying the characteristic size; one can solve the CSBM by(10)k=fge−gl,k≤γ,l≥lngf−lnγg,where γ is the acceptable absolute value of the inclination.


### 3.4.2. Relationship of CSBM and PJS

The CSBM was solved, as summarized in Table9, when the PJS was 10, 20, 30, 40, and 50 mm. The regression curves for the CSBM and PJS are shown in Figure 7.Table 9
Relationship between CSBM and PJS.
PJS (mm)1020304050Characteristic size of bulk modulus (mm)690.66683.52622.58570.2574.1Figure 7
Fitting curve of CSBM and PJS.Figure7 shows that, as the PJS increased from 10 to 50 mm, the CSBM decreased from 690.66 mm to 574.1 mm. Moreover, a linear relationship between CSBM and PJS was discerned, where the slope of the curve was found to be negative. Thus, according to the fitting curve, the following specific relationship was formalized:(11)Bs=−3.464s+732.144,where B(s) is the characteristic size of bulk modulus (unit: mm).(11) quantifies the CSBM, thereby providing a special relationship between the CSBM and PJS, which can be used for solving the CSBM in a two-dimensional plane. Note that, in field applications, the CSBM can be generally obtained when the PJS is measured.


### 3.4.3. Relationship of CBM and PJS

The value of CSBM was substituted into (5), and the characteristic bulk modulus (CBM) of rocks with different PJS values is described in Table 10.Table 10
Relationship between CBM and PJS.
PJS (mm)1020304050Characteristic bulk modulus (GPa)0.1150.1960.2480.2970.329Figure8 shows that the CBM increases gradually with the increase of the PJS, thereby exhibiting a linear relationship. On this basis, the special relation was obtained:(12)Kws=0.005s+0.078,where Kws is the CBM (unit: GPa).Figure 8
Fitting curve of CBM and PJS.(12) quantifies the CBM and provides a special relationship between the CBM and PJS. It can be used for solving the CBM on a two-dimensional plane. In field applications, the CBM can be obtained when the PJS is measured.


## 3.4.1. Derived Formula of CSBM

The size effect of K can be characterized by the characteristic size of the bulk modulus (CSBM). In particular, Liang et al. [20] previously provided the method for quantifying the characteristic size; one can solve the CSBM by(10)k=fge−gl,k≤γ,l≥lngf−lnγg,where γ is the acceptable absolute value of the inclination.


## 3.4.2. Relationship of CSBM and PJS

The CSBM was solved, as summarized in Table9, when the PJS was 10, 20, 30, 40, and 50 mm. The regression curves for the CSBM and PJS are shown in Figure 7.Table 9
Relationship between CSBM and PJS.
PJS (mm)1020304050Characteristic size of bulk modulus (mm)690.66683.52622.58570.2574.1Figure 7
Fitting curve of CSBM and PJS.Figure7 shows that, as the PJS increased from 10 to 50 mm, the CSBM decreased from 690.66 mm to 574.1 mm. Moreover, a linear relationship between CSBM and PJS was discerned, where the slope of the curve was found to be negative. Thus, according to the fitting curve, the following specific relationship was formalized:(11)Bs=−3.464s+732.144,where B(s) is the characteristic size of bulk modulus (unit: mm).(11) quantifies the CSBM, thereby providing a special relationship between the CSBM and PJS, which can be used for solving the CSBM in a two-dimensional plane. Note that, in field applications, the CSBM can be generally obtained when the PJS is measured.


## 3.4.3. Relationship of CBM and PJS

The value of CSBM was substituted into (5), and the characteristic bulk modulus (CBM) of rocks with different PJS values is described in Table 10.Table 10
Relationship between CBM and PJS.
PJS (mm)1020304050Characteristic bulk modulus (GPa)0.1150.1960.2480.2970.329Figure8 shows that the CBM increases gradually with the increase of the PJS, thereby exhibiting a linear relationship. On this basis, the special relation was obtained:(12)Kws=0.005s+0.078,where Kws is the CBM (unit: GPa).Figure 8
Fitting curve of CBM and PJS.(12) quantifies the CBM and provides a special relationship between the CBM and PJS. It can be used for solving the CBM on a two-dimensional plane. In field applications, the CBM can be obtained when the PJS is measured.


## 3.5. Verification Analysis

To evaluate the accuracy of our results, we analyzed the applicability of (5). The verification was conducted using the data of elastic moduli for different rock sizes from Figure 3.11 (page 44) in reference [21], as shown in Table 11. From these elastic modulus, the corresponding bulk modulus were calculated and summarized in Table 11. According to Table 11, we obtained a scatter plot and relationship curve of the bulk modulus and rock size (see Figure 9).Table 11
Elastic modulus and bulk modulus with different rock sizes.
Rock size2 m4 m8 m12 m16 mElastic modulus (GPa)11.388.067.667.397.29Bulk modulus (GPa)2.04841.45081.37881.33021.3122Figure 9
Fitting curves of bulk modulus with different sizes.The relationship between the bulk modulus and different sizes is shown in Figure9 and was formalized as follows:(13)Kl=1.34+4.37e−0.91l,where K(l) (GPa) is K of the rock when the rock size is l and l (m) is the rock size.The function type of (13) generally conformed to the mathematical model, introduced in (5). It can be therefore concluded that the numerical simulation and experimental conclusions were consistent. Moreover, the evaluation analysis indicates that the mathematical model, proposed in (5), is applicable for the solution of the bulk modulus with respect to the size.


## 4. Discussion

TheK values of rocks of different sizes varied with the PJS. Importantly, this study established the following four relationships: (1) K and PJS; (2) K and rock size; (3) CSBM and PJS; as well as (4) CBM and PJS. Only a few previous studies considered the size effect of rocks with PJS on K. Moreover, the effect of rock size changes on the K of rocks with PJS, and the effect of the changes in PJS on the size effect of K; both remained understudied.In this study, theK values of rocks with different sizes were obtained using the data from [21], and (15) was obtained, which confirmed the accuracy of (5). Moreover, it proved that the formula, proposed in the study, stands out with certain universality, thereby also confirming the high accuracy of the research results. However, some differences between the numerical and experimental results were discerned. In addition, the two-dimensional model is limited in simulating three-dimensional samples. Both these challenges can be alleviated in future studies.Overall, the specific relationships of rockK were obtained in this study. The information about the rock size effect and the PJS is essential as input for the setting of mechanical parameters such as K in rock engineering, which can prevent rock engineering disasters. When the PJS and rock size are available, rock K, CSBM, and CBM can be quickly obtained for a selected engineering site, thereby providing valuable data-driven guidelines for engineering purposes.


## 5. Conclusions

This study elucidated the size effect of rocks with a PJS onK using numerical simulations. The following conclusions can be drawn:(1)
The relationship between K and PJS is linear.(14)Ks=as+b.We also obtained the special relationship.(15)Ks=1.241sl−0.833+381.501l−1.272.Moreover, the relationship betweenK and rock size is exponential:(16)Kl=d+fe−gl.And we obtained the special relationship.(17)Kl=0.0083s+2.051e−1.048×10−4s+0.008l+0.005s+0.072.(2)
CSBM is fundamentally associated with PJS. This study obtained the following relationship:(18)Bs=−3.464s+732.144.(3)
The CBM was found to be related to PJS. This study provided the following relationship:(19)Kws=0.005s+0.078.



---

*Source: 8475726-2022-09-10.xml*