A SIMPLE METHOD OF INTERPRETATION OF RESISTIVITY SOUNDING DATA USING EXPONENTIAL APPROXIMATION OF THE KERNEL FUNCTION*

A simple unified equation of apparent resistivity for a general four-electrode array is developed. The main idea is the analytical integration of the Stefanescu expression for potential over a layered earth by writing an exponential approximation for the kernel function. Finally a matrix equation is developed to estimate the kernel function from observed apparent resistivity values. The general equation automatically reduces to the particular configuration once the electrode separations are modified suitably. Examples for Schlumberger and Wenner configurations are numerically calculated to estimate the precision of the method. Good results in a short execution time are obtained, irrespective of the shape of the apparent resistivity curve. Finally, the full interpretation of one theoretical resistivity curve and two field resistivity curves is demonstrated. The more stable ridge-regression estimation method is used in the identification of layer parameters from the kernel function.     

DIRECT INTERPRETATION OF TWO-ELECTRODE RESISTIVITY SOUNDINGS *

This paper describes the procedure for interpreting the apparent resistivity data measured with the two-electrode array directly with the help of kernel function. The calculation of kernel function from the observed resistivity curve is done by the method of decomposition. In the method of decomposition the resistivity curve is approximated by a sum of certain functions, whose choice is only restricted by the requirement that the contribution to the kernel function corresponding to them should be easily computable. A few such functions are classified. These, and the standard curves for corresponding kernel functions obtained by utilising an integral expression for two-electrode array expressing the kernel explicitly in terms of the apparent resistivity functions, are plotted on log-log scale. The determination of layer parameters, that is, the layer resistivities and thicknesses from the kernel function can be carried out by a method proposed by Pekeris (1940).     

A STUDY ON THE DIRECT INTERPRETATION OF RESISTIVITY SOUNDING DATA MEASURED BY WENNER ELECTRODE CONFIGURATION*

This paper describes certain procedures for deriving from the apparent resistivity data as measured by the Wenner electrode configuration two functions, known as the kernel and the associated kernel respectively, both of which are functions dependent on the layer resistivities and thicknesses. It is shown that the solution of the integral equation for the Wenner electrode configuration leads directly to the associated kernel, from which an integral expression expressing the kernel explicitly in terms of the apparent resistivity function can be derived. The kernel is related to the associated kernel by a simple functional equation where K₁(λ) is the kernel and B₁(λ) the associated kernel. Composite numerical quadrature formulas and also integration formulas based on partial approximation of the integrand by a parabolic arc within a small interval are developed for the calculation of the kernel and the associated kernel from apparent resistivity data. Both techniques of integration require knowledge of the values of the apparent resistivity function at points lying between the input data points. It is shown that such unknown values of the apparent resistivity function can satisfactorily be obtained by interpolation using the least-squares method. The least-squares method involves the approximation of the observed set of apparent resistivity data by orthogonal polynomials generated by Forsythe's method (Forsythe 1956). Values of the kernel and of the associated kernel obtained by numerical integration compare favourably with the corresponding theoretical values of these functions.     

INVERSE FILTER COEFFICIENTS FOR THE COMPUTATION OF APPARENT RESISTIVITY STANDARD CURVES FOR A HORIZONTALLY STRATIFIED EARTH*

In this paper a fast method is developed for computing apparent resistivity curves for known layer configurations. The method is based on the application of a linear filter to determine the apparent resistivity curve from, the kernel function.     

A FAST METHOD FOR DETERMINING THE LAYER DISTRIBUTION FROM THE RAISED KERNEL FUNCTION IN GEOELEGTRICAL SOUNDING*

In a previous publication (Koefoed 1968) a function called the “raised kernel function” has been introduced as an intermediate function in the interpretation of resistivity sounding data, and methods have been described both for the determination of the raised kernel function from the apparent resistivity function, and for the determination of the layer distribution from the raised kernel function. In the present paper a procedure is described by which the second step in this interpretation method–i.e. the determination of the layer distribution from the raised kernel function–is considerably accelerated. This gain in interpretation speed is attained by the use of a standard graph for a function which defines the reduction of the raised kernel function to a lower boundary plane.     

AN ANALYSIS OF EQUIVALENCE IN RESISTIVITY SOUNDING*

A mathematical analysis is given of the phenomenon of equivalence in resistivity sounding, which is based upon the properties of the raised kernel function. Analysis of this function instead of the apparent resistivity function is justified because, as has been shown in a previous publication, variations in the apparent resistivity function lead to variations in the raised kernel function with relative values of the same order of magnitude The expression for the raised kernel function is expanded into a Mac Laurin series. Equivalence can occur only if the second order term of this series is negligible. The coefficient of the first order term depends on the resistivity and the thickness of the layer under consideration. There is an infinite set of combinations of values for these two quantities, for which the coefficient of the first order term has the same value. All these combinations represent equivalent layer distributions.     

THE DIRECT INTERPRETATION OF RESISTIVITY OBSERVATIONS MADE WITH A WENNER ELECTRODE CONFIGURATION*

In a previous paper by the present author a method was developed for direct interpretation of resistivity observations made with a Schlumberger electrode configuration. This method consisted of two steps. The first of these was to derive the kernel function in the integral expression for the apparent resistivity from the observed data; the second step was to derive the resistivity stratification from this kernel function. The first of these two steps depends on the electrode configuration that has been used. In the present paper the above mentioned method is modified so as to make it apply to a Wenner electrode configuration. The procedure is indicated by which the method may be adapted to any other electrode configuration in which the distances between the electrodes are finite. The second step in the interpretation, i.e. the derivation of the resistivity stratification from the kernel function, is independent of the electrode configuration used, and therefore needs no further discussion in the present context.     

On the automatic interpretation of direct current resistivity soundings

A method for the automatic inversion of resistivity soundings is presented. The procedure consists of two main stages. First, application of linear filters which transforms the apparent resistivity curve into the kernel function, and vice versa. In the second stage the first and second derivatives of the kernel function are calculated and used in a second-order modified Newton-Raphson iterative fitting procedure. The model obtained is optimal in the least squares sense. The method has been tried on some field examples and produced resistivity models which show a good agreement with the geological well logs.     

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.     

LES COURBES DAR ZARROUK ET LES SONDAGES DIFFÉRENTIELS DE RÉSISTIVÉ*

Differential and difference soundings of resistivity for real, imaginary and complex values of the newly introduced parameters ν and μ are considered. The transverse and longitudinal differential soundings, known up to the present, correspond with the values ν= 1 and ν=— 1. The point of departure is the generalization of the concept of Dar Zarrouk parameters and the dimensional analogy between D.Z. resistivity and transverse and longitudinal differential resistivities. Properties of generalized D.Z. curves as well as the method of their construction are given. The examples of differential and difference curves are presented and the properties of these curves are discussed; particulary the effective spacings are determined. The problems of the realization of differential and difference soundings are considered; the bisymmetrical and asymmetrical arrays are proposed as a practical ones.     

THE USE OF FILTERED BESSEL FUNCTIONS IN DIRECT INTERPRETATION OF GEOELECTRICAL SOUNDINGS *

We start from the Hankel transform of Stefanescu's integral written in the convolutionintegral form suggested by Ghosh (1971). In this way it is possible to obtain the kernel function by the linear electric filter theory. Ghosh worked out the sets of filter coefficients in frequency domain and showed the very low content of high frequencies of apparent resistivity curves. Vertical soundings in the field measure a series of apparent resistivity values at a constant increment Δx of the logarithm of electrode spacing. Without loss of information we obtain the filter coefficient series by digital convolution of the Bessel function of exponential argument with sine function of the appropriate argument. With a series of forty-one values we obtain the kernel functions from the resistivity curves to an accuracy of better than 0.5%. With the digital method it is possible to calculate easily the filter coefficients for any electrode arrangement and any cut-off frequency.     

THE QUANTITATIVE INTERPRETATION OF DIPOLE SOUNDINGS BY MEANS OF THE RESISTIVITY TRANSFORM FUNCTION*

In this paper a theorem is demonstrated which allows—after the introduction of a suitable dipole kernel function or dipole resistivity transform function—to write the apparent resistivity function as an Hankel transformable integral expression. As a practical application of the theorem a procedure of quantitative interpretation of dipole soundings is suggested in which the dipole resistivity transform function obtained after inversion of the original dipole apparent resistivity data is used to control the goodness of the set of layering parameters which have been derived with our previous method of transformation of dipole sounding curves into equivalent Schlumberger diagrams.     

CONTRIBUTION DES MESURES DE L'ANISOTROPIE ÉLECTRIQUE À LA RECHERCHE DES AQUIFÈRES DE FRACTURE EN MILIEU CRISTALLIN AU BÉNIN1

A particular methodology adapted to crystalline formations with a thin weathered zone was developed for a village hydrological project, in Benin. A combination of electrical profiles, Schlumberger and square arrays, was able to locate the most fractured zones in the basement. We present some results obtained from theoretical models as well as from field data. The suggested methodology uses both measurements of resistivity and anisotropy. Strong anisotropy and low resistivity indicate the most productive hydrogeological areas.     

THE APPLICATION OF LINEAR FILTER THEORY TO THE DIRECT INTERPRETATION OF GEOELECTRICAL RESISTIVITY SOUNDING MEASUREMENTS*

Koefoed has given practical procedures of obtaining the layer parameters directly from the apparent resistivity sounding measurements by using the raised kernel function H(λ) as the intermediate step. However, it is felt that the first step of his method—namely the derivation of the H curve from the apparent resistivity curve—is relatively lengthy. In this paper a method is proposed of determining the resistivity transform T(λ), a function directly related to H(λ), from the resistivity field curve. It is shown that the apparent resistivity and the resistivity transform functions are linearily related to each other such that the principle of linear electric filter theory could be applied to obtain the latter from the former. Separate sets of filter coefficients have been worked out for the Schlumberger and the Wenner form of field procedures. The practical process of deriving the T curve simply amounts to running a weighted average of the sampled apparent resistivity field data with the pre-determined coefficients. The whole process could be graphically performed within an quarter of an hour with an accuracy of about 2%.     

随深度指数变化的视电阻率解

从横向均匀介质满足的基本方程出发,得到视电阻率核函数的一阶非线性微分方程,通过方程求解,并利用滤波系数法容易得到电阻率随深度任意变化的视电阻率问题。当各层介质电性结构随深度呈指数变化时,还可得到各层之间核函数的递推关系,这对实际介质的正反演问题都有重要意义和应用价值。     

NUMERICAL INTERPRETATION OF RESISTIVITY SOUNDINGS OVER HORIZONTAL BEDS*

The digital computer technique described for interpreting resistivity soundings over a horizontally stratified earth requires two steps. First, the kernel function is evaluated numerically from the inverse Hankel transform of the observed apparent resistivity curve. Special attention is given to the inversion of resistivity data recorded over a section with a resistant basement. The second step consists in the least-squares estimation of layer resistivities and thicknesses from the kernel function. For the case of S or T-equivalent beds only one layer-parameter can be obtained, either the longitudinal conductance, or the transverse resistance respectively. Two examples given in the paper show that a wide tolerance is permitted for Choosing the starting values of the layering parameters in the successive approximation procedure. Another important feature for practical applications is good convergence of the iterations. The method is probably best suited for interpreting profiles of electrical soundings with the purpose of mapping approximately horizontal interfaces at depth.     

RESISTIVITY SOUNDING ON A MULTI-LAYERED EARTH WITH TRANSITIONAL LAYERS. PART I: THEORY*

The electrical potential generated by a point source of current on the ground surface is studied for a multi-layered earth formed by layers alternatively characterized by a constant conductivity value and by conductivity varying linearly with depth. The problem is accounted for by solving a Laplace's differential equation for the uniform layers and a Poisson's differential equation for the transitional layers. Then, by a simple algorithm and by the introduction of a suitable kernel function, the general expression of the apparent resistivity for a Schlumberger array placed on the surface is obtained. Moreover some details are given for the solution of particular cases as 1) the presence of a infinitely resistive basement, 2) the absence of any one or more uniform layers, and 3) the absence of any one or more transitional layers. The new theory proves to be rather general, as it includes that for uniform layers with sharp boundaries as a particular case. Some mathematical properties of the kernel function are studied in view of the application of a direct system of quantitative interpretation. Two steps are considered for the solution of the direct problem: (i) The determination of the kernel function from the field measurements of the apparent resistivity. Owing to the identical mathematical formalism of the old with this new resistivity theory, the procedures there developed for the execution of the first step are here as well applicable without any change. Thus, some graphical and numerical procedures, already published, are recalled. (ii) The determination of the layer distribution from the kernel function. A recurrent procedure is proposed and studied in detail. This recurrent procedure follows the principle of the reduction to a lower boundary plane, as originally suggested by Koefoed for the old geoelectrical theory. Here the method differs mainly for the presence of reduction coefficients, which must be calculated each time when passing to a reduced earth section.     

CARACTERES DE ĽONDE ELECTROMAGNETIQUE DE SURFACE ENGENDREE PAR UN DIPOLE MAGNETIQUE *

The electromagnetic field radiated from a magnetic dipole lying on the ground is considered, in the extremely low frequency range (DC to 20 000 Hertz). Theoretical and experimental data are given on the characteristics of the surface wave (vanishing wave) generated at the air-ground interface, in the case of an homogeneous subsurface. The case of a subsurface with electrical resistivity varying with depth is considered. It is shown how the above-mentioned characteristics may be applied in the quantitative investigation of the electrical resistivity of the subsurface as a function of depth, in a method using the measurements on the ground of all the components of the radiated field: horizontal electrical components, vertical and horizontal magnetic components.     

A NUMERICAL METHOD TO COMPUTE THE RESISTIVITY TRANSFORM FROM WENNER SOUNDING DATA*

A numerical technique to compute the resistivity transform directly from the observed Wenner sounding data has been developed. In principle, the procedure is based on a decomposition method and consists of two steps: the first step determines a function that approximates the apparent resistivity data and the second step transforms this function into the corresponding kernel by an analytical operation. The proposed method is tested on some theoretical master curves. A high degree of precision is achieved with very little computer time. The applicability is shown on two field examples.     

水平层状介质视电阻率的高精度计算公式及其误差分析

为适应地电阻率法预报地震深入研究的需要,本文介绍了几种高精度计算水平层状介质视电阻率的公式及方法。对于MN→0的视电阻率计算,本文给出了J_1(x)滤波系数法的计算步骤,此外对二层、三层介质的传统级数求和法,给出了估算精度的公式;对于MN≠O的视电阻率计算,本文讨论了数值积分法和J_0(x)滤波系数法的计算步骤以及为提高精度所采取的技术措施,并且给出了一批用两种方法所得到的计算结果。