- Abstract
- Affiliation
- All
- Authors
- Book Series
- DOI
- EISBN
- EISSN
- Full Text
- GeoRef ID
- ISBN
- ISSN
- Issue
- Keyword (GeoRef Descriptor)
- Meeting Information
- Report #
- Title
- Volume

- Abstract
- Affiliation
- All
- Authors
- Book Series
- DOI
- EISBN
- EISSN
- Full Text
- GeoRef ID
- ISBN
- ISSN
- Issue
- Keyword (GeoRef Descriptor)
- Meeting Information
- Report #
- Title
- Volume

- Abstract
- Affiliation
- All
- Authors
- Book Series
- DOI
- EISBN
- EISSN
- Full Text
- GeoRef ID
- ISBN
- ISSN
- Issue
- Keyword (GeoRef Descriptor)
- Meeting Information
- Report #
- Title
- Volume

- Abstract
- Affiliation
- All
- Authors
- Book Series
- DOI
- EISBN
- EISSN
- Full Text
- GeoRef ID
- ISBN
- ISSN
- Issue
- Keyword (GeoRef Descriptor)
- Meeting Information
- Report #
- Title
- Volume

- Abstract
- Affiliation
- All
- Authors
- Book Series
- DOI
- EISBN
- EISSN
- Full Text
- GeoRef ID
- ISBN
- ISSN
- Issue
- Keyword (GeoRef Descriptor)
- Meeting Information
- Report #
- Title
- Volume

- Abstract
- Affiliation
- All
- Authors
- Book Series
- DOI
- EISBN
- EISSN
- Full Text
- GeoRef ID
- ISBN
- ISSN
- Issue
- Keyword (GeoRef Descriptor)
- Meeting Information
- Report #
- Title
- Volume

### NARROW

**Follow your search**

*Close Modal*

The front matter includes the title page, copyright page, contents, dedication, and preface.

Summary A simple exercise in inverting transient electromagnetic (TEM) data for the layering in a basin ran into difficulties when the inversion failed to converge for about half of the soundings. The failure was caused by the forward modeling program that created artifacts for models with thin, shallow conductive layers. The artifacts were difficult to identify visually on raw plots because of the large dynamic range of most TEM data. They did show up clearly in relative comparisons of the curves from two programs. This experience led to a broader study of the differences between the many layered-earth and 3-D TEM modeling codes available and of the various approximations used. This study has already resulted in improvements in several codes, which generated more confidence in their results.

#### Transient Diffusive Electromagnetic Field Computation—A Structured Approach Based on Reciprocity

Summary The reciprocity theorem for transient diffusive electromagnetic fields is taken as the point of departure for developing computational methods to model such fields. Mathematically, the theorem is representative of any weak formulation of the field problem. Physically, the theorem describes the interaction between (a discretized version of) the actual field and a suitably chosen computational state. The choice of the computational state determines which type of computational method results from the analysis. It is shown that the finite-element method, the integral-equation method, and the domain-integration method can be viewed as particular cases of discretization of the reciprocity relation. The local field representations of the electric- and the magnetic-field strengths in terms of edge-element expansion functions are worked out in some detail. The emphasis is on time-domain methods. The relationship with complex frequency-domain methods is indicated and used to symmetrize the basic field equations. This symmetrization expresses the correspondence that exists between transient electromagnetic wavefields in lossless media and transient diffusive electromagnetic fields in conductive media where the electric displacement-current contribution to the field can be neglected in the time window of observation. This aspect is also of importance in numerical modeling.

Summary Standard 3-D integral-equation methods using pulse expansion functions and point matching have difficulties with large contrasts in conductivity between the body and the host. We formulate a new expansion with rectangular building blocks and trilinear functions that more accurately models general current distributions inside the scattering body without increasing the number of unknowns. The formulation is in the frequency domain; to transform results to the time domain, we solve for the coefficients of an expansion of the decaying field in transients that are determined analytically by the behavior of an iterative solution of the integral equation. The method requires a modest number of frequencies and is stable and accurate even at late times.

Summary Numerical methods of electromagnetic (EM) modeling usually represent spatial and temporal variations of the field with piecewise polynomials of low order. We study the accuracy of several piecewise representations—constant, linear, quadratic, and cubic—using an integral equation to compute EM scattering from a thin wire. The actual scattered field of a uniform conductor should be holomorphic (i.e., have derivatives of all orders). Although holomorphism is impossible to achieve with piecewise polynomial expansions, we do ensure that the scattered electric field is continuous along the uniform conductor and has derivatives continuous to the order allowed by the polynomial basis. The higher-order bases dramatically improve the accuracy to which boundary conditions are satisfied. Higher-order expansions also reduce the degree to which the accuracy depends on the method used to solve the numerical equations. Moreover, convergence can be achieved with fewer discretization cells. Larger cell sizes, in turn, allow more accurate numerical evaluations of the integrals over the (singular) Green’s functions. We believe that these higher-order bases offer substantial improvements in 3-D numerical modeling of EM fields.

Summary. We describe a new method for solving Maxwell’s equations that couples scalar surface integral equations (SIEs) with an iterative implementation of the (vector) boundary conditions. The formulation can handle multiple homogeneous regions and, in computational efficiency, compares favorably with the more conventional volume integral method. The system equations are solved iteratively in three stages: (1) iterations to solve the coupled scalar SIE; (2) iterations to implement vector boundary conditions; and (3) iterations to calculate surface fields in three dimensions. Numerical results for a sphere excited by a plane wave and for an oblate spheroid in a half-space excited by a rectangular surface loop show that the method is stable for a wide range of conductivity contrasts and gives reasonable results for the electromagnetic response from the resistive up to the inductive limit.

Summary. Conventional methods for solving 3-D volume integral equations in electromagnetic (EM) modeling produce matrices that are ill-conditioned when the conductivity contrast between the target and the host is large or when the host resistivity itself is small. This problem is especially acute when pulse basis functions represent the scattering currents because the artificial discontinuities in the current create spurious charges in homogeneous conductive regions (where there are no physical charges). Basis functions with higher-order continuity do not entirely eliminate this problem. The field of these charges is amplified at high host resistivities by Green’s tensor and artificially attenuates the physical vortex currents. We present a new formulation that eliminates these spurious charges analytically by replacing the volume integral for the field caused by charges (on the boundary of homogeneous regions) with a surface integral, while retaining the volume integral for the field of the induced currents. This formulation is mathematically closer to the physics of EM induction and effectively removes a major source of error. Also, condition numbers of matrix systems without artificial charges are also much smaller than those of the conventional method.

#### Consistent Discretization of Electromagnetic Fields and Transient Modeling

Summary Differential forms provide an elegant formulation of electromagnetic field theory. Their geometric structure also leads to a self-consistent scheme for discretizing Maxwell’s equations in conducting media. This discretization, which is essentially a staggered-grid scheme, preserves differential operator identities, conservation laws, and physical boundary conditions. The self-consistent scheme has been implemented in a code to model transient electromagnetic fields in a half-space with a step-function current excitation. For small conductivity contrasts results from the self-consistent scheme agree very well with results from a 3-D integral-equation code and from a staggered-grid finite-difference code. For high contrasts, though, the results from the integral-equation code differ. The difference may be a result of improper discretization of the integral equation.

Summary Electrical anisotropy in the Earth is usually a scale effect, created by averaging over structures with a preferred orientation. Homogenization theory applied to typical geological structures can help to determine the range of values to be expected in electrically anisotropic formations. The staggered-grid finite-difference method for Maxwell’s equations can easily accommodate anisotropic regions, without a significant increase in computational load. A simple magnetotelluric model shows that, in elongated anisotropic structures, the tipper is no longer perpendicular to the strike direction; thus, the electric and magnetic polarizations are mixed.

Summary The standard staggered grid for Maxwell’s equations is awkward for anisotropic media because the different components of the electric and magnetic fields are located at different nodes. There is, however, a natural alternative that places all components of the electric field at each node of one grid and all components of the magnetic field at each node of a staggered grid. This staggering allows a conservative finite-difference approximation for Maxwell’s equations with arbitrary 3-D tensor electrical conductivity, magnetic permeability, and dielectric permittivity. An example of the time-domain solution using spectral Lanczos decomposition is considered.

#### Finite-Difference Modeling of 3-D EM Fields with Scalar and Vector Potentials

Summary Scalar and vector potentials provide a convenient formulation for finite-difference modeling of harmonic electromagnetic (EM) fields in 3-D media. Because the EM potentials are continuous everywhere, the finite-difference method can handle large contrasts in conductivity or dielectric constant. In the Coulomb gauge, the method is stable at low frequencies because the potentials approach those for a static (dc) formulation. The linear system of finite-difference equations can be solved iteratively with a biconjugate-gradient method that is diagonally preconditioned. An impedance boundary condition terminates the grid. Modeling of a crosswell EM survey with this method shows that 2-m electric dipoles radiating at 5 MHz can detect a small zone of dense nonaqueous-phase liquid (DNAPL) with high resistivity and low dielectric constant between wells 10 m apart. Because the responses to contrasts in dielectric constant and resistivity are distinctly different, both quantities could be estimated. The response changes significantly with the position of the DNAPL zone, the weakest response occurring when the zone is just below the water table.

Summary Accuracy in finite-difference (FD) modeling is closely related to the discretization scheme, whereas speed depends mainly on the equation solvers. We compare the accuracy of five FD discretization schemes for 3-D resistivity modeling. Three schemes yield good results: a method using volume-weighted averages from conductivities assigned to neighboring grid cells, a method that integrates over elemental volumes, and a resistivity network approach. Discretization by elemental volume leads to coupling coefficients that are similar to those derived from the volume-weight method. The coefficients only differ by a real factor. In the second section, the cumulative amount of numerical work as a measure of speed is compared for five different equation solvers with and without preconditioning. The most efficient equation solver for symmetric matrices is the preconditioned conjugate gradient method. General matrix solution methods for both symmetric and nonsymmetric matrices—such as ORTHOMIN and the methods of stabilized biconjugate gradients and squared conjugate gradients—also achieve satisfactory convergence rates.

Summary The electromagnetic (EM) inverse problem in geophysics consists of deducing the Earth’s conductivity from measurements of electric or magnetic fields near the surface. It requires fast and accurate forward modeling, a method of solving nonlinear equations, and a method of controlling nonuniqueness. Four modeling methods are considered for simulation of airborne EM data over a 3-D conductivity model: a one-dimensional (1-D) approximation, the 3-D Born approximation, a finite-element (FE) method with approximate boundary conditions, and an FE method with exact boundary conditions. The 1-D and Born approximations are fast, but not very accurate. The two FE methods give similar responses, and the method with approximate boundary conditions requires much less computer resources. Three methods are considered for solving the EM inverse problem by unconstrained nonlinear optimization: conjugate gradient (CG), quasi-Newton, and Gauss-Newton (GN) methods. For the simple models considered, the GN method required the least CPU time; the CG method, the most. An alternative to regularization for controlling nonuniqueness is joint inversion of several different surveys to produce a single model that fits all of the data. An example with airborne EM data and pole-pole dc resistivity data shows that such joint inversion can effectively reduce nonuniqueness in EM inversion.

Summary We propose a new algorithm for 3-D electromagnetic inversion that works with the magnetic-field integral equation (instead of the traditional electric-field integral equation). The forward and inverse integral equations are discretized by the finite element method; solution of the matrix system involves alternating conjugate-gradient and biconjugate-gradient iterations. The solution is regularized by a novel external annealing technique. A coupled-domain decomposition allows a very efficient (96%) implementation of the algorithm on massively parallel machines. Tests on both synthetic and field data from environmental sites yield good subsurface images in a reasonable amount of computational time.

Summary We develop an inverse scattering algorithm for imaging 3-D dielectric scatterers with microwaves. The algorithm uses the biconjugate gradient and fast Fourier transform method (BCG-FFT) for forward modeling and a conjugate gradient algorithm to solve the inverse problem, which is formulated as a nonlinear optimization problem. The BCG-FFT method reduces the computational complexity of forward modeling to order N log 2 N , where N is the number of unknowns used to represent the vector fields in 3-D space. The conjugate gradient method requires the computation of the Fréchet derivative (matrix of partial derivatives) operating on a vector of residuals, which is equivalent to solving another forward problem. Simulations show the efficiency of this algorithm, especially when the number of illumination angles is limited.

Summary Parameter estimation in geoelectromagnetics aims to obtain the most important parameters of a well-defined conductivity model of the Earth. These parameters are features of typical geological structures, such as depth and size of conductive or resistive targets, angle of dike inclination and its length, and conductivity of anomalous bodies. We develop this approach through regularized nonlinear optimization. We use finite differences of forward computations and Broyden’s updating formula to compute sensitivities (Frechet or partial derivatives) for each parameter. To estimate the optimal step length, we apply line search, with a simple and fast parabolic correction. Our inversion also includes Tikhonov’s regularization procedure. We use our method to study measurements of the magnetic fields from a conductive body excited by a loop source at the surface. Keeping the depth of the body constant, we estimate the horizontal coordinates of the body from three components of the magnetic field measured in a borehole. These measurements accurately determine the direction to the conductive target.

Summary The quasi-linear (QL) approximation replaces the (unknown) total field in the integral equation of electromagnetic (EM) scattering with a linear transformation of the primary field. This transformation involves the product of the primary field with a reflectivity tensor, which is assumed to vary slowly inside inhomogeneous regions and therefore can be determined numerically on a coarse grid by a simple optimization. The QL approximation predicts EM responses accurately over a wide range of frequencies for conductivity contrasts of more than 100 to 1 between the scatterer and the background medium. It also provides a fast-forward model for 3-D EM inversion. The inversion equation is linear with respect to a modified material property tensor, which is the product of the reflectivity tensor and the anomalous conductivity. We call the (regularized) solution of this equation a quasi-Born inversion . The material property tensor (obtained by inversion of the data) then is used to estimate the reflectivity tensor inside the inhomogeneous region and, in turn, the anomalous conductivity. Solution of the nonlinear inverse problem thus proceeds through a set of linear equations. In practice, we accomplish this inversion through gradient minimization of a cost function that measures the error in the equations and includes a regularization term. We use synthetic experiments with plane-wave and controlled sources to demonstrate the accuracy and speed of the method.

Summary Linearized, iterative inversion of electromagnetic data requires computation of partial derivatives (sensitivities) with respect to parameters of the model, e.g., the conductivities of cells. Methods based on the adjoint equation compute these sensitivities by integrating, over each cell, the scalar product of an adjoint electric field with the electric field produced by the forward modeling in the preceeding iteration. We approximate the partial derivative by computing the adjoint field in either a homogeneous or a layered half-space. Computation of the approximate adjoint field is significantly faster than that of the true adjoint field, with the relative efficiency increasing with the size of the problem. Approximate sensitivities compare well with exact values for sample controlled-source surveys in 2-D and 3-D models. We demonstrate that approximate sensitivities can be sufficiently accurate to drive an iterative algorithm by inverting synthetic magnetotelluric data. Approximate sensitivities should enable the solution of inverse problems larger than those now practical.

#### Linearized 3-D Electromagnetic Vector Wave Inversion

Summary Generalized (scalar) diffraction tomography is a linear inverse scattering method that can be extended to electromagnetic vector fields with complete polarization information. Its essential equation is a vector form of the Porter–Bojarski integral equation linearized in the material parameters through the Kirchhoff approximation to scattering. This vector equation can be inverted with dyadic algebra and the standard techniques of diffraction tomography using data from multiple frequencies (frequency diversity) or multiple angles of incidence (angle diversity). An algorithm using scattering at multiple frequencies to reconstruct perfectly conducting objects is discussed in detail and checked against synthetic data generated with the MAFIA code for an airplane model. The results are overwhelmingly superior to those obtained by scalar inversion.

Summary Theoretical inverse problem is the terminology (in the Russian literature) for a geophysical inverse problem in which the field is given by an explicit expression. This can arise, for example, when data are approximated by singular sources (monopoles, dipoles, etc.) in a half-space. This chapter derives explicit integrodifferential equations of theoretical inverse problems for 3-D electromagnetic fields satisfying the Helmholtz, telegraphic, and diffusion equations. Functional equations for a perfect conductor are presented, and some numerical examples are given.