New 3D parallel GILD electromagnetic modeling and nonlinear inversion using global magnetic integral and local differential equation
This paper describes a new 3D parallel GILD electromagnetic (EM) modeling and nonlinear inversion algorithm. The algorithm consists of: (a) a new magnetic integral equation instead of the electric integral equation to solve the electromagnetic forward modeling and inverse problem; (b) a collocation finite element method for solving the magnetic integral and a Galerkin finite element method for the magnetic differential equations; (c) a nonlinear regularizing optimization method to make the inversion stable and of high resolution; and (d) a new parallel 3D modeling and inversion using a global integral and local differential domain decomposition technique (GILD). The new 3D nonlinear electromagnetic inversion has been tested with synthetic data and field data. The authors obtained very good imaging for the synthetic data and reasonable subsurface EM imaging for the field data. The parallel algorithm has high parallel efficiency over 90% and can be a parallel solver for elliptic, parabolic, and hyperbolic modeling and inversion. The parallel GILD algorithm can be extended to develop a high resolution and large scale seismic and hydrology modeling and inversion in the massively parallel computer.
- Research Organization:
- Lawrence Berkeley National Lab., Earth Sciences Div., CA (United States)
- Sponsoring Organization:
- USDOE Office of Energy Research, Washington, DC (United States)
- DOE Contract Number:
- AC03-76SF00098
- OSTI ID:
- 334369
- Report Number(s):
- LBNL-42105; CONF-971111-; ON: DE99002422; TRN: AHC29914%%173
- Resource Relation:
- Conference: 67. annual meeting and expo of the Society of Exploration Geophysicists (SEG), Dallas, TX (United States), 2-7 Nov 1997; Other Information: PBD: Jul 1998
- Country of Publication:
- United States
- Language:
- English
Similar Records
A new parallel 3-D numerical modeling of the electromagnetic fields
Three-dimensional induced polarization data inversion for complex resistivity