Inversion of Magnetotelluric Data for a One-Dimensional Conductivity

Author(s):  
Kenneth P. Whittall ◽  
Douglas W. Oldenburg
Geophysics ◽  
1976 ◽  
Vol 41 (4) ◽  
pp. 766-770 ◽  
Author(s):  
F. E. M. Lilley

Observed magnetotelluric data are often transformed to the frequency domain and expressed as the relationship [Formula: see text]where [Formula: see text] [Formula: see text] and [Formula: see text] [Formula: see text] represent electric and magnetic components measured along two orthogonal axes (in this paper, for simplicity, to be north and east, respectively). The elements [Formula: see text] comprise the magnetotelluric impedance tensor, and they are generally complex due to phase differences between the electric and magnetic fields. All quantities in equation (1) are frequency dependent. For the special case of “two‐dimensional” geology (where structure can be described as having a certain strike direction along which it does not vary), [Formula: see text] with [Formula: see text]. For the special case of “one‐dimensional” geology (where structure varies with depth only, as if horizontally layered), [Formula: see text] and [Formula: see text].


1983 ◽  
Vol 88 (B3) ◽  
pp. 2407 ◽  
Author(s):  
Shimon Coen ◽  
Franchesca Quercia ◽  
Maria Mackiewicz

Geophysics ◽  
1977 ◽  
Vol 42 (5) ◽  
pp. 995-1005 ◽  
Author(s):  
W. R. Petrick ◽  
W. H. Pelton ◽  
S. H. Ward

Ridge regression inversion has been used to test the applicability of various one‐dimensional crustal models to the interpretation of deep Schlumberger sounding data from southern Africa (Van Zijl and Joubert, 1975). Four main models were investigated: a simple three‐layered earth, a layered earth with a transition zone exhibiting a linear decrease in log resistivity with depth, a similar earth with the transition zone determined by cubic splines, and a model having exponential resistivity behavior at depth. The last model corresponds to temperature‐dependent semiconduction through solid mineral grains (Brace, 1971). It was found that all of these models are capable of fitting the sounding data from southwestern Africa, while all except the semiconduction model fit the data from southeastern Africa. One is, thereby, immediately alerted to the problem of lack of resolution in Schlumberger sounding data where geologic control is not available. A major with the inversion of Schlumberger data alone is that accurate information is obtainable only for the resistivity‐thickness product of the resistive portion of the crust. On the other hand, magnetotelluric data, when available, tends to provide information on the thickness, but very little information on the true resistivity of the section. In order to resolve both resistivity and thickness it is possible to invert simultaneously Schlumberger and magnetotelluric (MT) data. Results obtained from the combined inversion of the African resistivity data and hypothetical MT data show that a considerable improvement in model resolution can be achieved using MT amplitude data even of poor accuracy from a relatively limited frequency range (0.1 to 100 Hz), whereas inclusion of MT phase information is of negligible additional benefit. Unfortunately, no significant test can be made, from data available at the time of our analysis, of the applicability of one‐dimensional inversion in a geologic circumstance which probably demands more dimension.


2018 ◽  
Vol 70 (1) ◽  
Author(s):  
Dennis Conway ◽  
Graham Heinson ◽  
Nigel Rees ◽  
Joseph Rugari

Geophysics ◽  
1986 ◽  
Vol 51 (8) ◽  
pp. 1603-1607 ◽  
Author(s):  
Kenneth P. Whittall

I present an algorithm for the one‐dimensional magnetotelluric inverse problem of finding conductivity as a function of depth in the earth. The algorithm uses linear programming to solve an integral form of a nonlinear Riccati equation. This iterative scheme sacrifices the efficiency of direct inversion for the overwhelming advantages of incorporating localized conductivity constraints. I use localized conductivity constraints in two ways to combat the nonuniqueness of the nonlinear inverse problem. First, I impose physical constraints derived from external sources to restrict the nonuniqueness and construct conductivity models that are closer to reality. Second, I impose arbitrary constraints in an effort to assess the extent of nonuniqueness and explore the range of acceptable profiles. The first technique enhances the reliability of an interpretation, and the second measures the plausibility of particular conductivity features.


2002 ◽  
Vol 359 (3-4) ◽  
pp. 371-380
Author(s):  
Kiyoshi Fuji-ta ◽  
Satoshi Fujiwara ◽  
Masahiro Ichiki ◽  
Yutaka Makino

1989 ◽  
Vol 20 (2) ◽  
pp. 335
Author(s):  
J.P. Cull ◽  
J.D. Gray

Magnetotelluric data obtained in the Eromanga Basin can be interpreted using one-dimensional models to describe plane layers consistent with geological mapping. Interpretations are based on the results of non-linear inversions generating a minimum least-squares error between the observations and the model. However there is no statistical justification for selecting highly complex starting models. In particular adequate solutions can be generated using models based on 2, 3 or 4 layers over basement; additional layers defining fine structure can only be retained using external geological constraints. Solutions based on random starting models suggest gradations in resistivity (1?4 ohm m) consistent with sediment compaction. Major discontinuities in all models (4?30 ohm m) are assumed to indicate a basement contact at depths of 7?8 km.


Sign in / Sign up

Export Citation Format

Share Document