Claims
- 1. A method of seismic analysis comprising acquiring seismic data, applying a temporal Fourier transform to the seismic data to yield transformed data, applying a spatial wavenumber estimation based on a parametric algorithm to the transformed data, and generating a wavenumber spectrum from the parametric algorithm.
- 2. The method of claim 1, wherein the second applying step comprises selecting a frequency span around a frequency F.
- 3. The method of claim 1, wherein the second applying step comprises generating a spatio-spectral correlation matrix.
- 4. The method of claim 3, wherein the second applying step further comprises decomposing the spatio-spectral correlation matrix to noise and signal subspace.
- 5. The method of claim 4, wherein the decomposing step is performed using an eigenvalue decomposition method.
- 6. The method of claim 5, where there is a magnitude threshold for eigenvalues.
- 7. The method of claim 1, wherein the parametric algorithm is a MUSIC algorithm,
- 8. The method of claim 1, wherein the parametric algorithm is Minimum Variance algorithm, Eigenvector algorithm, Maximum Likelihood algorithm, Pisarenco algorithm, an ARMA, AR or MA algorithm or Maximum Entropy algorithm.
- 9. The method of claim 1, wherein the second applying step comprises selecting a frequency span around a frequency F, generating a spatio-spectral correlation matrix, decomposing the spatio-spectral correlation matrix to noise and signal subspace, and wherein the parametric algorithm is a MUSIC algorithm, Minimum Variance algorithm, Eigenvector algorithm, Maximum Likelihood algorithm, Pisarenco algorithm, an ARMA, AR or MA algorithm or Maximum Entropy algorithm.
- 10. The method of claim 1, wherein the parametric algorithm is applied to data from equidistant spatial sampling.
- 11. The method of claim 1, wherein the parametric algorithm is applied to data from sparse spatial sampling.
- 12. The method of claim 11, wherein parameters passed to the parametric algorithm are selected from the group of order of Fourier transform, sampling frequency, distance between receivers in an x direction, distance between receivers in a y direction, number_of_sources, or f_span.
- 13. The method of claim 12, wherein the f_span passed as a parameter to the parametric algorithm determines a temporal frequency domain resolving the wavenumber spectrum.
- 14. The method of claim 12, wherein number_of_sources is passed to the parametric algorithm.
- 15. The method of claim 14, wherein the number_of_sources is not greater than the number of signals.
- 16. The method of claim 1, wherein seismic data is acquired from a one dimensional array of sensors.
- 17. The method of claim 16, wherein the data is fed into the parametric algorithm as an M×L matrix, wherein M is the number of signals in a one dimensional array and L is the number of samples per sensor.
- 18. The method of claim 16, wherein the generated wavenumber spectrum is a one-dimensional fk spectrum for each frequency span.
- 19. The method claim 18, comprising a further step of generating an fkk spectrum 3D volume.
- 20. The method of claim 1, wherein seismic data is acquired from a two dimensional array of sensors.
- 21. The method of claim 19, wherein the data is fed into the parametric algorithm as an MN×L matrix, wherein MN is the number of signals in a two dimensional array and L is the number of samples per sensor.
- 22. The method of claim 19, wherein the generated wavenumber spectrum is a two-dimensional fkk spectrum for each frequency span.
- 23. The method claim 21, comprising a further step of generating an fkk spectrum 3D volume.
- 24. The method of claim 1, wherein seismic data is acquired from a three dimensional array of sensors.
- 25. The method of claim 24, wherein the data is fed into the parametric algorithm as an MNO×L matrix, wherein MNO is the number of signals in a three dimensional array and L is the number of samples per sensor.
- 26. The method claim 24, wherein the generated wavenumber spectrum is a three-dimensional fkkk spectrum for each frequency span.
- 27. The method claim 26, comprising a further step of generating an fkkk spectrum 4D image.
- 28. The method of claim 1, comprising a further step of performing velocity analysis.
CROSS-REFERENCE TO RELATED APPLICATIONS
[0001] This application claims benefit of U.S. provisional patent application serial No. 60/368,122, filed Mar. 27, 2002, which is herein incorporated by reference.
Provisional Applications (1)
|
Number |
Date |
Country |
|
60368122 |
Mar 2002 |
US |