A 3D integrated saturation model was built for the Sierras Blancas Formation of the Neuquén Basin, Argentina. The saturation model was based on core, logs and seismic data. History match of reservoir pressure and well productivities were taken into account to accurately determine the gas in place and productive reservoir boundaries, specifically using 3D seismic water saturation in the gas condensate formation. The Sierras Blancas Formation is an eolian deposit. In tight, wet and diagenetically altered regions, the seismic inversion porosity and acoustic impedance based models were not adequate to describe the gas in place distribution. Further, the effective gas permeabilities in the tighter part of the reservoir are a strong function of the initial water saturation as evidenced by fewer condensate and water blocking problems of horizontal wells that navigated through low water saturation, high permeability regions. Any relationship between seismic impedance and porosity correlation degraded in areas affected by secondary diagenetic processes therefore necessitating the use of a saturation parameter. 30 vertical wells that had DT curves were selected based on their production and spatial location in order to establish a correlation between log saturation and seismic attributes. Seismic saturation cubes were generated by multiattribute seismic analysis and resampled into the simulation scale model. Log saturations were then co-kriged with the 3D seismic saturation. Water saturations obtained from the initialized simulation scale model were compared with the 2D saturation logs, the 3D seismic and the geological model scales. An objective function was defined to match the 3D seismic water saturation with the initialized simulation model water saturation. Model parameters were iterated until a satisfactory match with the initialized simulation model was obtained. By focusing the saturation match at the initialization stage, seismic, geological, petrophysical and SCAL models were ensured to be consistent prior to the full history match. Well history matching was consequently achieved much more simply and quickly. This paper presents a new detailed methodology of 3D pseudo-seismic water saturation generation, modeling and simulation used to accurately define OGIP, the productive boundaries of the reservoir, and to design trajectories for new horizontal wells.