Slow gas expansion in saturated natural porous media by gas injection and partitioning with non-aqueous phase liquids

The partitioning of volatile non-aqueous phase liquid (NAPL) compounds to a discontinuous gas phase can result in the expansion of that gas phase, and the resulting gas flow can significantly affect the mass transfer from NAPL source zones. This recently reported gas flow generated by the spontaneous expansion of a discontinuous gas phase has not been extensively characterized in the literature. This study measured the expansion rate of a single gas cluster in a 1.1 mm sand above a pool of trans-1,2-dichloroethene (tDCE) in small-scale flow cell experiments. To characterize the gas flow, gas injection experiments in three sizes of sand were conducted at very slow injection rates typical of gas flow rates produced by gas expansion due to NAPL partitioning. Gas cluster spontaneous expansion rates above a tDCE pool were found to be 0.34 ± 0.02 and 0.29 ± 0.01 mL/day in duplicate experiments, which is sufficiently slow to result in discontinuous gas flow in porous media with a grain size diameter greater than 0.02 mm. Measured capillary pressures during gas injection showed patterns consistent with discontinuous gas flow, and identified multiple fragmentation events and expansion by coalescence with trapped clusters. The combination of pressure data and light transmission images were used to identify fragmentation and obtain direct measurements of the critical cluster length (i.e. the length at which withdrawal of the gas phase from a pore space occurs) in quasi-two-dimensional porous media for the first time. The measured critical cluster lengths were 1.4–3.6, 3.2–6.0 and 2.8–6.5 cm in 1.1, 0.7 and 0.5 mm sands, respectively. These values agreed well with estimates of the critical cluster length made using previously reported equations, and parameters derived from the medium’s capillary pressure-saturation relationship.     

Pore-scale imbibition experiments in dry and prewetted porous media

Two-phase imbibition behavior of immiscible fluids was studied in dry and prewetted porous media using a laser-induced fluorescence technique. Imbibition was first investigated in two-dimensional (2-D) systems under conditions comparable to those for a study of drainage [Ovdat H, Berkowitz B. Pore-scale study of drainage displacement under combined capillary and gravity effects in index-matched porous media. Water Resources Research 2006;42:W06411. doi:10.1029/2005WR004553] in the capillary-dominated regime. The effect of initial wetting saturation (IWS) was then explored in 2-D and 3-D porous media under the combined effect of gravity, capillary and viscous forces, within and outside the capillary-dominated regime. Parameters that describe maximum vertical advance, volumetric fraction, total surface area and specific surface area of the invading fluid were used to quantify the behavior. Comparison of 2-D drainage and imbibition patterns demonstrates significant qualitative differences under analogous viscosity ratio, buoyancy number, and capillary number values. However, quantitative analyses show strong pore-scale similarities between these patterns. Invasion structures in 3-D, prewetted (IWS ≈ 8% of the pore volume) porous media are ramified, with lateral branching and regions containing trapped residual fluid. These structures are qualitatively and quantitatively different from the compact, branchless structures that develop in dry (IWS = 0) porous media.     

Micro-scale experimental investigation of the effect of flow rate on trapping in sandstone and carbonate rock samples

We present the results of a pore-scale experimental study of residual trapping in consolidated sandstone and carbonate rock samples under confining stress. We investigate how the changes in wetting phase flow rate impacts pore-scale distribution of fluids during imbibition in natural, water-wet porous media. We systematically study pore-scale trapping of the nonwetting phase as well as size and distribution of its disconnected globules. Seven sets of drainage-imbibition experiments were performed with brine and oil as the wetting and nonwetting phases, respectively. We utilized a two-phase miniature core-flooding apparatus integrated with an X-ray microtomography system to examine pore-scale fluid distributions in small Bentheimer sandstone (D = 4.9 mm and L = 13 mm) and Gambier limestone (D = 4.4 mm and L = 75 mm) core samples. The results show that with increase in capillary number, the residual oil saturation at the end of the imbibition reduces from 0.46 to 0.20 in Bemtheimer sandstone and from 0.46 to 0.28 in Gambier limestone. We use pore-scale displacement mechanisms, in-situ wettability characteristics, and pore size distribution information to explain the observed capillary desaturation trends. The reduction was believed to be caused by alteration of the order in which pore-scale displacements took place during imbibition. Furthermore, increase in capillary number produced significantly different pore-scale fluid distributions during imbibition. We explored the pore fluid occupancies and studied size and distribution of the trapped oil clusters during different imbibition experiments. The results clearly show that as the capillary number increases, imbibition produces smaller trapped oil globules. In other words, the volume of individual trapped oil globules decreased at higher brine flow rates. Finally, we observed that the pore space in the limestone sample was considerably altered through matrix dissolution at extremely high brine flow rates. This increased the sample porosity from 44% to 62% and permeability from 7.3 D to 80 D. Imbibition in the altered pore space produced lower residual oil saturation (from 0.28 to 0.22) and significantly different distribution of trapped oil globules.     

Chaotic advection at the pore scale: Mechanisms,upscaling and implications for macroscopic transport

The macroscopic spreading and mixing of solute plumes in saturated porous media is ultimately controlled by processes operating at the pore scale. Whilst the conventional picture of pore-scale mechanical dispersion and molecular diffusion leading to persistent hydrodynamic dispersion is well accepted, this paradigm is inherently two-dimensional (2D) in nature and neglects important three-dimensional (3D) phenomena. We discuss how the kinematics of steady 3D flow at the pore scale generate chaotic advection—involving exponential stretching and folding of fluid elements—the mechanisms by which it arises and implications of microscopic chaos for macroscopic dispersion and mixing. Prohibited in steady 2D flow due to topological constraints, these phenomena are ubiquitous due to the topological complexity inherent to all 3D porous media. Consequently 3D porous media flows generate profoundly different fluid deformation and mixing processes to those of 2D flow. The interplay of chaotic advection and broad transit time distributions can be incorporated into a continuous-time random walk (CTRW) framework to predict macroscopic solute mixing and spreading. We show how these results may be generalised to real porous architectures via a CTRW model of fluid deformation, leading to stochastic models of macroscopic dispersion and mixing which both honour the pore-scale kinematics and are directly conditioned on the pore-scale architecture.     

A transport phase diagram for pore-level correlated porous media

Transport in porous media is often characterized by the advection–dispersion equation, with the dispersion coefficient as the most important parameter that links the hydrodynamics to the transport processes. Morphological properties of any porous medium, such as pore size distribution, network topology, and correlation length control transport. In this study we explore the impact of correlation length on transport regime using pore-network modelling. Earlier direct simulation studies of dispersion in carbonate and sandstone rocks showed larger dispersion compared to granular homogenous sandpacks. However, in these studies, isolation of the impact of correlation length on transport regime was not possible due to the fundamentally different pore morphologies and pore-size distributions. Against this limitation, we simulate advection–dispersion transport for a wide range of Péclet numbers in unstructured irregular networks with “different” correlation lengths but “identical” pore size distributions and pore morphologies. Our simulation results show an increase in the magnitudes of the estimated dispersion coefficients in correlated networks compared to uncorrelated ones in the advection-controlled regime. The range of the Péclet numbers which dictate mixed advection–diffusion regime considerably reduces in the correlated networks. The findings emphasize the critical role of correlation length which is depicted in a conceptual transport phase diagram and the importance of accounting for the micro-scale correlation lengths into predictive stochastic pore-scale modelling.     

Tracer dispersal in sandy sediment porewater under enhanced physical forcing

This research (1) integrated a fluorescent dye injection and monitoring system for measuring the mixing of a fluorescent dye tracer (fluorescein) in permeable (sandy) sediments with a cabled ocean observatory, Kilo Nalu, Oahu, Hawaii, and (2) used this system to conduct remotely controlled in situ measurements of wave-enhanced porewater mixing in a physically well-characterized wave-dominated setting. Laboratory results indicated that the fiber-optic sensor is effective at measuring fluorescence-traced enhanced mixing in sandy sediments. Observed dye mixing, driven by paddle-generated surface waves in a laboratory wave channel was 2–3 orders of magnitude greater than molecular diffusion, and decreased with depth in the sediments. Similarly, field experiments used a remotely controlled submersible syringe pump for fluorescent tracer injection into sediment that was monitored with a fiber-optic sensor. These experiments were carried out at 10 m water depth, with mean wave heights of 0.3–1.5 m and periods of ∼15 s. At 15 cm sediment depth, transport rates of 0–5 cm h⁻¹ were measured, with maximum dispersion coefficients 2–3 orders of magnitude faster than sedimentary molecular diffusion. Hydrodynamic measurements collected simultaneously via Kilo Nalu correlated with porewater transport, with significant wave height and diffusion having the strongest covariation.     

Transport,sloughing and settling rates of estuarine macrophytes: Mechanisms and ecological implications

The study of plant-bound nutrient transport has been largely neglected in estuaries. Lately however, it has been shown that nutrients bound to macroalgae and seagrasses can constitute a major part of the nutrient transport in shallow tidal estuaries. Organic detritus in estuaries comes from various sources. This paper looks into the source of detritus from sloughing, and the transport behaviour of plant detritus under unidirectional flows. In order to determine the extent of the sloughing of macrophytes, the threshold current velocities for the traction and resuspension of the most common submersed macrophytes in Venice Lagoon (Ulva lactuca, Enteromorpha sp., Ceramium rubrum, Cladophora sp., and Chaetomorpha linum) were studied in a laboratory flume. It was found that all macrophytes subjected to flows of 1.5–3.0 cm s⁻¹ move initially as bed load. The threshold for suspension of the macrophyte tissue was at current speeds >3 cm s⁻¹. The exception was the filamentous macroalgae, C. linum, which moved as bed load at all current speeds. This implies that the advection of plant-bound nutrients in Venice Lagoon is widespread and takes place over virtually all stages of the tide.     

Rapid injection of particles and gas into non-fluidized granular material, and some volcanological implications

In diatremes and other volcanic vents, steep bodies of volcaniclastic material having differing properties (particle size distribution, proportion of lithic fragments, etc.) from those of the surrounding vent-filling volcaniclastic material are often found. It has been proposed that cylindrical or cone-shaped bodies result from the passage of “debris jets” generated after phreatomagmatic explosions or other discrete subterranean bursts. To learn more about such phenomena, we model experimentally the injection of gas-particulate dispersions through other particles. Analogue materials (glass beads or sand) and a finite amount of compressed air are used in the laboratory. The gas is made available by rapidly opening a valve—therefore the injection of gas and coloured particles into a granular host is a brief (<1 s), discrete event, comparable to what occurs in nature following subterranean explosions. The injection assumes a bubble shape while expanding and propagating upwards. In reaction, the upper part of the clastic host moves upward and outward above the ‘bubble’, forming a ‘dome’. The doming effect is much more pronounced for shallow injection depths (thin hosts), with dome angles reaching more than 45°. Significant surface doming is also observed for some full-scale subterranean blasts (e.g. buried nuclear explosions), so it is not an artefact of our setup. What happens next in the experiments depends on the depth of injection and the nature of the host material. With shallow injection into a permeable host (glass beads), the compressed air in the “bubble’ is able to diffuse rapidly through the roof. Meanwhile the coloured beads sediment into the transient cavity, which is also closing laterally because of inward-directed granular flow of the host. Depending on the initial gas pressure in the reservoir, the two-phase flow can “erupt” or not; non-erupting injections produce cylindrical bodies of coloured beads whereas erupting runs produce flaring upward or conical deposits. Changing the particle size of the host glass beads does not have a large effect under the size range investigated (100–200 to 300–400 μm). Doubling the host thickness (injection depth) requires a doubling of the initial gas pressure to produce similar phenomena. Such injections—whether erupting or wholly subterranean—provide a compelling explanation for the origin and characteristics of multiple cross-cutting bodies that have been documented for diatreme and other vent deposits.     

Micro-computed tomography pore-scale study of flow in porous media: Effect of voxel resolution

A fundamental understanding of flow in porous media at the pore-scale is necessary to be able to upscale average displacement processes from core to reservoir scale. The study of fluid flow in porous media at the pore-scale consists of two key procedures: Imaging - reconstruction of three-dimensional (3D) pore space images; and modelling such as with single and two-phase flow simulations with Lattice-Boltzmann (LB) or Pore-Network (PN) Modelling. Here we analyse pore-scale results to predict petrophysical properties such as porosity, single-phase permeability and multi-phase properties at different length scales. The fundamental issue is to understand the image resolution dependency of transport properties, in order to up-scale the flow physics from pore to core scale. In this work, we use a high resolution micro-computed tomography (micro-CT) scanner to image and reconstruct three dimensional pore-scale images of five sandstones (Bentheimer, Berea, Clashach, Doddington and Stainton) and five complex carbonates (Ketton, Estaillades, Middle Eastern sample 3, Middle Eastern sample 5 and Indiana Limestone 1) at four different voxel resolutions (4.4 µm, 6.2 µm, 8.3 µm and 10.2 µm), scanning the same physical field of view. Implementing three phase segmentation (macro-pore phase, intermediate phase and grain phase) on pore-scale images helps to understand the importance of connected macro-porosity in the fluid flow for the samples studied. We then compute the petrophysical properties for all the samples using PN and LB simulations in order to study the influence of voxel resolution on petrophysical properties. We then introduce a numerical coarsening scheme which is used to coarsen a high voxel resolution image (4.4 µm) to lower resolutions (6.2 µm, 8.3 µm and 10.2 µm) and study the impact of coarsening data on macroscopic and multi-phase properties. Numerical coarsening of high resolution data is found to be superior to using a lower resolution scan because it avoids the problem of partial volume effects and reduces the scaling effect by preserving the pore-space properties influencing the transport properties. This is evidently compared in this study by predicting several pore network properties such as number of pores and throats, average pore and throat radius and coordination number for both scan based analysis and numerical coarsened data.     

Comparing the behaviour of spherical beads and natural grains in bedload mixtures

It is common to use idealised materials to study the dynamics of granular transport in fluid flows, but the impact of this choice upon sediment behaviour has not been extensively explored. To tackle this research gap, two experiments were undertaken to explore the influence of a finer grain input to a channelized coarser granular flow driven by a shallow fluid flow. The first set of runs was undertaken using spherical glass beads, and the second set with natural fluvial sediment. The transport system approximates a narrow slice through the bedload at the bottom of a river. In the runs with natural grains, the infiltration of fine sediment into the bed was similar to the spherical glass beads, but with reduced infiltration capacity. We ascribe this behaviour to irregular and variable pore shapes and sizes in the natural material. The behaviour of the bedload in the natural material runs matched that of the bead runs only when the feed contained a high content of fines. When the feed contained a low content of fines the transport of natural grains was more complex, including the emergence of migrating collections of grains. However, the overall changes in bed and water slope due to the finer grain input were comparable in both sets of experiments. We conclude that artificial, idealised materials qualitatively represent sedimentary phenomena, but that quantitative differences in the outcomes must be expected. © 2020 John Wiley & Sons, Ltd.     

Grid convergence of variable-density flow simulations in discretely-fractured porous media

Grid convergence in space and time of variable-density flow in fractured-porous rock is systematically assessed. Convergence of the flow simulation is attained using both uniform and adaptive time-stepping. This contrasts to variable-density flow in unfractured porous media where grid convergence variable-density flow problems is almost never achieved. At high discretization levels, the number of fingers in fractured-porous rock is no longer influenced by spatial-temporal grid discretization, which is not the case in unfractured porous media. However, similar to unfractured porous media, the number of fingers in fractured-porous media varies at low discretization levels. Simulated convective pattern and penetration depth of the dense plume in fractured rock depend more on spatial discretization than on temporal discretization. The appropriate spatial-temporal grid is then used to examine some aspects of mixed convection in fractured-porous rock, characterized by the mixed convection number M. The critical mixed convection number M_c = 46 represents the transition between forced and free convection in fractured porous media, which is much higher than M_c = 1 in unfractured porous media. Thus, for mixed convective flow problems, the value of M_c is not a sufficient indicator to predict the convective mode (free convection-forced convection), and the presence of vertical fractures must be included in the prediction of convective flow modes.     

Direct simulations of two-phase flow on micro-CT images of porous media and upscaling of pore-scale forces

Pore-scale forces have a significant effect on the macroscopic behaviour of multiphase flow through porous media. This paper studies the effect of these forces using a new volume-of-fluid based finite volume method developed for simulating two-phase flow directly on micro-CT images of porous media. An analytical analysis of the relationship between the pore-scale forces and the Darcy-scale pressure drops is presented. We use this analysis to propose unambiguous definitions of Darcy-scale viscous pressure drops as the rate of energy dissipation per unit flow rate of each phase, and then use them to obtain the relative permeability curves. We show that this definition is consistent with conventional laboratory/field measurements by comparing our predictions with experimental relative permeability. We present single and two-phase flow simulations for primary oil injection followed by water injection on a sandpack and a Berea sandstone. The two-phase flow simulations are presented at different capillary numbers which cover the transition from capillary fingering at low capillary numbers to a more viscous fingering displacement pattern at higher capillary numbers, and the effect of capillary number on the relative permeability curves is investigated. Overall, this paper presents a new finite volume-based methodology for the detailed analysis of two-phase flow directly on micro-CT images of porous media and upscaling of the results to the Darcy scale.     

Dissolution of residual non-aqueous phase liquids in porous media: pore-scale mechanisms and mass transfer rates

Experiments designed to elucidate the pore-scale mechanisms of the dissolution of a residual non-aqueous phase liquid (NAPL), trapped in the form of ganglia within a porous medium, are discussed. These experiments were conducted using transparent glass micromodels with controlled pore geometry, so that the evolution of the size and shape of individual NAPL ganglia and, hence, the pore-scale mass transfer rates and mass transfer coefficients could be determined by image analysis. The micromodel design permitted reasonably accurate control of the pore water velocity, so that the mass transfer coefficients could be correlated in terms of a local (pore-scale) Peclet number. A simple mathematical model, incorporating convection and diffusion in a slit geometry was developed and used successfully to predict the observed mass transfer rates. For the case of non-wetting NAPL ganglia, water flow through the corners in the pore walls was seen to control the rate of NAPL dissolution, as recently postulated by Dillard and Blunt [Water Resour. Res. 36 (2000) 439–454]. Break-up of doublet non-wetting phase ganglia into singlet ganglia by snap-off in pore throats was also observed, confirming the interplay between capillarity and mass transfer. Additionally, the effect of wettability on dissolution mass transfer was demonstrated. Under conditions of preferential NAPL wettability, mass transfer from NAPL films covering the solid surfaces was seen to control the dissolution process. Supply of NAPL from the trapped ganglia to these films by capillary flow along pore corners was observed to result in a sequence of pore drainage events that increase the interfacial area for mass transfer. These observations provide new experimental evidence for the role of capillarity, wettability and corner flow on NAPL ganglia dissolution.     

Laboratory Study of Air Sparging: Air Flow Visualization

Laboratory flow visualization experiments, using glass beads as the porous medium, were conducted to study air sparging, an innovative technology for subsurface contaminant remediation. The purpose of these experiments was to observe how air flows through saturated porous media and to obtain a basic understanding of air plume formation and medium heterogeneity effects. The experiments indicate that air flow occurring in discrete, stable channels is the most probable flow behavior in medium to fine grained water saturated porous media and that medium heterogeneity plays an important role in the development of air channels. Several simulated scales of heterogeneities, from pore to field, have been studied. The results suggest that air channel formation is sensitive to the various scales of heterogeneities. Site-specific hydrogeologic settings have to be carefully reviewed before air sparging is applied to remediate sites contaminated by volatile organic compounds.     

Scouring of a granular bed by dam break: Experimental study and numerical simulation by a VOF-LPT coupling

The first part of the research reported here consists of an experimental campaign to study the scouring of a granular bed(glass beads, sand) induced by a dam break in an open channel. Two configurations are considered: with and without cylinders. In the second part of this study, the volume of fluid method coupled with the shear stress transport turbulent model and the lagrangian particle tracking method is used to simulate the local scour processes. The four-way coupling is realized by consider...     

Inviscid behaviour of fines-rich pyroclastic flows inferred from experiments on gas-particle mixtures

Experiments were carried out on granular flows generated by instantaneous release of gas-fluidised, bidisperse mixtures and propagating into a horizontal channel. The mixture consists of fine (< 100 μm) and coarse (> 100 μm) particles of same density, with corresponding grain size ratios of ∼ 2 to 9. Initial fluidisation of the mixture destroys the interparticle frictional contacts, and the flow behaviour then depends on the initial bed packing and on the timescale required to re-establish strong frictional contacts. At a fines mass fraction (α) below that of optimal packing (∼ 40%), the initial mixtures consist of a continuous network of coarse particles with fines in interstitial voids. Strong frictional contacts between the coarse particles are probably rapidly re-established and the flows steadily decelerate. Some internal friction reduction appears to occur as α and the grain size ratio increases, possibly due to particle rolling and the lower roughness of internal shear surfaces. Segregation only occurs at large grain size ratio due to dynamical sieving with fines concentrated at the flow base. In contrast, at α above that for optimal packing, the initial mixtures consist of coarse particles embedded in a matrix of fines. Flow velocities and run-outs are similar to that of the monodisperse fine end-member, thus showing that the coarse particles are transported passively within the matrix whatever their amount and grain size are. These flows propagate at constant height and velocity as inviscid fluid gravity currents, thus suggesting negligible interparticle friction. We have determined a Froude number of 2.61 ± 0.08 consistent with the dam-break model for fluid flows, and with no significant variation as a function of α, the grain size ratio, and the initial bed expansion. Very little segregation occurs, which suggests low intensity particle interactions during flow propagation and that active fluidisation is not taking place. Strong frictional contacts are only re-established in the final stages of emplacement and stop the flow motion. We infer that fines-rich (i.e. matrix-supported) pyroclastic flows propagate as inviscid fluid gravity currents for most of their emplacement, and this is consistent with some field data.     

Solute dilution under imbibition and drainage conditions in a heterogeneous structure: Modeling of a sand tank experiment

This study aims at modeling the transport of a conservative tracer in two dimensions, as experimentally observed in a strongly heterogeneous medium under conditions of variable water saturation during drainage and imbibition. Solute transport experiments were conducted in a sand tank containing an artificial packing of three quartz sands of different particle sizes. The packing was characterized by the presence of numerous homogeneous layers (0.5 × 5 × 5 cm) inclined at 45° and randomly distributed in a tank. Six different stationary flow conditions were sequentially established during imbibition and drainage. When a stationary flow regime was reached, several solute pulses were applied at different positions at the upper surface of the sand structure. The transport regime was studied by monitoring the tracer plumes injected as point-like pulses at the surface, as they travelled through the sand bedding.     

Exploring the nature of non-Fickian transport in laboratory experiments

Observations of non-Fickian transport in sandbox experiments [Levy M, Berkowitz B. Measurement and analysis of non-Fickian dispersion in heterogeneous porous media. J Contam Hydrol 2003;64:203–26] were analyzed previously using a power law tail ψ(t) ∼ t^−1−β with 0 < β < 2 for the spectrum of transition times comprising a tracer plume migration. For each sandbox medium a choice of β resulted in an excellent fit to the breakthrough curve (BTC) data, and the value of β decreased slowly with increasing flow velocity. Here, the data are reanalyzed with the full spectrum of ψ(t) gleaned from analytical calculations [Cortis A, Chen Y, Scher H, Berkowitz B. Quantitative characterization of pore-scale disorder effects on transport in “homogeneous” granular media. Phys Rev E 2004;10(70):041108. doi: 10.1103/PhysRevE.70.041108], numerical simulations [Bijeljic B, Blunt MJ. Pore-scale modeling and continuous time random walk analysis of dispersion in porous media. Water Resour Res 2006;42:W01202. doi: 10.1029/2005WR004578] and permeability fields [Di Donato G, Obi E-O, Blunt MJ. Anomalous transport in heterogeneous media demonstrated by streamline-based simulation. Geophys Res Lett 2003;30:1608–12s. doi: 10.1029/2003GL017196]. We represent the main features of the full spectrum of transition times with a truncated power law (TPL), ψ(t) ∼ (t₁ + t)^−1−βexp(−t/t₂), where t₁ and t₂ are the limits of the power law spectrum. An excellent fit to the entire BTC data set, including the changes in flow velocity, for each sandbox medium is obtained with a single set of values of t₁, β, t₂. The influence of the cutoff time t₂ is apparent even in the regime t < t₂. Significantly, we demonstrate that the previous apparent velocity dependence of β is a result of choosing a pure power law tail for ψ(t). The key is the change in the log–log slope of the TPL form of ψ(t) with a shifting observational time window caused by the change in the mean velocity. Hence, the use of the full spectrum of ψ(t) is not only necessary for the transition to Fickian behavior, but also to account for the dynamics of these laboratory observations of non-Fickian transport.     

Thellier palaeointensity experiments on Faroes flood basalts: technical aspects and geomagnetic implications

Based on palaeomagnetic and rock magnetic results, 108 samples from 27 flows from Palaeogene flood basalts of Faroe Islands were chosen for whole-rock Thellier palaeointensity experiments. Altogether 90 samples were rejected due to either chemical alterations or typical multidomain (MD) behaviour evidenced by pTRM-tails. AF pre-treatment was used to reduce the effect of MD grains on Thellier experiments. Only five flows (18 samples) yielded acceptable palaeointensity estimates, with flow mean VDMs ranging from 3.5 to 7.4×10²² A m². Modest selection criteria imposed on all published 5-160 Ma palaeointensity data left only 15 palaeomagnetic dipole moments: eight from whole-rock samples, six from submarine basaltic glass and one from single plagioclase crystals. More data are needed before the intriguing differences between results from different materials can be put into a geomagnetic context.