A method to predict tars and tar sands in petroleum exploration and production in a subsurface volume by building a viscoelastic property model of the volume using rock behavior including tars together with seismic inversion to derive a lithology volume identifying tars.