Mobility of daughter elements of 238U decay chain during leaching by In Situ Recovery (ISR): New insights from digital autoradiography

In highly permeable sedimentary rock formations, U extraction by in-situ leaching techniques (ISR - In-Situ Recovery) is generally considered to have a limited environmental impact at ground level. Significantly, this method of extraction produces neither mill tailings nor waste rocks. Underground, however, the outcome for ²³⁸U daughter elements in aquifers is not well known because of their trace concentrations in the host rocks. Thus, understanding the in-situ mobility of these elements remains a challenge. Two samples collected before and after six months of ISR experiments (Dulaan Uul, Mongolia) were studied with the help of a digital autoradiography technique (DA) of alpha particles, bulk alpha spectrometry, and complementary petrographic observation methods. These techniques demonstrate that before and after leaching, the radioactivity is concentrated in altered and microporous Fe-Ti oxides. Most of the daughter elements of U remain trapped in the rock after the leaching process. DA confirms that the alpha activity of the Fe-Ti oxides remains high after uranium leaching, and the initial secular equilibrium of the ²³⁸U series for ²³⁰Th to ²¹⁰Po daughter elements (including ²²⁶Ra) of the fresh rocks is maintained after leaching. While these findings should be confirmed by more systematic studies, they already identify potential mechanisms explaining why the U-daughter concentrations in leaching water are low.

Imaging connected porosity of crystalline rock by contrast agent-aided X-ray microtomography and scanning electron microscopy

We set out to study connected porosity of crystalline rock using X-ray microtomography and scanning electron microscopy with energy dispersive X-ray spectroscopy (SEM-EDS) with caesium chloride as a contrast agent. Caesium is an important radionuclide regarding the final deposition of nuclear waste and also forms dense phases that can be readily distinguished by X-ray microtomography and SEM-EDS. Six samples from two sites, Olkiluoto (Finland) and Grimsel (Switzerland), where transport properties of crystalline rock are being studied in situ, were investigated using X-ray microtomography and SEM-EDS. The samples were imaged with X-ray microtomography, immersed in a saturated caesium chloride (CsCl) solution for 141, 249 and 365 days and imaged again with X-ray microtomography. CsCl inside the samples was successfully detected with X-ray microtomography and it had completely penetrated all six samples. SEM-EDS elemental mapping was used to study the location of caesium in the samples in detail with quantitative mineral information. Precipitated CsCl was found in the connected pore space in Olkiluoto veined gneiss and in lesser amounts in Grimsel granodiorite. Only a very small amount of precipitated CsCl was observed in the Grimsel granodiorite samples. In Olkiluoto veined gneiss caesium was found in pinitised areas of cordierite grains. In the pinitised areas caesium was found in notable excess compared to chloride, possibly due to the combination of small pore size and negatively charged surfaces. In addition, elevated concentrations of caesium were found in kaolinite and sphalerite phases. The findings concerning the location of CsCl were congruent with X-ray microtomography.

Determination of Granites' Mineral Specific Porosities by PMMA Method and FESEM/EDAX

Over extended periods, long-lived radionuclides (RN) or activation products within geologic disposal sites may be released from the fuel and migrate to the geo/biosphere. In the bedrock, contaminants will be transported along fractures by advection and retarded by sorption on mineral surfaces and by molecular diffusion into stagnant pore water in the matrix along a connected system of pores and micro-fissures.

The objective of this paper was to determine the connective porosity and mineral-specific porosities for three granite samples by 14C methylmethacrylate (14C-PMMA) autoradiography. Scanning electron microscopy and energy-dispersive X-ray analyses (FESEM/EDAX) were performed in order to study the pore apertures of porous regions in greater detail and to identify the corresponding minerals. Finally, the porosity results were used to evaluate the diffusion coefficients of RNs from previous experiments which determined apparent diffusion coefficients for the main minerals in three granite samples by the Rutherford Backscattering technique.

The total porosity of the Grimsel granite (0.75%) was significantly higher than the porosities of the El Berrocal and Los Ratones granites (0.3%). The porosities of the Grimsel granite feldspars were two to three times higher than the porosities of the El Berrocal and Los Ratones granites' feldspars. However, there was no significant difference between the porosities of the dark minerals. A clear difference was found between the various quartz grains. Quartz crystals were nonporous in the El Berrocal and Los Ratones granites when measured by the PMMA method, but the quartz crystals in the Grimsel granite showed 0.5% intra granular porosity. The apparent diffusion coefficients calculated for uranium diffusion within Grimsel granite on different minerals were very similar (210-13 ± 0.5 m2/s), but differences within both Spanish granites were found from one mineral to another (9 ± 110-14 m2/s in feldspars and 4.5 ± 0.510-14 m2/s in quartz) - always presenting lower diffusion values than in the Grimsel granite.

Direct simulation of heterogeneous diffusion and inversion procedure applied to an out-diffusion experiment. Test case of Palmottu granite

An out-diffusion laboratory experiment using a non-reactive tracer was fitted using the Time Domain Diffusion (TDD) method. This rapid particle tracking method allows simulation of the heterogeneous diffusion based on pore-scale images and local values of diffusivities. The superimposed porosity and mineral 2D maps act as computation grids to condition diffusion pathways. We focused on a Palmottu granite sample, in which the connected pore space has a composite microstructure with cracks linking microporous minerals and is above the percolation threshold. Three main results were achieved: (i) When compared to the fitting obtained with one coefficient (best mean square residual R = 1.6 x 10(-2)), diffusion is shown to be suitably characterised with two coefficients related to cracks and microporous minerals (best R = 6.5 x 10(-4)), (ii) rather than imposing a local apparent diffusion coefficient D(a) independent of the local porosity Phi, a best fit is obtained by applying Archie's relationship D(a) = D(0) x G with G = Phi(m) to each pixel of the calculation grids (G is the geometry factor, D(0) is the diffusion coefficient in free fluid, and m is Archie's exponent), and (iii) the order of magnitude of the fitted diffusion coefficient or Archie's exponents (m=0 for microcracks and m=1.82 for microporous minerals) is physically realistic.

Attempt to model laboratory-scale diffusion and retardation data

Different approaches for measuring the interaction between radionuclides and rock matrix are needed to test the compatibility of experimental retardation parameters and transport models used in assessing the safety of the underground repositories for the spent nuclear fuel. In this work, the retardation of sodium, calcium and strontium was studied on mica gneiss, unaltered, moderately altered and strongly altered tonalite using dynamic fracture column method. In-diffusion of calcium into rock cubes was determined to predict retardation in columns. In-diffusion of calcium into moderately and strongly altered tonalite was interpreted using a numerical code FTRANS. The code was able to interprete in-diffusion of weakly sorbing calcium into the saturated porous matrix. Elution curves of calcium for the moderately and strongly altered tonalite fracture columns were explained adequately using FTRANS code and parameters obtained from in-diffusion calculations. In this paper, mass distribution ratio values of sodium, calcium and strontium for intact rock are compared to values, previously obtained for crushed rock from batch and crushed rock column experiments. Kd values obtained from fracture column experiments were one order of magnitude lower than Kd values from batch experiments.