Time-lapse simultaneous inversion of amplitudes and time shifts constrained by pre-computed input maps

A non-transitory computer-readable storage medium containing an executable code that when executed by a computer causes the computer to perform a method for seismic monitoring including a global inversion of multi-vintage seismic data, the method comprising:

• obtaining seismic data for vintages acquired over a subsurface structure at different times, each of the vintages including at least one seismic trace;
• generating an initial model M of the subsurface structure for all the vintages, the initial model containing initial values for geophysical properties of the subsurface structure;
• defining a cost function including a difference between a reference time shift map of seismic traces between two of the vintages and a calculated time shift map of the seismic traces between the two of the vintages, wherein the calculated time shift map is based on a proposed perturbation to at least one of the geophysical properties of the model M;
• and
• using simulated annealing to minimize the cost function and to produce updated current values of the geophysical properties of the subsurface structure,
• wherein changes in the geophysical properties between the vintages are used to monitor a gas and/or oil reservoir in the subsurface structure.