Abstarct: Magnetotelluric method (MT) is an electromagnetic (EM) method that uses natural EM fields, generated from Earth's magnetosphere for mapping the deep subsurface structures. This method with recording the vertical and horizontal components of magnetic and electric field in earth surface, achieved the electrical conductivity distribution in earth underground. High penetration depth of EM fields in MT method has made its wide application in deep target exploration such as geothermal and hydrocarbon resources. For this reason MT method has special status among the geophysical methods. Preliminary studies on geology of Bushli region that state in South West of Nir and hot springs spa Sareyn, showed relatively good potential for geothermal resources in this region that could be a good alternative to fossil fuels and future energy supply. In this study we aim to process 60 station of magnetotelluric data, modeling and interpretation of these data. At first, we reviewed 60 MT sounding time series data, and analyzed using the methodology of resistance to noise (Robust methods) and outlier elimination, reprocess the raw data to achieve high quality of the aperient resistivity and phase data at each desired frequency. For this purpose, we used specialized software such as Mapros for processing and WinGlixnk for 1-D and 2-D smooth inversion modeling. 1-D and 2-D modelling and electric vertical sections and also horizontal resistivity maps at different depths, was diagnosed the geothermal reservoir at a depth of 2500 to 3000 meters. Also, due to the same depth resistivity map, possible expansion tank, mainly to the south and partly in the East of the region. Probably eastern part of the region due to the many faults as well as being low height compared to neighboring areas gave rise to numerous hot springs in the area