The three-dimensional porosity simulation plays a pivotal role in the exploration and prospecting of oil fields. In this study, we introduce a hybrid methodology for subsismic resolution porosity estimation, combining algorithms based on fuzzy logic and geostatistics. We demonstrate how it is possible to use simulated facies association surfaces and structural surfaces as input parameters for the fuzzy logic algorithm to derive porosity trends, subsequently refined using simple kriging with magnetic resonance log obtained from wells. The selected simulation area is the Sapinhoá field, located in the Santos Basin, Brazil, covering the age interval from 123 Ma to 113 Ma. To statistically validate the results, we employed the Pearson correlation coefficient, along with graphs comparing the measured and simulated porosity in the wells. The obtained results showed satisfactory levels of accuracy in the evaluated wells. Moreover, in areas lacking well data, the simulated porosity volumes aligned with the trends modeled using fuzzy logic.