Philip I. Meldrum

Electrical resistivity tomography applied to geologic, hydrogeologic, and engineering investigations at a former waste-disposal site

A former dolerite quarry and landfill site was investigated using 2D and 3D electrical resistivity tomography (ERT), with the aims of determining buried quarry geometry, mapping bedrock contamination arising from the landfill, and characterizing site geology. Resistivity data were collected from a network of intersecting survey lines using a Wenner-based array configuration. Inversion of the data was carried out using 2D and 3D regularized least-squares optimization methods with robust (L1-norm) model constraints. For this site, where high resistivity contrasts were present, robust model constraints produced a more accurate recovery of subsurface structures when compared to the use of smooth (L2-norm) constraints. Integrated 3D spatial analysis of the ERT and conventional site investigation data proved in this case a highly effective means of characterizing the landfill and its environs. The 3D resistivity model was successfully used to confirm the position of the landfill boundaries, which appeared as electrically well-defined features that corresponded extremely closely to both historic maps and intrusive site investigation data. A potential zone of leachate migration from the landfill was identified from the electrical models; the location of this zone was consistent with the predicted direction of groundwater flow across the site. Unquarried areas of a dolerite sill were imaged as a resistive sheet-like feature, while the fault zone appeared in the 2D resistivity model as a dipping structure defined by contrasting bedrock resistivities.

Noninvasive monitoring of DNAPL migration through a saturated porous medium using electrical impedance tomography.

Electrical impedance tomography (EIT) was used to monitor the movement of a fluorinated hydrocarbon dense nonaqueous phase liquid (DNAPL) through a saturated porous medium within a laboratory column. Impedance measurements were made using a horizontal plane of 12 electrodes positioned at regular intervals around the centre of the column. A 2D inversion algorithm, which incorporated the cylindrical geometry of the column, was used to reconstruct resistivity and phase images from the measured data. Differential time-lapse images of DNAPL movement past the plane of electrodes were generated by the cell-by-cell subtraction of resistivity and phase baseline models from those associated with the DNAPL release stage of the experiment. The DNAPL pulse was clearly delineated as resistive anomalies in the differential time-lapse resistivity images. The spatial extent of the resistive anomalies indicated that in addition to vertical migration, some lateral spreading of the DNAPL had occurred. Residual contamination could be detected after quasi-static conditions were reestablished. Residual DNAPL saturation was estimated from the resistivity model data by applying Archies second equation.Despite significant measured phase changes due to DNAPL contamination, the differential phase images revealed only weak anomalies associated with DNAPL flow; these anomalies could be seen only in the initial stages of the experiment during peak flow through the plane of electrodes.

Automated monitoring of coastal aquifers with electrical resistivity tomography

An Automated time-Lapse Electrical Resistivity Tomography (ALERT) system has been developed for the long-term monitoring of coastal aquifers. This ALERT system has been permanently installed in the River Andarax, Almeria, Spain to monitor and manage the impact of climatic change and land-use practice on the underlying Quaternary aquifer. An electrode array, nearly 1.6 km long, has been buried below the normally dry riverbed with electrode take-outs at regular intervals of 10 m. The maximum depth of investigation is about 160 m below ground level. An unmanned, permanent control station, in a secure location, allows the aquifer to be interrogated remotely from the BGS office in the UK. Volumetric geoelectric images of the subsurface can be obtained ‘on demand’ or at regular intervals; thereby eliminating the need for expensive repeat surveys. The entire process from data capture to image on the office PC is fully automated and seamless. The ALERT technology can provide early warning of potential threats to vulnerable water systems such as over-exploitation, rising sea levels, anthropogenic pollutants and seawater intrusion. The electrical images obtained (in space and time) are interpreted in terms of the hydrogeologic features including the seawater-freshwater interface. The timely detection and imaging of groundwater changes can help to regulate pumping and irrigation schemes.

Fundamentals of the capacitive resistivity technique

Capacitive resistivity (CR) is an emerging geophysical technique designed to extend the scope of the conventional methodology of dc resistivity to environments where galvanic coupling is notoriously difficult to achieve — for example, across engineered structures (roads, pavements), hard rock, dry soil, or frozen ground. Conceptually, CR is based on a four-point array capacitively coupled to the ground. Under certain conditions, capacitive measurements of resistivity are equivalent to those obtained with the dc technique, thus making dc interpretation schemes applicable to CR data. The coupling properties of practical sensor realizations are shown to be a function of their geometrical arrangement. Separate bodies of theory are associated with two complementary but distinct sensor types: the capacitive-line antenna and the plate-wire combination. The use of plate-wire combinations results in localized coupling, which, in conjunction with a quasi-static (low-frequency) formulation of the transfer impedance,...

4D electrical resistivity tomography monitoring of soil moisture dynamics in an operational railway embankment

The internal moisture dynamics of an aged (> 100 years old) railway earthwork embankment, which is still in use, are investigated using 2D and 3D resistivity monitoring. A methodology was employed that included automated 3D ERT data capture and telemetric transfer with on-site power generation, the correction of resistivity models for seasonal temperature changes and the translation of subsurface resistivity distributions into moisture content based on petrophysical relationships developed for the embankment material. Visualization of the data as 2D sections, 3D tomograms and time series plots for different zones of the embankment enabled the development of seasonal wetting fronts within the embankment to be monitored at a high-spatial resolution and the respective distributions of moisture in the flanks, crest and toes of the embankment to be assessed. Although the embankment considered here is at no immediate risk of failure, the approach developed for this study is equally applicable to other more high-risk earthworks and natural slopes.

The Use of 3D Electrical Resistivity Tomography to Characterise Waste and Leachate Distribution within a Closed Landfill, Thriplow, UK

Electrical Resistivity Tomography (ERT) has been used to map the 3D spatial distribution of waste and leachate concentrations within a closed and unconfined landfill. Borehole sampling and a 2D ERT survey down-gradient of the landfill boundary had failed to detect a pollution plume. Accordingly, a 3D survey was undertaken to determine the pattern of leachate drainage within the waste so that a more refined contaminant transport model could be developed. A full 3D survey was undertaken by sub-dividing the landfill into a number of discrete rectangular blocks and acquiring data on multiple parallel lines. The line data were merged into a single x-y matrix file and then inverted using a 3D finite element algorithm. The results are presented as 3D volumetric tomograms to show the inferred waste distribution and leachate flow-paths. The resistivity models indicate that leachate has accumulated at several discrete localities within the landfill. The controlling mechanism appears to be the depth and geometry of ...

Quantification of reactivated landslide behaviour using acoustic emission monitoring

Slope failures world-wide cause many thousands of deaths each year and damage built environment infrastructure, costing billions of pounds to repair, resulting in thousands of people being made homeless and the breakdown of basic services such as water supply and transport. There is a clear need for affordable instrumentation that can provide an early warning of slope instability to enable the evacuation of vulnerable people and timely repair and maintenance of critical infrastructure. An approach, Assessment of Landslides using Acoustic Real-time Monitoring Systems (ALARMS) is described in the paper, and results of a field trial of sensors on an active landslide at Hollin Hill, North Yorkshire, UK, are described. Continuous and real-time monitoring of acoustic emission generated by the deforming slope has been compared to traditional inclinometer slope displacement measurements. Analysis of the results has established that there is a direct relationship between acoustic emission rate and displacement rate trends triggered by rainfall events. The technique has provided insight into reactivated slope movement kinematics.

3D ground model development for an active landslide in Lias mudrocks using geophysical, remote sensing and geotechnical methods

A ground model of an active and complex landslide system in instability prone Lias mudrocks of North Yorkshire, UK is developed through an integrated approach, utilising geophysical, geotechnical and remote sensing investigative methods. Surface geomorphology is mapped and interpreted using immersive 3D visualisation software to interpret airborne light detection and ranging data and aerial photographs. Subsurface structure is determined by core logging and 3D electrical resistivity tomography (ERT), which is deployed at two scales of resolution to provide a means of volumetrically characterising the subsurface expression of both site scale (tens of metres) geological structure, and finer (metre to sub-metre) scale earth-flow related structures. Petrophysical analysis of the borehole core samples is used to develop relationships between the electrical and physical formation properties, to aid calibration and interpretation of 3D ERT images. Results of the landslide investigation reveal that an integrated approach centred on volumetric geophysical imaging successfully achieves a detailed understanding of structure and lithology of a complex landslide system, which cannot be achieved through the use of remotely sensed data or discrete intrusive sampling alone.

Capacitive Resistivity Imaging with Towed Arrays

The capacitive resistivity (CR) technique is a generalization of the conventional DC resistivity method that facilitates measurements of electrical resistivity on engineered surfaces and highly resistive ground. The CR methodology allows the use of towed sensor arrays, thus enabling the rapid collection of high-resolution resistivity data. Under quasi-static conditions, CR data are equivalent to galvanic DC measurements so that CR datasets can be interpreted with conventional DC inversion algorithms. In this study, we demonstrate that the methodology and fundamental parameters of the CR technique facilitate spatial sampling at the centimeter scale for towing speeds of the order of 2.5 km/h. We argue that the dipole-dipole array is the most suitable geometry for dynamic CR measurements and present example data acquired with a prototype instrument using plate-wire combinations arranged in an equatorial geometry. The information content of raw CR data is found to be dominant over towing-induced noise and a direct comparison with DC profile data shows good agreement between both techniques. Based on these findings, we show that tomographic imaging is possible using datasets acquired with moving arrays. Closer investigation of the practical aspects of towed-array capacitive resistivity imaging (CRI) highlights the similarities with galvanic multi-electrode surveys, but the different geometric constraints and sampling regime of CRI give rise to advantages (high lateral resolution) as well as disadvantages (limitations in vertical resolution). We conclude our study with recommendations for practical CRI survey procedures and present a field example where we have successfully imaged a subsurface target in 3D. Towed-array CRI is found to provide equivalent (and in some ways superior) information about the shallow subsurface when compared to DC ERT (electrical resistivity tomography); it therefore widens the scope of electrical imaging surveys to environments where the conventional methodology would be impractical.

Optimization of Array Configurations and Panel Combinations for the Detection and Imaging of Abandoned Mineshafts using 3D Cross-Hole Electrical Resistivity Tomography

Cross-borehole electrical resistivity tomography was used to detect and image a concealed air-filled mineshaft at a greenfield test site. The measurement configurations and panel combinations were selected using a two-stage optimization process. An optimal set of array configurations was selected for each cross-borehole panel on the basis of the model resolution matrix. Subsequently, various combinations of panels were tested with synthetic and field data to determine the effects of coverage and data density on the resulting tomographic image. In the field trials, complicating factors were introduced by the use of resistive cement linings in the boreholes. A resistive feature was detected between the boreholes using a single panel and a 2.5D inversion, but the image quality was too poor to identify this as a mineshaft. A much-improved image was obtained using eight boreholes and eight panels with a full 3D inversion. Only four of these panels intersected the shaft. Crucially, the other panels provided coverage of outlying regions of the model, enabling the inversion algorithm to distinguish between the resistive effects of the borehole linings and the mineshaft.