Influence of the mode of matrix porosity determination on matrix diffusion calculations

The theoretical basis for matrix diffusion in fractured rocks and the methodology for the determination of diffusion coefficients in the laboratory are well established. One significant problem, however, remains in that it is difficult to quantify the degree of sample disturbance affecting the geometrical, geophysical and hydraulic properties of the rock matrix. A new technique, with in situ rock impregnation with resin, for examining the diffusion-accessible rock matrix has been developed and successfully adopted to the rock matrix behind a water-conducting fracture in host crystalline rocks at Nagra's Grimsel Test Site in Switzerland and JNC's Kamaishi In Situ Test Site in Japan. In line with the results of a large number of natural analogue and laboratory studies, the existence of an in situ interconnected pore network was substantiated. Matrix porosities determined on the laboratory samples from both the sites are 1.5-3 times higher than in situ values, irrespective of the technique applied. On the Grimsel granodiorite matrix, matrix porosity existing in situ and artefacts of stress release and physical disturbance, induced by sampling and sample preparation, were clearly distinguished, allowing in situ porosity to be quantified. Laboratory work with conventional techniques tends to overestimate the porosity of the rock matrix, hence leading to an overestimation of in situ matrix diffusion. The implications of these differences to a repository performance assessment are assessed with a couple of examples from existing assessments, and recommendations for future approaches to the examination of in situ matrix porosity are made.     

The assessment of electromigration as a new technique to study diffusion of radionuclides in clayey soils

For the performance assessment study of a geological disposal of High Level Waste (HLW) in clayey formations, migration studies are essential. For low permeability soils (clays), classical diffusion studies take a very long time. In order to reduce the experimental time, we propose an electrical field as driving force to accelerate the migration of ionic species. This paper reports the assessment of the electromigration technique as a powerful new and fast technique for migration studies. The apparent molecular diffusion coefficient can be derived by two independent methods using the migration parameters obtained from an electromigration experiment, namely the apparent dispersion coefficient and the apparent convection velocity. First, it can be calculated from the velocity of the migrating species by the Einstein relation. But, corrections are necessary for electroosmotic flow. The apparent electroosmotic mobility is experimentally determined as 2.2·10⁻⁹ m²/Vs. Second, it can be calculated from the relation between the apparent dispersion coefficient and the total apparent convection velocity. But it is necessary to know the dispersion length of the medium. The dispersion length for Boom Clay is experimentally determined as 8·10⁻⁵ m. Because of the serious reduction in time, it becomes possible to run series of experiments at different electrical fields to obtain averaged values for the apparent molecular diffusion coefficient according to the two methods. Experiments at different electrical fields have another advantage: the intercept of the linear relationship between the total apparent convection velocity and the apparent dispersion coefficient gives the apparent molecular diffusion coefficient. The apparent molecular diffusion coefficients obtained for ⁸⁵Sr, ¹³¹I and HTO are respectively 0.8·10⁻¹¹, 15·10⁻¹¹, and 24·10⁻¹¹ m²/s. These values are confirmed by pure diffusion experiments. The excellent agreement with the apparent molecular diffusion coefficients obtained by classical diffusion tests clearly demonstrates the feasibility of the electromigration technique for the determination of diffusion coefficients.     

Use of X-ray absorption imaging to examine heterogeneous diffusion in fractured crystalline rocks

Heterogeneous diffusion in different regions of a fractured granodiorite from Japan has been observed and measured through the use of X-ray absorption imaging. These regions include gouge-filled fractures, recrystallized fracture-filling material and hydrothermally altered matrix. With the X-ray absorption imaging technique, porosity, relative concentration, and relative mass of an iodine tracer were imaged in two dimensions with a sub-millimeter pixel size. Because portions of the samples analyzed have relatively low porosity values, imaging errors can potentially impact the results. For this reason, efforts were made to better understand and quantify this error. Based on the X-ray data, pore diffusion coefficients (Dp) for the different regions were estimated assuming a single diffusion rate and a lognormal multirate distribution of Dp. Results show Dp for the gouge-filled fractures are over an order of magnitude greater than those of the recrystallized fracture-filling material, which in turn is approximately two times greater than those for the altered matrix. The recrystallized fracture-filling material was found to exhibit the greatest degree of variability. The results of these experiments also provide evidence that diffusion from advective zones in fractures through the gouge-filled fractures and recrystallized fracture-filling material could increase the pore space available for matrix diffusion. This evidence is important for understanding the performance of potential nuclear waste repositories in crystalline rocks as diffusion is thought to be an important retardation mechanism for radionuclides.     

Field-scale effective matrix diffusion coefficient for fractured rock: results from literature survey

Matrix diffusion is an important mechanism for solute transport in fractured rock. We recently conducted a literature survey on the effective matrix diffusion coefficient, D_m^e, a key parameter for describing matrix diffusion processes at the field scale. Forty field tracer tests at 15 fractured geologic sites were surveyed and selected for the study, based on data availability and quality. Field-scale D_m^e values were calculated, either directly using data reported in the literature, or by reanalyzing the corresponding field tracer tests. The reanalysis was conducted for the selected tracer tests using analytic or semi-analytic solutions for tracer transport in linear, radial, or interwell flow fields. Surveyed data show that the scale factor of the effective matrix diffusion coefficient (defined as the ratio of D_m^e to the lab-scale matrix diffusion coefficient, D_m, of the same tracer) is generally larger than one, indicating that the effective matrix diffusion coefficient in the field is comparatively larger than the matrix diffusion coefficient at the rock-core scale. This larger value can be attributed to the many mass-transfer processes at different scales in naturally heterogeneous, fractured rock systems.Furthermore, we observed a moderate, on average trend toward systematic increase in the scale factor with observation scale. This trend suggests that the effective matrix diffusion coefficient is likely to be statistically scale-dependent. The scale-factor value ranges from 0.5 to 884 for observation scales from 5 to 2000 m. At a given scale, the scale factor varies by two orders of magnitude, reflecting the influence of differing degrees of fractured rock heterogeneity at different geologic sites. In addition, the surveyed data indicate that field-scale longitudinal dispersivity generally increases with observation scale, which is consistent with previous studies. The scale-dependent field-scale matrix diffusion coefficient (and dispersivity) may have significant implications for assessing long-term, large-scale radionuclide and contaminant transport events in fractured rock, both for nuclear waste disposal and contaminant remediation.     

The effects of matrix diffusion on radionuclide migration in rock column experiments

Rock column experiments were performed to examine the effects of matrix diffusion and hydrodynamic dispersion on the migration of radionuclides at the laboratory scale. Tritiated water and chloride transportation was studied in intact mica gneiss and in altered more porous tonalite columns with narrow flow channels. The column diffusion properties were estimated prior to water flow experiments using the gas diffusion method with helium as the tracer gas. The numerical compartment model for advection and dispersion, with and without matrix diffusion, was used to interpret the tracer transport in the columns. Matrix diffusion behavior was also distinguished from dominating hydrodynamic dispersion in rock column experiments at the slowest water flow rates.     

Migration of colloids in discretely fractured porous media: effect of colloidal matrix diffusion

Numerical simulations of colloid transport in discretely fractured porous media were performed to investigate the importance of matrix diffusion of colloids as well as the filtration and remobilization of colloidal particles in both the fractures and porous matrix. To achieve this objective a finite element numerical code entitled COLDIFF was developed. The processes that COLDIFF takes into account include advective-dispersive transport of colloids, filtration and remobilization of colloidal particles in both fractures and porous matrix, and diffusive interactions of colloids between the fractures and porous matrix. Three sets of simulations were conducted to examine the importance of parameters and processes controlling colloid migration. First, a sensitivity analysis was performed using a porous block containing a single fracture to determine the relative importance of various phenomenological coefficients on colloid transport. The primary result of the analysis showed that the porosity of the matrix and the process of colloid filtration in fractures play important roles in controlling colloid migration. Second, simulations were performed to replicate and examine the results of a laboratory column study using a fractured shale saprolite. Results of this analysis showed that the filtration of colloidal particles in the porous matrix can greatly affect the tailing of colloid concentrations after the colloid source was removed. Finally, field-scale simulations were performed to examine the effect of matrix porosity, fracture filtration and fracture remobilization on long-term colloid concentration and migration distance. The field scale simulations indicated that matrix diffusion and fracture filtration can significantly reduce colloid migration distance. Results of all three analyses indicated that in environments where porosity is relatively high and colloidal particles are small enough to diffuse out of fractures, the characteristics of the porous matrix that affect colloid transport become more important than those of the fracture network. Because the properties of the fracture network tend to have greater uncertainty due to difficulties in their measurement relative to those of the porous matrix, prediction uncertainties associated with colloid transport in discretely fractured porous media may be reduced.     

An interpretation of potential scale dependence of the effective matrix diffusion coefficient

Matrix diffusion is an important process for solute transport in fractured rock, and the matrix diffusion coefficient is a key parameter for describing this process. Previous studies have indicated that the effective matrix diffusion coefficient values, obtained from a large number of field tracer tests, are enhanced in comparison with local values and may increase with test scale. In this study, we have performed numerical experiments to investigate potential mechanisms behind possible scale-dependent behavior. The focus of the experiments is on solute transport in flow paths having geometries consistent with percolation theories and characterized by multiple local flow loops formed mainly by small-scale fractures. The water velocity distribution through a flow path was determined using discrete fracture network flow simulations, and solute transport was calculated using a previously derived impulse-response function and a particle-tracking scheme. Values for effective (or up-scaled) transport parameters were obtained by matching breakthrough curves from numerical experiments with an analytical solution for solute transport along a single fracture. Results indicate that a combination of local flow loops and the associated matrix diffusion process, together with scaling properties in flow path geometry, seems to be an important mechanism causing the observed scale dependence of the effective matrix diffusion coefficient (at a range of scales).     

Plant uptake of radionuclides in lysimeter experiments

The results of seven years lysimeter experiments to determine the uptake of 60Co, 137Cs and 226Ra into agricultural crops (endive, maize, wheat, mustard, sugarbeet, potato, Faba bean, rye grass) are described. The lysimeter consists of twelve monolithic soil profiles (four soil types and three replicates) and is located in Seibersdorf/Austria, a region with a pannonian climate (pronounced differences between hot and semi-arid summers and humid winter conditions, annual mean of precipitation: 517 mm, mean annual temperature: 9.8 degrees C). Besides soil-to-plant transfer factors (TF), fluxes were calculated taking into account biomass production and growth time. Total median values of TF's (dry matter basis) for the three radionuclides decreased from 226Ra (0.068 kg kg(-1)) to 137Cs (0.043 kg kg(-1)) and 60Co (0.018 kg kg(-1)); flux values exhibited the same ranking. The varying physical and chemical properties of the four experimental soils resulted in statistically significant differences in transfer factors or fluxes between the investigated soils for 137Cs and 226Ra, but not for 60Co. Differences in transfer between plant species and plant parts are distinct, with graminaceous species showing, on average, TF values 5.8 and 15 times lower than dicotyledonous species for 137Cs and 60Co, respectively. This pattern was not found for 226Ra. It can be concluded that 137Cs transfer is heavily influenced by soil characteristics, whilst the plant-specific factors are the main source of TF variability for 60Co. The variability of 226Ra transfer originates both from soil properties and plant species behaviour.     

Diffusion of radionuclides in concrete and concrete-bentonite systems

Various construction materials are under consideration for nuclear waste repositories. Two important materials are concrete and bentonite clay, which will act as mechanical barriers and prevent convective water flow. These barriers will also retard transport (diffusion controlled) of dissolved radionuclides by a combination of mechanical constraints and chemical interactions with the solid.An important issue is the possible change of the initial sodium bentonite into the calcium form due to interaction with calcium from the concrete. The initial leaching of concrete was studied using radioactive spiked concrete in contact with compacted bentonite.Measurement were made of the diffusion of Cs, Am and Pu into 5 different types of concrete in contact with pore water. The diffusivity measured for Cs agrees reasonably well with data found in the literature. No movement could be measured for Am and Pu (< 0.2 mm), even though the contact times were extremely long (2.5 and 5 yr, respectively). The diffusion of Na, Ca and Cs from concrete into bentonite was also measured.     

Matrix diffusion coefficients in volcanic rocks at the Nevada test site: influence of matrix porosity, matrix permeability, and fracture coating minerals

Diffusion cell experiments were conducted to measure nonsorbing solute matrix diffusion coefficients in forty-seven different volcanic rock matrix samples from eight different locations (with multiple depth intervals represented at several locations) at the Nevada Test Site. The solutes used in the experiments included bromide, iodide, pentafluorobenzoate (PFBA), and tritiated water ((3)HHO). The porosity and saturated permeability of most of the diffusion cell samples were measured to evaluate the correlation of these two variables with tracer matrix diffusion coefficients divided by the free-water diffusion coefficient (D(m)/D*). To investigate the influence of fracture coating minerals on matrix diffusion, ten of the diffusion cells represented paired samples from the same depth interval in which one sample contained a fracture surface with mineral coatings and the other sample consisted of only pure matrix. The log of (D(m)/D*) was found to be positively correlated with both the matrix porosity and the log of matrix permeability. A multiple linear regression analysis indicated that both parameters contributed significantly to the regression at the 95% confidence level. However, the log of the matrix diffusion coefficient was more highly-correlated with the log of matrix permeability than with matrix porosity, which suggests that matrix diffusion coefficients, like matrix permeabilities, have a greater dependence on the interconnectedness of matrix porosity than on the matrix porosity itself. The regression equation for the volcanic rocks was found to provide satisfactory predictions of log(D(m)/D*) for other types of rocks with similar ranges of matrix porosity and permeability as the volcanic rocks, but it did a poorer job predicting log(D(m)/D*) for rocks with lower porosities and/or permeabilities. The presence of mineral coatings on fracture walls did not appear to have a significant effect on matrix diffusion in the ten paired diffusion cell experiments.     

植物修复--大面积低剂量放射性污染的新治理技术

放射性污染是当今难以治理的环境污染问题,对于大面积低剂量的放射性污染,植物修复是种较适用的新型治理技术.系统阐述了植物修复技术的概念与内容,对放射性污染植物修复技术研究现状与进展、植物修复的优缺点进行了重点评述,并对这技术在未来放射性污染治理中的应用与研究提出了建议.     

Heterogeneous Fenton degradation of persistent organic pollutants using natural chalcopyrite: effect of water matrix and catalytic mechanism

Environmental Science and Pollution Research - Natural chalcopyrite was evaluated as heterogeneous Fenton catalyst. Catalytic performance was evaluated considering different systems, catalyst...     

A review of parameterizations for modelling dry deposition and scavenging of radionuclides

This article aims at reviewing the state-of-the-science parameterizations for modelling dry deposition and scavenging of atmospheric tracers, with a focus on radionuclides. These parameterizations are key components of the numerical models that are used for environmental forecast. We present detailed models and parameterizations. Both are characterized by many uncertainties.     

The deposition of radionuclides from Chernobyl to a forest in Belgium

In the course of a study performed on traditional pollutants, deposition of radionuclides from Chernobyl was observed at a forest in E Belgium. The main deposition flux occurred as a consequence of a 7.4 mm rainfall episode which took place on 4 May. The deposition behaviour of radionuclides is very similar to the behaviour of traditional pollutants with the noteworthy exception of radioactive caesium which seems to have been absorbed by tree foliage.     

The flux of radionuclides in flowback fluid from shale gas exploitation

This study considers the flux of radioactivity in flowback fluid from shale gas development in three areas: the Carboniferous, Bowland Shale, UK; the Silurian Shale, Poland; and the Carboniferous Barnett Shale, USA. The radioactive flux from these basins was estimated, given estimates of the number of wells developed or to be developed, the flowback volume per well and the concentration of K (potassium) and Ra (radium) in the flowback water. For comparative purposes, the range of concentration was itself considered within four scenarios for the concentration range of radioactive measured in each shale gas basin, the groundwater of the each shale gas basin, global groundwater and local surface water. The study found that (i) for the Barnett Shale and the Silurian Shale, Poland, the 1 % exceedance flux in flowback water was between seven and eight times that would be expected from local groundwater. However, for the Bowland Shale, UK, the 1 % exceedance flux (the flux that would only be expected to be exceeded 1 % of the time, i.e. a reasonable worst case scenario) in flowback water was 500 times that expected from local groundwater. (ii) In no scenario was the 1 % exceedance exposure greater than 1 mSv—the allowable annual exposure allowed for in the UK. (iii) The radioactive flux of per energy produced was lower for shale gas than for conventional oil and gas production, nuclear power production and electricity generated through burning coal.     

Modeling the long-term transport and accumulation of radionuclides in the landscape for derivation of dose conversion factors

To evaluate the radiological impact of potential releases to the biosphere from a geological repository for spent nuclear fuel, it is necessary to assess the long-term dynamics of the distribution of radionuclides in the environment. In this paper, we propose an approach for making prognoses of the distribution and fluxes of radionuclides released from the geosphere, in discharges of contaminated groundwater, to an evolving landscape. The biosphere changes during the temperate part (spanning approximately 20,000 years) of an interglacial period are handled by building biosphere models for the projected succession of situations. Radionuclide transport in the landscape is modeled dynamically with a series of interconnected radioecological models of those ecosystem types (sea, lake, running water, mire, agricultural land and forest) that occur at present, and are projected to occur in the future, in a candidate area for a geological repository in Sweden. The transformation between ecosystems is modeled as discrete events occurring every thousand years by substituting one model by another. Examples of predictions of the radionuclide distribution in the landscape are presented for several scenarios with discharge locations varying in time and space. The article also outlines an approach for estimating the exposure of man resulting from all possible reasonable uses of a potentially contaminated landscape, which was used for derivation of Landscape Dose Factors.     

煤矸石及其燃烧产物中天然放射性核素测定及组分分析

以朔州某煤矸石发电厂的煤矸石及其燃烧产物粉煤灰和炉渣为研究对象,采用高纯锗探测器和低本底多道γ能谱仪测定其放射性核素226Ra、232Th和40K,并分析比较了煤矸石中各组分在燃烧过程中的变化情况。通过对测定实验数据进行整理,经过燃烧反应后,煤矸石中的226Ra和232Th向燃烧产物粉煤灰中显著迁移,40K向燃烧产物炉渣中显著迁移。对实验样品进行SEM形貌扫描并作XRD图分析,并给出了煤矸石燃烧后226Ra、232Th和40K的富集因子,分别为3.67、2.08和2.30。且226Ra、232Th在燃烧产物粉煤灰和炉渣中的富集比例均约为7：3,40K的富集比例约为1：1。最后,依据国内外相关标准,分析了该地区煤矸石、粉煤灰和炉渣的放射性水平,对其适用范围进行了评估。     

Management of pikeperch migrating over management areas in a Baltic archipelago area

Pikeperch Sander lucioperca (L.) were tagged in an area where ten small independent management units regulate fisheries, in order to analyze the relevance of migrations for the management. A total of 465 fishes were tagged in connection with the spawning migration. The number of recaptures was 96. The distribution of recaptures between areas and that of total catches, estimated from questionnaire to housholds and, individuals in the survey area, were correlated. The comparison of tagging and recapture dates showed that the same pikeperch migrated to the spawning areas at the same time in two subsequent years, indicating the existence of individual spawning behavior. A modified yield/recruit-model demonstrated that in situations where the dispersal area of a pikeperch stock overlaps with several management areas, the risk of overfishing is very high. Management of this kind of fisheries needs cooperation between units as well as sound data on the movements of pikeperch to define the effects of fishing regulations.     

An axisymmetric diffusion experiment for the determination of diffusion and sorption coefficients of rock samples

Diffusion anisotropy is a critical property in predicting migration of substances in sedimentary formations with very low permeability. The diffusion anisotropy of sedimentary rocks has been evaluated mainly from laboratory diffusion experiments, in which the directional diffusivities are separately estimated by through-diffusion experiments using different rock samples, or concurrently by in-diffusion experiments in which only the tracer profile in a rock block is measured. To estimate the diffusion anisotropy from a single rock sample, this study proposes an axisymmetric diffusion test, in which tracer diffuses between a cylindrical rock sample and a surrounding solution reservoir. The tracer diffusion between the sample and reservoir can be monitored from the reservoir tracer concentrations, and the tracer profile could also be obtained after dismantling the sample. Semi-analytical solutions are derived for tracer concentrations in both the reservoir and sample, accounting for an anisotropic diffusion tensor of rank two as well as the dilution effects from sampling and replacement of reservoir solution. The transient and steady-state analyses were examined experimentally and numerically for different experimental configurations, but without the need for tracer profiling. These experimental configurations are tested for in- and out-diffusion experiments using Koetoi and Wakkanai mudstones and Shirahama sandstone, and are scrutinized by a numerical approach to identify favorable conditions for parameter estimation. The analysis reveals the difficulty in estimating diffusion anisotropy; test configurations are proposed for enhanced identifiability of diffusion anisotropy. Moreover, it is demonstrated that the axisymmetric diffusion test is efficient in obtaining the sorption parameter from both steady-state and transient data, and in determining the effective diffusion coefficient if isotropic diffusion is assumed. Moreover, measuring reservoir concentrations in an axisymmetric diffusion experiment coupled with tracer profiling may be a promising approach to estimate of diffusion anisotropy of sedimentary rocks.