1. A method of seismic data processing comprising: a) acquiring a first prestack seismic dataset and a second prestack seismic dataset wherein the second prestack seismic dataset has a time-lapse relationship to the first prestack seismic dataset; b) inverting the first prestack seismic dataset for determination of a first source wavelet data set and a first reflectivity parameters data set; c) obtaining a perturbation data set by using a computer to simultaneously minimize: i) an Lp-Norm between a portion of the first reflectivity parameters data set and reflectivity parameters for the second data set; ii) an Lp-Norm between the first source wavelet data set and a source wavelet for the second data set; iii) an Lp-Norm for data misfit between the first prestack seismic data set and the second prestack seismic data set.