A MATLAB PROGRAM FOR THE ANALYSIS AND INTERPRETATION OF TRANSIENT ELECTROMAGNETIC SOUNDING DATA


Published: Jan 25, 2010
Keywords:
Electromagnetic sounding TEM MATLAB geophysical interpretation
Andreas Tzanis
Abstract
Herein we present a software system, written in MATLAB, to interpret TEM sounding data. The program, dubbed maTEM, is designed to process, model and invert multiple soundings, either individually, or simultaneously along profiles. The latter capability allows for laterally constrained inversion, so as to generate pseudo-2D or 2D resistivity sections based on the program EM1DINV v2.13 by the Hydro-Geophysics Group of the University of Aarhus, Denmark. Using maTEM, the analyst may import and display data multiple data sets, denoise and smooth the data, perform approximate inversions, design 1-D model(s) graphically, perform forward modelling and inversion and generate/ update data base in which to store the results. Finally, the analyst may use the data base to create 2-D and 3-D displays of the geoelectric structure with built-in graphical functions. maTEM is highly modular so that additional functions can be added at any time, at minimal programming cost. Although the software presented herein is focused on the analysis of TEM data, the maTEM concept has been designed ready to incorporate additional electrical and EM geophysical sounding methods and to mutually constrained analysis of different geophysical data sets.
Article Details
  • Section
  • Geophysics
Downloads
Download data is not yet available.
References
Anderson, W.L., 1982a. Calculation of transient soundings for a coincident loop system (Program
TCOLOOP). USGS Open-File Report 82-378.
Anderson, W.L., 1982b. Nonlinear least-squares inversion of transient soundings for a coincident loop system
(Program NLSTCO), USGS Open-File Report 82-1064
Anderson, W.L., 1982c. Nonlinear least-squares inversion of transient soundings for a central induction
loop system (Program NLSTCI), USGS Open-File Report 82-1129
Auken E. and Christiansen A.V. 2004. Layered and laterally constrained 2D inversion of resistivity data.
Geophysics 69, 752–761.
XLIII, No 4 – 1951
Auken, E., Christiansen, A.V., Jacobsen, B.H., Foged, N. and Sørensen, K.I., 2005. Piecewise 1D laterally
constrained inversion of resistivity data, Geophysical Prospecting, 53, 497–506.
Barsukov, P.O., Fainberg, E.B. and Khabensky, E.O., 2007. Shallow investigations by TEM-FAST technique:
Methodology and examples, in V.V. Spichak (editor), Methods in Geochemistry and Geophysics,
Volume 40, 56-77.
Christiansen, A. V. and Auken, E., 2004: Optimizing a layered and laterally constrained 2D inversion of
resistivity data using Broyden’s update and 1D derivatives. Journal of Applied Geophysics, 56, 247-
Ingeman-Nielsen, T. and Baumgartner, F. 2006. CR1Dmod: A Matlab program to model 1D complex resistivity
effects in electrical and electromagnetic surveys. Computers & Geosciences, 32, 1411–1419
Karmis, P., 2003. Automated interpretation of transient electromagnetic soundings. Doctoral Dissertation,
N.K University of Athens (in Greek).
Tartaras, E., Zhdanov, M.S. Wada, K. Saito, A. and Hara, T., 2000. Fast Imaging of TDEM data based on
S-inversion. Journal of Applied Geophysics, 43, 15–32.
Wahba, G., 1990. Spline models for observational data. In: CBMS-NSF Regional Conference Series in Applied
Mathematics 59, Soc. Industrial and Applied Mathematics, Philadelphia, Pennsylvania (1990),
p. 169.