Next Article in Journal
Geochemistry, Mineralogy and Textural Properties of the Lower Globigerina Limestone Used in the Built Heritage
Next Article in Special Issue
Ex and In Situ Reactivity and Sorption of Selenium in Opalinus Clay in the Presence of a Selenium Reducing Microbial Community
Previous Article in Journal
A Mineralogical, Geochemical, and Geochronological Study of ‘Valencianite’ from La Valenciana Mine, Guanajuato, Mexico
Previous Article in Special Issue
Mineralogical Analysis of Bentonite from the ABM5 Heater Experiment at Äspö Hard Rock Laboratory, Sweden
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Swelling Pressure and Permeability of Compacted Bentonite from 10th Khutor Deposit (Russia)

by
Artur Yu. Meleshyn
1,
Sergey V. Zakusin
2,3 and
Victoria V. Krupskaya
2,4,*
1
Gesellschaft für Anlagen- und Reaktorsicherheit (GRS) gGmbH, 38122 Braunschweig, Germany
2
Institute of Ore Geology, Petrography, Mineralogy and Geochemistry, Russian Academy of Science (IGEM RAS), 119017 Moscow, Russia
3
Department of Engineering and Environmental Geology, Lomonosov Moscow State University, 119991 Moscow, Russia
4
Nuclear Safety Institute, Russian Academy of Science (IBRAE RAS), 115191 Moscow, Russia
*
Author to whom correspondence should be addressed.
Minerals 2021, 11(7), 742; https://doi.org/10.3390/min11070742
Submission received: 5 May 2021 / Revised: 21 June 2021 / Accepted: 28 June 2021 / Published: 7 July 2021

Abstract

:
Bentonites from the 10th Khutor deposit (Republic of Khakassia, Russia) are considered a potential buffer material for isolation of radioactive waste in the crystalline rocks of Yeniseyskiy site (Krasnoyarskiy region). This study presents the results of a series of permeameter experiments with bentonite compacted to dry densities of 1.4, 1.6, and 1.8 g/cm3, saturated and permeated by the artificial groundwater from Yeniseyskiy Site. Permeation was conducted at hydraulic gradients of 180–80,000 m/m to simulate potential hydraulic conditions in the early post-closure phase of a deep geological repository (DGR). The respective swelling pressures of 0.8 ± 0.3, 2.2 ± 0.6, and 6.3 ± 0.3 MPa and permeabilities of (27 ± 15) × 10−20, (3.4 ± 0.8) × 10−20, and (0.96 ± 0.26) × 10−20 m2 were observed for the hydraulic gradient of 2000 m/m, which is recommended for the determination of undisturbed swelling pressures and permeabilities in permeameter experiments. Upon incremental increases in the hydraulic gradient, swelling pressures at all densities and permeability at the density of 1.8 g/cm3 remained unchanged, whereas permeabilities at 1.4 and 1.6 g/cm3 decreased overall by a factor of approximately 5 and 1.7, respectively. Seepage-induced consolidation and/or reorganisation of bentonite microstructure are considered possible reasons for these decreases.

1. Introduction

In the Russian Federation, a concept is being developed to isolate class 1 and 2 high-level radioactive waste and spent nuclear fuel in crystalline rocks (Yeniseyskiy Site, Krasnoyarskiy region) at depths of about 500–700 m [1,2,3,4]. In this case, the safety function of radioactive waste disposal is provided by a passive multi-barrier system that ensures the isolation of radioactive waste from the biosphere [5,6,7].
The current Russian sealing concept of a deep geological radioactive waste disposal facility (DGR) in crystalline rocks foresees the use of compacted bentonite as a buffer material enclosing containers with radioactive waste [1,2,3,4] due to the stability of its properties and ability to prevent the mass transfer of radionuclides from advection and diffusion mechanisms [2,3,4,5,6,7]. The 10th Khutor bentonite deposit is located relatively close to the area of the deep geological disposal site and considered as a source for the potential buffer material in study [4] based on its favourable sealing properties and deposit resources. Swelling pressure and permeability (or proportional to that hydraulic conductivity) of compacted bentonite are key parameters in determining the sealing performance of bentonite-based geotechnical barriers. A direct determination of these parameters for a realistic buffer configuration under in situ conditions of DGR would require conducting sophisticated, long-term, and expensive experiments in an Underground Research Laboratory. Laboratory experiments with small specimens of compacted bentonite provide a feasible way to obtain estimates of the parameters necessary to conduct preliminary performance assessments of bentonite-based geotechnical barriers and facilitate the planning and optimisation in situ studies, possibly reducing their extent.
During the construction of geotechnical barriers, unsaturated bentonites will experience high fluid pressure differences—up to several MPa—which may result across the buffer in the course of establishing hydraulic conditions characteristic of the host formation. Therefore, fluid pressure differences of up to 8 MPa across the compacted specimens were imposed on the present study in order to test the sealing performance of a bentonite sample from the 10th Khutor deposit (described as “10X” or “10X bentonite”). Swelling pressure and hydraulic conductivity of 10X were estimated in a series of preliminary laboratory experiments with distilled water [4]. The present study provides an independent laboratory estimation of swelling pressure and permeability for 10X upon saturation and permeation with the synthetic groundwater of crystalline rocks at the Yeniseyskiy site.

2. Materials and Methods

2.1. Geological Site and Composition of 10X

The 10th Khutor deposit is located in the steppe area 8 km southwest of Chernogorsk city (Republic of Khakassia). It has reserves of about 3–4 million tons and belongs to a group of deposits with total reserves of approximately 13 million [8]. The deposit has a Carboniferous age (C2) expressed in several bentonite layers interbedded with clay rocks and carbonaceous layers [8,9,10]. Its genesis is considered volcanogenic-sedimentary.
In the process of studying bentonite powder with specified technical conditions, obtained characteristics were required to justify the use of clay material as a component of engineering safety barriers. This material was assigned a commercial name,10X.
The study of particle size distribution was conducted using laser deflation on an Analysette-22 (Fritsch GmbH, Idar-Oberstein, Germany) particle analyser in suspension for the measurement range of 0.1–300 μm in five different positions of the measuring cell and 100 measurements in each position. The results are shown in Table 1.
The mineral composition of the 10X bentonite [10,11] was estimated by X-ray diffraction using an Ultima-IV (Rigaku, Tokyo, Japan) X-ray diffractometer (Cu-Kα radiation, detector—D/Tex-Ultra, scan range 3–65° 2θ). The quantitative mineral analysis (Table 1) was conducted using the Rietveld method [12] with PROFEX GUI ver. 3.1.3, (Nicola Doeblin, Solothurn, Switzerland) software for BGMN [13].
The chemical composition (Table 1) was determined using X-ray fluorescence analysis using an Axios mAX (PANalytical, Almelo, the Netherlands) according to standard methods [11]. The hygroscopic water content of 8.94% was determined upon drying at 105 °C. The suction of the 10X bentonite was not determined.
The density of the solid phase was measured using the pycnometer method with kerosene following [14,15] and amounted to 2.8 g/cm3.
The value of the cation exchange capacity was determined from the adsorption of the Cu(II)-trien-complex using the spectrophotometric method [16,17,18] and amounted to 73.5 meq/100g. The composition of the exchangeable complex is mainly Na-Mg, and comprises Na—23.6, Ca—16.6, Mg—30.3, and K—1.2 meq/100g, measured with atomic absorption spectrometry (AAS) after repeated extraction with ammonium acetate solution.

2.2. Experiments and Conditions

For the research, triplicate experiments of 10X bentonite at the dry densities of 1.4, 1.6, and 1.8 g/cm3 were conducted using nine constant-head, rigid-wall permeameter cells made of titanium (Figure 1a). Compacted specimens were prepared by pressing directly into permeameter cells. Bentonite powder (Table 1a) was filled into a permeameter between two sintered steel filters (GKN SIKA-R 20 AX, a porosity of 43%, an average pore diameter of 23 µm, with an average permeability—measured in triplicate using three permeameter cells—of (3.4 ± 0.5) × 10−14 m2) and, after installing a load-transfer ram and a press flange above the top sintered filter, it was statically compacted (without adding any water) to a pellet specimen with a diameter of 5 cm, a height of 1 cm, and a dry density of 1.4, 1.6, and 1.8 g/cm3. The cap nut of the press flange was then tightened without adding axial stress. Based on the density of soil solids (2.80 g/cm3), porosity of specimens can be estimated to equal 0.50, 0.43, and 0.36, pore volumes to 9.8, 8.4, and 7.0 cm3, and saturation to 25%, 33%, and 45%, respectively. Respective bulk densities equal 2.02, 2.17, and 2.32 g/cm3.
The swelling pressure of the specimens in the permeameter cells was measured continuously by force transducers FKA613 (Ahlborn Mess- und Regelungstechnik GmbH, Holzkirchen, Germany). It should be noted that for the specimen in the ninth cell (at 1.8 g/cm3, the leftmost cell in Figure 1b), a behaviour strongly deviating from the other two specimens at 1.8 g/cm3 was observed, as its swelling pressure exceeded the measuring capacity of force transducers at 10.2 MPa within two days of starting the experiment. Two restarts and additional tests suggested a permeameter malfunction as the reason for this deviating behaviour; therefore, the experiment with the ninth cell was discontinued.
Since no reliable data currently exists on the composition of the groundwater of the Yeniseyskiy Site at the target depth, synthetic groundwater (hereafter referred to as SGW) was used in the experiments and developed based on available hydrogeological data and data in [19], which is characterized by a pH of 8.1, ionic strength of 6.5 mmol/L, and ionic composition: Na+—0.760, Mg2+—0.494, Ca2+—1.217, K+—0.115, Cl—2.549, SO42−—0.494, and HCO3—0.76 mmol/L.
Compacted bentonite specimens were allowed to saturate by absorbing SGW from burettes connected to inlets at cell bases to let trapped air escape through outlets at cell tops for 28 days, then to inlets at cell tops for 87 days to reach the most homogeneous saturation of bentonite possible (Figure 1b). The duration of 28 days for the first saturation stage was chosen based on a preceding study involving 15 different bentonites [20,21]. According to the experimental protocol established in that study, the second saturation stage was continued until the swelling pressure of the specimens was judged to level off. The fluid pressure difference across the specimens was equal to 0.018 MPa during this stage.
After disconnecting the burettes, permeation tests were started on the 133rd day (after a delay due to a malfunction and temporary unavailability of the pump) by injecting SGW through inlets at cell bases using a piston pump (BESTA HD 2-200, BESTA-Technik GmbH, Wilhelmsfeld, Germany) and a bypass (Figure 1b) to prevent uncontrolled fluid injection pressure surges observed in a preceding study discussed in [20,21]. Outlets at cell tops were connected to burettes for measuring volumes of percolating SGW used to calculate permeability. Swelling of compacted bentonites upon saturation was assumed to seal off possible flow pathways along the bentonite-cell interface with the solution percolating only through bentonite pores during permeability measurements, which was confirmed by the observed low permeabilities of specimens. The fluid pressure difference of 0.018 MPa maintained across the specimens for 132 days by the water column in the burettes and was then incrementally increased to 0.2, 0.5, 2.0, and 8.0 MPa using the pump. These fluid pressure differences correspond to hydraulic gradients of 180 up to 80,000 m/m across the compacted bentonite specimens. For the first stage of the permeation tests at a hydraulic gradient of 2000 m/m, levelling off of permeabilities while observing for the levelling off of swelling pressures, if reasonably possible, was used as a termination criterion. During this stage, the pore volume of specimens was replaced on average 13.7, 2.2, and 0.5 times at a dry density of 1.4, 1.6, and 1.8 g/cm3, respectively. For the following stages, a sufficient number of measurement points to capture the permeability trend was aimed, and the pore volume of specimens was replaced at a respective dry density of 1.4, 1.6, and 1.8 g/cm3 on average 0.9, 0.1, and 0.02 times at a hydraulic gradient of 5000 m/m up to 6.9, 4.0, and 1.8 times at a hydraulic gradient of 80,000 m/m. Room temperature and relative humidity during the entire experiment were 18.1 ± 0.4 °C and 48% ± 7%, respectively. A dynamic viscosity of 1.014 ± 0.009 µPa·s was determined for SGW from eight replicate measurements at room temperature.
Apart from the specimens for the swelling and permeation experiments described above, but from the same batch of 10X, further specimens were compacted at hygroscopic humidity for studies using the method of computed microtomography (μCT, using Yamato TDM1000H-2) and qualitative analyses of the microstructure by scanning electron microscopy (SEM) using LEO 1450 VP (Zeiss, Jena, Germany) with a guaranteed resolution of 5 nm). These specimens were prepared in the same manner as the swelling and permeation experiments with distilled water described in [4]. To study changes in the microstructure of 10X during swelling, SEM analyses were conducted with specimens after permeation with distilled water [4].

3. Results and Discussion

The µCT study of the structure of compacted 10X showed that the specimen at a density of 1.4 g/cm3 is characterized by an uneven distribution of moisture (Figure 2), and areas of non-compacted material are no more than 0.25 mm in size. The specimen at a density of 1.6 g/cm3 is characterized by greater homogeneity, although a small amount of a non-compacted material (up to 0.1 mm on average, some areas with 0.5 mm) and some moisture inhomogeneity are observed. When compacted to a density of 1.8 g/cm3, the obtained specimen does not contain any non-compacted parts, and no inhomogeneous distribution of moisture between the formed aggregates was observed. The appearance of cracks at the edges of the specimen at a density of 1.8 g/cm3 (Figure 2c) was attributed to a relaxation of the 10X’s fabric upon taking the specimen out of the compaction mold for the µCT study.
The initial stage of bentonite swelling upon saturation with solutions of low salinity [22,23] and—for some bentonites—with brines [20,21] is characterized by a “double-peak” pattern as a result of (i) swelling of aggregates (the first peak), (ii) their subsequent partial decomposition and a collapse of the inter-aggregate pores (the depression), and (iii) swelling of the reorganised aggregates (the second peak). This pattern is also observed for 10X specimens at the studied densities (graphs on left-side in Figure 3). The variation of this pattern with changing density further exemplifies that the collapse of the inter-aggregate pores is more strongly pronounced at the higher porosity so that the second peak is much lower than the first peak at 1.4 g/cm3. This effect diminishes with decreasing porosity so that the second peak becomes higher than the first peak within 1 to 4 days at 1.8 g/cm3. The higher second peak was observed previously for low- and high-salinity solutions ([20,21,22,23,24,25] and references therein), whereas the lower second peak was observed for high-salinity solutions [20,21]. By combining the latter observation with the present one, it can be suggested that the collapse of the inter-aggregate pores may depend on the electrostatic repulsion between diffuse-double layers (DDL) of single montmorillonite particles or their aggregates, as high salinity is known to depress their DDL [26] and increase pore volume available for their reorganisation accordingly.
It can be further seen in Figure 3 that despite the identical preparation protocol of specimens, the method used did not eliminate the inherent inhomogeneity of the microstructure of compacted bentonites. As a result of this inhomogeneity, measured swelling pressures and the time span of development for the second peak of the swelling pattern show significant variation at a given density and are otherwise similar to initial swelling patterns. This variation would not disappear (and even became larger at the densities of 1.4 and 1.6 g/cm3) upon switching the solution supply from cell bases to cell tops on the 28th day. Possible reasons for this behaviour can be better understood considering the solution volumes absorbed by specimens and channels of permeameter cells from burettes within the first 115 days of the experiment (Figure 4). Firstly, these data reveal that the bentonite saturation was not accomplished within this stage, as the solution uptake did not level off. This occurrence is a well-known drawback of permeameter tests performed without a back pressure high enough to ensure full saturation of compacted specimens [27]. Secondly, for specimens with estimated pore volumes of 7.0–9.8 cm3 at a given density, a variation in the absorbed solution volumes of up to 10.7 cm3 was observed, which clarifies that the air was still present in the channels of some permeameters or even in the specimens. Hence, a reliable determination of the solution uptake by specimens upon saturation from burettes in the experimental setup used seems only possible through a test termination, dismantling of permeameters, and weighing the specimens.
The subsequent application of the hydraulic gradient of 2000 m/m led to a significant gradual increase of swelling pressure in all specimens, presumably, as a result of displacement of residual air from specimens by the injected solution (as discussed below to permeabilities) and/or microstructural reorganisation (as discussed below). During this stage and the preceding saturation stage using burettes, the total pressure measured at cell tops was equal to the swelling pressure of specimens. Starting from the 246th day, however, it was equal to the sum of the swelling pressure of the specimens and the fluid pressure established at cell tops upon imposing the hydraulic gradient of ≥5000 m/m across the specimens. This can be seen in steep increases of the total pressure at top of all cells on the 246th, 261st (Figure 3), and 281st day (Figure 5) upon increasing the hydraulic gradient to 5000; 20,000, and 80,000 m/m, respectively. It also results from its steep decrease to the values measured before the 246th day upon the intermediate and final switching off the pump on the 294th and 297th day, respectively (Figure 5). Therefore, the swelling pressure of 10X compacted to the densities of 1.4, 1.6, and 1.8 g/cm3 can be concluded to equal to 0.8 ± 0.3, 2.2 ± 0.6, and 6.3 ± 0.3 MPa, respectively, measured at the hydraulic gradient of 2000 m/m (Table 2).
Based on these values, the fluid pressure at cell tops can be estimated to equal 1.3, 0.9, and 0.6 MPa for the fluid pressure of 2.0 MPa at cell base (hydraulic gradient of 20,000 m/m) and 6.8, 6.1, and 3.4 MPa for the fluid pressure of 8.0 MPa at cell base (hydraulic gradient of 80,000 m/m) at dry densities of 1.4, 1.6, and 1.8 g/cm3, respectively. The decrease of fluid pressure from cell bases to cell tops is in line with previous observations of the non-linear variation of the excess fluid pressure across specimens in constant-head experiments with compacted clays [28]. A decrease in these values with increasing density further indicates that the degree of non-linearity in this variation increases with the decreasing porosity of the specimens (0.50, 0.43, and 0.36, respectively).
In the preliminary 7-day experiments [4] with 10X and distilled water using permeameters similar to those of the present study, for dry densities in the range of 1.44–1.82 g/cm3, swelling pressures in the range of 0.52–9.39 MPa and a relationship of the form P = exp(6.4ρd − 9.403) between swelling pressure P (MPa) and dry density ρd (g/cm3) were obtained. At the hydraulic gradient of 2000 m/m in the present study, a somewhat different relationship P = exp(5.16ρd − 7.445) was obtained (Figure 6). The difference might be due to the use of only three data points for obtaining the relationship in the present study, or the difference between SGW and distilled water, or for other reasons, which cannot be identified in view of the limited information on parameters of the experiments and statistical weight of fitted data in the study [4]. In the study by Wang et al. [29], relationships of the form P = exp(6.75ρd − 8.634), P = exp(6.85ρd − 9.231), and P = exp(3.32ρd − 5.608) were provided for MX-80 (75–90% smectite), FEBEX (92% smectite), and Kunigel V1 (48% smectite) bentonites, respectively. Based on the smectite content of 73 mass% of 10X [4], it may be speculated that the value of the coefficient preceding the density ought to be in a range from 3.3 to 6.8. Although the values obtained in [4] and in the present study agree with this range, establishing values that are characteristic of 10X would require additional experimental data.
Permeabilities measured for 10X in the present study show a gradual decrease at all densities during the initial experimental stage at the hydraulic gradient of 2000 m/m (Figure 7). The gradual decrease of permeabilities suggests that the residual air, if any, was largely displaced from the specimens at the very beginning of this stage since air removal would otherwise lead to a non-transient increase of permeability [27], which was not observed later on. At densities of 1.4 and 1.6 g/cm3, a double linear decrease was observed with a steeper initial decrease during the first 10 days at 1.4 g/cm3 and 40 days at 1.6 g/cm3 followed by a flatter one. Shortly after further increasing the hydraulic gradients at these densities, considerably higher transient permeabilities were measured (data points on the vertical dotted lines in Figure 7, which should not be treated as outliers but rather attributed to the additional amount of water squeezed out of specimens by seepage-induced (local) consolidation under the pressure gradient [28,30]. According to [30], the initial transient response preceding the steady-state flow condition can take up to a few hours for clay specimens due to time-dependent changes in the volume or distribution of pore space in a specimen.
As a result of hydraulic gradient increases above 2000 m/m, steady-state permeability showed two major decreases by a factor of ~2–2.5 at a density of 1.4 g/cm3 and one major decrease by a factor of ~1.5 at 1.6 g/cm3, whereas it stabilized at ~6 × 10−21 m2 at 1.8 g/cm3 (Figure 7, Table 2). The seepage-induced consolidation, which reduces the specimen’s pore volume, was proposed as a reason for permeability decrease with increasing pressure gradient [28,30]. This reason cannot be discarded for the present experiment, as the observed dependence of the extent of permeability decrease on the density would be in line with a larger seepage-induced consolidation for less dense specimens. If this proposal is valid, a larger hydraulic gradient than those applied in the present study can be expected to cause a major decrease of steady-state permeability of 10X specimen at 1.8 g/cm3 as well.
Alternatively, the hydraulic conductivity of expansive materials was observed to decrease after full saturation due to reorganisation of microstructure in the course of a water potential re-equilibration and corresponding water redistribution between the inter-aggregate and intra-aggregate pores [31]. Accordingly, it can be assumed that such microstructural reorganisation may have occurred to a larger extent for 10X specimen at a density of 1.4 g/cm3 and to a lesser extent at 1.6 g/cm3, whereas it did not take place at 1.8 g/cm3, at least for hydraulic gradients ≥5000 m/m. This assumption would be in line with the decreasing pore volume available for a microstructural reorganisation with increasing density.
A change in the pore space during the SGW saturation was detected by scanning electron microscopy. The microstructure of the original compacted 10X is characterized by the predominance of fine interparticle micropores and large inter-micro-aggregate pores. During saturation with the SGW, the number and size of pores are decreased (Figure 8), inter-micro-aggregate and small inter-micro-aggregate micropores predominate, and an increase in the size and isometricity of large inter-micro-aggregate pores is also noted. There is also a decrease in the degree of orientation of structural elements (clay aggregates of different sizes) in the plane orthogonal to the pressing direction, acquired during the process of compacting the specimen, in which the microstructure becomes more uniform.
This uniformity in microstructure is explained by the processes of hydration of highly hydrophilic clay particles upon saturation of the specimens and accompanying swelling of microaggregates, closure of fine pores, and opening of inter-ultra-micro-aggregate (inter-micro-aggregate) micropores. As a result, the microstructure of the specimens after permeation becomes more homogeneous, intra- and inter-micro-aggregate pores become predominant since they open significantly during the swelling stage but are not affected by the fluid pressure during the permeation stage. This effect can also be associated with the swelling of the specimens during saturation and with the so-called “hinge-joint effect” of clay particles in a micro-aggregate during hydration [32].
The swelling of clays is often characterized by the effect of disordering the microstructure. According to the results of experiments and calculations given in the article [33], this effect is associated with the distribution of forces acting on particles and microaggregates, the result of which leads to the appearance of the hinge-joint effect. Thus, as a result of the wedging pressure of the hydration shells around the particles, their mutual repulsion occurs, while the forces of attraction at their edges prevent complete separation, resulting in the formation of a honeycomb structure with higher porosity and a lower degree of microstructure uniformity (Figure 8, sidebars).
A displacement of particles, which are not directly involved in the load-carrying of specimen’s fabric downward the flow direction, resulting in clogging of specimen’s pores was proposed as a further possible reason for seepage-induced changes in the permeability [30]. In the study by Al-Taie et al. [34], clogging of sintered filters by this mechanism was proposed as an alternative explanation for permeability decreases observed upon increases of hydraulic gradients above 1000 m/m in experiments with a smectite-rich clay and distilled water. However, in the infiltration tests of hydraulic gradients up to ~64,000 m/m [29], no density changes along the profile of a bentonite specimen, which would provide evidence of such displacement, were found. Furthermore, the initial permeability of compacted 10X was ~10−18 m2 at 1.4 g/cm3 (Figure 7), whereas that of sintered filters was (3.4 ± 0.5) × 10−14 m2, so that clogging of pore space would be more probable than that of sintered filters at 1.4 g/cm3 and even more so at 1.8 g/cm3. Since, no major decreases of steady-state permeability occurred at this density, clogging of pore space or sintered filters due to displacement of particles within the studied specimens was not a factor contributing significantly to the observed permeability decreases. Indeed, the permeability of uncleaned sintered filters after the experiment was measured to equal (1.9 ± 1.2) × 10−14, (1.5 ± 0.6) × 10−14, and (1.2 ± 0.2) × 10−14 m2 at a dry density of 1.4, 1.6, and 1.8 g/cm3, respectively, or a factor of 2 to 3 lower than that before the experiment. Thus, though some clogging of sintered filters occurred during the experiment, its contribution to the observed permeability decreases can be neglected.
In view of remaining uncertainty about the mechanism behind the permeability decreases observed at higher hydraulic gradients, permeability values of (27 ± 15) × 10−20, (3.4 ± 0.8) × 10−20, and (0.96 ± 0.26) × 10−20 m2 (Table 2) obtained at a hydraulic gradient of 2000 m/m were suggested as conservative estimates of permeabilities of 10X for SGW at a temperature of 18.1 ± 0.4 °C and dry densities of 1.4, 1.6, and 1.8 g/cm3, respectively. When assessing the performance of the bentonite in the early post-closure phase of DGR at hydraulic gradients higher than 2000 m/m, corresponding permeability estimates given in Table 2 may be used.
In the experiments with seawater, permeability κ (expressed in millidarcies, 9.869233 × 10−19 m2) of pure smectite and kaolinite was found to be related to bulk density ρbulk (g/cm3) by respective equations lnρbulk = −0.037lnκ + 0.27 and lnρbulk = −0.074lnκ + 0.40 [35]. Notably, as can be seen from the equations derived in the present study (Figure 9), the coefficient preceding ln κ at hydraulic gradients of 2000 and 5000 m/m (0.040 ± 0.005 and 0.035 ± 0.001, respectively) is close to the value obtained in study [35] for pure smectite, whereas at 80,000 m/m it is close to that for a mixture of 20% smectite and 80% kaolinite.
The equation ln ρbulk = −0.040lnκ(mD) + 0.47 (Figure 9), obtained in the present study at the hydraulic gradient of 2000 m/m, corresponds to ln κ(m2) = −13.4lnρdry − 38.3 (Figure 6), with κ being a power function of ρd. In this form it can be compared with the relationship ln κ(m2) = −9.03ρd − 14.778 obtained in the experiments with distilled water [4] where κ is an exponential function of ρd. In the semi-logarithmic plot (Figure 6), the latter relationship has the linear form whereas that used in the present study and the study by Mondol et al. [35] does not. In view of this difference, it appears that establishing the form of the relationship between permeability and density that is characteristic of 10X requires additional experimental data.
The comparison of the fitted curves in Figure 6 further reveals that permeabilities obtained in [4] at hydraulic gradients of 22,000 up to 92,000 m/m are about 2.5 times smaller than those in the present study obtained at a hydraulic gradient of 2000 m/m. Upon increasing the hydraulic gradient to 20,000 m/m in the present study, the permeability decreased by a factor of 2.5 at 1.4 g/cm3 and of 1.5 at 1.8 g/cm3 (Table 2), whereas at 1.6 g/cm3 it decreased by a factor of 1.7 only upon an increase of hydraulic gradient to 80,000 m/m. Therefore, the difference between the permeabilities measured in the present study and in [4] may be due to hydraulic gradients higher than 2000 m/m in the latter study. A contribution of other reasons to this difference cannot be excluded and needs to be verified in additional experiments.

4. Conclusions

Differences in the swelling pressures and permeabilities observed in the present study for 10X bentonite at given dry densities suggests that because of inherent inhomogeneity in the microstructure of compacted bentonites, at least duplicate, preferably triplicate, experiments should be conducted when using the present method to estimate swelling pressures or permeabilities. In this regard, a saturation of compacted bentonite specimens and measurements of their swelling pressure and permeability are recommended to be conducted at a hydraulic gradient of 2000 m/m. As higher hydraulic gradients cause modifications of compacted bentonite microstructure, which are still not well understood, the estimates obtained with those should be used with care, and their applicability may be restricted to the early post-closure phase of DGR.
To obtain a reliable estimate of the dependence of swelling pressure and permeability of compacted bentonite on its density, a series of triplicate experiments for at least 7–8 densities, preferably with an even distribution over a broader density range than that of 1.4–1.8 g/cm3 studied in work [4] and the present study, is needed. The dataset for 10X bentonite should be accordingly extended.

Author Contributions

Conceptualization, data curation, investigation, methodology, visualization, supervision, and writing—review and editing, A.Y.M.; investigation, methodology, visualization, S.V.Z. and V.V.K. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by the Federal Ministry of Economic Affairs and Energy (BMWi), represented by the Project Management Agency Karlsruhe (PTKA-WTE) within the frames of the project, “Sicherheitsanalytische Untersuchungen zu Endlagersystemen im Kristallin” (SUSE) with contract no. 02E11577B. Initial data, as well as some additional properties’ tests of bentonites and microstructure of compacted specimens, were supported by the Russian Science Foundation (project 16-17-10270).

Data Availability Statement

Not Applicable.

Acknowledgments

Swelling pressure and permeability results were obtained owing to the skillful lab work of Julia Gose of the geoscientific laboratory of the GRS. The authors thank Michail Chernov and Ruslan Kuznetsov (MSU, Geological Faculty) for providing data by the SEM and μCT. Experimental studies were partially performed using equipment acquired with the funding of the Moscow State University Development Program (X-ray Diffractometer Ultima-IV, Rigaku, and Scanning Electron Microscope LEO 1450VP, Carl Zeiss, X-Ray tomograph Yamato TDM1000, Yamato Scientific). We would like to thank two anonymous reviewers and Marvin Middlehoff (GRS) for constructive comments that helped to improve the original manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abramov, A.A.; Bolshov, L.A.; Dorofeeev, A.N.; Igin, I.M.; Kazakov, K.S.; Krasilnikov, V.Y.; Linge, I.I.; Trokhov, N.N.; Utkin, S.S. Underground research laboratory in the Nizhnekanskiy massif: Evolutionary design study. Radioact. Waste 2020, 1, 9–21. (In Russian) [Google Scholar] [CrossRef]
  2. Dorofeev, A.N.; Bolshov, L.A.; Linge, I.I.; Utkin, S.S.; Saveleva, E.A. Strategic master-plan for research demonstrating the safety of construction, operation and closure of a deep geological repository for radioactive waste. Radioact. Waste 2017, 1, 19–26. [Google Scholar]
  3. Krupskaya, V.V.; Biryukov, D.V.; Belousov, P.E.; Lekhov, V.A.; Romanchuk, A.Y.; Kalmykov, S.N. The use of natural clay materials to increase the nuclear and radiation safety level of nuclear legacy facilities. Radioact. Waste 2018, 2, 30–43. [Google Scholar]
  4. Krupskaya, V.V.; Zakusin, S.V.; Lekhov, V.A.; Dorzhieva, O.V.; Belousov, P.E.; Tyupina, E.A. Buffer properties of bentonite barrier systems for radioactive waste isolation in geological repository in the Nizhnekanskiy Massif. Radioact. Waste 2020, 1, 35–55. [Google Scholar] [CrossRef]
  5. Birgersson, M.; Hedström, M.; Karnland, O.; Sjöland, A. Bentonite buffer: Macroscopic performance from nanoscale properties. In Geological Repository Systems for Safe Disposal of Spent Nuclear Fuels and Radioactive Waste, 2nd ed.; Michael, J.A., Joonhong, A., Eds.; Woodhead Publishing: Duxford, UK, 2017; Volume 12, pp. 319–364. [Google Scholar] [CrossRef]
  6. Kaufhold, S.; Dohrmann, R. Distinguishing between more and less suitable bentonites for storage of high-level radioactive waste. Clay Miner. 2016, 51, 289–302. [Google Scholar] [CrossRef] [Green Version]
  7. Sellin, P.; Leupin, O.X. The use of clay as an engineered barrier in radioactive-waste management—A review. Clays Clay Miner. 2013, 61, 477–498. [Google Scholar] [CrossRef]
  8. Belousov, P.E.; Krupskaya, V.V. Bentonite clays of Russia and neighboring countries. Georesursy 2019, 21, 79–90. [Google Scholar] [CrossRef]
  9. Belousov, P.E.; Krupskaya, V.V.; Zakusin, S.V.; Zhigarev, V.V. Bentonite clays from 10th Khutor deposite: Features of genesis, composition and adsorption properties. RUDN J. Eng. Res. 2017, 18, 135–143. [Google Scholar] [CrossRef] [Green Version]
  10. Belousov, P.; Chupalenkov, N.; Zakusin, S.; Morozov, I.; Dorzhieva, O.; Chernov, M.; Krupskaya, V. 10th Khutor bentonite deposit (Russia): Geology, mineralogy and industrial application. Appl. Clay Sci. 2021, in press. [Google Scholar]
  11. Belousov, P.E.; Pokidko, B.V.; Zakusin, S.V.; Krupskaya, V.V. Quantitative methods for quantification of montmorillonite content in bentonite clays. Georesursy 2020, 22, 38–47. [Google Scholar] [CrossRef]
  12. Post, J.E.; Bish, D.L. Rietveld refinement of crystal structures using powder X-ray diffraction data. Rev. Miner. Geochem. 1989, 20, 277–308. [Google Scholar]
  13. Doebelin, N.; Kleeberg, R. Profex: A graphical user interface for the Rietveld refinement program BGMN. J. Appl. Crystallogr. 2015, 48, 1573–1580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. GOST 5180-2015 Soils. Laboratory Methods for Determination of Physical Characteristics; Standartinform: Moscow, Russia, 2019. [Google Scholar]
  15. ASTM D854-14. Standard Test Methods for Specific Gravity of Soil Solids by Water Pycnometer; ASTM International: West Conshohocken, PA, USA, 2014. [Google Scholar]
  16. Lorenz, P.; Meier, L.; Kahr, G. Determination of the cation exchange capacity (CEC) of clays minerals using the complexes of copper (II) ion with triethylenetetramine and tetraethylenepentamine. Clays Clay Miner. 1999, 47, 386–388. [Google Scholar]
  17. Dohrmann, R.; Genske, D.; Karnland, O.; Kaufhold, S.; Kiviranta, L.; Olsson, S.; Plötze, M.; Sandén, T.; Sellin, P.; Svensson, D.; et al. Interlaboratory CEC and exchangeable cation study of bentonite buffer materials: I. Cu(II)-triethylenetetramine method. Clays Clay Miner. 2012, 60, 162–175. [Google Scholar] [CrossRef]
  18. Stanjek, H.; Künkel, D. CEC determination with Cu-triethylenetetramine: Recommendations for improving reproducibility and accuracy. Clay Miner. 2016, 51, 1–17. [Google Scholar] [CrossRef]
  19. Rozov, K.B.; Rumynin, V.G.; Nikulenkov, A.M.; Leskova, P.G. Sorption of 137Cs, 90Sr, Se, 99Tc, 152(154)Eu, 239(240)Pu on fractured rocks of the Yeniseysky site (Nizhne-Kansky massif, Krasnoyarsk region, Russia). J. Environ. Radioact. 2018, 192, 513–523. [Google Scholar] [CrossRef]
  20. Meleshyn, A.Y. Swelling behaviour and permeability of compacted bentonites at high salinity as influenced by fluid pressure changes. Radioact. Waste 2020, 2, 109–119, (In Russian; the English version is in press). [Google Scholar] [CrossRef]
  21. Meleshyn, A. Mechanisms of Transformation of Bentonite Barriers (UMB); Report GRS-390 BMWi-FKZ 02 E 11344A.; GRS: Braunschweig, Germany, 2019; Available online: https://www.grs.de/publikationen/grs-390 (accessed on 30 June 2021).
  22. Yigzaw, Z.G.; Cuisinier, O.; Massat, L.; Masrouri, F. Role of different suction components on swelling behavior of compacted bentonites. Appl. Clay Sci. 2016, 120, 81–90. [Google Scholar] [CrossRef]
  23. Massat, L.; Cuisinier, O.; Bihannic, I.; Claret, F.; Pelletier, M.; Masrouri, F.; Gaboreau, S. Swelling pressure development and inter-aggregate porosity evolution upon hydration of a compacted swelling clay. Appl. Clay Sci. 2016, 124, 197–210. [Google Scholar] [CrossRef]
  24. Zhu, C.M.; Ye, W.M.; Chen, Y.G.; Chen, B.; Cui, Y.J. Influence of salt solutions on the swelling pressure and hydraulic conductivity of compacted GMZ01 bentonite. Eng. Geol. 2013, 166, 74–80. [Google Scholar] [CrossRef]
  25. Ye, W.M.; Wan, M.; Chen, B.; Chen, Y.G.; Cui, Y.J.; Wang, J. Temperature effects on the swelling pressure and saturated hydraulic conductivity of the compacted GMZ01 bentonite. Environ. Earth Sci. 2013, 68, 281–288. [Google Scholar] [CrossRef]
  26. Zhang, F.; Ye, W.M.; Wang, Q.; Chen, Y.G.; Chen, B. An insight into the swelling pressure of GMZ01 bentonite with consideration of salt solution effects. Eng. Geol. 2019, 251, 190–196. [Google Scholar] [CrossRef]
  27. Chapuis, R.P. Sand–bentonite liners: Predicting permeability from laboratory tests. Can. Geotech. J. 1990, 27, 47–57. [Google Scholar] [CrossRef]
  28. Pane, V.; Croce, P.; Znidarcic, D.; Ko, H.Y.; Olsen, H.W.; Schiffman, R.L. Effects of consolidation on permeability measurements for soft clay. Geotechnique 1983, 33, 67–72. [Google Scholar] [CrossRef]
  29. Wang, Q.; Tang, A.M.; Cui, Y.J.; Delage, P.; Gatmiri, B. Experimental study on the swelling behaviour of bentonite/claystone mixture. Eng. Geol. 2012, 124, 59–66. [Google Scholar] [CrossRef] [Green Version]
  30. Olsen, H.W.; Nichols, R.W.; Rice, T.L. Low gradient permeability measurements in a triaxial system. Geotechnique 1985, 35, 145–157. [Google Scholar] [CrossRef]
  31. Wang, Q.; Cui, Y.J.; Tang, A.M.; Barnichon, J.D.; Saba, S.; Ye, W.M. Hydraulic conductivity and microstructure changes of compacted bentonite/sand mixture during hydration. Eng. Geol. 2013, 164, 67–76. [Google Scholar] [CrossRef]
  32. Osipov, V.I.; Sokolov, V.N. Clays and Its Properties. In Composition, Structure and Properties Formation; GEOS: Moscow, Russia, 2013; p. 576. ISBN 978-5-89118-617-0. [Google Scholar]
  33. Osipov, V.I.; Babak, V.G. The nature and mechanism of clay swelling. Eng. Geol. 1987, 5, 18–27. [Google Scholar]
  34. Al-Taie, L.; Pusch, R.; Knutsson, S. Hydraulic properties of smectite rich clay controlled by hydraulic gradients and filter types. Appl. Clay Sci. 2014, 87, 73–80. [Google Scholar] [CrossRef]
  35. Mondol, N.H.; Bjørlykke, K.; Jahren, J. Experimental compaction of clays: Relationship between permeability and petrophysical properties in mudstones. Pet. Geosci. 2008, 14, 319–337. [Google Scholar] [CrossRef]
Figure 1. A construction drawing of experiments; (a); supply of SGW to inlets at cell bases by a piston pump; (b) permeability experiments (note that burettes connected to outlets at cell tops to collect percolating solutions).
Figure 1. A construction drawing of experiments; (a); supply of SGW to inlets at cell bases by a piston pump; (b) permeability experiments (note that burettes connected to outlets at cell tops to collect percolating solutions).
Minerals 11 00742 g001
Figure 2. μCT images of 10X specimens compacted to dry densities of (a) 1.4 g/cm3; (b) 1.6 g/cm3; (c) 1.8 g/cm3.
Figure 2. μCT images of 10X specimens compacted to dry densities of (a) 1.4 g/cm3; (b) 1.6 g/cm3; (c) 1.8 g/cm3.
Minerals 11 00742 g002
Figure 3. Total pressure at tops of permeameter cells for 10X specimens compacted to dry densities of 1.4 g/cm3 (cells A1–A3); 1.6 g/cm3 (cells B1–B3); and 1.8 g/cm3 (cells C2, C3), saturated and permeated by SGW. Graphs on left-side represent x- and y-axis close-ups of the data measured within the first four days of saturation. Vertical dotted lines mark onsets of solution supply by burettes through cell tops (28th day) and by the pump at hydraulic gradients of 2000 m/m (133rd day), 5000 m/m (246th day), and 20,000 m/m (261st day) (note that corresponding labels in the 2nd graph on right-side).
Figure 3. Total pressure at tops of permeameter cells for 10X specimens compacted to dry densities of 1.4 g/cm3 (cells A1–A3); 1.6 g/cm3 (cells B1–B3); and 1.8 g/cm3 (cells C2, C3), saturated and permeated by SGW. Graphs on left-side represent x- and y-axis close-ups of the data measured within the first four days of saturation. Vertical dotted lines mark onsets of solution supply by burettes through cell tops (28th day) and by the pump at hydraulic gradients of 2000 m/m (133rd day), 5000 m/m (246th day), and 20,000 m/m (261st day) (note that corresponding labels in the 2nd graph on right-side).
Minerals 11 00742 g003
Figure 4. Uptake of SGW by permeameter cells with 10X specimens compacted to dry densities of 1.4 g/cm3 (A1–A3); 1.6 g/cm3 (B1–B3); and 1.8 g/cm3 (C2, C3), upon saturation from burettes connected to inlets at cell bases for 27 days and cell tops afterwards.
Figure 4. Uptake of SGW by permeameter cells with 10X specimens compacted to dry densities of 1.4 g/cm3 (A1–A3); 1.6 g/cm3 (B1–B3); and 1.8 g/cm3 (C2, C3), upon saturation from burettes connected to inlets at cell bases for 27 days and cell tops afterwards.
Minerals 11 00742 g004
Figure 5. Total pressure at tops of permeameter cells for 10X specimens compacted to dry densities of 1.4 g/cm3 (A1–A3), 1.6 g/cm3 (B1–B3), and 1.8 g/cm3 (C2, C3) and permeated by SGW at a hydraulic gradient of 80,000 m/m. Real fluid injection pressure imposed by the pump is shown in the upper graph (note that the abrupt fluid injection pressure changes were caused by short switching off the pump before weekends to disconnect cells A1–A3, preventing an overflow of burettes connected to them and after weekends, reconnecting them).
Figure 5. Total pressure at tops of permeameter cells for 10X specimens compacted to dry densities of 1.4 g/cm3 (A1–A3), 1.6 g/cm3 (B1–B3), and 1.8 g/cm3 (C2, C3) and permeated by SGW at a hydraulic gradient of 80,000 m/m. Real fluid injection pressure imposed by the pump is shown in the upper graph (note that the abrupt fluid injection pressure changes were caused by short switching off the pump before weekends to disconnect cells A1–A3, preventing an overflow of burettes connected to them and after weekends, reconnecting them).
Minerals 11 00742 g005
Figure 6. Swelling pressure and permeability (m2) of 10X specimens permeated by SGW as a function of dry density ρd. Circles and solid lines represent measured values and fitted curves of the present study, whereas dashed lines represent fitted curves of the study [4] (hydraulic conductivities reported in [4] were converted to permeabilities for the purpose of comparison).
Figure 6. Swelling pressure and permeability (m2) of 10X specimens permeated by SGW as a function of dry density ρd. Circles and solid lines represent measured values and fitted curves of the present study, whereas dashed lines represent fitted curves of the study [4] (hydraulic conductivities reported in [4] were converted to permeabilities for the purpose of comparison).
Minerals 11 00742 g006
Figure 7. Permeability (m2) of 10X specimens compacted to dry densities of 1.4 g/cm3 (A1–A3); 1.6 g/cm3 (B1–B3); 1.8 g/cm3 (C2, C3), and permeated by SGW. Vertical dotted lines mark increases of hydraulic gradient to 5000 m/m (246th day); 20,000 m/m (261st day), and 80,000 m/m (281st day).
Figure 7. Permeability (m2) of 10X specimens compacted to dry densities of 1.4 g/cm3 (A1–A3); 1.6 g/cm3 (B1–B3); 1.8 g/cm3 (C2, C3), and permeated by SGW. Vertical dotted lines mark increases of hydraulic gradient to 5000 m/m (246th day); 20,000 m/m (261st day), and 80,000 m/m (281st day).
Minerals 11 00742 g007
Figure 8. SEM images of the compacted 10X, fraction <0.25 mm, density 1.65 g/cm3, at 32,000×: (a) natural micro-aggregates of clay particles after an air-drying; (b) partially swelled micro-aggregate of clay particles. Sidebars on the top corners: structural model of clay swelling with different moisture (W) content: (a) W ≤ Wmg; (b) Wp < W < Wsw (Wmg—maximal hygroscopic moisture content; Wp—moisture of the plastic limit; Wsw—moisture of swelling). 1—clay particle; 2—adsorbed water; 3—osmotic water.
Figure 8. SEM images of the compacted 10X, fraction <0.25 mm, density 1.65 g/cm3, at 32,000×: (a) natural micro-aggregates of clay particles after an air-drying; (b) partially swelled micro-aggregate of clay particles. Sidebars on the top corners: structural model of clay swelling with different moisture (W) content: (a) W ≤ Wmg; (b) Wp < W < Wsw (Wmg—maximal hygroscopic moisture content; Wp—moisture of the plastic limit; Wsw—moisture of swelling). 1—clay particle; 2—adsorbed water; 3—osmotic water.
Minerals 11 00742 g008
Figure 9. Permeability (mD) of 10X permeated by SGW at hydraulic gradients of 2000–80,000 m/m as a function of bulk density. Circles and solid lines represent measured values and fitted curves.
Figure 9. Permeability (mD) of 10X permeated by SGW at hydraulic gradients of 2000–80,000 m/m as a function of bulk density. Circles and solid lines represent measured values and fitted curves.
Minerals 11 00742 g009
Table 1. Basic characteristics of bentonite powder 10X.
Table 1. Basic characteristics of bentonite powder 10X.
a—Particle Size Distribution
Mods, µmMedian, µmParticle Size Distribution, %
0.1–1 µm1–10 µm10–25 µm
1.792.9015.672.611.8
b—Mineral Composition (%)
Smectite
(Montmorillonite)
ChloriteQuartzFeldspars
(Microcline and Albite)
CalciteAnataseGypsum
73.00.913.910.80.70.30.4
c—Chemical Composition
Na2OMgOAl2O3SiO2K2OCaOTiO2Fe2O3FeOMnOLOI (105–1000 °C)
0.1072.6718.8451.870.486.580.703.370.330.017.64
Table 2. Total pressures at the top (σ) and permeabilities (κ) averaged over the last two weeks of measurement at a given hydraulic gradient (iH) of 10X specimens compacted to dry densities of 1.4, 1.6, and 1.8 g/cm3; saturated and permeated by SGW.
Table 2. Total pressures at the top (σ) and permeabilities (κ) averaged over the last two weeks of measurement at a given hydraulic gradient (iH) of 10X specimens compacted to dry densities of 1.4, 1.6, and 1.8 g/cm3; saturated and permeated by SGW.
iH, m/mσ, MPaκ, 10−20 m2
1.4 g/cm31.6 g/cm31.8 g/cm31.4 g/cm31.6 g/cm31.8 g/cm3
1800.7 ± 0.21.5 ± 0.65.3 ± 0.9---
20000.8 ± 0.32.2 ± 0.66.3 ± 0.327 ± 153.4 ± 0.80.96 ± 0.26
50000.9 ± 0.32.3 ± 0.66.4 ± 0.325 ± 143.2 ± 0.40.54 ± 0.10
20,0002.1 ± 0.23.1 ± 0.66.9 ± 0.311 ± 42.9 ± 0.40.68 ± 0.02
80,0007.6 ± 0.18.3 ± 0.89.7 ± 0.35.4 ± 1.52.0 ± 0.20.67 ± 0.07
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Meleshyn, A.Y.; Zakusin, S.V.; Krupskaya, V.V. Swelling Pressure and Permeability of Compacted Bentonite from 10th Khutor Deposit (Russia). Minerals 2021, 11, 742. https://doi.org/10.3390/min11070742

AMA Style

Meleshyn AY, Zakusin SV, Krupskaya VV. Swelling Pressure and Permeability of Compacted Bentonite from 10th Khutor Deposit (Russia). Minerals. 2021; 11(7):742. https://doi.org/10.3390/min11070742

Chicago/Turabian Style

Meleshyn, Artur Yu., Sergey V. Zakusin, and Victoria V. Krupskaya. 2021. "Swelling Pressure and Permeability of Compacted Bentonite from 10th Khutor Deposit (Russia)" Minerals 11, no. 7: 742. https://doi.org/10.3390/min11070742

APA Style

Meleshyn, A. Y., Zakusin, S. V., & Krupskaya, V. V. (2021). Swelling Pressure and Permeability of Compacted Bentonite from 10th Khutor Deposit (Russia). Minerals, 11(7), 742. https://doi.org/10.3390/min11070742

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