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

College of Safety Science and Engineering, Liaoning Technical University, Huludao, Liaoning 125105, China School of Safety and Emergency Management Engineering, Taiyuan University of Technology, Taiyuan, Shanxi 030024, China Key Laboratory of Mine­ermodynamic Disasters and Control, Ministry of Education, Liaoning Technical University, Huludao, Liaoning 125105, China School of Civil Engineering, Shaoxing University, Shaoxing, Zhejiang 312000, China Key Laboratory of Rock Mechanics and Geohazards of Zhejiang Province, Shaoxing, Zhejiang 312000, China


Introduction
e rock bulk modulus re ects the ability of a rock to resist deformation under external loads, using which one can characterize the mechanical characteristics of rocks. e bulk modulus is a comprehensive response of the rock mineral composition, uid, pore, and structure under the action of the original internal environment, such as the formation pressure and temperature [1]. ese unique characteristics make the bulk modulus an important parameter to study. e joints in a rock fundamentally a ect its bulk modulus. For instance, parallel joints can exist in sedimentary rocks, thereby signi cantly a ecting the bulk modulus [2]. At the same time, the number of parallel joint spacings (PJS) in the rock causes a di erence in the bulk modulus. For instance, Zhao et al. [3] have used the nite element method to study the bulk modulus of rock.
Moreover, some scholars have studied the e ects of di erent factors on bulk modulus. For instance, Blake and Faulkner [4] have conducted seepage tests on granite to examine the e ect 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. ey 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. e 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. e 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.
is study established 12 numerical models to investigate the influence of the PJS and rock size on rock K. 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.

Numerical Simulation Plans
is 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].
is study uses RFPA as the simulation software. e boundary conditions and rock mechanical parameters used in the numerical simulation are summarized in referred to [15]. e 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°.

Stress-Strain Curve Analysis.
e research analysis included plotting of the stress-strain curves in the plans 1 to 7 ( Figure 1). Furthermore, the stress-strain curves in plans 8 to 12 were plotted as well ( Figure 2). Figure 1 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. e 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. Figure 2 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.
e elastic modulus fundamentally reflects the proportional relationship between the stress and strain of a material during the elastic deformation stage. Figures 1 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.
Bulk modulus is a relatively stable material constant. Fundamentally, there is a relationship between bulk modulus K, elastic modulus E, and Poisson's ratio ]: ). 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 3 summarizes the statistical data, which revealed no clear relationship between Poisson's ratio, rock size, and PJS, and the data obtained are relatively discrete.     l=100 mm×100 mm l=200 mm×200 mm l=400 mm×400 mm l=600 mm×600 mm l=800 mm×800 mm l=1000 mm×1000 mm l=1200 mm×1200 mm
As seen in Figure 3, when the rock size was fixed and unchanged, K was affected by the PJS and increased with the increasing PJS. is 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. is 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. e fitting formula from Table 5 was further used. e relationship indicated that K and PJS exhibited a linear relationship, and the mathematical model for K and PJS was formalized as where K(s) is K when PJS is s (GPa) and s is PJS (mm); a and b are the parameters. e values of the parameters a and b from Table 5 are listed in Table 6 according to (1). e values from Table 6 were used to draw the fitting curves of the rock size and the parameters a and b (see Figure 4). Figure 4 shows that the parameters a and b exhibited a power function relationship to s. us, we established the following formulae: From equations (1)-(3), we obtained a special relational formula for K and PJS: (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.

Influence of the Size Effect of K.
A data point plot of K and the rock size under each PJS was drawn according to the values from Table 4. e corresponding curve was fitted, as shown in Figure 5. Figure 5 shows that when the PJS was fixed and unchanged, K was affected by the rock size and decreased with the increase in size. is 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. is 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. e fitting formula in Table 7 was used to infer an exponential relationship between K and the rock size. e mathematical model for K and rock size was proposed: where K(l) is K at the size l and the unit is GPa; d, f, and g are parameters. e values of the parameters d, 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. Figure 6 shows that the parameters d, f, and g all were linearly associated to s, thereby laying the foundation for the following formulae: g � 1.048 × 10 − 4 s + 0.008.
Furthermore, equations (5)-(8) were used to obtain a special relational formula for K and the rock size: (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. us, K can be obtained when PJS is known and when the rock size is determined.

Derived Formula of CSBM.
e 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 where c is the acceptable absolute value of the inclination.

Relationship of CSBM and PJS.
e CSBM was solved, as summarized in Table 9, when the PJS was 10, 20, 30, 40, and 50 mm. e regression curves for the CSBM and PJS are shown in Figure 7. Figure 7 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 6 Advances in Civil Engineering negative. us, according to the fitting curve, the following specific relationship was formalized: 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.

Relationship of CBM and PJS.
e value of CSBM was substituted into (5), and the characteristic bulk modulus (CBM) of rocks with different PJS values is described in Table 10.
l=100 mm discrete point l=200 mm discrete point l=400 mm discrete point l=600 mm discrete point l=800 mm discrete point l=1000 mm discrete point l=1200 mm discrete point l=100 mm curve fitting l=200 mm curve fitting l=400 mm curve fitting l=600 mm curve fitting l=800 mm curve fitting l=1000 mm curve fitting l=1200 mm curve fitting   Advances in Civil Engineering Figure 8 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: where K w (s) is the CBM (unit: GPa). (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.

Verification Analysis.
To evaluate the accuracy of our results, we analyzed the applicability of (5). e 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). e relationship between the bulk modulus and different sizes is shown in Figure 9 and was formalized as follows: where K(l) (GPa) is K of the rock when the rock size is l and l (m) is the rock size. e 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.

Discussion
e K 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, the K 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 rock K were obtained in this study. e 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.

Conclusions
is study elucidated the size effect of rocks with a PJS on K using numerical simulations. e following conclusions can be drawn: (1) e relationship between K and PJS is linear.    We also obtained the special relationship. K(s) � 1.241sl − 0.833 + 381.501l − 1.272 . (15) Moreover, the relationship between K and rock size is exponential: And we obtained the special relationship.
(3) e CBM was found to be related to PJS. is study provided the following relationship: K w (s) � 0.005s + 0.078.
Data Availability e data supporting the findings of this study are available from the corresponding author upon reasonable request.

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

Authors' Contributions
G Hu contributed to conceptualization, software, and funding acquisition. T Zhao contributed to methodology, data curation, investigation, writing, and original draft preparation. T Wang contributed to data curation and formal analysis. H Zhang was responsible for funding acquisition and formal analysis.