Patent ID: 8811116

Claim:
A method of performing a geophysical survey, comprising: using a plurality of vibratory seismic sources to generate a series of two or more seismic sweeps by activating a first set of one or more of the plurality of vibratory seismic sources to produce a first seismic sweep at time T 0 and activating a second set of one or more of the plurality of vibratory seismic sources to produce a second seismic sweep at time T 1 , wherein: each of the first and the second seismic sweeps comprises a sweep period and a listening time L; a timing of the first and the second seismic sweeps are configured such that T 1 <T 0 +S+L; the first seismic sweep has a minimum frequency Fmin and a maximum frequency Fmax; the first seismic sweep overlaps in time with the second seismic sweep; and harmonics of the second sweep overlap with a fundamental of the first sweep in an overlap region in the time-frequency domain; determining a ground force generated by at least one of the first and the second seismic sweeps; using a plurality of receivers to measure seismic signals generated by the first and the second seismic sweeps; using the determined ground force to determine properties of the harmonics of the second sweep at one or more frequencies in the overlap region; using the measured seismic signals for a plurality of frequencies adjacent to the overlap region to determine weightings applicable to the harmonics of the second sweep, wherein the weightings applicable to the harmonics of the second sweep are determined from the following equation: ρ ⁡ ( f 1 ) = ρ GF ⁡ ( f 1 ) ⁢ ρ ⁡ ( f 2 ) ρ GF ⁡ ( f 2 ) , where: ρ is a harmonic weighting for a frequency; f 1 is a frequency at which overlap between the harmonics of the second sweep and the fundamental of the first sweep in the time-frequency domain occurs, f 2 is a frequency at which there is no overlap between the harmonics of the second sweep and the fundamental of the first sweep in the time-frequency domain, and wherein f 2 is between the minimum frequency Fmin of the first seismic sweep and the maximum frequency Fmax of the first seismic sweep, and superscript GF denotes that the harmonic weighting is derived from ground-force measurements; and using the determined properties of the harmonics of the second sweep and the determined weightings applicable to the harmonics of the second sweep to attenuate harmonic noise from the measured seismic signals and produce harmonic noise attenuated seismic data.