Post-mineralization activities may cause difficulties in the process of ore modeling in porphyry deposits. Sungun, NW Iran, is one of the porphyry copper deposits, in which dyke intrusions have made ore modeling more complicated than expected. Among different kinds of dykes, two types were chosen and the consequent geostatistical analyses were applied on. In this study, simple directional variograms were used for extracting relevant information from dyke systems on which the Sequential Indicator Simulations were applied consequently. One hundred realizations were produced on the simulation grid considering the anisotropy characteristics and E-type map has been provided averaging all realizations. Moreover, a binary state between dyke and non-dyke environments was produced putting threshold on the E-type grid node values to discriminate the ore from barren dykes. Hole-effect models were fitted to the empirical variograms perpendicular to the dyke strike. Dimensional information was elicited from these models and the results were compared with the previously carried out geological investigations, and finally a good numerical match was found between these two sources of information.