The fourier transform is sometimes denoted by the operator fand its inverse by f1, so that. There is also an inverse fourier transform that mathematically synthesizes the original function from its frequency domain representation. Antileakage fourier transform for seismic data regularization. For fourier transform purposes, sampling is modeled as a product between st and a dirac comb function. The frequency quantum is the change in frequency that is represented by two sequential bucket numbers. As a consequence, energy from one fourier coefficient leaks onto other coefficients. However, the computational cost dra matically increases with the sampling dimensions. The antileakage fourier transform alft is a regularization method using an iterative procedure for computing the spectrum of irregularly sampled data. Introduction to fourier transforms fourier transform as a limit of the fourier series inverse fourier transform. Breve introduzione al mondo dellelaborazione dei segnali in ambito digitale.
This method uses discrete fourier transform dft, although for regular data with missing traces or gaps it can use the fast fourier transform fft zwartjes and sacchi 2007. Ep2103959a2 method for interpolating seismic data by anti. However, this is a very sharpened window that causes broad leakage. Matching pursuit fourier interpolation using priors derived.
Transform each trace from time domain to frequency domain using the fast fourier transform fft generate a data series using the available traces for each frequency a frequency slice use the antileakage lssa to regularize each frequency slice transform each trace from frequency domain to time domain using the inverse fft. Anti leakage fourier transform the first documented solution to the leakage problem was performed by hogbom 1974 as a solution within astrophysics. Even if the window contains integer multiples of a period erroneous spectral components occur. Summary the antileakage fourier transform alft is a regularization method using an iterative procedure for computing the spectrum of irregularly. Discrete fourier transform dft is the basic part of various algorithms of signal processing in many fields of science and technology. The antileakage fourier transform algorithm can overcome the spectral leakage problem and handles aliased data. Antileakage fourier transform among fourier based seismic interpolation methods, alft is an iterative algorithm that acts on each timefrequency slice separately, and reconstructs a temporal frequency spectrum along spatial dimensions. Many advantages of regularization in the fourier domain come with some serious problems.
These algorithms are established mainly on the basis of fourier transform. The application of fourier analysis to forecasting the. This is an iterative algorithm, which exploits the fact that leakage. One of the popular seismic interpolation techniques is antileakage fourier transform alft which was introduced by xu and pham. Half range fourier sine or cosine series a half range fourier sine or cosine series is a series in which only sine terms or. The reverse transform can then be used either to output the data on a regular grid for full regularisation or to. Pdf antileakage fourier transform for seismic data regularization. Transform each trace from time domain to frequency domain using the fast fourier transform fft generate a data series using the available traces for each frequency a frequency slice use the anti leakage lssa to regularize each frequency slice transform each trace from frequency domain to time domain using the inverse fft. Seismic data regularization with the antialias antileakage. Journal of geophysics and engineering, volume, number 1. Chapter 3 fourier transforms of distributions questions 1 how do we transform a function f. In this talk first several aspects of the alft are illustrated using simple synthetic signals.
Then, the maximum fourier component is selected and transformed back to the irregular grid. Michel schonewille, andreas klaedtke, alan vigner, john brittan and tony martin present a solution to a problem with regularly sampled data found in the use of the. One of the methods for 5d fourier based seismic interpolation is the antileakage fourier transfrom alft. Pdf antileakage fourier transform for seismic data. The inverse fourier transform for linearsystems we saw that it is convenient to represent a signal fx as a sum of scaled and shifted sinusoids. For example, you can effectively acquire timedomain signals, measure the frequency content, and convert the results to realworld units and displays as shown on traditional benchtop. Dec 22, 2010 the antileakage fourier transform algorithm can overcome the spectral leakage problem and handles aliased data. Spectral leakage by douglas lyon abstract this paper is part 4 in a series of papers about the discrete fourier transform dft and the inverse discrete fourier transform idft.
Chapter 1 the fourier transform institute for mathematics. Consequently, spectrum leakage appears in frequency domain. Antileakage fourier transform alft and its application for seismic data reconstruction article pdf available in journal of the earth and space physics 394. D032 seismic data regularization with antileakage fourier transform sheng xu and don pham abstract 1 veritas dcg inc. Spectral leakage applies to all forms of dft, including the fft. The application of fourier analysis to forecasting the inbound call time series of a call centre bruce g. Any other type of operation creates new frequency components that may be referred to as spectral leakage in the broadest sense. Antileakage ndft and its application in data reconstruction ieee. Antileakage fourier transform for seismic data regularization free. However, the computational cost dramatically increases with the sampling dimensions. Digital signal processing algorithms are commonly used to obtain radio spectrum estimates based on measurements. Antileakage fourier transform we have developed an algorithm that performs a forward irregular fourier transform in all five data dimensions simultaneously to build a representation of the data in the frequency domain. Seismic data regularization with antileakage fourier transform. Deriving fourier transform from fourier series consider a periodic signal f with period t.
In the antileakage fourier transform, we first solve for the most energetic fourier coefficient, assuming that it causes the most severe leakage. For analysis of signals with dft, the length of discrete signal must be finite so signal to be analyzed must be divided to some windows. Equation 3 can also be used to obtain the frequency quantum. Town park drive houston tx 77072 usa in the theory of fourier reconstruction from discrete seismic data it aims to estimate the spatial frequency content on an irregularly sampled grid. In the antileakage fourier transform, we first solve for the most energetic. Mar 09, 2015 leakage effect occuring when a discrete fourier transform is performed over a hamming time window. Antileakage fourier transform for seismic data regularization in higher dimensions sheng xu1, yu zhang1, and gilles lambare2 abstract wideazimuth seismic data sets are generally acquired more sparsely than narrowazimuth seismic data sets. D032 seismic data regularization with anti leakage fourier transform sheng xu and don pham abstract 1 veritas dcg inc. To generalize it to higher dimensions, we propose an area weighting scheme to accurately estimate the fourier components. The resulting transform pairs are shown below to a common horizontal scale. Imaging subsalt structures under complex overburdens is a chal. Antileakage fourier transform for seismic data regularization in. The antileakage fourier transform algo rithm can overcome the spectral leakage problem and han dles aliased data. Leakage effect occuring when a discrete fourier transform is performed over a hamming time window.
Resolution and leakage of fourier transform concepts and. The anti leakage fourier transform reconstruction xu et al. Jan 29, 2016 one of the popular seismic interpolation techniques is anti leakage fourier transform alft which was introduced by xu and pham. In practice, the measured signals are limited in time and the fft calculates the frequency transform over a certain number of discrete frequencies called bins. Leakage effect of a discrete fourier transform hamming window. The antileakage fourier transform reconstruction xu et al. The antileakage fourier transform alft is a fourier regularization method. Antileakage leastsquares spectral analysis for data. Antileakage fourier transform the first documented solution to the leakage problem was performed by hogbom 1974 as a solution within astrophysics. The general principle of fourier regularization is to compute the spectrum of irregularly sampled data, and then use an inverse transform to reconstruct data at new locations. Thus the ripples of the sinc function will show up as leakage far from any spectral peak that is not perfectly periodic in the ffts length. Seismic data regularization with the antialias antileakage fourier.
Despite the simplicity of the antileakage fourier transform alft, and its ability to regularize pure nonuniform data sets and handle higher dimension seismic data, alft still remains less appealing because of its iterative algorithm and so its costly computation. An adaptation of this can give a robust solution to fourier regularisation, as discussed by xu and pham 2004a, xu and pham 2004b and schonewille 2009, and has been named the antileakage fourier transform due to its ability to suppress spectral leakage. Bellc answ police assistance line, tuggerah, nsw 2259, email. An estimated frequencywavenumber spectrum is generated by applying a first anti leakage fourier transform method to aliased frequency components in temporaltransformed seismic data and applying a second anti leakage fourier transform method to unaliased frequency components in the temporaltransformed seismic data. And the properties of the discrete fourier transform, as well see, the differences in the properties between this and the fourier transforms weve talked about in previous lectures are related very closely to the implied periodicity in the sequence or in the discrete fourier. Despite the simplicity of the anti leakage fourier transform alft, and its ability to regularize pure nonuniform data sets and handle higher dimension seismic data, alft still remains less appealing because of its iterative algorithm and so its costly computation. Mar 23, 2016 dal dominio del tempo a quello delle frequenze. For each iteration a discrete fourier transform is performed. The key to resolving the spectral leakage is to reduce the leakages among fourier coefficients in the original data before the calculation. An estimated frequencywavenumber spectrum is generated by applying a first antileakage fourier transform method to aliased frequency components in temporaltransformed seismic data and applying a second antileakage fourier transform method to unaliased frequency components in the temporaltransformed seismic data. Anti leakage fourier transform we have developed an algorithm that performs a forward irregular fourier transform in all five data dimensions simultaneously to build a representation of the data in the frequency domain. Fourier series corresponding to an even function, only cosine terms and possibly a constant which we shall consider a cosine term can be present.
In reality, signals are of timelimited nature and nothing can be known about the signal beyond the measured interval. Both resolution and leakage depend on the shape of the data window. We investigate the nonorthogonality of the fourier basis on an irregularly sampled grid and propose a technique called antileakage fourier transform to overcome the spectral leakage. Us20090231956a1 method for interpolating seismic data by. This brings new challenges to seismic data regularization algo. Antialiasing, antileakage fourier transform earthdoc. The fourier transform is not limited to functions of time, but the domain of the original function is commonly referred to as the time domain.
Excessive highfrequency leakage causes ringing and can accentuate the. Since each of the rectangular pulses on the right has a fourier transform given by 2 sin ww, the convolution property tells us that the triangular function will have a fourier transform given by the square of 2 sin ww. Despite the simplicity of the antileakage fourier transform alft, and its ability to regularize pure nonuniform data sets and. Pdf seismic data regularization, which spatially transforms irregularly sampled. Pham, 2004, seismic data regularization with antileakage. Unlike resolution, leakage can be reduced by increased time domain sampling because of the increased separation between spectral replicas. This is a good point to illustrate a property of transform pairs. Fourier transform has many applications in physics and engineering such as analysis of lti systems, radar, astronomy, signal processing etc. To generalize it to higher dimensions, we proposeanareaweightingschemetoaccuratelyestimatethe fourier components. A first, standard anti leakage fourier transform method is applied to unaliased frequency components of temporaltransformed seismic data and, a second, nonstandard anti leakage fourier transform method is applied to aliased frequency components of the temporaltransformed seismic data seismic data. Sampling, for instance, produces leakage, which we call aliases of the original spectral component. In this paper, we present an algorithm, called antileakage fourier transform, for seismic data reconstruction from an irregularly sampled grid to a regular grid that overcomes these difficulties. The anti leakage fourier transform alft is a regularization method using an iterative procedure for computing the spectrum of irregularly sampled data.
Michel schonewille, andreas klaedtke, alan vigner, john brittan and tony martin present a solution to a problem with regularly sampled data found in the use of the anti alias anti leakage fourier transform alft. Im also overjoyed that my method gives identical results to its much speedier cousin, the fast fourier transform fft, because the slow fourier transform sft has a number of very significant advantages over the fft. Its discrete fourier transform is likewise a finite length sequence. An acceleration method for the antileakage parabolic radon. Consider this fourier transform pair for a small t and large t, say t 1 and t 5. Returning to case of the spectrum analyzer, the fourier transform, leakage causes a conflict of interests. To improve spectral resolution of irregular sampling signals, this paper presents an antileakage nonuniform discrete fourier transform ndft. Radio spectrum estimates using windowed data and the discrete. Note that the transform of a rectangular window is a sinc function sinxx, which has infinite width. This wellknown phenomenon is called spectral leakage. This interspectral mixing is called leakage and is distinct from the loss of resolution that occurs within the spectral replica. Figure 1 schematically describes the main areas of seis mic processing in which interpolation and regularization.
1514 606 944 439 700 1338 820 984 1303 764 364 1407 259 825 1146 15 526 172 968 842 525 611 461 834 913 1455 608 662 900 181 1280 788 1244 1049 813 1282 820 1564 1268 1092 504 937 666 448 644 390