Skip to Main Content
Book Chapter

PART IV: 3-D EM and Parallel Computers

Published:
January 01, 1999

Summary

A numerical method has been developed to simulate the electromagnetic (EM) response of a 3-D earth to a dipole source at frequencies ranging from 100 Hz to 100 MHz. The problem is formulated in the frequency domain with a modified vector Helmholtz equation for the scattered electric fields. The differential equation is approximated on a staggered finite-difference grid, giving a sparse complex symmetric matrix equation. The system is solved by a preconditioned quasi-minimum–residual method.

Dirichlet boundary conditions are imposed at the edges of the mesh by setting the tangential electric fields equal to zero. At frequencies less than 1 MHz, grid stretching reduces reflections off the grid boundaries. At higher frequencies, absorbing boundary conditions are imposed by making the stretching parameters of the modified vector Helmholtz equation complex to introduce loss at the boundaries.

Iterative solution to the nonlinear 3-D EM inverse problem proceeds by linearized model updates using the method of conjugate gradients. Full wave equation modeling is employed to compute model sensitivities and predicted data in the frequency domain with the 3-D finite-difference algorithm. Both the forward and inverse solutions are implemented on a massively parallel computing platform which allows forward and inverse models with millions and ten of thousands of parameters, respectively.

You do not currently have access to this article.

Figures & Tables

Contents

Geophysical Developments Series

Three-Dimensional Electromagnetics

Society of Exploration Geophysicists
Volume
7
ISBN electronic:
9781560802154
Publication date:
January 01, 1999

GeoRef

References

Related

Citing Books via

Close Modal
This Feature Is Available To Subscribers Only

Sign In or Create an Account

Close Modal
Close Modal