Magnetotelluric (MT) is a passive electromagnetic (EM) method which measure natural variations of electric and magnetic vector fields at the Earth surface to map subsurface electrical conductivity/resistivity structure. In this study, we obtained numerical solution of three-dimensional (3-D) MT using vector finite element method by solving second order Maxwell differential equation describing diffusion of plane wave through the conductive earth. Rather than the nodes of the element, the edges of the element is used as a vector basis to overcome the occurrence of nonphysical solutions that usually faced by scalar (node based) finite element method. Electric vector fields formulation was used and the resulting system of equation was solved using direct solution method to obtain the electric vector field distribution throughout the earth resistivity model structure. The resulting MT response functions was verified with 1-D layered Earth and 3-D2 COMMEMI outcropping structure. Good agreement is achieved for both structure models.

1.
Avdeev
,
D.
, “
Three-dimensional electromagnetic modelling and inversion from theory to application
”,
Surveys in Geophysics.
,
26
,
767
799
(
2005
).
2.
Börner
,
U. R.
, “
Numerical modelling in geo-electromagnetics: advances and challenges
”,
Surveys in Geophysics.
,
31
,
225
245
(
2010
).
3.
Cagniard
,
L.
, “
Basic theory of the magneto-telluric method of geophysical prospecting
”,
Geophysics
,
18
,
605
635
(
1953
).
4.
Farquharson
,
C. G.
, and
Miensopust
,
M. P.
, “
Three-dimensional finite-element modeling of magnetotelluric data with a divergence correction
”,
Journal of Applied Geophysics
,
75
,
699
710
(
2011
).
5.
Jin
,
J. M.
, “
The finite element method in electromagnetics
”,
John Wiley and Sons, Inc
. (
1993
).
6.
Mackie
,
R. L.
,
Madden
,
T. R.
, and
Wannamaker
P. E.
, 3-D “
Magnetotelluric modeling using difference equations – theory and comparisons to integral equation solutions
”,
Geophysics
,
58
,
215
226
(
1993
).
7.
Mogi
,
T.
, “
Three-dimensional modeling of magnetotelluric data using finite element method
”,
Journal of Applied Geophysics
,
35
,
185
189
(
1996
).
8.
Sam
,
C. T.
and
Hohmann
,
G. W.
, “
Integral equation modeling of tree-dimensional magnetotelluric response
”,
Geophysics
,
46
,
182
197
(
1981
).
9.
Schenk
,
O.
and
Gärtner
,
K.
, “
Solving unsymmetric sparse systems of linear equations with PARDISO
”,
Future Generation of Computer System
,
20
,
475
487
(
2004
).
10.
Siripunvaraporn
,
W.
,
Egbert
,
G.
, and
Lenbury
,
Y.
, “
Numerical accuracy of magnetotelluric modeling: A comparison of finite difference approximations
”,
Earth Planets Space
,
54
,
721
725
(
2002
).
11.
Wannamaker
,
P. E.
, “
Advances in three-dimensional magnetotelluric modeling using integral equations
”,
Geophysics
,
56
,
1716
1728
(
1991
).
12.
Zhdanov
,
M. S.
,
Varentsov
,
I. M.
,
Weaver
,
J. T.
,
Golubev
,
N. G.
, and
Krylov
,
V. A.
, “
Methods for modeling electromagnetic fields Results from COMMEMI –the international project on the comparison of modeling methods for electromagnetic induction
”,
Journal of Applied Geophysics
,
37
,
133
271
(
1997
).
This content is only available via PDF.
You do not currently have access to this content.