Using 2D inversion of magnetic resonance soundings to locate a water-filled karst conduit

A new methodology for magnetic resonance sounding (MRS) data acquisition and interpretation was developed for locating water-filled karst cavities. This methodology was used to investigate the Ouysse karst system in the Poumeyssens shaft in the Causse de Gramat (France). A new 2D numerical MRS response model was designed for improved accuracy over the previous 1D MRS approach. A special survey performed by cave divers confirmed the accuracy of the MRS results. Field results demonstrated that in favourable conditions (a low EM noise environment and a relatively shallow, large target) the MRS method, used with a coincident transmitter/receiver loop, can be an effective tool for locating a water-filled karst conduit. It was shown numerically that because an a priori orientation of the MRS profile with the karst conduit is used in the inversion scheme (perpendicular for instance), any error in this assumption introduces an additional error in locating the karst. However, the resulting error is within acceptable limits when the deviation is less than 30°. The MRS results were compared with an electrical resistivity tomography (ERT) survey. It was found that in Poumeyssens, ERT is not able to locate the water-filled karst. On the other hand, ERT provides additional information about heterogeneities in the limestone.     

EQUIVALENCE IN ELECTROMAGNETIC (FREQUENCY) SOUNDING*

An analysis of the equivalence phenomenon, analogous to that encountered in resistivity sounding, has been made for electromagnetic (frequency) sounding with the following systems: horizontal coplanar coils, perpendicular coils, vertical coplanar coils, and vertical coaxial coils. Cases analyzed are three-layer H and K type resistivity distributions. The theoretical responses for the EM sounding systems have been computed by the digital linear filter method using short filters. An analysis has been made of the resolution by the four EM systems considered f the equivalence effect. It is concluded that from the equivalence point of view electromagnetic (frequency) sounding has relatively better resolution compared to resistivity sounding for the H type three-layer cases. For the K type three-layer cases electromagnetic (frequency) sounding provides a serious problem in resolution similar to that encountered in resistivity sounding. Empirical relationships have been established for the H and K type cases respectively in order to numerically quantify the equivalence involved in electromagnetic (frequency) sounding.     

Challenges of using electrical resistivity method to locate karst conduits—A field case in the Inner Bluegrass Region, Kentucky

Conduits serve as major pathways for groundwater flow in karst aquifers. Locating them from the surface, however, is one of the most challenging tasks in karst research. Geophysical methods are often deployed to help locate voids by mapping variations of physical properties of the subsurface. Conduits can cause significant contrasts of some physical properties that can be detected; other subsurface features such as water-bearing fractures often yield similar contrasts, which are difficult to distinguish from the effects of the conduits. This study used electrical resistivity method to search for an unmapped karst conduit that recharges Royal Spring in the Inner Bluegrass karst region, Kentucky, USA. Three types of resistivity techniques (surface 2D survey, quasi-3D survey, and time-lapse survey) were used to map and characterize resistivity anomalies. Some of the major anomalies were selected as drilling targets to verify the existence of the conduits. Drilling near an anomaly identified by an electrical resistivity profile resulted in successful penetration of a major water-filled conduit. The drilling results also suggest that, in this study area, low resistivity anomalies in general are associated with water-bearing features. However, differences in the anomaly signals between the water-filled conduit and other water-bearing features such as water-filled fracture zones were undistinguishable. The electrical resistivity method is useful in conduit detection by providing potential drilling targets. Knowledge of geology and hydrogeology about the site and professional judgment also played important roles in locating the major conduit.     

The Application of Geotomography in Engineering

It seems that electromagnetic tomography is quite suitable to be used in metal and groundwaterexploration. especial for limestone areas where the absorption coefficient of electromagnetic wave is smalland geological body has a large property contrast with country rock it should be pointed out that there exitsome uncertainties in geological interpretation of EM tomogram, corresponding laboratory tests for rockphysical property are always necessary.Resistivity tomography provides us a new useful…     

Extending electromagnetic methods to map coastal pore water salinities

The feasibility of mapping pore water salinity based on surface electromagnetic (EM) methods over land and shallow marine water is examined in a coastal wetland on Tampa Bay, Florida. Forward models predict that useful information on seabed conductivity can be obtained through <1.5 m of saline water, using floating EM-31 and EM-34 instruments from Geonics Ltd. The EM-31 functioned as predicted when compared against resistivity soundings and pore water samples and proved valuable for profiling in otherwise inaccessible terrain due to its relatively small size. Experiments with the EM-34 in marine water, however, did not reproduce the theoretical instrument response. The most effective technique for predicting pore water conductivities based on EM data entailed (1) computing formation factors from resistivity surveys and pore water samples at representative sites and (2) combining these formation factors with onshore and offshore EM-31 readings for broader spatial coverage. This method proved successful for imaging zones of elevated pore water conductivities/salinities associated with mangrove forests, presumably caused by salt water exclusion by mangrove roots. These zones extend 5 to 10 m seaward from mangrove trunks fringing Tampa Bay. Modeling indicates that EM-31 measurements lack the resolution necessary to image the subtle pore water conductivity variations expected in association with diffuse submarine ground water discharge of fresher water in the marine water of Tampa Bay. The technique has potential for locating high-contrast zones and other pore water salinity anomalies in areas not accessible to conventional marine- or land-based resistivity arrays and hence may be useful for studies of coastal-wetland ecosystems.     

Delineation of a quick clay zone at Smørgrav,Norway, with electromagnetic methods under geotechnical constraints

In many coastal areas of North America and Scandinavia, post-glacial clay sediments have emerged above sea level due to iso-static uplift. These clays are often destabilised by fresh water leaching and transformed to so-called quick clays as at the investigated area at Smørgrav, Norway. Slight mechanical disturbances of these materials may trigger landslides. Since the leaching increases the electrical resistivity of quick clay as compared to normal marine clay, the application of electromagnetic (EM) methods is of particular interest in the study of quick clay structures.For the first time, single and joint inversions of direct-current resistivity (DCR), radiomagnetotelluric (RMT) and controlled-source audiomagnetotelluric (CSAMT) data were applied to delineate a zone of quick clay. The resulting 2-D models of electrical resistivity correlate excellently with previously published data from a ground conductivity metre and resistivity logs from two resistivity cone penetration tests (RCPT) into marine clay and quick clay. The RCPT log into the central part of the quick clay identifies the electrical resistivity of the quick clay structure to lie between 10 and 80 Ω m. In combination with the 2-D inversion models, it becomes possible to delineate the vertical and horizontal extent of the quick clay zone. As compared to the inversions of single data sets, the joint inversion model exhibits sharper resistivity contrasts and its resistivity values are more characteristic of the expected geology. In our preferred joint inversion model, there is a clear demarcation between dry soil, marine clay, quick clay and bedrock, which consists of alum shale and limestone.     

Detection of alteration at the Millennium uranium deposit in the Athabasca Basin: a comparison of data from two airborne electromagnetic systems with ground resistivity data

The Millennium uranium deposit is located within the Athabasca Basin of northern Saskatchewan. The basement rocks, comprised primarily of paleo‐Proterozoic gneisses, are electrically resistive. However, the deposit is associated with highly conductive graphitic metasediments that are intercalated with the gneisses. An unconformity separates the basement rocks from the overlying, horizontally stratified, Proterozoic sandstones of the Athabasca Group (which are also highly resistive). The strike extents of the graphitic metasedimenary packages are extensive and therefore electromagnetic (EM) survey techniques are successful at identifying these zones but do not identify the specific locations where they are enriched in uranium. Through drilling it has been noted that hydrothermal processes associated with mineralization has altered the rocks in the vicinity of the deposits, which should in theory result in a resistivity low. A significant resistivity low has been mapped coincident with the Millennium deposit using ground resistivity survey techniques. However, a comparison of the airborne EM and ground resistivity results reveals that the two data sets have imaged different features. The resistive‐limit (on‐time) windows of the MEGATEM data show conductive features corresponding to lakes located to the west and south of the deposit. The late‐time windows show a feature to the east of the deposit, interpreted as being associated with the east‐dipping graphitic basement conductors (similar to that observed in historical ground EM data collected in this area). The early‐time TEMPEST windows (delay times less than 0.2 ms) show a broad resistivity low located at approximately the same location as where the alteration has been identified through drilling. Modelling the data is not easy but a response that decays prior to 0.3 ms is consistent with 500 Ωm material in the sandstone, a resistivity value close to the lower limit with respect to the hydrothermally altered Athabasca group sediments in this area. The MEGATEM system does not see a conductive zone over the alteration as clearly but the high signal‐to‐noise ratio in the late‐time MEGATEM data means that the conductive material at a greater depth is more coherently imaged.     

The competing impacts of geology and groundwater on electrical resistivity around Wrightsville Beach,NC

Measurements of electrical resistivity have long been used to find freshwater resources below the earth's surface. Recently, offshore resistivity and electromagnetic techniques have been used to map occurrences of submarine groundwater originating from the offshore extension of terrestrial aquifers. In many cases, observations of a high resistivity (low conductivity) anomaly in the seafloor are sufficient to suggest the presence of fresh (and less conductive) pore waters. Data from offshore Wrightsville Beach, NC show highly variable resistivity structure, with moderately high resistivity at depths of ～20 m subsurface that is at least in part caused by lithologic complexity in an underlying limestone unit, the Castle Hayne. These offshore results suggest caution in the interpretation of resistivity anomalies simply in terms of groundwater volumes. In contrast, low onshore resistivities show evidence for intrusion of saltwater into the subsurface beneath the beach, adjacent to areas of pumping for water supply.     

ELECTROMAGNETIC PROSPECTING FOR GROUNDWATER IN PRECAMBRIAN TERRAINS IN THE REPUBLIC OF UPPER VOLTA*

The “Autorité des amenegements des valées des Voltas (AVV)” is establishing new rural settlements in the Volta valleys. First, a survey of available water supplies is performed. Economic aquifers in Precambrian terrains are deep (15–50 m) and usually occur in fractured zones accompanying faults. Such zones can be identified on aerial photographs, but their precise location on the ground is virtually impossible by visual means. Because of the small size of the aquifers, a location error of 5 m can make the difference between a productive well and a dry hole. Traditionally, resistivity profiling has been used as the means of locating the fractured zones in the field. Our studies suggest that the task can be performed faster, cheaper and more accurately by VLF and EM methods. Because of the limited choice of transmitting stations reccivable in Upper Volta, the VLF method is not sufficiently sensitive to detect conductors with a strike between 45° and 105°. The results obtained with a multifrequency, horizontal-loop EM (HLEM) system were satisfactory in all investigated areas. During the 1980 field season, 35 target areas were surveyed. Of the 24 holes drilled so far, 23 are productive. The weathered layer is a source of distinctive HLEM anomalies, which are characteristic of the underlying rocks. Therefore, different interpretational procedures had to be developed for granitic and volcano-sedimentary areas. Despite the high background level of in-phase and quadrature components, which varied with thickness and conductivity of the weathered layer, aquifers could be detected at a depth greater than 30 m. Attempts were made to interpret the HLEM results quantitatively using two models: a three-layer medium and a valley discontinuity. The latter model is more realistic, but more scale modelling will have to be performed to permit development of viable interpretational procedures. Meanwhile, phasor diagrams based on drilling and resistivity sounding data can be used to estimate the aquifer depth.     

The applicability of earth resistivity methods for saline interface definition

Direct current resistivity traversing has been used to characterise the nature of the saline interface at Te Horo on the Kapiti Coast in New Zealand. The results show that the interface in the vicinity of the settlement, which relies on bores for potable water, has intruded inland 10 m further than in undeveloped areas. Resistivity traversing has been particularly successful in defining subsurface areas of higher salinity by providing a two-dimensional image of the bulk resistivity structure. The results of the resistivity surveys are supported by bore water chemistry, which show evidence of saltwater mixing. Bores on beachfront properties have concentrations of up to 1% seawater. A resistivity formation factor has been derived to allow the pore fluid resistivity to be estimated for future coastal surveys. The results also illustrate the problems associated with using standard models to predict the location of the saline interface when small amounts of diffusive mixing occur.     

Electromagnetic mapping of electrical conductivity beneath the Columbia basalts1

Sedimentary rocks beneath the Columbia River Basalt Group are recognized as having potential for oil and gas production, but the overlying layered basalts effectively mask seismic reflections from the underlying sediments. Four electromagnetic (EM) methods have been applied on profiles crossing Boylston Ridge, a typical east–west trending anticline of the Yakima Fold Belt, in an attempt to map the resistivity interface between the basalts and the sediments and to map variations in structure and resistivity within the sediments. The EM surveys detected strong variations in resistivity within the basalts, and in particular the continuous magnetotelluric array profiling (EMAP) revealed resistivity lows beneath the surface anticlines. These low resistivity zones probably coincide with fracturing in the core of the anticlines and they appear to correlate well with similar zones of low seismic velocity observed on a nearby seismic profile. The controlled-source EM surveys (in-loop transient, long-offset transient, and variable-offset frequency-domain) were designed in anticipation of relatively uniform high resistivity basalts, and were found to have been seriously distorted by the intrabasalt conductors discovered in the field. In particular, the resistivity sections derived from 1D inversions were found to be inconsistent and misleading. The EMAP survey provided the most information about the subsurface resistivity distribution, and was certainly the most cost-effective. However, both controlled-source and EMAP surveys call for accurate 2D or 3D inversion to accommodate the geological objectives of this project.     

APPLICATION OF THE MULTIFREQUENCY HORIZONTAL-LOOP EM METHOD IN OVERBURDEN INVESTIGATIONS1

In recent years, geophysical methods (shallow seismic, electromagnetic, resistivity, ground penetrating radar) have been increasingly applied to overburden investigations. Their effectiveness has been found to depend significantly on local geological conditions. Compared with advanced seismic techniques, EM methods are faster and hence more cost-effective, but they have not been considered sufficiently accurate. Analysis is carried out of data obtained with the multifrequency horizontal-loop method (HLEM) in northeastern Ontario, where the overburden consists of Quaternary glacial and glaciolacustrine sediments. Surveying along 1-6 km long profiles permitted recognition of bedrock inhomogeneities and selection of sites suitable for HLEM data interpretation using the layered model. Phasor diagrams and computer inversion based on the ridge regression technique were used to interpret HLEM soundings obtained at eight frequencies. Interpreted layer resistivities and thicknesses were correlated with the results of Rotasonic drilling at 70 sites. Relatively accurate estimates of overburden thickness (within 10%) could be obtained in about 80% of the cases. Nine examples of HLEM soundings are given and discussed: three each of one-, two- and three-layer situations. An appropriate interpretation model cannot be selected simply by minimizing the rms error or by analysing the parameter resolution matrix. Frequently, the most effective way of evaluating a solution is to consider whether resistivity values determined by inversion fit any of the ranges determined by statistical analyses of sediment resistivities. A previously published study of electrical properties of Quaternary sediments indicated that resistivities of clay, till and sand are stable within a fairly large area, such as the one under investigation. While the application of HLEM methods to mapping of Quaternary sediments can be considered a success, interpretation of EM data in regions covered by glacial sediments is more difficult than in weathered terrains, where near-surface layering is more predictable. The problem of equivalence causes non-uniqueness in interpretation. Thickness equivalence, which results in poor resistivity estimates, was found to affect areas convered by sand and till. Conductance equivalence caused poor resolution of thickness and resistivity for thin clay layers (less than 10 m).     

DETECTION OF SOLUTION CHANNELS IN LIMESTONE BY ELECTRICAL RESISTIVITY METHOD*

The results of geophysical surveys carried out in the Kopili Hydel Project site in the United North Cachar and Mikir Hills district, Assam, to locate and trace underground caverns in limestone are presented. Electrical resistivity observations over known limestone caverns showed significant anomalies (high resistivity). In the course of the systematic electrical resistivity survey conducted at the eastern fringe of the Umrong basin to study the competency of the proposed reservoir, a number of indications of possible caverns in limestone were obtained and a few of them have been verified by drilling. An analysis of the electrical resistivity data also suggests that the deep seated caverns are aligned along an axis trending NNW-SSE, which may represent a zone of weakness– presumably a fault–that has influenced the formation of the caverns in preferred direction. Drilling over a few of the geophysical indications have corroborated the presence of a fault along the eastern fringe of the Umrong basin.     

航空电磁勘查技术发展现状及展望

航空电磁作为一种高效的地球物理勘查技术手段,其发展在国外(加拿大、澳大利亚等国家)已十分成熟.然而,在我国该项技术仍处于发展当中,在国内目前尚未形成具有实际探测能力的航空电磁系统和解释手段.这一现状严重影响了我国对地形地质条件复杂区域(比如广大西部地区)矿产资源勘查的需求.本文旨在通过系统介绍航空电磁勘查技术中的基础理论、关键技术、仪器系统、数据处理、解释及应用,并对未来我国航空电磁勘查技术的发展提出建议,使读者了解该技术未来发展方向和研究热点,以期该项技术在我国得到快速发展并获得广泛应用.     

Pitfalls of 1D inversion of small-loop electromagnetic data for detecting man-made objects

In order to interpret field data from small-loop electromagnetic (EM) instruments with fixed source–receiver separation, 1D inversion method is commonly used due to its efficiency with regard to computation costs. This application of 1D inversion is based on the assumption that small-offset broadband EM signals are insensitive to lateral resistivity variation. However, this assumption can be false when isolated conductive bodies such as man-made objects are embedded in the earth. Thus, we need to clarify the applicability of the 1D inversion method for small-loop EM data. In order to systematically analyze this conventional inversion approach, we developed a 2D EM inversion algorithm and verified this algorithm with a synthetic EM data set. 1D and 2D inversions were applied to synthetic and field EM data sets. The comparison of these inversion results shows that the resistivity distribution of the subsurface constructed by the 1D inversion approach can be distorted when the earth contains man-made objects, because they induce drastic variation of the resistivity distribution. By analyzing the integrated sensitivity of the small-loop EM method, we found that this pitfall of 1D inversion may be caused by the considerable sensitivity of the small-loop EM responses to lateral resistivity variation. However, the application of our 2D inversion algorithm to synthetic and field EM data sets demonstrate that the pitfall of 1D inversion due to man-made objects can be successfully alleviated. Thus, 2D EM inversion is strongly recommended for detecting conductive isolated bodies, such as man-made objects, whereas this approach may not always be essential for interpreting the EM field data.     

利用核磁共振方法探查岩溶水

核磁共振（ＮＭＲ）找水方法是目前唯一直接找水的物探新方法,该方法可以探查各种类型的地下水。本文仅论述了在岩溶发育区常规电阻率方法找水时遇到的地质干扰。ＮＭＲ找水方法的原理,仪器的组成,并以实例说明了ＮＭＲ找水方法探查岩溶水的效果。     

电偶源频率电磁测深二维阻抗视电阻率计算的源效应校正法

研讨了频率域电磁法中不同源装置的大地电磁测深、线源频率电磁测深和偶极源频率电磁测深阻抗视电阻率的源效应影响特征。在唯象分析的基础上，提出了几种电磁测深法阻抗视电阻率的相互换算法──源效应校正法（大地电磁测深二维ＴＥ极化视电阻率和其它两种电磁法的赤道装置二维阻抗视电阻车）。模型试验表明，利用这一源效应校正法可以由大地电磁二维视电阻率近似地计算出线源频率电磁二维阻抗视电阻率。这一方法被尝试应用于由线源频率电磁二维阻抗视电阻率估算偶极源频率电磁二维阻抗视电阻率。     

Time‐lapse electrical resistivity imaging of solute transport in a karst conduit

The use of electrical resistivity surveys to locate karst conduits has shown mixed success. However, time‐lapse electrical resistivity imaging combined with salt injection improves conduit detection and can yield valuable insight into solute transport behaviour. We present a proof of concept above a known karst conduit in the Kentucky Horse Park (Lexington, Kentucky). A salt tracer solution was injected into a karst window over a 45‐min interval, and repeat resistivity surveys were collected every 20 min along a 125‐m transect near a monitoring well approximately 750 m downgradient from the injection site. In situ fluid conductivity measurements in the well peaked at approximately 25% of the initial value about 3 h after salt injection. Time‐lapse electrical resistivity inversions show two broad zones at the approximate conduit depth where resistivity decreased and then recovered in general agreement with in situ measurements. Combined salt injection and electrical resistivity imaging are a promising tool for locating karst conduits. The method is also useful for gaining insight into conduit geometry and could be expanded to include multiple electrical resistivity transects. Copyright © 2015 John Wiley & Sons, Ltd.     

Forward modeling of marine DC resistivity method for a layered anisotropic earth

Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell’s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic (EM) fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.     

A Method to Remove Electromagnetic Coupling from Induced Polarization Data for an "Exponential" Earth Model

--The electromagnetic (EM) coupling effect in induced polarization (IP) data is an important problem. In many works it has been computed only considering homogeneous or layered earth models with discretely uniform conductivity. In this study, an algorithm has been developed to compute the EM coupling effect in IP data measured on the earth, whose conductivity varies (increases or decreases) exponentially with depth. The EM coupling effects for Percent Frequency Effect (PFE) and phase data are computed for a dipole-dipole array with different separations, however the method can be applied to any electrode array. The results obtained for the cases of increasing and decreasing conductivity as a function of depth indicate that the EM coupling effect strongly depends on the subsurface resistivity and the dipole length. Here an "exponential" earth model is considered to remove EM coupling from the IP data in frequency and phase domain. For this purpose, first, the region of pseudo-section is divided into segments, and within each segment a typical average apparent resistivity („_a) curve is constructed. An exponential conductivity model is fitted to average „_a data. The conductivity model is then used to compute EM responses. Next the data are corrected for the EM coupling contribution. This decoupling process is applied to field data from a galenite-pyrite mineralization area at the Dolluk site, in western Turkey. The results from the decoupling method developed here are compared with other techniques.