Next Article in Journal
Energy Action Mechanism of Reinforced Sandstone under Triaxial Cyclic Loading and Unloading
Next Article in Special Issue
Deformation Mechanism and Control of In-Situ Assembling Caisson Technology in Soft Soil Area under Field Measurement and Numerical Simulation
Previous Article in Journal
Modulation of Structure and Optical Property of Nitrogen-Incorporated VO2 (M1) Thin Films by Polyvinyl Pyrrolidone
Previous Article in Special Issue
Experimental Investigation of Pre-Flawed Rocks under Dynamic Loading: Insights from Fracturing Characteristics and Energy Evolution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Model of Stress-Damage-Permeability Relationship of Weakly Cemented Rocks under Triaxial Compressive Conditions

1
School of Mines, China University of Mining & Technology, Xuzhou 221116, China
2
School of Resources and Earth Science, China University of Mining & Technology, Xuzhou 221116, China
3
School of Mines, Saint Petersburg Mining University, 199106 Saint Petersburg, Russia
*
Author to whom correspondence should be addressed.
Materials 2023, 16(1), 210; https://doi.org/10.3390/ma16010210
Submission received: 26 November 2022 / Revised: 23 December 2022 / Accepted: 23 December 2022 / Published: 26 December 2022

Abstract

:
To unravel the permeability variation mechanism of weakly cemented rocks (WCR), the paper conducted triaxial permeability tests on weakly cemented sandstones (WCS) collected from the Jurassic formation in northwest China. The paper identified the correlation of WCS permeability versus porosity, cementation structure, and mineral composition, further developing a model to characterize the WCS stress–damage–permeability relationship. The research indicated that the WCS permeability was initially high due to the naturally high porosity, large pore diameter, and loose particle cementation, thus favoring a significant decline as pore convergence in the compaction stage. In the residual stage, kaolinite and montmorillonite minerals disintegrated into water and narrowed fractures, causing a slight permeability increase from the initial to the maximum and residual stages. The WCS matrix fracturing was phenomenologically accompanied by clay mineral disintegration. By assuming that the matrix can be compressed, jointed, and fractured, the paper defined a damage variable D and accordingly developed a stress–damage–permeability relationship model that incorporated matrix compression, jointing, and fracturing. The model can describe the WCS permeability regime regarding the high initial permeability and slight difference of the maximum and residual permeabilities versus the initial.

1. Introduction

Weakly cemented rocks (WCR) widely exist in nature and contribute to the impact of geotechnical projects on the localized water system and ecological environment. Understanding WCR permeability responses to different stress conditions helps provide guidance for operating hydropower, underground tunneling, and mine projects [1,2,3,4]. In northwest China, the rock units in Jurassic and Cretaceous formations are often weakly cemented, generally featuring young sedimentation, low strength, and loose cementation. Such properties cause weakly cemented rocks to be sensitive to engineering perturbations [5,6]. Moreover, the arid and semi-arid climate, water resource scarcity, and eco-environmental fragility in northwest China further aggravate the impact of WCR permeability increase on groundwater depressurization and thus environmental deterioration [7,8]. Therefore, understanding the WCR permeability variation mechanism is of great significance for both geotechnical practice and eco-environmental protection.
The current progress of weakly cemented rocks is mainly in rock properties and permeability [9,10,11]. In rock mechanics, Sharma and Wang et al. conducted laboratorial tests and quantified the correlation of bulk density and P-wave velocity versus axial compressive strength and elastic modulus of weakly cemented mudstone (WCM) and sandstone (WCS) [12,13]. Konstantinou et al. analyzed the relationship of WCS tensile strength with fracture patterns, surface features, and cementing degree [14]. Li and Zhao et al. studied the stress–strain, acoustic emission, and energy evolution characteristics of weakly cemented rocks under Brazilian test conditions, by which their relationship was discussed [15,16]. Liu et al. conducted triaxial compression tests and analyzed the impact of water content and confining pressure on fracture pattern and stress–strain curves of WCM and WCS samples [17]. Tommasi and Ye et al. proposed that WCR mechanical behavior is correlated to physical and structural properties, further analyzing the impact of pore abundance, mineral composition, and matrix-pore continuity on rock strength [18,19]. Regarding WCR hydraulics, Nguyen and Ruistuen et al. analyzed the impact of the stress path on WCS permeability and identified that the permeability had a significant decline when the shear strength was exceeded, which was attributed to particle friction and pore tortuosity increase [20,21]. You et al. and Liu et al. analyzed the impact of temperature, confining pressure, and water content on WCR permeability and developed a fractional quaternion model to describe WCR creep behavior in three continuous stages [22,23]. Liu et al. analyzed WCS permeability variation with porosity and mercury injection rate, by which the pore–throat ratio was considered a major indicator of permeability mutation [24]. Lyu et al. conducted laboratorial tests and analyzed WCS permeability features with varying particle diameters, revealing that the permeability coefficient was positively correlated to particle diameter in the same confining pressure and hydraulic pressure condition, between which the pore water pressure had greater impacts [25].
Previous research has investigated WCR strength properties and tried qualitative analysis on the permeability variation from physical perspectives, consolidating the basis for this paper. According to knowledge from geotechnical practices, the microstructural damage induced by stress environment change should be the critical contribution to WCR permeability variation. In this context, the paper collected WCR samples from the Jurassic formation in Ili mine fields, China and conducted complete stress–strain permeability tests using Rock Top multi-physics platform. The paper further studied the impact of porosity feature, cementation structure, and mineral composition on permeability variation. A WCR representative element volume (REV) model was constructed by considering physical mechanics and permeability property, by which, coupled with defining a damage variable D, a WCR stress–damage–permeability relationship model was developed. This research can provide theoretical reference for geotechnical operation and eco-environmental protection in a weakly cemented rock environment.

2. Materials

2.1. Sampling

The research first collected WCS samples from the Middle Jurassic Xishanyao Formation in Ili No. 4 Mine in Ili, Xinjiang Province, China. The downhole sampling adopted pump-suction reverse circulation drilling technology to maintain sample intactness to the maximum extent. The sampling position was 138.55 to 147.34 m below the topography, as shown in Figure 1.

2.2. Rock Physics

As indicated by some studies, rock mass permeability is closely correlated to the porosity, cementing form, and sediment percentage [26,27,28,29]. Before laboratorial permeability tests, the paper conducted mercury intrusion tests, scanning electron microscope (SEM, Tescan, Brno, Southern Moravia, Czech Republic) imaging, and X-ray diffraction (XRD, Bruker, Billerica, MA, USA) measurements to obtain WCS porosity, cementation structure, and mineral composition data.
Mercury intrusion tests measured WCR porosity of 14.78%. The mercury injection rate of WCR was fast, indicating that the proportion of large-diameter pores in rocks was high. At the same time, the mercury ejection rate of WCR was slow, indicating that the pore connectivity of WCR was poor and the pores were mainly closed or partially closed. Rock pores characterize medium diameters: 1.8–10 μm diameter dominated, taking 33.53% of the total amount, followed by 10–50 μm diameter by 30.43%. WCR was mainly composed of intergranular middle pores [30].The mercury injection–ejection curves and porosity distribution are shown in Figure 2a,b, respectively. Figure 2c shows SEM imaging results, indicating that the pores were extensive and partially filled with fragments, compared to which the mineral particles had greater diameters about 25–60 μm and loosely distributed with explicit boundaries. The XRD measurement obtained WCS mineral compositions, as shown in Figure 2d, revealing that quartz, kaolinite, and montmorillonite took 90% of the total mineral content, among which the quartz and kaolinite minerals were about 34.7% and 33.2%, almost equivalent. As indicated by the XRD results, WCS samples were rich in clay minerals including kaolinite, montmorillonite, and illite, with a percentage of about 55.9%.

3. Complete Stress–Strain Permeability Tests

3.1. Methods

Permeability tests were conducted using the Rock Top multi-physics platform. The platform included a controller system for parameter and servo-control configuration, a hydraulic loading system for executing loading and unloading operations along axial and radius directions, and a monitoring system for sample encapsulation and data (displacement, pressure, and flux) monitoring. The maximum axial stress, maximum confining pressure, and maximum water pressure of the loading system were 500 MPa, 60 MPa, and 50 MPa, respectively. In order to accurately measure the axial deformation and radial strain, two Linear Variable Displacement Transducer (LVDT) displacement sensors were placed in parallel on both sides of the sample with a range of 12 mm and an accuracy of 0.001 mm [31]. Figure 3 exhibits the platform and its three components. The virgin cores were standardized into cylinders with a 100 mm height and 50 mm diameter; the parallelism tolerance of both ends was controlled within 0.1 mm and the diameter difference of both ends were within 0.2 mm. The tests adopted the displacement control method. The samples were saturated using pure water and then tested by keeping the hydraulic pressure difference and loading rate constant. Parameters of the permeability tests as shown in Table 1.

3.2. Results

Figure 4 and Figure 5 show WCS stress–strain–permeability curves and sample failure scenarios, respectively.
Figure 4 shows that WCS permeability fluctuated throughout the complete stress–strain stages. The permeability first declined due to rock matrix compression, followed by a further but slight decline to the minimum in the elastic stage. Since rock yielding, the permeability increased and reached the maximum in the middle failure stage. There was again a permeability decline in the later failure stage, after which the permeability was almost unchanged in the residual stage. Of note is that the peak permeability was behind the peak deviatoric stress. An explicit feature of weakly cemented rocks was that the initial permeability (k0) was relatively high and had a slight difference against the maximum (kmax) and residual permeability (kres); for example, under the 5 MPa confining pressure condition, the ratio kmax/k0 and kres/k0 was only 3.34 and 2.85, respectively. Table 2 lists the three permeability thresholds.
From Figure 5, it can be seen that the WCS sample failure was mainly characterized by shear fracturing; along fractures there was explicit disintegration of clay minerals.

3.3. Relationship between WCR Permeability and Physical Properties

The permeability variation indicated that WCR permeability was closely correlated to the porosity, cementation structure, and mineral composition. The high initial permeability was directly dependent on the high virgin porosity, large pore diameter, and loose connection of mineral particles. The permeability significantly decreased in the compaction stage, practically resulting from the loading-induced compression on extensive pores and thus flow pathways. In the failure stage, the maximum permeability was slightly different than the initial permeability since the abundant clay minerals such as kaolinite and montmorillonite expanded due to encountering water and partially narrowed the fluid flow pathways. As a result, the maximum permeability kept a low level. The residual stage experienced fracture propagation and closure and clay mineral disintegration, thus favoring an almost constant permeability due to a dynamic balance between fresh and existing fractures. In addition, clay minerals took up to 50%, being a major part of the WCR matrix. The clay mineral disintegration suggested that the WCR matrix fractured under the triaxial compression condition.

4. Stress–Damage–Permeability Relationship Model

The above analysis indicated that the high porosity, weak cementation structure, and abundant clay minerals caused WCRs to characterize distinctive permeability regimes. Considering that permeability variation is a reflection of rock material damage induced by pressurization, the paper introduced a damage variable D to help quantify WCS permeability variation. On this basis, a stress–damage–permeability relationship model was then developed.

4.1. Model Construction

Based on WCS physics and permeability features identified above and existing studies [32], the paper placed the following hypotheses [33,34,35,36].
(1)
Weakly cemented rocks were composed of matrix and fractures; the matrix contained natural pores that can be squeezed and closed in the compaction stage and allow matrix materials to be jointed.
(2)
The matrix can fracture in response to triaxial compression, further affecting the effective aperture of fractures and thus permeabilities.
The weakly cemented rocks were described using representative element volume (REV) units as shown in Figure 6.
The REV was the smallest micro element structure that represents the deformation of loaded coal rock mass. Accordingly, the porosity of REV units can be expressed via fracture aperture and matrix geometry [37]:
ϕ f ϕ f 0 = a 0 b 0 b a = ( 1 + Δ b b 0 ) a 0 a
where ϕ f and ϕ f 0 are respectively the current and initial porosity of fractures, a 0 and b 0 are the initial geometry of matrix and pores, a and b are the current geometry of matrix and pores, and Δ b represents the variation of pore geometry.
Existing studies on WCR porosity variation generally placed less focus on the impact of matrix deformation on permeability. In this research, the laboratorial tests revealed that weakly cemented rocks were of high porosity and large pore diameter; the large pores can be geometrically compressed in the compaction stage, and WCR samples well degraded in water. By implication, the paper assumed that the fractured matrix can joint in compressive and fluid flow conditions.
According to this hypothesis, the WCR pore geometry variation Δ b was considered composed of matrix geometry change b M and fracture aperture change b F , with proportion of γ M and γ F respectively. Therefore,
{ Δ b = b M + b F γ M = b M / Δ b γ F = b F / Δ b
where b M and b F are respectively the contribution of matrix and fracture to pore geometry variation, and γ M and γ F are ratios of the respective contribution.
The volumes of fracture and matrix are:
{ V F = a 2 b F ; V F 0 = a 0 2 b F 0 V M = a 2 b M ; V F 0 = a 0 2 b M 0
where V F and V M are volumes of fracture and matrix, with initial values expressed as V F 0 and V M 0 , respectively.
The fracture volume and matrix volume have another expression in natural strain and engineering strain [38]:
{ V F = V F 0 exp ( Δ σ e K F ) V M = V M 0 ( 1 Δ σ e K M )
where Δ σ e is effective stress increment, K F is fracture bulk modulus, and K M is matrix bulk modulus. Considering that the matrix was difficult to compress, this work focused on the influence of K F on the WCR permeability. Meanwhile, the determination of K F was mainly obtained by fitting the stress–strain and permeability curves of WCR.
Combining Equations (2)–(4) obtained the pore geometry variation Δ b :
Δ b = ( a 0 a ) 2 [ b F 0 exp ( Δ σ e K F ) + b M 0 ( 1 Δ σ e K M ) ]
The matrix itself was rather rigid but can deform with the contribution of matrix jointing. Consider that the matrix geometry variation obeyed a function of strain ε e , and
b M = b M 0 a 0 a ( 1 f ( ε e ) ) 2
where ε e is matrix strain and f ( ε e ) the matrix geometry change, following f ( ε e ) = exp ( ε e ) .
Substituting Equation (6) into Equation (5) obtained:
Δ b = ( a 0 a ) 2 [ b F 0 exp ( Δ σ e K F ) + b M 0 ( 1 f ( ε e ) ) 2 ]
Incorporating Equations (1), (2), and (7), with cubic low, obtained WCR permeability expression:
k = k 0 ( a 0 a ) 3 { 1 + ( a 0 a ) 2 [ γ F exp ( Δ σ e K F ) + ( 1 γ F ) ( 1 f ( ε e ) ) 2 ] } 3
where k 0 and k are respectively the initial and current permeability.
With continuous pressurization, WCR units damaged and then fractured as the damage gradually accumulated to a certain extent. The fracture initiation and development were accompanied by the matrix geometry variation. By assuming the number of natural fractures and fractures corresponding to the WCR unit damage was respectively N 0 and N , the matrix geometry variation can be written as:
a 0 a = N N 0 = N 0 + Δ N N 0 = 1 + Δ N N 0 = 1 + D ς
where Δ N refers to the fracture number increment, D is damage variable, and ς is the coefficient of matrix geometry variation.
Substituting Equation (9) into Equation (8) obtained a stress–damage–permeability relationship:
k = k 0 ( 1 + D ς ) 3 _ { 1 + ( 1 + D ς ) 2 _ [ γ F exp ( Δ σ e K F ) _ + ( 1 γ F ) ( 1 f ( ε e ) ) 2 _ ] } 3                     
where items ①, ②, and ③ represent the impact of fracture number, aperture convergence, and matrix joint on permeability, respectively.

4.2. Damage Variable Determination

Assume the damage variable D obeying Weibull distribution [39]:
P ( F * ) = m 0 F 0 ( F * F 0 ) m 0 1 exp [ ( F * F 0 ) m 0 ]
where F * represents microstructure strength, such as deformation modulus, tensile strength, and cohesion, m0 is coefficient of average degree, and F0 is mean value of the strength parameter F * .
By Drucker–Prager failure criterion, the microstructure strength F * can be expressed as:
{ F * = f ( σ * ) = α 0 I 1 + J 2 1 / 2 α 0 = sin φ 9 + 3 sin 2 φ I 1 = σ 1 + σ 2 + σ 3 J 2 = 1 6 [ ( σ 1 σ 2 ) 2 + ( σ 2 σ 3 ) 2 + ( σ 1 σ 3 ) 2 ]
where φ is angle of friction, and σ 1 , σ 2 , and σ 3 refer to three effective principal stress. In this test, the second and third effective principal stresses ( σ 2 and σ 3 ) were consisted equivalent. I 1 and J 2 are respectively the first and second invariants of effective stress tensor.
Considering that the damage degree varied with positions, it was obtained via elastic damage theory [33,40]:
{ σ = E ε ( 1 C n D ) D = 1 exp [ ( F * F 0 ) m 0 ]
where E is deformation modulus, and C n is proportionality coefficient of the damage, reflecting the residual strength of WCR material.
To better reflect WCR mechanics and hydraulics regarding high initial permeability, rapid decline in the compaction stage, and slight change in the residual stage, the overall permeability variation was segmented into compaction elasticity, yielding damage, failure, and residual stages, corresponded to the complete stress–strain profile. Combining Equations (11)–(13) obtained a WCR statistical damage model:
σ 1 = { E ε 1 + 2 μ σ 3 ε 1 ε ela E ( ε 1 ε ela ) ( 1 C n + C n exp ( E ( ε 1 ε ela ) 3 ( σ 1 2 μ σ 3 ) ( sin φ ( σ 1 σ ela + 2 σ 3 ) 3 + sin 2 φ + σ 1 σ 3 ) + σ 1 σ ela σ 3 F 0 ) m 0 ) + σ ela + 2 μ σ 3    ε ela ε 1 ε max E ( ε max ε 1 ) 2 + 2 μ σ 3    ε max < ε 1 ε res   σ res       ε 1 > ε res
where ε ela is the strain by the end of the elastic stage, ε max is the peak strain, ε res is the strain at the start of the residual stage, and σ ela measures the axial stress when axial strain is ε ela .
Fitting the stress–strain data of each damage stage and residual coefficient determined from residual strength obtained m0, F0, ζ, and further the damage variable D of the WCR stress–strain curves under different confining pressures. Table 3 lists the curve fitting results.

4.3. Model Verification

To verify the WCR stress–damage–permeability relationship model, the statistical damage parameters summarized in Table 3 were further used to determine the damage variable D. Then, the WCR permeability variation under different confinement conditions can be calculated by substituting D into Equation (10). The comparison between calculation results and measurement is shown in Figure 7.
Figure 7 suggests a good agreement between theoretical and laboratorial test results, with an error of 12% in deviatoric stress and 15% in permeability, Table 4. This result indicates that the stress–damage–permeability relationship model (Equation (10)) can well describe WCR permeability regime in complete stress–strain stages.

5. Discussion

The WCRs in Jurassic and Cretaceous generally characterize young sedimentation, low strength, and loose cementation. In this work, the Jurassic weakly cemented sandstone is selected from Ili No. 4 Coal Mine in Xinjiang Province, China. The sample selection is very typical and representative. Therefore, the authors believe that the results of permeability test and the established stress–damage–seepage relationship model of WCR based on the typical samples can express its general permeability characteristics. Of course, in the future, the WCRs in more different areas should be tested and fitted to further verify the validity and universal applicability of the established stress–damage–seepage relationship model.
It can be seen from Figure 7c that when the confining pressure is 5.0 MPa, the fitting error of the deviatoric stress and permeability of weakly cemented sandstone is more than 10%. This may be due to the fact that when the confining pressure is 5.0 MPa, the weakly cemented sandstone has obvious elastic brittle failure characteristics (as shown in Figure 7c), and its stress–strain and permeability curves have mutation points (or distortion points). The existence of mutation points destroys the continuity of the curves, which makes its fitting error larger. In future studies, more weakly cemented samples under high confining pressure should be tested and fitted to further explain the above phenomenon.

6. Conclusions

Weakly cemented rocks showed permeability fluctuation in the complete stress–strain stages, including a decrease in compaction and elastic segments, an increase in yield and failure segments, and an again decrease to a leveling-off value in the ultimate residual segment. Such a permeability regime was physically correlated to WCR porosity, cementation structure, and mineral composition properties. WCR material had a high initial permeability due to the naturally high porosity, large pore diameter, and loose cementation structure, thus favoring a rapid permeability decline due to pore compression. Since the yield stage, WCR clay minerals such as kaolinite and montmorillonite decomposed into water and narrowed the fractures, causing that the maximum and residual permeabilities were slightly greater than the initial. This minor difference indicated that WCR materials had a limited permeability responses in triaxial confinement conditions, which helps WCR strata maintain a good water-resisting functionality.
XRD measurements revealed that WCR materials were rich in clay minerals. Kaolinite, montmorillonite, and illite minerals took more than 50%, constituting a major part of the WCR matrix. Evident clay mineral disintegration was observed with sample fracturing, implying matrix fracturing in response to triaxial compression.
By assuming that WCR matrix can be compressed, jointed, and fractured, the paper defined a damage variable D and further developed a stress–damage–permeability relationship model. Then, a D determination method was established according to Weibull distribution and the Drucker–Prager principle. The relationship model was verified with errors in deviatoric stress and permeability within 15% with respect to laboratorial measurements, suggesting that the model can represent the permeability regime of weakly cemented rocks in the complete stress–strain stages.

Author Contributions

Conceptualization, D.Z. and W.L.; data curation, S.L. and Z.F.; investigation, T.L. and S.L.; methodology, G.F.; validation, S.Z., T.L. and Z.F.; writing—original draft; S.Z.; writing—review and editing; G.F., D.Z. and W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 52204161; 51974291), the Fundamental Research Funds for the Central Universities (Grant No. 2022QN1008; 2021ZDPY0226), the Jiangsu Funding Program for Excellent Postdoctoral Talent (Grant No. 2022ZB511), the Shanxi Province Unveils Bidding Project (Grant No. 20201101009).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to the anonymous reviewers for their helpful comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Grgic, D.; Al Sahyouni, F.; Golfier, F.; Moumni, M.; Schoumacker, L. Evolution of gas permeability of rock salt under different loading conditions and implications on the underground hydrogen storage in salt caverns. Rock Mech. Rock Eng. 2021, 55, 691–714. [Google Scholar] [CrossRef]
  2. Rutqvist, J.; Borgesson, L.; Chijimatsu, M.; Hernelind, J.; Jing, L.; Kobayashi, A.; Nguyen, S. Modeling of damage, permeability changes and pressure responses during excavation of the TSX tunnel in granitic rock at URL, Canada. Environ. Geol. 2009, 57, 1263–1274. [Google Scholar] [CrossRef] [Green Version]
  3. Akbarimehr, D.; Aflaki, E. Site investigation and use of artificial neural networks to predict rock permeability at the Siazakh Dam, Iran. Q. J. Eng. Geol. Hydrogeol. 2019, 52, 230–239. [Google Scholar] [CrossRef]
  4. Khanal, M.; Guo, H.; Adhikary, D. 3D Numerical Study of Underground Coal Mining Induced Strata Deformation and Subsequent Permeability Change. Geotech. Geol. Eng. 2019, 37, 235–249. [Google Scholar] [CrossRef]
  5. Zhao, Z.; Lv, X.; Wang, W.; Tan, Y. Damage evolution of bi-body model composed of weakly cemented soft rock and coal considering different interface effect. Springerplus 2016, 5, 292. [Google Scholar] [CrossRef] [Green Version]
  6. Zhang, S.; Fan, G.; Zhang, D.; Li, Q. Physical simulation research on evolution laws of clay aquifuge stability during slice mining. Environ. Earth Sci. 2018, 77, 278. [Google Scholar] [CrossRef]
  7. Fan, G.; Chen, M.; Zhang, D.; Wang, Z.; Zhang, S.; Zhang, C.; Li, Q.; Cao, B. Experimental study on the permeability of weakly cemented rock under different stress states in triaxial compression tests. Geofluids 2018, 2018, 9035654. [Google Scholar] [CrossRef]
  8. Zhang, S.; Fan, G.; Zhang, D.; Li, S.; Fan, Y.; Luo, T. Impacts of longwall mining speeds on permeability of weakly cemented strata and subsurface watertable: A case study. Geomat. Nat. Hazards Risk 2021, 12, 3063–3088. [Google Scholar] [CrossRef]
  9. Dong, J.; Tu, C.; Lee, W.; Jheng, Y. Effects of hydraulic conductivity/strength anisotropy on the stability of stratified, poorly cemented rock slopes. Comput. Geotech. 2012, 40, 147–159. [Google Scholar] [CrossRef]
  10. Holtzman, R. Micromechanical model of weakly-cemented sediments. Int. J. Numer. Anal. Methods Geomech. 2012, 36, 944–958. [Google Scholar] [CrossRef]
  11. Yeh, P.; Lee, K.Z.; Chang, K. 3D Effects of permeability and strength anisotropy on the stability of weakly cemented rock slopes subjected to rainfall infiltration. Eng. Geol. 2020, 266, 105459. [Google Scholar] [CrossRef]
  12. Wang, Z.; Li, W.; Wang, Q.; Liu, S.; Hu, Y.; Fan, K. Relationships between the petrographic, physical and mechanical characteristics of sedimentary rocks in Jurassic weakly cemented strata. Environ. Earth Sci. 2019, 78, 131. [Google Scholar] [CrossRef]
  13. Sharma, M.S.R.; O’Regan, M.; Baxter, C.D.P.; Moran, K.; Vaziri, H.; Narayanasamy, R. Empirical relationship between strength and geophysical properties for weakly cemented formations. J. Pet. Sci. Eng. 2010, 72, 134–142. [Google Scholar] [CrossRef]
  14. Konstantinou, C.; Biscontin, G.; Logothetis, F. Tensile strength of artificially cemented sandstone generated via microbially induced carbonate precipitation. Materials 2021, 14, 4735. [Google Scholar] [CrossRef]
  15. Zhao, Y.; Liu, B. Deformation field and acoustic emission characteristics of weakly cemented rock under brazilian splitting test. Nat. Resour. Res. 2021, 30, 1925–1939. [Google Scholar] [CrossRef]
  16. Li, P.; Li, Z.; Liu, B.; Teng, T.; Guo, J. Experimental investigation on the tensile properties of weakly cemented sandstone in China Shendong mining area. Therm. Sci. 2020, 24, 3987–3994. [Google Scholar] [CrossRef]
  17. Liu, H.; Zhang, D.; Zhao, H.; Chi, M.; Yu, W. Behavior of weakly cemented rock with different moisture contents under various tri-axial loading stats. Energies 2019, 12, 1563. [Google Scholar] [CrossRef] [Green Version]
  18. Yu, W.; Li, K.; Liu, Z.; An, B.; Wang, P.; Wu, H. Mechanical characteristics and deformation control of surrounding rock in weakly cemented siltstone. Environ. Earth Sci. 2021, 80, 337. [Google Scholar] [CrossRef]
  19. Tommasi, P.; Verrucci, L.; Rotonda, T. Mechanical properties of a weak pyroclastic rock and their relationship with microstructure. Can. Geotech. J. 2015, 52, 211–223. [Google Scholar] [CrossRef]
  20. Ruistuen, H.; Teufel, L.W.; Rhett, D. Influence of reservoir stress path on deformation and permeability of weakly cemented sandstone reservoirs. SPE Reserv. Eval. Eng. 1999, 2, 266–272. [Google Scholar] [CrossRef]
  21. Nguyen, V.H.; Gland, N.; Dautriat, J.; David, C.; Wassermann, J.; Guélard, J. Compaction, permeability evolution and stress path effects in unconsolidated sand and weakly consolidated sandstone. Int. J. Rock Mech. Min. Sci. 2014, 67, 226–239. [Google Scholar] [CrossRef]
  22. You, S.; Ji, H.; Wang, T.; Song, Z. Thermal and mechanical coupling effects on permeability of weakly cemented sandstone. Emerg. Mater. Res. 2018, 7, 100–108. [Google Scholar] [CrossRef]
  23. Liu, J.; Jiang, H.; Meng, B.; Wang, L.; Yang, J.; Zhang, X. A four-element fractional creep model of weakly cemented soft rock. Bull. Eng. Geol. Environ. 2020, 79, 5569–5584. [Google Scholar] [CrossRef]
  24. Liu, Q.; Sun, Y.; Li, J. Experimental study on seepage characteristics of Jurassic weakly cemented sandstone under water-rock interaction. Geofluids 2020, 2020, 854368. [Google Scholar] [CrossRef]
  25. Lyu, X.; Zhao, Z.; Wang, X.; Wang, W. Study on the Permeability of Weakly Cemented Sandstones. Geofluids 2019, 2019, 8310128. [Google Scholar] [CrossRef] [Green Version]
  26. Han, T.; Best, A.I.; Sothcott, J.; North, L.J.; Macgregor, L.M. Relationships among low frequency (2 Hz) electrical resistivity, porosity, clay content and permeability in reservoir sandstones. J. Appl. Geophys. 2015, 112, 279–289. [Google Scholar] [CrossRef]
  27. Chao, Z.; Ma, G.; He, K.; Wang, M. Investigating low-permeability sandstone based on physical experiments and predictive modeling. Undergr. Space 2021, 6, 364–378. [Google Scholar] [CrossRef]
  28. Singh, V.K.; Kumar, D.; Kashyap, P.S.; Singh, P.K.; Kumar, A.; Singh, S.K. Modelling of soil permeability using different data driven algorithms based on physical properties of soil. J. Hydrol. 2020, 580, 124223. [Google Scholar] [CrossRef]
  29. Liu, X.; Xu, M.; Wang, K. Mechanism of permeability evolution for reservoir sandstone with different physical properties. Geofluids 2018, 2018, 5327895. [Google Scholar] [CrossRef] [Green Version]
  30. Saki, M.; Siahpoush, S.; Khaz’ali, A.R. A new generalized equation for estimation of sandstone and carbonate permeability from mercury intrusion porosimetry data. J. Pet. Explor. Prod. Technol. 2020, 10, 2637–2644. [Google Scholar] [CrossRef]
  31. Chen, L.; Zhang, D.; Zhang, W.; Guo, J.; Yao, N.; Fan, G.; Zhang, S.; Wang, X.; Wu, P. Experimental Investigation on Post-Peak Permeability Evolution Law of Saturated Sandstone under Various Cyclic Loading-Unloading and Confining Pressure. Water 2022, 14, 1773. [Google Scholar] [CrossRef]
  32. Greco, O.D.; Ferrero, A.M.; Oggeri, C. Experimental and analytical interpretation of the behaviour of laboratory tests on composite specimens. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1993, 30, 1539–1543. [Google Scholar] [CrossRef]
  33. Fan, Z.; Fan, G.; Zhang, D. Permeability and energy evolution characteristics of heterogeneous coal and rock mass. Nat. Resour. Res. 2021, 30, 4493–4514. [Google Scholar] [CrossRef]
  34. Li, L.; Tu, W.; Zhou, Z.; Shi, S.; Zhang, M.; Chen, Y. Dynamic unloading instability mechanism of underground cavern based on seepage-damage coupling. KSCE J. Civ. Eng. 2020, 24, 1620–1631. [Google Scholar] [CrossRef]
  35. Jiang, A.; Yang, X.; Xu, M.; Jiang, T. Coupled hydrologic-mechanical-damage analysis and its application to diversion tunnels of hydropower station. Adv. Civ. Eng. 2021, 2021, 8341528. [Google Scholar] [CrossRef]
  36. Wang, S.; Yang, T.; Zhang, Z.; Sun, Z. Unsaturated seepage-stress-damage coupling and dynamic analysis of stability on discrete fractured rock slope. Environ. Earth Sci. 2021, 80, 660. [Google Scholar] [CrossRef]
  37. Fan, Z.; Fan, G.; Zhang, D.; Zhang, L.; Zhang, S.; Liang, S.; Yu, W. Optimal injection timing and gas mixture proportion for enhancing coalbed methane recovery. Energy 2021, 222, 119880. [Google Scholar] [CrossRef]
  38. Qiao, T.; Zhang, G.; Hu, S.; Zhao, S. Research on constitutive model of microcrack compaction section in brittle rock. J. Saf. Sci. Technol. 2017, 13, 128–135. [Google Scholar] [CrossRef]
  39. Pan, X.; Guo, W.; Wu, S.; Chu, J. An experimental approach for determination of the Weibull homogeneity index of rock or rock-like materials. Acta Geotech. 2020, 15, 375–391. [Google Scholar] [CrossRef]
  40. Xiao, W.; Zhang, D.; Wang, X.; Yang, H.; Wang, X.; Wang, C. Research on Microscopic Fracture Morphology and Damage Constitutive Model of Red Sandstone Under Seepage Pressure. Nat. Resour. Res. 2020, 29, 3335–3350. [Google Scholar] [CrossRef]
Figure 1. Borehole profile and WCS sampling position.
Figure 1. Borehole profile and WCS sampling position.
Materials 16 00210 g001
Figure 2. Physical properties of weakly cemented rocks: (a) mercury injection–ejection curves; (b) pore diameter distribution; (c) microscopic structures of pores and minerals via SEM imaging; (d) mineral composition via XRD measurement.
Figure 2. Physical properties of weakly cemented rocks: (a) mercury injection–ejection curves; (b) pore diameter distribution; (c) microscopic structures of pores and minerals via SEM imaging; (d) mineral composition via XRD measurement.
Materials 16 00210 g002
Figure 3. Rock Top multi-physics platform: (a) physical picture; (b) schematic diagram.
Figure 3. Rock Top multi-physics platform: (a) physical picture; (b) schematic diagram.
Materials 16 00210 g003
Figure 4. WCS stress–strain–permeability curves under different confining pressures: (a) confining pressure 2.0 MPa; (b) confining pressure 3.5 MPa; (c) confining pressure 5.0 MPa.
Figure 4. WCS stress–strain–permeability curves under different confining pressures: (a) confining pressure 2.0 MPa; (b) confining pressure 3.5 MPa; (c) confining pressure 5.0 MPa.
Materials 16 00210 g004
Figure 5. WCS sample failure.
Figure 5. WCS sample failure.
Materials 16 00210 g005
Figure 6. REV description of weakly cemented rocks.
Figure 6. REV description of weakly cemented rocks.
Materials 16 00210 g006
Figure 7. WCR permeability variation curves under different confinement conditions: (a) confining pressure 2.0 MPa; (b) confining pressure 3.5 MPa; (c) confining pressure 5.0 MPa.
Figure 7. WCR permeability variation curves under different confinement conditions: (a) confining pressure 2.0 MPa; (b) confining pressure 3.5 MPa; (c) confining pressure 5.0 MPa.
Materials 16 00210 g007
Table 1. Parameters of the permeability tests.
Table 1. Parameters of the permeability tests.
ProjectsConfining Pressure (MPa)Osmotic Pressure by Both Ends (MPa)Hydraulic Difference (MPa)Loading Controlled byLoading Rate (mm/min)
12.0Top: 0.5
Bottom: 1.5
1.0Displacement0.03
23.5
35.0
Table 2. WCS permeability thresholds (initial, maximum, and residual) in the complete stress–strain stages.
Table 2. WCS permeability thresholds (initial, maximum, and residual) in the complete stress–strain stages.
LithologyConfining Pressure (MPa)k0 (m2)kmax (m2)kres (m2)Ratio
kmax/k0kres/k0kmax/kres
Sandstone2.01.59 × 10−172.91 × 10−172.61 × 10−171.831.641.11
3.52.91 × 10−186.50 × 10−186.06 × 10−182.232.081.07
5.08.75 × 10−192.92 × 10−182.49 × 10−183.342.851.17
Table 3. WCR damage parameters.
Table 3. WCR damage parameters.
Lithology Confining Pressure (MPa)Statistical Parameters
m0F0ζ
Sandstone2.03.597.021.74
3.53.749.091.52
5.01.7414.471.02
Table 4. Error of theoretical WCR deviatoric stress and permeability versus laboratorial measurements.
Table 4. Error of theoretical WCR deviatoric stress and permeability versus laboratorial measurements.
LithologyConfining Pressure (MPa)Mean Relative Error (%)
Deviatoric StressPermeability
Sandstone2.07.247.13
3.55.596.32
5.011.3914.79
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, S.; Fan, G.; Zhang, D.; Li, W.; Luo, T.; Liang, S.; Fan, Z. A Model of Stress-Damage-Permeability Relationship of Weakly Cemented Rocks under Triaxial Compressive Conditions. Materials 2023, 16, 210. https://doi.org/10.3390/ma16010210

AMA Style

Zhang S, Fan G, Zhang D, Li W, Luo T, Liang S, Fan Z. A Model of Stress-Damage-Permeability Relationship of Weakly Cemented Rocks under Triaxial Compressive Conditions. Materials. 2023; 16(1):210. https://doi.org/10.3390/ma16010210

Chicago/Turabian Style

Zhang, Shizhong, Gangwei Fan, Dongsheng Zhang, Wenping Li, Tao Luo, Shuaishuai Liang, and Zhanglei Fan. 2023. "A Model of Stress-Damage-Permeability Relationship of Weakly Cemented Rocks under Triaxial Compressive Conditions" Materials 16, no. 1: 210. https://doi.org/10.3390/ma16010210

APA Style

Zhang, S., Fan, G., Zhang, D., Li, W., Luo, T., Liang, S., & Fan, Z. (2023). A Model of Stress-Damage-Permeability Relationship of Weakly Cemented Rocks under Triaxial Compressive Conditions. Materials, 16(1), 210. https://doi.org/10.3390/ma16010210

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop