3D inversion of the time-domain electromagnetic data for exploration of submarine hydrothermal deposits using the GEMTIP model

Author(s):  
Masashi Endo ◽  
Leif Cox ◽  
David Sunwall ◽  
Michael Zhdanov ◽  
Eiichi Asakawa
Geophysics ◽  
2013 ◽  
Vol 78 (1) ◽  
pp. E47-E57 ◽  
Author(s):  
Douglas W. Oldenburg ◽  
Eldad Haber ◽  
Roman Shekhtman

We present a 3D inversion methodology for multisource time-domain electromagnetic data. The forward model consists of Maxwell’s equations in time where the permeability is fixed but electrical conductivity can be highly discontinuous. The goal of the inversion is to recover the conductivity-given measurements of the electric and/or magnetic fields. The availability of matrix-factorization software and high-performance computing has allowed us to solve the 3D time domain EM problem using direct solvers. This is particularly advantageous when data from many transmitters and over many decades are available. We first formulate Maxwell’s equations in terms of the magnetic field, [Formula: see text]. The problem is then discretized using a finite volume technique in space and backward Euler in time. The forward operator is symmetric positive definite and a Cholesky decomposition can be performed with the work distributed over an array of processors. The forward modeling is quickly carried out using the factored operator. Time savings are considerable and they make 3D inversion of large ground or airborne data sets feasible. This is illustrated by using synthetic examples and by inverting a multisource UTEM field data set acquired at San Nicolás, which is a massive sulfide deposit in Mexico.


2013 ◽  
Vol 2013 (1) ◽  
pp. 1-4
Author(s):  
Christoph Schwarzbach ◽  
Elliot Holtham ◽  
Eldad Haber

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. E287-E295 ◽  
Author(s):  
Aihua Weng ◽  
Dajun Li ◽  
Yue Yang ◽  
Sirui Li ◽  
Jianping Li ◽  
...  

Wide applications of time-domain electromagnetic (TEM) data require 3D inversion. A possible strategy is to use the developed 3D inversion algorithms in frequency-domain (FD) electromagnetic (EM) methods. Thus, the key of the strategy is how to transform the time-domain ([Formula: see text]-[Formula: see text]) EM signal into the FD. An inversion algorithm has been developed to transform the [Formula: see text]-[Formula: see text] signal into a corresponding FD response. In this method, a step-off current is presumed. Under this assumption, the Fourier transform relating the EM FD response to the [Formula: see text]-[Formula: see text] signal becomes a sine or cosine transformation. Using the polynomial approximation method, the transformation turns into a linear equation. From a set of [Formula: see text]-[Formula: see text] signals, FD responses could be obtained by solving these linear equations in the least-squares sense. To reduce the nonuniqueness of the solution, and enhance the solution stability, an additional smoothness constraint on the FD response is imposed, thus converting the minimization problem into a regularization inversion problem. The algorithm is applied to synthetic and field vertical magnetic data in the in-loop TEM surveying mode. The numerical results show that in the entire audio-frequency range, the relative errors between the inversed and theoretical FD responses of the real and imaginary parts are almost all less than 1%, with the largest discrepancy of 5% occurring at high frequencies. There are two significances behind our work: First, the possibility of accurately transforming [Formula: see text]-[Formula: see text] response into FD response in audio-frequency range is coming into true, thereby (from the mathematical perspective) implementing the equivalence between the responses of the EM method in the time domain and the FD. Second, the algorithm provides a new approach to interpret TEM data in 3D mode by using developed 3D FEM inversion techniques.


Geophysics ◽  
2014 ◽  
Vol 79 (6) ◽  
pp. E303-E314 ◽  
Author(s):  
David Marchant ◽  
Eldad Haber ◽  
Douglas W. Oldenburg

Understanding the effects of induced-polarization (IP) effects on time-domain electromagnetic data requires the ability to simulate common survey techniques when taking chargeability into account. Most existing techniques preform this modeling in the frequency domain prior to transforming their results to the time domain. Even though this technique can allow for chargeable material to be easily incorporated, its application for some problems can be computationally limiting. We developed a new technique for forward modeling the time-domain electromagnetic response of chargeable materials in three dimensions. The frequency dependence of Ohms’ law translates to an ordinary differential equation when considered in the time domain. The system of ordinary-partial differential equations was then discretized using an implicit time-stepping algorithm, that yielded absolute stability. This approach allowed us to operate directly in the time domain and avoid frequency to time-domain transformations. Although this approach can be applied directly to materials exhibiting Debye dispersions, other Cole-Cole dispersions resulted in fractional derivatives in time. To overcome this difficulty, Padé approximations were used to represent the frequency dependence as a rational series of integer order terms. The resulting method was then simplified to generate a reduced time-domain model that can be used to forward model the IP decay curves in the absence of any electromagnetic coupling. We found numerical examples in which the method produced accurate results. The potential application of the method was demonstrated by modeling the full time-domain electromagnetic response of a gradient array IP survey, and the occurrence of negative transients in airborne time-domain electromagnetic data.


2021 ◽  
Vol 11 (5) ◽  
pp. 2060 ◽  
Author(s):  
Alexander Parshin ◽  
Ayur Bashkeev ◽  
Yuriy Davidenko ◽  
Marina Persova ◽  
Sergey Iakovlev ◽  
...  

Nowadays in solving geological problems, the technologies of UAV-geophysics, primarily magnetic and gamma surveys, are being increasingly used. However, for the formation of the classical triad of airborne geophysics methods in the UAV version, there was not enough technology for UAV-electromagnetic sounding, which would allow studying the geological environment at depths of tens and hundreds of meters with high detail. This article describes apparently the first technology of UAV-electromagnetic sounding in the time domain (TDEM, TEM), implemented as an unmanned system based on a light multi-rotor UAV. A measuring system with an inductive sensor—an analogue of a 20 × 20 or 50 × 50 m receiving loop is towed by a UAV, and a galvanically grounded power transmitter is on the ground and connected to a pulse generator. The survey is carried out along a network of parallel lines at low altitude with a terrain draping at a speed of 7–8 m/s, the maximum distance of the UAV’s departure from the transmitter line can reach several kilometers, thus the created technology is optimal for performing detailed areal electromagnetic soundings in areas of several square kilometers. The results of the use of the unmanned system (UAS) in real conditions of the mountainous regions of Eastern Siberia are presented. Based on the obtained data, the sensitivity of the system was simulated and it was shown that the developed technology allows one to collect informative data and create geophysical sections and maps of electrical resistivity in various geological situations. According to the authors, the emergence of UAV-TEM systems in the near future will significantly affect the practice of geophysical work, as it was earlier with UAV-magnetic prospecting and gamma-ray survey.


Sign in / Sign up

Export Citation Format

Share Document