Patent EA-003760-B1 (EA)
Title: Method and apparatus for doppler smear correction in marine seismology measurements
Inventor: KENNETH PAUL ALLEN
IPC: G01V1/36
Filed: 2000-06-16
Abstract: 1. A method for correcting seismic data for distortion caused by the motion of a seismic source and the motion of a seismic receiver, comprising: correcting the seismic data for motion of the receiver to obtain fixed-receiver seismic data; correlating the fixed-receiver seismic data with a frequency sweep signal to obtain correlated seismic data; determining a seismic-velocity model from the correlated seismic data; and computing a dilation function for a chosen dip slice from the seismic velocity model. 2. The method of claim 1, wherein step of determining a seismic-velocity model includes: calculating root-mean-square velocities for reflected seismic waves; converting the root-mean-square velocities into interval velocities; and determining interval depths through which the reflected seismic waves traveled at corresponding interval velocities. 3. The method of claim 2, wherein step of computing a dilation function includes: determining a dilation and a travel time of seismic waves reflected from the seismic source to the seismic receiver for each diffractor in an array of diffractors assumed to be present in the seismic-velocity model; segregating the diffractors into sets sharing a common receiver dip slice; and determining a dilation versus time function from the dilation and travel times of diffractors in a selected set. 4. The method of claim 1, further comprising: computing a forward F-K transform of the correlated seismic data to produce F-K transformed data; selecting a constant dip slice of said F-K transformed data; computing an inverse F-K transform of said selected dip slice of said F-K transformed data to produce correlated seismic data for the selected dip slice; and correcting the correlated seismic data for the selected dip slice to substantially eliminate dilation as indicated by the computed dilation function from step of computing a dilation function. 5. The method of claim 4, wherein step of correcting the correlated seismic data includes: dilating a frequency sweep signal by an amount indicated from the dilation function to obtain a dilated frequency sweep signal; correlating the frequency sweep signal with the dilated frequency sweep signal to determine a phase error; applying a phase correction filter to the correlated seismic data to substantially eliminate the phase error. 6. The method of claim 4, further comprising: repeating selection a constant dip slice steps, computing an inverse F-K transform, correcting the correlated seismic data to obtain dilation-compensated seismic data for a plurality of constant dip slices; and combining the dilation-compensated seismic data for the dip slices to obtain seismic data which has been corrected for source and receiver motion. 7. The method of claim 1, further comprising: computing a forward F-K transform of the fixed-receiver seismic data to produce F-K transformed data; selecting a constant dip slice of said F-K transformed data; computing an inverse F-K transform of said constant dip slice of said F-K transformed data to produce fixed-receiver seismic data for the chosen dip slice; and correcting the fixed-receiver seismic data for the chosen dip slice to substantially eliminate dilation as indicated by the computed dilation function from step of computing a dilation function. 8. The method of claim 7, wherein step of correcting the correlated seismic data includes: dilating a frequency sweep signal by an amount indicated from the dilation function to obtain a dilated frequency signal; and correlating the fixed-receiver seismic data with the dilated frequency sweep signal to produce dilation-compensated correlated seismic data. 9. The method of claim 7, further comprising: repeating selection a constant dip slice steps, computing an inverse F-K transform, correcting the correlated seismic data to obtain dilation-compensated seismic data for a plurality of constant dip slices; and combining the dilation-compensated seismic data for the dip slices to obtain seismic data which has been corrected for source and receiver motion. 10. The method of claim 7 wherein said step of selecting a constant dip slice of said F-K transformed data includes selecting F-K transformed data corresponding to a predetermined range of dip angles. 11. A seismic measurement system, comprising: a user input device; a display device; a data storage unit; a processing unit coupled to the user input device, the display device, and the data storage unit, wherein the processing unit is configured to create a seismic velocity model from measured seismic data, and is configured to use the seismic velocity model to calculate a dilation function indicative of distortion caused by motion of a seismic source. 12. The seismic measurement system of claim 11, wherein said processing unit is further configured to correct the measured seismic data for the distortion caused by motion of the seismic source by computing an F-K transform of the measured seismic data, partitioning the transformed seismic data into dip slices, inverse transforming the individual dip slices, correcting the inverse transformed dip slice data, and adding together corresponding time values of the corrected dip slice data. 13. The seismic measurement system of claim 12, wherein the processing unit corrects the inverse transformed dip slice data by applying a series of filters, wherein the filters are determined by dilating a frequency sweep signal as indicated by the dilation function. 14. The seismic measurement system of claim 13, wherein the dilation function is divided into time segments and a filter is determined for each time segment and applied to a corresponding time segment of the inverse-transformed dip slice data. 15. A seismic system, comprising: at least one seismic source configured to transmit a frequency sweep signal; at least one seismic receiver configured to provide measured seismic signals; and a seismic measurement and processing system which includes: a user input device; a display device; a data storage unit; a source and receiver interface unit coupled to said seismic source and seismic receiver; a processing unit coupled to the user input device, the display device, and the data storage unit, wherein the processing unit is configured to create a seismic velocity model from measured seismic data, and is configured to use the seismic velocity model to calculate a dilation function indicative of distortion caused by motion of a seismic source. 16. The seismic measurement system of claim 15, wherein said processing unit is further configured to correct the measured seismic data for the distortion caused by motion of the seismic source by computing an F-K transform of the measured seismic data, partitioning the transformed seismic data into dip slices, inverse transforming the individual dip slices, correcting the inverse transformed dip slice data, and adding together corresponding time values of the corrected dip slice data. 17. The seismic measurement system of claim 16, wherein the processing unit corrects the inverse transformed dip slice data by applying a series of filters, wherein the filters are determined by dilating a frequency sweep signal as indicated by the dilation function. 18. The seismic measurement system of claim 17, wherein the dilation function is divided into time segments and a filter is determined for each time segment and applied to a corresponding time segment of the inverse-transformed dip slice data.