A new method for geoelectrical investigations underwater1

A new electrical method is proposed for determining the apparent resistivity of multi-earth layers located underwater. The method is based on direct current geoelectric sounding principles. A layered earth model is used to simulate the stratigraphic target. The measurement array is of pole-pole type; it is located underwater and is orientated vertically. This particular electrode configuration is very useful when conventional electrical methods cannot be used, especially if the water depth becomes very important. The calculated apparent resistivity shows a substantial quality increase in the measured signal caused by the underwater targets, from which little or no response is measured using conventional surface electrode methods. In practice, however, different factors such as water stratification, underwater streams or meteorological conditions complicate the interpretation of the field results. A case study is presented, where field surveys carried out on Lake Geneva were interpreted using the calculated apparent resistivity master-curves.     

Integrated geophysical and chemical study of saline water intrusion

Surface geophysical surveys provide an effective way to image the subsurface and the ground water zone without a large number of observation wells. DC resistivity sounding generally identifies the subsurface formations-the aquifer zone as well as the formations saturated with saline/brackish water. However, the method has serious ambiguities in distinguishing the geological formations of similar resistivities such as saline sand and saline clay, or water quality such as fresh or saline, in a low resistivity formation. In order to minimize the ambiguity and ascertain the efficacy of data integration techniques in ground water and saline contamination studies, a combined geophysical survey and periodic chemical analysis of ground water were carried out employing DC resistivity profiling, resistivity sounding, and shallow seismic refraction methods. By constraining resistivity interpretation with inputs from seismic refraction and chemical analysis, the data integration study proved to be a powerful method for identification of the subsurface formations, ground water zones, the subsurface saline/brackish water zones, and the probable mode and cause of saline water intrusion in an inland aquifer. A case study presented here illustrates these principles. Resistivity sounding alone had earlier failed to identify the different formations in the saline environment. Data integration and resistivity interpretation constrained by water quality analysis led to a new concept of minimum resistivity for ground water-bearing zones, which is the optimum value of resistivity of a subsurface formation in an area below which ground water contained in it is saline/brackish and unsuitable for drinking.     

ELECTRICAL SOUNDING OF A HALF-SPACE WITH A MONOTONIC CONTINUOUS VARIATION OF THE RESISTIVITY WITH DEPTH*

The interpretation of electrical sounding data for a subsurface with monotonic continuous variation of the resistivity with depth is becoming increasingly necessary. The contribution of this article is the derivation of the solution for the Wenner and the Schlumberger apparent resistivity functions for a resistivity varying as a real power of a linear positive function of the depth. The interpretation of sounding data in these cases can be used to estimate the variation of the porosity or the salt content of the pore water with depth.     

On Electrical Conductivity at Depth in the Shanuch Copper–Nickel and the Aginskoe Gold-Ore Fields,Kamchatka

This paper is concerned with the method of interpretation and with the results of magnetotelluric sounding (MTS) along the line connecting the Shanuch copper–nickel and the Aginskoe gold ore deposit. The interpretation revealed that the upper section divides into two blocks according to the electrical resistivity. The western block has a high resistivity, which is due to the metamorphic and igneous basement rocks that are exposed at the ground surface. The eastern block has a low resistivity, which is due to the presence of a conductive sedimentary–volcanogenic cover involving a wide network of discontinuities. A low resistivity layer was identified at depths of 30–60 km. The deep-seated faults that have been identified by interpretation in the eastern block penetrate through the lithosphere as far as the conductive layer; fluids can move along these faults upward into the upper crust where at lower temperatures they release minerals whose accumulations make up ore deposits of nonferrous and noble metals.     

Assessment of aquifer vulnerability to industrial waste water using resistivity measurements. A case study, along El-Gharbyia main drain, Nile Delta, Egypt

1D resistivity sounding and 2D resistivity imaging surveys were integrated with geological and hydrochemical data to assess the aquifer vulnerability and saltwater intrusion in the north of Nile Delta, Egypt. In the present study, the El-Gharbyia main drain was considered as a case study to map the sand bodies within the upper silt and clay aquitard. Twenty Schlumberger soundings and six 2D dipole-dipole profiles were executed along one profile close to the western side of the main drain. In addition, 14 groundwater samples and 4 surface water samples from the main drain were chemically analyzed to obtain the major and trace elements concentrations.The results from the resistivity and hydrochemical data were used to assess the protection of the groundwater aquifer and the potential risk of groundwater pollution. The inverted resistivities and thicknesses of the layers above the aquifer layer were used to estimate the integrated electrical conductivity (IEC) that can be used for quantification of aquifer vulnerability. According to the aquifer vulnerability assessment of an underlying sand aquifer, the southern part of the area is characterized by high vulnerability zone with slightly fresh to brackish groundwater and resistivity values of 11-23 Ω.m below the clay cap. The resistivity sections exhibit some sand bodies within the clay cap that lead to increase the recharging of surface waste water (650 mg/l salinity) and flushing the upper part of underlying saltwater aquifer. The region in the north has saltwater with resistivity less than 6 Ω.m and local vulnerable zones within the clay cap. The inverted 2D dipole-dipole profiles in the vulnerable zones, in combination with drilling information have allowed the identification of subsoil structure around the main drain that is highly affected by waste water.     

Pole-dipole resistivity sounding technique for shallow investigations in hard rock areas

Theory of the pole-dipole resistivity sounding technique and its application in the hard rock areas for shallow groundwater exploration is presented in this paper. The different components of electric field produced by the point source of current, situated over the ground surface, are measured by a dipole placed at a large distance from the source. The theory of the method is rather simple, suggesting two configurations, namely radial pole-dipole and axial pole-dipole. Theoretical expressions derived for the apparent resistivity over layered Earth are directly related to the Schlumberger apparent resistivity, whereas expressions for geometrical factor for pole-dipole and Schlumberger configuration are different. The proposed technique has been tested in actual field conditions having different rock types. A few examples are presented along with Schlumberger sounding curves which confirm the applicability of the proposed sounding technique.     

坑道直流电阻率测深异常研究

电阻率测深数据可利用抗道的四个面和四个角得到。但由于受到坑道本身和体积效应影响,坑道内电阻率测深资料解释的困难很大。本文用有限单元法对坑道内的直流电阻率测深进行数值模拟计算。首先给出坑道中正常电位的取值,并利用有限单元异常电位计算法对坑道内电阻率测深进行正演计算。结果表明,在异常体断面尺寸、埋深及所用电极距与坑道截面尺寸相差不大时,坑道面上的测深断面反映的基本是本坑道面外侧地质体的信息,对其他位置的地质体没有明显反映,坑道角断面则反映的是相邻两个坑道面外侧地质体的综合信息。     

Analysis and application of a non-conventional underwater geoelectrical method in Lake Geneva,Switzerland

The electrical method presented is used for determining the resistivity of lake-bottom sediments and is based on the d.c. electrical sounding principles. The electrode array, called the fishing rod (FR), is of pole-pole type and is orientated vertically on a line perpendicular to the surface of the water. The technique is used for mapping resistivity anomalies located deep underwater. This paper presents an analysis of the resolution capabilities of the FR method and the results of a case study carried out in Lake Geneva, where measurements were interpreted using a one-dimensional (1D) multilayer earth model. The analysis of the uncertainty in the model parameters of a 1D multilayer earth model is carried out using the covariance matrix of the linearized inversion problem. The results of the analyses show that when the thickness and resistivity of the water layer is known, the resistivity of the sediment layer is well determined under most circumstances. The thickness of the sediment layer is well determined when resistivity contrasts are not too low. In Lake Geneva the FR method has been used to study an old depression with a resistive channel. This application shows the efficiency of the method compared with conventional electrical methods, where water depth becomes a limiting factor. The use of an automated iterative inversion scheme in this particular case is advantageous, as a joint interpretation of the three different data sets measured with the FR method can be carried out. Finally, the result of the inversion is compared with the trial-and-error interpretations of a previous study.     

水下激发极化测深异常特征研究

利用有限单元法对水下激发极化测深进行数值模拟以研究其异常特征和适用条件.水下的电阻率测深和激电测深均可对水底地质体进行近距离探测,提高探测分辨率.而水底地形对视电阻率测深的影响很大.水底极化体有限元正演模拟的结果表明,水底地形对激发极化测深结果没有影响.采用水下激发极化法进行地质体探测是可行、有效的.通过计算不同水体和水底岩石电阻率对视极化率幅值的影响发现,水体电阻率是影响水下极化率测深的主要因素.由于海水的电阻率过低,观测视极化率异常微弱而不适合开展水下激电工作.     

Theoretical study of slope effects in resistivity surveys and applications

When electric soundings are made over an irregular terrain, topographic effects can influence the values of apparent resistivity and lead to erroneous 1D interpretation. A 3D finite-element method has been applied to study the topographical effect of a slope on Schlumberger soundings parallel to the strike. When the resistivity survey is performed at the top of the slope, the apparent resistivity values can be two times higher than in the flat-earth case, depending on the angle (α) and height (H) of the slope, and on the distance (X) between the sounding and the slope top. The results are presented as nondimensional curves which can be used for evaluating topographic anomalies for any value of the parameters α, H and X. It is numerically shown that the topographic effects can be removed from measurements on horizontally layered structures with an irregular earth surface. Real measurements were performed in different geological conditions over an irregular terrain. The correction method based on the nondimensional curves has been applied to the data and has enabled the determination of the correct layered ground configuration using 1D interpretation.     

Automated 1D interpretation of resistivity soundings by simultaneous use of the direct and iterative methods*

Noise contamination of measured data greatly affects the final results of inversion. Three types of noise source — random and systematic errors and the uncertainties due to the inadequacy of the mathematical model in representing the actual physical conditions — are discussed in the framework of resistivity sounding data. Two methods are proposed for describing these uncertainties. The first possibility is to smooth the measured data by a combination of simple fitting functions that satisfies the ‘1D smoothness’ criteria and consequently simulates the behaviour of a 1D Schlumberger apparent resistivity curve. The second method is to derive weight coefficients from the differences between the measured and the smoothed data sets. Both methods are carried out under the control of the interpreter. The relative merits and drawbacks of the direct and iterative interpretation methods used for the estimation of the parameters of the layered earth model are summarized. Two variants of the combination of these methods are presented to obtain more powerful and automatic interpretation schemes. In the sequential interpretation, an initial guess supplied by the direct method is improved by the iterative method to obtain a reasonable fit between the measured data and the model response. In the simultaneous interpretation, the successive application of the direct and iterative methods is carried out, starting from the first branch of the apparent resistivity curve. The operation is then shifted to subsequent branches that represent the deeper parts of the geoelectric section. This is similar to the data acquisition applied in direct current sounding in which the depth penetration is increased by expanding the current electrode spacings. The proposed sequential and simultaneous interpretation algorithms require minimum aids and efforts of the interpreter.     

ELECTRICAL SOUNDING OF A HALF SPACE WHOSE RESISTIVITY OR ITS INVERSE FUNCTION VARIES LINEARLY WITH DEPTH*

The solution for the potential distribution about a point source of current placed at the surface of a continuous half-space is obtained for two cases: (1) the resistivity increases linearly with depth; (2) the conductivity increases linearly with depth. In each case, an expression for the apparent resistivity is established and master curves are presented for both the Wenner and the Schlumberger configurations. The results can be used in the interpretation of electrical sounding data in specified geologic situations. Furthermore, they may be used as a first step in the development of solutions for the more complex electric sounding problems.     

降雨对昌黎台地电阻率的干扰分析

昌黎台地电阻率受降雨影响明显, 降雨量越大, 地电阻率下降幅度越大, 两者存在良好的线性关系。 基于昌黎台电测深资料反演的水平层状均匀介质结果, 建立地表薄层电性结构变化模型, 通过影响系数理论定性分析降雨对地电阻率的即时效应, 研究结果显示降雨主要是通过改变地表介质含水量和电性结构造成地电阻率变化的。 利用褶积滤波法定量模拟降雨对地电阻率的短期影响过程, 可以较好去除降雨对地电阻率的即时影响和滞后响应。     

MATRIX METHOD FOR THE TRANSFORMATION OF RESISTIVITY SOUNDING DATA OF ONE ELECTRODE CONFIGURATION TO THAT OF ANOTHER CONFIGURATION1

Matrix equations are derived to transform the resistivity sounding data obtained in one type of a four-electrode array to the corresponding resistivity sounding data that would be obtained using a different four-electrode array. These expressions are based primarily on recent work in which we have established a linear relation between the apparent resistivity and the kernel function by using a powerful exponential approximation for the kernel function. It is shown that the resistivity sounding data of two different four-electrode arrays have a linear relation through an essentially non-singular matrix operator and, as such, one is derivable from the other for a one-dimensional model and it can also be extended to two-dimensions. Some numerical examples considering synthetic data are presented which demonstrates the efficiency of the method in such transformations. Two published field examples are also considered for transformation giving a reliable interpretation.     

STUDYING THE INTERCONNECTION OF GROUND AND SURFACE WATERS IN REPRESENTATIVE BASINS UNDER THE PROGRAMME OF THE INTERNATIONAL HYDROLOGICAL DECADE

Abstract

Four geoelectrical soundings were measured with a combination of Schlumberger and azimuthal or equatorial dipole electrode arrays on a Carboniferous limestone basin of the Condroz area, Belgium. The measuring technique is briefly outlined as well as the interpretation procedure, which follows a closed-loop scheme with control of calculated model curves. Some general problems of interpretation of geoelectrical sounding curves are tackled, as far as they have a practical bearing on the treatment of Condroz soundings.

The problem of determining the very high resistivity of limestone is approached through ARCHIE's formula, an empirical relation between the bulk rock resistivity, the porosity and the electrolyte resistivity. An evaluation of the latter two parameters, combined with electrical horizontal conductance measurements directly made on resistivity sounding curves, offers a possibility for fast determination of the total water storage in a limestone aquifer. Such storage determinations could be applied whenever an aquifer shows up as a conductive layer interbedded between two highly resistant layers (e.g. nonsaturated limestone and compact, non-fractured limestone).     

USE OF THE RESISTIVITY METHOD IN GEOLOGICAL MAPPING—CASE HISTORIES FROM RANIGANJ COALFIELD,INDIA*

The effectiveness of the electrical resistivity method has been studied using various configurations with different spacings over the Southern Boundary Fault in the northwestern part, and across a dolerite dyke named “Salma dyke” in the central part, of the Raniganj Coalfield, India. It has been observed that the delineation of the fault and the dyke was made possible under differential tropical weathering conditions. Geoelectric sections across the fault and the dyke have been prepared on the basis of Schlumberger sounding results. In profiling, Wenner, two-electrode, half-Schlumberger and part of Al-Chalabi's configurations were applied over the Southern Boundary Fault. Across the dyke, alpha-, beta-, and gamma-Wenner, Schlumberger, half-Schlumberger and two-electrode configurations were employed. Azi-muthal two-electrode sounding was also studied over the fault. The interpretation of the results of apparent resistivity profiles across the Southern Boundary Fault suggests that the Wenner and two-electrode configurations possess certain diagnostic features which help in mapping a single lithological contact, provided sufficient resistivity contrast exists. Although Schlumberger configuration seems to be quite suitable, other configurations may also be usefully employed over the dyke.     

SIMPLE DATA PROCESSING OF TRIPOTENTIAL APPARENT RESISTIVITY MEASUREMENTS AS AN AID TO THE INTERPRETATION OF SUBSURFACE STRUCTURE*

The advantages of the Wenner tripotential method (Carpenter 1955) for apparent resistivity profiling are described and two new data processing techniques introduced as an aid to the interpretation of apparent resistivity sections (pseudo-sections). These techniques were developed from model data computed using a two-dimensional finite difference method. Oscillatory components present in anomalies on tripotential profiles and related to electrode spacing are shown to be effectively removed by linear filtering that also simplifies their form and aids recognition. Furthermore, it is shown that the ratio of the beta- and gamma-apparent resistivities is a good indicator of resistivity variation, and is particularly sensitive to lateral change. Model data indicates that, over a wide range of conditions, enough subsurface information can be obtained by inspection of tripotential resistivity and ratio profiles, and from space sections to make possible a useful—and sometimes semi-quantitative—interpretation. A rationale for the general interpretation of tripotential data is developed. Field data are described from an area of weathered granite basement in Nigeria. A model of the subsurface is developed using parameters derived from the processed observations. The observed and calculated apparent resistivity space sections are very similar.     

AN INTERPRETATION SYSTEM FOR GEO-ELECTRICAL SOUNDING GRAPHS *

For the two and three layer cases geo-electrical sounding graphs can be rapidly and accurately evaluated by comparing them with an adequate set of standard model graphs. The variety of model graphs required is reasonably limited and the use of a computer is unnecessary for this type of interpretation. For more than three layers a compilation of model graphs is not possible, because the variety of curves required in practice increases immensely. To evaluate a measured graph under these conditions, a model graph is calculated by computer for an approximately calculated resistivity profile which is determined, for example, by means of the auxiliary point methods. This model graph is then compared with the measured curve, and from the deviation between the curves a new resistivity profile is derived, the model graph of which is calculated for another comparison procedure, etc. This type of interpretation, although exact, is very inconvenient and time-consuming, because there is no simple method by which an improved resistivity profile can be derived from the deviations between a model graph and a measured graph. The aim of this paper is, on the one hand, to give a simple interpretation method, suitable for use during field work, for multi-layer geo-electrical sounding graphs, and, on the other hand, to indicate an automatic evaluation procedure based on these principles, suitable for use by digital computer. This interpretation system is based on the resolution of the kernel function of Stefanescu's integral into partial fractions. The system consists of a calculation method for an arbitrary multi-layer case and a highly accurate approximation method for determining those partial fractions which are important for interpretation. The partial fractions are found by fitting three-layer graphs to a measured curve. Using the roots and coefficients of these partial fractions and simple equations derived from the kernel function of Stefanescu's integral, the thicknesses and resistivities of layers may be directly calculated for successively increasing depths. The system also provides a simple method for the approximative construction of model graphs.     

Correlation of vertical electrical sounding and electrical borehole log data for groundwater exploration

We have correlated the longitudinal unit conductance CL obtained from interpreted vertical electrical sounding data with the formation resistivity R_t and the formation resistivity factor F, obtained by carrying out electrical borehole logging. Interpreted geophysical data of eleven soundings and two electrical borehole log records are used for the analysis. The geophysical data used were acquired in a sedimentary basin. The study area is called Lower Maner Basin located in the province of Andhra Pradesh, India. Vertical electrical soundings were carried out using a Schlumberger configuration with half current electrode separation varying from 600–1000 m. For logging the two boreholes, a Widco logger‐model 3200 PLS was used. True formation resistivity R_t was calculated from a resistivity log. Formation resistivity factor F was also calculated at various depths using R_t values. An appreciable inverse relation exists between the correlated parameters. The borehole resistivity R_t and the formation resistivity factor F decrease with the increase in the longitudinal unit conductance CL. We have shown the use of such a relation in computing borehole resistivity R_t and formation resistivity factor F at sites that posses only vertical electrical sounding data, with a fair degree of accuracy. Validation of the correlation is satisfactory. Scope for updating the correlation is discussed. Significance and applications of the relation for exploration of groundwater, namely to update the vertical electrical sounding data interpretation by translating the vertical electrical sounding data into electrical borehole log parameters, to facilitate correlations studies and to estimate the porosity (φ), permeability (K) and water saturation S_w of water bearing zones are discussed.     

Some geoelectrical investigations for potential groundwater zones in part of Azamgarh area of U.P.

Summary The results of some electrical sounding measurements have been used for defining potential groundwater zones in part of Azamgarh district of Uttar Pradesh, India. The area covered is a part of gangetic alluvial tract. Aquifers consist of lenticular sand masses distributed within the alluvium. The ideal geoelectrical conditions, generally desired for the interpretation of resistivity curves, lack in the area. Among the number of layers encountered in lithologs, only a few are delineated by quantitative interpretation of the sounding data. A study of the total transverse resistance and longitudinal conductance for a particular depth of the alluvial column has helped in defining potential groundwater zones.