Rapid and accurate characterization of time-lapse seismic data is important to enable operational adjustments to be made and provide a guide for future drilling. A facies-based Bayesian inversion offers some advantages over traditional simultaneous prestack inversion, primarily avoiding the laborious construction of low-frequency models. To implement a facies-based inversion method, we adjust the model parameterization to be the ratio of monitor to baseline elastic properties. With this parameterization, the set of facies is reduced to those corresponding to specific production scenarios (production facies) that characterize expected subsurface changes between the monitor and baseline acquisitions. Production facies’ elastic properties are generally modelled through rock physics relationships. The inversion operates on the difference of the angle stacks directly and hence requires properly calibrated and registered baseline and monitor data. The result is a rapid workflow that can image changes in elastic properties accurately. We demonstrate the technique on a synthetic example, and also on field data from an operating oil sands thermal recovery project in Alberta, Canada.