Method of interpreting seismic data
First Claim
1. In seismic exploration of a subsurface formation wherein seismic source are activited at spaced locations along a line of exploration and the amplitude of reflections are recorded as a function of time along said lines to produce a series of seismic traces, a method of compensating the amplitude spectrum of a selected one of said seismic traces for the effects of non-white reflectivity of said subsurface formation comprising:
- whitening effective reflectivity of said amplitude spectrum of said selected seismic trace;
reducing the effects of reflectivity leakage on said amplitude spectrum of said selected seismic trace, said step of reducing the effects of reflectivity leakage further including the steps of transforming said amplitude spectrum from the frequency domain into the frequency domain, selecting a first cepstral window length, applying first cepstral window to said frequency domain amplitude spectrum, and transforming said amplitude spectrum from said frequency domain to said frequency domain;
compensating said amplitude spectrum of said selected seismic trace to approximate that of an analog seismic wavelet; and
smoothing said compensated amplitude spectrum to eliminate rapidly varying components, said step of smoothing said compensated amplitude spectrum further including the steps of transforming said amplitude from the frequency domain into the frequency domain, selecting a second cepstral window length, applying said second cepstral window to said frequency domain amplitude spectrum, and transforming said amplitude spectrum from said frequency domain to said time domain.
2 Assignments
0 Petitions
Accused Products
Abstract
Method to process seismic data to yield useful information regarding the subsurface area explored. Seismic data which is acquired from the seismic exploration is autocorrelated and the amplitude spectrums estimated, stacked and averaged to reduce the effects of nonwhite reflectivity upon the data. The amplitude spectrum is cepstral windowed a first time to reduce the effects of reflectivity leakage on the amplitude spectrum, amortized at both low and high frequencies so that the amplitude spectrum will approximate the analog wavelet and cepstral windowed a second time to eliminate any remaining rapidly varying components of the amplitude spectrum. The resultant amplitude spectrum is then used to estimate the phase spectrum for the wavelet and a two-sided least squares filter is then designed to shape the estimated seismic wavelet into a zero-phase, low sidelobe output wavelet.
-
Citations
23 Claims
-
1. In seismic exploration of a subsurface formation wherein seismic source are activited at spaced locations along a line of exploration and the amplitude of reflections are recorded as a function of time along said lines to produce a series of seismic traces, a method of compensating the amplitude spectrum of a selected one of said seismic traces for the effects of non-white reflectivity of said subsurface formation comprising:
-
whitening effective reflectivity of said amplitude spectrum of said selected seismic trace; reducing the effects of reflectivity leakage on said amplitude spectrum of said selected seismic trace, said step of reducing the effects of reflectivity leakage further including the steps of transforming said amplitude spectrum from the frequency domain into the frequency domain, selecting a first cepstral window length, applying first cepstral window to said frequency domain amplitude spectrum, and transforming said amplitude spectrum from said frequency domain to said frequency domain; compensating said amplitude spectrum of said selected seismic trace to approximate that of an analog seismic wavelet; and smoothing said compensated amplitude spectrum to eliminate rapidly varying components, said step of smoothing said compensated amplitude spectrum further including the steps of transforming said amplitude from the frequency domain into the frequency domain, selecting a second cepstral window length, applying said second cepstral window to said frequency domain amplitude spectrum, and transforming said amplitude spectrum from said frequency domain to said time domain. - View Dependent Claims (2, 3, 4, 5, 6, 7)
-
-
8. A method of processing geophysical data to make it more useful in interpreting subsurface characterstics of the earth comprising:
-
acquiring at least one seismic trace containing information related to the subsurface characteristics of the earth; determining an amplitude spectrum for each of said at least one seismic trace; whitening effective reflectivity of each of said at least one seismic trace; compensating said amplitude spectrum of said at least one seismic trace to approximate an analog seismic wavelet; estimating a phase spectrum for said seismic wavelet from said compensated seismic wavelet; determining a two sided, least squares filter from said estimated phase spectrum; and determining an output seismic wavelet from said compensated amplitude spectrum and said estimated phase spectrum by applying said two sided, least squares filter to said compensated amplitude spectrum in the time domain; wherein said seismic wavelet is shaped into a zero-phase, low side-lobe output wavelet. - View Dependent Claims (9, 10, 11, 12, 13, 14, 15, 16, 17, 18)
-
-
19. In seismic exploration wherein seismic sources are activated at spaced locations along a line of exploration and where the amplitude of reflections are recorded as a function of time along said line to produce a series of seismic traces, a method of interpreting said seismic traces to provide information related to subsurface characteristics of the earth comprising the steps of:
-
determining an amplitude spectrum for each of said series of seismic traces; whitening effective reflectivity of said series of seismic traces; reducing the effects of relectivity leakage on said amplitude spectrums of said series of seismic traces, said step of reducing the effects of reflectivity leakage further including the stpes of transforming said amplitude spectrums from the frequency domain into the frequency domain, selecting a first cepstral window length, applying said first cepstral window to said frequency domain amplitude spectrums, and transforming said amplitude spectrums from said frequency domain to said frequency domain; compensating said amplitude spectrums of said series of seismic traces to approximate an analog seismic wavelet, said step of compensating said amplitude spectrum of said series of seismic traces to approximate an analog seismic wavelet further including the steps of amortizing low frequencies of said amplitude spectrums to remove the effect of low frequency noise such as ground roll from the estimate of the wavelet amplitude spectrum and amortizing high frequencies of said amplitude spectrum to remove the effect of high frequency noise from the estimate of the wavelet amplitude spectrum; smoothing said compensated amplitude spectrums to eliminate rapidly varying components, said step of smoothing said compensated amplitude spectrums further including the steps of transforming said amplitude spectrums from the frequency domain into the frequency domain, selecting a second cepstral window length, applying said second cepstral window to said frequency domain amplitude spectrums, and transforming said amplitude spectrums from said frequency domain to said time domain; estimating a phase spectrum for said seismic wavelet from said smoothed compensated seismic wavelet; determining a two sided, least squares filter from said estimated phase spectrum; and determining an output seismic wavelet from said smoothed compensatedd amplitude spectrum said estimated phase spectrum by applying said two sided, least squares filter to said smoothed compensated amplitude spectrum in the time domain wherein said seismic wavelet is shaped into a zero-phase, low side-lobe wavelet. - View Dependent Claims (20, 21, 22, 23)
-
Specification