Geological Repository Research Articles

Undrained triaxial compression tests on normally consolidated bentonite considering temperature/confining pressure dependency

In the geological repository of high-level radioactive wastes (HLRW), one of the most important issues is the stability of artificial barrier basically using bentonite. However, due to the exist gaps between bentonite block and surround rock mass, the bentonite can swell relatively freely, whose constraining condition is quite different from the normal swelling pressure test. As an extreme case, the bentonite may engage completely into the surrounding ground and become a complete normal consolidated state. Therefore, in order to investigate the thermo-mechanical behavior of saturated bentonite at normally consolidated state, a newly proposed static compaction method for preparing specimen of saturated normally consolidated bentonite was proposed first. Then, a series of undrained triaxial compression tests were conducted under different confining pressure and temperature, by which the thermo-mechanical behavior of bentonite, used as an artificial barrier of geological repository of HLRW, is investigated systematically. The tests results reveal that the specimen prepared by the proposed method is confirmed to have sufficient saturation (> 0.92), meanwhile the prepared specimen is at slightly overconsolidated state that roughly equal to a K0-line consolidation. The deformation of the bentonite in undrained triaxial compression tests changes from plastic to brittle as the temperature increases. Simulation based on a sophisticated constitutive model considering water compressibility at high confining pressure was also discussed. The results are useful for assessing the mechanical behavior of bentonite after hundreds years of geological repository of HLRW.

Thermo-mechanical modelling of spalling around the deposition boreholes in an underground nuclear waste repository during its thermal phase

This paper presents a three-dimensional numerical analysis of multiple fracture growth leading to the development of excavation disturbed zones and spalling around deposition boreholes in a geological disposal facility. The development of fracture patterns is simulated with the Imperial College Geomechanics Toolkit, a finite-element based simulator that can model the simultaneous nucleation, growth, and coalescence of multiple fractures in quasi-brittle rock. In these simulations, fractures develop due to the stress concentrations around the borehole wall, caused by the local in situ stresses, and due to the thermal stresses caused by the radioactive decay of the waste. Fracture patterns, and the extent of the spalled zone, are computed after the borehole drilling, heating, and cooling stages, at the Forsmark repository site in Sweden. The effect of temperature on the nucleation and growth of spalling fractures, as well as on the reactivation of pre-existing fractures, is assessed qualitatively, by comparing fracture patterns, and quantitatively, in terms of the maximum spalling depth, width, and increase in the total fractured surface area. Overall, the simulations presented herein indicate that thermal spalling will increase the depths (away from the borehole) and angular widths of the spalled zone, but is not likely to lead to major increases in fracture aperture, and concomitant increases in hydraulic transmissivity and permeability of the spalled zone, above that which has already been caused by mechanical spalling.

Visualization of hydraulic fracturing in compacted bentonite: The roles of dry density, water content, and pressurization rate

Deep geological repository is typically situated at depths ranging from several hundred to 1000 m below ground, making bentonite engineered barrier potentially vulnerable to high water pressure and even inducing hydraulic fracturing. This study conducted injection tests on compacted GMZ (Gaomiaozi) bentonite with a self-developed visualization set-up. The objective was to unveil the roles of dry density, water content, and pressurization rate in hydraulic fracturing from the perspective of fracturing macro-morphological dynamics and breakthrough characteristics. Moreover, the relationships between breakthrough characteristics and microstructure were examined by MIP (mercury intrusion porosimetry) analysis. Results showed that the fracturing dynamics were characterized by three stages: hydration, cracking, and fracturing stages. Compared to water content and pressurization rate, dry density exerted more pronounced effects on these stages. Increasing dry density can lead to an expansion of circular hydration zone, a more complex cracking network, and a change in fracturing patterns from long and clear to short and fuzzy. In terms of breakthrough characteristics, the breakthrough pressure was positively correlated with dry density and negatively correlated with water content. Interestingly, there is a good and unique logarithmic correlation between the breakthrough pressure and the ratio eM/em of inter-aggregate void ratio and intra-aggregate void ratio, regardless of dry density and water content. Within a certain range (i.e. 200-50 kPa/min), breakthrough pressure showed slight dependency on pressurization rate. Nevertheless, an extremely low pressurization rate of 20 kPa/min caused a transition for the specimen from quasi-brittle to plastic state owning to more water infiltration, thereby hindering fracture initiation and propagation.

The value of simplified models of radionuclide transport for the safety assessment of nuclear waste repositories: A benchmark study

In order to assess sites for a deep geological repository for storing high-level nuclear waste safely in Germany, various numerical models and tools will be in use. For their interaction within one workflow, their reproducibility, and reliability version-controlled open-source solutions and careful documentation of model setups, results and verifications are of special value. However, spatially fully resolved models including all relevant physical and chemical processes are neither computationally feasible for large domains nor is the data typically available to parameterize such models. Thus, simplified models are crucial for the pre-assessment of possible sites to narrow down the list of suitable candidates for which detailed site investigations and fully resolved models will be done at a later stage. Still, the accuracy of these simplified models is of importance as the pre-assessment of suitable sites will be based on them. In this study, we compare the modelling capabilities of TransPyREnd, a one-dimensional transport code based on finite differences, specifically developed for the fast estimation of radionuclide transport by the German federal company for radioactive waste disposal (BGE), with OpenGeoSys, which is a modelling platform based on finite elements in up to three spatial dimensions. Both codes are used in the site selection procedure for the German nuclear waste repository. The comparison of the model results of TransPyREnd and OpenGeoSys is augmented by comparisons with an analytical solution for a homogeneous material. For the purpose of numerical benchmarking, we consider a geological profile located in southern Germany as an example where the hypothetical repository is located in a clay-stone formation. TransPyREnd and OpenGeoSys yield overall similar results. However, both codes use different discretizations which impact is the highest for strongly sorbing compounds, while the difference gets negligible for less sorbing and more diffusive compounds as higher diffusion tends to blur the initial conditions. Overall, the OpenGeoSys model is more exact whereas the TransPyREnd model has considerable faster run times. We found in our example, that significant substance amounts are only leaving the host rock formation, if apparent diffusion is high, for which case both codes give similar results, while relative differences are considerable for strongly sorbing compounds. However, in the latter case no significant substance amount of radionuclides leaves the host-rock formation, thus deeming the differences in the model results minor for the overall safety assessments of sites.

Oxic corrosion model for KAERI Reference disposal system via O2 consumption reactions and mixed-potential theory

The durability of copper (Cu) canisters against corrosion is critical for the licensing of deep geological repositories. Assessing oxic corrosion, a primary degradation mechanism, is essential for ensuring the reliability of such repositories. Due to the complex interactions influencing oxic corrosion, a comprehensive model is necessary for evaluating Cu canister corrosion. This study develops a model for the KAERI Reference Disposal System (KRS), incorporating mixed-potential theory with key O2 consumption reactions, including Cu corrosion, Cu(I) oxidation, FeS2 oxidation, aerobic microbial activity, and O2 dissolution and consumption. Simulation of 11 scenarios revealed that the representative KRS case would experience a maximum corrosion depth of 9.3 μm on the Cu canister after 2.3 years due to oxic corrosion, under conditions that are unfavorable for the initiation of pitting corrosion. These results suggest that oxic corrosion is not a threat to the isolation of spent nuclear fuels in KRS.

Sedimentary basins in the southeastern part of the Laptev sea and the adjacent area of the Siberian platform

Background. The sedimentary cover of the Leno-Anabar trough and the adjacent Anabar-Khatanga saddle show signs of oil and gas bearing capacity, which indicates the presence of generation and accumulation hydrocarbon systems. In addition, the offshore area of the Laptev plate and the adjacent Siberian platform is promising in terms of oil and gas potential, since the sedimentary basins, including hydrocarbon systems, are characterized by a transboundary (onshore/offshore) position. Until now, the offshore and onshore parts of this area have been investigated separately. Such an approach significantly reduces the estimation accuracy of hydrocarbon potential. This determines the relevance of developing a unified regional geological model of the junction zone, which may serve as the basis for further oil and gas geological studies.Aim. Construction of a geological model to trace the development of sedimentary basins within the junction area of the Laptev plate and the Siberian platform.Materials and methods. Collection and analysis of published materials, as well as data from geological repositories related to geological, geophysical, and geochemical investigations of the sedimentary cover in the construction of a geoinformation database; development of a numerical geological model of the studied area; carrying out paleogeographic reconstructions with identification of the main depocenters and domains of consistent subsidence delineation; assessment of sedimentation rates of the main sedimentary complexes.Results. A unified geological model of the junction area of the Laptev plate and the Siberian Platform was created, enabling investigation of the sedimentary basins and delineation of the essential hydrocarbon system elements.

No signs of microbial-influenced corrosion of cast iron and copper in bentonite microcosms after 400 days

High-level radioactive waste needs to be safely stored for a long time in a deep geological repository by using a multi-barrier system. In this system, suitable barrier materials are selected that ideally show long-term stability to prevent early radionuclide release into the biosphere. In this study, different container matals (copper and cast iron) and pore water compositions (Opalinus Clay pore water and saline cap rock solution) were combined with Bavarian bentonite in static batch experiments to investigate microbial-influenced corrosion. The increasing concentration of iron and copper in the solution as well as detected corrosion products on the metal surface are indicative of anaerobic corrosion of the respective metals during an incubation of 400 days at 37 °C. However, although the intrinsic microbial bentonite community was stimulated with either lactate or H2, an acceleration of cast iron- and copper corrosion did not occur. Furthermore, neither corrosive bacteria nor conventional bacterial corrosion products, such as metal sulfides, were detected in any of the analyzed samples. The analyses of geochemical parameters (e.g. ferrous iron-, iron-, copper- and potassium concentrations as well as redox potentials) showed significant changes in some cast iron- and copper-containing setups, but these changes did not correlate with the microbial community structure in the respective microcosms, as confirmed by statistical analyses. Hence, the analyzed Bavarian bentonite (type B25) showed no significant contribution to cast iron and copper corrosion under the applied conditions after 400 days of incubation. From this perspective, bentonite B25 could be a suitable candidate as a geotechnical barrier in future repositories.

Analytical Prediction of Temperature Distribution in a Deep Geological Nuclear Waste Repository

Temperature distribution plays a crucial role in the safety performance assessment and thermal dimensioning design in a deep geological repository for disposing high-level waste. In this study, a two-dimensional axisymmetric model of a single container for heat transfer was created. The fully analytical solution to temperature distribution in the repository was derived by utilizing the methods of separation of variables, impulse theorem, and Fourier transform. The fully analytical solution was validated by comparing with the existing semi-analytical solution and line heat source solution. The temperature change in the near field around the container was analyzed using the present solution, and the influences of different parameters on the container surface temperature were investigated. Furthermore, the proposed fully analytical solution was used to predict the results of the in situ test. The findings indicate that the temperature in the buffer layer rapidly increases and reaches its peak value within the first 2 years, then gradually decreases thereafter with time. The thickness of the bentonite pellet layer had a greater effect on the container surface temperature than that of the bentonite block layer. A comparison between the fully analytical solution and the results of the in situ heating test demonstrated that the proposed fully analytical solution can accurately predict the temperature variations in the in situ heating test.

Comparative study of hydro-mechanical behaviors of compacted bentonite powder and granular bentonite

In the deep geological disposal repository of high-level radioactive waste, buffer/backfill materials typically consist of compacted bentonite block and granular bentonite. As these materials undergo a long-term hydration, it is anticipated that the two forms of bentonite materials (i.e. compacted bentonite powder (CBP) and granular bentonite (GB)) are expected to exhibit differing hydro-mechanical behaviors due to the differences in their structures. This work aims to investigate the differences in swelling pressure and compressibility through a series of swelling pressure tests, compression tests and mercury intrusion porosimetry (MIP) tests. The experimental results demonstrated that swelling pressure curves of the CBP specimens showed higher first peak values and more pronounced collapse than those of the GB specimens at a given dry density, regardless of vapor-water hydration or liquid-water hydration. The final swelling pressures of the two materials were similar at the same dry density, suggesting an independent correlation between swelling pressure and dry density. At the high suction range, the compression curves exhibited an obvious bi-linear pattern for the CBP specimens and a significant nonlinearity for the GB specimens. Meanwhile, the CBP specimens presented higher pre-consolidation pressures and larger compression indices than the GB specimens at a given suction. As suction decreased, the compression curves of the two materials gradually approached each other and their differences were reduced accordingly. After reaching saturation, a good consistency between them was observed whether for final swelling pressure or compressibility. Pore structure analysis revealed that the two materials both presented an initially double structure, and their differences were primarily manifested at the macrostructural level. Eventually, the differences in swelling pressure or compression curves of the two materials were well interpreted by combining microstructural evolutions.

Kinematic and paleostress analyses in the KURT site to reconstruct the long-term structural evolution to apply for deep geological repository of nuclear waste

Kinematic and paleostress analyses in the KURT site to reconstruct the long-term structural evolution to apply for deep geological repository of nuclear waste

Hybrid Twins Modeling of a High-Level Radioactive Waste Cell Demonstrator for Long-Term Temperature Monitoring and Forecasting.

Monitoring a deep geological repository for radioactive waste during the operational phases relies on a combination of fit-for-purpose numerical simulations and online sensor measurements, both producing complementary massive data, which can then be compared to predict reliable and integrated information (e.g., in a digital twin) reflecting the actual physical evolution of the installation over the long term (i.e., a century), the ultimate objective being to assess that the repository components/processes are effectively following the expected trajectory towards the closure phase. Data prediction involves using historical data and statistical methods to forecast future outcomes, but it faces challenges such as data quality issues, the complexity of real-world data, and the difficulty in balancing model complexity. Feature selection, overfitting, and the interpretability of complex models further contribute to the complexity. Data reconciliation involves aligning model with in situ data, but a major challenge is to create models capturing all the complexity of the real world, encompassing dynamic variables, as well as the residual and complex near-field effects on measurements (e.g., sensors coupling). This difficulty can result in residual discrepancies between simulated and real data, highlighting the challenge of accurately estimating real-world intricacies within predictive models during the reconciliation process. The paper delves into these challenges for complex and instrumented systems (multi-scale, multi-physics, and multi-media), discussing practical applications of machine and deep learning methods in the case study of thermal loading monitoring of a high-level waste (HLW) cell demonstrator (called ALC1605) implemented at Andra's underground research laboratory.

Predicting bentonite swelling pressure: optimized XGBoost versus neural networks

The swelling pressure of bentonite and bentonite mixtures is critical in designing barrier systems for deep geological radioactive waste repositories. Accurately predicting the maximum swelling pressure is essential for ensuring these systems' long-term stability and sealing characteristics. In this study, we developed a constrained machine learning model based on the extreme gradient boosting (XGBoost) algorithm tuned with grey wolf optimization (GWO) to determine the maximum swelling pressure of bentonite and bentonite mixtures. A dataset containing 305 experimental data points was compiled, including relevant soil properties such as montmorillonite content, liquid limit, plastic limit, plasticity index, initial water content, and soil dry density. The GWO-XGBoost model, incorporating a penalty term in the loss function, achieved an R2 value of 0.9832 and an RMSE of 0.5248 MPa in the testing phase, outperforming feed-forward and cascade-forward neural network models. The feature importance analysis revealed that dry density and montmorillonite content were the most influential factors in predicting maximum swelling pressure. While the developed model demonstrates high accuracy and reliability, it may have limitations in capturing extreme values due to the complex nature of bentonite swelling behavior. The proposed approach provides a valuable tool for predicting the maximum swelling pressure of bentonite-based materials under various conditions, supporting the design and analysis of effective barrier systems in geotechnical engineering applications.

A CSS surface based THM bounding constitutive model for volumetric behavior of bentonite

Development of the constitutive model for bentonite under coupled thermo-hydro-mechanical (THM) condition is of great significance for the construction and safety assessment of deep geological disposal repositories for high-level radioactive waste (HLW). In this work, a new temperature-suction-mean net stress (T-s-p) space with the conception of critical saturated state (CSS) surface was defined to represent the actual stress state of bentonite under coupled THM condition. Then, based on the CSS surface, a THM constitutive model was proposed for describing the volumetric behavior of compacted bentonite. Under the THM model framework, two bounding surfaces were proposed to describe the elastoplastic volume changes induced by mechanical response of skeleton and hydration of montmorillonite, respectively. The model responses upon some typical THM paths were simulated and discussed to reveal the performance of the proposed constitutive model for bentonite. Finally, the proposed model was validated by simulating by several volume change tests carried out on different bentonites. The results confirmed that the proposed model shows more advantages in describing the THM volumetric behavior.

Gas migration at the granite–bentonite interface under semirigid boundary conditions in the context of high‐level radioactive waste disposal

AbstractThe corrosion of waste canisters in the deep geological disposal facilities (GDFs) for high‐level radioactive waste (HLRW) can generate gas, which escapes from the engineered barrier system through the interfaces between the bentonite buffer blocks and the host rock and those between the bentonite blocks. In this study, a series of water infiltration and gas breakthrough experiments were conducted on granite and on granite–bentonite specimens with smooth and grooved interfaces. On this basis, this study presents new insights and a quantitative assessment of the impact of the interface between clay and host rock on gas transport. As the results show, the water permeability values from water infiltration tests on granite and granite–bentonite samples (10−19–10−20 m2) are found to be slightly higher than that of bentonite. The gas permeability of the mock‐up samples with smooth interfaces is one order of magnitude larger than that of the mock‐up with grooved interfaces. The gas results of breakthrough pressures for the granite and the granite–bentonite mock‐up samples are significantly lower than that of bentonite. The results highlight the potential existence of preferential gas migration channels between the rock and bentonite buffer that require further considerations in safety assessment.

Numerical assessment of the time-dependent behaviour of a tunnel with a yielding support system for a potential repository in claystone

For the disposal of heat-generating radioactive waste and spent fuel in claystone formations, support structures for underground excavations are essential for the operational safety of geological disposal facilities (GDF). Claystone formations of moderate strength, at considerable depths, are characterised by their squeezing, creeping, and, sometimes, swelling behaviour that results in continuous tunnel convergence over time. Under these conditions, a rigid support system can be subjected to very high loads, requiring a considerable thickness and/or high-performance concretes. Therefore, yielding support systems are being currently investigated as a promising alternative for GDF. This work presents a numerical study of the behaviour of a tunnel in a claystone formation with a yielding support system. The use of a compressible mortar between the rock and the lining is assumed, as a means of reducing loads and mitigating the effects of creep deformations. A key aspect of the analyses is that the host rock is characterised by a constitutive model that includes a number of features that are relevant for the satisfactory description of the hydro-mechanical behaviour of stiff clayey materials, such as mechanical anisotropy, creep, strain softening, and its ability to simulate localised deformations through a nonlocal regularisation. An elastoplastic constitutive model was also developed to represent the behaviour of the compressible mortar. Results provide relevant insights into the performance of the adopted yielding support system, particularly regarding the effect of time-dependent deformations and the additional relaxation of the rock on the fractured zone near the excavation.

Archaeological analogues for the lifetime prediction of the canister for spent nuclear fuel in the deep geological repository: Part I. Composition of pore water and corrosion products

AbstractThis work focuses on archaeological analogues as a way of obtaining long‐term data for assessing the lifetime of a spent nuclear fuel canister. An analysis of the environment around the excavated artefacts is presented, followed by geochemical modelling of the likely corrosion products and subsequent comparison with real phases in the corrosion layer estimated by X‐ray diffraction. A total of 16 archaeological sites with fine clay soils and the potential for long‐term flooding were evaluated. Although the soil pore solution contained high levels of silicates, no such phases with iron in the corrosion products were confirmed. Carbonate corrosion products, typical of bentonite environments, were also not observed. Although oxygen access is very limited in the environments of all sites, even low concentrations shifted the equilibrium of corrosion products formed in favour of (hydroxy)oxides.

Value of abstraction in performance assessment – When is a higher level of detail necessary?

In this study, different approaches in performance assessment (PA) of the long-term safety of a repository for radioactive waste were examined. This investigation was carried out as part of the DECOVALEX-2023 project, an international collaborative effort for research and model comparison. One specific task of the DECOVALEX-2023 project was the Salt Performance Assessment Modelling task (Salt PA), which aimed at comparing various models and methods employed in the performance assessment of deep geological repositories in salt. In the context of the Salt PA task, three distinct teams from SNL (United States), Quintessa Ltd (United Kingdom), and GRS (Germany) examined the consequences of employing different levels of abstractions when modelling the repository's geometry and implementing various features and processes, using the example of a simple hypothetical repository structure in domal salt. Each team applied their own tools: PFLOTRAN (SNL), QPAC (Quintessa) and LOPOS (GRS). These differ essentially regarding numerical concept and degree of detail in the representation of the underlying physical processes. The discussion focused on when simplifications can be appropriately applied and what consequences result from them. Furthermore, it was explored when and if a higher level of fidelity in geometry or physical processes is required.

Unveiling the geoheritage, cultural geomorphology and geotourism potential of Zanskar region, NW Himalaya, India

The Zanskar region, nestled in the heart of the Tethyan Himalaya, represents a significant geological and cultural repository that has remained relatively unexplored by the global tourism industry. This study investigates Zanskar's geoheritage and cultural geomorphology, highlighting its importance for geotourism and advocating for the implementation of sustainable development to preserve its fragile ecosystem. Flowing northward from the Higher Himalayas, the Zanskar River, one of the largest northeast-flowing tributary of the Upper Indus River Basin (UIRB), traverses the complex, strongly folded, and thrusted Zanskar Range via a steep, narrow 60-kilometer-long gorge. Zanskar has consistently garnered attention from geoscientists and tourists owing to its unique geographic location and remoteness. Many areas and routes in Zanskar are challenging to reach, but it stands out as a prime location where glaciers can be accessed by road, making it an ideal destination for studying geoscience. Geologically, the study area is part of the Zanskar Range, with its boundaries defined by the Ladakh Range to the northeast and the Great Himalayan Range to the southwest. The Indus Tsangpo Suture Zone, the Tethyan Sedimentary Sequence, and the Higher Himalayan Crystalline Sequence are exposed from north to south. Zanskar serves as a pristine field museum for studying young Himalayan orogenic tectonics, featuring prominent exposures such as the Zanskar Shear Zone (ZSZ)/South Tibetan Detachment System (STDS) and the structural outcrop of Zangla, showcasing a variety of folding types. These folds signify different phases of deformation resulting from the Indo-Eurasian collision. Moreover, it serves as a prime sequence of Quaternary Glacial Periods (QGP), providing an ideal natural laboratory and an open, practical classroom for geoscientists and Earth science students. We propose ten geoheritage sites in Zanskar, which encompass cultural, geomorphological, and geological attractions with potential for geotourism in this region.

Thermal effects on the strain rate-dependent behavior of highly compacted GMZ01 bentonite

Investigation of thermal effects on the strain rate-dependent properties of compacted bentonite is crucial for the long-term safety assessment of deep geological repository for disposal of high-level radioactive waste. In the present work, cylindrical GMZ01 bentonite specimens were compacted with suction-controlled by the vapor equilibrium technique. Then, a series of temperature- and suction-controlled stepwise constant rate strain (CRS) tests was performed and the rate-dependent compressibility behavior of the highly compacted GMZ01 bentonite was investigated. The plastic compressibility parameter λ, the elastic compressibility parameter κ, the yield stress p0, as well as the viscous parameter α were determined. Results indicate that λ, κ and α decrease and p0 increases as suction increases. Upon heating, parameters λ, α and p0 decrease. It is also found that p0 increases linearly with increasing CRS in a double-logarithm coordinate. Based on the experimental results, a viscosity parameter α(s, T) was fitted to capture the effects of suction s and temperature T on the relationship between yield stress and strain rate. Then, an elastic-thermo-viscoplastic model for unsaturated soils was developed to describe the thermal effects on the rate-dependent behavior of highly compacted GMZ01 bentonite. Validation showed that the calculated results agreed well to the measured ones.

Mineralogical and geochemical composition of a cementitious grout and its evolution during interaction with water

In the present study, the chemical composition, mineralogy, and mechanisms of alteration of a cementitious grout based on a CEM III/C with addition of smectite, hydrotalcite, and silica fume, are studied using a combination of chemical and physical methods. This material was designed in the context of geological repository of radioactive wastes, with a twofold aim: first, to fill the technical voids left by drilling operations at the interface between the geological formation and the disposal galleries. Second, to neutralize a potential acidic transient due to pyrite oxidation, and to create an environment that favors low corrosion rates of carbon steels. The grout is mainly composed of calcium silicate hydrates having a Ca/Si ratio of ~0.8, incorporating Al in the bridging site of the Si chains (C-A-S-H), and accounting for 29–36 wt.% of the sample. It also contains silica fume (38–48 wt.%), smectite with interlayer Na (11–17 wt.%), hydrotalcite with interlayer CO32− (3–4 wt.%), and lower amounts of portlandite, calcite, and possibly gibbsite and gypsum. Upon alteration by water in a flow-through reactor, the main modifications affecting the sample are calcite and gypsum dissolution, hence releasing aqueous Ca2+ that is adsorbed in smectite interlayer by replacing Na+, and stoichiometric C-A-S-H dissolution. The evolution of solution chemistry and of the solid phase composition are reproduced successfully using a thermokinetic model.