Open access

Claystone as a Potential Host Rock for Nuclear Waste Storage

Written By

Károly Lázár and Zoltán Máthé

Submitted: 06 November 2011 Published: 12 September 2012

DOI: 10.5772/48123

From the Edited Volume

Clay Minerals in Nature - Their Characterization, Modification and Application

Edited by Marta Valaškova and Gražyna Simha Martynkova

Chapter metrics overview

3,112 Chapter Downloads

View Full Metrics

1. Introduction

Nuclear energy is widely used for production of electricity. There are several advantages of this type of generation of electric power, but severe drawbacks also emerge. These difficulties should be eliminated for appropriate large-scale utilization. One of these difficulties is the formation of long life-time isotopes during the nuclear reactions (fission and capture of neutrons). The life-time of some of these products reaches even the 105-106 year range. There are procedures, by which these long life time isotopes can be extracted or eliminated (reprocessing and transmutation), but finally, some isotopes are still remaining with significantly decreased activity with life-times in ca. thousand years range.

The reliable and safe deposition of these products should be elaborated. The long-term geological disposal is considered as a suitable option for the isolation of these isotopes from the biosphere during their long life-time. Properties of various types of host rocks were evaluated from this aspect, crystalline (granitic), clayey or even salt types among them.

The clayey host rocks exhibit some andvantages in this respect, namely they may have significant capacity for sorption and ion exhange, both processes considerably may retard the migration of dangerous components. Some of these formations have been characterised in very details – even underground laboratories were established in them for their in situ characterisation (Callovo-Oxfordian formation in France, Opalinus Clay in Switzerland, and Boom Clay in Belgium).

The issue is about to be also addressed in Hungary. C.a. 40 % of electricity is procuded in a nuclear power plant since the 1980’s. A significant amount of collected spent fuel has been stored in temporary repositories since then. The procedure for final deposition should be elaborated and implemented within the next 20 – 30 years.

A preliminary screening had been performed, the Boda Claystone Formation has emerged as a potential media for deposition. This formation is an extensive one – samples collected from different parts of it have been characterised. An underground research facility had even been established in it in a depth of 1050 m from 1994 – 1999. A detailed characterization of the formation has been performed based on the data collected from results of measurements performed in this facility.

The intention with the recent chapter is to provide an illustration and to present some aspects of the evaluation on the potential possibilities of the application of clayey geological media for storage of wastes of nuclear origin. Following this short introduction in the second part an overview is given on the types and most important long life-time isotopes, and on the ongoing processes which might occur during the several thousand years of operation time of the waste storage facility. In the third part a short account is presented on the properties of clay minerals which influence their behaviour as an isolation media. Results of some specific evaluations of clays (Callovo-Oxfordian, Opalinus) are also presented in this part to provide an overview and comparison. In the fourth part a more detailed description and characterisation of Boda Claystone is presented. Its formation, components, characterisation of clay minerals is described in more detail. Special attention is devoted to measurements correlated to isolation properties of this media against migration of long life-time isotopes.

Advertisement

2. Nuclear wastes

The production of nuclear wastes is connected primarily to generation of electricity in nuclear power plants. Other sources of nuclear wastes are common as well. For example those radionuclides which had been used for therapeutical or industrial applications. These latter radioisotopes can be generated by neutron or proton irradiation of certain stable isotopes, they usually have short half life-times (less than ca. 30 year), and their amount is comparatively small. The nuclear wastes usally are classified to two main groups: to low-level and to high-level wastes. The criteria for the distinction are the half-life and the amounts of the radionuclides present in them. The main components in the low-level wastes have short half life-times, they can be stored in ground surface facilities, their operation time is expected to elapse for a few hundred years. In contrast, the radioisotopes in the high level nuclear wastes have considerably longer half life-time (up to the range of 105 – 106 year). Construction of underground facilities is considered for the final disposal as a reliable mean to maintain an appropriate isolation from the biosphere for the long time interval necessary for the decay of the critical isotopes [1]. One option for the disposal of high level wastes can be the utilisation of clays as isolating media. Prior to discussing some properties of clays it is worth to discuss some properties of the high level wastes in more details with particular respect to their possible interactions with clayey minerals occuring during the long operation time of the disposal facility.

2.1. High level nuclear waste

In conventional nuclear reactors the energy is produced by splitting 235U nuclei with neutron irradiation. In a conventional fuel type the 235U content is enriched to ca. 4 % from the natural 0.6 – 0.8 per cent in the ores. The rest of the fuel is 238U (both isotopes in oxide form). The fission of the heavy U nucleus results in formation of two novel nuclei, and a few neutrons are also emitted which may propagate the process. A small difference exists between the masses of the starting and final products. The corresponding energy equivalent of this difference is the generated heat in the final gross balance. A variety of nuclei are produced in the fission, a number of them is instable and will be stabilized in later processes by emitting or radiation. There are several instable nuclei among these fission products which have considerably long half life-time (Table 1.) Further on, in the nuclear reactors not only fission but capture of neutrons may also proceed. For instance 238U may capture one neutron and the with simultaneous radiation 239Np forms. This nucleus is also instable and emits one electron by forming 239Pu. This pruduct is still instable but is has a much longer life-time. Plutonium can also be used as a reactor fuel, but for splitting of 239Pu nuclei other type of reactor conditions are optimal [2].

IsotopeHalf-timeRel % energy
106 yearkeV
99Tc0.2116.1292
93Zr1.535.490
135Cs2.36.9210
129I15.70.8198

Table 1.

Some typical long life-time fission isotopes produced in nuclear reactors

A large variety of other short half life-time isotopes are also formed, they influence the working conditions of the reactor. After a certain period of the usage (3-4 years) the fuel rods are removed from the reactor. They can be stored afterwards without any further utilization. They are highly radiating, a great amount of heat is still generated in them (without exposed them to neutron irradiation). Thus, after operation they are usually kept in a water pool, then in other temporary storage facilities. But, as was shown above, isotopes with long half life-time are present in them, their final safe disposal should be provided in some manner. This type of the straight utilization of the nuclear fuel is called open fuel cycle.

Another option is the reprocessing. The used fuel still contains 235U amounting to ca. ~ 1 %, and it is also worth to extract the newly formed transuranic isotopes (eg. 239Pu) which can be utilized as a nuclear fuel (although in modified type of reactors). During the various steps of the reprocession the amounts of the long life-time isotopes are significantly decreased, but they are still present in significant amounts at the end of the procedure (Figure 1). This version of utilization is named closed fuel cycle. It is seen that the activity is rather high, even after c.a. one million years of storage the tera-Becquerel level is still maintained.

The transmutation is also a further possibility to reduce the amounts of the radiounuclides present in the high level waste. Specific controlled nuclear reactions can be initiated by irradiating different target isotopes with accelerated electrons, protons, etc. and the energy of stabilisation from the ustable target nuclei can also be extracted. This process is very sophisticated, should be performed in different manner with each isotopes separately, thus its practical application cannot be expected in the near future.

Figure 1.

Activity of high level waste remaining after processing of 1 ton spent fuel in dependence of time (reprinted with permission from [1].)

This brief overview can be summarised with the conclusions that long life-time isotopes are byproducts in generation of electricity in nuclear power plants. The life time of some of these isotopes may expand to even 106 years. The amounts of this isotopes can be decreased by various procedures but they cannot be fully eliminated.

Due to their long life-time these waste isotopes should be isolated from the biosphere. By commonly accepted recent considerations a perspective solution might be the disposal in geological media, i.e. host rocks.

2.2. Disposal and subsequent possible processes in host rocks

Two types of barriers are usually mentioned with respect to the geological disposal, the engineered and the natural ones. The engineered barrier comprises the immediate capsulation (metal cannister), the backfill material (high sorption capacity porous material – e.g. bentonite, etc.) and the strengthened wall of the shaft (made of concrete) in which the cannisters are situated. The natural barrier is the geological media. In the first few hundred years the engineered barrier should sustain the isolation from the long half-life isotopes. During this period the components of this barrier will probably lose their ability for protection due corrosion and other hydrothermal processes, Thus the role of the natural barrier will probably gain emphasised importance later on [3].

Several processes may take place in the working cycle of a disposal site. It should be taken into account that the extent of decaying nuclei is large, thus, a significant amount of heat is generated simultaneously. Further on, most of the design schemes consider deposition sites below the ground level, thus presence some amount of water cannot be excluded either. In addition, due to radiation effects, free oxidizing or reducing radicals may also form, (e.g. OH from water).

The situation can be illustrated for example with the case of neptunium on a Pourbaix (Eh-pH) diagram [4]. At the closure of the disposal site the environment is exposed to air, the pH of the ground water (in equilibrium with the CO2 of the air) is c.a. 5.5 (position A in Figure 2.). After the closure the site is isolated, the pH approaches to the mean of the rock (usually it is between 7 and 8). As the time elapses, the concrete sealing around the engineered barrier starts to erode, the pH in the close environment may increase to ~10. Later, very slowly the system may approach the state charcterising generally the whole isolating rock host (Eh -0.1 V, pH 7.5 can be estimated, see Figure 2).

Figure 2.

Schematic representation of the change of conditions on the Eh-pH diagram of Np during the life-time of a disposal site. Processes start with the construction, from a stage marked with A. (the base Eh-pH diagram is taken from [4])

It is seen in Figure 2 that Np is present as a cationic species in NpO2+ form at the start, which can be sorbed in aqueous media on the components of the rock minerals. However, the speciation changes with the increase of pH, the anionic Np(OH)5- component will be dominant later. This component probably will migrate easily due to its negative charge. Having migrated for several hundred years NpO(OH)2 component may be formed at the end, which is neutral, thus it will precipitate.

The simplified example may already provide an impression that the speciation and the pH dependence may be different for the various long half life-time isotopes. The real situation is far more sophisticated: the temperature may reach 80 – 100 °C (depending on the time elapsed, and the depth), the hydrostatic pressure may be in the 50 – 80 bar range (depending on the depth), various interlinked geochemical processes may occur under these conditions which influence the spreading of isotopes in the host rock (for further details see e.g. in [5]).

2.3. Simple practical approaches to estimate the speed of migration

The speed of migration of the various long half life-time components is a crucial factor for the evaluation of the properties of the perspective host rocks. In the first approach the migration of isotopes can be considered as a cyclic repetition of frequent sorption/desorption steps in a slowly flowing aqueous media. The strength of sorption can simply be characterised by the distribution coefficient, Kd.

Kd=I0IeIeVmE1

where I-s are concentration (activity) values of the studied isotope, I0 in the starting solution, Ie in equilibrium, after having kept the m mass of sample in contact with V volume of solution.

The Kd values of various rock samples can easily be obtained, the concentration change of the particular ion should only be determined experimentally. Consequently, Kd values are widely used. Further on, the speed of the migration relative to the slowly flowing media can also be estimated by using the Kd values:

vradvgw=11+(1εε)ρKd,E2

where vrad is the velocity of migration of the given radionuclide, vgw is the velocity of the flow of the ground water, is the porosity, is the density and Kd is the distribution coefficient of the isotope in question.

Equation (2) has been widely used since Kd values can easily be measured. However it should be used with particular precaution. Namely, it should be taken into account that:

  • The structure of the powdered rock is different from the real one – usually not all the mineral components are contacted with the aqueous medium in the real case.

  • The mass ratios of the liquid solution related to the amount of the solid rock are distinctly different in the Kd determinations and in the real situation. In the former instance the liquid component is applied in large excess (at least ten-fold), whereas the conditions in the reality are the opposite, usually only a few per cent of pore water is present in the rocks.

  • Equation (2) provides a relative, retention velocity. It may occur that the flow is practically zero in the rocky media, whilst the simultaneous diffusion takes place still. Under these conditions Eq. (2) cannot be applied at all.

There are other simple experimental methods by which more reliable determination of the crucial constants can be provided. Rock samples preserved in conditions closer to the natural ones are suitable for this type of studies, obtained for example from bore cores. One simple type of these more reliable measurements is the so-called „break-through” type. Two compartments are separated with the studied borecore sample in the measuring cell. One compartment contains the very dilute solution of the studied tracer isotope (in ground water), and the break-through (i.e. the appearance) and the rate of further increase of its concentration is followed in the ground water in the other, opposite compartment. If the trace isotopes are used in really low concentrations the experimental conditions may actually be close to the real, natural ones. In other words the applied small concentration gradient hardly modifies the natural conditions. On the other side, however, rather long time intervals are necessary to perform these studies (several hundred days). From these measurements the effective diffusion constant (Deff) can directly be determined (Eq. 3)

C'(t)C0=ADeffV'LtαAL6V',E3

where C0 is the original concentration (activity) of the isotope in a compartment at the start of the experiment, C’(t) is the concentration measured in the opposite compartment, after having the break-through started, A is the cross section of the borecore sample, L is the thickness of the borecore, V’ is the volume of the compartment. is the so called retention factor, its value is related to the delay of the break-through. Deff is proportional to the slope of the increase of concentration of the studied isotope with time in the originally non-active compartment.

Deff is one of the most appropriate parameters which can be used for the characterisation of rocks from the point of view of migration of isotopes expected to proceed in them.

Advertisement

3. Application of clayrocks for the disposal

Clay rocks are widely considered as appropriate media for final waste disposal site. Clay minerals have several advantageous properties for the isolation. They have open, layered structure in which they easily ad- or chemisorb cations.

3.1. Recent investigations on clayrocks

Due to the advantageous properties of clayrocks significant amount of efforts is devoted to investigate them as optional disposal sites. For instance, at the Nuclear Energy Agency (NEA) a working group has been formed devoted to study the scientific bases for stability and buffering capacity of deep geological waste management systems („Clay Club”). This Club organised several conferences, and the results have been published in workshop proceedings [6]. Other specific conferences have also been devoted to this topic (see e.g. [7]).

At present, three clayrock types have been characterised in very detail in Europe. Namely, the Callovo-Oxfordian formation in Bure, France, the Opalinus Clay in Mont Terri, Switzerland, and the Boom Clay in Mol, Belgium [8]. Some basic properties and some characteristic minerals of these clays are summarized and compared in Table 2 [9]. Properties of Boda Claystone are also included for later comparisons.

ClayrockAgeDepthMax.
temp.
Pore
water
Organic matterClay fractionAlbiteCalcite
106 yearm°C%%%%%
Boom30 - 3620016~ 252 - 335 - 651 - 2~ 1
Callovo-Oxfordian15550040~ 7~ 1~ 40~ 1~ 35
Opalinus18050085~ 7~ 1.545 - 60~ 14 - 8
Boda2501000220~ 2-35 - 504 - 83 - 5

Table 2.

Certain characteristic data of clayrock types considered as perspective host [9].

Principally, each rock types contain clays in large proportions. Hovewer, in other properties the clays exhibit significant differences. Their ages, the maximal temperatures experienced during formation, water contents, characteristic minerals present in them etc. are distinctly different, attesting their unique history of formation for each of them.

The safety aspects for the construction of disposal sites have also been discussed and evaluated in detail (see eg. for the Callovo-Oxfordian in [10], for the Oplalinus Clay in [11]).

Advertisement

4. Boda Claystone

There had been considerations with respect to the future disposal sites for high level nuclear wastes in Hungary, too. Boda Claystone emerged as potential host. A brief comparison of data shown in the previous Table 2 shows that each clayrock has distincly different properties – each of them should be characterized specifically. Thus, Boda Claystone was characterised in the course of various programmes as well. Some results of these characterisations are summarized below.

4.1. General characterisation

The Upper Permian sedimentary sequence of the Boda Claystone Formation (BCF) is located in Western Mecsek Mountains, southern Transdanubia, SW Hungary. The Mecsek Mts. is part of the Tisza Megaunit comprising the basement of the south-eastern half of the Pannonian Basin [12] (Figure 3). The geological map of BCF is shown in Figure 4.

Following the Variscan orthogeny the continental sedimentation in the Mecsek Mts. began in the Early Permian and lasted until the Lower Triassic. The BCF is part of this about 2000-4000 m thick siliciclastic sequence (continental red beds); its transition with underlying Cserdi Formation is conformable, while the boundary with overlying strata of the Kővágószőlős Formation is usually sharp (Figure 5) [13].

Figure 3.

Overview of plate tectonic units in the Pannonian Basin (after [12]).

Figure 4.

Geological map with depth contour of the top of the Boda Claystone Formation and locations of studied objects (boreholes: Bat-14, Ib-4, Alfa Drift) after [13]

Figure 5.

Idealised lithological column of Boda Claystone Formation [13]

On the basis of data of boreholes and geological mapping BCF is known in an area of 150 km2. Within this area approximately 15 km2 can be found on the surface in W Mecsek Mountains (see Figure 4). Two distribution areas of BCF are known in W Mecsek Mountains: i/ perianticlinal structure of the W Mecsek Mountains; ii/ so called Gorica block. In the Gorica block outcrop of BCF is not known; in this block only the borehole Ib-4 recovers sequence of BCF in significant thickness (between 494,2 and 709 m) (Figure 4).

On the basis of the deep drillings total thickness of BCF is estimated to be about 700-900 m (in perianticlinal structure of the W Mecsek Mountains) whereas according to our knowledge its thickness is smaller in the Gorica block (about 350 m). The sediments of the BCF are dominantly red and reddish brown in color, reflecting the dominantly oxidizing nature of the depositional and early diagenetic environments [14-18] (Figure 6).

Figure 6.

Reddish brown albitic analcite-bearing claystone with interbedded aleurolite (Al) and dolomite (D) layers (Ib-4, 591-596 m).

The BCF was deposited in a shallow-water lacustrine environment (playa mudflat, playa lake), under semi-arid to arid climatic conditions. According to our present-day knowledge the middle thickest unit of the BCF has only one reduce interbedding (greyish black albitic claystone containing pyrite and finaly disseminated organic matter), its tickness is about 3-4 m. However, several reduce thin layers (green, greenish-gray claystone, siltstone) can be observed in its lower and upper transitional zones (Figure 5).

On the basis of mineralogical investigation with X-ray diffraction (XRD), differential thermal analysis (DTA), and electron microscopy (EMP), the main rock-forming minerals of the BCF are: clay minerals (absolute dominant are illite-muscovite and chlorite; smectite, kaolinite, vermiculite and mixed-layer clay minerals were identified in considerable amounts), authigenic albite, quartz, carbonate minerals (calcite and dolomite) and hematite [16-19]. In addition, barite, anhydrite, authigenic K-feldspar and detrital constituents (muscovite, biotite, chlorite, zircon, rutile, apatite, ilmenite, Ca-bearing plagioclase) were always also identified in trace amounts. The authigenic albite is present as albite cement (typical of all rock types of formation), albite and carbonate-lined disseminated irregular white voids (typical of albitic claystone) and albite replacement of detrital feldspars in sandstone beds [16-17, 19].

Six main rock types of BCF can be defined based on mineralogical, geochemical and textural considerations, namely, albitic claystone, albitolite, „true” siltstone, dolomite interbeddings, sandstone, and conglomerata [13, 15-17, 19,]. Their mineral compositions are shown in Table 3.

Rock typesclay minerals (wt %)authigene albite
(wt %)
quartz
(wt %)
carbonates (wt %)hematite (wt %)
albitic claystone20-5020-505-10107-10
albitolite,<25"/>50<10105-6
„true” siltstoneapprox. 10"/>35"/>25approx. 105
dolomite interbeddings1030-40535-505
sandstone5-traces25-40 *20-305-205-traces

Table 3.

Mineral compositions of the main rock types in Boda Claystone Formation

This mineralogical composition is typical of in perianticlinal structure of the W Mecsek Mountains. The principally dominant rock type of the formation is albitic claystone.

BCF recovered from borehole Ib-4 (Gorica block) differs in its mineralogical composition. This succession of BCF contains abundant analcime in significant amounts in addition to the above minerals, in a range between 8 and 25 wt %; (typical examples are shown in the table below).

Samples
Ib-4
clay mineralsauthigene albitanalcimequartzcarbonateshematite
wt %
527,2 m4281213186
538,7 m3471628105
560,64 m4316203108

Table 4.

Mineral compositions of typical samples taken from Ib-4 borehole in Gorica block

In the Gorica block the same rock types can be defined based on mineralogical, geochemical and textural considerations as in the perianticlinal structure of the W Mecsek Mountains.

The Boda Claystone Formation is divided into three main sections [13]; for the characterization of the units see Figure 5. In the Middle Cretaceous on the basis of thickness of overlying strata in perianticlinal structure of West Mecsek Mts. the BCF was located at least at 3,5 to 4,0 km burial depth. Data determined in area of the perianticlinal structure of the W Mecsek Mountains point to late or deep diagenesis, max. 200-250 °C according to illite and chlorite crystallinity as well as vitrinite reflectance [17]. Higher illite and chlorite crystallinity is determined in core samples of the deep drilling Ib-4 (Gorica block), suggesting that BCF in Gorica block underwent lower grade diagenesis.

There were various stages of exploration and detailed characterisations of BCF. For a period even an Underground Research Laboratory was established and had been maintained in a depth of 1050 m below ground level (1994 – 1998). The results of studies performed there are collected in various reports. The data had been analysed from the aspects of criteria for establishing a waste disposal site, the main conclusions of them are compiled in a Digest [20]. The most important conclusions were presented also at one of the „Clay Club”s conferences [21].

4.2. Specific characterisations - samples

In a further stage of investigations selected Boda Claystone samples were used for measurements performed by applying directly radionuclides to obtain informations on the isolation properties of the rock against long half life-time isotopes.

Sorption and diffusion studies were performed by using fission isotopes and one of the most relavant acitnoide element, uranium. Among others samples from Delta-9, Bat-14 and IB-4 borecores were selected.

Delta-9 deep drilling is located in Alfa Drift at the 1050 m depth below ground level, it is nearly a horizontal borehole; Bat-14 was drilled near the outcrop of BCF, here the formation is covered by 25,5 m thick Quaternary sediments, and Ib-4 borehole starts also from the surface in the Gorica block, and the studied samples were collected from a distant depth 570 m (for their locations see Figure 4). Some specific studies on these samples are presented in the following paragraphs, the results of measurements with radionuclides are presented afterwards.

4.2.1. High resolution SEM with electron diffraction

These measurements were performed on samples from borehole IB-4 (570 m). A SEM/Electrondiffraction image is shown in Figure 7 from a small area of the sample. The corresponding general mineral composition is shown in Table 5.

Figure 7 illustrates the structure of the sample. Very small crystals originated from the secondary processes form the consolidated clay rock. It is also shown that within small distances the composition may change significantly (see e.g the lower right corner – it is more abundant in analcime than the rest). Parallel porosity measurements were also performed from the same regions. It was found that the porosity at the analcime rich bottom right corner is ca. 1.5 %, whereas in the upper part is 3.5 %. In contrast the portion of calcite is larger in the top regions

Figure 7.

Mineral composition obtained from a SEM & electron diffraction image of a 5 x 2.5 mm section from IB-4 (570 m) sample (IS mel: mixed illite-smectite. The image was recorded in ERM, Poitiers, France)

MineralContent (%)
Illite -
muscovite
33
Chlorite2
Smectite2
Analcime23
Albite12
Quartz11
Calcite9
Dolomite1
Hematite7

Table 5.

Average mineral composition of the IB-4 (570 m) sample determined by quantitative XRD analysis

4.2.2. Mössbauer spectroscopy

Oxidation and coordination states of iron can conveniently be studied by Mössbauer spectroscopy. The method can also be used to analyse iron-bearing minerals. (Corresponding Mössbauer parameters of over 400 iron containing minerals are collected in the handbook [22].) Boda Claystone contains also iron-bearing minerals in considerable amounts (hematite, clay minerals), thus samples can also be analysed, and some processes occured in the formation of the rock can be traced.

The ferric iron oxide, hematite exhibits a magnetically split characteristic sextet. Iron in the clay minerals displays two-line doublets in the spectra. Doublets of Fe2+ and Fe3+ states can clearly be distinguished, their positions are different. Two examples are shown below, in Figures 8 and 9. The effect of weathering can be tracked in Figure 8 in which spectra of samples collected from layers close to the surface (below the 25m Quaternary sediments) are shown for Bat-14 borehole. The very right peak in the spectra belongs to the Fe2+ ions in chlorite mineral [22]. It is seen that the intensity of this peak inreases in correspondence with the sampling depth. I.e. clay minerals in the top layers are much abundant in Fe3+ whereas in the deeper layers Fe2+ starts to dominate in clay minerals. On the other side, the amount of hematite (marked with H in the right side spectra) practically is the same in each sample. Thus it can be concluded that Fe2+ Fe3+ oxidation had taken place in the layers close to the surface, in other word the process of weathering is reflected in the Fe2+/Fe3+ change.

The conditions controlling the formation of minerals can be deduced from specra of Figure 9. Boda Claystone is described in the previous parts of this chapter as a rock which was formed essentially under oxidizing conditions. However a few exceptional strata also exist.

Figure 8.

Mössbauer spectra on samples collected from borehole BAT-14. The main iron bearing mineral is hematite (left spectra). Effect of weathering is shown in the right side spectra, approaching to the surface Fe3+/Fe ratio increases in the clay mineral (chlorite)

Figure 9.

Crossing a reducing zone at 83.7 m in Delta-9 borehole. Iron is present only in Fe2+ form in chlorite. Hematite starts to appear in other distances.

For example Delta-9 borehole crosses a few meter thick layer in which reducing conditions controlled the formation of the minerals. This layer is located ca. at the 83 – 84 m distance from the commencement of the borehole (Note: the borehole was drilled at the 1050 m depth.) The overwhelming part of iron is present in Fe2+ form in chlorite – hematite is not present at all in the top spectrum of Figure 9. By moving apart a few meters from this layer hematite appears and its presence becames prevailing by moving even farther.

Both examples show that iron ions located in the layered clay minerals are sensitive to the change of redox conditons and they can be used as indicators to monitor the alterations in them. Further description and interpretation of these measurements can be found in [23].

4.2.3. Reduction/oxidation studies from voltammetry of microparticles

The reduction/oxidation properties of rocks and minerals can also be studied by voltammetric measurements. Even small amounts of them can be analyzed by voltammetry of microparticles. Samples shown in Figure 9 were also characterised at the LCPME – CNRS, Nancy, France). For illustraton, Figure 10 displays some cyclic voltammograms obtained on samples collected from Delta 9 borehole.

The results of voltammeric measurements are in good correspondence with the Mössbauer analysis as Figure 10 displays. High anodic current (A1 peak at 0.65 V – relative to standard calomel electrode) could be detected which can be attributed to extended oxidation of Fe2+ ions on the sample formed under reducing conditions (collected from 83.7 m position). On the other two samples formed under oxidative environment reduction of iron ions ions can be observed (C1 peak probably corresponds to reduction of ions in the dissolved ions, C2 at –0.45 V to the reduction of ferric ions in the clays and C3 at -0.74 V is for ferric ions in hematite). More extended description of the method, samples and measurements can be found in [24].

Figure 10.

Cyclic voltammograms obtained on samples collected from Delta-9 borehole. The numbers in the top left corner represent the relative amounts of iron bearing minerals extracted from Mössbauer spectra of Figure 9. A1 stands for anodic, C1, C2 and C3 for cathodic peaks, Ei=0 is the zero potential. (The voltammograms were measured at LCPME Nancy, France)

4.3. Studies with long life-time radionuclides

4.3.1. Sorption measurements

In the first stage of investigations sorption of radionuclides, mostly fission products, were investigated on various types of samples. Namely Kd values were determined for 125I, 137Cs, 60Co and 85Sr on various Boda Claystone samples collected from boreholes drilled in the 1050 m depth exploratory tunnel. Measurements were performed with 2 g amounts of crushed samples. Radionuclides were added in trace amounts to 20 ml volume of ground water (10-5 – 10-3 mol/L). Kd values were determined by using Eq. (1). There was a scatter in the results depending of the samples, for 125I 0.2 -2.0, for 85Sr 60 – 120, for 137Cs 600 – 5000 and for 60Co 1200 – 15000 values were determined. Concentration dependence of the sorption was also measured by diluting the tracer isotopes with inactive (natural) isotopes. Sorption of caesium showed a regular behaviour, the constants of a non-linear Freundlich-type isotherm could also be determined. Similar results could not be obtained with Sr2+ and Co2+, both ions have low solubility in the applied ground water, they easily form precipitates (SrCO3, and Co(OH)2). Further details and results are described in [25].

It should be mentioned that similar sorption measurements were also performed on samples originated from another facies of Boda Siltstone, rich in analcime component originated from IB-4 borehole (see Figure 7 and Table 5). In these measurements non-radioactive natural isotopes were applied with higher concentrations, their amounts were determined by atomic absorption spectroscopy measurements. The results are in good correspondence with those obtained with radioisotopes [26]. The obtained Kd values might be applied for the further evaluation of the isolation properties of the rock against migration of radionuclides by using Eq. (2).

In relation with Eq. (2) it was already pointed out that principally it decribes the relative retardation of the velocity of migration of radionuclides compared to the velocity of hydraulic flow. The equation is primarily valid for cations, tending to adsorb on negatively charged clay minerals in the aqueous media. Negatively charged species hardly sorb as is reflected for example with the mentioned low Kd value of 125I. Further on, Boda Claystone has small porosity (1-3 %) and the water permeability also lies in an extremely low value region (10–20 – 10-23 m2) [27]. Thus, the velocity if hydraulic flow can also be practically negligible. Consequently, there are at least two reasons that the velocity of spreading of negatively charged anionic species cannot be estimated by using Eq. (2). One reason is that the vgw term it close to zero in the equation, the other reason is that anionic species do not sorb on clay minerals.

In order to obtain reliable information on the migration of negatively charged anionic species another option should be applied, namely to perform direct diffusion measurements.

4.3.2. Measurement of the diffusion coefficient of anionic species and HTO

To obtain approximately reliable information on the diffusion of radionuclides it is advised to use compact rock samples with their preserved consolidated structure. Various arrangements can be used. One option is to use a simple measuring arrangement, the break-through cells as mentioned in Section 2.3. The effective diffusion coefficients can be experimentally determined by using Eq (3).

It is also possible to collect information on the permeability and self diffusion of water itself by using a particular radiotracer, that is tritiated water, HTO.

Break-through measurements were carried out with HTO and with anionic long life-time fission product tracers (99TcO4- and 125I) and additionally, with H14CO3-. The compartments of c.a. 90 volume in the break-through cells were filled with ground water. The studied 8 mm thick borecore disc samples were obtained from the 570 m depth in the IB-4 borehole. Ground water in one side of the compartments was separately spiked with tracer amounts of the respective isotopes (10 MBq HTO, 6.2 MBq 99TcO4- and 8 MBq H14CO3-, respectively). These radionuclides emit -radiation, the activities in small amounts of samples (10 microliter) taken regularly from the compartments were determined liquid scintillation detection. As an illustration, the measurements performed in a series with H14CO3- are shown in Figure 11 (next page).

The top part of the figure shows that after a minor decrease the activity (i.e. the concentration of H14CO3-) is constant. The bottom part displays the break-through, the appearance and the increase of the concentration of the radionuclide in the ground water in the compartment on the other side of the sample disc. Note that only c.a. 2 % of the total amount of HCO3- penetrates through the 8 mm thick disc within a year. The delay in the commencement of the break through indicates the interaction between the solute and solid phases. The intercept on the horizontal axis is c.a. 220 days, this is the value of (retention coefficient) in Eq. (3).

Figure 11.

The relative activities of 14CO32- in the active (top) and originally inactive (bottom) compartments of the break-through cell at pH = 12 on two samples (marked with squares and triangles).

The delay in the break-through can probably be attributed to the exchange of radiocarbon to the carbonates in minerals of solid phase. Namely, the mesurements presented in Figure 11 were performed at pH=12. At this pH the 2 HCO3- CO32- + H2O + CO2 equilibrium is shifted to the right, ie. CO32- ions prevail over the hydrocarbonates, and some of the 14CO32- ions may be exchanged from the solution to the solid phase with carbonates in calcite (dolomite). The main data extracted from the measurements and evaluated in correspondence with Eq (3) are collected in Table 6. Further experimental details, description of samples and measurements performed with the -emitting isotopes can be found in [28].

Migrating
species
Deff
10-12 m2 s-1
HTO14~ 0
H14CO3-~ 1
~ 1
~ 0 a
1.5 b
99TcO4-~ 1~ 0
125I-~ 8 - 14 c0.3 – 1.3 c

Table 6.

Effective diffusion coefficients and retention factors ( in Eq. 3) of neutral (HTO) and anionic migrating radionuclides

As for investigations related to the diffusion of iodide anion it should be mentioned that several more extended measurements have also been performed by using -emitting 125I isotope (35.5 keV with half life-time only 60.2 days). The -radiation allows a more easy detection of the radiating isotope, for example presence and acitivity of 125I in solid samples can directly be measured. The migration behaviour is expected to be the very same as the long half life time fission product 129I (1.7 107 years).

Two additional types of diffusion measurements were performed with 125I. The in-diffusion was measured from spiked groundwater into borecore samples, and concentrations of iodine were determined in slices cut from the end of the solid borecore. In another arrangement the „in situ” conditions were modeled by placing the borecore samples under 100 bar pressure at 50 °C during the period of the measurements. These conditions correspond to the values of the ambient pressure and temperature in a rock located in 1000 m depth (characterising e.g. the former exploratory alpha drift). Two sets samples were used for these measurements, a few one was kept under these conditions for 118, another few for 340 days, the iodine contents were measured in the solid slices of sample borecore afterwards.

The results of the three types of measurements of diffusion coefficient of iodide (the break-through, the in-diffusion and the „in situ”) were in good correspondence with similar values. Further experimental details and evaluation are described in [29].

As for the general evaluation of measurements with anions, data of Table 6 are in good correlation with the expectations. Namely, water and the anionic iodide exhibit high mobility, with similar Deff values. This reflects the lack of interactions with clays. Hydrocarbonate and technetate anions move more slowly. In the first instance some interaction (exchange) may take place between carbonate ions in the solution and in the solid phase. Technetate may move slower due to the larger radius of the hydrated anion. Further on, it is worth comparing these values with the self diffusion coefficients of ions in strong aqueous electrolytes. For this end, the Bruggeman relation should be recalled which establishes connection between the measured effective diffusion coefficient (Deff) and the self diffusion constant (Ds) through the porosity () [30]:

Deff=ε3Ds.E4

The self diffusion constant for ions in aqueous electrolytes is 1.7 10-9 m2 s-1 [31]. Taking into account the 1.5 – 3 % porosity, the measured Deff values correlate well to the self-diffusion constant values characterising the electrolytes in aqueous soulutions. The similarity of the compared values reflects that there is no significant interaction between the anionic component and the minerals of the clay – they diffuse through the rock without any interaction with the pore walls (provided the rock capacity factors are close to zero).

As an upper limit for the velocity of the migration of long life-time isotopes the obtained 10-12 m2 s-1 can be considered. (The cations migrate more slowly since they interact and may became sorbed on the minerals.) As a very rough approximation the

xDefftE5

distance (radius) – time relation can be considered within which the radionuclides spread during their life-time. Inserting the long half life-times values into this approaching equation the order of magnitude of migration distance can be estimated. For illustration, for 99TcO4- (half life-time 2.1 105 years) ca. 1.4 m is the distance, for the longer lifetime 129I (1.7 107 years) ~ 30 m is the distance of spreading within one half lifetime period in Boda Claystone media. These distance values should be considered with precaution. First, only one half lifetime period was considered, the amount of radiation emitting isotopes will be only halved within this period. The second point is our starting precondition, namely the hydraulic flow was neglected. Thus, in real cases the distances of migration can be significanly longer.

Some properties of the other clayrock types considered for final disposal for nuclear wastes were also briefly mentioned in previous section 3.1 (Boom Clay, Opalinus Clay, and Callovo-Oxfordian argillite). Here, at the end of the present section dealing with diffusion constant measurements on Boda Claystone samples, a reference should also be given to an overview, which provides a detailed compilation and description of recent diffusion studies performed in the framework of the EC FUNMIG IP project on all the mentioned four types of clayrocks [32]. In general diffusion data obtained on Boda Claystone suit well within those collected on the other three types of rocks which had been extensively studied and characterized in the project.

4.3.3. Break-through studies with uranium

In the preceeding sections mostly sorption and diffusion properties of fission products were discussed on Boda Claystone samples. To consider the behaviour and interactions of the long life-time actinide components with minerals in Boda Claystone the major component of the high level nuclear wastes, uranium was also investigated in break-through experiments. Uranium is similar to neptunium in respect with the occurence of various oxidation states (as was shown for Np in Figure 2.) Depending on the Eh – pH conditions uranium may be stabilized either in cationic, or in anionic form or as stable neutral precipitate. The U-C-O-H diagram is shown – it can be applied for carbonate containing rocks and in correspondence hydrocarbonate contaning ground water, which is the case for Boda Claystone. Stability regions of iron are also shown and separated in black dashed lines (Figure 12 - based on [4]).

The life cycle of the disposal site is also shown in blue (similarly to Np in Fig 2). At the start uranium (or a part of it) may be dissolved in form uranyl ions (UO22+) into the ground water. This species is a cation and prefers to ad- or chemisorb on the clay minerals. After having the site closed the potential decreases, the pH increases – neutral UO2CO3 may precipitate. In a subsequent stage with the further increase of pH negatively charged carbonates can be formed which tend to migrate in the aqueous phase. In the final stage with the decrease of potential the reduction prevails and stable UO2 forms from the various migrating U(VI) uranyl carbonate species (shown as the orange encircled area in Figure 12). In the encircled region Fe2+ Fe3+ redox processes may also proceed, maybe they can be coupled with the mentioned U(VI) U(IV) reduction.

In the corresponding break-through experiments the ground water was saturated with uranyl acetate (using uranium with natural isotope composition, dominant 238U with ca. 1% 235U). The changes of activities of uranium were measured in the compartments of the break-through cells. Furthermore, the distribution of uranium along the path of migration in the bore core samples was also determined by high sensitivity laser ablation ion-coupled mass spectrometry (LA-ICP-MS) after finishing the experiments and dismounting the break-through cells. The change of concentration (activity) of uranyl ions in dependence of time is shown in Figure 13. Note the long time interval, the samples were kept in the cells for ca. five years.

Figure 12.

The Eh-pH Pourbaix-diagram of uranium. The dashed line represents the possible changes during the life cycle of the repository, the encircled region might be the final stage.

Primarily the depletion of uranium is seen in the starting solution in Figure 13. In the very first period a sudden drop takes place, c.a. 20 % of uranium disappears from the solution, most probably due to precipitation of carbonates via the HCO3- CO32- equilibrium in the solution, where almost all the carbonates are consumed. After having this fast process finished the amount of uranium had been stabilised for c.a. one year at the 70 – 80 % level of the original concentration. Later on the starting solution is considerably depleted in uranium. Uranium is apparently trapped in the solid borecore disc, simultaneous measurements on liquid samples does not show any sign of through diffusion. Negligible increase of the amount of uranium is detected in samples taken from the opposite compartment of break-through cell (bottom curve in the Figure).

In correspondence, the LA-ICP-MS measurements directly prove the event of trapping. A typical distribution of concentration (intensity) of uranium is shown in Figure 14, which displays the distribution of uranium along the migration path, crossing the borecore dics. (Notice the logarithmic scale on the vertical intensity/concentration axis.)

Figure 13.

Changes of concentations of uranium in the two compartments of a break-through cell during ca. 5 years. Top curve is the original solution, the bottom curve is for the opposite compartment.

Figure 14.

and 235U contents measured in the sample after finishing the 5 year break-through experiment. The left side of sample contacted the uranyl solution. Note the logarithmic scale on the vertical axis.

The figure illustrates the sensitivity of the method, even the minor uranium component, 235U (a few ppb in the sample, ~ 1 % in the natural uranium component), can convincingly be measured. It is seen that the characteristic longest migration distance was only c.a. 4 mm during the five years duration. (The horizontal sections between 4 and 10 mm distances refer to the 238U and 235U contents present originally in the sample before the start of measurements.) At the first approximation the drops of concentrations/intensities of 238U and 235U can be represented by straight lines in logarithmic scale, thus the decrease can be described some power dependence of distance covered by the migration. Further experimental details of the measurements are described in [33].

A rough estimation can be performed by using the square-root approximation mentioned with Eg (5). Inserting 4 mm and 1900 days into Eq (5), a nominal Deff 10-13 m2 s-1 can be obtained. This value is less with two orders of magnitude than those obtained for the mobile iodine and HTO, and is less with one order of magnitude obtained for the less mobile hydrocarbonate and technetate (shown in Table 6). Thus this comparison with neutral and anionic species can also be considered as an indirect proof the measurements with uranium. The slow migration of uranium can probably be attributed to more than one single factor. Beside chemisorption of uranyl ions on clay minerals interaction with (hydro)carbonate anions may take place with forming partly immoblie carbonate species (see the Eh-pH diagram in Figure 12). Fe2+/Fe3+ ratios were measured by Mössbauer spectroscopy in the surface layers of the borecore samples after having the uranium break-through measurements finished. They were the same as at the start. In consequence the Fe2+ Fe3+ process coupled with U(VI) UO2 reduction, as mentioned in the introductory part of this section, does not play role [33].

Advertisement

5. Summary

Certain claystones have been considered as appropriate media for long time isolation of high level nuclear waste depositories. Various interactions between the migrating anionic and/or cationic species of radionuclides with the constituents of minerals may take place. For example, the open, layered structure of clay minerals is advantageous for this purpose resulting in high ad- and chemisorption capacities to retard migration of long life-time radioisotopes. Characteristic related properties of different clayrock formations considered as possible hosts in Europe are briefly mentioned and compared. Boda Claystone is discussed in more detail, diagenesis, lithology and mineralogical properties are described. Results of various specific measurements on claystone samples are presented. Particular attention is devoted to sorption and diffusion measurements performed with anionic fission product and actinide radionuclides. From these results some characteristic distances which can be traversed by these radionuclides in the clay rock during the life time of the depository site are estimated.

Acknowledgement

The research leading to these results has received funding from the European Union’s European Atomic Energy Community’s (EURATOM) sixth and seventh Framework Programmes under grant agreements FP6-516514 (FUNMIG) and FP7-212287 (ReCosy), respectively. The authors are thankful for the fruitful informations and knowledge provided by the opportunity with the participation in these projects. Special thanks are due to the colleagues who participated the work in its different stages (J. Megyeri, P. Mell., T. Szarvas, and late L. Riess, at the Institute of Isotopes), and particular thanks are also due to colleagues who contributed to the characterisation of samples by analysing them with specific and unique methods (J.-C. Parneix and M. Perdicakis). The courtesy of PURAM for providing the samples for analysis is also appreciated.

References

  1. 1. Ojovan M.I, Lee W.E,2005An Inroduction to Nuclear Waste Immobilisation, Elsevier, 310 p.
  2. 2. Lieser K H1997Nuclear and Radiochemistry: Fundamentals and Applications,VCH, Wiley, Weinheim, 460 p.
  3. 3. Extrapolation of Short Term Observations to Time Periods Relevant to the Isolation of LongLived Radioactive Waste, IAEA-TECDOC-1177, IAEA, Vienna,2000p.
  4. 4. Brookins D G,1988Eh-pH Diagrams for GeochemistrySpringer, Berlin 176 p.
  5. 5. Brookins D G1984Geochemical Aspects of Radioactive Waste DisposalSpringer, Berlin. 347 p.
  6. 6. Stability and Buffering Capacity of the Geosphere for the Long-term Isolation of Radioactive Waste (Application to Argillaceous Media), „Clay Club” Workshop Proceedings, Braunschweig, Germany, Dec. 2003, OECD 2004, NEA 5303, 244 p
  7. 7. Clay in Natural and Engineered Barriers for Radioactive Waste Confinement, Ed: Aranyossy J-F.Special issue of Physics and Chemistry of the Earth, 322007p.
  8. 8. The Use of Scientific and Technical Results from Underground Research Laboratory Investigations for the Geological Disposal of Radioactive Waste, IAEA-TECDOC-1243, IAEA, Vienna2001p.
  9. 9. AltmannS.2009RTD Component 3, in: Fundamental Processes of Radionuclide Migration, FUNMIG Workshop Proceedings, (Ed. Buckau G et al.) Forschungszentrum Karlsruhe, Wissenschaftliche Berichte FZKA 74614582
  10. 10. DossierArgile, Safety Evaluation of a Geological Repository (2005ANDRA 782 p
  11. 11. Project Opalinus Clay, Safety Report, NAGRA Technical Report 02-05,2002p.
  12. 12. HaasJ.PéróCs.2004Mesosoic evolution of the Tisza Mega-unit. International Journal of Earth Sciences 93297313
  13. 13. KonrádGy.SebeK.HalászA.BabinszkiE.2010Sedimentology of a Permian playa lake: the Boda Claystone Formation, HungaryGeologos162741
  14. 14. JámborÁ.1964Lower Permian formations of the Mecsek Mountains. Manuscript, Mecsek Ore Environment Company, Pécs (In Hungarian.)
  15. 15. BarabásA.Barabás-StuhlÁ.1998Stratigraphy of the Permian formations in the Mecsek Mountains and its surroundings. In: Stratigraphy of Geological Formations of Hungary. Mol Plc. and Hung. Geol. Institute, Budapest, 187215In Hungarian.)
  16. 16. Summary Report of the Site Characterisation Program of the Boda Siltstone Formation19984Ed: Máthé Z., Manuscript, Mecsek Ore Environment Co., Pécs. 76 p.
  17. 17. ÁrkaiP.BaloghK.DeményA.FórizsI.NagyG.MáthéZ. (2000), Composition, diagenetic and post-diagenetic alterations of a possible radioactive waste repository site: the Boda Albitic Claystone Formation, southern Hungary. Acta Geologica Hungarica, 43 351378 .
  18. 18. -VargaR.SzakmányA.GyRaucsik. B.MáthéZ.2005Chemical composition, provenance and early diagenetic processes of playa lake deposits from the Boda Siltstone Formation (Upper Permian), SW Hungary. Acta Geol. Hungarica 484968
  19. 19. -VargaR.RaucsikA.SzakmányB.GyMáthé. Z.2006Mineralogical, petrological and geochemical characteristics of the siliciclastic rock types of Boda Siltstone Formation. Bulletin of the Hungarian Geological Society, 136201232
  20. 20. Digest on the Results of the Short-term Characterisation of the Boda Claystone Formation, Ed.Kovács L., Mecsekérc-Puram, Pécs-Paks, (1999p.
  21. 21. SzűcsI.CsicsákJ.ÓváriÁ.KovácsL.NagyZ. (2004) Confinement performance of Boda Claystone Formation, Hungary, in: Stability and Buffering Capacity of the Geosphere for the Long-term Isolation of Radioactive Waste (Application to Argillaceous Media), „Clay Club” Workshop Proceedings, OECD 2004, NEA 5303. 209224 .
  22. 22. Mösssbauer Mineral Handbook1998Eds: Stevens J.G. et al., Mössbauer Effect Data Center, 527 p.
  23. 23. LázárK.MáthéZ.FöldváriM.2010Various redox conditions in Boda Claystone as reflected in the change of Fe2+/Fe3+ ratio in clay minerals, Journal of Physics: Conference series 217: 012053.
  24. 24. PerdicakisM.XuY. L.LázárK.MáthéZ.RouillardL.2012Voltammetric characterization of Boda Albitic Claystone: Comparison with Mössbauer spectroscopy data, Electroanalysis, submitted.
  25. 25. MellP.MegyeriJ.RiessL.MáthéZ.CsicsákJ.LázárK.2006Sorption of Co, Cs, Sr and I onto argillaceous rock as studied by radiotracers, Journal of Radioanalytical and Nuclear Chemistry, 268405410
  26. 26. SiposP.NémethT.MáthéZ.2010Preliminary results on the Co, Sr and Cs sorption properties of the analcime-containing rock type of the Boda Siltstone Formation, Central European Geology, 536778
  27. 27. FedorF.HámosG.JobbikA.MáthéZ.SomodiG.SzűcsI.2008Laboratory pressure pulse decay permeability measurement of Boda Claystone, Mecsek Mts., SW Hungary, Physics and Chemistry of the Earth 33: S45S53.
  28. 28. LázárK.MegyeriJ.Parneix-CJ.MáthéZ.SzarvasT.2009Diffusion of anionic species (99TcO4-, H14CO3-) and HTO in Boda Claystone borecore samples, in: Fundamental Processes of Radionuclide Migration, FUNMIG Workshop Proceedings, (Ed. Buckau G et al.) Forschungszentrum Karlsruhe, Wissenschaftliche Berichte FZKA 7461199204
  29. 29. MellP.MegyeriJ.RiessL.MáthéZ.HámosG.LázárK.2006Diffusion of Sr, Cs, Co and I in argillaceous rock as studied by radiotracers, Journal of Radioanalytical and Nuclear Chemistry, 268411417
  30. 30. Rose D A, (1963) Water movement in porous materials, British Journal of Applied Physics, 14 256262
  31. 31. ErdeyGrúz. T.SchayG.1962Theoretical Physical Chemistry, Tankönyvkiadó, Budapest, 3p. (in Hungarian)
  32. 32. AltmannS.TournassatC.GoutelardF.Parneix-CJ.GimmiT.MaesN.2012Diffusion-driven transport in clayrock formations, 27463478
  33. 33. LázárK.MegyeriJ.MácsikZs.SzélesÉ.MáthéZ.2011Migration of uranyl ions in Boda Claystone samples, in: Redox Phenomena Controlling Systems, RECOSY 3rd Annual Workshop Proceedings, (Ed. Altmaier M, et al.) KIT Scientific Reports 76039197

Written By

Károly Lázár and Zoltán Máthé

Submitted: 06 November 2011 Published: 12 September 2012