Direct inversion of one-dimensional magnetotelluric data

1983 ◽  
Vol 88 (B3) ◽  
pp. 2407 ◽  
Author(s):  
Shimon Coen ◽  
Franchesca Quercia ◽  
Maria Mackiewicz
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.


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


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

2020 ◽  
Vol 222 (2) ◽  
pp. 1023-1033
Author(s):  
Max Moorkamp ◽  
Anna Avdeeva

SUMMARY We present a new approach that allows for the inversion of quantities derived from the observed data using non-diagonal data covariance matrices. For example, we can invert approximations of apparent resistivity and phase instead of magnetotelluric impedance using this methodology. Compared to the direct inversion of these derived quantities, the proposed methodology has two advantages: (i) If an inversion algorithm allows for the specification of a full data covariance matrix, users can invert for arbitrary derived quantities by specifying the appropriate covariance matrix instead of having to rely on the inversion code to have implemented this feature. (ii) It is fully compatible with the assumptions of least-squares optimization and thus avoids potential issues with bias when inverting quantities that are nonlinear functions of the original data, We discuss the theory of this approach and show an example using magnetotelluric data. However, the same method can be applied to other types of geophysical data, for example gravity gradient measurements.


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

Sign in / Sign up

Export Citation Format

Share Document