Premium
A MAN/COMPUTER INTERPRETATION SYSTEM FOR RESISTIVITY SOUNDINGS OVER A HORIZONTALLY STRAFIFIED EARTH *
Author(s) -
JOHANSEN H. K.
Publication year - 1977
Publication title -
geophysical prospecting
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.735
H-Index - 79
eISSN - 1365-2478
pISSN - 0016-8025
DOI - 10.1111/j.1365-2478.1977.tb01196.x
Subject(s) - singular value decomposition , logarithm , eigenvalues and eigenvectors , least squares function approximation , sequence (biology) , mathematics , minification , interpretation (philosophy) , ellipsoid , value (mathematics) , decomposition , geology , algorithm , mathematical analysis , computer science , mathematical optimization , statistics , chemistry , geodesy , physics , biochemistry , quantum mechanics , estimator , programming language , organic chemistry
A bstract The proposed system works as follows:1 By a trial‐and‐error procedure using a graphic display terminal a geologically relevant layer sequence with parameters (ρ j , d j ) is adjusted to yield roughly the measured curve. 2 The resulting layer sequence is used as starting model for an iterative least squares procedure with singular value decomposition. Minimization of the sum of the squares of the logarithmic differences between measured and calculated values with respect to the logarithms of the resistivities and thicknesses as parameters linearizes the problem to a great extent, with two important implications:a) a considerable increase in speed (the number of iterations goes down), thus making it cheap to achieve the optimum solution; b) the confidence surfaces in parameter space are well approximated by the hyper‐ellipsoids defined by the eigenvalues and eigenvectors of the normal equations.Since these are known from the singular value decomposition we do in fact know all possible solutions compatible with the measured curve and the geological concept.3 It is possible to “freeze” any combination of parameters at predetermined values. Thus extra knowledge and/or hypotheses are easily incorporated and can be tested by rerunning step (2). The overall computing time for a practical case is of the order of 10 sec on a CDC 6400.