Could the processed seismic wavelet be simpler than we think?

Geophysics ◽  
1991 ◽  
Vol 56 (5) ◽  
pp. 681-690 ◽  
Author(s):  
N. S. Neidell

J. P. Lindsey, (1988) in a clearly written short piece, opens an old question which concerns the analytic properties of seismic wavelets. This well conceived study concludes that most of the roots of a seismic wavelet as expressed by its z transform representation lie on or are very near the unit circle. The present discussion does not seek to characterize the form of all seismic wavelets, but only many if not most of those which have been processed with deconvolutions or “inversion” type operators to have reduced length, broadened bandwidth, and some desirable phase property. For such wavelets, despite the diversity by which they are obtained, remarkably simple operations having very few parameters can be extremely effective. As a case in point, constant‐phase rotations appear to carry such wavelets to zero‐phase symmetric form to a very good approximation. I start with empirical attributes which appear to characterize most processed seismic wavelets. Such wavelets tend to be of 40–100 ms duration with a smooth and unimodal amplitude spectrum of “peak” or “central” frequency between 15 and 30 Hz. The amplitude spectrum itself is further largely concentrated at frequencies between 5 and 55 Hz. A z transform root structure having essentially all of its roots only on the unit circle and on the real axis seems able to characterize all of the observed attributes rather well. This structure will be termed the band‐limiting root approximation (BLRA) and describes the attributes I seek to explain which are not as readily understood from alternative descriptions of the wavelets. Since the class of wavelets we address is obtained by a variety of means, and because the differences in character are at best subtle according to interpretive criteria, my justification is heuristic. The BLRA wavelet structure can be represented with remarkably few parameters (typically fewer than five). Of these few parameters, two relate to the frequency distribution. Such a formalism should be exceptionally useful for designing seismic techniques which seek to extract interpretive information based on properties of the wavelet.

Geophysics ◽  
1997 ◽  
Vol 62 (1) ◽  
pp. 288-290 ◽  
Author(s):  
Richard E. Duren ◽  
E. Clark Trantham

A controlled‐phase acquisition and processing methodology for our company has been described by Trantham (1994). He pointed out that it is careful attention to wavelet phase that leads to improved well ties and a more geologically accurate seismic image. In addition, we prefer zero‐phase wavelets on our seismic sections. For a given amplitude spectrum they have the simplest shape and the highest peak; further, the peak occurs at the reflection time of the event. This alignment is important since the seismic wavelet generally broadens with increasing depth with a zero‐phase wavelet remaining symmetrical about the event time. Our experience has been that a true zero‐phase section can be tied over the entire length of a synthetic trace without having to slide the synthetic trace to tie different time zones.


Geophysics ◽  
1977 ◽  
Vol 42 (7) ◽  
pp. 1369-1383 ◽  
Author(s):  
A. J. Berkhout

Detailed comparison between borehole data and seismic data has taught that, in general, conventional seismic inverse filtering is not effective enough to produce desirable deconvolution results, i.e., seismic sections with broad‐band zero‐phase wavelets. Application of conventional seismic reverse filters has the advantage that very little information is needed from the user. However, as is shown in this paper, the phase spectra of these filters may be seriously in error, even if the seismic wavelet has the minimum‐phase property. In wavelet deconvolution the phase spectrum of the filter is correct, provided a good estimate of the seismic wavelet is available. In this paper, wavelet deconvolution is compared with Wiener filtering. The main conclusions are illustrated by examples.


1969 ◽  
Vol 59 (3) ◽  
pp. 1247-1258
Author(s):  
Abou-Bakr K. Ibrahim

abstract The amplitude spectrum obtained from Haskell's matrix formulation for body waves travelling through a horizontally layered crustal model shows a sequence of minima and maxima. It is known that multiple reflections within the crustal layers produce constructive and destructive interferences, which are shown as maxima and minima in the amplitude spectrum. Analysis of the minima in the amplitude spectra, which correspond to zero phase in the phase spectra, enables us to determine the thickness of the crust, provided the ratio of wave velocity in the crust to velocity under the Moho is known.


2014 ◽  
Vol 644-650 ◽  
pp. 4334-4337
Author(s):  
Bin Wang

In this paper, frequency spectrum of infrasound signal collected was analyzed based on MATLAB. The infrasound power level of change with frequency distribution and the infrasound harmonic of each frequency component of the sound pressure amplitude can be made up judgment visually by frequency spectrum analysis and amplitude spectrum analysis, and thus we can obtain different conditions of infrasound properties scientifically and quantitatively.


Geophysics ◽  
1991 ◽  
Vol 56 (7) ◽  
pp. 1058-1060 ◽  
Author(s):  
A. E. Barnes

Robertson and Nogami (1984) have shown that the instantaneous frequency at the peak of a zero‐phase Ricker wavelet is exactly equal to that wavelet’s average Fourier spectral frequency weighted by its amplitude spectrum. Bodine (1986) gave an example which shows this is also true for constant‐phase bandpass wavelets. Here I prove that this holds for any constant‐phase wavelet. I then develop an equation expressing this quantity as a function of propagation time through an attenuating medium. A corresponding equation is derived for the amplitude of the envelope peak. Taken together, these may aid in the analysis of seismic data as suggested by Robertson and Nogami (1984), Bodine (1986), and Robertson and Fisher (1988).


Geophysics ◽  
1993 ◽  
Vol 58 (1) ◽  
pp. 39-46 ◽  
Author(s):  
Ralph W. Knapp

The suite of a wavelet is defined as being all wavelets that share a common amplitude spectrum and total energy but differ in phase spectra. Within a suite there are also classes of wavelets. A wavelet class has a common amplitude envelope and energy distribution. As such, it includes all wavelets that differ by only a constant‐angle phase shift. Of all wavelets within suite, the zero‐phase wavelet has the minimum energy envelope width; its energy is confined to minimum time dispersion. Therefore, the zero‐phase wavelet has maximum resolving power within the suite. Because a zero‐phase wavelet shares its amplitude envelope with a class of wavelets that differ by only a constant phase shift, all wavelets of the class also have maximum resolving power within the suite. The most familiar of these is the quadrature‐phase wavelet (90‐degree phase shift). Use of the complex trace results in an evaluation of the total energy, both potential and kinetic, of the wavelet signal. Assuming the wavelet signal is the output of a velocity geophone, partial energy represents only kinetic energy. Total energy better represents wavelet energy propagating through the earth. Use of partial energy (real signal only) applies a bias that favors the zero‐phase wavelets with respect to others of its class despite identical energy distribution. This bias is corrected when the wavelet envelope is used in the evaluation rather than wavelet trace amplitude. On a wiggle‐trace seismic section (amplitude display) a zero‐phase wavelet maintains a detectability advantage in the presence of noise because of a slightly greater amplitude; however, the advantage is lost in complex trace sections (energy displays) because both reflection strength and instantaneous frequency are independent of a constant phase shift in the wavelet. These sections are identical whether the wavelet is zero‐phase, quadrature‐phase or any other constant phase value, i.e., a wavelet within the zero‐phase class. (This does not imply that reflection strength sections should replace wiggle trace ones, only that they have advantages in the solution of some problems.)


Geophysics ◽  
2008 ◽  
Vol 73 (2) ◽  
pp. V11-V18 ◽  
Author(s):  
Mirko van der Baan

Phase mismatches sometimes occur between final processed sections and zero-phase synthetics based on well logs, despite best efforts for controlled-phase acquisition and processing. The latter are often based on deterministic corrections derived from field measurements and physical laws. A statistical analysis of the data can reveal whether a time-varying nonzero phase is present. This assumes that the data should be white with respect to all statistical orders after proper deterministic corrections have been applied. Kurtosis maximization by constant phase rotation is a statistical method that can reveal the phase of a seismic wavelet. It is robust enough to detect time-varying phase changes. Phase-only corrections can then be applied by means of a time-varying phase rotation. Alternatively, amplitude and phase deconvolution can be achieved using time-varying Wiener filtering. Time-varying wavelet extraction and deconvolution can also be used as a data-driven alternative to amplitude-only inverse-[Formula: see text] deconvolution.


Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. A75-A80 ◽  
Author(s):  
Mirko van der Baan ◽  
Sergey Fomel

Phase mismatches sometimes occur between final processed seismic sections and zero-phase synthetics based on well logs — despite best efforts for controlled-phase acquisition and processing. Statistical estimation of the phase of a seismic wavelet is feasible using kurtosis maximization by constant-phase rotation, even if the phase is nonstationary. We cast the phase-estimation problem into an optimization framework to improve the stability of an earlier method based on a piecewise-stationarity assumption. After estimation, we achieve space-and-time-varying zero-phasing by phase rotation.


Geophysics ◽  
1983 ◽  
Vol 48 (11) ◽  
pp. 1468-1485 ◽  
Author(s):  
Dushan B. Jovanovich ◽  
Roger D. Sumner ◽  
Sharon L. Akins‐Easterlin

Detailed lithologic interpretation of seismic sections and/or pseudo‐sonic logs generated from seismic data requires that the seismic trace can be modeled as a reflection series convolved with a zero‐phase broadband wavelet. Ghosting and marine signature deconvolution processing is a prerequisite for assuring that the seismic wavelet on a marine CDP section will be zero phase. A deterministic approach to deconvolution is centered around the concept of abandoning the purely statistical method of wavelet estimation and actually measuring the seismic wavelet. A proper signature recording for marine data is, therefore, a crucial component of deterministic deconvolution. Another important element in the deterministic deconvolution sequence is the application of a deghosting filter to remove near‐surface reflections. Proper application of a deghosting filter significantly improves the correlation between log synthetics and the seismic trace. It has been found that statistical deconvolution schemes, because of the number of statistical hypotheses required to produce a deconvolution filter, produce residual wavelets that are highly variable in character and whose average phases cover the entire phase spectrum, modulo 2π. Examples of a Gulf Coast marine line which was shot with Aquapulse™, air gun, and Maxipulse™ sources by the RV Hollis Hedberg are presented to demonstrate the differences between statistical and deterministic deconvolution processing sequences. It will be shown, using sonic logs from wells adjacent to the seismic line, that the deterministic deconvolution sections for all three sources are close to zero phase while the statistical deconvolution sections have residual average phase errors between 180 and 270 degrees. The deterministic deconvolution sections have a high degree of correlation among themselves and to the wells adjacent to the line, while the statistical deconvolution sections correlate poorly to each other and to the wells. Synthetic seismograms and their impedance logs, and the seismic sections and their corresponding pseudo‐sonic logs, are used to demonstrate how deconvolution influences lithologic interpretation. ™Western Geophysics Co.


Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 637-647 ◽  
Author(s):  
Milton J. Porsani ◽  
Bjorn Ursin

We describe a new algorithm for mixed‐phase deconvolution. It is valid only for pulses whose Z-transform has no zeros on the unit circle. That is, the amplitude spectrum cannot be zero for any frequency. Using the Z-transform of a discrete‐time signal, and assuming that the signal has α zeros inside the unit circle, the inverse of its minimum‐delay component may be estimated by solving the extended Yule‐Walker (EYW) system of equations with the lag α of the autocorrelation function (ACF) on diagonal of the coefficient matrix. This property of the solution of EYW equations is exploited to derive mixed‐phase inverse filters and their corresponding mixed‐phase pulses. For different values of α, a suite of inverse filters is generated using the same ACF. To choose the best decomposition and its corresponding mixed‐phase inverse filter, we have used the value of α which gives the maximum value of the Lp norm of the filtered signal. The optimal value of α does not seem to be very sensitive to the choice of norm as long as p > 2. In the numerical examples, we have used p = 5. The mixed‐phase deconvolution filter performs better than minimum‐phase deconvolution on the synthetic and real data examples shown.


Sign in / Sign up

Export Citation Format

Share Document