Patent EA-002508-B1 (EA)
Title: Adaptive seismic noise and interference attenuation method
Inventor: OZBEK ALI
IPC: G01V1/36
Filed: 1999-05-18
Abstract: 1. Method for filtering noise from discrete noisy seismic signals, said method comprising the steps of: receiving signals representing seismic energy reflected and/or refracted from an earth structure using a plurality of seismic receivers; and filtering received signals using an at least partially adaptive filter that attenuates signals that lack predetermined propagation characteristics: said filtering step comprising the steps of: defining at least two independent sets of conditions with a first set defining a desired response, and a second set, defining the propagation characteristics of signals to be preserved wherein the first set is not affected by the second set; and adapting filter coefficients of said filter subject to said at least two independent sets of conditions so as to optimize the filter output for signals that lack the propagation characteristics of signals to be preserved. 2. The method of claim 1 wherein the step of adapting comprises optimizing a cost function representing the output power of the filter and the white noise gain of the filter. 3. The method of claim 2 wherein the cost function represents a weighted sum of the output power of the filter and the white noise gain of the filter. 4. The method of claim 3 wherein the relative weight of the output power and the "white noise gain" in the cost function is adjustable. 5. The method of claim 1 wherein the propagation characteristics of a signal is determined by its temporal and spatial spectral content. 6. The method of claim 1 wherein the propagation characteristics of a signal is determined by its wavenumber and frequency. 7. The method of claim 1 wherein the propagation characteristics of a signal is defined as a region or as regions in the frequency-wavenumber domain. 8. The method of claim 1 wherein the propagation characteristics of a signal includes a plurality directions in the frequency-wavenumber domain. 9. The method of claim 1 wherein the conditions define a constraint matrix comprising at least two mutually orthogonal subspaces with a first subspace defining the desired response of the filter and a second subspace defining the propagation characteristics of the signals to be preserved with said matrix being applied during the adaptation process of filter coefficients. 10. The method of claim 1 wherein the conditions define a constraint matrix using a quiescent response vector and principle eigenvectors of a correlation matrix made of steering vectors defining at least one region of the frequency-wavenumber domain wherein the signals are preserved. 11. The method of claim 1 wherein the conditions are defined such that signals with pre-determined propagation characteristics are suppressed by the filter. 12. The method of claim 1, wherein said filtering step further comprises defining further sets of conditions with one set of conditions forcing the desired response of the filter to zero and one set of conditions defining the propagation characteristics of signals to be suppressed by the filter. 13. The method of claim 1 wherein the filter comprises M temporally local filters, said filters forming a filter bank, and M being a number equal to or larger than two. 14. The method of claim 13 including the step of multiplying M filtered estimates with temporal window functions. 15. The method of claim 14 wherein the temporal window functions are characterized by the requirement that only adjoining windows overlap. 16. The method of claim 14 wherein the cost function is minimized using the approximation that the sum, weighted by window functions, of the output of adjacent filters of the M filters is equal when applied to the same signal in time regions where said window functions overlap. 17. The method of claim 1 comprising the further step of dividing filter coefficients into a fixed part and a adaptive part. 18. The method of claim 1 comprising the further step of splitting the signals into at least two frequency bands, filtering at least one of said at least two frequency band, and, before further processing, recombining said frequency bands to regain noise attenuated seismic signals. 19. The method of claim 1 wherein the discrete noisy signals used as input are recordings from individual seismic sensors prior to any group forming techniques. 20. The method of claim 1 comprising the further steps of positioning at least one seismic source and a plurality of seismic receivers at appropriate positions in a land, marine or transitional zone environment; activating said at least one source to propagate energy through subterranean formations, and using said receivers to measure energy as single sensor recordings, converting said single sensor recordings into discrete seismic signals; and transmitting said signals as input to the filter.