Method of focusing neural network approximators for solving a nonlinear inverse problem of geophysics

Author(s):  
I.E. Obornev ◽  
M.I. Shimelevich ◽  
E.A. Obornev ◽  
E.A. Rodionov
Geophysics ◽  
1994 ◽  
Vol 59 (9) ◽  
pp. 1327-1341 ◽  
Author(s):  
Douglas W. Oldenburg ◽  
Yaoguo Li

We develop three methods to invert induced polarization (IP) data. The foundation for our algorithms is an assumption that the ultimate effect of chargeability is to alter the effective conductivity when current is applied. This assumption, which was first put forth by Siegel and has been routinely adopted in the literature, permits the IP responses to be numerically modeled by carrying out two forward modelings using a DC resistivity algorithm. The intimate connection between DC and IP data means that inversion of IP data is a two‐step process. First, the DC potentials are inverted to recover a background conductivity. The distribution of chargeability can then be found by using any one of the three following techniques: (1) linearizing the IP data equation and solving a linear inverse problem, (2) manipulating the conductivities obtained after performing two DC resistivity inversions, and (3) solving a nonlinear inverse problem. Our procedure for performing the inversion is to divide the earth into rectangular prisms and to assume that the conductivity σ and chargeability η are constant in each cell. To emulate complicated earth structure we allow many cells, usually far more than there are data. The inverse problem, which has many solutions, is then solved as a problem in optimization theory. A model objective function is designed, and a “model” (either the distribution of σ or η)is sought that minimizes the objective function subject to adequately fitting the data. Generalized subspace methodologies are used to solve both inverse problems, and positivity constraints are included. The IP inversion procedures we design are generic and can be applied to 1-D, 2-D, or 3-D earth models and with any configuration of current and potential electrodes. We illustrate our methods by inverting synthetic DC/IP data taken over a 2-D earth structure and by inverting dipole‐dipole data taken in Quebec.


Author(s):  
Lukas Hecker ◽  
Rebekka Rupprecht ◽  
Ludger Tebartz van Elst ◽  
Juergen Kornmeier

AbstractEEG and MEG are well-established non-invasive methods in neuroscientific research and clinical diagnostics. Both methods provide a high temporal but low spatial resolution of brain activity. In order to gain insight about the spatial dynamics of the M/EEG one has to solve the inverse problem, which means that more than one configuration of neural sources can evoke one and the same distribution of EEG activity on the scalp. Artificial neural networks have been previously used successfully to find either one or two dipoles sources. These approaches, however, have never solved the inverse problem in a distributed dipole model with more than two dipole sources. We present ConvDip, a novel convolutional neural network (CNN) architecture that solves the EEG inverse problem in a distributed dipole model based on simulated EEG data. We show that (1) ConvDip learned to produce inverse solutions from a single time point of EEG data and (2) outperforms state-of-the-art methods (eLORETA and LCMV beamforming) on all focused performance measures. (3) It is more flexible when dealing with varying number of sources, produces less ghost sources and misses less real sources than the comparison methods. (4) It produces plausible inverse solutions for real-world EEG recordings and needs less than 40 ms for a single forward pass. Our results qualify ConvDip as an efficient and easy-to-apply novel method for source localization in EEG and MEG data, with high relevance for clinical applications, e.g. in epileptology and real time applications.


Sign in / Sign up

Export Citation Format

Share Document