G eothermal exploration involves a systematic study of subsurface conditions to identify areas of high geothermal potential, including surface and subsurface methods. This study combines surface and subsurface analysis in order to determine the distribution of reservoir zone and its relation to the permeability zone of the study area. The fault and fracture density method is used for zoning the value of the lineament density on the surface of a geothermal prospect area that is related to the emergence of geothermal manifestations on the surface. The magnetotelluric (MT) geophysical method has become the method of choice for exploration of subsurface geothermal resources. This method can describe the resistivity values of the main components of the geothermal system, making it effective for describing subsurface geometry including detecting reservoir zones. The microeaethquake (MEQ) method is used to detect subsurface weak zones related to geological structures, this is important in selecting permeability pathways from geological structures found in geothermal prospect areas because of their relation to the reservoir zone as a pathway for fluid to migrate from the geothermal reservoir. The results of the analysis of the three methods show that the geothermal reservoir zones are distributed in the central and southeastern parts at a depth of approximately 1000 m.