Geothermal energy is the energy that has been obtained from the heated fluid (hot water or steam) within depth of the earth. In a geothermal system heat stored in rocks and lava located deep in the earth is transferred to the surface through a carrier fluid. Currently due to the growth in population and industry, energy shortage and the crisis of global warming that has been created by fossil fuels, geothermal energy has become an attractive energy source in the world. Geothermal reservoirs are usually located at a depth range of 1 to 5 km, so to efficiently utilize such resources an advanced prospecting method is needed to detect these deep geologic structures. The TEM (Transient Electro Magnetic) and MT (Magnetotelluric) methods are commonly used in geothermal exploration. Because it can image the low resistivity, low permeability smectite clay that caps most geothermal reservoirs, magnetotellurics (MT) is commonly used to help target geothermal wells and assess resource capacity. The methods give valuable information where the resistivity structure of the geothermal systems may be connected to temperature and other attributes that are of interest. In our country there are also the areas with high geothermal gradient that with the studies done, about 18 promising area for the utilization of this energy is detected. Sabalan geothermal field is a high potential geothermal system in Ardabil that So far, much research has been done to explore and assess its thermal potential. History feasibility of operations in this area dates back to 1350 years that the first time by an Italian company, several studies conducted in this area and Sabalan region as a first priority to the exploration of geothermal energy for power generation and other industrial applications, was introduced in the country. In 1998, the Magnetotelluric data has been collected in area in order to determine the resistivity of the superficial and deep structure by a New Zealand company. In this study modeling of data collected in 1998 is performed in order to investigate and determine the major components of geothermal system and Status of resistivity in the upper layers of geothermal reservoir. For this purpose, 11 profiles from the data approximately 20 km length in two orthogonal and aligned with the major geological structures in the southern of Meshkinshahr (mean fault in Moil area) are selected and data modeling on approximately 13 sites along each of the profiles is performed. Results of one-dimensional inverse modeling, two-dimensional inverse modeling and the final three dimensional model indicate the presence of two zones around the valley and one zone in the northeast of moil Valley with resistivity less than 11 ohm-m in the elevation range of 788 to 1588 meter above sea level, which is consistent with previous studies and geological information in the area. Finally, based on the product area of ??this region in the conductivity value as criteria rated as areas with high potential, in order of priority three optimal points were proposed.