Abstarct: The find and exploration of mineral deposits requires the identification of target areas. Geophysical methods are used as indirect methods and with low cost and time for find and exploration of mineral deposits. Induction polarization (IP) is one of the most commonly used methods for the identification of sulfide and mextal minerals, especially dispersed sulfides. According to the purpose of the study, determining the zone of copper mineralization in the Sonajeel region of East Azarbaijan province, IP data modeling has been selected as the main method in this research. Sonajeel is located 17 kilometers from the city of Harris in the province of East Azarbaijan. The main stone units of the area from the old to the new are: volcanic and volcanoclastic rocks, Sonageel porphyry stock, Incheh granitoid stock, and Okuzdaghi volcanic rocks. To access the results closer to reality And a precise interpretation, the combination of the IP method and resistivity and magnetic methods are used. IP data and resistive data were collected using a pole-dipole array. The reason for choosing this array is the high penetration depth (relative to dipole - dipole array) and its lower sensitivity to telluronic and electromagnetic coupling. IP data and resistive data were collected during 6 parallel profiles in the west-east direction with 30 meter electrode intervals. Two-dimensional inverse modeling of IP and resistive data was performed using Res2dinv and ZondRes2d softwares. According to the results of data modeling, the sites with high chargeability and high resistivity values, Indicates the position of mineralization. hence, it can be concluded that copper mineralization is porphyry and scattered. In continuation, in order to determine more precisely the mineralization range, The three-dimensional inverse modeling of data was done using the Res3dinv software. magnetic collected in this region due to accompaniment, mineralization of sulfide with mineral magnetite and the presence of clay in the study area was carried out as an indirect method. Magnetometric data collection was carried out In length 19 parallel profiles with 20 m station spacing and profiles distances 50 m (excluding profile 19, located 30 meters in profile 18). After the collected data , the correction, processing and application of different magnetic filters are used, to determine the anomalys., Including filters that apply on the magnetic data for the boundary estimates, separation of a local magnetic field From regional are applied to the data; Can to filteres, RTP, Upward, Downward, Vertical derivatives, Analytical signals, Trend, Horizontal derivative, Cited. Due to the effect of geologic structure on the intensity of the magnetic field, three-dimensional modeling of the data was investigated. The process of magnetic data modeling is modeled by Mag3D software. By reviewing the results . By reviewing the results, can Found out to the presence of magnetite in the region. According to the results of modeling and interpretation of IP findings, resistivity, magnetism, and also drilling and geological information, the direct relation between copper and magnetism mineralization has been proved. In such a way that regions with high magnetic field strengths, In sections and IP and resistive models, it also shows high values. By comparing and interpreting the magnetic results, their Assessment with geological data, The existence probabilit of magnetite cany and, consequently, copper mineralization in the area proved to be .