THE USE OF TWO-ELECTRODE AND SCHLUMBERGER FILTERS FOR COMPUTING RESISTIVITY AND EM SOUNDING CURVES*

Different sets of filter coefficients for the linear filter technique for the computations of resistivity and EM sounding curves are evaluated for several electrode and coil configurations. Instead of this procedure, the two-electrode filter can be used for computations of Wenner, Schlumberger, and dipole—dipole apparent resistivity model curves by defining convolutional expressions which contain the new input functions in terms of the resistivity transform function. Similarly, the Schlumberger filter performs the computations of dipole—dipole apparent resistivity model curves. The Wenner, Schlumberger, and dipole—dipole filter functions are defined in terms of the two-electrode filter using the new convolutional expressions. A relationship between the Schlumberger and dipole—dipole filter functions is given. The above arguments are adopted for the computations of EM sounding curves. It is shown that the EM filter for the horizontal coplanar loop system (which is identical to the two-electrode filter) performs the computations of the mutual coupling ratios for perpendicular, vertical coplanar, and vertical coaxial loop systems. In the same way, the Schlumberger filter can be used to compute vertical coaxial sounding curves. The corresponding input functions are defined in terms of the EM kernel for all convolutional expressions presented. After these considerations, integral expressions of the mutual coupling ratios involving zero-order Bessel function are derived. The mutual coupling ratio for the vertical coaxial loop system is given in the same form as the mutual coupling ratio for the vertical coplanar loop system.     

ON THE TRANSFORMATION OF DIPOLE TO SCHLUMBERGER SOUNDING CURVES *

A method is described to transform a dipole sounding curve, obtained with any one of the common dipole arrays over a horizontally layered earth, to the form of a Schlumberger sounding curve. Starting from the general expression which relates the dipole apparent resistivity to the Schlumberger apparent resistivity and its derivative with respect to the spacing, it is possible with some approximations to derive an easy numerical computation procedure in order to perform the transformation. The applicability of the method is discussed briefly.     

TRANSFORMATION OF DIPOLE TO SCHLUMBERGER SOUNDING CURVES BY MEANS OF DIGITAL LINEAR FILTERS*

Linear relationship between dipole and Schlumberger sounding resistivities leads to the use of digital filters to transform the former to the latter. This transformation is of importance from the viewpoint that Schlumberger interpretational techniques and know-how could then be applied to the pseudo-Schlumberger field curve. Filters for this transformation are presented for the radial, perpendicular, and parallel (30°) dipole method. The characteristics of these filters are similar to the ones for transforming dipole data to the kernel and are favourable in that they do not amplify noise. A sampling interval of (In 10) /6 has been used in determining the filter yielding good accuracy. Like previous filters the present one is handy and fast in operation.     

TRANSFORMATION OF SCHLUMBERGER APPARENT RESISTIVITY TO DIPOLE APPARENT RESISTIVITY OVER LAYERED EARTH BY THE APPLICATION OF DIGITAL LINEAR FILTERS *

A simple extension of our previous work in which digital filters were developed to transform dipole resistivity measurements over layered earth to Schlumberger ones leads us to the development of filters for transforming the latter to the former. As in the previous work we use a sampling interval of 1/6 In 10 in designing the filters that are both accurate and fast in operation.     

DIGITAL LINEAR FILTER FOR COMPUTING TYPE CURVES FOR THE TWO-ELECTRODE SYSTEM OF RESISTIVITY SOUNDING*

In this paper a technique for computing type curves for the two-electrode resistivity soundings is presented. It is shown that the apparent resistivity due to the system can be represented by a convolution integral. Thus, it is possible to apply the principle of digital linear filtering and compute the desired type-curves. The filter function required for the purpose is found to be identical with that used to compute the EM sounding curves for the two coplanar horizontal loop system. It is further shown that from the two-electrode apparent-resistivity expression one can easily derive the apparent resistivities for other configurations. A comparison of depths of investigation for various systems reveals that the two-electrode system has greater depth of investigation than other conventional systems. This is also supported by the field example presented in the end, which illustrates the relative performance of the two-electrode system vis-a-vis the Wenner system.     

FINITE-DIFFERENCE EVALUATION OF APPARENT RESISTIVITY CURVES*

The problem of numerical evaluation of apparent resistivity curves is treated by finite difference modeling. The models proposed are set up in cylindrical coordinates and yield the potential field due to a point source located in a radially symmetric environment. The Schlumberger configuration, widely used for surface measurements, is emphasized. However, the treatment is equally applicable to other similar situations such as the computation of synthetic electric logs when the resistivity of the borehole fluid is different from that of the surrounding uniform or stratified medium. Moreover, the individual layers may not necessarily be isotropic. The medium under investigation is discretized by using a very coarse system of horizontal and vertical grid lines whose distance from the source increases logarithmically; consequently, the physical dimensions of the medium can be made “infinite” without affecting the numerical size of the model. Finer features such as a thin but anomalously resistive or conductive bed which would ordinarily be missed in coarse discretization are accurately taken into account, since the calculations are done in terms of the Dar Zarrouk parameters derived from the exact resistivity distribution of the model. This enables one to compute the potential field by inverting a small sparse matrix. When the medium comprises only a few layers, the efficiency of the finite-difference model is comparable to that of the known analytical methods; for more complicated structures, however, the finite-difference model becomes more efficient. The accuracy of finite-difference results is demonstrated by comparing them with the corresponding analytically obtained data.     

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.     

LOW-PASS FILTERING OF NOISY SCHLUMBERGER SOUNDING CURVES PART 1: THEORY*

A contribution is given to the solution of the problem of filtering noise-degraded Schlumberger sounding curves. It is shown that the transformation to the pole-pole system is actually a smoothing operation that filters high-frequency noise. In the case of residual noise contamination in the transformed pole-pole curve, it is demonstrated that a subsequent application of a conventional rectangular low-pass filter, with cut-off frequency not less than the right-hand frequency limit of the main message pass-band, may satisfactorily solve the problem by leaving a pole-pole curve available for interpretation. An attempt is also made to understand the essential peculiarities of the pole-pole system as far as penetration depth, resolving power and selectivity power are concerned.     

电阻率测深的数字解释

本文主要介绍了应用积分变换的方法和采样定理将视电阻率ρs曲线作线性滤波,得出一新的电阻率转换函数T′曲线,然后,以层参数（各层的电阻率和厚度）算得的T用最优化数值方法在DJS-6型电子计算机上与其进行自动拟合,以达到解释电测深曲线的目的。 文中简述了戈什（Ghosh）提出的对ρs作线性滤波的原理,介绍了与国外不同的取样间距和滤波系数的确定以及阻尼最小二乘法和变尺度最优化法的计算框图和应用,最后附有实例和简要的讨论。     

LOW-PASS FILTERING OF NOISY FIELD SCHLUMBERGER SOUNDING CURVES PART II: APPLICATION*

The basic principles of the application of the linear system theory for smoothing noise-degraded d.c. geoelectrical sounding curves were recently established by Patella. A field Schlumberger sounding is presented to demonstrate first their application and validity. To achieve this purpose, firstly it is pointed out that the required smoothing or low-pass filtering can be considered as an intrinsic property of the transformation of original Schlumberger sounding curves into pole-pole (two-electrode) curves. Then we sketch a numerical algorithm to perform the transformation, opportunely modified from a known procedure for transforming dipole diagrams into Schlumberger ones. Finally we show a field example with the double aim of demonstrating (i) the high quality of the low-pass filtering, and (ii) the reliability of the transformed pole-pole curve as far as quantitative interpretation is concerned.     

电阻率测深的数字解释

本文主要介绍了应用积分变换的方法和采样定理将视电阻率ρ_s曲线作线性滤波,得出一新的电阻率转换函数T′曲线,然后,以层参数（各层的电阻率和厚度）算得的T用最优化数值方法在DJS-6型电子计算机上与其进行自动拟合,以达到解释电测深曲线的目的。 文中简述了戈什（Ghosh）提出的对ρ_s作线性滤波的原理,介绍了与国外不同的取样间距和滤波系数的确定以及阻尼最小二乘法和变尺度最优化法的计算框图和应用,最后附有实例和简要的讨论。     

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 NUMERICAL METHOD OF CALCULATING THE KERNEL FUNCTION FROM SCHLUMBERGER APPARENT RESISTIVITY DATA*

A method to calculate the resistivity transform of Schlumberger VES curves has been developed. It consists in approximating the field apparent resistivity data by utilizing a linear combination of simple functions, which must satisfy the following requirements: (i) they must be suitable for fitting the resistivity data; (ii) once the fitting function has been obtained they allow the kernel to be determined in an analytic way. The fitting operation is carried out by the least mean squares method, which also accomplishes a useful smoothing of the field curve (and therefore a partial noise filtering). It gives the possibility of assigning different weights to the apparent resistivity values to be approximated according to their different reliability. For several examples (theoretical resistivity curves in order to estimate the precision of the method and with field data to verify the practicality) yield good results with short execution time independent of shape the apparent resistivity curve.     

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.     

DEPTH OF INVESTIGATION IN WENNER,THREE-ELECTRODE AND DIPOLE-DIPOLE DC RESISTIVITY METHODS*

Using a method and definition given earlier (Roy and Apparao, 1971), this paper computes the depths of investigation in homogeneous ground for (a) the Wenner α, β and γ configurations, (b) the three electrode system and (c) the dipole-dipole arrangements when the dipole lengths are not infinitesimally small. The results for (a) and (b) have been summarised in a table, while those for (c) are shown as contour diagrams. In all the dipolar arrangements examined in this paper, except the equatorial, the depth of investigation decreases (and the vertical resolution increases) with increase in any or both of the dipole lengths. For the equatorial set up, this decrease (or increase) is very small.     

电偶源频率电磁测深中的静位移现象及视电阻率曲线的畸变

以电偶源频率电磁测深 (FEMS)二维数字模拟为例 ,从理论上简略讨论了波区定义的 5种视电阻率的变化特征及静位移现象 ,定性分析了表面局部水平非均匀体对视电阻率曲线畸变的影响 .     

长偏移距瞬变电磁测深甚晚期响应及视电阻率的数值计算

文中讨论了提高长偏移距瞬变电磁测深甚晚期响应计算精度及相应的全区视电阻率的计算方法。在频率域 ,利用直接数值积分结合连分式展开提高低频段响应的计算精度 ,从而提高时间域甚晚期响应的计算精度。同时 ,利用连分式渐进展开 ,将均匀半空间电阻率表示为该模型长偏移距瞬变电磁测深响应的反函数 ,利用迭代方法求出适合全时间段的全区视电阻率。数值结果表明 ,文中的方法可有效地提高长偏移距瞬变电磁测深 (甚 )晚期数据的模拟与解释精度     

UNCERTAIN RESISTIVITY SOUNDING AND EQUIVALENT MODELS*

Problems of the resistivity sounding method may be reparametrized by means of Dar Zarrouk—or generalized D.Z.—transformations. A complete set of equivalent models is subdivided into three classes: initial models, intermediate models, and all other models. Equivalent reparametrized models are subdivided into three analogous classes. There is no one-to-one mapping of intermediate primary and intermediate reparametrized models. Non-trivial equivalent models are generated from initial models via intermediate reparametrized models. A method is developed to find equivalent models relative to the uncertainty of resistivity soundings. There are different forms of approximate equivalence laws, apart from those commonly accepted. The equivalence depends considerably on the “flatness” of the sounding curve. For mean flatness, the electric models exhibit the (hp^±2)-equivalence rather than the (hp^±1)-equivalence.     

TRANSFORMATION OF RESISTIVITY SOUNDING MEASUREMENTS OBTAINED IN ONE ELECTRODE CONFIGURATION TO ANOTHER CONFIGURATION BY MEANS OF DIGITAL LINEAR FILTERING*

The technique of digital linear filtering is used for transformation of apparent resistivity data from one electrode configuration into another. Usually filter spectra are determined via the discrete Fourier transforms of input and output functions: the filter characteristic is the quotient of the spectra of the output function and input function. In this paper, the transformation of the apparent resistivities is presented for four electrode configurations (Wenner, the two-electrode, Schlumberger, and dipole configurations). In our method, there is no need to use the discrete Fourier transform of the input and output functions in order to determine the filter spectrum for converting apparent resistivity in one electrode configuration to any other configuration. Sine responses for determination of the derivative of apparent resistivities are given in analytical form. If the filter spectrum for converting the apparent resistivity to the resistivity transform for one electrode configuration is known, the filter spectra for transforming the apparent resistivity to the resistivity transform for any electrode configurations can be calculated by using newly derived expressions.