The study aims to implement a high-resolution Extended Elastic Impedance (EEI) inversion to estimate the petrophysical properties (e.g., porosity, saturation and volume of shale) from seismic and well log data. The inversion resolves the pitfall of basic EEI inversion in inverting below-tuning seismic data. The resolution, dimensionality and absolute value of basic EEI inversion are improved by employing stochastic perturbation constrained by integrated energy spectra attribute in a Bayesian Markov Chain Monte Carlo framework. A general regression neural network (GRNN) is trained to learn and memorize the relationship between the stochastically perturbed EEI and the associated well petrophysical log data. The trained GRNN is then used to predict the petrophysical properties of any given stochastic processed EEI. The proposed inversion was successfully conducted to invert the volume of shale, porosity and water saturation of a 4.0 m thick gas sand reservoir in Sarawak Basin, Malaysia. The three petrophysical geobodies were successfully built using the discovery wells cut-off values, showing that the inverted petrophysical properties satisfactorily reconstruct the well petrophysical logs with sufficient resolution and an accurate absolute value at the well site and are laterally conformable with seismic data. Inversion provides reliable petrophysical properties prediction that potentially helps further reservoir development for the study field.