Biogeochemical processes in a clay formation in situ experiment: Part F - Reactive transport modelling

Reactive transport modelling was used to simulate solute transport, thermodynamic reactions, ion exchange and biodegradation in the Porewater Chemistry (PC) experiment at the Mont Terri Rock Laboratory. Simulations show that the most important chemical processes controlling the fluid composition within the borehole and the surrounding formation during the experiment are ion exchange, biodegradation and dissolution/precipitation reactions involving pyrite and carbonate minerals. In contrast, thermodynamic mineral dissolution/precipitation reactions involving alumo-silicate minerals have little impact on the fluid composition on the time-scale of the experiment. With the accurate description of the initial chemical condition in the formation in combination with kinetic formulations describing the different stages of bacterial activities, it has been possible to reproduce the evolution of important system parameters, such as the pH, redox potential, total organic C. dissolved inorganic C and SO(4) concentration. Leaching of glycerol from the pH-electrode may be the primary source of organic material that initiated bacterial growth, which caused the chemical perturbation in the borehole. Results from these simulations are consistent with data from the over-coring and demonstrate that the Opalinus Clay has a high buffering capacity in terms of chemical perturbations caused by bacterial activity. This buffering capacity can be attributed to the carbonate system as well as to the reactivity of clay surfaces.

Hoisting C Structures Into Clay In Device Drivers

This project addresses the unreliability of operating system code, in particular in device drivers. Device driver software is the interface between the operating system and the device's hardware. Device drivers are written in low level code, making them difficult to understand. Almost all device drivers are written in the programming language C which allows for direct manipulation of memory. Due to the complexity of manual movement of data, most mistakes in operating systems occur in device driver code. The programming language Clay can be used to check device driver code at compile-time. Clay does most of its error checking statically to minimize the overhead of run-time checks in order to stay competitive with C's performance time. The Clay compiler can detect a lot more types of errors than the C compiler like buffer overflows, kernel stack overflows, NULL pointer uses, freed memory uses, and aliasing errors. Clay code that successfully compiles is guaranteed to run without failing on errors that Clay can detect. Even though C is unsafe, currently most device drivers are written in it. Not only are device drivers the part of the operating system most likely to fail, they also are the largest part of the operating system. As rewriting every existing device driver in Clay by hand would be impractical, this thesis is part of a project to automate translation of existing drivers from C to Clay. Although C and Clay both allow low level manipulation of data and fill the same niche for developing low level code, they have different syntax, type systems, and paradigms. This paper explores how C can be translated into Clay. It identifies what part of C device drivers cannot be translated into Clay and what information drivers in Clay will require that C cannot provide. It also explains how these translations will occur by explaining how each C structure is represented in the compiler and how these structures are changed to represent a Clay structure.

Microscopic and Spectroscopic Methods for Probing the Clay Water Interface

Clay minerals have a fundamental importance in many processes in soils and sediments such as the bioavailability of nutrients, water retention, the adsorption of common pollutants, and the formation of an impermeable barrier upon swelling. Many of the properties of clay minerals are due to the unique environment present at the clay mineral/water interface. Traditional techniques such as X-ray diffraction (XRD) and absorption isotherms have provided a wealth of information about this interface but have suffered from limitations. The methods and results presented herein are designed to yield new experimental information about the clay mineral/water interface.A new method of studying the swelling dynamics of clay minerals was developed using in situ atomic force microscopy (AFM). The preliminary results presented here demonstrate that this technique allows one to study individual clay mineral unit layers, explore the natural heterogeneities of samples, and monitor swelling dynamics of clay minerals in real time. Cation exchange experiments were conducted monitoring the swelling change of individual nontronite quasi-crystals as the chemical composition of the surrounding environment was manipulated several times. A proof of concept study has shown that the changes in swelling are from the exchange of interlayer cations and not from the mechanical force of replacing the solution in the fluid cell. A series of attenuated total internal reflection Fourier transform infrared spectroscopy (ATR-FTIR) experiments were performed to gain a better understanding of the organization of water within the interlayer region of two Fe-bearing clay minerals. These experiments made use of the Subtractive Kramers-Kronig (SKK) Transform and the calculation of difference spectra to obtain information about interfacial water hidden within the absorption bands of bulk water. The results indicate that the reduction of structural iron disrupts the organization of water around a strongly hydrated cation such as sodium as the cation transitions from an outer-sphere complex with the mineral surface to an inner-sphere complex. In the case of a less strongly hydrated cation such as potassium, reduction of structural iron actually increases the ordering of water molecules at the mineral surface. These effects were only noticed with the reduction of iron in the tetrahedral sheet close to the basal surface where the increased charge density is localized closer to the cations in the interlayer.

Translational diffusion of water in compacted clay systems

The translational diffusion of water in compacted clays at a high hydration level has been investigated by quasielastic neutron scattering at a time-of-flight spectrometer FOCUS (SINQ). Four compacted clays with systematic structural differences have been studied: Na-montmorillonite, Na-illite, kaolinite and pyrophyllite. The QENS experiments were performed using two different incident wavelengths in order to access a larger Q range and verify the data analysis. The translational diffusion coefficient for water in Na-montmorillonite and Na-illite are lower than those for bulk water, whereas the preliminary results for kaolinite and pyrophyllite show larger diffusion coefficient.

Characterisation of Gas Transport Properties of the Opalinus Clay, a Potential Host Rock Formation for Radioactive Waste Disposal

The Opalinus Clay in Northern Switzerland has been identified as a potential host rock formation for the disposal of radioactive waste. Comprehensive understanding of gas transport processes through this low-permeability formation forms a key issue in the assessment of repository performance. Field investigations and laboratory experiments suggest an intrinsic permeability of the Opalinus Clay in the order of 10(-20) to 10(-21) m(2) and a moderate anisotropy ratio < 10. Porosity depends on clay content and burial depth; values of similar to 0.12 are reported for the region of interest. Porosimetry indicates that about 10-30 of voids can be classed as macropores, corresponding to an equivalent pore radius > 25 nm. The determined entry pressures are in the range of 0.4-10 MPa and exhibit a marked dependence on intrinsic permeability. Both in situ gas tests and gas permeameter tests on drillcores demonstrate that gas transport through the rock is accompanied by porewater displacement, suggesting that classical flow concepts of immiscible displacement in porous media can be applied when the gas entry pressure (i.e. capillary threshold pressure) is less than the minimum principal stress acting within the rock. Essentially, the pore space accessible to gas flow is restricted to the network of connected macropores, which implies a very low degree of desaturation of the rock during the gas imbibition process. At elevated gas pressures (i.e. when gas pressure approaches the level of total stress that acts on the rock body), evidence was seen for dilatancy controlled gas transport mechanisms. Further field experiments were aimed at creating extended tensile fractures with high fracture transmissivity (hydro- or gasfracs). The test results lead to the conclusion that gas fracturing can be largely ruled out as a risk for post-closure repository performance.

Transport of 234U in the Opalinus Clay on centimetre to decimetre scales

The Opalinus Clay formation in North Switzerland is a potential host rock for a deep underground radioactive waste repository. The distribution of U-238, U-234 and Th-230 was studied in rock samples of the Opalinus Clay from an exploratory borehole at Benken (Canton of Zurich) using MC-ICP-MS. The aim of U-234 was to assess the in situ, long-term migration behaviour in this rock. Very low hydraulic conductivities of the Opalinus Clay, reducing potential of the pore water and its chemical equilibrium with the host rock are expected to render both U-238 and Th-230 immobile. If U is heterogeneously distributed in the Opalinus Clay, gradients in the supply of U-234 from the rock matrix to the pore water by the decay of U-238 will be established. Diffusive redistribution separates U-234 from its immobile parent U-238 resulting in bulk rock U-234/U-238 activity disequilibria. These may provide a means of estimating the mobility of U-234 in the rock if the diffusion rate of U-234 is significant compared to its decay rate. Sampling was carried out on two scales. Drilling of cm-spaced samples from the drill-core was done to study mobility over short distances and elucidate possible small-scale lithological control. Homogenized 25-cm-long portions of a 2-m-long drill-core section were prepared to provide information on transport over a longer distance. Variations in U and/or Th content on the cm-scale between clays and carbonate-sandy layers are revealed by beta-scanning, which shows that the (dominant) clay is richer in both elements. Samples were digested using aqua regia followed by total HF dissolution, yielding two fractions. in all studied samples U was found to be concentrated in the HF digestion fraction. It has a high U/Th ratio and a study by SEM-EDS points to sub-mu m up to several mu m in size zircon grains as the main U-rich phase. This fraction consistently has U-234/U-238 activity ratios below unity. The minute zircon grains constitute the major reservoir of U in the rock and act as constant rate suppliers of U-234 into the rock matrix and the pore water. The aqua regia leach fraction was found to be enriched in Th, and complementary to the HF fraction, having U-234/U-238 activity ratios above unity. It is believed that these U activity ratios reflect the surplus of having U-234 delivered from the zircon grains. Some cm-spaced samples show bulk rock U-234/U-238 activity ratios that are markedly out of equilibrium. In most of them a striking negative correlation between the total U content and the bulk rock U-234/U-238 activity ratios is observed. This is interpreted to indicate net U-234 transfer from regions of higher supply of U-234 towards those of lower supply which is, in most cases, equivalent to transfer from clayey towards carbonate/sandy portions of the rock. In contrast, the 25 cm averaged samples all have uniform bulk rock U-234/U-238 activity ratios in equilibrium, indicating U immobility in the last 1-1.5 Ma on this spatial scale. It is concluded that the small-scale lithological variations which govern U spatial distribution in the Opalinus Clay are the major factor determining U-234 in situ supply rates, regulating its diffusive fluxes and controlling the observed bulk rock U-234/U-238 activity ratios. A simple box-model is presented to simulate the measured bulk rock U-234/U-238 activity ratios and to give an additional insight into the studied system. (C) 2008 Elsevier Ltd. All rights reserved.

In-situ diffusion of HTO, 22Na+, Cs+ and I- in Opalinus Clay at the Mont Terri underground rock laboratory

The diffusion properties of the Opalinus Clay were studied in the underground research laboratory at Mont Terri (Canton Jura, Switzerland) and the results were compared with diffusion data measured in the laboratory on small-scale samples. The diffusion of HTO, Na-22(+), Cs+ and I- were investigated for a period of 10 months. The diffusion equipment used in the field experiment was designed in such a way that a solution of tracers was circulated through a sintered metal screen placed at the end of a borehole drilled in the formation. The concentration decrease caused by the diffusion of tracers into the rock could be followed with time and allowed first estimations of the effective diffusion coefficient. After 10 months, the diffusion zone was over-cored and the tracer profiles measured. From these profiles, effective diffusion coefficients and rock capacity factors Could be extracted by applying a two-dimensional transport model including diffusion and sorption. The simulations were done with the reactive transport code CRUNCH. In addition, results obtained from through-diffusion experiments oil small-sized samples with HTO, Cl-36(-) and Na-22(+) are presented and compared with the in situ data. In all cases. excellent agreement between the two data sets exists. Results for Cs+ indicated five times higher diffusion rates relative to HTO. Corresponding laboratory diffusion measurements are still lacking. However. our Cs+ data are in qualitative agreement wish through-diffusion data for Callovo-Oxfordian argillite rock samples. which also indicate significantly higher effective diffusivities for Cs+ relative to HTO.

Diffusion of HTO, Br-, I-, Cs+, Sr-85(2+) and Co-60(2+) in a clay formation: Results and modelling from an in situ experiment in Opalinus Clay

The migration of radioactive and chemical contaminants in clay materials and argillaceous host rocks is characterised by diffusion and retention processes. Valuable information on such processes can be gained by combining diffusion studies at laboratory scale with field migration tests. In this work, the outcome of a multi-tracer in situ migration test performed in the Opalinus Clay formation in the Mont Terri underground rock laboratory (Switzerland) is presented. Thus, 1.16 x 10(5) Bq/L of HTO, 3.96 x 10(3) Bq/L of Sr-85, 6.29 x 10(2) Bq/L of Co-60, 2.01 x 10(-3) mol/L Cs, 9.10 x 10(-4) mol/L I and 1.04 x 10(-3) mol/L Br were injected into the borehole. The decrease of the radioisotope concentrations in the borehole was monitored using in situ gamma-spectrometry. The other tracers were analyzed with state-of-the-art laboratory procedures after sampling of small water aliquots from the reservoir. The diffusion experiment was carried out over a period of one year after which the interval section was overcored and analyzed. Based on the experimental data from the tracer evolution in the borehole and the tracer profiles in the rock, the diffusion of tracers was modelled with the numerical code CRUNCH. The results obtained for HTO (H-3), I- and Br- confirm previous lab and in situ diffusion data. Anionic fluxes into the formation were smaller compared to HTO because of anion exclusion effects. The migration of the cations Sr-85(2+), Cs+ and Co-60(2+) was found to be governed by both diffusion and sorption processes. For Sr-85(2+), the slightly higher diffusivity relative to HTO and the low sorption value are consistent with laboratory diffusion measurements on small-scale samples. In the case of Cs+, the numerically deduced high diffusivity and the Freundlich-type sorption behaviour is also supported by ongoing laboratory data. For Co, no laboratory diffusion data were yet available for comparison; however, the modelled data suggests that Co-60(2+) sorption was weaker than would be expected from available batch sorption data. Overall, the results demonstrate the feasibility of the experimental setup for obtaining high-quality diffusion data for conservative and sorbing tracers. (C) 2007 Elsevier Ltd. All rights reserved.