Reservoir characterization of sand-shale sequences has always challenged geoscientists due to the presence of anisotropy in the form of shale lenses or shale layers. Water saturation and volume of shale are among the fundamental reservoir properties of interest for sand-shale intervals, and relate to the amount of fluid content and accumulating potentials of such media. This paper suggests an integrated workflow using synthetic data for the characterization of shaley-sand media based on anisotropic rock physics (T-matrix approximation) and seismic reflectivity modelling. A Bayesian inversion scheme for estimating reservoir parameters from amplitude vs. offset (AVO) data was used to obtain the information about uncertainties as well as their most likely values. The results from our workflow give reliable estimates of water saturation from AVO data at small uncertainties, provided background sand porosity values and isotropic overburden properties are known. For volume of shale, the proposed workflow provides reasonable estimates even when larger uncertainties are present in AVO data.