Seismic reservoir characterization becomes challenging when reservoir thickness goes beyond the limits of seismic resolution. Geostatistical inversion techniques are being considered to overcome the resolution limitations of conventional inversion methods and to provide an intuitive understanding of subsurface uncertainty. Geostatistical inversion was applied on a highly compartmentalized area of Sapphire gas field, offshore Nile Delta, Egypt, with the aim of understanding the distribution of thin sands and their impact on reservoir connectivity. The integration of high-resolution well data with seismic partial-angle-stack volumes into geostatistical inversion has resulted in multiple elastic property realizations at the desired resolution. The multitude of inverted elastic properties are analyzed to improve reservoir characterization and reflect the inversion nonuniqueness. These property realizations are then classified into facies probability cubes and ranked based on pay sand volumes to quantify the volumetric uncertainty in static reservoir modeling. Stochastic connectivity analysis was also applied on facies models to assess the possible connected volumes. Sand connectivity analysis showed that the connected pay sand volume derived from the posterior mean of property realizations, which is analogous to deterministic inversion, is much smaller than the volumes generated by any high-frequency realization. This observation supports the role of thin interbed reservoirs in facilitating connectivity between the main sand units.