US5051963A
Expired / Lapsed
US
⏱️ 0.0 years remaining
Process for cancelling stationary sinusoidal noise in seismic data
Assignee: Mobil Oil AS
•
Filing: 1990-07-09
•
Granted: 1991-09-24
•
Expires: 2010-07-09
Abstract
Abstract An improved method for cancelling a stationary sinusoid corrupting a data trace is described. The frequency of the corrupting sinusoid is determined precisely and a Wiener-Levinson least squares filter is designed, using a reference sinusoid at precisely the determined frequency and correlated with the corrupted data trace. The filter is applied to the reference trace yielding a noise estimate which can then be subtracted from the data trace.
Inventors
- 👤 Andrew F. Linville, Jr.
- 👤 Robert A. Meek
Bibliographic Details
Application Number
—
Jurisdiction
US
Data Source
google_patents
Last Synchronized
2026-10-24 15:25
CPC Classifications
G; G01; G01V; G01V1/00; G01V1/28; G01V1/36; G01V1/364; G01V1/366; G01V2210/00; G01V2210/30; G01V2210/32
Claims Analysis
Total Claims: 20
Claim 1
Independent Claim
109 words
A method of removing stationary sinusoidal noise corrupting a seismic data trace from the trace, comprising the steps of: determining the precise frequency F o of the stationary sinusoidal noise; computing a reference sinusoidal trace R of frequency F o ; designing a filter F for shaping the reference sinusoidal trace to correspond precisely to the sinusoidal noise corrupting the data, wherein the amplitude and phase of the reference trace are adjusted to precisely match those of the sinusoid corrupting the data trace; convolving the filter F thus designed with the reference trace to yield a noise estimate R'; and subtracting the noise estimate R' from the data trace.
Claim 2
Independent Claim
109 words
A method of removing stationary sinusoidal noise corrupting a seismic data trace from the trace, comprising the steps of: determining the precise frequency F o of the stationary sinusoidal noise; computing a reference sinusoidal trace R of frequency F o ; designing a filter F for shaping the reference sinusoidal trace to correspond precisely to the sinusoidal noise corrupting the data, wherein the amplitude and phase of the reference trace are adjusted to precisely match those of the sinusoid corrupting the data trace; convolving the filter F thus designed with the reference trace to yield a noise estimate R'; and subtracting the noise estimate R' from the data trace.
Claim 3
Dependent Claim
54 words
The method of claim 1 wherein said filter design step involves the steps of: computing the autocorrelation R×R of the reference sinusoidal trace R; computing the crosscorrelation D×R of the data trace D with the reference trace R; and designing the filter F using the Wiener-Levinson algorithm, using the calculated autocorrelation and crosscorrelation values.
Claim 4
Dependent Claim
103 words
The method of claim 1 wherein said frequency F o of the stationary sinusoidal noise is determined by the following steps: (a) making an initial estimate of the frequency F o ; (b) computing a discrete Fourier transform with respect to the data trace for each of a plurality of frequencies closely spaced about F o ; (c) selecting the frequency with respect to which the value of the discrete Fourier transform is maximized as a new F o ; and (d) repeating said steps (b) and (c), iteratively reducing the spacing of the plurality of frequencies in successive performances of step (b).
Claim 5
Dependent Claim
44 words
The method of claim 3, wherein said steps (b) and (c) are repeated, iteratively reducing the spacing of the plurality of frequencies, until the spacing of the frequencies with respect to which the discrete Fourier transforms are calculated is less than a predetermined value.
Claim 6
Dependent Claim
40 words
The method of claim 1 wherein the step of determining the precise frequency F o is carried out with respect to a portion of the seismic trace which is considered to be substantially entirely sinusoidal noise from a particular source.
Claim 7
Dependent Claim
51 words
The method of claim 1 wherein prior to subtraction of the noise estimate R' from a particular trace, the amplitude of the discrete Fourier transform corresponding to frequency F o is calculated and compared to a threshold value, to determine whether removal of the noise from that particular trace is necessary.
Claim 8
Dependent Claim
14 words
The method of claim 1 wherein said filter is a Wiener-Levinson least squares filter.
Claim 9
Dependent Claim
62 words
The method of claim 7 wherein said Wiener-Levinson filter is designed by successively calculating pairs of filter operators of increasingly greater order on either side of a determined lag of maximum correlation, calculating the prediction error of each of the two filter operators, and selecting the filter operator having the lesser error as the input to the subsequent stage of filter design.
Claim 10
Dependent Claim
32 words
The method of claim 1 wherein extended portions are added to the values of the reference sinusoid prior to said filter design step, and are separated therefrom prior to said subtraction step.
Claim 11
Independent Claim
105 words
A method of removing a stationary sinusoid corrupting a data signal, comprising the steps of: determining the precise frequency F o of the stationary sinusoidal noise; computing a reference sinusoidal trace R of frequency F o ; designing a filter F for shaping the reference sinusoidal trace to correspond precisely to the sinusoidal noise corrupting the data, wherein the amplitude and phase of the reference trace are adjusted to precisely match those of the sinusoid corrupting the data trace; convolving the filter F thus designed with the reference trace to yield a noise estimate R'; and subtracting the noise estimate R' from the data trace.
Claim 12
Independent Claim
105 words
A method of removing a stationary sinusoid corrupting a data signal, comprising the steps of: determining the precise frequency F o of the stationary sinusoidal noise; computing a reference sinusoidal trace R of frequency F o ; designing a filter F for shaping the reference sinusoidal trace to correspond precisely to the sinusoidal noise corrupting the data, wherein the amplitude and phase of the reference trace are adjusted to precisely match those of the sinusoid corrupting the data trace; convolving the filter F thus designed with the reference trace to yield a noise estimate R'; and subtracting the noise estimate R' from the data trace.
Claim 13
Dependent Claim
53 words
The method of claim 10 wherein said filter design step comprises the steps of: computing the autocorrelation R×R of the reference sinusoidal trace R; computing the crosscorrelation D×R of the data trace D with the reference trace R; and designing the filter using the Wiener-Levinson algorithm, employing the calculated autocorrelation and crosscorrelation values.
Claim 14
Dependent Claim
100 words
The method of claim 10 wherein said frequency F o of the stationary sinusoidal noise is determined by the following steps: (a) making an initial estimate of the frequency F o ; (b) computing a discrete Fourier transform with respect to the data trace for each of a plurality of frequencies closely spaced about F o ; (c) selecting the frequency with respect to which the discrete Fourier transform is maximized as a new F o ; and (d) repeating said steps (b) and (c), iteratively reducing the spacing of the plurality of frequencies in successive performance of step (b).
Claim 15
Dependent Claim
44 words
The method of claim 10, wherein said steps (b) and (c) are repeated, iteratively reducing the spacing of the plurality of frequencies, until the spacing of the frequencies with respect to which the discrete Fourier transforms are calculated is less than a predetermined value.
Claim 16
Dependent Claim
40 words
The method of claim 10 wherein the step of determining the precise frequency F o is carried out with respect to a portion of the seismic trace which is considered to be substantially entirely sinusoidal noise from a particular source.
Claim 17
Dependent Claim
51 words
The method of claim 10 wherein prior to subtraction of the noise estimate R' from a particular trace, the amplitude of the discrete Fourier transform corresponding to frequency F o is calculated and compared to a threshold value, to determine whether removal of the noise from that particular trace is necessary.
Claim 18
Dependent Claim
14 words
The method of claim 10 wherein said filter is a Wiener-Levinson least squares filter.
Claim 19
Dependent Claim
59 words
The method of claim 16 wherein said Wiener-Levinson filter is designed by successively calculating pairs of filter operators of increasingly greater order on either side of a determined lag of maximum correlation, calculating the prediction error of the filter operators, and selecting the filter operator having the lesser error as the input to the subsequent stage of filter design.
Claim 20
Dependent Claim
32 words
The method of claim 10 wherein extended portions are added to the values of the reference sinusoid prior to said filter design step, and are separated therefrom prior to said subtraction step.
References Cited (Backward Citations)
USPTO Annuity & Maintenance Windows
| Annuity Stage | Window Opens | Payment Deadline | Surcharge Period Until | Estimated Fee | Status |
|---|---|---|---|---|---|
| 3.5 Years (Year 4) | 1994-09-23 | 1995-03-25 | 1995-09-24 | $2000.00 | Window Closed / Lapsed |
| 7.5 Years (Year 8) | 1998-09-23 | 1999-03-25 | 1999-09-24 | $3760.00 | Window Closed / Lapsed |
| 11.5 Years (Year 12) | 2002-09-23 | 2003-03-25 | 2003-09-24 | $7700.00 | Window Closed / Lapsed |