The invention is a method for inferring the saturation and pressure change of a reservoir by combining the information from 4D (time-lapse) seismic and time lag data volumes (7) derived from the 4D seismic, well logs (4), and reservoir simulation results (when simulator results are available) and featuring one or more 4D well ties (1) for a quantitative 4D interpretation. The inventive method uses model-based inversion incorporating rock physics (2) at well locations (5), and is statistical-based (6) away from wells. The method thus allows integration (8) of rock physics model and reservoir simulation and honors 4D seismic change.