Three-dimensional magnetotelluric inversion using a heterogeneous smoothness-constrained least-squares method

Tateyuki Negi, Hideki Mizunaga, Koichi Asamori, Koji Umeda

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

This paper presents a fast algorithm for electromagnetic data inversion to three-dimensional (3D) resistivity models. The algorithm is distinctive for the level of accuracy it attains while bypassing the sensitivity matrix update. A common sensitivity matrix for homogeneous half-space is used in all iterations. Instead of updating the sensitivity matrix, the smoothness filter coefficients at each model element are updated, based on the spatial variations in resistivity in the model derived from the latest iteration. This substitution is expected not only to reduce the computation time required for large-scale inversions, such as those for 3D surveys, but also to allow the resolution of sharp boundaries in resistivity structures. Our algorithm was applied to 3D magnetotelluric inversion in order to confirm its effectiveness. Using synthetic examples under several conditions, we demonstrated that the method can reduce the number of forward calculations required to reduce data misfits to noise level, and that the method is robust for constructing target models even with sharp boundaries without generating fatally false resistivity structures or boundaries under noisy conditions.

Original languageEnglish
Pages (from-to)145-155
Number of pages11
JournalExploration Geophysics
Volume44
Issue number3
DOIs
Publication statusPublished - Sep 16 2013

Fingerprint

least squares method
electrical resistivity
inversions
matrix
iteration
sensitivity
matrices
data inversion
half space
half spaces
substitution
spatial variation
substitutes
electromagnetism
filter
filters
inversion
coefficients
method

All Science Journal Classification (ASJC) codes

  • Geophysics
  • Geology

Cite this

Three-dimensional magnetotelluric inversion using a heterogeneous smoothness-constrained least-squares method. / Negi, Tateyuki; Mizunaga, Hideki; Asamori, Koichi; Umeda, Koji.

In: Exploration Geophysics, Vol. 44, No. 3, 16.09.2013, p. 145-155.

Research output: Contribution to journalArticle

@article{3d4355a7761d4acea23c87f47caeca6a,
title = "Three-dimensional magnetotelluric inversion using a heterogeneous smoothness-constrained least-squares method",
abstract = "This paper presents a fast algorithm for electromagnetic data inversion to three-dimensional (3D) resistivity models. The algorithm is distinctive for the level of accuracy it attains while bypassing the sensitivity matrix update. A common sensitivity matrix for homogeneous half-space is used in all iterations. Instead of updating the sensitivity matrix, the smoothness filter coefficients at each model element are updated, based on the spatial variations in resistivity in the model derived from the latest iteration. This substitution is expected not only to reduce the computation time required for large-scale inversions, such as those for 3D surveys, but also to allow the resolution of sharp boundaries in resistivity structures. Our algorithm was applied to 3D magnetotelluric inversion in order to confirm its effectiveness. Using synthetic examples under several conditions, we demonstrated that the method can reduce the number of forward calculations required to reduce data misfits to noise level, and that the method is robust for constructing target models even with sharp boundaries without generating fatally false resistivity structures or boundaries under noisy conditions.",
author = "Tateyuki Negi and Hideki Mizunaga and Koichi Asamori and Koji Umeda",
year = "2013",
month = "9",
day = "16",
doi = "10.1071/EG13026",
language = "English",
volume = "44",
pages = "145--155",
journal = "Exploration Geophysics",
issn = "0071-3473",
publisher = "Consultants Bureau",
number = "3",

}

TY - JOUR

T1 - Three-dimensional magnetotelluric inversion using a heterogeneous smoothness-constrained least-squares method

AU - Negi, Tateyuki

AU - Mizunaga, Hideki

AU - Asamori, Koichi

AU - Umeda, Koji

PY - 2013/9/16

Y1 - 2013/9/16

N2 - This paper presents a fast algorithm for electromagnetic data inversion to three-dimensional (3D) resistivity models. The algorithm is distinctive for the level of accuracy it attains while bypassing the sensitivity matrix update. A common sensitivity matrix for homogeneous half-space is used in all iterations. Instead of updating the sensitivity matrix, the smoothness filter coefficients at each model element are updated, based on the spatial variations in resistivity in the model derived from the latest iteration. This substitution is expected not only to reduce the computation time required for large-scale inversions, such as those for 3D surveys, but also to allow the resolution of sharp boundaries in resistivity structures. Our algorithm was applied to 3D magnetotelluric inversion in order to confirm its effectiveness. Using synthetic examples under several conditions, we demonstrated that the method can reduce the number of forward calculations required to reduce data misfits to noise level, and that the method is robust for constructing target models even with sharp boundaries without generating fatally false resistivity structures or boundaries under noisy conditions.

AB - This paper presents a fast algorithm for electromagnetic data inversion to three-dimensional (3D) resistivity models. The algorithm is distinctive for the level of accuracy it attains while bypassing the sensitivity matrix update. A common sensitivity matrix for homogeneous half-space is used in all iterations. Instead of updating the sensitivity matrix, the smoothness filter coefficients at each model element are updated, based on the spatial variations in resistivity in the model derived from the latest iteration. This substitution is expected not only to reduce the computation time required for large-scale inversions, such as those for 3D surveys, but also to allow the resolution of sharp boundaries in resistivity structures. Our algorithm was applied to 3D magnetotelluric inversion in order to confirm its effectiveness. Using synthetic examples under several conditions, we demonstrated that the method can reduce the number of forward calculations required to reduce data misfits to noise level, and that the method is robust for constructing target models even with sharp boundaries without generating fatally false resistivity structures or boundaries under noisy conditions.

UR - http://www.scopus.com/inward/record.url?scp=84883737606&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84883737606&partnerID=8YFLogxK

U2 - 10.1071/EG13026

DO - 10.1071/EG13026

M3 - Article

AN - SCOPUS:84883737606

VL - 44

SP - 145

EP - 155

JO - Exploration Geophysics

JF - Exploration Geophysics

SN - 0071-3473

IS - 3

ER -