Inversion of potential field data in Fourier transform domain

Geophysics ◽  
1985 ◽  
Vol 50 (4) ◽  
pp. 685-691 ◽  
Author(s):  
J. C. Mareschal

A relationship is derived between the Fourier transform of a potential field at the Earth’s surface and the transform of the inducing source distribution. The Fourier transform of the field is the Laplace transform of the source distribution spectrum when the Laplace transform variable p is equal to the wavenumber. This relationship can be used to determine all possible source distributions compatible with the data. The solution is the superposition of a particular solution to an inhomogeneous problem and of the general solution to the homogeneous problem (i.e., for which the field vanishes at the surface). Source distribution can be expanded into a set of known functions; coefficients of the expansion are determined by solving a system of linear equations. Physical constraints can be introduced to restrict the variation range of the coefficients of expansion. Two examples are presented to illustrate the method: a synthetic gravity profile and a heat flow profile are inverted to determine density or heat source distributions compatible with the data.

Author(s):  
Michael J Corinthios

In this paper, the author uses his recently proposed complex variable generalized distribution theory to expand the domains of existence of bilateral Laplace and z transforms, as well as a whole new class of related transforms. A vast expansion of the domains of existence of bilateral Laplace and z transforms and continuous-time and discrete-time Hilbert, Hartley and Mellin transforms, as well as transforms of multidimensional functions and sequences are obtained. It is noted that the Fourier transform and its applications have advanced by leaps and bounds during the last century, thanks to the introduction of the theory of distributions and, in particular, the concept of the Dirac-delta impulse. Meanwhile, however, the truly two-sided ‘bilateral’ Laplace and z transforms, which are more general than Fourier, remained at a standstill incapable of transforming the most basic of functions. In fact, they were reduced by half to one-sided transforms and received no more than a passing reference in the literature. It is shown that the newly proposed generalized distributions expand the domains of existence and application of Laplace and z transforms similar to and even more extensively than the expansion of the domain of Fourier transform that resulted from the introduction, nearly a century ago, of the theory of distributions and the Dirac-delta impulse. It is also shown that the new generalized distributions put an end to an anomaly that still exists today, which meant that for a large class of basic functions, the Fourier transform exists while the more general Laplace and z transforms do not. The anomaly further manifests itself in the fact that even for the one-sided causal functions, such as the Heaviside unit step function u ( t ) and the sinusoid sin βtu ( t ), the Laplace transform does not exist on the j ω -axis, and the Fourier transform which does exist cannot be deduced thereof by the substitution s =j ω in the Laplace transform, which by definition it should. The extended generalized transforms are well defined for a large class of functions ranging from the most basic to highly complex fast-rising exponential ones that have so far had no transform. Among basic applications, the solution of partial differential equations using the extended generalized transforms is provided. This paper clearly presents and articulates the significant impact of extending the domains of Laplace and z transforms on a large family of related transforms, after nearly a century during which bilateral Laplace and z transforms of even the most basic of functions were undefined, and the domains of definition of related transforms such as Hilbert, Hartley and Mellin transforms were confined to a fraction of the space they can now occupy.


Filomat ◽  
2020 ◽  
Vol 34 (9) ◽  
pp. 2869-2876
Author(s):  
H.M. Srivastava ◽  
Mohammad Masjed-Jamei ◽  
Rabia Aktaş

This article deals with a general class of differential equations and two general classes of integral equations. By using the Laplace transform and the Fourier transform, analytical solutions are derived for each of these classes of differential and integral equations. Some illustrative examples and particular cases are also considered. The various analytical solutions presented in this article are potentially useful in solving the corresponding simpler differential and integral equations.


Author(s):  
A Michel ◽  
J-P Boy

Summary Long term deformations strongly depend on the Earth model and its rheological parameters, and in particular its viscosity. We give the general theory and the numerical scheme to compute them for any spherically non rotating isotropic Earth model with linear rheology, either elastic or viscoelastic. Although the Laplace transform is classically used to compute viscoelastic deformation, we choose here instead, to implement the integration with the Fourier transform in order to take advantage of the Fast Fourier Transform algorithm and avoid some of the Laplace transform mathematical difficulties. We describe the methodology to calculate deformations induced by several geophysical signals regardless of whether they are periodic or not, especially by choosing an adapted time sampling for the Fourier transform. As examples, we investigate the sensitivity of the displacements due to long period solid Earth tides, Glacial Isostatic Adjustment (GIA), and present-day ice melting, to anelastic parameters of the mantle. We find that the effects of anelasticity are important for long period deformation and relatively low values of viscosities for both Maxwell and Burgers models. We show that slight modifications in the rheological models could significantly change the amplitude of deformation but also affect the spatial and temporal pattern of the signal to a lesser extent. Especially, we highlight the importance of the mantle anelasticity in the low degrees deformation due to present-day ice melting and encourage its inclusion in future models.


1991 ◽  
Vol 45 (1) ◽  
pp. 59-70 ◽  
Author(s):  
Riccardo Croci

The purpose of this paper is to derive the asymptotic solutions to a class of inhomogeneous integral equations that reduce to algebraic equations when a parameter ε goes to zero (the kernel becoming proportional to a Dirac δ function). This class includes the integral equations obtained from the system of Vlasov and Poisson equations for the Fourier transform in space and the Laplace transform in time of the electrostatic potential, when the equilibrium magnetic field is uniform and the equilibrium plasma density depends on εx, with the co-ordinate z being the direction of the magnetic field. In this case the inhomogeneous term is given by the initial conditions and possibly by sources, and the Laplace-transform variable ω is the eigenvalue parameter.


2019 ◽  
Vol 23 (Suppl. 6) ◽  
pp. 1929-1934
Author(s):  
Resat Yilmazer ◽  
Neslihan Demirel

In this study, we obtain the solution of a local fractional Schrodinger equation. The solution is obtained by the implementation of the Laplace transform and Fourier transform in closed form in terms of the Mittag-Leffler function.


1968 ◽  
Vol 8 (3) ◽  
pp. 557-567 ◽  
Author(s):  
M. J. Englefield

Laplace transform techniques for solving differential equations do not seem to have been directly applied to the Schrödinger equation in quantum mechanics. This may be because the Laplace transform of a wave function, in contrast to the Fourier transform, has no direct physical significance. However, this paper will show that scattering phase shifts and bound state energies can be determined from the singularities of the Laplace transform of the wave function. The Laplace transform method can thereby simplify calculations if the potential allows a straightforward solution of the transformed Schrödinger equation. Suitable cases are the Coulomb, oscillator and exponential potentials and the Yamaguchi separable non-local potential.


2017 ◽  
Vol 120 (2) ◽  
pp. 225 ◽  
Author(s):  
Marcus Carlsson ◽  
Jens Wittsten

We revisit the classical problem of when a given function, which is analytic in the upper half plane $\mathbb{C} _+$, can be written as the Fourier transform of a function or distribution with support on a half axis $(-\infty ,b]$, $b\in \mathbb{R} $. We derive slight improvements of the classical Paley-Wiener-Schwartz Theorem, as well as softer conditions for verifying membership in classical function spaces such as $H^p(\mathbb{C} _+)$.


Geophysics ◽  
1990 ◽  
Vol 55 (5) ◽  
pp. 549-555 ◽  
Author(s):  
Mark Pilkington ◽  
W. E. S. Urquhart

Most existing techniques for potential field data enhancement and interpretation require data on a horizontal plane. Hence, when observations are made on an irregular surface, reduction to a horizontal plane is necessary. To effect this reduction, an equivalent source distribution that models the observed field is computed on a mirror image of the observation surface. This irregular mirror image surface is then replaced by a horizontal plane and the effect of the equivalent sources is computed on the required horizontal level. This calculated field approximates the field reduced to a horizontal plane. The good quality of this approximation is demonstrated by two‐dimensional synthetic data examples in which the maximum errors occur in areas of steep topographic gradients and increased magnetic field intensity. The approach is also applied to a portion of a helicopter‐borne aeromagnetic survey from the Gaspé region in Quebec, Canada, where the results are a horizontal shifting of anomaly maxima of up to 150 m and changes in anomaly amplitudes of up to 100 nT.


Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. J9-J25 ◽  
Author(s):  
Mark Pilkington ◽  
Olivier Boulanger

The continuation of potential field data from one irregular surface to another, not always horizontal, is often a necessary component within the data processing and interpretation stream. The most common requirement is to reduce field values (or some related component or derivative) to a horizontal plane, to facilitate further quantitative processing. Methods available to continue data comprise two main approaches. The first (source-based) involves calculating a source distribution that produces a fit to the data and can be used to calculate the field at any other point above. The second (field-based) requires no source determinations and deals with only fields but may involve calculating the field on some intermediate surface. Nine different continuation methods were compared (four source based and five field based) through synthetic tests and on real data from a helicopter-borne survey in Yukon, Canada. The preferred methods of Guspi and Hansen are those that do not involve any theoretical or geometric approximations and involve intermediate calculations on a plane or surface close to the observation surface. The Guspi approach is faster, based on using frequency-domain processing, but the Hansen method uses equivalent sources close enough to and consistently below the observation surface so that no low-pass filtering needs to be used.


Geophysics ◽  
1977 ◽  
Vol 42 (7) ◽  
pp. 1384-1393 ◽  
Author(s):  
Anas M. Abo‐Zena

For an elastic material with an infinite circular cylindrical hole, the exact solution due to a pressure on a finite length of the cylinder is obtained as a function of the Laplace transform parameter on time and Fourier transform parameter on the z-coordinate (the axis of the cylinder). The applied pressure is a function of the time and the position z. Numerical inversion of the Laplace and Fourier transforms are required to determine the field quantities in the time and space parameters. In the far field, the inverse Fourier transform can be obtained by an asymptotic expansion. It remains to obtain the inverse Laplace transform numerically. We have found that for cylinders whose radius is small compared with the smallest wavelength of interest, an analytical solution can be obtained. Graphical results for the cases of instantaneous explosion and progression of the detonation with constant velocity are given. In both cases an exponential decay of the explosion pressure is assumed.


Sign in / Sign up

Export Citation Format

Share Document