P1Freely available online through the publisher-supported open access option.


Corrosion, water radiolysis and microbial degradation will result in the generation of gas within repositories designed for the geological disposal of high-level radioactive waste. It is therefore crucial in the design of such facilities that the relevant mechanisms allowing gas migration through repository materials, both engineered barriers and clay-based candidate host rocks, are correctly identified. In Belgium, the Boom Clay represents a candidate host material for which the advective gas breakthrough characteristics and transport properties have been extensively tested and are well defined by numerous studies. The Boom Clay displays a significant capacity for self-sealing and both laboratory and field tests indicate that advective gas transport occurs not by visco-capillary flow, but instead through the formation of pressure-induced dilatant pathways. In this study, we present results from a gas injection test designed to demonstrate the presence of these features by injecting nanoparticulate tracers with helium gas into a sample of Boom Clay. The results provide conclusive evidence for the formation of transient, dilatant gas pathways within a candidate clay-based host rock. This technique provides a novel diagnostic tool for the identification of processes governing multi-phase flow, supporting robust long-term assessments of repository performance.


A full understanding of the migration behaviour of gas generated during canister corrosion is of fundamental importance to the development of geological disposal facilities for radioactive waste. There is significant evidence that in some scenarios the rate of gas production may exceed the diffusional capacity of the engineered barriers and clay-based candidate host rocks (Weetjens and Sillen, 2006; Ortiz et al., 2002; Wikramaratna et al., 1993). If this occurs, capillary restrictions (Aziz and Settari, 1979) on the movement of gas will result in the build-up of pressure to a critical value. To accurately predict the movement of gas through argillaceous materials (Neretnieks, 1984; Kreis, 1991; Askarieh et al., 2000; Ekeroth et al., 2006; Smart et al., 2006), it is first necessary to define the correct conceptual model that best represents the empirical data. Four such models for gas flow exist: (1) gas movement by solution and/or diffusion governed by Henry's and Fick's laws, respectively; (2) gas migration by visco-capillary flow, governed by a generalized form of Darcy's Law; (3) gas flow through dilatant pathways, the permeability of which is dependent on an interplay between local gas pressure and the effective stress state; and (4) gas flow along macro fractures similar in form to those observed in hydrofracture activities during reservoir stimulation, where fracture initiation occurs when the gas pressure exceeds the sum of the minor principle stress and tensile strength. For engineering problems such as radioactive waste disposal, underground gas storage and carbon dioxide sequestration, most interest is focussed on understanding the processes and mechanisms governing the advective transport of gas (mechanisms 2–4).

Past laboratory and field tests indicate that, in the case of the Boom Clay, advective gas flow occurs through the creation of transient, localized pathways, rather than standard visco-capillary flow (Horseman and Harrington, 1994; Ortiz et al., 1996, 2002; Harrington and Horseman, 1997, 1999). These pathways have been shown to form in a number of argillaceous materials and are pressure-induced features (Horseman and Harrington, 1994, 1999; Gallé and Tanai, 1998; Horseman et al., 1999; Ortiz et al., 2002; Autio et al., 2006; Angeli, 2009), the aperture of which is a function of the internal gas pressure and structural constraints within the clay. However, the presence and spatial distribution of such features has not yet been directly imaged.

The purpose of this study is not to define the (already well characterized) breakthrough properties of Boom Clay, but to provide the first direct microscale evidence for this mechanism of flow. To help differentiate between competing theories and processes, a novel technique has been developed utilizing the injection of nanoparticles in a gaseous stream, to provide new insight into the formation, geometry and distribution of conductive pathways. In the results section we present images showing localized aggregates of gold and hydrophobic titanium oxide (TiO2) nanoparticles, successfully injected into the clay, and draping of the host fabric around nanoparticle aggregates; clearly demonstrating the principle of plastic self-sealing which is commonly attributed to such materials. This method and the resulting observations represent a completely new approach to imaging ephemeral fluid flow pathways, which will aid investigations into the scale and nature of gas migration through low-permeability clays and mudrocks. Such observations are of direct relevance to a broad range of geoscientific problems, from the disposal of heat-emitting radioactive waste to the geological storage of carbon dioxide.


A well-preserved sample of Boom Clay, deposited in the north European sedimentary province during the middle Oligocene, (Tertiary period, ~30 Ma) was taken from drilling CG R13U in the HADES underground research laboratory (located at Mol, in northeast Belgium) at a depth of 178.6 m TAW1. This was sub-sampled (to minimize perturbation) and a number of cylindrical specimens (50 mm diameter × 25 mm length) were manufactured using standard geotechnical techniques for poorly consolidated clays (dry-cutting with a stainless steel ring-former with a sharp, leading edge). Gallé and Tanai (1998) demonstrated that gas entry pressure is dependent upon the initial saturation of the clay. Prior to testing the sample saturation was estimated (assuming a grain density of 2680 kg m−3) to be ~93.1%. However, given the plasticity of Boom Clay and its capacity for self-sealing (Davies, 2005), the application of confining pressure is likely to affect the effective saturation under test conditions. Assuming ideal gas behaviour, the degree of saturation of a compressible clay-like material under confinement, Srtc, can be estimated as a function of the void ratio and saturation at atmospheric conditions, eat and Swat, respectively, as follows (Horseman and Harrington, 1997):  
where Rp is the ratio of atmospheric pressure to the air pressure in the sample when subjected to confinement, w is the gravimetric water content and Gs is the average specific gravity of the mineral solids. Under the conditions applied in this experiment, and assuming that as much as 80% of the confining pressure is carried through the sample fabric rather than the pores, this leads to an estimated effective saturation of 99%, even before allowing for any additional dissolution of gas into solution. As such, it seems unlikely that a fully interconnected pathway of gas was present within the sample prior to gas injection. This assertion is borne out by observations presented in the results section of this paper. It should be noted that, due to the introduction of additional mass (injected particles) into the sample (and remaining at the sample interface), it was not possible to derive the post-test saturation of the sample.

A 50:50 mixture (by volume) of gold and hydrophobic coated TiO2 nanoparticles (MK Impex, Canada: at <100 nm and <50 nm, respectively) was selected as a tracer for the gas injection experiment. Gold nanoparticles were chosen for their very high backscattered electron coefficient (due to density and large atomic number) and natural absence within the Boom Clay. As such, the gold can be readily differentiated by scanning electron microscopy (SEM) techniques, even in nanoparticle form. In the initial development of this technique, gold nanoparticles alone were used as a tracer, but their tendency to aggregate, led to poor dispersion in the clay. In the experiment presented here, mixing of the gold nanoparticles with additional TiO2 nanoparticles (50:50 by volume) allowed aggregation to be reduced and dispersion improved.

The test sample was sandwiched between two stainless-steel platens, with the nanoparticle mixture preloaded into a recess located within the injection platen. The sample was isolated from the confining fluid (deionized water) using a flexible Teflon sheath and an isotropic confinement of 2.2 MPa was applied to the clay. Injection and confining pressures were controlled using two high-precision ISCO syringe pumps (following the technique outlined by Harrington et al., 2009). As gas could potentially leak past the syringe seal, the injection pump was charged with water which was used to displace helium into the sample from a pre-charged vessel. This approach also ensures that the helium is water-saturated in advance, negating the possibility that it may desiccate the clay.

The mass flux of gas per unit time, t, into the clay during gas injection testing was derived using a data reduction algorithm based on the ideal gas law. At any point in time, the total volume of the injection system Vs is:  
where Vw and Vg are the volume of water and gas respectively. As the injection pump advances (i.e. water is injected into the interface vessel), the volume of gas reduces by an equal amount as compression of the gas occurs:  
By applying the ideal gas law, where P is the gas pressure, n is the number of moles of gas in the system, R is the universal gas constant and T is ambient temperature, the number of moles of gas injected into the clay is:  

Hence, by multiplying the ideal gas law by the molar volume of gas (taken as 0.02241m3 mol−1) defined at standard temperature and pressure conditions (STP2), the equivalent volumetric flow rate of gas (at STP) injected into the clay can be obtained.

At day zero, helium injection was initiated through one platen, at a constant rate of 1650 μl hr−1, and the opposing filter was maintained at atmospheric pressure. This generated an increasing differential pressure across the sample with time (Fig. 1), whilst the resulting gas inflow was continuously monitored (Fig. 2). This phase of gas testing continued for just over six days, during the course of which no clear evidence for gas breakthrough was observed. Once the applied gas pressure reached a value of 2.37 MPa, the upstream pressure was reduced to 2.1 MPa and the injection pump was switched to maintain this at a constant value. As such, the test history is non-ideal, in that the applied gas pressure was allowed to briefly exceed the confining pressure (for a period of 3 hours) by a maximum value of 0.27 MPa. However, any impact on sample behaviour is unlikely for a number of reasons: (1) a simple force balance analysis demonstrates that the force applied through the area of the injection filter during this time is negligible in comparison to the force applied to the sample by the confining pressure transferred by the injection platen; (2) the sample sheath was locked into position, providing some additional resistance to tensile strain; (3) there was no evidence for inflow into the sample at this time (or for the following three days), as might be expected if the sample was no longer gas-tight.

Following this event, constant gas pressure was maintained for an additional three and a half days, at which time (day 9.3) a notable increase in inflow was detected, indicating gas entry into the specimen. This inflow continued until day twelve when the pump ran out of fluid and the test was terminated. Such ‘time-dependent’ breakthrough behaviour has previously been observed in other argillaceous materials (Harrington et al., 2012) and it has been suggested that this ‘delay’ is related to the requirement for the generation of a sufficient network of dilatant pathways prior to significant gas migration.

Upon test completion, the core was examined using a combination of optical and electron microscopy for evidence of gold and TiO2 nanoparticle penetration. This was facilitated by submerging the sample in liquid nitrogen, thereby embrittling the clay. It was then split axially (on its side to prevent contamination), providing a longitudinal profile through the core, perpendicular to the gas injection point. The freshly fractured surfaces were photographed and visually examined under a binocular optical microscope. Higher resolution observations were then undertaken using a FEI QUANTA 600 environmental scanning electron microscope (ESEM) equipped with an Oxford Instruments INCA Energy 450 energy-dispersive X-ray microanalysis (EDXA) system with a 50 mm2 X-Max silicon-drift detector. The fracture surfaces were examined directly in the ESEM under low vacuum (0.98 Torr H2O), without a conductive coating in both secondary (SE) and backscattered electron imaging (BSE) modes, using a 10 kV beam. The EDXA was used to aid nanoparticle identification and to map the distribution of gold and TiO2 particles on the fractured clay surface.

Experimental observations

Flow observations and post-test sample imaging indicate that a network of conductive gas pathways was successfully established across the specimen during testing. Flux across the sample was maintained for sufficient time to promote the ingress of nanoparticles. Upon completion of the injection stage the sample was removed from the apparatus. Visual and optical microscope observations showed that the freshly-fractured core surface displayed a very diffuse, pale bluish white ‘bloom’, extending from immediately beneath the nanoparticle injection cake to a distance of up to 20 mm into the sample, as a ‘tongue-like’ plume with a generally narrowing tail running approximately parallel to the core axis (Fig. 3a). Small diffuse lobes were also observed to spread laterally from the margins of the main plume, parallel to the bedding lamination (the core axis was perpendicular to bedding lamination). Elemental mapping showed that the plume of bluish white discolouration observed on the fractured core surface corresponded precisely with areas of high TiO2 concentration (Fig. 3b). High-resolution SEM observations showed that the fractured clay surface in these areas was coated with a very fine, diffuse film of TiO2 nanoparticles combined with large ‘fluffy’ aggregates of TiO2 nanoparticles (Fig. 4a). The EDXA elemental mapping also showed that the TiO2 nanoparticle-aggregates were concentrated along a series of small, interconnected, sub-horizontal, braided, ‘step-like’ features or channels in the broken core surface. These correspond to where the core has broken along a network of very fine, anastomosing sub-vertical, laterally discontinuous microfractures (orientated sub-parallel to the core axis) and their intersection with the bedding lamination.

In contrast to the TiO2 nanoparticles, elemental mapping showed that gold nanoparticles were only very sparsely distributed (and barely detectable) on the fractured core surface, except within the nanoparticle cake and where relatively large aggregates of gold nanoparticles had penetrated into sandy laminae just beneath the point of gas injection. Although improved dispersion of the gold nanoparticles was achieved through using a mixture of gold and TiO2 nanoparticles, the analyses revealed that the gold nanoparticles still clumped together in much coarser aggregates than the TiO2 nanoparticles. The gold nanoparticle aggregates were usually seen to be enclosed within a mass of finer TiO2 nanoparticles on the fractured clay surface. However, nanoparticle aggregates and isolated nanoparticles of gold were also seen to be trapped and enclosed within the clay matrix of the rock (Fig. 4b). The gold nanoparticle aggregates typically varied between <0.5 and 4 μm in diameter, but close to the gas injection interface reached up to 10 μm in size. The aggregate diameter generally decreased with increasing distance from the gas injection interface, but traces of the gold nanoparticles could be found up to 20 mm from the point of gas injection.

Linear ‘trails’ (up to 300 μm long) of gold and TiO2 nanoparticles were observed trapped within the clay matrix in the fractured clay surface (Figs 5 and 6). These appear to delineate the probable traces of sealed pathways intersecting with the freshly broken core surface being observed under the ESEM. Larger open microfractures were also observed, usually developed parallel to bedding. These are quite different to the healed microfractures, and were largely devoid of gold and TiO2 nanoparticles (Fig. 6), and attributed to post-experimental artefacts produced by sample shrinkage during drying in the ESEM instrument. Similar nanoparticle aggregates were rarely seen to be present along the grain boundaries of quartz and feldspar silt particles (Fig. 7). These observations indicate that the nanoparticles have been injected along transient pathways, through the clay matrix and around grain boundaries of coarser detrital siliciclastic silt and sand grains, which were induced and dilated during gas injection. It would appear that these features then sealed, trapping any nanoparticles that were being transported in the gas stream. The size of the nanoparticle aggregates, typically between 0.5 and 4 μm (though up to 10 μm diameter), provides an indication of the minimum aperture of the dilatant pathways induced in the clay during gas injection.


This study has successfully demonstrated that gold and TiO2 nanoparticles can be introduced into a gas stream and used as a tracer for studying gas transport pathways through clay or mudrock. This is the first time that the location of gas pathways within low permeability mudrock or clays have been directly observed, and provides a new insight into the mechanism of gas migration through these rocks. These observations highlight the potential importance of dilatancy as a fundamental process in the advective movement of gas through such materials. These observations support previous indirect measurements which indicated the existence of such features (Horseman and Harrington, 1994; Ortiz et al., 1997, 2002; Rodwell, 2000; Harrington and Horseman, 1999; Van Geet et al., 2008; Angeli et al., 2009). The SEM observations show that relatively large aggregates of gold nanoparticles, generally between 0.5 to 4 μm in diameter, and up to 10 μm, have been injected with the gas stream, migrating through pressure induced pathways for distances up to 20 mm from the point of injection. This continuous network of pathways is approximately two orders of magnitude larger than the average pathways dimensions for undisturbed clay based on mercury intrusion porosimetry (Dehandschutter et al., 2004). The size distribution of aggregates entrained within the clay matrix clearly indicates variable pathway apertures. Current understanding of such features is that they collapse as gas pressure decays to a critical threshold (Horseman and Harrington, 1994; Ortiz et al., 1996; Harrington and Horseman, 1999, 2003; Rodwell, 2000). Gold (and TiO2) nanoparticles are trapped within the clay matrix, forming inclusion trails. These ‘strands’ of particulates delineate the traces of formerly conductive pathways, which have then collapsed, providing conclusive evidence of pathway closure. These pathways appear to be ephemeral features opening and closing with the passage of gas. The capacity to differentiate process (i.e. viscocapillary flow vs. pathway dilation) is fundamental when attempting to accurately describe the multi-phase flow behaviour of argillaceous media, particularly when predicting the long-term behaviour of such materials used for the disposal of heat emitting radioactive waste or as caprocks for the disposal of CO2.


Studies on nanoparticle migration were funded as a component of the BGS Geosphere Containment and Development of Capability projects (part of NERC's National Capability Programme). This work is dedicated to Dr Steve Horseman, who pioneered research on the movement of gas through clay-based host rocks and engineered barriers. He sadly died before nanoparticle injection was available as a viable method for the direct imagining of dilatant pathways, which supports his long-held belief for the importance of such features in multi-phase flow in clay-rich materials. In addition, the authors wish to thank Norbert Maes and Geert Volckaert (Studiecentrum voor Kernenergie, SCK-CEN) for their long support of this research area and Ondraf Niras for the provision of test material and review comments from both organizations. This paper is published with the permission of the Director, British Geological Survey (NERC).

TAW stands for Tweede Algemene Waterpassing which is the Belgian national datum.
Defined as 273.15 K and 101.325 kPa.