Version française abrégée
1 Introduction
L'Agence nationale pour la gestion des déchets radioactifs (ANDRA) construit un laboratoire de recherche souterrain sur la commune de Bure (Nord-Est de la France), pour étudier la faisabilité d'un stockage profond de déchets radioactifs au sein de la formation argileuse du Callovo-Oxfordien (COX, ). Pour connaître la durabilité chimique des matériaux introduits sur le site (verres, aciers, bétons, barrières ouvragées argileuses...) et la mobilité des radioéléments dans cet environnement, il est nécessaire de connaître la chimie des eaux de la formation. Pour répondre à cet objectif, l'ANDRA a carotté un grand nombre d'échantillons au niveau du laboratoire et dans son environnement régional. Ces carottes ont permis d'obtenir un grand nombre de données géochimiques utilisables pour définir la chimie des eaux interstitielles de la formation, à une échelle régionale. La littérature scientifique dans ce domaine [2,26] montre que la définition de cette chimie est très difficile. L'extraction de l'eau interstitielle par pressage [7,8], ou l'échantillonnage des eaux en forages spéciaux à chambre sous N2 [17,29,32] génèrent des artéfacts modifiant cette chimie dans des proportions variables, l'oxydation des échantillons étant l'artéfact le plus difficile à maîtriser [6,8]. Dans le cadre du partenariat BRGM–ANDRA, une méthode a été développée, d'abord au sein du laboratoire de recherche du mont Terri [24], puis appliquée à Bure. À partir d'un jeu complet de données géochimiques, il est possible de construire des modèles d'équilibre eau/argilite sans ajustement de paramètres. Dans les modèles antérieurs, Bradbury et Baeyens [4] étaient obligés d'estimer la , Beaucaire et al. [3] négligeaient les argiles 2:1 et les réactions d'échange d'ions, ou encore Motellier et al. [21] faisaient l'impasse sur les teneurs en sulfates.
2 Situation géologique
La zone étudiée par l'ANDRA [13] concerne un bloc structural homogène ceinturé par des accidents structuraux (Fig. 1). Le COX imperméable est recouvert par l'Oxfordien calcaire et le Kimméridgien, qui présentent quelques niveaux producteurs en eau. Il surmonte le Dogger peu productif en eau, le Lias, qui est imperméable et les formations triasiques qui contiennent des saumures.
3 Acquisition des données expérimentales
Les paramètres utilisés pour construire les modèles sont les suivants : (1) une identification exhaustive des minéraux et une qualification de leur état cristallin (Tableau 1) ; (2) une qualification de l'état réduit du système, compte tenu de la présence de pyrite et d'une capacité réductrice de la roche en moyenne de 4,7 meq g−1 (milliéquivalents d'électrons par gramme de solide sec) ; (3) une mesure de la teneur en eau (Tableau 2) ; (4) une mesure des surfaces spécifiques mouillées par l'eau, pour calculer la répartition entre eau libre et eau liée tenant compte des phénomènes d'exclusion anionique [18,23] existant dans les bancs argileux de la formation (Tableau 3) ; (5) des expériences de lessivage pour connaître les teneurs en anions dissous (Tableau 4, Fig. 2) ; (6) des mesures de cations adsorbés (Fig. 4) ; (7) l'acquisition des constantes d'échange ionique spécifiques de l'argilite (Tableau 6, Fig. 5) ; (8) la mesure des pressions partielles en CO2 (Fig. 6) ; (9) la mesure des du CO2 (Fig. 7). Les protocoles expérimentaux correspondant peuvent être trouvés dans les références [12–15].
Synthesis of the mineralogical data gathered for the Callovian–Oxfordian formation
Synthèse des données minéralogiques obtenues pour la formation du Callovo-Oxfordien
Mineral | Abundance | Depth at the URL location | Crystallinity | Insertion in the modeling |
Illite | +++ | ever present | good, Fe-rich | yes |
Interlayered I/SR0 | +++ | >493.4 m | good | yes |
Interlayered I/SR1 | +++ | <489.2 m | good | yes |
Glauconite | trace | no systematic | ? | no |
Kaolinite | + | <453.7 m | good | can be tested |
Montmorillonite | ++ | special level at 425.03–425.08 | good | can be tested |
Chlorite (chamosite) | trace | ever present | various Fe-content | yes |
Biotite | trace | no systematic | chloritised | no |
Muscovite/phengite | trace | no systematic | detrital | no |
Microcline/orthoclase | trace | no systematic | good, detrital | can be tested |
Albite/plagioclase | very rare | very rare | detrital, alteration figures | no |
Quartz | +++ | ever present | good, detrital figures | yes |
Chalcedony | trace | no systematic | good | can be tested |
Calcite | +++ | ever present | good, several crystallizations | yes |
Aragonite | trace | no systematic | good | no |
Ca-rich dolomite | ++ | ever present | good | yes |
Ankerite/siderite | ++ | ever present | good | yes |
Celestite | trace to ++ | ever present | good | yes |
Pyrite | + | ever present | polymorph, good | yes |
Sphalerite/chalcopyrite | very rare | very rare | ? | no |
Native sulfur | very rare | very rare | very small particle | no |
Rutile/anatase | trace | ever present | good | no |
Apatite/collophane | trace | no systematic | ? | no |
Mean values of mass water content (%) in the Callovian–Oxfordian formation at the regional scale
Valeurs moyennes (%) des teneurs massiques en eau de la formation du Callovo-Oxfordien à l'échelle régionale
Borehole | Water content (%) | Standard error (%) |
EST312 NE | 8.04 | 1.34 |
EST205 URL | 7.03 | 1.80 |
EST212 URL | 6.62 | 0.92 |
EST207 URL | 6.73 | 1.00 |
EST322 SW | 5.54 | 2.09 |
EST342 SW | 5.57 | 1.98 |
Results of H2O specific surface and calculation of the % of free water
Résultats de surface spécifique à l'eau et calcul du % d'eau libre
Est205 Sample & depth (m) | H2O Specific Surface (m2/g) H2O size: 10.6 Å | Total water content (g/100 g) (dry base) | Free water (g/100 g) (dry base) | Free water (%) |
K100 424m | 30.0 | 6.13 | 3.79 | 61.8 |
K108 447m | 35.6 | 6.28 | 3.50 | 55.8 |
K113 462m | 36.0 | 6.18 | 3.34 | 54.1 |
K119 477m | 56.1 | 8.98 | 4.60 | 51.3 |
Concentrations of Cl−, Br−, in the total of the porosity (TW) and in the geochemical porosity (GP) calculated from leaching experiments
Concentrations de Cl−, Br−, dans la porosité totale (TW) ou dans la porosité géochimique (GP) calculées à partir des expériences de lessivage
EST205 Sample | Depth (m) | Cl− (mol l−1) | Br− (mol l−1) | ||||
TW | GP | TW | GP | TW | GP | ||
K100 424m | 423.68−424.48 | 1.79×10−2 | 2.90×10−2 | 3.2×10−5 | 5.1×10−5 | 5.17×10−2 | 8.35×10−2 |
K119 477m | 476.57−476.97 | 1.54×10−2 | 3.01×10−2 | 3.4×10−5 | 6.5×10−5 | 3.22×10−2 | 6.28×10−2 |
Callovian–Oxfordian formation (fraction <2 μm). Gaines–Thomas constants for a global type of exchange sites
Formation du Callovo-Oxfordien (fraction <2 μm). Constantes de Gaines–Thomas pour un type global de sites d'échange
Model 1 type of sites | Constant | Standard deviation | |
H/Na | 13.0 | 2.63 | 1.11 |
K/Na | 2.78 | 0.12 | 0.44 |
NH4/Na | 2.98 | 0.14 | 0.47 |
Ca/Na | 3.17 | 0.10 | 0.50 |
Mg/Na | 2.41 | 0.09 | 0.38 |
Sr/Na | 4.75 | 0.19 | 0.68 |
4 Résultats et modélisation
Les observations montrent une grande homogénéité dans la nature des minéraux rencontrés à l'échelle régionale. Les constantes d'échange ionique correspondent à celles attendues pour des argiles de types illite et smectite [1,9,16,27]. La composition isotopique du CO2 montre que la concentration de ce gaz dans l'eau de la formation est régulée de manière constante par l'équilibre entre la solution et les minéraux de la roche, comme cela est prévu par Coudrain-Ribstein et al. [5]. En revanche, les teneurs en Cl− chutent de 100 à 15 mmol l−1 en suivant un axe NE–SW (du forage EST312 au forage EST322) (Fig. 3). Cette variation s'accompagne d'une évolution des cations adsorbés (Fig. 4). Compte tenu du profil de concentration en Cl− du Dogger, on peut émettre l'hypothèse d'une entrée de solutions salées provenant du Trias par un système de failles situé au nord-est du forage EST312 et diffusant horizontalement au sein du Dogger (les gradients hydrauliques, assez mal connus, paraissent très faibles dans cette formation). Le Dogger influence nettement la chimie des eaux du COX.
Le modèle d'équilibre thermodynamique entre la solution, la minéralogie, les surfaces argileuses, et la utilisant le code PHREEQC [22] et la base de données llnl.dat [33] a permis de calculer plusieurs compositions d'eau, compte tenu d'hypothèses sur les valeurs de mesurées. Le modèle de Coudrain-Ribstein et al. [5], ne fixant pas la , mais régulant ce paramètre par un assemblage de minéraux en équilibre avec la solution, a été testé et donne de bons résultats (Tableau 7). Appliquée à une échelle régionale, la modélisation a permis de calculer la variabilité de la composition chimique des eaux (Fig. 8).
Model results (PHREEQC, llnl.dat). Various test cases of porewater chemistry: Cl− concentration in the geochemical porosity, ‘A’ simulation use the maximum measured value and ‘B’ a value in the best analytical range, ‘Chlo’ corresponds to test cases when the is not fixed. In this model, is determined by the equilibrium between calcite, dolomite, and Mg-chlorite
Résultats de modélisations (PHREEQC, llnl.dat) : Différents cas tests pour la chimie de l'eau interstitielle : concentration en Cl− dans la porosité géochimique, la simulation ‘A’ utilise la valeur mesurée de maximum, ‘B’ utilise la valeur de obtenue dans les meilleures conditions analytiques, ‘Chlo’ correspond à un cas test, où la n'est pas fixée. Dans ce dernier modèle, la est déterminée par équilibre avec la calcite, la dolomite et une chlorite-Mg
Model | A | B | Chlo |
Porewater chemistry | mol/kgw | mol/kgw | mol/kgw |
pH | 7.10 | 7.28 | 7.00 |
−163 | −176 | −156 | |
Alk | 1.97×10−3 | 1.29×10−3 | 2.46×10−3 |
Ionic strength | 1.16×10−1 | 1.16×10−1 | 1.17×10−1 |
Cl | 3.01×10−2 | 3.01×10−2 | 3.01×10−2 |
S | 3.39×10−2 | 3.40×10−2 | 3.39×10−2 |
Na | 3.21×10−2 | 3.20×10−2 | 3.21×10−2 |
K | 7.09×10−3 | 7.07×10−3 | 7.10×10−3 |
Ca | 1.49×10−2 | 1.48×10−2 | 1.50×10−2 |
Mg | 1.41×10−2 | 1.41×10−2 | 1.42×10−2 |
Sr | 1.12×10−3 | 1.12×10−3 | 1.12×10−3 |
Si | 9.41×10−5 | 9.43×10−5 | 9.40×10−5 |
Al | 7.36×10−9 | 8.62×10−9 | 6.85×10−9 |
Fe | 2.14×10−4 | 9.40×10−5 | 3.32×10−4 |
saturation index (si) | |||
−1.96 | |||
si illite | 0 | 0 | 0 |
si montmor-Ca | −0.177 | −0.143 | −0.195 |
si microcline | 0.5581 | 0.652 | 0.509 |
si calcite | 0 | 0 | 0 |
si dolomite | 0 | 0 | 0 |
si Fe-chlorite | 0 | 0 | 0 |
si quartz | 0 | 0 | 0 |
si siderite | 0.0906 | −0.235 | 0.261 |
si pyrite | 0 | 0 | 0 |
si celestite | 0 | 0 | 0 |
si gypsum | −0.0085 | −0.0086 | −0.0085 |
si Mg-chlorite | 0.85 | 2.48 | 0 |
ion exchange | |||
SrX2 | 4% | 4% | 4% |
KX | 8% | 8% | 8% |
NaX | 14% | 14% | 14% |
MgX2 | 28% | 28% | 28% |
CaX2 | 45% | 45% | 45% |
5 Conclusion
La méthode employée permet de définir la chimie des eaux sans estimation de paramètres. Cependant, différents modèles peuvent être testés compte tenu des informations géochimiques recueillies. Cela permet de donner une gamme de compositions relativement restreinte pour la chimie des eaux au niveau du laboratoire souterrain de Bure. Les variations de teneurs en Na et Cl suivant un axe NE–SW au sein de la formation sont imputables au Dogger. Il conviendrait de construire un modèle 2D rendant compte de cette dynamique de transport.
1 Introduction
1.1 Context of the study
The French Agency for Nuclear Waste Management (ANDRA) is building an underground research laboratory (URL) in Bure (northeastern France) to study the disposal of radioactive waste in a deep clay formation. The main gallery is located in the Callovian–Oxfordian formation (COX) at an approximate depth of 490 m. To assess the durability of materials (glass, stainless steel, concrete, engineered clay barrier, etc.) and determine the mobility of the radionuclides in the host rock formation, it is necessary to define the chemistry of the COX porewater. ANDRA has therefore drilled a series of boreholes in the vicinity of the URL (Fig. 1). Core drilling provides a vast amount of geochemical data. This paper presents a method enabling the modelling of the porewater chemistry. This method is applied at a regional scale and enables a proper description of the chemical conditions in the formation.
1.2 Studying the geochemistry of clay systems
Studying the water chemistry of an argillaceous formation is a challenging task. The recovery of undisturbed water samples by classical methods is almost impossible due to the low permeability and low water content. Recently developed techniques for water extraction (squeezing, leaching, special borehole with gas control) give valuable results that can be used to characterize the water chemistry. However, each one of these techniques introduces experimental and analytical artefacts. A review of studies of the porewater chemistry of mudrocks and extraction techniques can be found in [2,26]. The chemistry of porewater in clay systems is difficult to define for several reasons, including the very low porewater content and its strong interactions with clay minerals. The physical interaction between water and minerals means that either high pressure must be used to recover water by squeezing, or long-term sampling must be done to recover water in special boreholes in URLs. The high pressures used for squeezing techniques can dissolve some minerals. Even for conservative anions like Cl, squeezing techniques generate artefacts [7,8]. Sampling in special boreholes requires preventive measures. As the porewater slowly fills the borehole chamber (specifically drilled and/or stored under N2 to minimize the oxidation) [17,29], it will degas and bacterial development can perturb the system [32]. The reducing conditions that prevail in such systems could be another major source of concerns. Oxidation might significantly modify the mineralogy and the porewater quality [6,8]. On the contrary, growing of anaerobic bacteria can strongly reduce the system [32].
1.3 Development of the THERMOAR method
In a preliminary step to the URL porewater recovering experiments, a strategy to characterize and model the geochemical system was developed within the framework of the THERMOAR project through an ANDRA–BRGM partnership. The method was first developed for Opalinus Clay (Switzerland, [24]). The first samples from Bure came from the EST205 borehole drilled on the URL site where the laboratory's auxiliary shaft has later been dug. Many additional samples from regional boreholes drilled later by ANDRA have been added to the original dataset. Most of the geochemical data obtained by ANDRA on the formation have been used to define the porewater chemistry at a regional scale.
The construction of the thermodynamic model enables us to define the porewater chemistry and also to describe reaction pathways. With this tool, we are able to model perturbations (alkaline solution, oxidation, microbial disturbance, etc.). In the THERMOAR method, data for a complete set of physical and chemical parameters were obtained to constrain the model. The originality of this approach is that no parameters are estimated and all major mineral phases are taken into account. In the past, clay porewater models have generally prescribed a hypothetical CO2 partial pressure (e.g., the model developed by Bradbury and Baeyens for the Opalinus Clay formation (Switzerland) [4]). Beaucaire et al. [3] disregard 2:1 clay and ion exchange reactions in a model for the Boom Clay porewater. Another model for the Callovian–Oxfordian formation developed by Motellier et al. [21] includes ion exchanges for three types of sorption sites, but calculates the aqueous speciation partially and does not take into consideration the content.
After a short description of the geology of the area, the data are presented and used to constrain the models. As several options are possible, the results of a series of models are presented and the most realistic model for the URL was applied to the evolving chemical conditions observed northeast and southwest of the URL site.
2 Geological setting
The geology of the site is described in detail in [13]. The eastern Paris Basin sedimentary formations slightly dip (1 to 2°) to the west within a simple monocline related to the subsidence of the central part of the basin. Except for major regional fault zones (the Bray, Vittel and Metz faults) and the Gondrecourt and Marne faults delimiting the area investigated by ANDRA (Fig. 1), there is very little evidence of deformation.
At Bure, the Callovian–Oxfordian (COX) argillite (conventionally referred to as C2) is about 130 m thick. Five lithostratigraphic units have been identified and labelled C2a, C2b1, C2b2, C2c and C2d, from base to top. The sedimentary series is primarily clayey in the bottom two thirds, and becomes progressively siltier, and then increasingly calcareous towards the top. A remarkable sedimentary surface is associated with the maximum clay zone (MCZ, 489–486 m): the maximum flooding surface (MFS) of the Callovian–Oxfordian cycle. This cycle marks the transition from transgression at the base to a regression, which led to the installation of the Oxfordian carbonate platforms. The top of the clayey formation is located at a depth of around 417 m, overlapped by 243 m of Middle and Late Oxfordian limestone, followed by 175 m of Kimmeridgian marl and limestone, capped by the Barrois limestone of Kimmeridgian age. The clayey formation is separated from the Dogger limestone by a regional palaeo-erosion surface. Under the Dogger limestone or marl (Bathonian, Bajocian, Aalenian), the Lias is marl-rich and capped the Triassic formations (sandstone and evaporite).
Some levels in the Kimmeridgian, Oxfordian and Dogger limestone are permeable. In the last two formations, water productivity is low. The hydraulic gradient in the Dogger is known as low. The COX is considered as impermeable. The Lias is impermeable or an aquitard formation. The Triassic is known to contain brines.
3 Materials and methods
3.1 Origin of the core samples
The location of the boreholes is shown in Fig. 1. Boreholes EST205, EST207, EST212 concern the URL and additional data were obtained on samples from the laboratory auxiliary access shaft (PAX). The distance from the URL to boreholes MSE101, EST312, EST322 is 12.5, 15.5 and 9 km, respectively. They are located, respectively, northwest, northeast and southwest of the URL. All of these boreholes and shafts are located in the same structural block. Borehole EST342 is located beyond the ‘Marne graben’ and the ‘Joinville graben’ structural accidents. The depths of boreholes in Figs. 3 and 6 are normalized to depths corresponding to borehole EST205 (URL site). As the thicknesses are not identical at all locations, the correspondence of depths between boreholes is not perfect.
Borehole EST205 was drilled with an oil-based mud consisting of 80% (in volume) oil (gas–oil containing 0.4% aromatics) and 20% brine (CaCl2, 350 g l−1). The oil and brine mixture forms an emulsion whose high Ca content prevents swelling of the clay minerals. The addition of CaCO3 to the mud helps to obtain densities above 1.1, necessary for coring in clay formations. This type of drilling mud helped prevent the mechanical evolution of the core samples and the swelling of clay minerals. However, the brine imposed a very high Cl gradient, which contaminated the cores by diffusion. An immediate treatment of the cores eliminates the edges, which preserves the chemistry of the samples. This was done for a limited number of borehole EST205 samples. The authors do not recommend the use of Cl data published in [13] which are too high because of the contamination by the oil-based mud. (The same problem exists for borehole EST211.) The samples from the shafts were air drilled without mud. The mud used for boreholes EST212, EST207, EST342, EST312, EST322, MSE101 was a common mixture of bentonite and polymers. No major contamination was detected for this last series of boreholes.
3.2 Physical and geochemical parameters
The parameters used in the model are the following:
- • mineral identification and crystallinity of minerals;
- • total reduced capacity of the rock (TRC);
- • water content and argillite texture;
- • anion leaching;
- • cation-exchange capacity (CEC) and adsorbed cation population;
- • ion-exchange constants (, , , , );
- • CO2(g) partial pressure and isotopic signature of the gas.
The analytical techniques used to identify and characterize the minerals, to measure the TRC, to determine the water content and the textural properties, the leaching of anions, the cation exchange capacity and the adsorbed cation population are detailed in [13]. The experimental protocol used to determine ion exchange constants has been published in [12]. The method used to obtain the and the isotopic signature of the gas is presented in [15].
4 Experimental and modelling
4.1 Experimental data
4.1.1 Identification and crystallinity of minerals
An accurate modelling of the water–rock interactions takes into account all minerals present, even at trace levels, and requires a complete mineralogical analysis of the samples. Because several minerals can control the aqueous speciation of an element, it was also necessary to determine some crystallinity criteria. Indeed, any mineral having alteration figures or an unbalanced structural formula is not considered at equilibrium in the model. Table 1 lists all of the minerals found in the Callovian–Oxfordian formation with their relative abundance, distribution, and information on their crystallinity. Table 1 also indicates whether the phase must be taken into account in the thermodynamic model.
4.1.2 Total reduced capacity of the rock (TRC)
The TRC is an index of the redox potential of the rock. It represents its ‘oxidizability’ and is expressed in term of mole of available electrons. The mean value for the formation is . These values are similar to those measured in the Opalinus Clay Formation, [30]. TRC values do not vary a lot in the clay formation. The lower values are mainly due to calcareous beds (1.6 mequiv g−1). TRC is linked to the organic matter content and mainly to sulphide minerals and to a lesser extent to siderite. Conditions in the formation are clearly reducing. In the model, the redox potential can be imposed by pyrite.
4.1.3 Water content and argillite texture
The Callovian–Oxfordian formation has a very low water content. Water content ranges from 6% to 9% (mass water content dry base). Low contents (1.5%) correspond to levels rich in fine-grained carbonate cement. There is a very good correlation between high carbonate content () and very low water content (). At the regional scale, variations in water content seem to be insignificant (Table 2).
The water in the clay is of four different types [18]:
- • interlayer water within the sheet structure of clay minerals. This water is strongly adsorbed;
- • another strongly adsorbed water at the clay mineral surface is the (structured) water participating in the hydration shells of the cations (which locally balance the mineral surface charges);
- • loosely held water in the diffuse-layer;
- • free water.
For geochemical purposes, it is convenient to distinguish only between ‘bound’ water (which includes interlayer, strongly adsorbed and some diffuse layer waters) and ‘free’ water (which includes, not only Horseman's free water [18], but also some of the diffuse layer and surface-sorbed waters as well) [23]. Note that the concept of geochemical porosity has been developed for clay to take into account anionic repulsion due to the negatively charged surface of the clay. Anionic repulsion limits the access of anions to part of the porosity. Geochemical porosity is the porosity that is accessible to mobile anions (such as Cl or Br, for example). In a carbonaceous rock, anionic repulsion is limited and anions have access to all of the porosity. This distinction between geochemical porosity and total porosity for argillaceous rock is approximate, however practical. Thus, it is possible to quantify the two types of water by measuring the specific wetted surface using water sorption isotherms () [10]. For the argillite, the ‘bound’ water is estimated to be about 2 to 3 water layers thick. Knowing the and the thickness of the water layer, the volume of ‘bound’ water is estimated to be about 50% of the total water volume (Table 3). The ionic concentrations can be calculated either with the total water content (TW) or only with the ‘geochemical porosity’ (GP). The advantages of the ‘geochemical porosity’ concept will be demonstrated in the next paragraph.
4.1.4 Leaching of anions
Leaching of ions is a classic technique to obtain Cl− and Br− (and, in some cases, I− and ) contents in clay porewater. Information concerning cations, aqueous silica, carbonates and fluorides is generally disturbed by the dissolution of minerals. When plotted, concentrations of leached conservative anions give straight lines as a function of the solid/liquid ratio of the batch systems. The dissolution of a mineral gives a hyperbole (Fig. 2), and the oxidation of rock samples produces scattered values.
The behaviour of during the leaching experiments with a regression coefficient of 0.997 shows a perfect conservative character. The experiments seem to be undisturbed by pyrite oxidation. The Sr concentration value indicates that celestite is not significantly dissolved during leaching. The concentrations obtained by leaching represent probably the maximum values present in the porewater. Examples of concentrations of Cl−, Br−, in the total water (TW) and in the geochemical porosity (GP) are given in Table 4 for borehole EST205 (external part of the samples immediately cut after drilling to avoid contamination).
An interesting observation was made during the calculation of Cl− concentrations in TW and GP. If we know the Cl− concentration in mg kg−1 of rock, we can calculate the Cl− concentration in the porewater. For the samples from the auxiliary shaft, the mean concentration of Cl− in the total porosity is . If a distinction is made between the argillaceous layers and the carbonaceous layers (water content ), two mean TW concentration values are observed (Table 5). The Cl− concentration in the carbonaceous layers is about twice as great as in the argillaceous layers. If we consider that the concept of geochemical porosity is applicable only to the argillaceous layers, the Cl− concentration values in the two media are very similar. The Cl− content in the formation is rather homogeneous at the URL site and a value of can be used in the model.
Comparison of calculated concentrations of Cl in the calcareous and argillaceous layers as a function of the accessible porosity
Comparaison des concentrations en Cl calculées au sein des bancs argileux ou calcaires en fonction de l'accessibilité de la porosité
Rock | Anions accessible Porosity (geochemical porosity concept) | Cl− concentration in the total porosity (TW) (mol l−1) | Cl− concentration in the geochemical porosity (GP) (mol l−1) |
Carbonates | 100% | 3.48×10−2±1.96×10−2 | 3.48×10−2±1.96×10−2 |
Argillites | 50% | 1.65×10−2±5.80×10−3 | 3.61×10−2±1.68×10−2 |
The Cl− contents were determined for all of the boreholes (Fig. 3). Cl− concentrations were calculated in the geochemical porosity for the argillaceous layers and in the total porosity for the calcareous layers. Dogger porewater was sampled in the regional boreholes and the Cl− concentrations are given in Fig. 3, which shows the strong influence of this formation on the Cl− content in the Callovian formation. The Cl− content varies with depth in borehole EST312. A straight line can be drawn from the Dogger at 110 mmol l−1 to the top of the formation at 40 mmol l−1. For borehole MSE101 and at the URL site, there is less vertical variation between 50–60 and 30 mmol l−1. Although MSE101 is located northeast of the URL, no Cl− variations can be detected between the two locations. For the two boreholes located southwest of the site, the vertical variation is negligible (EST 322: , EST342: ). There is, therefore, a regional variation with a decreasing concentration of Cl− from northeast to southwest.
One possible explanation for this might be that saline water migrates from the Triassic formation via a fault system located northeast of borehole EST312. Since there is a Cl− concentration gradient in the Dogger formation, diffusion or a very slow advection might explain the horizontal migration of Cl−. The Dogger clearly influences the chemistry of porewater in the Callovian–Oxfordian formation.
4.1.5 Cation exchange capacity (CEC) and adsorbed cationic population
At the URL (EST205/PAX data), the mean CEC of the Callovian–Oxfordian formation is . The lowest values are measured in the more calcareous beds (2.1 meq/100 g) and the highest values in the zone (25.1 meq/100 g at 463.74 m). The (less rich in smectite) has a very stable average value of .
As Bradbury and Baeyens [4] pointed out, the adsorbed cation population is a reflection of the solution in equilibrium with the mineral assemblage. Exchange reactions are very rapid and adsorbed cations are immediately in equilibrium with the porewater. The proportion of exchangeable cations is, therefore, very useful information. At the URL site, the distribution of exchangeable cations slowly evolves. In addition to the 4% of exchangeable Sr2+, the adsorption sites are taken up mainly by of Ca, of Mg, of Na, and of K.
Exchangeable Sr was not systematically determined at the regional scale. Only the concentrations of exchangeable K, Na, Mg, and Ca will, therefore, be discussed. No significant vertical variations were observed for the different boreholes. This is also true for EST312, which presents the major Cl evolution. From northeast to southwest, there are large variations in the exchangeable-cation distribution (Fig. 4). CaX2 values are identical in borehole EST312 and in the URL (EST205/PAX), and MgX2 increases to the detriment of alkaline ions. For example, there is only half as much KX. Southwest of the URL site at borehole EST322, the proportion of MgX2 is identical to that at the URL, there is a greater proportion of CaX2 and only one third as much NaX. Borehole EST342, located beyond the Marne and Joinville fault systems, has a greater proportion of NaX, to the detriment of MgX2.
Considering the ion exchange populations of the regional boreholes, we can conclude that, for the COX porewater:
- • the proportion of Na decreases from northeast to southwest following the decrease of the Cl− content;
- • the proportion of K also decreases, in a lesser extent;
- • Ca and Mg concentrations are stable considering the electroneutrality.
4.1.6 Ion-exchange constants (KH/Na, KK/Na, KGT NH4/Na, KCa/Na, KMg/Na, KSr/Na)
The work of many authors [1,4,16,21,24,25,28,31] shows that there is a large consensus that ion-exchange reactions are essential if we wish to model the porewater chemistry of clay-rich formations. We therefore determined a complete set of ion-exchange constants (, , , , , ) on the <2-μm clay fractions using the Gaines–Thomas [11] approach. As opposed to [19,21], which use H+, our reference cation is Na+. This gives constants that can be compared to data presented in the other studies after formalism adaptation [9,27].
For a detailed description of the following experimental method, see [12].
Batch experiments were done using the purified clay fraction () from a sample (K119 at 476 m) from borehole EST205. The clay (0.5 g) was equilibrated (72 h) with a solution (25 ml) of NaxA1−xCl (A = K or ) or NaxB0.5−xCl (B = Ca, Mg, Sr). The discrete variations of x enabled a progression of the ion-exchange reactions. After separation (centrifuging and filtration), this first solution was analyzed and the equivalent fractions of aqueous species were determined. The clay was re-dispersed (24 h) in a concentrated solution of cobaltihexamine (4 g l−1). After separation (centrifuging and filtration), this second solution was analyzed and the equivalent fractions of adsorbed species were calculated. During this calculation, the molar quantity of aqueous species remaining in the residual porewater (after the first centrifugation) was corrected. The equivalent fractions of aqueous and adsorbed Na and A (or B) were used to construct isotherms (Fig. 5).
Using the Gaines–Thomas approach [11], it is possible to calculate for the exchange sites , the selectivity coefficient (Eq. (1)) for an exchange between two cations M and N of charges and , respectively. and (Eq. (2)) are the equivalent fractions of the adsorbed species. In Eq. (2), is the exchange capacity of the sites.
(1) |
(2) |
Isotherm data are processed from Eq. (1) using a non-linear regression [20]. By modelling the Na/K isotherm, we were able to identify illite-type and smectite-type sites. For the Opalinus Clay, a similar fit was done with two types of sites [12]. However, in this study, the global approach considering the ion exchange sites as a whole was tested for the modelling of porewater in neutral conditions. Therefore, in this paper, only the global constants are given in Table 6.
4.1.7 CO2(g) partial pressure and isotopic signature of the gas
The most original component of the THERMOAR project is the direct analysis of CO2 partial pressure on core samples [14]. The core samples were placed in a rigid container immediately after they had been drilled and the edges sawed off. Then the container was filled with an inert gas (N2, He, or Ar). The results are slightly different, depending on the inert gas used. The best results were obtained when the initial total pressure in the container was sufficiently high . The temperature and pressure in the cell were monitored continuously. The container was placed in a ventilated cupboard thermostated at a temperature of 25 °C. The gases released from the container were analyzed by gas chromatography until the stabilized. At this point, analyses of CO2 isotopes were done [15].
The mean values for the reference borehole EST205 are between 3.10 and 7.17 mbar (). These results are in the range of data compiled by Coudrain-Ribstein et al. [5] for sedimentary confined aquifers (1 to 100 mbar at 25 °C). Fig. 6 presents the data obtained at a regional scale in the best analytical conditions (). Values range between 1.4 and 4.1 mbar . No variation as a function of depth or at the regional scale can be observed. This can be a good indication of an internal control of this parameter by the mineralogy, the nature of which is constant throughout the formation. The isotopic composition of the gas confirms this. It shows (Fig. 7) that the CO2 has an original isotopic composition independent of the mantle, the atmosphere, organic matter or bacterial activity. Considering the gas values, we were able to determine the isotopic composition of the bicarbonate ions at equilibrium with the gas. The values are between −2 and . The signatures of the bicarbonates and of the carbonate cements of the clay rocks (mainly calcite and some dolomite and siderite) are similar. This proves that the is controlled by the mineralogical assemblage (carbonates) and not by degradation of organic matter or diffusion of gas coming from other geological formations.
4.2 Modelling porewater chemistry at a regional scale
4.2.1 Introduction
In the Callovian–Oxfordian formation, the nature of the mineralogy and the regulation mechanisms can be considered as invariable at a regional scale. As the nature of the clay is also stable, the exchange selectivity coefficients can be used for the entire formation. The regional variation of the porewater chemistry was observed for Cl− and the ion exchange population. The model was initially constructed to define the porewater chemistry at the URL and, in a second step, was adapted to regional variations.
4.2.2 Constraints to model the porewater at the URL
Due to the stability of the mineralogy and the geochemical data, the porewater chemistry was considered to be constant in the clay formation at the URL location. All calculations were done at thermodynamic equilibrium. Due to the very low permeability of the clay formation and the probable residence time of the water, it is realistic to suppose that thermodynamic equilibrium can be reached in these geological media. The numerical model of the porewater chemistry involves two steps. First, seawater is diluted by pure water, and then the resulting water is equilibrated with the mineralogical assemblage and the exchange sites (X) of the clay. The dilution is defined so as to obtain the Cl− concentration determined in Section 4.1.4. The entire mineral assemblage (Section 4.1.1) controls the following elements: e−, S, C, Si, Ca, Mg, Na, K, Sr, Fe, Al. Concerning the mineralogy, some simplifications have been considered: pure phases (dolomite and siderite) and end-members for chlorite solid solution (daphnite [Fe] and clinochlore [Mg]) have been introduced in the modelling, in place of the observed minerals. Some geochemical subsystems can be distinguished in the model. Na is only determined by ion exchange (NaX). A subsystem Fe, S, e−, Sr, accounts for daphnite (Fe-chlorite), celestite, pyrite, and ion exchange (SrX2). A subsystem Si, Al, K is determined by quartz, illite, and ion exchange (KX). A subsystem C, Ca, Mg, is controlled by , calcite, and dolomite. The ion exchange equations for Ca and Mg are considered in the equation system. However, they are not the determining factors for the water chemistry. They can only adjust the distribution of cations on the exchanger. The pH is calculated by iteration to obtain the electroneutrality of the solution.
The calculations are performed with the PHREEQC code [22] using the llnl.dat database [33]. Three cases are presented in Table 7. The simulations use the Cl− concentrations determined in the geochemical porosity. ‘A’ simulation uses the maximum measured value and ‘B’ simulation a value in the best analytical range (Section 4.1.7). ‘Chlo’ corresponds to test cases where is not fixed. Mg-chlorite was added to the phase assemblage as in the model of Coudrain-Ribstein et al. [5]. In this model, is not controlled externally and is determined by the equilibrium between calcite, dolomite, and the Mg-chlorite.
4.2.3 Results: Composition of the porewater at the URL
Fig. 8 shows that the solutions are relatively concentrated, with an ionic strength of . The pH values are between 7.00 and 7.28. The values are low . The waters have a Cl/SO4 facies. Na is twice as concentrated as Ca and Mg and four times as concentrated as K. The alkalinity is moderate, between 1 and 3 mmol kg−1w. The Al and Si contents agree with a system controlled by silicates and aluminosilicates. The Fe contents (determined by the equilibrium with daphnite) are quite high. However, the model calculates a quasi-equilibrium with siderite, which has been found in places in the formation. The value calculated with the ‘Chlo’ model is similar to the maximum measured values . This is a good indication that the model converges toward a representative chemistry. Another possible verification is the proportion of adsorbed cations that are accounted for by model. Table 7 shows that the results are close to the measured values. On the contrary, if microcline is tested in the model instead of illite, the percentage of adsorbed cations is significantly different. This indicates that this phase does not participate in the regulation of the porewater chemistry.
4.2.4 Modelling porewater at the regional scale
To model the porewater at a regional scale, the ‘B’ model described in Section 4.2.2 was adjusted using the Cl− concentrations obtained for boreholes EST312 and EST322 (Section 4.1.4). The corresponding adsorbed ion concentrations are considered in the model (Section 4.1.5). Fig. 8 shows the evolution of the porewater chemistry, which is related mainly to a decrease in salinity (NaCl) from northeast to southwest in the structural block in which the URL is located. The decrease is of about one order of magnitude for the Cl and Na concentrations. The other concentrations are determined by the mineralogy and do not evolved significantly.
5 Conclusions
The THERMOAR method enables a detailed characterization of the geochemistry of the COX. The dataset is complete and a model can therefore be developed without any estimated parameters. The different models presented in this paper for the ANDRA URL give an indication of the range of variation that can be calculated for the porewater chemistry. At the regional scale, the Cl− contents change as well as the adsorbed cations. It was also possible to model the porewater chemistry for the different regional boreholes. The challenge is now to model the transport of a saline porewater originating in the northeast and migrating horizontally in the Dogger Formation and then vertically and probably slowly in the Callovian–Oxfordian formation.
The porewater chemistry model will also enable us to interpret future data obtained experimentally, directly in the ANDRA underground research laboratory.
Acknowledgement
This work was funded within the framework of the ANDRA–BRGM partnership (THERMOAR project 1998–2005). We thank H. Rebours (ANDRA/EST) and the ERM team (Poitiers) for their contributions. The reviews of Dr Pierre de Cannière (CEN/SCK, Belgique) and of Dr Gilles Berger (CNRS, France) have been very helpful.