Madan Bozorg is an active copper mine located in NE Iran, which is a part of the very wide copper mineralization zone named Miami-Sabzevar copper belt. The main goal of this research work is the 3D model construction of the induced polarization (IP) and resistivity (Rs) data with quantifying the uncertainties using geostatistical methods and drilling. Four profiles were designed and surveyed using the CRSP array based on the boreholes. The data obtained was processed, 2D sections of IP and Rs were prepared for each profile by inverting the data, and these sections were evaluated by some exploratory boreholes in the studied area. Based on the geostatistical methods, 3D block models were constructed for the 2D IP and Rs data, and the uncertainties in the prepared models were obtained. The mineralization location was determined according to the geophysical detected anomalies. In order to check the models, some locations were proposed for drilling in the cases that the borehole data was unavailable. The drilling results indicated a high correlation between the identified anomalies from the models and mineralization in the boreholes. The results obtained show that it is possible to construct 3D models from surveyed 2D IP & Rs data with an acceptable error level. In this way, the suggested omitted drilling locations were optimized so that more potentials could be obtained for copper exploration by the least number of boreholes