Rapid least‐squares inversion of apparent resistivity pseudosections by a quasi‐Newton method1

Geophysical Prospecting - Tập 44 Số 1 - Trang 131-152 - 1996
M.H. Loke1, R. D. Barker2
1School of Physics, Universiti Sains Malaysia, 11800 Penang, Malaysia. Forrnerly at
2School of Earth Sciences, The University of Birmingham, Edgbaston, Birmingham B15 2TT, Kingdom.

Tóm tắt

Abstract

A fast inversion technique for the interpretation of data from resistivity tomography surveys has been developed for operation on a microcomputer. This technique is based on the smoothness‐constrained least‐squares method and it produces a two‐dimensional subsurface model from the apparent resistivity pseudosection. In the first iteration, a homogeneous earth model is used as the starting model for which the apparent resistivity partial derivative values can be calculated analytically. For subsequent iterations, a quasi‐Newton method is used to estimate the partial derivatives which reduces the computer time and memory space required by about eight and twelve times, respectively, compared to the conventional least‐squares method. Tests with a variety of computer models and data from field surveys show that this technique is insensitive to random noise and converges rapidly. This technique takes about one minute to invert a single data set on an 80486DX microcomputer.

Từ khóa


Tài liệu tham khảo

AndrewsR.1993.The impact of sewage sludge application of nitrate leaching from arable land on the unconfined chalk aquifer of East Anglia England.PhD thesis The University of Birmingham.

10.3997/1365-2397.1992004

10.1090/S0025-5718-1965-0198670-6

Broyden C.G., 1972, Numerical Methods for Unconstrained Optimization, 87

Burden R.L., 1981, Numerical Analysis

10.1190/1.1442813

Dennis J.E., 1983, Numerical Methods for Unconstrained Obtimization and Nonlinear Equations

10.1111/j.1365-2478.1979.tb00961.x

10.1190/1.1440762

George A., 1981, Computer Solution of Large Sparse Positive Definite Systems

Golub G.H., 1989, Matrix Computations

10.1016/0926-9851(93)90005-J

10.3997/1365-2397.1990008

HallofP.G.1957.On the interpretation of resistivity and induced polarization measurements.PhD thesis M.I.T.

Keller G.V., 1966, Electrical Methods in Geophysical Prospecting

Lawson C.L., 1974, Solving Last‐squares Problems

10.1111/j.1365-2478.1984.tb00726.x

LokeM.H.1994.The inversion of two‐dimensional apparent resistivity data.PhD thesis The University of Birmingham.

10.1190/1.1443900

10.1007/BF01436248

10.1111/j.1365-2478.1990.tb01859.x

10.2307/2005323

10.1190/1.1443128

Press W.H., 1988, Numerical Recipes in C

10.1190/1.1442649

10.1111/j.1365-2478.1992.tb00536.x

Scales L.E., 1985, Intoduction to Non‐linear Optimization, 10.1007/978-1-349-17741-7

Schwarz H.R., 1973, Numerical Analysis of Symmetric Matrices

Sheriff R.E., 1991, Encyclopedic Dictionary of Exploration Geophysics

10.1109/TGRS.1984.350575

10.1190/1.1441578

TurnbullJ.1936.A field and interpretation technique for resistivity surveying over two‐dimensional structures.PhD thesis. The University of Birmingham.

10.1016/B978-0-12-141050-6.50030-1