US4884247A
Expired / Lapsed
US
⏱️ 0.0 years remaining
Method of processing geophysical data to compensate for earth filter attenuation
Assignee: ExxonMobil Oil Corp
•
Filing: 1987-10-13
•
Granted: 1989-11-28
•
Expires: 2007-10-13
Abstract
Abstract Method to process seismic data to yield useful information regarding the subsurface area explored. Generated seismic energy travels into a subsurface formation where it is distorted and attenuated. The distorted and attenuated seismic pulse is detected by a geophone or other seismic receiver and recorded. The received seismic traces are first compensated for the attenuation effects of the earth filter and the compensated seismic traces are then used to produce a close approximation of a zero-phase, band-limited output wavelet of the reflectivity of the geophysical formation under exploration.
Inventors
- 👤 Mohamed T. Hadidi
- 👤 Andrejs Jurkevics
- 👤 Andrew F. Linville, Jr.
Bibliographic Details
Application Number
—
Jurisdiction
US
Data Source
google_patents
Last Synchronized
2026-10-29 07:33
CPC Classifications
G; G01; G01V; G01V1/00; G01V1/28; G01V1/36; G01V1/364; G01V1/366; G01V2210/00; G01V2210/50; G01V2210/58
Claims Analysis
Total Claims: 41
Claim 1
Independent Claim
89 words
A method of compensating for the effects of intrinsic and apparent attenuation of a seismic signal transmitted through a geophysical formation comprising the steps of: determining an antisymmetric component of the cepstrum of an inverse earth filter; determining a symmetric component of the cepstrum of said inverse earth filter modified to effect gain limiting; combining said modified symmetric component and said antisymmetric compenasted to determine said inverse earth filter of said geophysical formation; and applying said inverse earth filter to said seismic signal to determine a compensated seismic signal.
Claim 2
Independent Claim
89 words
A method of compensating for the effects of intrinsic and apparent attenuation of a seismic signal transmitted through a geophysical formation comprising the steps of: determining an antisymmetric component of the cepstrum of an inverse earth filter; determining a symmetric component of the cepstrum of said inverse earth filter modified to effect gain limiting; combining said modified symmetric component and said antisymmetric compenasted to determine said inverse earth filter of said geophysical formation; and applying said inverse earth filter to said seismic signal to determine a compensated seismic signal.
Claim 3
Dependent Claim
68 words
The method of claim 1 wherein the step of determining said modified symmetric component of the cepstrum of the inverse earth filter further comprises the steps of: determining a real component of the logarithm of the frequency domain inverse earth filter; modifying said real component to stabilize the inverse earth filter; and taking the inverse Fourier transform of said modified real component to determine said modified symmetric component.
Claim 4
Dependent Claim
55 words
The method of claim 1 wherein the step of determining said antisymmetric component of said cepstrum of the inverse earth filter further comprises the steps of: determining an imaginary component of the logarithm of the frequency domain inverse earth filter; and taking the inverse Fourier transform of said imaginary component to determine said antisymmetric component.
Claim 5
Dependent Claim
15 words
The method of claim 1 further comprising the step of deconvolving said compensated seismic signal.
Claim 6
Dependent Claim
29 words
The method of claim 4 wherein the step of deconvolving said compensated seismic signal further comprises the step of deconvolving said compensated seismic signal according to spiking deconvolution methods.
Claim 7
Dependent Claim
29 words
The method of claim 4 wherein the step of deconvolving said compensated seismic signal further comprises the step of deconvolving said compensated seismic signal according to flatiron deconvolution methods.
Claim 8
Dependent Claim
74 words
The method of claim 4 further comprising the steps of: whitening the effective reflectivity of the amplitude spectrum of said compensated seismic signal; reducing the effects of reflectivity leakage on said amplitude spectrum of said compensated seismic signal; smoothing said amplitude spectrum of said compensated seismic signal to eliminate rapidly varying components of said amplitude spectrum; and compensating said amplitude spectrum of said compensated seismic signal to approximate that of an analog seismic wavelet.
Claim 9
Dependent Claim
174 words
The method according to claim 4 further comprising the steps of: whitening effective reflectivity of said amplitude spectrum of said seismic signal; reducing the effects of reflectivity leakage on said amplitude spectrum of said seismic signal, said step of reducing the effects of reflectivity leakage further including the steps of transforming said amplitude spectrum from the frequency domain into the quefrency domain, selecting a first cepstral window length, applying said first cepstral window to said quefrency domain amplitude spectrum, and transforming said amplitude spectrum from said quefrency domain to said frequency domain; compensating said amplitude spectrum of said seismic signal 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 spectrum from the frequency domain into the quefrency domain, selecting a second cepstral window length, applying said second cepstral window to said quefrency domain amplitude spectrum, and transforming said amplitude spectrum from said quefrency domain to said time domain.
Claim 10
Dependent Claim
64 words
The method according to claim 1 wherein the step of applying said inverse earth filter to said seismic signal to determine said compensated seismic signal further comprises the steps of: repeatedly applying said combination to said seismic signal to produce a series of traces; and summing said series of traces produced by said repeated applications and said seismic signal to determine said compensated signal.
Claim 11
Dependent Claim
41 words
The method according to claim 9 wherein the step of repeatedly applying said combination to said seismic signal to produce a series of traces further comprises the step of convolving the Taylor series expansion of said combination with said seismic signal.
Claim 12
Independent Claim
111 words
A method of compensating for the effects of intrinsic and apparent attenuation of a series of seismic signals transmitted through a geophysical formation comprising the steps of: selecting a first seismic signal; determining an antisymmetric component of the cepstrum of the inverse earth filter; determining a symmetric component of the cepstrum of the inverse earth filter modified to effect gain limiting; combining said symmetric and antisymmetric components; repeatedly applying said combination to said first seismic signal; combining partial results produced by said repeated applications to determine a first compensated seismic signal; selecting a next seismic signal; and repeating said selecting, applying, and combining steps until all seismic signals have been compensated.
Claim 13
Independent Claim
111 words
A method of compensating for the effects of intrinsic and apparent attenuation of a series of seismic signals transmitted through a geophysical formation comprising the steps of: selecting a first seismic signal; determining an antisymmetric component of the cepstrum of the inverse earth filter; determining a symmetric component of the cepstrum of the inverse earth filter modified to effect gain limiting; combining said symmetric and antisymmetric components; repeatedly applying said combination to said first seismic signal; combining partial results produced by said repeated applications to determine a first compensated seismic signal; selecting a next seismic signal; and repeating said selecting, applying, and combining steps until all seismic signals have been compensated.
Claim 14
Dependent Claim
68 words
The method of claim 11 wherein the step of determining said modified symmetric component of the cepstrum of the inverse earth filter further comprises the steps of: determining a real component of the logarithm of the frequency domain inverse earth filter; modifying said real component to stabilize the inverse earth filter; and taking the inverse Fourier transform of said modified real component to determine said modified symmetric component.
Claim 15
Dependent Claim
52 words
The method of claim 11 wherein the step of determining said antisymmetric component of said inverse earth filter further comprises the steps of: determining an imaginary component of the logarithm of the frequency domain inverse earth filter; and taking the inverse Fourier transform of said imaginary component to determine said antisymmetric component.
Claim 16
Dependent Claim
15 words
The method of claim 11 further comprising the step of deconvolving said compensated seismic signal.
Claim 17
Dependent Claim
29 words
The method of claim 11 wherein the step of deconvolving said compensated seismic signal further comprises the step of deconvolving said compensated seismic signal according to spiking deconvolution methods.
Claim 18
Dependent Claim
29 words
The method of claim 11 wherein the step of deconvolving said compensated seismic signal further comprises the step of deconvolving said compensated seismic signal according to flatiron deconvolution methods.
Claim 19
Dependent Claim
74 words
The method of claim 11 further comprising the steps of: whitening the effective reflectivity of the amplitude spectrum of said compensated seismic signal; reducing the effects of reflectivity leakage on said amplitude spectrum of said compensated seismic signal; smoothing said amplitude spectrum of said compensated seismic signal to eliminate rapidly varying components of said amplitude spectrum; and compensating said amplitude spectrum of said compensated seismic signal to approximate that of an analog seismic wavelet.
Claim 20
Dependent Claim
36 words
The method according to claim 11 wherein the step of repeatedly applying said combination to said first seismic signal further comprises the step of convolving the Taylor series expansion of said combination with said seismic signal.
Claim 21
Dependent Claim
174 words
The method according to claim 11 further comprising the steps of: whitening effective reflectivity of said amplitude spectrum of said seismic signal; reducing the effects of reflectivity leakage on said amplitude spectrum of said seismic signal, said step of reducing the effects of reflectivity leakage further including the steps of transforming said amplitude spectrum from the frequency domain into the quefrency domain, selecting a first cepstral window length, applying said first cepstral window to said quefrency domain amplitude spectrum, and transforming said amplitude spectrum from said quefrency domain to said frequency domain; compensating said amplitude spectrum of said seismic signal to approximate that of an analog seismic wavelet; and smoothing said compensated amplitude spectrum to elimate rapidly varying components, said step of smoothing said compensated amplitude spectrum further including the steps of transforming said amplitude spectrum from the frequency domain into the quefrency domain, selecting a second cepstral window length, applying said second cepstral window to said quefrency domain amplitude spectrum, and transforming said amplitude spectrum from said quefrency domain to said time domain.
Claim 22
Independent Claim
155 words
A method of compensating the amplitude spectrum of a seismic trace for the effects of earth filtering, nonwhite reflectivity of the subsurface formation being explored and sampling, comprising: determining an antisymmetric component of the cepstrum of an inverse earth filter; determining a symmetric component of the cepstrum of said inverse earth filter modified to effect gain limiting; combining said modified symmetric component and said antisymmetric component to determine said inverse earth filter of said geophysical formation; applying said inverse earth filter to said seismic signal to determine a compensated seismic signal; whitening effective reflectivity of the amplitude spectrum of said compensated seismic signal; reducing the effects of reflectivity leakage on said amplitude spectrum of said compensated seismic signal; compensating said amplitude spectrum of said seismic signal to approximate that of an analog seismic wavelet; and smoothing said compensated amplitude spectrum of said compensated seismic signal to eliminate rapidly varying components of said compensated amplitude spectrum.
Claim 23
Independent Claim
155 words
A method of compensating the amplitude spectrum of a seismic trace for the effects of earth filtering, nonwhite reflectivity of the subsurface formation being explored and sampling, comprising: determining an antisymmetric component of the cepstrum of an inverse earth filter; determining a symmetric component of the cepstrum of said inverse earth filter modified to effect gain limiting; combining said modified symmetric component and said antisymmetric component to determine said inverse earth filter of said geophysical formation; applying said inverse earth filter to said seismic signal to determine a compensated seismic signal; whitening effective reflectivity of the amplitude spectrum of said compensated seismic signal; reducing the effects of reflectivity leakage on said amplitude spectrum of said compensated seismic signal; compensating said amplitude spectrum of said seismic signal to approximate that of an analog seismic wavelet; and smoothing said compensated amplitude spectrum of said compensated seismic signal to eliminate rapidly varying components of said compensated amplitude spectrum.
Claim 24
Dependent Claim
106 words
The method of claim 20 wherein the step of whitening the effective reflectivity of said amplitude spectrum of said seismic trace further comprises the steps of: selecting a first seismic trace; determining an autocorrelation function for said first seismic trace; determining an amplitude spectrum for said first seismic trace from said autocorrelation function; stacking said amplitude spectrum onto said stack; selecting a next seismic trace; repeating the steps of determining the autocorrelation function, determining the amplitude spectrum, and stacking said amplitude spectrum for said next seismic trace until all seismic traces have been stacked; and averaging said stacked amplitude spectrum to yield an average amplitude spectrum.
Claim 25
Dependent Claim
75 words
The method of claim 20 wherein the step of compensating said amplitude spectrum to approximate an analog seismic wavelet further comprises the steps of: amortizing low frequencies of said amplitude spectrum to remove the effect of low frequency noise such as ground roll from the estimated 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.
Claim 26
Dependent Claim
65 words
The method of claim 22 wherein the step of amortizing low frequencies of said amplitude spectrum to eliminate low frequency ground roll noise further comprises the steps of: selecting a low frequency value at which to commence low end amortization; and applying the slope of the low cut filter to the portion of said amplitude spectrum having a frequency less than said low frequency value.
Claim 27
Dependent Claim
107 words
The method of claim 22 wherein the step of amortizing high frequencies of said amplitude spectrum to remove the effect of high frequency noise from the estimate of the wavelet amplitude spectrum further comprises the steps of: selecting a first and second high frequency value to commence high end amortization; applying said characteristics of the earth filter to the portion of the ampltidue spectrum between said first and said second high frequency values; and applying the slope of the high cut filter and the slope of the earth filter to the portion of said amplitude spectrum between said second high frequency value and an upper limit frequency.
Claim 28
Dependent Claim
72 words
The method according to claim 20 wherein the step of reducing the effects of reflectivity leakage on said amplitude spectrum of said compensated seismic signal further comprises the steps of transforming said amplitude spectrum from the frequency domain into the quefrency domain, selecting a first cepstral window length, applying said first cepstral window to said quefrency domain amplitude spectrum, and transforming said amplitude spectrum from said quefrency domain to said frequency domain.
Claim 29
Dependent Claim
78 words
The method according to claim 25 where in the step of smoothing said compensated amplitude spectrum of said compensated seismic signal to eliminate rapidly varying components of said compensated amplitude spectrum further comprises the steps of transforming said amplitude spectrum from the frequency domain into the quefrency domain, selecting a second cepstral window length, applying said second cepstral window to said quefrency domain amplitude spectrum, and transforming said amplitude spectrum from said quefrency domain to said time domain.
Claim 30
Independent Claim
218 words
A method of processing geophysical data to make it more useful in interpreting the subsurface characteristics of the earth comprising: acquiring at least one seismic trace containing information related to the subsurface characteristics of the earth; determining an antisymmetric component of the cepstrum of an inverse earth filter; determining a symmetric component of the cepstrum of said inverse earth filter; combining said modified symmetric component and said antisymmetric component to determine said inverse earth filter of said geophysical formation; applying said inverse earth filter to said at least one seismic trace to determine at least one compensated seismic trace; determining an amplitude spectrum for said at least one compensated trace; whitening effective reflectivity of each of said at least one compensated 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 zerophase, low sidelobe output wavelet by said two sided, least squares filter.
Claim 31
Independent Claim
218 words
A method of processing geophysical data to make it more useful in interpreting the subsurface characteristics of the earth comprising: acquiring at least one seismic trace containing information related to the subsurface characteristics of the earth; determining an antisymmetric component of the cepstrum of an inverse earth filter; determining a symmetric component of the cepstrum of said inverse earth filter; combining said modified symmetric component and said antisymmetric component to determine said inverse earth filter of said geophysical formation; applying said inverse earth filter to said at least one seismic trace to determine at least one compensated seismic trace; determining an amplitude spectrum for said at least one compensated trace; whitening effective reflectivity of each of said at least one compensated 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 zerophase, low sidelobe output wavelet by said two sided, least squares filter.
Claim 32
Dependent Claim
69 words
The method of claim 27 further comprising the steps of: determining the autocorrelation function for each of said at least one seismic trace; determining an amplitude spectrum for each of said at least one seismic trace from said autocorrelation function; stacking said amplitude spectrums for said at least one seismic trace; and determining an average amplitude spectrum for said at least one seismic trace from said stacked amplitude spectrum.
Claim 33
Dependent Claim
79 words
The method of claim 27 wherein the step of compensating said amplitude spectrum of said at least one seismic trace to approximate an analog seismic wavelet further comprises: amortizing low frequencies of said amplitude spectrum to remove the effect of low frequency noise such as ground roll from the estimate of the wavelet amplitude spectrums; and amortizing high frequencies of said amplitude spectrum to remove the effect of high frequency noise from the estimate of the wavelet amplitude spectrum.
Claim 34
Dependent Claim
71 words
The method of claim 29 wherein the step of amortizing said low frequencies of said amplitude spectrum to remove the effect of low frequency noise such as ground roll further comprises the steps of: selecting a low frequency value at which to commence low end amortization; and applying the slope of the low cut filter to the portion of said amplitude spectrum having a frequency less than said low frequency value.
Claim 35
Dependent Claim
107 words
The method of claim 29 wherein the step of amortizing the high frequencies of said amplitude spectrum remove the effect of high frequency noise from the estimate of the wavelet amplitude spectrum further comprises the steps of: selecting a first and second high frequency value to commence high end amortization; applying said characteristics of the earth filter to the portion of the amplitude spectrum between said first and said second high frequency values; and applying the slope of the high cut filter and the slope of the earth filter to the portion of said amplitude spectrum between said second high frequency value and an upper limit frequency.
Claim 36
Dependent Claim
23 words
The method of claim 27 further comprising the step of selecting a portion of said at least one compensated seismic trace for interpretation.
Claim 37
Dependent Claim
33 words
The method of claim 27 wherein the step of compensating said amplitude spectrum to approximate an analog wavelet further comprises the step of reducing the effects of reflectivity leakage on said amplitude spectrum.
Claim 38
Dependent Claim
66 words
The method of claim 33 wherein the step of reducing the effects of reflectivity leakage on said amplitude spectrum further comprises the steps of: transforming said amplitude spectrum from the frequency domain into the quefrency domain; selecting a first cepstral window length; applying said first cepstral window to said quefrency domain amplitude spectrum; and transforming said amplitude spectrums from said quefrency domain to said time domain.
Claim 39
Dependent Claim
23 words
The method of claim 33 further comprising the step of smoothing said amplitude spectrum to eliminate rapidly varying components of said amplitude spectrum.
Claim 40
Dependent Claim
97 words
The method of claim 35 wherein the step of smoothing said amplitude spectrum to eliminate rapidly varying components of said amplitude spectrum further comprises the steps of transforming said amplitude spectrums from the frequency domain into the quefrency domain; doubling all valves in the first half period of said transformed amplitude spectrum greater than zero; eliminating all valves in the second half period of said transformed amplitude spectrum; selecting a second cepstral window length; applying said second cepstral window to said transformed amplitude spectrum; and transforming said amplitude spectrum from said quefrency domain to said time domain.
Claim 41
Dependent Claim
69 words
The method of claim 35 wherein the step of smoothing said amplitude spectrum to eliminate rapidly varying components of said amplitude spectrum further comprises the steps of: transforming said amplitude spectrum from the frequency domain into the quefrency domain; selecting a second cepstral window length; applying said second cepstral window to said quefrency domain amplitude spectrum; and transforming said amplitude spectrum from said quefrency domain to said time domain.
References Cited (Backward Citations)
US3396365A
Cited
US3581274A
Cited
US3634668A
Cited
US3662108A
Cited
US4052889A
Cited
US4164021A
Cited
US4209672A
Cited
US4213183A
Cited
US4310721A
Cited
US4545039A
Cited
US4545250A
Cited
US4608673A
Cited
US4655228A
Cited
US4780859A
Cited
USPTO Annuity & Maintenance Windows
| Annuity Stage | Window Opens | Payment Deadline | Surcharge Period Until | Estimated Fee | Status |
|---|---|---|---|---|---|
| 3.5 Years (Year 4) | 1992-11-27 | 1993-05-29 | 1993-11-28 | $2000.00 | Window Closed / Lapsed |
| 7.5 Years (Year 8) | 1996-11-27 | 1997-05-29 | 1997-11-28 | $3760.00 | Window Closed / Lapsed |
| 11.5 Years (Year 12) | 2000-11-27 | 2001-05-29 | 2001-11-28 | $7700.00 | Window Closed / Lapsed |