A challenging situation in reservoir characterization arises when reservoirs occur below the vertical resolution of seismic, such as the 6 m gas-prone sands of the heterogenous Khadro Formation. This formation is the Zamzama gas field's secondary reservoir and needs detailed evaluation for optimum production of the field as the deeper primary reservoir of Pab sandstone is facing the problem of early water encroachment in its producing wells. The characterization of thin heterogeneous gas-bearing facies is hampered by the insensitivity of well-bore recording tools in precisely capturing elastic properties. Consequently, petro-elastic models (PEMs) of identified facies are generated that produce consistent elastic responses and discriminate facies in their true elastic domains. The resolution of seismic elastic properties is enhanced to illuminate the thin gas-bearing sands by adopting a Bayesian stochastic inversion process in collaboration with PEMs modeled elastic responses on a fine-scale stratigraphic grid. In addition, the Bayesian framework is used to estimate the litho-facies probability cubes to mitigate the drilling risks by integrating highly sampled posterior elastic volumes, modeled elastic logs, and geologic information of identified litho-facies. The outcomes included prospect identification at the drilled well locations that need to be perforated, while new potential zones are present for additional wells.