Three‐dimensional wideband electromagnetic modeling on massively parallel computers

Abstract

A method is presented for modeling the wideband, frequency domain electromagnetic (EM) response of a three‐dimensional (3‐D) earth to dipole sources operating at frequencies where EM diffusion dominates the response (less than 100 kHz) up into the range where propagation dominates (greater than 10 MHz). The scheme employs the modified form of the vector Helmholtz equation for the scattered electric fields to model variations in electrical conductivity, dielectric permitivity and magnetic permeability. The use of the modified form of the Helmholtz equation allows for perfectly matched layer ( PML) absorbing boundary conditions to be employed through the use of complex grid stretching. Applying the finite difference operator to the modified Helmholtz equation produces a linear system of equations for which the matrix is sparse and complex symmetrical. The solution is obtained using either the biconjugate gradient (BICG) or quasi‐minimum residual (QMR) methods with preconditioning; in general we employ the QMR method with Jacobi scaling preconditioning due to stability. In order to simulate larger, more realistic models than has been previously possible, the scheme has been modified to run on massively parallel (MP) computer architectures. Execution on the 1840‐processor Intel Paragon has indicated a maximum model size of 280 × 260 × 200 cells with a maximum flop rate of 14.7 Gflops. Three different geologic models are simulated to demonstrate the use of the code for frequencies ranging from 100 Hz to 30 MHz and for different source types and polarizations. The simulations show that the scheme is correctly able to model the air‐earth interface and the jump in the electric and magnetic fields normal to discontinuities. For frequencies greater than 10 MHz, complex grid stretching must be employed to incorporate absorbing boundaries while below this normal (real) grid stretching can be employed.

A hybrid finite-difference and integral-equation method for modeling and inversion of marine controlled-source electromagnetic data

“…The implementation of the FD method developed in this paper follows that of Newman and Alumbaugh (1995) and Alumbaugh et al (1996). The method solves Maxwell's equations in the frequency domain based on an FD scheme on a staggered grid and uses the anomalous field formulation with the total field being decomposed into background E b and anomalous E a fields.…”

“…The conductivities where the electric fields are located are represented by a weighted average of conductivities of the four adjoining cells based on Ampere's law (Wang and Hohmann, 1993;Alumbaugh et al, 1996).…”

A review of block Krylov subspace methods for multisource electromagnetic modelling

“…Large-scale, massively parallel schemes typically use a domain decomposition paradigm and parallelize solver (Alumbaugh et al 1996). The achieved speedup slowly degrades with increasing number of cores, because of the overhead due to message passing.…”

Electromagnetic inversion for environmental site characterization: Data quality versus image resolution

“…The roughness or regularization matrix (W) consists of a fmite difference approxiination to the Laplacian operator, the observed data are represented by the vector d, and the predicted data arising from the model m 1 at the i'th iteration of the inversion procedure are denoted by These values are calculated using the fmite differences scheme given in Alumbaugh et al (1996). The data weighting matrix D is diagonal and consists of individual estimates of the data noise.…”

“…To demonstrate the problems caused by noisy in-phase data, a simulation of the CTP survey was calculated for two O.1S/m blocks (9m wide by Sm tail) separated by 9m and buried at 3m depth in a 0.OLS/m half-space using the scheme described in Alumbaugh et al (1996). A data sampling interval of 3m was used with a source-receiver separation of lOm, and following the specifications of the MaxMin system eight frequencies were employed.…”