Calibration, Conversion, and Quantitative Multi-Layer Inversion of Multi-Coil Rigid-Boom Electromagnetic Induction Data.

Multi-coil electromagnetic induction (EMI) systems induce magnetic fields below and above the subsurface. The resulting magnetic field is measured at multiple coils increasingly separated from the transmitter in a rigid boom. This field relates to the subsurface apparent electrical conductivity (σa), and σa represents an average value for the depth range investigated with a specific coil separation and orientation. Multi-coil EMI data can be inverted to obtain layered bulk electrical conductivity models. However, above-ground stationary influences alter the signal and the inversion results can be unreliable. This study proposes an improved data processing chain, including EMI data calibration, conversion, and inversion. For the calibration of σa, three direct current resistivity techniques are compared: Electrical resistivity tomography with Dipole-Dipole and Schlumberger electrode arrays and vertical electrical soundings. All three methods obtained robust calibration results. The Dipole-Dipole-based calibration proved stable upon testing on different soil types. To further improve accuracy, we propose a non-linear exact EMI conversion to convert the magnetic field to σa. The complete processing workflow provides accurate and quantitative EMI data and the inversions reliable estimates of the intrinsic electrical conductivities. This improves the ability to combine EMI with, e.g., remote sensing, and the use of EMI for monitoring purposes.

Methane gas transport in unconfined aquifers: A numerical sensitivity study of a controlled release experiment at CFB Borden.

Subsurface leakage of natural gas from petroleum wells can impact freshwater aquifers. Accurate prediction of gas migration in the subsurface will depend on knowledge of permeability, porosity, and flow system conditions. A series of two-dimensional numerical multi-phase flow simulations (CFbio) were conducted to investigate the role of multi-phase parameters (relative permeability and air entry pressure), flow system conditions (intrinsic permeability, anisotropy, and groundwater velocity), and geometric properties (layer thickness and layer lateral continuity) on the flow of gas-phase methane emanating from two variable-rate point sources in an unconfined sandy aquifer. Numerical simulations showed that for a homogeneous, weakly anisotropic aquifer, gas migrates almost exclusively vertically due to buoyancy, before venting to the vadose zone and atmosphere. As vertical migration became restricted through increased anisotropy, inclusion of lower-permeable layers, or increased horizontal groundwater velocity, an increase in the lateral component of gas migration was observed. This led to the formation of a broader lateral migration of the gas-phase plume and establishment of variably distributed vertical preferential flow paths, ultimately resulting in increased gas retention in the aquifer with relatively less methane reaching the vadose zone or atmosphere. The inclusion of a thin layer with moderately lower permeability (1-2 orders of magnitude) and increased air entry pressure was used to depict a fine-grained sand lens within a uniform aquifer. This subtle feature led to the formation of thin gas pools extending up- and down-gradient beneath the lens, allowing methane to travel much farther and faster than by groundwater advection alone, which is consistent with field observations during the experiment. In all scenarios investigated gas-phase methane was shown to migrate predominantly vertically due to buoyancy, until the aquitard permeability was <30% of the aquifer permeability. Our modelling demonstrates that even subtle permeability contrasts, together with capillary pressure changes demarcating grain-scale bedding, will lead to extensive lateral free-phase gas migration, and development of a more extensive and complex zone of impacted aquifer than presupposed.

Monitoring the evolution and migration of a methane gas plume in an unconfined sandy aquifer using time-lapse GPR and ERT.

Fugitive methane (CH4) leakage associated with conventional and unconventional petroleum development (e.g., shale gas) may pose significant risks to shallow groundwater. While the potential threat of stray (CH4) gas in aquifers has been acknowledged, few studies have examined the nature of its migration and fate in a shallow groundwater flow system. This study examines the geophysical responses observed from surface during a 72day field-scale simulated CH4 leak in an unconfined sandy aquifer at Canadian Forces Base Borden, Canada, to better understand the transient behaviour of fugitive CH4 gas in the subsurface. Time-lapse ground-penetrating radar (GPR) and electrical resistivity tomography (ERT) were used to monitor the distribution and migration of the gas-phase and assess any impacts to groundwater hydrochemistry. Geophysical measurements captured the transient formation of a CH4 gas plume emanating from the injector, which was accompanied by an increase in total dissolved gas pressure (PTDG). Subsequent reductions in PTDG were accompanied by reduced bulk resistivity around the injector along with an increase in the GPR reflectivity along horizontal bedding reflectors farther downgradient. Repeat temporal GPR reflection profiling identified three events with major peaks in reflectivity, interpreted to represent episodic lateral CH4 gas release events into the aquifer. Here, a gradual increase in PTDG near the injector caused a sudden lateral breakthrough of gas in the direction of groundwater flow, causing free-phase CH4 to migrate much farther than anticipated based on groundwater advection. CH4 accumulated along subtle permeability boundaries demarcated by grain-scale bedding within the aquifer characteristic of numerous Borden-aquifer multi-phase flow experiments. Diminishing reflectivity over a period of days to weeks suggests buoyancy-driven migration to the vadose zone and/or CH4 dissolution into groundwater. Lateral and vertical CH4 migration was primarily governed by subtle, yet measurable heterogeneity and anisotropy in the aquifer.

Estimating water content in an active landfill with the aid of GPR.

Landfill gas (LFG) receives a great deal of attention due to both negative and positive environmental impacts, global warming and a green energy source, respectively. However, predicting the quantity of LFG generated at a given landfill, whether active or closed is difficult due to the heterogeneities present in waste, and the lack of accurate in situ waste parameters like water content. Accordingly, ground penetrating radar (GPR) was evaluated as a tool for estimating in situ water content. Due to the large degree of subsurface heterogeneity and the electrically conductive clay cap covering landfills, both of which affect the transmission of the electromagnetic pulses, there is much scepticism concerning the use of GPR to quantify in situ water content within a municipal landfill. Two landfills were studied. The first landfill was used to develop the measurement protocols, while the second landfill provided a means of confirming these protocols. GPR measurements were initially completed using the surface GPR approach, but the lack of success led to the use of borehole (BH) GPR. Both zero offset profiling (ZOP) and multiple offset gathers (MOG) modes were tried, with the results indicating that BH GPR using the ZOP mode is the most simple and efficient method to measure in situ water content. The best results were obtained at a separation distance of 2m, where higher the water content, smaller the effective separation distance. However, an increase in water content did appear to increase the accuracy of the GPR measurements. For the effective separation distance of 2m at both landfills, the difference between GPR and lab measured water contents were reasonable at 33.9% for the drier landfill and 18.1% for the wetter landfill. Infiltration experiments also showed the potential to measure small increases in water content.

Complex dielectric response of ellipsoidal particles with surface conduction.

Both particle shape and surface phenomena significantly affect the effective complex dielectric properties of colloidal systems. The treatment of particle shape has generally relied on the extrapolation from the solution of the spherical case proposed by O'Konski [J. Chem. Phys. 64, 605 (1960)] that treats ellipsoidal particles possessing surface conductivity as equivalent homogeneous anisotropic ellipsoids with bulk conduction. To test this approach, we have performed a rigorous analysis of the complex dielectric response of an ellipsoidal particle with surface conductivity using the generalization of the O'Konski boundary conditions to an ellipsoidal shape. The resulting closed-form solution obtained shows that surface conduction effects are represented by an equivalent inhomogeneous anisotropic ellipsoid. For the case of a spheroidal particle, the principle axes of the effective dielectric permittivity tensor of the equivalent particle are aligned with its geometrical principal axes; the effective permittivity varies in the direction of the unique spheroidal axis. In addition, numerical results indicate that the product of the surface area to volume ratio and the specific surface conductivity completely characterizes the effect of the surface phenomena on the response of spheroidal particles with a given shape. Numerical simulations show that spherical and prolate spheroidal particles exhibit a progressive dielectric enhancement while more disklike oblate spheroidal particles undergo an initial dielectric suppression followed by a subsequent enhancement with increasing surface conduction. A comparison of our model predictions with those obtained using the O'Konski approximation revealed significant differences in the magnitude of the low-frequency dielectric enhancement and relaxation frequency for ellipsoidal particle suspensions.

Long-term ground penetrating radar monitoring of a small volume DNAPL release in a natural groundwater flow field.

An earlier field experiment at Canadian Forces Base Borden by Brewster and Annan [Geophysics 59 (1994) 1211] clearly demonstrated the capability of ground penetrating radar (GPR) reflection profiling to detect and monitor the formation of DNAPL layers in the subsurface. Their experiment involved a large volume release (770 L) of tetrachloroethylene into a portion of the sand aquifer that was hydraulically isolated from groundwater flow by sheet pile walls. In this study, we evaluated the ability of GPR profiling to detect and monitor much smaller volume releases (50 L). No subsurface confining structure was used in this experiment; hence, the DNAPL impacted zone was subjected to the natural groundwater flow regime. This condition allowed us to geophysically monitor the DNAPL mass loss over a 66 month period. Reflectivity variations on the GPR profiles were used to infer the presence and evolution of the solvent layers. GPR imaging found significant reflectivity increases due to solvent layer formation during the two week period immediately after the release. These results demonstrated the capacity of GPR profiling for the detection and monitoring of lesser volume DNAPL releases that are more representative of small-scale industrial spills. The GPR imaged solvent layers subsequently reduced in both areal extent and reflectivity after 29 months and almost completely disappeared by the end of the 66 month monitoring period. Total DNAPL mass estimates based on GPR profiling data indicated that the solvent mass was reduced to 34%-36% of its maximum value after 29 months; only 4%-9% of the solvent mass remained in the study area after 66 months. These results are consistent with independent hydrogeological estimates of remaining DNAPL mass based on the downgradient monitoring of the dissolved solvent phase. Hence, we have concluded that the long-term GPR reflectivity changes of the DNAPL layers are likely the result from the dissolution of chlorinated solvents residing in those layers. The long-term monitoring results demonstrated that GPR profiling is a promising non-invasive method for use at DNAPL contaminated sites in sandy aquifers where temporal information about immiscible contaminant mass depletion due to either natural flow or remediation is needed. However, our results also indicated that the GPR signature of older DNAPL impacted zones may not differ greatly from the uncontaminated background if significant mass reduction due to dissolution has occurred.