Obtaining quantitative information from 4D seismic data about changes in the saturation and pressure in a reservoir due to production is very important for improved hydrocarbon recovery and management of the field. However, in heterogeneous reservoirs where the dynamic behavior varies spatially below the seismic resolution, conventional 4D interpretation methods fail to track quantitative changes in individual lithologies. The lack of information about the particular lithology to which a given 4D effect is connected can be detrimental to a general understanding of the reservoir and the forecast.