Surface instability of a homogeneous half-space coupled with a finite number of laminas

1995 ◽  
Vol 76 (3) ◽  
pp. 2343-2347
Author(s):  
V. N. Chekhov ◽  
V. S. Stukotilov
Keyword(s):  
Finite Number ◽  
Half Space ◽  
Surface Instability ◽  
Homogeneous Half Space

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
Keyword(s):  
Half Space ◽  
Born Approximation ◽  
High Efficiency ◽  
Three Dimensions ◽  
Background Model ◽  
Desktop Computer ◽  
Near Surface ◽  
Homogeneous Half Space ◽  
Integral Equation Formulation ◽  
Forward Operator

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.

Numerical seismograms obtained using ω-κ integrals, for a point P source in a vertically inhomogeneous anelastic model

1996 ◽  
Vol 86 (3) ◽  
pp. 750-760
Author(s):  
F. Abramovici ◽  
L. H. T. Le ◽  
E. R. Kanasewich
Keyword(s):  
Half Space ◽  
Homogeneous Layer ◽  
Variable Step Size ◽  
Step Size ◽  
Adaptive Process ◽  
Cpu Time ◽  
Homogeneous Half Space ◽  
Linear Layer ◽  
Variable Step ◽  
Viscoelastic Layers

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.

Surface-wave dispersion computations

1970 ◽  
Vol 60 (2) ◽  
pp. 321-344 ◽  
Author(s):  
Fred Schwab ◽  
Leon Knopoff
Keyword(s):  
Half Space ◽  
Rayleigh Waves ◽  
Wave Dispersion ◽  
Reduction Technique ◽  
Surface Wave Dispersion ◽  
Single Precision ◽  
Homogeneous Half Space ◽  
Full Control ◽  
Significant Figures ◽  
Rayleigh Wave Dispersion

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.

scholarly journals Rayleigh-type waves on a coated elastic half-space with a clamped surface

2019 ◽  
Vol 377 (2156) ◽  
pp. 20190111 ◽  
Author(s):  
J. Kaplunov ◽  
D. Prikazchikov ◽  
L. Sultanova
Keyword(s):  
Half Space ◽  
Numerical Solutions ◽  
Substrate Surface ◽  
Numerical Data ◽  
Elastic Half Space ◽  
Asymptotic Results ◽  
Homogeneous Half Space ◽  
Upper Face ◽  
Localized Waves ◽  
Perturbed Equation

Elastodynamics of a half-space coated by a thin soft layer with a clamped upper face is considered. The focus is on the analysis of localized waves that do not exist on a clamped homogeneous half-space. Non-traditional effective boundary conditions along the substrate surface incorporating the effect of the coating are derived using a long-wave high-frequency procedure. The derived conditions are implemented within the framework of the earlier developed specialized formulation for surface waves, resulting in a perturbation of the shortened equation of surface motion in the form of an integral or pseudo-differential operator. Non-uniform asymptotic formula for the speeds of the sought for Rayleigh-type waves, failing near zero frequency and the thickness resonances of a layer with both clamped faces, follow from the aforementioned perturbed equation. Asymptotic results are compared with the numerical solutions of the full dispersion relation for a clamped coated half-space. A similarity with Love-type waves proves to be useful for interpreting numerical data. This article is part of the theme issue ‘Modelling of dynamic phenomena and localization in structured media (part 1)’.

Computation of Green’s tensor integrals for three‐dimensional electromagnetic problems using fast Hankel transforms

Geophysics ◽  
1984 ◽  
Vol 49 (10) ◽  
pp. 1754-1759 ◽  
Author(s):  
Walter L. Anderson
Keyword(s):  
Half Space ◽  
Three Dimensional ◽  
Hankel Transform ◽  
Electromagnetic Modeling ◽  
Electromagnetic Problems ◽  
Hankel Transforms ◽  
Homogeneous Half Space ◽  
Green's Tensor ◽  
Green’S Tensor ◽  
The Many

A new method is presented that rapidly evaluates the many Green’s tensor integrals encountered in three‐dimensional electromagnetic modeling using an integral equation. Application of a fast Hankel transform (FHT) algorithm (Anderson, 1982) is the basis for the new solution, where efficient and accurate computation of Hankel transforms are obtained by related and lagged convolutions (linear digital filtering). The FHT algorithm is briefly reviewed and compared to earlier convolution algorithms written by the author. The homogeneous and layered half‐space cases for the Green’s tensor integrals are presented in a form so that the FHT can be easily applied in practice. Computer timing runs comparing the FHT to conventional direct convolution methods are discussed, where the FHT’s performance was about 6 times faster for a homogeneous half‐space, and about 108 times faster for a five‐layer half‐space. Subsequent interpolation after the FHT is called is required to compute specific values of the tensor integrals at selected transform arguments; however, due to the relatively small lagged convolution interval used (same as the digital filter’s), a simple and fast interpolation is sufficient (e.g., by cubic splines).

Sign in / Sign up

Export Citation Format

Share Document