TRAITEMENT AUTOMATIQUE DES SONDAGES ELECTRIQUES AUTOMATIC PROCESSING OF ELECTRICAL SOUNDINGS*

I. Introduction In this section the problem is stated, its physical and mathematical difficulties are indicated, and the way the authors try to overcome them are briefly outlined. Made up of a few measurements of limited accuracy, an electrical sounding does not define a unique solution for the variation of the earth resistivities, even in the case of an isotropic horizontal layering. Interpretation (i.e. the determination of the true resistivities and thicknesses of the ground-layers) requires, therefore, additional information drawn from various more or less reliable geological or other geophysical sources. The introduction of such information into an automatic processing is rather difficult; hence the authors developped a two-stage procedure:

• a) the field measurements are automatically processed, without loss of information, into more easily usable data;
• b) some additional information is then introduced, permitting the determination of several geologically conceivable solutions.

The final interpretation remains with the geophysicist who has to adjust the results of the processing to all the specific conditions of his actual problem. II. Principles of the procedure In this section the fundamental idea of the procedure is given as well as an outline of its successive stages. Since the early thirties, geophysicists have been working on direct methods of interpreting E.S. related to a tabular ground (sequence of parallel, homogeneous, isotropic layers of thicknesses h_i and resistivities ρ_i). They generally started by calculating the Stefanesco (or a similar) kernel function, from the integral equation of the apparent resistivity: where r is the distance between the current source and the observation point, S₀ the Stefanesco function, ρ(z) the resistivity as a function of the depth z, J₁ the Bessel function of order 1 and λ the integration variable. Thicknesses and resistivities had then to be deduced from S₀ step by step. Unfortunately, it is difficult to perform automatically this type of procedure due to the rapid accumulation of the errors which originate in the experimental data that may lead to physically impossible results (e.g. negative thicknesses or resistivities) (II. 1). The authors start from a different integral representation of the apparent resistivity: where K₁ is the modified Bessel function of order I. Using dimensionless variables t = r/2h₀ and y(t)=ζ (r)/ρ₁ and subdividing the earth into layers of equal thicknesses h₀ (highest common factor of the thicknesses h_i), ø becomes an even periodic function (period 2π) and the integral takes the form: The advantage of this representation is due to the fact that its kernel ø (function of the resistivities of the layers), if positive or null, always yields a sequence of positive resistivities for all values of θ and thus a solution which is surely convenient physically, if not geologically (II.3). Besides, it can be proved that ø(θ) is the Fourier transform of the sequence of the electric images of the current source in the successive interfaces (II.4). Thus, the main steps of the procedure are: a) determination of a non-negative periodic, even function ø(θ) which satisfies in the best way the integral equation of apparent resistivity for the points where measurements were made; b) a Fourier transform gives the electric images from which, c) the resistivities are obtained. This sequence of resistivities is called the “comprehensive solution”; it includes all the information contained in the original E.S. diagram, even if its too great detail has no practical significance. Simplification of the comprehensive solution leads to geologically conceivable distributions (h, ρ) called “particular solutions”. The smoothing is carried out through the Dar-Zarrouk curve (Maillet 1947) which shows the variations of parameters (transverse resistance R_i= h_i.ρ_i–as function of the longitudinal conductance C_i=h_i/ρ_i) well suited to reflect the laws of electrical prospecting (principles of equivalence and suppression). Comprehensive and particular solutions help the geophysicist in making the final interpretation (II.5). III. Computing methods In this section the mathematical operations involved in processing the data are outlined. The function ø(θ) is given by an integral equation; but taking into account the small number and the limited accuracy of the measurements, the determination of ø(θ) is performed by minimising the mean square of the weighted relative differences between the measured and the calculated apparent resistivities: minimum with inequalities as constraints: where t_l are the values of t for the sequence of measured resistivities and p_l are the weights chosen according to their estimated accuracy. When the integral in the above expression is conveniently replaced by a finite sum, the problem of minimization becomes one known as quadratic programming. Moreover, the geophysicist may, if it is considered to be necessary, impose that the automatic solution keep close to a given distribution (h, ρ) (resulting for instance from a preliminary interpretation). If φ(θ) is the ø-function corresponding to the fixed distribution, the quantity to minimize takes the form: where: The images are then calculated by Fourier transformation (III.2) and the resistivities are derived from the images through an algorithm almost identical to a procedure used in seismic prospecting (determination of the transmission coefficients) (III.3). As for the presentation of the results, resorting to the Dar-Zarrouk curve permits: a) to get a diagram somewhat similar to the E.S. curve (bilogarithmic scales coordinates: cumulative R and C) that is an already “smoothed” diagram where deeper layers show up less than superficial ones and b) to simplify the comprehensive solution. In fact, in arithmetic scales (R versus C) the Dar-Zarrouk curve consists of a many-sided polygonal contour which múst be replaced by an “equivalent” contour having a smaller number of sides. Though manually possible, this operation is automatically performed and additional constraints (e.g. geological information concerning thicknesses and resistivities) can be introduced at this stage. At present, the constraint used is the number of layers (III.4). Each solution (comprehensive and particular) is checked against the original data by calculating the E.S. diagrams corresponding to the distributions (thickness, resistivity) proposed. If the discrepancies are too large, the process is resumed (III.5). IV. Examples Several examples illustrate the procedure (IV). The first ones concern calculated E.S. diagrams, i.e. curves devoid of experimental errors and corresponding to a known distribution of resistivities and thicknesses (IV. 1). Example I shows how an E.S. curve is sampled. Several distributions (thickness, resistivity) were found: one is similar to, others differ from, the original one, although all E.S. diagrams are alike and characteristic parameters (transverse resistance of resistive layers and longitudinal conductance of conductive layers) are well determined. Additional informations must be introduced by the interpreter to remove the indeterminacy (IV.1.1). Examples 2 and 3 illustrate the principles of equivalence and suppression and give an idea of the sensitivity of the process, which seems accurate enough to make a correct distinction between calculated E.S. whose difference is less than what might be considered as significant in field curves (IV. 1.2 and IV. 1.3). The following example (number 4) concerns a multy-layer case which cannot be correctly approximated by a much smaller number of layers. It indicates that the result of the processing reflects correctly the trend of the changes in resistivity with depth but that, without additional information, several equally satisfactory solutions can be obtained (IV. 1.4). A second series of examples illustrates how the process behaves in presence of different kinds of errors on the original data (IV.2). A few anomalous points inserted into a series of accurate values of resistivities cause no problem, since the automatic processing practically replaces the wrong values (example 5) by what they should be had the E.S. diagram not been wilfully disturbed (IV.2.1). However, the procedure becomes less able to make a correct distinction, as the number of erroneous points increases. Weights must then be introduced, in order to determine the tolerance acceptable at each point as a function of its supposed accuracy. Example 6 shows how the weighting system used works (IV.2.2). The foregoing examples concern E.S. which include anomalous points that might have been caused by erroneous measurements. Geological effects (dipping layers for instance) while continuing to give smooth curves might introduce anomalous curvatures in an E.S. Example 7 indicates that in such a case the automatic processing gives distributions (thicknesses, resistivities) whose E.S. diagrams differ from the original curve only where curvatures exceed the limit corresponding to a horizontal stratification (IV.2.3). Numerous field diagrams have been processed (IV. 3). A first case (example 8) illustrates the various stages of the operation, chiefly the sampling of the E.S. (choice of the left cross, the weights and the resistivity of the substratum) and the selection of a solution, adapted from the automatic results (IV.3.1). The following examples (Nrs 9 and 10) show that electrical prospecting for deep seated layers can be usefully guided by the automatic processing of the E.S., even when difficult field conditions give original curves of low accuracy. A bore-hole proved the automatic solution proposed for E.S. no 10, slightly modified by the interpreter, to be correct.     

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%.     

AN INTERPRETATION THEORY FOR INDUCED POLARIZATION VERTICAL SOUNDINGS (TIME-DOMAIN)*

In this paper it is shown how one may obtain a generalized Ohm's law which relates the induced polarization electric field to the steady-state current density through the introduction of a fictitious resistivity defined as the product of the chargeability and the resistivity of a given medium. The potential generated by the induced polarization is calculated at any point in a layered earth by the same procedure as used for calculating the potential due to a point source of direct current. On the basis of the definition of the apparent chargeability m_a, the expressions of m_a for different stratigraphie situations are obtained, provided the IP measurements are carried out on surface with an appropriate AMNB array. These expressions may be used to plot master curves for IP vertical soundings. Finally some field experiments over sedimentary formations and the quantitative interpretation procedure are reported.     

ZERO—PHASE FORWARD FILTERS FOR RESISTIVITY SOUNDING*

Forward filters to transform the apparent resistivity function over a layered half-space into the resistivity transform have been derived for a number of sample intervals. The filters have no apparent Gibbs' oscillations and hence require no phase shift. In addition, the end points of the filter were modified to compensate for truncation. The filters were tested on simulated ascending and descending two-layer cases. As expected, “dense” filters with sample spacing of In (10)/6 or smaller performed very well. However, even “sparse” filters with spacing of In (10)/2 and a total of nine coefficients have peak errors of less than 5% for p₁:p₂ ratios of 10^–6 to 10⁶. If a peak error of 5.5% is acceptable, then an even sparser filter with only seven coefficients at a spacing of 3 In (10)/5 may be used.     

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.     

AN INTERACTIVE COMPUTER/GRAPHIC-DISPLAY-TERMINAL SYSTEM FOR INTERPRETATION OF RESISTIVITY SOUNDINGS*

A fast computer-procedure giving the apparent resistivity curve as well as the partial derivatives with respect to the layer-parameters is presented. It is based on the linear filter method developed by D. P. Ghosh in 1971. The sampling frequency is 10 points per decade, and 3 decades are covered. The maximum relative error is less than 10^?3, and in most cases orders of magnitude smaller. The computation time on a CDC 6400 for one curve given in 30 points ranges linearly from .17s for a two-layer case to .36s for a ten-layer case. The procedure is used to plot master curves interactively on a graphic display terminal (Tektronix 4010) connected to the CDC 6400. By trial-and-error adjustments a set of layer-parameters is found, giving essentially the measured curve.     

NEW WAYS FOR THE INTERPRETATION OF GEOELECTRICAL RESISTIVITY MEASUREMENTS IN THE SEARCH FOR AND DELIMITATION OF AQUIFERS

Abstract

Geoelectric resistivity measurements by means of direct current for solving hydrogeological problems have become increasingly significant in recent years. Measurements on the surface according to the four-point-method (SCHLUMBERGER or WENNER arrangement) result in so-called “apparent” resistivities ? α as a function of the electrode distance L. The evaluation of these measuring data, here in form of sounding graphs ? α(L/2), consists of the determination of true resistivities as a function of the depth z. Since a direct computation of ? (z) from α (L/2) is not possible in practice, theoretically computed master curves constitute the essential auxiliary means for the evaluation.

New simplified calculation techniques allow to establish accurately computed master curves for an underground consisting of more than three layers. By means of such standard graphs special problems of hydrogeology can quantitatively be solved by applying geoelectrical methods. The procedure is demonstrated on hand of complicated cases of aquifers devided into several storeys.     

RESISTIVITY MEASUREMENTS IN VALLEYS WITH ELLIPTIC CROSS-SECTION*

In this paper an idealized valley of a semi-elliptic cross-section is considered. For a Schlumberger configuration on the axis, sets of master curves are calculated for the ratio of semi-axis a/b= 1, 2, 3, 4, 5 in corresponding to various resistivity ratios of surrounding rocks and valley sediments. For small resistivity ratios, these model curves have the shape of three-layer curves for horizontal bedding and are often equivalent to them within the accuracy of measurements. The axial ratio a/b considerably affects the depth determination of valley sediments. In the special case of a circular cross-section (a/b= 1) the influence of the position of the electrodes on the sounding curve is studied in more detail. The application of the master curves in practice shows that the influence of the specific shape of the valley on soundings should not be neglected. In general, the valleys have a greater “true” depth than can be seen from the interpretation of a sounding by master curves corresponding to the horizontal bedding.     

GEOELECTRICAL MODEL CALCULATIONS FOR TWO-DIMENSIONAL RESISTIVITY DISTRIBUTIONS*

For the calculation of geoelectrical model curves for a two-dimensional resistivity distribution, the potential equation is transformed by means of a Fourier cosine transform into a two-dimensional Helmholtz equation containing the separation parameter λ. The numerical solution of this equation for different values of λ for an irregular grid is obtained using the method of finite differences combined with the method of overrelaxation. The method by which derivatives are replaced by finite differences turned out to be very important, especially for high resistivity contrasts. After testing several methods designed to deal with any type of resistivity distribution, a method of discretization similar to that used by Brewitt—Taylor and Weaver (1976) for magnetotelluric modeling for H polarization was found the best. Examples are given of model curves for Schlumberger soundings over a vertical fault covered by overburden. The incorrect use of horizontal-layer models leads to erroneous interpretations that are more complex than the real subsurface situations.     

STACKING FILTERS AND THEIR CHARACTERIZATION IN THE (f—k) DOMAIN *

The implementation of a stacking filter involves the filtering of each trace with an individual filter and the subsequent summing of all outputs. The actual position of a trace in space as well as certain simultaneous shifts of traces and filter components in time do not influence the process. The resulting output is consequently invariant to various arbitrary coordinate transformations. For a certain useful class of ensembles of non-linear moveout arrival times for signals a particular transformation can be found which transforms a given ensemble into one consisting only of straight lines. It is thus possible to reduce, for instance, the analysis of a stacking filter designed for hyperbola-like moveout curves to the analysis of a velocity filter with linear moveout curves. As the (f—k) transform is a very useful concept to describe a velocity filter, it can consequently be applied to characterize a stacking filter in regard to its performance on input signals with non-linear moveout.     

MATHEMATICAL MODELS OF CONDUCTING ORE BODIES FOR DIRECT CURRENT ELECTRICAL PROSPECTING *

The authors generalize a method expounded in a previous paper (1971, Geoph. Prosp. 18, 786-799) to the case of a local conductivity σ(M) of the infinite medium satisfying the relation where the R_i's are the distances from the point M to n fixed points S_i (i= 1,. n), k is a positive real constant and C_i, C_ii are constants ensuring the condition α > O. The sub-surface conductivity distributions (half-spaces) complying with (1) provide a wide variety of conducting structures, which can fit quite successfully the rather complicated distributions of conductivity occurring in natural ore bodies. An exact algebraic calculation of the apparent resistivity for these grounds, valid for any dc electrical prospecting devices (Wenner, Schlumberger, dipole, etc.) leads to a set of simultaneous linear equations, with a matrix which is invariant with respect to the position of the quadrupole being used. This greatly simplifies the numerical computation. We also present some examples of cross sections for the real and apparent resistivity obtained by this method.     

THE DETERMINATION OF FILTER COEFFICIENTS FOR THE COMPUTATION OF STANDARD CURVES FOR DIPOLE RESISTIVITY SOUNDING OVER LAYERED EARTH BY LINEAR DIGITAL FILTERING *

The technique of linear digital filtering developed for the computation of standard curves for conventional resistivity and electromagnetic depth soundings is applied to the determination of filter coefficients for the computation of dipole curves from the resistivity transform function by convolution. In designing the filter function from which the coefficients are derived, a sampling interval shorter than the one used in the earlier work on resistivity sounding is found to be necessary. The performance of the filter sets is tested and found to be highly accurate. The method is also simple and very fast in application.     

THE QUANTITATIVE ANALYSIS OF FREQUENCY-DOMAIN INDUCED POLARIZATION SOUNDINGS OVER HORIZONTAL BEDS*

Following up our recent study of an indirect procedure for the practical determination of the maximum frequency-effect, defined as fe = 1 ? pρ_∞/ρdc with ρ_∞ the resistivity at infinite frequency, we show at first how, through the Laplace transform theory, ρ_∞ can be related to stationary field vectors in the simple form of Ohm's law. Then applying the equation of continuity for stationary currents with a suitable set of boundary conditions, we derive the integral expression of the apparent resistivity at infinite frequency ρ_∞,a in the case of a horizontally layered earth. Finally, from the definition of the maximum apparent frequency-effect, analytical expressions of fe_α are obtained for both Schlumberger and dipole arrays placed on the surface of the multi-layered earth section in the most general situation of vertical changes in induced polarization together with dc resistivity variations not at the same interfaces. Direct interpretation procedures are suggested for obtaining the layering parameters directly from the analysis of the sounding curves.     

ESSAIS D'AMELIORATION DU TRAITEMENT ET DE L'INTERPRETATION DES ENREGISTREMENTS MAGNETOTELLURIQUES*

The MT method has proved to be a useful tool for the resolution of some geophysical problems not well adapted to seismic prospecting. To meet the accuracy needed by such cases, the original field and processing techniques has to be somewhat improved. The processing of the natural EM field (with two electrical and three magnetic components recorded) is made in the frequency domain, and the usual MT parameters are calculated. However, the dispersion of the results is given a great importance and the final plots (apparent resistivity vs. frequency, azimuths of the main directions vs. frequency, etc.) illustrate the distribution of the calculated values. Average curves can then be selected and studied by means of the known master curves. However, if the earth can be assumed to consist of parallel layers, a computer program derives from the apparent resistivity diagram an automatic solution consisting of a great number of layers with the same ratio of thickness to square root of resistivity. This comprehensive solution is then simplified through an auxiliary diagram called the MT Dar Zarrouk curve and a set of geologically consistent solutions can be found.     

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.     

用双侧向测井响应求取地层参数的Newton－SVD反演方法

针对测井中经常遇到的完全非均匀的地层模型，利用深、浅双侧向测井的视电阻率曲线作为约束条件，给出了反演侵入半径ｒ_ｉ、原状地层电阻率Ｒ_ｔ、上（下）围岩电阻率Ｒ_{ｓｕ（Ｒｓｄ）的Ｎｅｗｔｏｎ－ＳＶＤ反演方法．数值模拟表明这一方法是可行的，经实际资料处理．说明该方法实用．}     

Electric Potential and Fréchet Derivatives for a Uniform Anisotropic Medium with a Tilted Axis of Symmetry

In this paper we develop analytic solutions for the electric potential, current density and Fréchet derivatives at any interior point within a 3-D transversely isotropic medium having a tilted axis of symmetry. The current electrode is assumed to be on the surface of the Earth and the plane of stratification given arbitrary strike and dip. Profiles can be computed for any azimuth. The equipotentials exhibit an elliptical pattern and are not orthogonal to the current density vectors, which are strongly angle dependent. Current density reaches its maximum value in a direction parallel to the longitudinal conductivity direction. Illustrative examples of the Fréchet derivatives are given for the 2.5-D problem, in which the profile is taken perpendicular to strike. All three derivatives of the Green’s function with respect to longitudinal conductivity, transverse resistivity and dip angle of the symmetry axis (dG/dσ_l, dG/dσ_t, dG/dθ₀) show a strongly asymmetric pattern compared to the isotropic case. The patterns are aligned in the direction of the tilt angle. Such sensitivity patterns are useful in real-time experimental design as well as in the fast inversion of resistivity data collected over an anisotropic earth.     

Application of a conceptual method for separating runoff components in daily hydrographs in Kimakia forest catchments,Kenya

This paper highlights the use of a conceptual method for separating runoff components in daily hydrographs, contrary to the traditionally used graphical method of separation. In the conceptual method, the components, viz. surface flow, interflow and baseflow, are regarded as high, medium and low frequency signals and their separation is done using the principle of a recursive digital filter commonly used in signal analysis and processing. It requires estimates of the direct runoff (β_d) and surface runoff (β_s) filter parameters which are obtained by a least‐squares procedure involving baseflow and interflow indices based on graphical and recursive digital filter estimation techniques. The method thus circumvents the subjective element associated with the graphical procedure of hydrograph separation, in which case the eye approximation and/or one's skill at plotting is the prime basis for the whole analysis. The analysis based on three forest catchments in Kimakia, Kenya, East Africa, revealed that β_d=K_b and β_s=K_i , where K_b and K_i are the baseflow and interflow recession constants. Copyright © 1999 John Wiley & Sons, Ltd.     

A compact representation of magnetotelluric responses for two-layer models

In this paper the locations where ρ_app = ρ₁ and ? = π/4 and where these parameters reach an extreme value in two-layer magnetotelluric (MT) sounding curves are summarized in an extremely compact form. The key parameters over two-layer models with conductivities σ₁, σ₂ and upper layer thickness h are the real S and α, where S is the conductivity contrast and α is the distance between the observation site and the conductivity interface, normalized to the half skindepth in the first layer. If the impedance components, various resistivity definitions ( ρ_{Re Z}, ρ_{Im Z} and ρ_|Z|, based on different parts of the complex impedance Z ) and the magnetotelluric phase ? are derived as a function of S and α, then the conditions for the apparent resistivity ρ_app and the phase ? are that they either satisfy ρ_app = ρ₁ and ? = π/4 or attain extreme values which can be given in terms of simple algebraic equations between S and α. All equations are valid for observation sites at any depth 0 ≤ z ≤ h in the first layer. The set of equations, presented in a tabular form, may make it possible to determine a layer boundary from the short period part of the sounding curves, in particular the ρ_{Re Z} and the ?_MT curves.     

FUNDAMENTAL FUNCTIONS FOR HORIZONTALLY STRATIFIED EARTH*

The potential distribution and the wave propagation in a horizontally stratified earth is considered and the analogy of the mathematical expression for seismic transfer function, electromagnetic and electric kernel functions, and magnetotelluric input impedance is discussed. Although these specific functions are conveniently treated by a separate expression in each method, it is indicated that the function for seismic and electromagnetic methods is mathematically the same with a change in the physical meaning of the variables from one method to the other. Similarly, the identity of the mathematical expressions of the resistivity kernel function and magnetotelluric input impedance is noticed. In each method a specific geophysical function depends on the thickness and the physical properties of the various layers. Every specific function involves two interdependent fundamental functions, that is P_n and Q_n, or P_n and P*_n, having different physical meaning for different methods. Specific functions are expressible as a ratio P_n/Q_n or P*_n/P_n. Fundamental functions may be reduced to polynomials. The fundamental polynomials Q*_n and P*_n describing the horizontally stratified media are a system of polynomials orthogonal on the unit circle, of first and second order, respectively. The interpretation of geophysical problems corresponds to the identification of the parameters of a system of fundamental orthogonal polynomials. The theorems of orthogonal polynomials are applied to the solution of identification problems. A formula for calculating theoretical curves and direct resistivity interpretation is proposed for the case of arbitrary resistivity of the substratum. The basic equation for synthetic seismograms is reformulated in appendix A. In appendix B a method is indicated for the conversion of the seismic transfer function from arbitrary to perfectly reflective substratum.