The sandy-clayey hydrocarbon reservoirs of the Upper Paleocene and Lower Eocene located to the north of Veracruz State, Mexico, present highly complex geological and petrophysical characteristics. These reservoirs, which consist of sandstone and shale bodies within a depth interval ranging from 500 to 2000 m, were characterized statistically by means of fractal modeling and geostatistical tools. For 14 wells within an area of study of approximately 6 km², various geophysical well logs were initially edited and further analyzed to establish a correlation between logs and core data. The fractal modeling based on the R/S (rescaled range) methodology and the interpolation method by successive random additions were used to generate pseudo-well logs between observed wells. The application of geostatistical tools, sequential Gaussian simulation and exponential model variograms contributed to estimate the spatial distribution of petrophysical properties such as effective porosity (PHIE), permeability (K) and shale volume (VSH). From the analysis and correlation of the information generated in the present study, it can be said, from a general point of view, that the results not only are correlated with already reported information but also provide significant characterization elements that would be hardly obtained by means of conventional techniques.