skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: New advances in three-dimensional controlled-sourceelectromagnetic inversion

Abstract

New techniques for improving both the computational andimaging performance of the three dimensional (3D) electromagnetic inverseproblem are presented. A non-linear conjugate gradient algorithm is theframework of the inversion scheme. Full wave equation modelling forcontrolled sources is utilized for data simulation along with anefficient gradient computation approach for the model update. Improvingthe modelling efficiency of the 3D finite difference method involves theseparation of the potentially large modelling mesh, defining the set ofmodel parameters, from the computational finite difference meshes usedfor field simulation. Grid spacings and thus overall grid sizes can bereduced and optimized according to source frequencies and source-receiveroffsets of a given input data set. Further computational efficiency isobtained by combining different levels of parallelization. While theparallel scheme allows for an arbitrarily large number of parallel tasks,the relative amount of message passing is kept constant. Imageenhancement is achieved by model parameter transformation functions,which enforce bounded conductivity parameters and thus prevent parameterovershoots. Further, a remedy for treating distorted data within theinversion process is presented. Data distortions simulated here includepositioning errors and a highly conductive overburden, hiding the desiredtarget signal. The methods are demonstrated using both synthetic andfield data.

Authors:
;
Publication Date:
Research Org.:
Ernest Orlando Lawrence Berkeley NationalLaboratory, Berkeley, CA (US)
Sponsoring Org.:
USDOE Director. Office of Science. Basic EnergySciences
OSTI Identifier:
925536
Report Number(s):
LBNL-63010
Journal ID: ISSN 0956-540X; R&D Project: 468113; BnR: KC0303010; TRN: US200809%%793
DOE Contract Number:
DE-AC02-05CH11231
Resource Type:
Journal Article
Resource Relation:
Journal Name: Geophysical Journal International; Journal Volume: 172; Journal Issue: 2; Related Information: Journal Publication Date: 02/2008
Country of Publication:
United States
Language:
English
Subject:
54; ALGORITHMS; EFFICIENCY; FINITE DIFFERENCE METHOD; OVERBURDEN; PERFORMANCE; POSITIONING; SIMULATION; TARGETS; TRANSFORMATIONS; WAVE EQUATIONS

Citation Formats

Commer, Michael, and Newman, Gregory A. New advances in three-dimensional controlled-sourceelectromagnetic inversion. United States: N. p., 2007. Web.
Commer, Michael, & Newman, Gregory A. New advances in three-dimensional controlled-sourceelectromagnetic inversion. United States.
Commer, Michael, and Newman, Gregory A. Sat . "New advances in three-dimensional controlled-sourceelectromagnetic inversion". United States. doi:. https://www.osti.gov/servlets/purl/925536.
@article{osti_925536,
title = {New advances in three-dimensional controlled-sourceelectromagnetic inversion},
author = {Commer, Michael and Newman, Gregory A.},
abstractNote = {New techniques for improving both the computational andimaging performance of the three dimensional (3D) electromagnetic inverseproblem are presented. A non-linear conjugate gradient algorithm is theframework of the inversion scheme. Full wave equation modelling forcontrolled sources is utilized for data simulation along with anefficient gradient computation approach for the model update. Improvingthe modelling efficiency of the 3D finite difference method involves theseparation of the potentially large modelling mesh, defining the set ofmodel parameters, from the computational finite difference meshes usedfor field simulation. Grid spacings and thus overall grid sizes can bereduced and optimized according to source frequencies and source-receiveroffsets of a given input data set. Further computational efficiency isobtained by combining different levels of parallelization. While theparallel scheme allows for an arbitrarily large number of parallel tasks,the relative amount of message passing is kept constant. Imageenhancement is achieved by model parameter transformation functions,which enforce bounded conductivity parameters and thus prevent parameterovershoots. Further, a remedy for treating distorted data within theinversion process is presented. Data distortions simulated here includepositioning errors and a highly conductive overburden, hiding the desiredtarget signal. The methods are demonstrated using both synthetic andfield data.},
doi = {},
journal = {Geophysical Journal International},
number = 2,
volume = 172,
place = {United States},
year = {Sat May 19 00:00:00 EDT 2007},
month = {Sat May 19 00:00:00 EDT 2007}
}
  • Inversion of transient electromagnetic (TEM) data sets to image the subsurface three-dimensional (3-D) electrical conductivity and magnetic permeability properties can be done directly in the time domain. The technique, first introduced by Wang et al. (1994) for causal and diffusive electromagnetic fields and subsequently implemented by Zhdanov and Portniaguine (1997) in the framework of iterative migration, is based upon imaging methods originally developed for seismic wavefields (Claerbout, 1971; Tarantola, 1984). In this paper we advance the original derivations of Wang et al. (1994) and Zhdanov and Portniaguine (1997) to treat non-causal TEM fields, as well as correct a flaw inmore » the theory for treatment of magnetic field data. Our 3D imaging scheme is based on a conjugate-gradient search for the minimum of an error functional involving EM measurements governed by Maxwell's equations without displacement currents. Treatment for magnetic field, voltage (time derivative of the magnetic field) and electric field data are given. The functional can be computed by propagating the data errors back into the model in reverse time along with a DC field, sourced by the integrated data errors over the measurement time range. By correlating these fields, including the time-integrated back-propagated fields, with the corresponding incident field and its initial value at each image point, efficient computational forms for the gradients are developed. The forms of the gradients allow for additional efficiencies when voltage and electric field data are inverted. In such instances the combined data errors can be back-propagated jointly, significantly reducing the computation time required to solve the inverse problem. The inversion algorithm is applied to the long offset transient electromagnetic measurement (LOTEM) configuration thereby demonstrating its capability in inverting non-causal field measurements of electric field and voltage, sourced by a grounded wire, over complex structures. Findings also show that migration, without iteration or preconditioning, is not an effective imaging strategy; reconstructions at the first inversion iteration bear little resemblance to simple or complex test models.« less
  • The growing use of the controlled-source electromagnetic method (CSEM) and magnetotellurics (MT) for exploration applications has been driving the development of data acquisition technologies, and three-dimensional (3-D) modeling and imaging techniques. However, targeting increasingly complex geological environments also further enhances the problems inherent in large-scale inversion, such as non-uniqueness and resolution issues. In this paper, we report on two techniques to mitigate these problems. We use 3-D joint CSEM and MT inversion to improve the model resolution. To avoid the suppression of the resolution capacities of one data type, and thus to balance the use of inherent, and ideally complementarymore » information content, different data reweighting schemes are proposed. Further, a hybrid model parameterization approach is presented, where traditional cell-based model parameters are used simultaneously within a parametric inversion. The idea is to limit the non-uniqueness problem, typical for 3-D imaging problems, in order to allow for a more focusing inversion. The methods are demonstrated using synthetic data generated from models with a strong practical relevance.« less
  • Linear regularization has been applied to the HL-2A infrared imaging bolometer to reconstruct local plasma emission with one-dimensional (1D) and three-dimensional (3D) modeling under the assumption of toroidal symmetry. In the 3D modeling, a new method to calculate the detector point response function is introduced. This method can be adapted to an arbitrarily shaped pinhole. With the full 3D treatment of the detector geometry, up to 50% of the mean-squared error is reduced compared with the 1D modeling. This is attributed to the effects of finite detector size being taken into account in the 3D modeling. Meanwhile, the number ofmore » the bolometer pixels has been optimized to 20 × 20 by making a trade-off between the number of bolometer pixels and the sensitivity of the system. The plasma radiated power density distributions have been calculated as a demonstration using 1D modeling and 3D modeling, respectively.« less
  • The estimation of Q values and/or source corner frequencies f{sub c} from single-station narrow-band recordings of microearthquake spectra is a strongly nonunique problem. This is due to the fact that the spectra can be equally well fitted with low-Q/high-f{sub c} or a high-Q/low f{sub c} spectral models. Here, a method is proposed to constrain this ambiguity by inverting a set of microearthquake spectra for a three-dimensional Q model structure and model source parameters seismic moment (M{sub o}) and corner frequency (f{sub c}) simultaneously. The inversion of whole path Q can be stated as a linear problem in the attenuation operatormore » l* and solved using a tomographic reconstruction of the three-dimensional Q structure. This Q structure is then used as a geometrical constraint for a nonlinear Marquardt-Levenberg inversion of M{sub o} and f{sub c} and a new Q value. The method has been tested on a set of approximately 2,800 P wave spectra from the recordings of 635 microearthquakes from the Kaoiki seismic zone in Hawaii which were recorded at up to six stations. The hypocenters are distributed within a volume of approximately 18 {times} 18 {times} 15 km (depth). The Q model uncertainties have been estimated on the basis of several different tests: Self-consistency, constraining the corner frequencies, and additionally splitting the data set. The standard deviation of the final Q model which used a grid size of 1.5 {times} 1.5 {times} 2.0 km (depth) was less than 3% for the depth range 0-5 km, less than 5% between 5 and 7 km, and 7% between 7 and 9 km. For the given data set and depths below 3-5 km, the method is believed to be able to resolve the model dependent attenuation structure on a scale down to 1-2 km with a resolution of a few percent.« less