3.1 Test results of mudstone-like material
The physical, mechanical, and hydraulic properties of 16 groups of specimens are presented in Table 8. The results suggest that the material's mineral composition was the same as that of natural Cretaceous mudstone, and the synthesis process was a reaction between minerals and water that did not cause chemical changes in the sample, which has been confirmed by the SEM tests. Additionally, the material’s water content (\(\:\omega\:\)), density (\(\:\rho\:\)), uniaxial compressive strength (UCS), elastic modulus(E), cohesion and angle of internal friction (c and \(\:\phi\:\)), and disintegration time (t) are mainly distributed in the ranges of 5.65–12.59%, 1.62–2.08 g/cm3, 0.09–3.59 MPa, 24.71–505.93 MPa, 97.43–495.32 kPa, 6.28–64.26° and 21–90 mi. While the distributions of natural Cretaceous mudstones listed in Section 2.1 are 2.73–10.53%, 1.97–2.16 g/cm3, 1.21–5.79 MPa, 131–687 MPa, 450–980 kPa, 9.26–27.82° and 10–40 min respectively. Figure 5 illustrates a detailed comparison of the range of distributions of these parameters between the two materials. It shows that the physical, mechanical and hydraulic parameter distributions of the newly-prepared mudstone-like materials have an apparent overlap with those of the natural mudstone, indicating that this new material can better meet the requirements of simulating natural mudstone. Furthermore, according to the demands of specific projects, we can flexibly adjust the mix proportions of base materials to customize the mudstone-like materials for stress-seepage coupling model test under certain similarity scales.
Table 8
Physical and mechanical properties test results of the mudstone-like material.
Group | | Major mineral content/% | Physical property | strength parameters | t/min |
Gypsum | Kaolinite | Rectorite | Quartz | \(\:\omega\:\)/% | \(\:\rho\:\)/(g.cm-3) | UCS/ MPa | E/MPa | c/kPa | \(\:\phi\:\)/(\(\:^\circ\:\)) |
1 | 12.80 | 1.95 | 1.40 | 68.45 | 8.70 | 1.62 | 0.19 | 96.51 | 189.43 | 27.56 | 78 |
2 | 9.60 | 3.9 | 2.81 | 56.61 | 10.87 | 1.72 | 0.15 | 42.13 | 139.76 | 42.83 | 43 |
3 | 6.40 | 5.86 | 4.21 | 50.24 | 9.80 | 1.89 | 0.09 | 32.01 | 217.88 | 20.06 | 50 |
4 | 3.20 | 7.81 | 5.62 | 53.01 | 9.41 | 1.82 | 0.12 | 24.71 | 97.43 | 32.78 | 21 |
5 | 19.20 | 2.93 | 2.11 | 51.15 | 11.01 | 1.68 | 1.55 | 144.51 | 322.76 | 54.28 | 90 |
6 | 14.40 | 5.86 | 4.21 | 57.17 | 9.02 | 1.62 | 0.90 | 70.12 | 209.29 | 54.51 | 75 |
7 | 9.60 | 8.78 | 6.32 | 23.96 | 5.65 | 2.08 | 2.23 | 187.19 | 129.54 | 44.32 | 54 |
8 | 4.80 | 11.71 | 8.43 | 39.95 | 6.48 | 2.04 | 1.46 | 113.43 | 167.65 | 30.41 | 46 |
9 | 25.60 | 3.9 | 2.81 | 42.84 | 12.35 | 1.66 | 1.98 | 255.13 | 417.49 | 29.34 | 85 |
10 | 19.20 | 7.81 | 5.62 | 31.52 | 8.67 | 1.82 | 2.61 | 319.33 | 492.75 | 6.28 | 73 |
11 | 12.80 | 11.71 | 8.43 | 50.67 | 8.51 | 1.74 | 2.01 | 176.82 | 220.58 | 37.60 | 61 |
12 | 6.40 | 15.62 | 11.24 | 41.07 | 7.76 | 1.91 | 2.21 | 207.51 | 166.38 | 41.80 | 58 |
13 | 32.00 | 4.88 | 3.51 | 21.69 | 12.59 | 1.83 | 3.59 | 505.93 | 442.67 | 64.26 | 82 |
14 | 24.00 | 9.76 | 7.02 | 25.67 | 11.81 | 1.83 | 2.95 | 389.98 | 495.32 | 54.65 | 75 |
15 | 16.00 | 14.64 | 10.53 | 38.27 | 9.85 | 1.73 | 2.21 | 217.83 | 226.23 | 38.75 | 68 |
16 | 8.00 | 19.52 | 14.04 | 41.67 | 11.20 | 1.82 | 1.73 | 146.90 | 213.93 | 9.69 | 62 |
In conclusion, the mudstone–like materials in this study are consistent well with the natural Cretaceous mudstones in terms of mineral composition, water content (ω), density (ρ), uniaxial compressive strength (UCS), elastic modulus(E), cohesion and angle of internal friction (c and φ), and disintegration time (t), proving that the prepare method of this kind of mudstone–like materials is quite practicable.
3.2 Sensitivity Analysis
The aforementioned experimental studies have demonstrated that the newly synthetic mudstone-like materials can accurately simulate the disintegration property of natural rock and instability mechanisms of coal mine roadways under long-term water-induced strength deterioration. However, the complexity of the in-situ geotechnical engineering disaster mechanism may be beyond one's imagination. Therefore, to improve its applicability to more practical geotechnical engineering problems, sensitivity analyses were carried out to identify the main controlling factors affecting the physical, mechanical and hydraulic properties of this mudstone-like material.
3.2.1 Analysis of factors affecting the density
Based on orthogonal experimental results, Fig. 6 presents the range R and impacts of 4 factors on density of the specimens. It can be observed that the influence degree of each factor on the density of mudstone-like materials is B ≈ C > D > A. This implies that the clay minerals and iron sand content exert a more significant influence on the density of the materials than that of the binder and powder sand content. Specifically, the density tends to increase linearly with increasing clay-gypsum ratio and iron sand content, while linearly decreasing with the increase of powder sand content. In addition, the density was least affected by the aggregate-binder ratio, fluctuating with changes in this influencing factor.
3.2.2 Analysis of factors affecting UCS
Based on orthogonal experimental results, Fig. 7 exhibits the range R and impacts of 4 factors on UCS of the specimens. We can obtain that the sensitivity extent of each factor on the UCS is A > C > D ≈ B, suggesting that the aggregate-binder ratio exerts a more significant influence on the UCS of mudstone-like materials than that of the other three factors. Specifically, UCS grows rapidly with decreasing aggregate-binder ratio but decreases slightly with increasing clay mixtures and powder content. In addition, UCS increases gradually with increasing iron sand content, but the increase is clearly smaller than that of UCS with decreasing aggregate-binder ratio.
3.2.3 Analysis of factors affecting E
Based on orthogonal experimental results, Fig. 8 illustrates the range R and impacts of 4 factors on the elasticity modulus (E) of the specimens. It can be seen that the sensitivity rank of the four factors on the E is A > C > B > D, demonstrating that the aggregate-binder ratio plays a more dominant role on the E of mudstone-like materials than that of the other three factors. Specifically, the E increases sharply and linearly with decreasing aggregate-binder ratio but decreases slightly with increasing clay and powder content. Moreover, with increasing iron sand content, the E value also increases, but significantly less than the E that increases with decreasing aggregate-binder ratios.
3.2.4 Analysis of factors affecting cohesion
Based on orthogonal experimental results, Fig. 9 shows the range R and impacts of 4 factors on cohesion of the specimens. It can be concluded that the influence degree of each factor on the cohesion of mudstone-like materials is A ≈ B > C > D, which reveals that both the aggregate-binder ratio and clay-gypsum ratio exert a more obvious influence on the cohesion of mudstone-like materials than that of the content of iron sand and powder sand. However, the trend of the effect of Factor A and B on cohesion is diametrically opposed. Specifically, the cohesion rises quickly with decreasing of aggregate-binder ratio, whereas it decreases rapidly with increasing clay minerals content, demonstrating a clear deterioration effect. Additionally, iron and chalk content had a small effect on cohesion, and there was no significant linear relationship between the two factors and cohesion.
3.2.5 Analysis of factors affecting internal friction angle
Based on orthogonal experimental results, Fig. 10 presents the range R and impacts of 4 factors on internal friction angle of the specimens. The conclusion can be drawn that the order of sensitivity of each factor on the internal friction angle is D > A > B > C, implying that the powder sand content exerts a more significant influence on the internal friction angle of mudstone-like materials than that of the other three factors. Nevertheless, the internal friction angle changed irregularly as the powder sand content increased, indicating no obvious relationship between them, and the same holds true for Factor A and C. Moreover, with increasing clay content, the internal friction angle tended to decrease linearly, which could be explained by the addition of clay minerals weakening the frictional characteristics of the mudstone-like materials.
3.2.6 Analysis of factors affecting disintegration time
Based on orthogonal experimental results, Fig. 11 depicts the range R and impacts of 4 factors on the complete disintegration consumption time (t) of the specimens. The result shows that the sensitivity rank of the four factors on the complete disintegration consumption time is B > A > C > D, suggesting that the clay-gypsum ratio plays a more major part in the disintegration time of mudstone-like materials than that of the other factors. Specifically, the disintegration time decreases rapidly and linearly with increasing clay-gypsum ratio while increases moderately with decreasing aggregate-binder ratio. This finding suggests that the presence of clay mineral additions, especially swelling clay minerals, leads to the disintegration of the material in water. Moreover, with increasing iron and powder sand content, the disintegration time fluctuated within a small range, indicating that the parameter is less sensitive to these two factors.
In summary, the findings of the range analysis showed that the aggregate-binder ratio significantly impacted the UCS, E and cohesion, while the clay-gypsum ratio content strongly influenced the density and disintegration time of mudstone-like materials. From the mineralogical point of view, gypsum significantly enhances the material's mechanical properties, i.e., as the aggregate-binder ratio decreases, the UCS, E, and cohesion of the material show a significant increase. In contrast, clay minerals have a remarkable deterioration effect on the material's mechanical properties, i.e., as its content rises, all the strength parameters display the same linear decreasing trend. Besides, the disintegration property of the materials also increased linearly with the increase in clay minerals. These findings imply that the mineralogy of the materials has a significant effect on the material's macroscopic physic-mechanical parameters.
3.3 microstructural analysis based on micromineralogy
The sensitivity analyses of similar materials described above indicate that clay minerals significantly impact the macroscopic physical properties (including disintegration) and macroscopic mechanical properties of the materials. In order to reveal the mechanism of this influence, we carried out micromorphological observations based on micromineralogy and quantitative microstructural analyses based on Avizo software on four samples with different gypsum-kaolinite ratios (referred to as G: K).
3.3.1 qualitive microstructural analysis
Figure 11 shows the micro-morphological images of four different gypsum-kaolinite ratios for samples S1-S4 (see Fig. 12(a1) -(d1), with the yellow dashed box as the selected ROI region), the zoomed-in images of the corresponding ROI regions (see Fig. 12(a2) -(d2)) and the results of pore segmentation of the (a1) -(d1) images based on the dynamic threshold segmentation method of Avizo (see Fig. 12(a3) -(d3)).
Figures 12(a1), (a2), and (a3) show the typical crystal cementation microstructural morphology of the S1 sample and the corresponding pore segmentation results, respectively. From (a1) and (a2), it can be seen that the S1 sample is a typical crystalline cementation, with a large number of needle-like gypsum crystals interlaced with each other, which effectively bonded the quartz sand particles together to form a dense micromechanical structure. Its pore segmentation results show that when G: K = 4: 0, the material, the pores are mostly uniformly distributed isolated tiny pores; although the number of pores is relatively large, the pore size and the total area of the pores are small.
Figure 12(b1), (b2) shows the microstructural morphology of the S2 sample (G: K = 3:1). Due to the presence of a small amount of kaolinite minerals (near-pristine flake minerals adhering to the surface of the gypsum crystals, blue-circled area in Fig. 12(b2)), large isolated pores begin to grow in the vicinity of the kaolinite minerals, attributed to the water-loss and shrinkage properties of kaolinite. The micro-mechanical structure formed by gypsum crystals starts to deteriorate in S2 sample compared to S1. Its pore segmentation results show (see Fig. 12(b3)) that the overall uniformly distributed isolated tiny pores are still dominant within the S2 sample compared to S1, but the generation of larger diameter pores results in a moderate increase in the average pore size and total pore area of the S2 sample.
Figure 12 (c1) and (c2) show the microstructural morphology characteristics of the S3 sample (G: K = 1:3). The results show that a large number of cemented associations are distributed inside the sample, consisting of larger-sized cement formed by the absorption and fusion of the dominant kaolinite minerals and a small proportion of needle-like gypsum crystals. Compared to S2, larger-sized water-loss shrinkage pores begin to be formed and interconnected, resulting in a significant deterioration of the microstructure within the material. The pore segmentation results of the S3 sample showed (see Fig. 12(c3)) that larger-sized interconnected pores were distributed non-uniformly within the S3 sample compared to S2, resulting in a significant increase in the average pore size and total pore area within the material.
Figure 12(d1) -(d2) show the typical microstructural morphology of loose argillaceous cementation for S4 sample (G: K = 0:4). Compared with S3, not only the shrinkage pores inside the material further expanded and connected to form larger-size pores, but the large-size cracks began to generate and connected with the surrounding pores and formed a complex pore-crack network with higher connectivity, resulting in an extremely loose microstructure of the material. The pore segmentation results of S3 (see Fig. 12(d3)) showed that the internal pores of the material are mostly non-uniformly distributed and connected large-sized types, leading to a high average diameter and total area of pore cleavage (although the number of pores is significantly lower).
3.3.2 Quantitative microstructural analysis
Understanding the variation of pore size distribution with clay content in geotechnical materials is important for revealing the deterioration mechanism of clay minerals. The frequency distribution of pore number and area of the four samples in Fig. 12(a3) -(d3) were detailed quantified by using the Label Analysis module in Avizo, and the results are shown in Fig. 12. Overall, the total pore counts of the four samples from S1-S4 were 1640, 1407, 1218, and 323, respectively, which decreased significantly with the increase of kaolinite content. Meanwhile, the total pore areas of the S1-S4 samples were 2265.13, 2540.15, 3239.52, and 3827.69 µm2, respectively, exhibiting a noticeable increase. As can be seen from the pore size distribution (PSD), most of the pores of all the samples were distributed in the range of 0–3 µm, indicating that a large number of micropores were distributed within all four samples. However, the incorporation of kaolinite significantly changed the PSD characteristics. Specifically, the pore number, number frequency, and area frequency within 0–3 µm decreased from 1573, 95.91%, and 52.75% of S1 samples to 283, 87.68%, and 7.14% of S4. And the corresponding parameters of pores > 10 µm in equivalent diameter increased from 0 in S1 and S2 to 5, 0.41%, 17.25% in S3, and 8, 3.10%, 74.19% in S4. In addition, the number, number frequency, and area frequency of 3–10 µm pores in S3 (96,7.88%,66.57%) also showed a significant increase compared to S1 (66, 4.08%, 47.25%) and S2 (89,6.32%, 60.63%), indicating that tiny pores gradually expanded and merged into larger-sized ones, which can be regarded as a transitional pattern of pore structure in S4 samples.
To further investigate the change in other characteristics of pores with increasing kaolinite content, the Volume Fraction and Fractal Dimension modules in Avizo were used to calculate the porosity and fractal dimension of the specimens, respectively, as shown in Fig. 14. below. It indicated that the samples' equivalent diameter, porosity, and fractal dimension increased to different degrees with the increase of kaolinite content. That is, the equivalent diameter and porosity of the samples increased exponentially from 0.93 µm and 11.82% to 1.87 µm and 20.57%, respectively, while the fractal dimension increased linearly from 1.639 to 1.675. Analyzed in conjunction with Fig. 13(b)-(d) as well as Fig. 12, this can be attributed to the strong development of 10–15 µm pores in the S3 sample and > 10 µm pores in the S4 sample. These large pores are interconnected with other pores to form a more complex pore network with greater connectivity. In conclusion, the incorporation of kaolinite obviously changed the pore structure characteristics and strongly weakened the cementation strength of the material, transforming the material from a dense crystalline cementation to a porous clay cementation, showing an obvious deterioration effect.
In conclusion, with the increase of kaolinite content, the PSD of the pores within the samples, as well as the micro-mechanical parameters, such as porosity and fractal dimension, changed significantly, causing the rock microstructure to gradually change from a homogeneous and dense crystalline cement to an anisotropic and loose argillaceous cement. Combined with the conclusion of the sensitivity analysis in Section 3.2, we can conclude that the presence of clay minerals completely changes the mineralogical properties and micro-mechanical properties of the material, which in turn causes the degradation of the macro-mechanical properties of the material and the enhancement of water sensitivity.