3-D Magnetotelluric Exploration of Tenerife Geothermal System (Canary Islands,Spain)

The resistivity structure of the Tenerife geothermal system has been determined by the 3-D inversion of data from different magnetotelluric surveys. In this paper, the ocean and topography effects on the magnetotelluric data were investigated by constructing a 3-D conceptual geoelectrical model of the island. The study showed that these effects should be taken into account in order to obtain a reliable subsurface model of the island. Data from 148 sites were used during three-dimensional inversion. The most interesting feature in the final geoelectrical model of the geothermal system is a low resistivity structure (<10 Ωm) above the resistive core of the system. The low resistivity structure has been interpreted as a hydrothermal clay alteration cap typically generated in the conventional geothermal systems. The resistivity model has been correlated with a recent seismic velocity model, showing that a low resistivity structure surrounds an area with high P wave velocity and medium–high resistivity. This medium–high resistivity area can be associated with a slowly solidified magma and, therefore, with a hotter part of the geothermal system.     

Rapid inversion of data from 2D resistivity surveys with electrode displacements

Resistivity monitoring surveys are used to detect temporal changes in the subsurface using repeated measurements over the same site. The positions of the electrodes are typically measured at the start of the survey program and possibly at occasional later times. In areas with unstable ground, such as landslide‐prone slopes, the positions of the electrodes can be displaced by ground movements. If this occurs at times when the positions of the electrodes are not directly measured, they have to be estimated. This can be done by interpolation or, as in recent developments, from the resistivity data using new inverse methods. The smoothness‐constrained least squares optimisation method can be modified to include the electrode positions as additional unknown parameters. The Jacobian matrices with the sensitivity of the apparent resistivity measurements to changes in the electrode positions are then required by the optimisation method. In this paper, a fast adjoint‐equation method is used to calculate the Jacobian matrices required by the least squares method to reduce the calculation time. In areas with large near‐surface resistivity contrasts, the inversion routine sometimes cannot accurately distinguish between electrode displacements and subsurface resistivity variations. To overcome this problem, the model for the initial time‐lapse dataset (with accurately known electrode positions) is used as the starting model for the inversion of the later‐time dataset. This greatly improves the accuracy of the estimated electrode positions compared to the use of a homogeneous half‐space starting model. In areas where the movement of the electrodes is expected to occur in a fixed direction, the method of transformations can be used to include this information as an additional constraint in the optimisation routine.     

The cone penetration test and 2D imaging resistivity as tools to simulate the distribution of hydrocarbons in soil

The purpose of geophysical electrical surveys is to determine the subsurface resistivity distribution by making measurements on the ground surface. From these measurements, the true resistivity of the subsurface can be estimated. The ground resistivity is related to various geological parameters, such as the mineral and fluid content, porosity and degree of water saturation in the rock. Electrical resistivity surveys have been used for many decades in hydrogeological, mining and geotechnical investigations. More recently, they have been used for environmental surveys. To obtain a more accurate subsurface model than is possible with a simple 1-D model, a more complex model must be used. In a 2-D model, the resistivity values are allowed to vary in one horizontal direction (usually referred to as the x direction) but are assumed to be constant in the other horizontal (the y) direction. A more realistic model would be a fully 3-D model where the resistivity values are allowed to change in all three directions. In this research, a simulation of the cone penetration test and 2D imaging resistivity are used as tools to simulate the distribution of hydrocarbons in soil.     

非平坦地形条件下电阻率三维反演

本文实现了非平坦地形条件下电阻率三维反演,讨论了几种消除反演中地形影响的方法. 结果表明,只有将地形直接带入反演算法中,进行带地形电阻率三维反演才能有效消除地形影响及其对反演结果的偏差,得到与地下电性结构相符的反演结果.     

Two-dimensional resistivity imaging and modelling in areas of complex geology

A system is described for the automatic measurement of electrical resistivity pseudo-sections. This comprises a linear array of up to 32 electrodes connected through a multicore cable to a computer controlled switching module and a resistivity meter. The processing of the measured sections to produce two-dimensional true resistivity images of the subsurface is briefly described. Some account is given of the capabilities and limitations of the technique. This is illustrated by a series of computed constant separation traverses for models of simple subsurface structures. Examples of processed images derived from sections measured in areas of relatively complex geology follow, a comparison being made of the interpretations obtained using an automatic imaging method and a manual iterative approach. It is concluded that with the equipment and software so far developed, in areas of modest subsurface geological complexity where some control is available and where the structures are essentially two-dimensional, then good approximations to the true geoelectric sections can be obtained down to depths of between 100 and 200 m.     

The inversion of data from very large three-dimensional electrical resistivity tomography mobile surveys

New developments in mobile resistivity meter instrumentation have made it possible to survey large areas with dense data coverage. The mobile system usually has a limited number of electrodes attached to a cable that is pulled along behind an operator so that a large area can be covered within a short time. Such surveys can produce three-dimensional datasets with hundreds of thousands of electrodes positions and data points. Similarly, the inverse model used to interpret the data can have several hundred thousand cells. It is impractical to model such large datasets within a reasonable time on microcomputers used by many small companies employing standard inversion techniques. We describe a model segmentation technique that subdivides the finite-element mesh used to calculate the apparent resistivity and Jacobian matrix values into a number of smaller meshes. A fast technique that optimizes the calculation of the Jacobian matrix values for multi-channel systems was also developed. A one-dimensional wavelet transform method was then used to compress the storage of the Jacobian matrix, in turn reducing the computer time and memory required to solve the least-squares optimization equation to determine the inverse model resistivity values. The new techniques reduce the calculation time and memory required by more than 80% while producing models that differ by less than 1% from that obtained using the standard inversion technique with a single mesh. We present results using a synthetic model and a field dataset that illustrates the effectiveness of the proposed techniques.     

Practical techniques for 3D resistivity surveys and data inversion1

Techniques to reduce the time needed to carry out 3D resistivity surveys with a moderate number (25 to 100) of electrodes and the computing time required to interpret the data have been developed. The electrodes in a 3D survey are normally arranged in a square grid and the pole-pole array is used to make the potential measurements. The number of measurements required can be reduced to about one-third of the maximum possible number without seriously degrading the resolution of the resulting inversion model by making measurements along the horizontal, vertical and 45° diagonal rows of electrodes passing through the current electrode. The smoothness-constrained least-squares inversion method is used for the data interpretation. The computing time required by this technique can be greatly reduced by using a homogeneous half-space as the starting model so that the Jacobian matrix of partial derivatives can be calculated analytically. A quasi-Newton updating method is then used to estimate the partial derivatives for subsequent iterations. This inversion technique has been tested on synthetic and field data where a satisfactory model is obtained using a modest amount of computer time. On an 80486DX2/66 microcomputer, it takes about 20 minutes to invert the data from a 7 by 7 electrode survey grid. using the techniques described below, 3D resistivity surveys and data inversion can be carried out using commercially available field equipment and an inexpensive microcomputer.     

基于非结构网格的电阻率三维带地形反演

地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探.     

A Coastal Aquifer Study Using Magnetotelluric and Gravity Methods in Abo Zenema, Egypt

Magnetotelluric (MT) soundings and gravity methods were employed to study the deep freshwater aquifer in the area north of Abo Zenema city on the eastern side of the Gulf of Suez, Egypt. Seven MT sites and 48 gravity stations were surveyed along northeast–southwest profiles as close as possible to a line perpendicular to the coast of the Gulf of Suez. The MT survey was conducted using high and low frequencies to investigate shallow and deep areas, respectively. One-dimensional inversion was conducted using a heuristic inversion scheme of the Bostick algorithm. The MT data were also inverted with a 2-D smooth model inversion routine using the nonlinear conjugate gradient method to infer variation in vertical and lateral resistivity inside the Earth. A 100-Ohm-m homogeneous half-space initial model was used to invert the TE mode data only. Then, the inverted model obtained from the TE mode data was used as an initial model for inversion of the TM mode data. The inverted model thus obtained from the TM mode data inversion was used as an initial model for the inversion of the joint TE and TM responses. Two-dimensional (2-D) forward modeling of the gravity data was conducted using the 2-D polygon method of Talwani’s algorithm for an arbitrarily shaped body and was based on the subsurface information from the MT survey and the available information about the geological structure of the study area. This method enabled us to obtain the basement structure of the coastal aquifer in the study area. The results from the analysis and the interpretation of MT and gravity data were used to detect and delineate the groundwater coastal aquifer in the study area.     

Illuminating and optimising a three‐dimensional model of an oil shale seam and its volume distribution using the transient electromagnetic induction method,central part of Jordan

The oil shale exploration program in Jordan is undertaking great activity in the domain of applied geophysical methods to evaluate bitumen‐bearing rock. In the study area, the bituminous marl or oil shale exhibits a rock type dominated by lithofacies layers composed of chalky limestone, marls, clayey marls, and phosphatic marls. The study aims to present enhancements for oil shale seam detection using progressive interpretation from a one‐dimensional inversion to a three‐dimensional modelling and inversion of ground‐based transient electromagnetic data at an area of stressed geological layers. The geophysical survey combined 58 transient electromagnetic sites to produce geoelectrical structures at different depth slices, and cross sections were used to characterise the horizon of the most likely sites for mining oil shale. The results show valuable information on the thickness of the oil shale seam at 3.7 Ωm, which is correlated to the geoelectrical layer between 2‐ and 4 ms transient time delays, and at depths ranging between 85 and 105 m. The 300 m penetrated depth of the transient electromagnetic soundings allows the resolution of the main geological units at narrow resistivity contrast and the distinction of the main geological structures that constrain the detection of the oil shale seam. This geoelectrical layer at different depth slices illustrates a localised oil shale setting and can be spatially correlated with an area bounded by fold and fault systems. Also, three‐dimensional modelling and inversion for synthetic and experimental data are introduced at the faulted area. The results show the limitations of oil shale imaging at a depth exceeding 130 m, which depends on the near‐surface resistivity layer, the low resistivity contrast of the main lithological units, and the degree of geological detail achieved at a suitable model's misfit value.     

Delineation of layer boundaries from smooth models obtained from the geoelectrical sounding data inversion

Two techniques have been presented for the delineation of boundaries from smooth models obtained by smooth inversion techniques of geoelectrical sounding data, such as straightforward inversion scheme, Occam’s and Zohdy’s methods. The smooth model consists of a large number of equally spaced layers, wherein the real geological boundaries are missing. The techniques proposed here suppress the geologically irrelevant boundaries and support the real structural boundaries present in the geoelectrical data. In the first technique, solution of linear inverse problem is improved iteratively through weighted minimum norm inverse, the weight being taken from the current solution. The technique is referred as Iterative Straightforward Inversion Scheme. The second method is analytical, based on the application of smoothing filter, referred in the literature as edge-preserving smoothing. A few examples of theoretical magnetotelluric, dc resistivity and field sounding data have been presented to demonstrate the capabilities of the techniques. The methodologies also reduce the conspicuous oscillations in the smooth solutions caused by the conversion of sharp boundaries to the smooth ones.     

Focused imaging of electrical resistivity data in archaeological prospecting

Electrical resistivity mapping and electrical resistivity profiling are powerful instruments for investigating archaeological structures. Interpretation of geoelectrical data is complicated by near-surface anomalies and the characteristics of the applied electrode arrays. Averaging Wenner α and Wenner β data as an alternative method of focused imaging is presented to overcome these problems. The mechanism of focused imaging is explained using the sensitivity distribution of the combined arrays. Various methods of imaging geoelectrical data are examined with synthetic and field data. In electrical resistivity mapping, inversion of the data is unnecessary when using focused imaging. In electrical resistivity profiling, focused imaging gives a first idea about the subsurface resistivity distribution without achieving the quality obtainable by inversion.     

Pseudo-3D constrained inversion of transient electromagnetic data for a polarizable SMS hydrothermal system in the Deep Sea

Seafloor massive sulfide (SMS) deposits are generated by high-temperature hydrothermal systems. Their precious resources have attracted global interest. A number of investigations with controlled-source electromagnetic (CSEM) methods have been implemented in recent years. There are three major problems with SMS surveying using EM methods. First, SMS imaging techniques for hydrothermal systems have a limited range. Simulations and applications have validated only simple layered models. Second, their inversion efficiencies must be improved further. Laterally constrained inversions and spatially constrained inversions are usually used to map geological structures. However, choosing their suitable weighting parameters is inefficient. Third, the effects of induced polarization (IP) on ore deposits are not considered in such inversions. A non-polarizable model is unable to accurately depict a polarizable model. To resolve these problems, an advanced strategy is used to improve the efficiency of the pseudo-3D inversion process. The proposed imaging method has the ability to map complex 3D geoelectrical structures, and therefore, it can both obtain information regarding surface ore deposits and distinguish between active and inactive hydrothermal systems. However, this method can also be used to depict the distributions of alteration zones and buried deposits. Furthermore, the influences of IP on the inversion are discussed with respect to the Cole- Cole model, and it is shown that the effects of IP on polarizable deposits cannot be ignored during the inversion.     

AIRBORNE RESISTIVITY MAPPING USING A MULTIFREQUENCY ELECTROMAGNETIC SYSTEM*

Airborne electromagnetic methods are most commonly used in mineral exploration. However, new developments, such as multifrequency capability and digital on-board field recording, as well as improvements in instrumentation resulting in high signal-to-noise ratios in recorded data, have made their application in geological mapping possible. A three-frequency airborne EM survey carried out over an area northwest of Timmins, Ontario, was interpreted in terms of thickness and resistivity of the layers of a two-layer earth section. Since both in-phase and quadrature components are measured, this provides six independent parameters at each point in space. Based on prior geological information and a preliminary interpretation of the field records, two two-layer models of the subsurface seemed to be appropriate for most of the survey area. An automatic computerized interpretation procedure was devised to interpret the field data at each point in terms of thickness and resistivity parameters of those two models. When the geology is more complex, the data do not fit the models and no interpretations are made. Two maps illustrating the variation of the resistivity and the thicknesses of the layers were constructed from the interpreted data. These maps agree with the known geological information about the distribution of glacial clay in the area. Areas where the layered models do not fit are known to be areas where the geology is complex with a large number of dykes and other lateral inhomogeneities. The study shows that multifrequency airborne EM surveys can be very useful in geological mapping over inaccessible terrain and can significantly help the mapping geologist where outcrops are scarce.     

Laterally constrained inversion of helicopter-borne frequency-domain electromagnetic data

Helicopter-borne frequency-domain electromagnetic (HEM) surveys are used for fast high-resolution, three-dimensional resistivity mapping. Standard interpretation tools are often based on layered earth inversion procedures which, in general, explain the HEM data sufficiently. As a HEM system is moved while measuring, noise on the data is a common problem. Generally, noisy data will be smoothed prior to inversion using appropriate low-pass filters and consequently information may be lost.For the first time the laterally constrained inversion (LCI) technique has been applied to HEM data combined with the automatic generation of dynamic starting models. The latter is important because it takes the penetration depth of the electromagnetic fields, which can heavily vary in survey areas with different geological settings, into account. The LCI technique, which has been applied to diverse airborne and ground geophysical data sets, has proven to be able to improve the HEM inversion results of layered earth structures. Although single-site 1-D inversion is generally faster and — in case of strong lateral resistivity variations — more flexible, LCI produces resistivity — depth sections which are nearly identical to those derived from noise-free data.The LCI results are compared with standard single-site Marquardt–Levenberg inversion procedures on the basis of synthetic data as well as field data. The model chosen for the generation of synthetic data represents a layered earth structure having an inhomogeneous top layer in order to study the influence of shallow resistivity variations on the resolution of deep horizontal conductors in one-dimensional inversion results. The field data example comprises a wide resistivity range in a sedimentary as well as hard-rock environment.If a sufficient resistivity contrast between air and subsurface exists, the LCI technique is also very useful in correcting for incorrect system altitude measurements by using the altitude as a constrained inversion parameter.     

An audio-magnetotelluric investigation in Terceira Island (Azores)

Ten audio-magnetotelluric soundings have been carried out along a profile crossing the Serra do Cume caldera in the eastern part of the Terceira Island (Azores). The main objectives of this investigation were to detect geoelectrical features related with tectonic structures and to characterize regional hydrological and hydrothermal aspects mainly those related to geothermal fluid dynamics.Three-dimensional numerical investigation showed that the data acquired at periods shorter than 1 s are not significantly affected by ocean effect. The data was analysed using the Smith's decomposition method in order to investigate possible distortions caused by superficial structures and to estimate a global regional strike. The results suggest that in general the soundings were not distorted. A regional N55°W strike was chosen for the two-dimensional data inversion.The low-resistivity zones (10–30 ohm-m) displayed in the central part of the 2-D geoelectrical model have been interpreted as caused by hydrothermal circulation. The low-resistivity anomalies at the ends of the profile might be attributed to alteration zones with interaction of seawater intrusion. High-resistivity (> 300 ohm-m) values have been related with less permeable zones in the SW of Cinco Picos and Guilherme Moniz caldera walls.     

Road evaluation with ground penetrating radar

This paper provides a status report of the Ground Penetrating Radar (GPR) highway applications based on studies conducted in both Scandinavia and the USA. After several years of research local transportation agencies are now beginning to implement GPR technology for both network and project level surveys. This paper summarizes the principles of operation of both ground-coupled and air-launched GPR systems together with a discussion of both signal processing and data interpretation techniques. In the area of subgrade soil evaluation GPR techniques have been used to nondestructively identify soil type, to estimate the thickness of overburden and to evaluate the compressibility and frost susceptibility of subgrade soil. In road structure surveys, GPR has been used to measure layer thickness, to detect subsurface defects and to evaluate base course quality. In quality control surveys, GPR techniques have been used for thickness measurements, to estimate air void content of asphalt surfaces and to detect mix segregation. Future developments are described where the technique has great potential in assisting pavement engineers with their new pavement designs and in determining the optimal repair strategies for deteriorated roadways.     

Exploring seismic inversion methodologies for non‐stationary geological environments: a benchmark study between deterministic and geostatistical seismic inversion

This paper presents a comparison between subsurface impedance models derived from different deterministic and geostatistical seismic inversion methodologies applied to a challenging synthetic dataset. Geostatistical seismic inversion methodologies nowadays are common place in both industry and academia, contrasting with traditional deterministic seismic inversion methodologies that are becoming less used as part of the geo‐modelling workflow. While the first set of techniques allows the simultaneous inference of the best‐fit inverse model along with the spatial uncertainty of the subsurface elastic property of interest, the second family of inverse methodology has proven results in correctly predicting the subsurface elastic properties of interest with comparatively less computational cost. We present herein the results of a benchmark study performed over a realistic three‐dimensional non‐stationary synthetic dataset in order to assess the performance and convergence of different deterministic and geostatistical seismic inverse methodologies. We also compare and discuss the impact of the inversion parameterisation over the exploration of the model parameter space. The results show that the chosen seismic inversion methodology should always be dependent on the type and quantity of the available data, both seismic and well‐log, and the complexity of the geological environment versus the assumptions behind each inversion technique. The assessment of the model parameter space shows that the initial guess of traditional deterministic seismic inversion methodologies is of high importance since it will determine the location of the best‐fit inverse solution.     

Non-invasive Geophysical Surveys in Search of the Roman Temple of Augustus Under the Cathedral of Tarragona (Catalonia,Spain): A Case Study

An integrated geophysical survey has been conducted at the Tarragona’s Cathedral (Catalonia, NE Spain) with the aim to confirm the potential occurrence of archaeological remains of the Roman Temple dedicated to the Emperor Augustus. Many hypotheses have been proposed about its possible location, the last ones regarding the inner part of the Cathedral, which is one of the most renowned temples of Spain (twelfth century) evolving from Romanesque to Gothic styles. A geophysical project including electrical resistivity tomography (ERT) and ground probing radar (GPR) was planned over 1 year considering the administrative and logistic difficulties of such a project inside a cathedral of religious veneration. Finally, both ERT and GPR have been conducted during a week of intensive overnight surveys that provided detailed information on subsurface existing structures. The ERT method has been applied using different techniques and arrays, ranging from standard Wenner–Schlumberger 2D sections to full 3D electrical imaging with the advanced Maximum Yield Grid array. Electrical resistivity data were recorded extensively, making available many thousands of apparent resistivity data to obtain a complete 3D image after a full inversion. In conclusion, some significant buried structures have been revealed providing conclusive information for archaeologists. GPR results provided additional information about shallowest structures. The geophysical results were clear enough to persuade religious authorities and archaeologists to conduct selected excavations in the most promising areas that confirmed the interpretation of geophysical data. In conclusion, the significant buried structures revealed by geophysical methods under the cathedral were confirmed by archaeological digging as the basement of the impressive Roman Temple that headed the Provincial Forum of Tarraco, seat of the Concilium of Hispania Citerior Province.     

Delineating coastal groundwater discharge processes in a wetland area by means of electrical resistivity imaging, 224Ra and 222Rn

Coastal groundwater discharge (CGD) plays an important role in coastal hydrogeological systems as they are a water resource that needs to be managed, particularly in wetland areas. Despite its importance, identifying and monitoring CGD often presents physical and logistical constraints, restraining the application of more traditional submarine groundwater discharge surveying techniques. Here we investigate the capability of electrical resistivity imaging (ERI) in the Peníscola wetland (Mediterranean coast, Spain). ERI surveying made it possible to identify and delineate an ascending regional groundwater flow of thermal and Ra‐enriched groundwater converging with local flows and seawater intrusion. The continuous inputs of Ra‐rich groundwater have induced high activities of Ra isotopes and ²²²Rn into the marsh area, becoming among the highest previously reported in wetlands and coastal lagoons. Geoelectrical imaging enabled inferring focused upward discharging areas, leaking from the aquifer roof through a confining unit and culminating as spring pools nourishing the wetland system. Forward modelling over idealized subsurface configurations, borehole datasets, potentiometric records from standpipe piezometers, petrophysical analysis, and four natural and independent tracers (²²⁴Ra, ²²²Rn, temperature and salinity) permitted assessing the geoelectrical model and a derived hydrogeological pattern. The research highlights the potential of ERI to improve hydrogeological characterization of subsurface processes in complex contexts, with different converging flows. Additionally, a hydrogeological conceptual model for a groundwater‐fed coastal wetland was proposed, based on the integration of surveying datasets. Copyright © 2013 John Wiley & Sons, Ltd.