Research Article  Open Access
Yong He, Keneng Zhang, Dongyu Wu, "Experimental and Modeling Study of Soil Water Retention Curves of Compacted Bentonite considering Salt Solution Effects", Geofluids, vol. 2019, Article ID 4508603, 11 pages, 2019. https://doi.org/10.1155/2019/4508603
Experimental and Modeling Study of Soil Water Retention Curves of Compacted Bentonite considering Salt Solution Effects
Abstract
Due to its low hydraulic conductivity, high swelling capacity, and good adsorption properties, the Gaomiaozi (GMZ) bentonite has been selected as potential buffer/backfill materials for construction of engineered barriers in the deep geological repository for disposal of highlevel nuclear waste (HLW) in China. Investigation of salt solution effects on the water retention properties of compacted bentonite is of great importance in the context of geological disposal of HLW based on the multibarrier concept. In this study, amended specimens were obtained through a spray of different concentrations of salt solutions to reach target salt contents, respectively. With employment of the vapor phase technique for suction control, water retention tests were conducted on densely compacted GMZ bentonite (1.7 Mg/m^{3}) with different salt contents under confined conditions. Corresponding soil water retention curves (SWRCs) were obtained. Analysis indicates that, for a given suction, the measured water content of GMZ bentonite specimen increases as the salt content (or pore fluid concentration) increases. The influencing rate depends on suction. For lower suctions (lower than 38 MPa), the water retention capacity increases as the salt content increases, while for higher suctions (higher than 38 MPa), the influence can be negligible. Based on the Fredlund and Xing (1994) equation, a soil water retention model was proposed for simulation of the SWRCs of compacted GMZ bentonite with consideration of pore fluid chemistry. Parameters were analyzed and determined with consideration of influences of the pore fluid concentration. Verification indicates that the SWRCs simulated by the proposed model are well agreed with the measured ones.
1. Introduction
Compacted bentonite has been recommended as a buffer/backfill material in many countries for engineered barriers in the deep geological repository for disposing highlevel nuclear waste (HLW) [1–9]. The compacted bentonite will undergo a complex—highly coupled thermal (T), hydro (H), mechanical (M), and chemical (C)—process during the longterm operation of a repository [10–12], i.e., simultaneous heating from decay of nuclear waste in the canister, infiltration of groundwater from the surrounding geological formations, possible leakage of nuclides from the canister, and swelling pressure of bentonite on hydration. To assess these coupled processes of the engineered barrier system, a number of investigations have been conducted [1, 3, 9, 13–15] and related modeling approaches [16–18] have been proposed.
The resaturation of the initially unsaturated repositories is predicted to last hundreds of years after the closure of a repository [19, 20]. During this period, it is inevitable for the saturation/water content of the compacted bentonite to increase by groundwater ingress or decrease by evaporation because of high temperature (50210°C, [21]) originating from decay heat. Thus, the compacted bentonite will experience cyclic wetting/drying processes accompanied by cyclic salinizationdesalination, which could be enhanced by the coupling effects of the sanity groundwater intrusion and the possible water evaporation, or the ingress of fresh groundwater. Therefore, the salt content or fluid concentration in the compacted bentonite will also undergo cyclic processes, which results in alterations of properties of the compacted bentonite [22].
The soil water retention curve (SWRC) has been widely used to evaluate unsaturated soil properties, including saturated and unsaturated permeability functions, shear strength parameters, and volume change [15, 23]. Meanwhile, in unsaturated soils, modeling of water distribution and flow needs knowledge of SWRC, which also plays a key role in water management and predicting transportation of solute and contaminant in the environment [24]. As an important hydraulic characteristic, water retention is closely correlated with the size and connectivity of soil pore space; therefore, it is strongly influenced by the texture and structure of the soil and other components such as soil organic matter [25]. Hence, it is important to have a reasonably accurate and more practical characterization of the SWRC.
It is generally agreed that direct measurement of SWRC in the laboratory requires hard work and is timeconsuming [26]. Therefore, many empirical equations were developed for description of the SWRC in the past decades [24, 27–30]. For compacted bentonite, modified water retention models were also developed with consideration of temperature and dry density effects [31, 32]. These equations are composed of two or three fitting parameters, which are evaluated by using either laboratory results at different moisture contents or statistical relationships between the suction and other basic soil properties.
Lots of factors can influence the SWRC of geomaterials, including soil type, initial water content, mineralogy, density, temperature, texture, stress history, method of compaction, and net confining stress. In this regard, various investigations [33–36] have been conducted. Meanwhile, experimental results also indicated that the salt concentration of infiltrating groundwater can influence the water retention properties of compacted bentonite/expansive clay specimens when used as buffer/backfill materials in the HLW repository [37–39]. For wetting or drying with a given concentration pore fluid, Fredlund and Xing [28] assumed that the maximum osmotic suction in less active soils is approximately 1.5 MPa. It is worth noting that osmotic suction changes may be of importance in some cases. For instance, the effect of osmotic suction is significant when the soil salt content is in a high level [40]. Thyagaraj and Rao [35] also found that the osmotic suction, derived by salt water, should not be ignored in the suction measurement of an Indian bentonite. This observation could be confirmed by Mata et al. [38] with measurement of the water retention curve of MX80 bentonite/sand mixture, in which different concentrations of salt solutions were used to speed up the saturation process. In Europe, the Boom Clay, which is a 100 m thick marine clay, presents favorable properties to limit and delay the migration of the leached radionuclides over extended periods of time. Based on the measurement of suction of Boom Clay containing NaNO_{3} solution using the chilledmirror psychrometer (WP4 Dewpoint Potential Meter) and investigation of the corresponding changes of the microstructure of the specimen tested with the MIP (mercury intrusion porosimetry) technique, Mokni [41] stated that the total suction of specimens saturated with salt solution was higher compared with that of specimens saturated with distilled water, while the matric suction in the salt solutionsaturated group was smaller than that in the distilled watersaturated group.
China’s deep geological disposal program for HLW began in the 1980s. Beishan of Gansu has been selected as the preferred site for construction of repository in China [23, 37]. Bentonite from Gaomiaozi (GMZ), 300 km northwest of Beijing, has been reckoned as the preferred “buffer/backfill material” for building engineering barriers in deep geological repositories [37].
Efforts have been made to thoroughly understand the water retention property of GMZ bentonite in recent years [30, 42–46]. Chen et al. [42] measured the SWRCs of the densely compacted GMZ bentonite at room temperature (20°C) under confined and free swelling conditions using the vapor phase technique and the osmotic technique. Results show that influence of confining conditions on water retention depends on suction, and it becomes less significant as the suction increases from 4 .2MPa. For considering temperature effects, Ye et al. [46] measured SWRCs of GMZ bentonite at different temperatures and reported that temperature effects on the soilwater characteristics depend strongly on the constraint conditions. For unconfined compacted GMZ bentonite specimens, the effect of temperature on the hysteresis behaviour is related to the suction level: a negative correlation between hysteresis behaviour and temperature. Based on these and a comprehensive review of literatures [27, 28, 31, 47], Wan [30] extended the Fredlund and Xing [28] equation for description of water retention curves at room temperature (20°C) to a novel one that can consider temperature effects.
According to in situ investigations in Beishan, the average dissolved solid (TDS) concentration of groundwater is in the range of 3 g/L to 12 g/L, with a maximum up to 50 g/L [48, 49]. Apparently, the influence of salt solution on water retention of GMZ bentonite becomes a key issue to be reckoned with and of significance to China’s HLW disposal project.
In this work, GMZ bentonite powder with different salt contents was obtained after being sprayed with different concentrations of NaCl solutions. The specimens were prepared by compacting the amended bentonite powder to target dimensions and dry density. Suctioncontrolled water retention tests were carried on with the compacted amended GMZ bentonite specimens under confined conditions. Based on these, a modified soil water retention model for the compacted GMZ bentonite with consideration of salt solution concentration effects was established and verified.
2. Materials and Methodology
2.1. Materials and Specimen Preparation
The material tested in this work is GMZ bentonite. The GMZ bentonite was sampled in the northern Chinese Inner Mongolia autonomous region, 300 km northwest from Beijing. There are 160 million tons with 120 million tons of Nabentonite reserves in the deposit (Figure 1), and the mine area is about 72 km^{2}. The GMZ bentonite used in this work belongs to a Natype bentonite. The bulk chemical components of the sample were analyzed using Xray fluorescence spectrometry as follows (in mass): SiO_{2} 67.43%, Al_{2}O_{3} 14.20%, Fe_{2}O_{3} 2.40%, Na_{2}O 1.75%, CaO 1.13%, K_{2}O 0.73%, FeO 0.29%, TiO_{2} 0.12%, MgO 0.10%, P_{2}O_{5} 0.02%, and MnO 0.02%. In addition, GMZ bentonite contains some lanthanides, such as La, Ce, and Nd [50, 51]. Relevant properties of this bentonite are listed in Table 1 [14, 50, 51]. A high cation exchange capacity and adsorption ability can be identified. The XRD test results (Figure 2) show that the characteristic diffraction peaks of the plane at 6.03°, 20.0°, and 35.0° in 2 values and the corresponding basal spacing are 14.70, 4.44, and 2.50 Å, respectively, indicating that the main compound in the specimens tested in the present work is montmorillonite (M). The XRD pattern also indicates the existence of quartz (Q: 26.8°), cristobalite (Cri: 22.2°), and albite (Al: 28.1°) [52].

For most investigations involving chemical effects on the hydromechanical behaviours of soils in the laboratory, various methods were developed for specimen preparation: (1) mixing the soil powder with desired solution, (2) immersing a soil specimen in a solution, and (3) percolating a soil specimen by solution [53]. Mixing the GMZ bentonite powder with desired solutions was employed for the specimen preparation in this work.
Following this aspect, a pretreating course was conducted. The GMZ bentonite powder with an initial water content of 11.7% was dried in an oven at 105°C for 24 h. After that, distilled water and 0.5mol/L, 1.0mol/L, and 1.5mol/L NaCl solutions were sprayed on the bentonite powder to achieve a target water content of 20%, respectively. The mass ratio of salt and bentonite was defined as (mol/g), and corresponding salt contents of the amended bentonite powder are listed in Table 2. After being selfcured in sealed boxes for about 60 h, the amended bentonite powder was dried in an oven at 55°C [54] for another 48 h. After being cooled on a metal plate, the dried specimen was transferred to a mortar and grinded for avoiding particle clusters. Then, the grinded bentonite powder was suctioncontrolled to target suction (CaCl_{2}, 139 MPa) using the vapor phase method at a room temperature of 20°C. After the equilibrium was reached, the amended bentonite powder was compacted into cylindrical specimen with a constant speed at 0.1 mm/min to a target dimension of 6 mm in height, a diameter of 20 mm, and a dry density of 1.70 Mg/m^{3}. The specifications of the compacted specimens used in this work are listed in Table 2.

It should be noted that, as the total salt quantity in the specimen is kept constant during the wetting/drying processes, the current salt concentration of the pore fluid can be determined based on an assumption that there were no chemical interactions between fluid phase and solid phase [55]: where is the target water content (20%), is the concentration of salt solutions (0 .5mol/L, 1 .0mol/L, and 1 .5mol/L) used for the preparation of amended GMZ bentonite powder, and is the current water content.
2.2. Methodology
2.2.1. Suction Control
In this work, the vapor phase technique was employed for suction control (Figure 3), in which total suction was applied. For application of the vapor phase method, the confined compacted GMZ bentonite specimens were placed in a desiccator, which was filled with a target saturated salt solution. Saturated salt solutions and their corresponding controlling suctions are listed in Table 3 [56, 57]. The specimen was regularly weighted until the equilibrium was reached. For each given suction, this process took about 2 to 4 weeks.
(a)
(b)
(c)

The wetting course from 139 MPa to 4.2 MPa, including 6 suctions listed in Table 3, was followed for each of the 4 specimens listed in Table 2.
2.2.2. Determination of Water Content
The gravimetric water content of each specimen was measured by drying for 24 h at 105°C in an oven. Based on this, the corresponding volumetric water content can be calculated by [40] where is the degree of saturation and is the relative density ( in Table 1).
For determination of the influence of pore fluid chemistry on the saturated (“zero air voids” [40]) water content, was assumed as the saturated water content corresponding to matric suction equal to 0 [31]. Following this aspect, specimens with different salt contents were inundated in distilled water. This saturation process lasted about 68 days [58]. After the equilibrium was reached, water content was measured and the saturated volumetric water content was calculated using equation (2).
3. Results
3.1. Influence of Salt Content on SWRC
The measured SWRCs of confined compacted GMZ bentonite with different salt contents in wetting path are shown in Figure 4.
According to Figure 4, the measured volumetric water content decreases with the increasing suction. For a given suction, the measured water content of specimens increases as the salt content (or pore fluid concentration) increases. Furthermore, the influencing rate depends on suction. For lower suctions (lower than 38 MPa), water retention capacity increases with the increase in salt content, while for higher suctions (higher than 38 MPa), the influence was negligible. It can be inferred that, for low suctions, the high saturation induced an increase in macropore size [38, 59], which was further enhanced due to the increase in salt content. However, this phenomenon was faded as the suction increased, which may attribute to the decrease in water content.
3.2. Influence of Pore Fluid Concentration on Saturated Water Content
For the confined compacted GMZ bentonite at a dry density of 1.7 Mg/m^{3} with different salt contents, measured saturated water content θ_{s} versus pore fluid concentration () is presented in Figure 5.
It is obvious that water content θ_{s} increases with the increase in pore fluid salt concentration, while the increasing rate decreases as the salt concentration increases and becomes stable for the fluid salt concentration higher than 0.735 mol/L. Possible explanations to pore fluid chemistry effects on the water content could be that the increase in the concentration of pore fluid results in the decrease in the diffuse double layer of bentonite [60]; i.e., a high water salinity of the pore fluid could induce the contraction of clay aggregates and thus increase the void volume between the aggregates (interaggregate or macropores) [59].
4. A Modified SWRC Model with Consideration of Salt Solution Effects
4.1. Equation for SWRC
As mentioned before, due to difficulties and labor costs for measuring SWRC, many empirical models for SWRC have been developed. Actually, among all the models available, the commonly used ones are those proposed by van Genuchten [29], Fredlund and Xing [28], and Houston et al. [61], while the Fredlund and Xing model appears to be the most appropriate one for fitting the experimental data [62].
The Fredlund and Xing model [28] can be expressed as follows: where is the suction and is a correction function that can be defined as where is the suction corresponding to the residual water content.
Parameter , a suction that is closely related with the SWRC’s inflection point, can be considered as the air entry value when parameter is small, even if the parameter value is usually bigger than that of air entry [28]. Parameter makes the curve towards the region of higher soil suction as its value increases but has no effect on the overall shape of curves [63]. Parameter , which represents the curve asymmetry around the inflection point, is associated with the residual water content. Parameter affects the slope of SWRC in the desaturation region. The value of is associated with pore size distribution, which always shows a positive relationship with the uniformity of the pore sizes in the soil [63, 64].
4.2. A Modified SWRC Model with Consideration of Pore Fluid Concentration Effects
4.2.1. Determination of Parameters
(1) Determination of with Consideration of Pore Fluid Concentration. To consider the pore fluid concentration effects in equation (3), influence of the pore fluid concentration on the saturated water content should be included.
It can be observed from Figure 5 that the measured relationship between the water content () and the pore fluid concentration () looks like an asymptotic line. Therefore, a secondorder equation (5) was proposed for fitting the saturated water content and pore fluid concentration at the saturated state.
Integrating equation (5) from , to and (where , which was calculated by equation (2) using the gravimetric water content measured at matric suction 0 MPa, and , which was assumed as the maximum value of volumetric water content of specimens tested), yields, where is the pore fluid concentration of chemical components corresponding to matric suction MPa, is a coefficient related to salt concentration, and is the saturated water content at the current chemical concentration. The values of and saturated water content can be calculated using equations (2) and (6).
With salt concentrations 0.01, 0.384, 0.735, and 2.0 mol/L and the corresponding measured saturated water contents in Figure 5, and can be fitted using equation (6) and substituted back into equation (6), which gives
With equation (7), the relationship between and can be depicted in Figure 5. Figure 5 indicates that the secondorder equation (5) proposed in this work can satisfactorily describe the saturated water content with consideration of pore fluid chemistry.
(2) Determination of with Consideration of Pore Fluid Concentration. To consider the pore fluid concentration effects in equation (3), the effects of the pore fluid concentration on the air entry value (parameter ) should be analyzed as well.
In this respect, Karagunduz et al. [65] once proposed an approximate linear relationship (equation (8)) between the air entry valuerelated parameter and the chemical concentration of the pore fluid. where is the pore fluid concentration corresponding to matric suction MPa, is a concentrationrelated coefficient, and is the air entry valuerelated parameter of specimen hydrated with distilled water.
The fitted relationship between the air entry valuerelated parameter and the pore fluid concentration for confined compacted GMZ bentonite is presented in Figure 6. Results show that a very good linear relationship can be observed.
With the relationship between the air entry valuerelated parameter and the pore fluid concentration in Figure 6, parameter can be fitted using equation (8) with a value of 0.7734 for confined compacted GMZ bentonite.
(3) Determination of with Consideration of Pore Fluid Concentration. In equation (3), reflects the slope of the SWRC in the desaturation zone [30], which is determined by the pore size distribution of the specimen.
The concentration of infiltration pore fluid will influence pore size distribution, which results in variation of the value and the shape of SWRC. This phenomenon has been observed and traditionally explained by means of the diffuse doublelayer theory (DDL) [66]. Thickness of the diffuse double layer decreases with increasing concentration of electrolytes or their valence in the pore fluid, which results in the contraction of clay aggregates locally and the corresponding increase in connecting pores in the clay. Therefore, although the use of hypotheses assumed by the doublelayer theory is still restrictive and merely for engineering purposes [60, 67], it is still helpful and reasonable to interpret most of the relations between pore fluid chemistry and macroscopic behaviour of clays [41, 68].
The influence of pore fluid concentration on water retention properties of clay has succeeded to attract sufficient attention by researchers [38, 39]. Among the literatures, the work of Ravi and Rao [39] provides a noteworthy contribution in this area. The measurement of SWRCs of bentonite/sand mixture specimens with different dry densities, respectively, in the distilled water group and NaCl solution group, was reported in the literature. Results showed that the concentration of infiltration pore fluid has a significant effect on SWRCs. On the basis of the measured data, the relationship between the slope value of SWRCs of bentonite/sandmixed specimens and the corresponding concentrations of pore fluid was fitted and depicted, as shown in Figure 7.
Based on the fitted results in Figure 7, the relationship between the slope values and the concentrations of pore fluid could be assumed as where is a concentrationrelated coefficient to be fitted, is the pore fluid concentration corresponding to matric suction MPa, and is the slope value of SWRC of specimens hydrated with distilled water.
With SWRCs of confined compacted GMZ bentonite specimens with dry density of 1.7 Mg/m^{3} tested using distilled water and different concentrations of NaCl solutions in this work, the relationship between the slope values of SWRCs and the corresponding pore fluid concentrations is shown in Figure 7. Parameter was fitted with a value of 1.113. Therefore, the relationship between the slope values of SWRCs and pore fluid concentrations for confined compacted GMZ bentonite can be expressed as
It can be observed that the value of the concentrationrelated parameter of 1.5 Mg/m^{3} () bentonite/sand mixture specimens is larger than that of 1.7 Mg/m^{3} GMZ bentonite specimen (). This observation indicates that the influence of pore fluid concentrations on the slope values of SWRCs depends on dry density. A possible explanation to this observation could be that the specimens with lower dry density have more pore space for pore fluid intrusion, which enhanced the increase in the interaggregate or macropores induced by the contraction of clay aggregates [59]. However, for compacted bentonite with higher dry density, less porefluid intrusion is likely to be insufficient to form enough DDL in a small interlayer space [37, 69].
Comparison results indicate that the proposed exponential function (equation (9)) can satisfactorily describe the relationship of the slope values of SWRCs and the pore fluid concentrations.
4.2.2. A Modified SWRC Model with Consideration of Pore Fluid Chemistry
Substituting equations (6), (8), and (10) into equation (3), equation (11) can be obtained. where , , , and can be calculated using equations (4), (6), (8), and (10), respectively.
It should be noted that, considering that the influence of the pore fluid concentration on the residual water content of compacted clays is insignificant (Figure 4), the influence of pore fluid chemistry on parameter is safely neglected in this model.
Deduced from the Fredlund and Xing [28] SWRC equation, equation (11) is the modified model that can be used for predicting the SWRCs of confined compacted GMZ bentonite with consideration of pore fluid chemistry.
It should be noted that, in equation (11), for calculation of using equation (4), the residual suction should be determined. Commonly, was practically determined or assumed that the maximum suction is applied, which is corresponding to a lower enough water content [30]. For instance, Phayak [70] chose the suction corresponding to 3.51% water content as the residual suction . Wan [30] assumed the maximum suction of 309 MPa applied at 20°C to be the residual suction corresponding to a water content of about 5% for compacted GMZ bentonite. In this study, the maximum suction of 139 MPa applied is assumed to be the residual suction of GMZ bentonite, which is corresponding to a water content around 8%.
4.2.3. Verification
Figure 8 shows the measured SWRCs of the compacted GMZ bentonite hydrated with different salt concentrations (, , , and ) and distilled water. Based on the measured SWRCs of , , and , as well as equations (1), (2), and (11), parameters , , , and in equation (11) for the compacted GMZ bentonite with 1.7 Mg/m^{3} dry density can be determined: , , , and . Substituting , , , and into equation (11) gives where
Equation (12) is the modified model for describing the SWRCs of the confined compacted GMZ bentonite taking account of the pore fluid concentration.
With equation (12), the SWRC of specimen was calculated. According to the contrast between the calculated results and measurement data (Figure 8), the proposed model can well describe the effects of the pore fluid chemistry process on the SWRCs of confined compacted GMZ bentonite.
5. Conclusions
In this work, the SWRCs of confined compacted GMZ bentonite with different salt contents were attained by a series of suctioncontrolled tests. The effects of pore fluid chemistry on the water retention properties of GMZ bentonite were analyzed. A modified model considering the pore fluid chemistry effects for describing the SWRCs of confined compacted GMZ bentonite was proposed and verified. The main conclusions can be drawn as follows.
For the compacted GMZ bentonite, the effect of pore fluid chemistry on the water retention properties depends on suction. For a given suction, the measured water content of specimens increases with the increase in pore fluid concentration. Furthermore, the influencing rate depends on suction.
With the Fredlund and Xing equation, parameters including the air entry valuerelated parameter , the saturated water content , and the slope value of SWRCs were determined with consideration of the influence of the pore fluid concentration. Based on these, a modified model with consideration of pore fluid chemistry effects was proposed to describe the SWRCs of the compacted GMZ bentonite. The contrast of the calculated results and the experimental data demonstrates that the modified model can satisfactorily describe the water retention properties of compacted GMZ bentonite with consideration of pore fluid chemistry effects.
It should be noted that the proposed model only verified with test results of confined compacted GMZ bentonite. Further investigations should be done in the coming work for extending its applications to related geomaterials.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors are grateful to the National Natural Science Foundation of China (Project Nos. 41807253, 41672271, and 41527801) for the financial support. The authors also wish to acknowledge the support of the European Commission via the Marie Curie IRSES project GREAT—Geotechnical and geological Responses to climate change: Exchanging Approaches and Technologies on a worldwide scale (FP7PEOPLE2013IRSES612665). The authors are also grateful to the Research Fund Program of the Key Laboratory of Geotechnical and Underground Engineering (Tongji University), Ministry of Education (Project No. KLETJGEB1803).
References
 ENRESA, FullScale Engineered Barriers Experiment for Deep Geological Repository for HighLevel Radioactive Waste in Crystalline Rock (FEBEX Project), EUR 19147, Nuclear Science and Technology Series, Luxembourg, 2000.
 Y. He, Y. G. Chen, and W. M. Ye, “Equilibrium, kinetic, and thermodynamic studies of adsorption of Sr(II) from aqueous solution onto GMZ bentonite,” Environmental Earth Sciences, vol. 75, no. 9, p. 807, 2016. View at: Publisher Site  Google Scholar
 JNC, H12 Project to Establish Technical Basis for HLW Disposal in Japan, Support Report 2, Japan Nuclear Cycle Development Institute, 1999.
 M. J. Kim, S. R. Lee, J. S. Jeon, and S. Yoon, “Sensitivity analysis of bentonite buffer peak temperature in a highlevel waste repository,” Annals of Nuclear Energy, vol. 123, pp. 190–199, 2019. View at: Publisher Site  Google Scholar
 D. Marcial, P. Delage, and Y. J. Cui, “On the high stress compression of bentonites,” Canadian Geotechnical Journal, vol. 39, no. 4, pp. 812–820, 2002. View at: Publisher Site  Google Scholar
 P. L. Martín, J. M. Barcala, and F. Huertas, “Largescale and longterm coupled thermohydromechanic experiments with bentonite: the FEBEX mockup test,” Journal of Iberian Geology, vol. 32, no. 2, pp. 259–282, 2006. View at: Google Scholar
 R. Pusch and R. Weston, “Superior techniques for disposal of highly radioactive waste (HLW),” Environmental Earth Sciences, vol. 73, no. 9, pp. 5219–5231, 2015. View at: Publisher Site  Google Scholar
 S. M. Rao and K. Ravi, “Influence of initial degree of saturation on swell pressures of compacted Barmer bentonite specimens,” Annals of Nuclear Energy, vol. 80, pp. 303–311, 2015. View at: Publisher Site  Google Scholar
 SKB, “SR97 Waste, repository design and sites,” Tech. Rep., SKB Technical Report TR9908, Stockholm, Sweden, 1999. View at: Google Scholar
 A. Dutt, M. S. Saini, T. N. Singh, A. K. Verma, and R. K. Bajpai, “Analysis of thermohydrologicmechanical impact of repository for highlevel radioactive waste in clay host formation: an Indian reference disposal system,” Environmental Earth Sciences, vol. 66, no. 8, pp. 2327–2341, 2012. View at: Publisher Site  Google Scholar
 A. P. Fraser Harris, C. I. McDermott, O. Kolditz, and R. S. Haszeldine, “Modelling groundwater flow changes due to thermal effects of radioactive waste disposal at a hypothetical repository site near Sellafield, UK,” Environmental Earth Sciences, vol. 74, no. 2, pp. 1589–1602, 2015. View at: Publisher Site  Google Scholar
 M. Xie, S. Bauer, O. Kolditz, T. Nowak, and H. Shao, “Numerical simulation of reactive processes in an experiment with partially saturated bentonite,” Journal of Contaminant Hydrology, vol. 83, no. 12, pp. 122–147, 2006. View at: Publisher Site  Google Scholar
 J. Kasbohm, R. Pusch, L. NguyenThanh, and T. HoangMinh, “Labscale performance of selected expandable clays under HLW repository conditions,” Environmental Earth Sciences, vol. 69, no. 8, pp. 2569–2579, 2013. View at: Publisher Site  Google Scholar
 W. M. Ye, M. Wan, B. Chen, Y. G. Chen, Y. J. Cui, and J. Wang, “Temperature effects on the unsaturated permeability of the densely compacted GMZ01 bentonite under confined conditions,” Engineering Geology, vol. 126, pp. 1–7, 2012. View at: Publisher Site  Google Scholar
 W. M. Ye, M. Wan, B. Chen, Y. G. Chen, Y. J. Cui, and J. Wang, “Temperature effects on the swelling pressure and saturated hydraulic conductivity of the compacted GMZ01 bentonite,” Environmental Earth Sciences, vol. 68, no. 1, pp. 281–288, 2013. View at: Publisher Site  Google Scholar
 P. J. Cleall, T. A. Melhuish, and H. R. Thomas, “Modelling the threedimensional behaviour of a prototype nuclear waste repository,” Engineering Geology, vol. 85, no. 12, pp. 212–220, 2006. View at: Publisher Site  Google Scholar
 A. Lloret and M. V. Villar, “Advances on the knowledge of the thermohydromechanical behaviour of heavily compacted “FEBEX” bentonite,” Physics and Chemistry of the Earth, Parts A/B/C, vol. 32, no. 814, pp. 701–715, 2007. View at: Publisher Site  Google Scholar
 J. Rutqvist, D. Barr, J. T. Birkholzer et al., “Results from an international simulation study on coupled thermal, hydrological, and mechanical processes near geological nuclear waste repositories,” Nuclear Technology, vol. 163, no. 1, pp. 101–109, 2008. View at: Publisher Site  Google Scholar
 L. H. Johnson, D. M. LeNeveu, D. W. Shoesmith et al., The Disposal of Canada’s Nuclear Fuel Waste: the Vault Model for Postclosure Assessment, AECL Research Document No. AECL10714, COG938, Atomic Energy of Canada Limited, Pinawa, Manitoba, Canada, 1994.
 R. Senger, P. Marschall, and S. Finsterle, “Investigation of twophase flow phenomena associated with corrosion in an SF/HLW repository in Opalinus Clay, Switzerland,” Physics and Chemistry of the Earth, Parts A/B/C, vol. 33, pp. S317–S326, 2008. View at: Publisher Site  Google Scholar
 S. GarcíaGarcía, M. Jonsson, and S. Wold, “Temperature effect on the stability of bentonite colloids in water,” Journal of Colloid and Interface Science, vol. 298, no. 2, pp. 694–705, 2006. View at: Publisher Site  Google Scholar
 E. Castellanos, M. V. Villar, E. Romero, A. Lloret, and A. Gens, “Chemical impact on the hydromechanical behaviour of highdensity FEBEX bentonite,” Physics and Chemistry of the Earth, Parts A/B/C, vol. 33, pp. S516–S526, 2008. View at: Publisher Site  Google Scholar
 W. M. Ye, Y. J. Cui, L. X. Qian, and B. Chen, “An experimental study of the water transfer through confined compacted GMZ bentonite,” Engineering Geology, vol. 108, no. 34, pp. 169–176, 2009. View at: Publisher Site  Google Scholar
 F. J. Leij, W. B. Russell, and S. M. Lesch, “Closedform expressions for water retention and conductivity data,” Ground Water, vol. 35, no. 5, pp. 848–858, 1997. View at: Publisher Site  Google Scholar
 M. Tuller and D. Or, “Retention of water in soil and the soil water characteristic curve,” Encyclopedia of Soils in the Environment, vol. 4, pp. 278–289, 2004. View at: Google Scholar
 A. Hildenbrand, S. Schlömer, and B. M. Krooss, “Gas breakthrough experiments on finegrained sedimentary rocks,” Geofluids, vol. 2, no. 1, 23 pages, 2002. View at: Publisher Site  Google Scholar
 R. H. Brooks and A. T. Corey, Hydraulic Properties of Porous Media, Hydrol. Pap. 3, Colorado State University, Fort Collins, CO, USA, 1964.
 D. G. Fredlund and A. Xing, “Equations for the soilwater characteristic curve,” Canadian Geotechnical Journal, vol. 31, no. 4, pp. 521–532, 1994. View at: Publisher Site  Google Scholar
 M. T. H. van Genuchten, “A closedform equation for predicting the hydraulic conductivity of unsaturated soils,” Soil Science Society of America Journal, vol. 44, no. 5, p. 892, 1980. View at: Publisher Site  Google Scholar
 M. Wan, Study on SoilWater Characteristics and Permeability of Highly Compacted GMZ Bentonite with Temperature Control, [Ph.D. Thesis], Tongji University, 2010.
 A. Jacinto, M. Villar, R. Gomezespina, and A. Ledesma, “Adaptation of the van Genuchten expression to the effects of temperature and density for compacted bentonites,” Applied Clay Science, vol. 42, no. 34, pp. 575–582, 2009. View at: Publisher Site  Google Scholar
 M. V. Villar, R. GómezEspina, and A. Lloret, “Experimental investigation into temperature effect on hydromechanical behaviours of bentonite,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 2, no. 1, pp. 71–78, 2010. View at: Google Scholar
 N. Lu and W. J. Likos, Unsaturated Soil Mechanics, Wiley, New York, NY, USA, 2004.
 T. M. Thu, H. Rahardjo, and E. C. Leong, “Soilwater characteristic curve and consolidation behavior for a compacted silt,” Canadian Geotechnical Journal, vol. 44, no. 3, pp. 266–275, 2007. View at: Publisher Site  Google Scholar
 T. Thyagaraj and S. M. Rao, “Influence of osmotic suction on the soilwater characteristic curves of compacted expansive clay,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 136, no. 12, pp. 1695–1702, 2010. View at: Publisher Site  Google Scholar
 M. V. Villar and A. Lloret, “Influence of temperature on the hydromechanical behaviour of a compacted bentonite,” Applied Clay Science, vol. 26, no. 14, pp. 337–350, 2004. View at: Publisher Site  Google Scholar
 Y. He, Y.G. Chen, W.M. Ye, B. Chen, and Y.J. Cui, “Influence of salt concentration on volume shrinkage and water retention characteristics of compacted GMZ bentonite,” Environmental Earth Sciences, vol. 75, no. 6, p. 535, 2016. View at: Publisher Site  Google Scholar
 C. Mata, E. Romero, and A. Ledesma, “Hydrochemical effects on water retention in bentonitesand mixtures,” in Proceeding of the 3rd International Conference on Unsaturated soil, J. F. T. Jucá, T. M. P. Campos, and F. A. M. Marinho, Eds., vol. 1, pp. 283–288, Recife, Brazil, 2002. View at: Google Scholar
 K. Ravi and S. M. Rao, “Influence of infiltration of sodium chloride solutions on SWCC of compacted bentonite–sand specimens,” Geotechnical and Geological Engineering, vol. 31, no. 4, pp. 1291–1303, 2013. View at: Publisher Site  Google Scholar
 D. G. Fredlund and H. Rahardjo, Soil Mechanics for Unsaturated Soils, Wiley, New York, NY, USA, 1993. View at: Publisher Site
 N. Mokni, Deformation and Flow Driven by Osmotic Processes in Porous Materials, [Ph.D. Thesis], Department of Geotechnical Engineering and Geoscience (ETCG) Technical University of Catalunya, UPC, Barcelona, Spain, 2011.
 B. Chen, L. X. Qian, W. M. Ye, Y. J. Cui, and J. Wang, “Soil water characteristic curves of Gaomiaozi bentonite,” Chinese Journal of Rock Mechanics and Engineering, vol. 25, no. 4, pp. 788–794, 2006. View at: Google Scholar
 Y. He, W. M. Ye, Y. G. Chen, B. Chen, B. Ye, and Y. J. Cui, “Influence of pore fluid concentration on water retention properties of compacted GMZ01 bentonite,” Applied Clay Science, vol. 129, pp. 131–141, 2016. View at: Publisher Site  Google Scholar
 W. M. Ye, N. C. Borrell, J. Y. Zhu, B. Chen, and Y. G. Chen, “Advances on the investigation of the hydraulic behavior of compacted GMZ bentonite,” Engineering Geology, vol. 169, pp. 41–49, 2014. View at: Publisher Site  Google Scholar
 W. M. Ye, Y. G. Chen, B. Chen, Q. Wang, and J. Wang, “Advances on the knowledge of the buffer/backfill properties of heavilycompacted GMZ bentonite,” Engineering Geology, vol. 116, no. 12, pp. 12–20, 2010. View at: Publisher Site  Google Scholar
 W. M. Ye, M. Wan, B. Chen, Y. G. Chen, Y. J. Cui, and J. Wang, “Effect of temperature on soilwater characteristics and hysteresis of compacted Gaomiaozi bentonite,” Journal of Central South University of Technology, vol. 16, no. 5, pp. 821–826, 2009. View at: Publisher Site  Google Scholar
 M. Sanchez, ThermoHydroMechanical Coupled Analyses in Low Permeability Media, [Ph.D. Thesis], Universitat Politecnica de Catalunya, Barcelona, Spain, 2004.
 Y. H. Guo, T. X. Yang, and S. F. Liu, “Hydrogeological characteristics of Beishan preselected area, Gansu province for China’s highlevel radioactive waste repository,” Uranium Geology, vol. 17, no. 3, pp. 184–189, 2001. View at: Google Scholar
 J. Wang, L. Chen, R. Su, and X. Zhao, “The Beishan underground research laboratory for geological disposal of highlevel radioactive waste in China: planning, site selection, site characterization and in situ tests,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 10, no. 3, pp. 411–435, 2018. View at: Publisher Site  Google Scholar
 Z. J. Wen, “Selection and basic properties of China’s buffer materials for high level radioactive waste repository,” Acta Petrologica et Mineralogica, vol. 24, no. 6, pp. 583–586, 2005. View at: Google Scholar
 Z. J. Wen, “Selection and basic properties of the buffer material for highlevel radioactive waste repository in China,” Acta Geologica Sinica  English Edition, vol. 82, no. 5, pp. 1050–1055, 2008. View at: Publisher Site  Google Scholar
 Y. He, W. M. Ye, Y. G. Chen, and Y. J. Cui, “Effects of K^{+} solutions on swelling behavior of compacted GMZ bentonite,” Engineering Geology, vol. 249, pp. 241–248, 2019. View at: Publisher Site  Google Scholar
 X. P. Nguyen, Y. J. Cui, A. M. Tang, Y. F. Deng, X. L. Li, and L. Wouters, “Effects of pore water chemical composition on the hydromechanical behavior of natural stiff clays,” Engineering Geology, vol. 166, pp. 52–64, 2013. View at: Publisher Site  Google Scholar
 O. Karnland, S. Olsson, U. Nilsson, and P. Sellin, “Experimentally determined swelling pressures and geochemical interactions of compacted Wyoming bentonite with highly alkaline solutions,” Physics and Chemistry of the Earth, Parts A/B/C, vol. 32, no. 17, pp. 275–286, 2007. View at: Publisher Site  Google Scholar
 C. Mata, Hydraulic Behavior of Bentonite Based Mixtures in Engineered Barriers: the Backfill and Plug Test at the ÄSPÖ HRL (Sweden), [Ph.D. Thesis], Technical University of Catalonia Department of Geotechnical Engineering and Geosciences School of Civil Engineering, Barcelona, Spain, 2003.
 P. Delage, M. D. Howat, and Y. J. Cui, “The relationship between suction and swelling properties in a heavily compacted unsaturated clay,” Engineering Geology, vol. 50, no. 12, pp. 31–48, 1998. View at: Publisher Site  Google Scholar
 A. M. Tang and Y. J. Cui, “Controlling suction by the vapour equilibrium technique at different temperatures and its application in determining the water retention properties of MX80 clay,” Canadian Geotechnical Journal, vol. 42, no. 1, pp. 287–296, 2005. View at: Publisher Site  Google Scholar
 C. M. Zhu, W. M. Ye, Y. G. Chen, B. Chen, and Y. J. Cui, “Influence of salt solutions on the swelling pressure and hydraulic conductivity of compacted GMZ01 bentonite,” Engineering Geology, vol. 166, pp. 74–80, 2013. View at: Publisher Site  Google Scholar
 G. Musso, E. Romero, and G. D. Vecchia, “Doublestructure effects on the chemohydromechanical behaviour of a compacted active clay,” Géotechnique, vol. 63, no. 3, pp. 206–220, 2013. View at: Publisher Site  Google Scholar
 T. A. Hueckel, “Water–mineral interaction in hygromechanics of clays exposed to environmental loads: a mixture theory approach,” Canadian Geotechnical Journal, vol. 29, no. 6, pp. 1071–1086, 1992. View at: Publisher Site  Google Scholar
 W. N. Houston, H. B. Dye, C. E. Zapata, Y. Y. Perera, and A. Harraz, “Determination of SWCC using one point suction measurement and standard curves,” in Geotechnical Special Publication, pp. 1482–1493, Reston, VA, USA, 147 edition, 2006. View at: Google Scholar
 S. Nam, M. Gutierrez, P. Diplas et al., “Comparison of testing techniques and models for establishing the SWCC of riverbank soils,” Engineering Geology, vol. 110, no. 12, pp. 1–10, 2010. View at: Publisher Site  Google Scholar
 G. L. S. Babu, J. Peter, M. Mukesh, and E. Gartung, “Significance of soil suction and soil water characteristic curve parameters,” Geotechnical Testing Journal, vol. 28, no. 1, pp. 102–107, 2005. View at: Publisher Site  Google Scholar
 W. C. Sillers, D. G. Fredlund, and H. Rahardjo, “Mathematical attributes of some soil—water characteristic curve models,” in Unsaturated Soil Concepts and Their Application in Geotechnical Practice, pp. 243–283, Springer, Dordrecht, Netherlands, 2001. View at: Publisher Site  Google Scholar
 A. Karagunduz, K. D. Pennell, and M. H. Young, “Influence of a nonionic surfactant on the water retention properties of unsaturated soils,” Soil Science Society of America Journal, vol. 65, no. 5, pp. 1392–1399, 2001. View at: Publisher Site  Google Scholar
 J. K. Mitchell, Fundamentals of Soil Behavior, John Wiley and Sons, New York, NY, USA, 2nd edition, 1993.
 Y. He, Y. J. Cui, W. M. Ye, and N. Conil, “Effects of wettingdrying cycles on the air permeability of compacted Téguline clay,” Engineering Geology, vol. 228, pp. 173–179, 2017. View at: Publisher Site  Google Scholar
 H. Komine and N. Ogata, “Prediction for swelling characteristics of compacted bentonite,” Canadian Geotechnical Journal, vol. 33, no. 1, pp. 11–22, 1996. View at: Publisher Site  Google Scholar
 R. Pusch and R. N. Yong, Microstructure of Smectite Clays and Engineering Performance, Taylor & Francis, New York, NY, USA, 2006.
 T. Phayak, Experimental Investigations on SmallStrain Stiffness Properties of Partially Saturated Soils via Resonant Column and Bender Element Testing, [Ph.D. Thesis], University of Texas at Arlington, 2006.
Copyright
Copyright © 2019 Yong He et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.