Inversion of slingram electromagnetic induction data using a Born approximation

Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. E201-E212 ◽  
Author(s):  
Jochen Kamm ◽  
Michael Becken ◽  
Laust B. Pedersen

We present an efficient approximate inversion scheme for near-surface loop-loop EM induction data (slingram) that can be applied to obtain 2D or 3D models on a normal desktop computer. Our approach is derived from a volume integral equation formulation with an arbitrarily conductive homogeneous half-space as a background model. The measurements are not required to fulfill the low induction number condition (low frequency and conductivity). The high efficiency of the method is achieved by invoking the Born approximation around a half-space background. The Born approximation renders the forward operator linear. The choice of a homogeneous half-space yields closed form expressions for the required electromagnetic normal fields. It also yields a translationally invariant forward operator, i.e., a highly redundant Jacobian. In connection with the application of a matrix-free conjugate gradient method, this allows for very low memory requirements during the inversion, even in three dimensions. As a consequence of the Born approximation, strong conductive deviations from the background model are underestimated. Highly resistive anomalies are in principle overestimated, but at the same time difficult to resolve with induction methods. In the case of extreme contrasts, our forward model may fail in simultaneously explaining all the data collected. We applied the method to EM34 data from a profile that has been extensively studied with other electromagnetic methods and compare the results. Then, we invert three conductivity maps from the same area in a 3D inversion.

Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1769-1778 ◽  
Author(s):  
Sándor Szalai ◽  
László Szarka ◽  
Erno˝ Prácser ◽  
Frank Bosch ◽  
Imre Müller ◽  
...  

The term “null array” is introduced for those electrode configurations where the measured potential difference is zero above a homogeneous half‐space when using a measuring dipole M0N0. Different types of null arrays (three‐electrode, Schlumberger, and dipole axial/equatorial null arrays) and their corresponding traditional arrays are studied. It was shown in a field study carried out in a karstified limestone area covered by thin sediments that it is possible to obtain geologically meaningful results with null‐array techniques. The main features of the null‐array data are as follows. (1) Null‐array data appear to be more spatially variable than the classical data. The spatial variability provides information about the presence of karstic fractures in the subsurface; (2) The null‐array anomalies caused by nearly vertical karstic fractures in the limestone basement do not decay with depth as quickly as the classical array anomalies. (3) The strike direction of the fractures is much less ambiguous than that found by using classical arrays. Nevertheless, the depth variation of the basement is more reliably observed in geoelectric anomalies obtained using traditional arrays. Therefore a joint use of classical arrays and their corresponding null methods is recommended, because the combined methods provide more information about the subsurface structure.


Geophysics ◽  
1989 ◽  
Vol 54 (3) ◽  
pp. 350-358 ◽  
Author(s):  
G. Nolet ◽  
R. Sleeman ◽  
V. Nijhof ◽  
B. L. N. Kennett

We present a simple algorithm for computing the acoustic response of a layered structure containing three‐dimensional (3-D) irregularities, using a locked‐mode approach and the Born approximation. The effects of anelasticity are incorporated by use of Rayleigh’s principle. The method is particularly attractive at somewhat larger offsets, but computations for near‐source offsets are stable as well, due to the introduction of anelastic damping. Calculations can be done on small minicomputers. The algorithm developed in this paper can be used to calculate the response of complicated models in three dimensions. It is more efficient than any other method whenever many sources are involved. The results are useful for modeling, as well as for generating test signals for data processing with realistic, model‐induced “noise.” Also, this approach provides an alternative to 2-D finite‐difference calculations that is efficient enough for application to large‐scale inverse problems. The method is illustrated by application to a simple 3-D structure in a layered medium.


1978 ◽  
Vol 45 (4) ◽  
pp. 812-816 ◽  
Author(s):  
B. S. Berger ◽  
B. Alabi

A solution has been derived for the Navier equations in orthogonal cylindrical curvilinear coordinates in which the axial variable, X3, is suppressed through a Fourier transform. The necessary coordinate transformation may be found either analytically or numerically for given geometries. The finite-difference forms of the mapped Navier equations and boundary conditions are solved in a rectangular region in the curvilinear coordinaties. Numerical results are given for the half space with various surface shapes and boundary conditions in two and three dimensions.


1996 ◽  
Vol 86 (3) ◽  
pp. 750-760
Author(s):  
F. Abramovici ◽  
L. H. T. Le ◽  
E. R. Kanasewich

Abstract This article presents some numerical experiments in using a computer program for calculating the displacements due to a P source in a vertically inhomogeneous structure, based on the Fourier-Bessel representation. The structure may contain homogeneous, inhomogeneous, elastic, or viscoelastic layers. The source may act in any type of sublayer or in the half-space. Synthetic results for the simple case of a homogeneous layer overlaying a homogeneous half-space compare favorably with computations based on the Cagniard method. Numerical seismograms for an elastic layer having velocities and density varying linearly with depth were computed by integrating numerically the governing differential systems and compared with results based on the Haskell model of splitting the linear layer in homogeneous sublayers. Even an adaptive process with a variable step size based on the Haskell model has a poorer performance on the accuracy-cpu time scale than numerical integration.


1970 ◽  
Vol 60 (2) ◽  
pp. 321-344 ◽  
Author(s):  
Fred Schwab ◽  
Leon Knopoff

abstract Fundamental-mode Love- and Rayleigh-wave dispersion computations for multilayered, perfectly-elastic media were studied. The speed of these computations was improved, and the accuracy brought under full control. With sixteen decimal digits employed in these computations, fifteen significant-figure accuracy was found possible with Love waves and twelve to thirteen figure accuracy with Rayleigh waves. In order to ensure that the computed dispersion is correct to a specified accuracy, say σ significant figures, (σ + 1)/4 wavelengths of layered structure must be retained above a homogeneous half-space. To this accuracy, the homogeneous half-space is a sufficient model of the true layering it replaces. Using this result, it was possible to refine the usual layer-reduction technique so as to ensure retention of the specified accuracy while employing reduction. With this reduction technique in effect, and with σ specified below single-precision accuracy, the program can be run entirely in single precision; the specified accuracy is maintained without overflow or loss-of-precision problems being encountered during calculations.


2019 ◽  
Vol 879 ◽  
pp. 808-833 ◽  
Author(s):  
B. J. Walker ◽  
K. Ishimoto ◽  
H. Gadêlha ◽  
E. A. Gaffney

We present a generalisation of efficient numerical frameworks for modelling fluid–filament interactions via the discretisation of a recently developed, non-local integral equation formulation to incorporate regularised Stokeslets with half-space boundary conditions, as motivated by the importance of confining geometries in many applications. We proceed to utilise this framework to examine the drag on slender inextensible filaments moving near a boundary, firstly with a relatively simple example, evaluating the accuracy of resistive force theories near boundaries using regularised Stokeslet segments. This highlights that resistive force theories do not accurately quantify filament dynamics in a range of circumstances, even with analytical corrections for the boundary. However, there is the notable and important exception of movement in a plane parallel to the boundary, where accuracy is maintained. In particular, this justifies the judicious use of resistive force theories in examining the mechanics of filaments and monoflagellate microswimmers with planar flagellar patterns moving parallel to boundaries. We proceed to apply the numerical framework developed here to consider how filament elastohydrodynamics can impact drag near a boundary, analysing in detail the complex responses of a passive cantilevered filament to an oscillatory flow. In particular, we document the emergence of an asymmetric periodic beating in passive filaments in particular parameter regimes, which are remarkably similar to the power and reverse strokes exhibited by motile$9+2$cilia. Furthermore, these changes in the morphology of the filament beating, arising from the fluid–structure interactions, also induce a significant increase in the hydrodynamic drag of the filament.


Sign in / Sign up

Export Citation Format

Share Document