A computer implemented method of correcting time lapse seismic data for overburden and recording variations present in the data, the data being obtained in a succession of seismic surveys over a period of time from an area of interest regarding a subsurface reservoir beneath an overburden, the method comprising the computer implemented steps of: forming a measure of average amplitude of seismic energy at the reservoir in a trace in the data of a selected one of the time lapse surveys; forming a measure of corrected reservoir amplitude of seismic energy at the reservoir in the trace in the data of the selected one of the time lapse surveys; forming a correction scalar measure indicating the ratio of the measure of corrected reservoir amplitude to the average amplitude at the reservoir in the trace of the selected one of the time lapse surveys; applying the correction scalar to the trace of the selected one of the time lapse surveys to form a corrected amplitude seismic trace; and forming a record of the corrected amplitude seismic trace; wherein forming a measure of corrected reservoir amplitude of seismic energy comprises: performing a transformation of the trace from a time domain measure to a frequency domain measure to obtain an amplitude spectrum of the seismic energy at the reservoir as a function of frequency; performing a transformation of the trace from a time domain measure to a frequency domain measure to obtain an amplitude spectrum of the seismic energy at the overburden as a function of frequency; forming an average of a ratio of the amplitude spectrum at the reservoir to the amplitude spectrum at the overburden for a select range of frequencies; or, alternatively, wherein forming a measure of corrected reservoir amplitude of seismic energy comprises: performing a time frequency transform of the data in the trace; forming an average trace from a selected range of frequencies in the time frequency transform; forming a measure of average amplitude of the seismic energy at the overburden; forming a measure of average amplitude of the seismic energy at the reservoir; and forming a ratio of the average amplitude of the seismic energy at the reservoir to the measure of average amplitude of the seismic energy at the overburden; forming a measure of average amplitude of seismic energy at the reservoir in a trace in the data of a selected one of the time lapse surveys; forming a measure of corrected reservoir amplitude of seismic energy at the reservoir in the trace in the data of the selected one of the time lapse surveys; forming a correction scalar measure indicating the ratio of the measure of corrected reservoir amplitude to the average amplitude at the reservoir in the trace of the selected one of the time lapse surveys; applying the correction scalar to the trace of the selected one of the time lapse surveys to form a corrected amplitude seismic trace; and forming a record of the corrected amplitude seismic trace; wherein forming a measure of corrected reservoir amplitude of seismic energy comprises: performing a transformation of the trace from a time domain measure to a frequency domain measure to obtain an amplitude spectrum of the seismic energy at the reservoir as a function of frequency; performing a transformation of the trace from a time domain measure to a frequency domain measure to obtain an amplitude spectrum of the seismic energy at the overburden as a function of frequency; forming an average of a ratio of the amplitude spectrum at the reservoir to the amplitude spectrum at the overburden for a select range of frequencies; or, alternatively, wherein forming a measure of corrected reservoir amplitude of seismic energy comprises: performing a time frequency transform of the data in the trace; forming an average trace from a selected range of frequencies in the time frequency transform; forming a measure of average amplitude of the seismic energy at the overburden; forming a measure of average amplitude of the seismic energy at the reservoir; and forming a ratio of the average amplitude of the seismic energy at the reservoir to the measure of average amplitude of the seismic energy at the overburden; performing a transformation of the trace from a time domain measure to a frequency domain measure to obtain an amplitude spectrum of the seismic energy at the reservoir as a function of frequency; performing a transformation of the trace from a time domain measure to a frequency domain measure to obtain an amplitude spectrum of the seismic energy at the overburden as a function of frequency; forming an average of a ratio of the amplitude spectrum at the reservoir to the amplitude spectrum at the overburden for a select range of frequencies; or, alternatively, wherein forming a measure of corrected reservoir amplitude of seismic energy comprises: performing a time frequency transform of the data in the trace; forming an average trace from a selected range of frequencies in the time frequency transform; forming a measure of average amplitude of the seismic energy at the overburden; forming a measure of average amplitude of the seismic energy at the reservoir; and forming a ratio of the average amplitude of the seismic energy at the reservoir to the measure of average amplitude of the seismic energy at the overburden;