Finite element resistivity modelling for three-dimensional structures with arbitrary anisotropy

We present a three-dimensional finite element algorithm for direct current resistivity modelling. The standard Fortran code allows for nearly arbitrary conductivity structures including general anisotropy. The problem is formulated in terms of secondary potentials where mixed boundary conditions are...

Full description

Saved in:
Bibliographic Details
Published in:Physics of the earth and planetary interiors 2005-05, Vol.150 (1), p.15-27
Main Authors: Li, Yuguo, Spitzer, Klaus
Format: Article
Language:English
Subjects:
Electrical anisotropy
Finite elements
Singularity removal
Three-dimensional dc resistivity modelling
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We present a three-dimensional finite element algorithm for direct current resistivity modelling. The standard Fortran code allows for nearly arbitrary conductivity structures including general anisotropy. The problem is formulated in terms of secondary potentials where mixed boundary conditions are incorporated. Also in case of anisotropy, this type of boundary condition is superior to the Dirichlet type. We have verified the finite element method using an anisotropic two-layered earth, whose analytical solutions are available. For this simple model, the algorithm achieves high accuracy. The relative deviation between numerical and analytical solution is less than 1.2%. Due to the lack of further analytical references, the responses of three representative model types serve as a cross-check for plausibility and prove the operativeness of the code.
ISSN:0031-9201
1872-7395
DOI:10.1016/j.pepi.2004.08.014